Lines Matching refs:p_th
85 PetscReal p_th; /* Thermodynamic pressure */ member
225 const PetscReal p_th = PetscRealPart(constants[P_TH]); in f0_conduct_quadratic_v() local
226 const PetscReal rho = p_th / (t + X[0] + X[1] + 1.); in f0_conduct_quadratic_v()
237 const PetscReal p_th = PetscRealPart(constants[P_TH]); in f0_conduct_quadratic_q() local
239 …f0[0] += p_th * (S + 2. * t + 3. * X[0] * X[0] - 2. * X[0] * X[1] + X[1] * X[1]) / PetscSqr(t + X[… in f0_conduct_quadratic_q()
246 const PetscReal p_th = PetscRealPart(constants[P_TH]); in f0_conduct_quadratic_w() local
248 …f0[0] -= c_p * p_th * (S + 2. * t + 3. * X[0] * X[0] - 2. * X[0] * X[1] + X[1] * X[1]) / (t + X[0]… in f0_conduct_quadratic_w()
570 const PetscReal p_th = PetscRealPart(constants[P_TH]); in f0_conduct_pipe_v() local
572 const PetscReal rho = p_th / (X[1] * (1. - X[1]) + T_in); in f0_conduct_pipe_v()
686 const PetscReal p_th = PetscRealPart(constants[P_TH]); in f0_conduct_pipe_wiggly_v() local
689 const PetscReal rho = p_th / (X[1] * (1. - X[1]) + T_in); in f0_conduct_pipe_wiggly_v()
733 const PetscReal p_th = PetscRealPart(constants[P_TH]); in f0_conduct_q() local
737 f0[0] += -u_t[uOff[TEMP]] * S * p_th / PetscSqr(u[uOff[TEMP]]); in f0_conduct_q()
740 for (d = 0; d < dim; ++d) f0[0] += p_th / u[uOff[TEMP]] * u_x[uOff_x[VEL] + d * dim + d]; in f0_conduct_q()
743 …for (d = 0; d < dim; ++d) f0[0] -= p_th / (u[uOff[TEMP]] * u[uOff[TEMP]]) * u[uOff[VEL] + d] * u_x… in f0_conduct_q()
768 const PetscReal p_th = PetscRealPart(constants[P_TH]); in f0_conduct_v() local
769 const PetscReal rho = p_th / PetscRealPart(u[uOff[TEMP]]); in f0_conduct_v()
843 const PetscReal p_th = PetscRealPart(constants[P_TH]); in f0_conduct_w() local
847 f0[0] = c_p * S * p_th / u[uOff[TEMP]] * u_t[uOff[TEMP]]; in f0_conduct_w()
850 …for (d = 0; d < dim; ++d) f0[0] += c_p * p_th / u[uOff[TEMP]] * u[uOff[VEL] + d] * u_x[uOff_x[TEMP… in f0_conduct_w()
910 const PetscReal p_th = PetscRealPart(constants[P_TH]); in g0_conduct_qu() local
914 for (d = 0; d < dim; ++d) g0[d] = -p_th / PetscSqr(u[uOff[TEMP]]) * u_x[uOff_x[TEMP] + d]; in g0_conduct_qu()
919 const PetscReal p_th = PetscRealPart(constants[P_TH]); in g1_conduct_qu() local
923 for (d = 0; d < dim; ++d) g1[d * dim + d] = p_th / u[uOff[TEMP]]; in g1_conduct_qu()
929 const PetscReal p_th = PetscRealPart(constants[P_TH]); in g0_conduct_qT() local
933 g0[0] -= S * p_th / PetscSqr(u[uOff[TEMP]]) * u_tShift; in g0_conduct_qT()
935 g0[0] += 2.0 * S * p_th / PetscPowScalarInt(u[uOff[TEMP]], 3) * u_t[uOff[TEMP]]; in g0_conduct_qT()
937 …for (d = 0; d < dim; ++d) g0[0] += p_th / PetscSqr(u[uOff[TEMP]]) * (-u_x[uOff_x[VEL] + d * dim + … in g0_conduct_qT()
942 const PetscReal p_th = PetscRealPart(constants[P_TH]); in g1_conduct_qT() local
946 for (d = 0; d < dim; ++d) g1[d] = -p_th / PetscSqr(u[uOff[TEMP]]) * u[uOff[VEL] + d]; in g1_conduct_qT()
973 const PetscReal p_th = PetscRealPart(constants[P_TH]); in g0_conduct_vT() local
979 for (d = 0; d < dim; ++d) g0[d] -= p_th * S / PetscSqr(u[uOff[TEMP]]) * u_t[uOff[VEL] + d]; in g0_conduct_vT()
983 …for (d = 0; d < dim; ++d) g0[c] -= p_th / PetscSqr(u[uOff[TEMP]]) * u[uOff[VEL] + d] * u_x[uOff_x[… in g0_conduct_vT()
987 g0[gdir] -= p_th / PetscSqr(u[uOff[TEMP]] * F); in g0_conduct_vT()
993 const PetscReal p_th = PetscRealPart(constants[P_TH]); in g0_conduct_vu() local
998 for (d = 0; d < dim; ++d) g0[d * dim + d] = S * p_th / u[uOff[TEMP]] * u_tShift; in g0_conduct_vu()
1002 … for (d = 0; d < dim; ++d) g0[c * Nc + d] += p_th / u[uOff[TEMP]] * u_x[uOff_x[VEL] + c * Nc + d]; in g0_conduct_vu()
1008 const PetscReal p_th = PetscRealPart(constants[P_TH]); in g1_conduct_vu() local
1017 if (c == d) g1[(c * NcJ + d) * dim + e] += p_th / u[uOff[TEMP]] * u[uOff[VEL] + e]; in g1_conduct_vu()
1075 const PetscReal p_th = PetscRealPart(constants[P_TH]); in g0_conduct_wu() local
1080 for (d = 0; d < dim; ++d) g0[d] = c_p * p_th / u[uOff[TEMP]] * u_x[uOff_x[TEMP] + d]; in g0_conduct_wu()
1086 const PetscReal p_th = PetscRealPart(constants[P_TH]); in g0_conduct_wT() local
1091 g0[0] += c_p * S * p_th / u[uOff[TEMP]] * u_tShift; in g0_conduct_wT()
1093 g0[0] -= c_p * S * p_th / PetscSqr(u[uOff[TEMP]]) * u_t[uOff[TEMP]]; in g0_conduct_wT()
1095 …for (d = 0; d < dim; ++d) g0[0] -= c_p * p_th / PetscSqr(u[uOff[TEMP]]) * u[uOff[VEL] + d] * u_x[u… in g0_conduct_wT()
1100 const PetscReal p_th = PetscRealPart(constants[P_TH]); in g1_conduct_wT() local
1105 for (d = 0; d < dim; ++d) g1[d] += c_p * p_th / u[uOff[TEMP]] * u[uOff[VEL] + d]; in g1_conduct_wT()
1157 PetscCall(PetscBagRegisterReal(bag, &p->p_th, 1.0, "p_th", "Thermodynamic pressure")); in SetupParameters()
1398 constants[P_TH] = param->p_th; in SetupProblem()