Verify calculation of the fluid properties.
31 const Real thermal_exp = 2.14E-4;
34 const Real bulk_modulus = 2.0E9;
35 const Real thermal_cond = 0.6;
36 const Real entropy = 300.0;
37 const Real visc = 1.0E-3;
38 const Real density0 = 1000.0;
39 const Real pp_coef = 0.0;
41 const Real tol = REL_TOL_CONSISTENCY;
42 const Real large_tol = 10 *
tol;
48 Real h =
cv * T + p / _fp->rho_from_p_T(p, T);
50 ABS_TEST(_fp->beta_from_p_T(p, T), thermal_exp,
tol);
51 ABS_TEST(_fp->cp_from_p_T(p, T),
cp,
tol);
52 ABS_TEST(_fp->cv_from_p_T(p, T),
cv,
tol);
53 ABS_TEST(_fp->c_from_p_T(p, T), std::sqrt(bulk_modulus / _fp->rho_from_p_T(p, T)),
tol);
54 ABS_TEST(_fp->k_from_p_T(p, T), thermal_cond,
tol);
55 ABS_TEST(_fp->k_from_p_T(p, T), thermal_cond,
tol);
56 ABS_TEST(_fp->s_from_p_T(p, T), entropy,
tol);
57 ABS_TEST(_fp->rho_from_p_T(p, T), density0 * std::exp(p / bulk_modulus - thermal_exp * T),
tol);
58 ABS_TEST(_fp->e_from_p_T(p, T),
cv * T,
tol);
59 ABS_TEST(_fp->e_from_p_rho(p, 1. /
v),
cv * T,
tol);
60 ABS_TEST(_fp->e_from_v_h(
v, h),
cv * T, large_tol);
61 ABS_TEST(_fp->mu_from_p_T(p, T), visc,
tol);
62 ABS_TEST(_fp->h_from_p_T(p, T), h,
tol);
63 ABS_TEST(_fp2->h_from_p_T(p, T),
cv * T + p * pp_coef / _fp2->rho_from_p_T(p, T),
tol);
64 ABS_TEST(_fp->cp_from_v_e(
v, e),
cp,
tol);
65 ABS_TEST(_fp->cv_from_v_e(
v, e),
cv,
tol);
67 _fp->c_from_v_e(
v, e),
68 std::sqrt(bulk_modulus / _fp->rho_from_p_T(_fp->p_from_v_e(
v, e), _fp->T_from_v_e(
v, e))),
70 ABS_TEST(_fp->k_from_v_e(
v, e), thermal_cond,
tol);
71 ABS_TEST(_fp->s_from_v_e(
v, e), entropy,
tol);
72 ABS_TEST(_fp->s_from_h_p(_fp->h_from_p_T(p, T), p), entropy,
tol);
73 ABS_TEST(_fp->T_from_v_e(
v, e), T,
tol);
74 ABS_TEST(_fp->T_from_p_h(p, _fp->h_from_p_T(p, T)), T,
tol);
75 ABS_TEST(_fp->T_from_p_rho(p, 1. /
v), T,
tol);
76 ABS_TEST(_fp->p_from_v_e(
v, e),
77 bulk_modulus * (thermal_exp * _fp->T_from_v_e(
v, e) + std::log(1 / (
v * density0))),
79 ABS_TEST(_fp->mu_from_v_e(
v, e), visc,
tol);
84 v = 6.249991432791718E-4;
85 h =
cv * T + p / _fp->rho_from_p_T(p, T);
87 ABS_TEST(_fp->beta_from_p_T(p, T), thermal_exp,
tol);
88 ABS_TEST(_fp->cp_from_p_T(p, T),
cp,
tol);
89 ABS_TEST(_fp->cv_from_p_T(p, T),
cv,
tol);
90 ABS_TEST(_fp->c_from_p_T(p, T), std::sqrt(bulk_modulus / _fp->rho_from_p_T(p, T)),
tol);
91 ABS_TEST(_fp->k_from_p_T(p, T), thermal_cond,
tol);
92 ABS_TEST(_fp->k_from_p_T(p, T), thermal_cond,
tol);
93 ABS_TEST(_fp->s_from_p_T(p, T), entropy,
tol);
94 ABS_TEST(_fp->rho_from_p_T(p, T), density0 * std::exp(p / bulk_modulus - thermal_exp * T),
tol);
95 ABS_TEST(_fp->e_from_p_T(p, T),
cv * T,
tol);
96 ABS_TEST(_fp->e_from_p_rho(p, 1. /
v),
cv * T, large_tol);
97 ABS_TEST(_fp->e_from_v_h(
v, h),
cv * T, large_tol);
98 ABS_TEST(_fp->mu_from_p_T(p, T), visc,
tol);
99 ABS_TEST(_fp->h_from_p_T(p, T), h,
tol);
100 ABS_TEST(_fp2->h_from_p_T(p, T),
cv * T + p * pp_coef / _fp2->rho_from_p_T(p, T),
tol);
101 ABS_TEST(_fp->cp_from_v_e(
v, e),
cp,
tol);
102 ABS_TEST(_fp->cv_from_v_e(
v, e),
cv,
tol);
104 _fp->c_from_v_e(
v, e),
105 std::sqrt(bulk_modulus / _fp->rho_from_p_T(_fp->p_from_v_e(
v, e), _fp->T_from_v_e(
v, e))),
107 ABS_TEST(_fp->k_from_v_e(
v, e), thermal_cond,
tol);
108 ABS_TEST(_fp->s_from_v_e(
v, e), entropy,
tol);
109 ABS_TEST(_fp->s_from_h_p(_fp->h_from_p_T(p, T), p), entropy,
tol);
110 ABS_TEST(_fp->T_from_v_e(
v, e), T,
tol);
111 ABS_TEST(_fp->T_from_p_h(p, _fp->h_from_p_T(p, T)), T,
tol);
112 ABS_TEST(_fp->T_from_p_rho(p, 1. /
v), T,
tol);
113 ABS_TEST(_fp->p_from_v_e(
v, e),
114 bulk_modulus * (thermal_exp * _fp->T_from_v_e(
v, e) + std::log(1 / (
v * density0))),
116 ABS_TEST(_fp->mu_from_v_e(
v, e), visc,
tol);
static const std::string cv
static const std::string cp
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
static const std::string v