23 ABS_TEST(_fp->molarMass(), 2.01588e-3, REL_TOL_SAVED_VALUE);
31 ABS_TEST(_fp->criticalPressure(), 1.315e6, REL_TOL_SAVED_VALUE);
32 ABS_TEST(_fp->criticalTemperature(), 33.19, REL_TOL_SAVED_VALUE);
33 ABS_TEST(_fp->criticalDensity(), 31.26226704, REL_TOL_SAVED_VALUE);
41 ABS_TEST(_fp->triplePointPressure(), 7.7e3, REL_TOL_SAVED_VALUE);
42 ABS_TEST(_fp->triplePointTemperature(), 13.952, REL_TOL_SAVED_VALUE);
51 const Real tol = 10.0 * REL_TOL_EXTERNAL_VALUE;
54 REL_TEST(_fp->vaporPressure(14.0), 7.541e3,
tol);
55 REL_TEST(_fp->vaporPressure(30.0), 0.80432e6,
tol);
65 const Real tol = REL_TOL_EXTERNAL_VALUE;
66 const std::vector<Real> hc = _fp->henryCoefficients();
68 REL_TEST(hc[0], -4.73284,
tol);
69 REL_TEST(hc[1], 6.08954,
tol);
70 REL_TEST(hc[2], 6.06066,
tol);
81 const Real tol = REL_TOL_EXTERNAL_VALUE;
83 REL_TEST(_fp->k_from_rho_T(0.80844, 298.15), 186.97e-3,
tol);
84 REL_TEST(_fp->k_from_rho_T(30.0, 35.0), 71.854e-3,
tol);
85 REL_TEST(_fp->k_from_rho_T(15.879, 400.0), 248.6e-3,
tol);
87 REL_TEST(_fp->k_from_p_T(30.0e6, 400.0), 248.6e-3,
tol);
98 const Real tol = 10.0 * REL_TOL_EXTERNAL_VALUE;
100 REL_TEST(_fp->mu_from_rho_T(0.0, 298.15), 8.8997e-6,
tol);
101 REL_TEST(_fp->mu_from_rho_T(0.68911, 350.0), 9.9645e-6,
tol);
102 REL_TEST(_fp->mu_from_rho_T(6.5764, 350.0), 1.0147e-5,
tol);
104 REL_TEST(_fp->mu_from_p_T(1.0e6, 350.0), 9.9645e-6,
tol);
105 REL_TEST(_fp->mu_from_p_T(10.0e6, 350.0), 1.0147e-5,
tol);
118 const Real tol = REL_TOL_EXTERNAL_VALUE;
120 REL_TEST(_fp->rho_from_p_T(p, T), 0.86069,
tol);
121 REL_TEST(_fp->h_from_p_T(p, T), 3676.1e3,
tol);
122 REL_TEST(_fp->e_from_p_T(p, T), 2514.2e3,
tol);
123 REL_TEST(_fp->s_from_p_T(p, T), 43.025e3,
tol);
124 REL_TEST(_fp->cp_from_p_T(p, T), 14.26e3,
tol);
125 REL_TEST(_fp->cv_from_p_T(p, T), 10.112e3,
tol);
126 REL_TEST(_fp->c_from_p_T(p, T), 1283.9,
tol);
130 REL_TEST(_fp->rho_from_p_T(p, T), 0.48302,
tol);
131 REL_TEST(_fp->h_from_p_T(p, T), 6856.5e3,
tol);
132 REL_TEST(_fp->e_from_p_T(p, T), 4786.2e3,
tol);
133 REL_TEST(_fp->s_from_p_T(p, T), 51.401e3,
tol);
134 REL_TEST(_fp->cp_from_p_T(p, T), 14.52e3,
tol);
135 REL_TEST(_fp->cv_from_p_T(p, T), 10.392e3,
tol);
136 REL_TEST(_fp->c_from_p_T(p, T), 1704.1,
tol);
140 REL_TEST(_fp->rho_from_p_T(p, T), 4.6658,
tol);
141 REL_TEST(_fp->h_from_p_T(p, T), 6923.4e3,
tol);
142 REL_TEST(_fp->e_from_p_T(p, T), 4780.1e3,
tol);
143 REL_TEST(_fp->s_from_p_T(p, T), 41.892e3,
tol);
144 REL_TEST(_fp->cp_from_p_T(p, T), 14.584e3,
tol);
145 REL_TEST(_fp->cv_from_p_T(p, T), 10.436e3,
tol);
146 REL_TEST(_fp->c_from_p_T(p, T), 1764.3,
tol);
155 const Real tol = REL_TOL_DERIVATIVE;
157 const Real p = 1.0e6;
160 DERIV_TEST(_fp->rho_from_p_T, p, T,
tol);
161 DERIV_TEST(_fp->mu_from_p_T, p, T,
tol);
162 DERIV_TEST(_fp->e_from_p_T, p, T,
tol);
163 DERIV_TEST(_fp->h_from_p_T, p, T,
tol);
164 DERIV_TEST(_fp->k_from_p_T, p, T,
tol);
168 const Real drho = 1.0e-4;
169 Real rho, drho_dp, drho_dT;
170 _fp->rho_from_p_T(p, T, rho, drho_dp, drho_dT);
173 (_fp->mu_from_rho_T(rho + drho, T) - _fp->mu_from_rho_T(rho - drho, T)) / (2.0 * drho);
174 Real mu = 0.0, dmu_drho = 0.0, dmu_dT = 0.0;
175 _fp->mu_from_rho_T(rho, T, drho_dT,
mu, dmu_drho, dmu_dT);
177 ABS_TEST(
mu, _fp->mu_from_rho_T(rho, T), REL_TOL_CONSISTENCY);
178 REL_TEST(dmu_drho, dmu_drho_fd,
tol);
182 const Real dp = 1.0e1;
183 const Real dT = 1.0e-4;
184 _fp->rho_from_p_T(p, T, rho, drho_dp, drho_dT);
185 _fp->mu_from_rho_T(rho, T, drho_dT,
mu, dmu_drho, dmu_dT);
187 Real dmu_dT_fd = (_fp->mu_from_rho_T(_fp->rho_from_p_T(p, T + dT), T + dT) -
188 _fp->mu_from_rho_T(_fp->rho_from_p_T(p, T - dT), T - dT)) /
191 REL_TEST(dmu_dT, dmu_dT_fd,
tol);
193 Real dmu_dp_fd = (_fp->mu_from_p_T(p + dp, T) - _fp->mu_from_p_T(p - dp, T)) / (2.0 * dp);
195 _fp->mu_from_p_T(p, T,
mu, dmu_dp, dmu_dT);
197 ABS_TEST(
mu, _fp->mu_from_p_T(p, T), REL_TOL_CONSISTENCY);
198 REL_TEST(dmu_dp, dmu_dp_fd,
tol);
200 _fp->mu_from_p_T(p, T,
mu, dmu_dp, dmu_dT);
201 dmu_dT_fd = (_fp->mu_from_p_T(p, T + dT) - _fp->mu_from_p_T(p, T - dT)) / (2.0 * dT);
203 REL_TEST(dmu_dT, dmu_dT_fd,
tol);
212 const Real p = 1.0e6;
213 const Real T = 300.0;
214 const Real tol = REL_TOL_CONSISTENCY;
217 Real rho, drho_dp, drho_dT;
218 _fp->rho_from_p_T(p, T, rho, drho_dp, drho_dT);
220 _fp->mu_from_p_T(p, T,
mu, dmu_dp, dmu_dT);
221 Real e, de_dp, de_dT;
222 _fp->e_from_p_T(p, T, e, de_dp, de_dT);
225 Real rho2, drho2_dp, drho2_dT, mu2, dmu2_dp, dmu2_dT, e2, de2_dp, de2_dT;
226 _fp->rho_mu_from_p_T(p, T, rho2, mu2);
228 ABS_TEST(rho, rho2,
tol);
229 ABS_TEST(
mu, mu2,
tol);
231 _fp->rho_mu_from_p_T(p, T, rho2, drho2_dp, drho2_dT, mu2, dmu2_dp, dmu2_dT);
232 ABS_TEST(rho, rho2,
tol);
233 ABS_TEST(drho_dp, drho2_dp,
tol);
234 ABS_TEST(drho_dT, drho2_dT,
tol);
235 ABS_TEST(
mu, mu2,
tol);
236 ABS_TEST(dmu_dp, dmu2_dp,
tol);
237 ABS_TEST(dmu_dT, dmu2_dT,
tol);
239 _fp->rho_e_from_p_T(p, T, rho2, drho2_dp, drho2_dT, e2, de2_dp, de2_dT);
240 ABS_TEST(rho, rho2,
tol);
241 ABS_TEST(drho_dp, drho2_dp,
tol);
242 ABS_TEST(drho_dT, drho2_dT,
tol);
243 ABS_TEST(e, e2,
tol);
244 ABS_TEST(de_dp, de2_dp,
tol);
245 ABS_TEST(de_dT, de2_dT,
tol);
static const std::string mu
TEST_F(HydrogenFluidPropertiesTest, fluidName)
Test that the fluid name is correctly returned.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real