TemperaturePressureFunctionFluidPropertiesTest.C File Reference

Go to the source code of this file.

Functions
	TEST_F (TemperaturePressureFunctionFluidPropertiesTest, fluidName)
	Test that the fluid name is correctly returned.
	TEST_F (TemperaturePressureFunctionFluidPropertiesTest, properties_constant_cv)
	Verify calculation of the fluid properties with a constant cv.
	TEST_F (TemperaturePressureFunctionFluidPropertiesTest, properties_function_cp)
	Verify calculation of the fluid properties with a function cp.
	TEST_F (TemperaturePressureFunctionFluidPropertiesTest, derivatives)
	Verify calculation of the derivatives by comparing with finite differences.
	TEST_F (TemperaturePressureFunctionFluidPropertiesTest, combined)
	Verify that the methods that return multiple properties in one call return identical values as the individual methods.
	TEST_F (TemperaturePressureFunctionFluidPropertiesTest, heatCapacityDifference)
	Verify the thermodynamic relation cp - cv = T alpha^2 / (rho kappa_T), with alpha = -(1/rho) drho/dT the thermal expansion coefficient and kappa_T = (1/rho) drho/dp > 0 the isothermal compressibility, so that cp >= cv.
	TEST_F (TemperaturePressureFunctionFluidPropertiesTest, derivativesAtZeroArguments)
	Verify that the finite-difference derivatives remain finite at zero pressure and zero internal energy, where a purely relative perturbation vanishes.

Function Documentation

TEST_F() [1/7]

TEST_F	(	TemperaturePressureFunctionFluidPropertiesTest	,
		combined
	)

Verify that the methods that return multiple properties in one call return identical values as the individual methods.

Definition at line 292 of file TemperaturePressureFunctionFluidPropertiesTest.C.

{
  const Real p = 1.0e6;
  const Real T = 300.0;
  const Real tol = REL_TOL_CONSISTENCY;
 
  // Single property methods
  Real rho, drho_dp, drho_dT;
  _fp->rho_from_p_T(p, T, rho, drho_dp, drho_dT);
  Real mu, dmu_dp, dmu_dT;
  _fp->mu_from_p_T(p, T, mu, dmu_dp, dmu_dT);
  Real e, de_dp, de_dT;
  _fp->e_from_p_T(p, T, e, de_dp, de_dT);
 
  // Combined property methods
  Real rho2, drho2_dp, drho2_dT, mu2, dmu2_dp, dmu2_dT, e2, de2_dp, de2_dT;
  _fp->rho_mu_from_p_T(p, T, rho2, mu2);
 
  ABS_TEST(rho, rho2, tol);
  ABS_TEST(mu, mu2, tol);
 
  _fp->rho_mu_from_p_T(p, T, rho2, drho2_dp, drho2_dT, mu2, dmu2_dp, dmu2_dT);
  ABS_TEST(rho, rho2, tol);
  ABS_TEST(drho_dp, drho2_dp, tol);
  ABS_TEST(drho_dT, drho2_dT, tol);
  ABS_TEST(mu, mu2, tol);
  ABS_TEST(dmu_dp, dmu2_dp, tol);
  ABS_TEST(dmu_dT, dmu2_dT, tol);
 
  _fp->rho_e_from_p_T(p, T, rho2, drho2_dp, drho2_dT, e2, de2_dp, de2_dT);
  ABS_TEST(rho, rho2, tol);
  ABS_TEST(drho_dp, drho2_dp, tol);
  ABS_TEST(drho_dT, drho2_dT, tol);
  ABS_TEST(e, e2, tol);
  ABS_TEST(de_dp, de2_dp, tol);
  ABS_TEST(de_dT, de2_dT, tol);
}

TEST_F() [2/7]

TEST_F	(	TemperaturePressureFunctionFluidPropertiesTest	,
		derivatives
	)

Verify calculation of the derivatives by comparing with finite differences.

Definition at line 216 of file TemperaturePressureFunctionFluidPropertiesTest.C.

{
  const Real tol = REL_TOL_DERIVATIVE;
  // Finite difference just does not get much
  const Real large_tol = 1000 * tol;
 
  // p, T and v, e are not consistent here but that's ok
  Real p = 1.0E7;
  Real T = 10.0;
  Real e = 8.372E5;
  Real v = 1.25E-3;
 
  AD_DERIV_TEST(_fp->rho_from_p_T, p, T, tol);
  DERIV_TEST(_fp->rho_from_p_T, p, T, tol);
  DERIV_TEST(_fp->v_from_p_T, p, T, tol);
  DERIV_TEST(_fp->mu_from_p_T, p, T, tol);
  DERIV_TEST(_fp->e_from_p_T, p, T, tol);
  DERIV_TEST(_fp->h_from_p_T, p, T, tol);
  DERIV_TEST(_fp->k_from_p_T, p, T, tol);
  DERIV_TEST(_fp->cp_from_p_T, p, T, 4 * large_tol); // uses finite differencing
  DERIV_TEST(_fp->cv_from_p_T, p, T, tol);
  DERIV_TEST(_fp->cp_from_v_e, v, e, large_tol); // uses finite differencing
  DERIV_TEST(_fp->cv_from_v_e, v, e, tol);
 
  e = 80150.0458189499477;
  v = 0.00021621621621621621;
 
  AD_DERIV_TEST(_fp_cp->rho_from_p_T, p, T, tol);
  DERIV_TEST(_fp_cp->rho_from_p_T, p, T, tol);
  DERIV_TEST(_fp_cp->v_from_p_T, p, T, tol);
  DERIV_TEST(_fp_cp->mu_from_p_T, p, T, tol);
  DERIV_TEST(_fp_cp->e_from_p_T, p, T, tol);
  DERIV_TEST(_fp_cp->h_from_p_T, p, T, tol);
  DERIV_TEST(_fp_cp->k_from_p_T, p, T, tol);
  DERIV_TEST(_fp_cp->cp_from_p_T, p, T, large_tol); // uses finite differencing
  DERIV_TEST(_fp_cp->cv_from_p_T, p, T, large_tol);
  DERIV_TEST(_fp_cp->cp_from_v_e, v, e, large_tol); // uses finite differencing
  DERIV_TEST(_fp_cp->cv_from_v_e, v, e, large_tol);
 
  p = 5.0E7;
  T = 190.0;
  e = 1.6744E6;
  v = 6.25E-4;
 
  AD_DERIV_TEST(_fp->rho_from_p_T, p, T, tol);
  DERIV_TEST(_fp->rho_from_p_T, p, T, tol);
  DERIV_TEST(_fp->v_from_p_T, p, T, tol);
  DERIV_TEST(_fp->mu_from_p_T, p, T, tol);
  DERIV_TEST(_fp->e_from_p_T, p, T, tol);
  DERIV_TEST(_fp->h_from_p_T, p, T, tol);
  DERIV_TEST(_fp->k_from_p_T, p, T, tol);
  DERIV_TEST(_fp->cp_from_p_T, p, T, 2 * large_tol); // uses finite differencing
  DERIV_TEST(_fp->cv_from_p_T, p, T, tol);
  DERIV_TEST(_fp->cp_from_v_e, v, e, 2 * large_tol); // uses finite differencing
  DERIV_TEST(_fp->cv_from_v_e, v, e, tol);
 
  e = 5374151.12329934724;
  v = 5.5944055944055945e-05;
 
  AD_DERIV_TEST(_fp_cp->rho_from_p_T, p, T, tol);
  DERIV_TEST(_fp_cp->rho_from_p_T, p, T, tol);
  DERIV_TEST(_fp_cp->v_from_p_T, p, T, tol);
  DERIV_TEST(_fp_cp->mu_from_p_T, p, T, tol);
  DERIV_TEST(_fp_cp->e_from_p_T, p, T, tol);
  DERIV_TEST(_fp_cp->h_from_p_T, p, T, tol);
  DERIV_TEST(_fp_cp->k_from_p_T, p, T, tol);
  DERIV_TEST(_fp_cp->cp_from_p_T, p, T, large_tol); // uses finite differencing
  DERIV_TEST(_fp_cp->cv_from_p_T, p, T, large_tol);
  DERIV_TEST(_fp_cp->cp_from_v_e, v, e, large_tol); // uses finite differencing
  DERIV_TEST(_fp_cp->cv_from_v_e, v, e, large_tol);
}

TEST_F() [3/7]

TEST_F	(	TemperaturePressureFunctionFluidPropertiesTest	,
		derivativesAtZeroArguments
	)

Verify that the finite-difference derivatives remain finite at zero pressure and zero internal energy, where a purely relative perturbation vanishes.

Definition at line 355 of file TemperaturePressureFunctionFluidPropertiesTest.C.

{
  // e_ref = 0 at T_ref = 0, so the internal energy vanishes at T = 0 for the constant-cv object
  const Real T = 250.0;
  Real cp, dcp_dp, dcp_dT;
  _fp->cp_from_p_T(0.0, T, cp, dcp_dp, dcp_dT);
  EXPECT_TRUE(std::isfinite(dcp_dp));
  Real cv, dcv_dp, dcv_dT;
  _fp_cp->cv_from_p_T(0.0, T, cv, dcv_dp, dcv_dT);
  EXPECT_TRUE(std::isfinite(dcv_dp));
  Real e, de_dp, de_dT;
  _fp_cp->e_from_p_T(0.0, T, e, de_dp, de_dT);
  EXPECT_TRUE(std::isfinite(de_dp));
 
  // The constant-cv object reaches e = 0 at T = T_ref = 0; use a pressure that keeps the
  // density positive
  const Real p = 1e6;
  const Real v = _fp->v_from_p_T(p, T);
  Real cp_ve, dcp_dv, dcp_de;
  _fp->cp_from_v_e(v, 0.0, cp_ve, dcp_dv, dcp_de);
  EXPECT_TRUE(std::isfinite(dcp_de));
}

TEST_F() [4/7]

TEST_F	(	TemperaturePressureFunctionFluidPropertiesTest	,
		fluidName
	)

Test that the fluid name is correctly returned.

Definition at line 16 of file TemperaturePressureFunctionFluidPropertiesTest.C.

{
  EXPECT_EQ(_fp->fluidName(), "TemperaturePressureFunctionFluidProperties");
}

TEST_F() [5/7]

TEST_F	(	TemperaturePressureFunctionFluidPropertiesTest	,
		heatCapacityDifference
	)

Verify the thermodynamic relation cp - cv = T alpha^2 / (rho kappa_T), with alpha = -(1/rho) drho/dT the thermal expansion coefficient and kappa_T = (1/rho) drho/dp > 0 the isothermal compressibility, so that cp >= cv.

Definition at line 335 of file TemperaturePressureFunctionFluidPropertiesTest.C.

{
  const Real p = 8.56E7;
  const Real T = 200.0;
  const Real rho = 1400 + 2.5 * T + 32e-5 * p;
  const Real alpha = -2.5 / rho;
  const Real kappa_T = 32e-5 / rho;
  const Real cp_minus_cv = T * alpha * alpha / (rho * kappa_T);
  EXPECT_GT(cp_minus_cv, 0.0);
 
  // constant cv
  REL_TEST(_fp->cp_from_p_T(p, T), _fp->cv_from_p_T(p, T) + cp_minus_cv, REL_TOL_CONSISTENCY);
  // function cp
  REL_TEST(_fp_cp->cv_from_p_T(p, T), _fp_cp->cp_from_p_T(p, T) - cp_minus_cv, REL_TOL_CONSISTENCY);
}

TEST_F() [6/7]

TEST_F	(	TemperaturePressureFunctionFluidPropertiesTest	,
		properties_constant_cv
	)

Verify calculation of the fluid properties with a constant cv.

Definition at line 24 of file TemperaturePressureFunctionFluidPropertiesTest.C.

{
  const Real cv = 4186.0;
 
  const Real tol = REL_TOL_CONSISTENCY;
  const Real large_tol = 100 * tol;
 
  Real p = 8.56E7;
  Real T = 200.0;
  Real e = cv * T;
  Real v = 1 / (1400 + 2.5 * T + 32e-5 * p);
  Real h = cv * T + p / _fp->rho_from_p_T(p, T);
 
  // See header for expressions
  Real thermal_cond = 14 + 2e-2 * T + 3e-5 * p;
  Real visc = 1e-3 + 1e-5 * T - 3e-9 * p;
  Real density = 1 / v;
  Real alpha = -1. / density * 2.5;
  Real beta = 1. / density * 32e-5;
  Real cp = cv + MathUtils::pow(alpha, 2) * T / density / beta;
 
  // Testing the properties with a constant cv
  ABS_TEST(_fp->cp_from_p_T(p, T), cp, tol);
  ABS_TEST(_fp->cv_from_p_T(p, T), cv, tol);
  ABS_TEST(_fp->k_from_p_T(p, T), thermal_cond, tol);
  ABS_TEST(_fp->k_from_p_T(p, T), thermal_cond, tol);
  ABS_TEST(_fp->rho_from_p_T(p, T), density, tol);
  ABS_TEST(_fp->v_from_p_T(p, T), 1 / density, tol);
  ABS_TEST(_fp->e_from_p_T(p, T), cv * T, tol);
  ABS_TEST(_fp->e_from_p_rho(p, 1. / v),
           cv * T,
           large_tol); // uses a Newton solve for variable set inversion
  ABS_TEST(_fp->mu_from_p_T(p, T), visc, tol);
  ABS_TEST(_fp->h_from_p_T(p, T), h, tol);
  ABS_TEST(_fp->cp_from_v_e(v, e), cp, tol);
  ABS_TEST(_fp->cv_from_v_e(v, e), cv, tol);
  ABS_TEST(_fp->k_from_v_e(v, e), thermal_cond, tol);
  ABS_TEST(_fp->T_from_v_e(v, e), T, tol);
  ABS_TEST(_fp->T_from_p_h(p, _fp->h_from_p_T(p, T)), T, tol);
  ABS_TEST(_fp->T_from_p_rho(p, 1. / v), T, tol);
  ABS_TEST(
      _fp->p_from_v_e(v, e), p, 1.5 * large_tol); // uses a Newton solve for variable set inversion
  ABS_TEST(_fp->mu_from_v_e(v, e), visc, tol);
 
  p = 1.06841E9;
  T = 300.0;
  e = cv * T;
  v = 1 / (1400 + 2.5 * T + 32e-5 * p);
  h = cv * T + p / _fp->rho_from_p_T(p, T);
 
  // See header for expressions
  thermal_cond = 14 + 2e-2 * T + 3e-5 * p;
  visc = 1e-3 + 1e-5 * T - 3e-9 * p;
  density = 1 / v;
  alpha = -1. / density * 2.5;
  beta = 1. / density * 32e-5;
  cp = cv + MathUtils::pow(alpha, 2) * T / density / beta;
 
  ABS_TEST(_fp->cp_from_p_T(p, T), cp, tol);
  ABS_TEST(_fp->cv_from_p_T(p, T), cv, tol);
  ABS_TEST(_fp->k_from_p_T(p, T), thermal_cond, tol);
  ABS_TEST(_fp->k_from_p_T(p, T), thermal_cond, tol);
  ABS_TEST(_fp->rho_from_p_T(p, T), density, tol);
  ABS_TEST(_fp->v_from_p_T(p, T), 1 / density, tol);
  ABS_TEST(_fp->e_from_p_T(p, T), cv * T, tol);
  ABS_TEST(_fp->e_from_p_rho(p, 1. / v),
           cv * T,
           large_tol); // uses a Newton solve for variable set inversion
  ABS_TEST(_fp->mu_from_p_T(p, T), visc, tol);
  ABS_TEST(_fp->h_from_p_T(p, T), h, tol);
  ABS_TEST(_fp->cp_from_v_e(v, e), cp, tol);
  ABS_TEST(_fp->cv_from_v_e(v, e), cv, tol);
  ABS_TEST(_fp->k_from_v_e(v, e), thermal_cond, tol);
  ABS_TEST(_fp->T_from_v_e(v, e), T, tol);
  ABS_TEST(_fp->T_from_p_h(p, _fp->h_from_p_T(p, T)), T, tol);
  ABS_TEST(_fp->T_from_p_rho(p, 1. / v), T, tol);
  ABS_TEST(_fp->p_from_v_e(v, e), p, large_tol); // uses a Newton solve for variable set inversion
  ABS_TEST(_fp->mu_from_v_e(v, e), visc, tol);
}

TEST_F() [7/7]

TEST_F	(	TemperaturePressureFunctionFluidPropertiesTest	,
		properties_function_cp
	)

Verify calculation of the fluid properties with a function cp.

Definition at line 107 of file TemperaturePressureFunctionFluidPropertiesTest.C.

{
  const Real tol = REL_TOL_CONSISTENCY;
  // properties computed using numerical integral evaluation
  const Real large_tol = 1e-8;
  // properties using both numerical integration and Newton's method
  const Real newton_tol = 1e-4;
 
  {
    Real p = 8.56E7;
    Real T = 200.0;
    Real cp = 3000. + 3 * T + 5e-4 * p;
    Real e = 9219999.5394894015;
    Real v = 1 / (1400 + 2.5 * T + 32e-5 * p);
    Real h = e + p * v;
 
    // See header for expressions
    Real thermal_cond = 14 + 2e-2 * T + 3e-5 * p;
    Real visc = 1e-3 + 1e-5 * T - 3e-9 * p;
    Real density = 1 / v;
    Real alpha = -1. / density * 2.5;
    Real beta = 1. / density * 32e-5;
    Real cv = cp - MathUtils::pow(alpha, 2) * T / density / beta;
 
    ABS_TEST(_fp_cp->cp_from_p_T(p, T), cp, tol);
    ABS_TEST(_fp_cp->cv_from_p_T(p, T), cv, tol);
    ABS_TEST(_fp_cp->k_from_p_T(p, T), thermal_cond, tol);
    ABS_TEST(_fp_cp->k_from_p_T(p, T), thermal_cond, tol);
    ABS_TEST(_fp_cp->rho_from_p_T(p, T), density, tol);
    ABS_TEST(_fp_cp->v_from_p_T(p, T), 1 / density, tol);
    ABS_TEST(_fp_cp->e_from_p_T(p, T), e, tol);
    ABS_TEST(_fp_cp->e_from_p_rho(p, 1. / v), e, 10 * large_tol);
    ABS_TEST(_fp_cp->mu_from_p_T(p, T), visc, tol);
    ABS_TEST(_fp_cp->h_from_p_T(p, T), h, tol);
    ABS_TEST(_fp_cp->cp_from_v_e(v, e), cp, large_tol);
    ABS_TEST(_fp_cp->cv_from_v_e(v, e), cv, large_tol);
    ABS_TEST(_fp_cp->k_from_v_e(v, e), thermal_cond, large_tol);
    ABS_TEST(_fp_cp->T_from_v_e(v, e), T, large_tol);
    ABS_TEST(_fp_cp->T_from_p_h(p, _fp_cp->h_from_p_T(p, T)), T, large_tol);
    ABS_TEST(_fp_cp->T_from_p_rho(p, 1. / v), T, large_tol);
    ABS_TEST(_fp_cp->p_from_v_e(v, e), p, newton_tol);
    ABS_TEST(_fp_cp->mu_from_v_e(v, e), visc, large_tol);
  }
 
  {
    Real p = 1.06841E8;
    Real T = 300.0;
    Real cp = 3000. + 3 * T + 5e-4 * p;
    Real e = 17061149.3251282759;
    Real v = 1 / (1400 + 2.5 * T + 32e-5 * p);
    Real h = e + p * v;
 
    // See header for expressions
    Real thermal_cond = 14 + 2e-2 * T + 3e-5 * p;
    Real visc = 1e-3 + 1e-5 * T - 3e-9 * p;
    Real density = 1 / v;
    Real alpha = -1. / density * 2.5;
    Real beta = 1. / density * 32e-5;
    Real cv = cp - MathUtils::pow(alpha, 2) * T / density / beta;
 
    ABS_TEST(_fp_cp->cp_from_p_T(p, T), cp, tol);
    ABS_TEST(_fp_cp->cv_from_p_T(p, T), cv, tol);
    ABS_TEST(_fp_cp->k_from_p_T(p, T), thermal_cond, tol);
    ABS_TEST(_fp_cp->k_from_p_T(p, T), thermal_cond, tol);
    ABS_TEST(_fp_cp->rho_from_p_T(p, T), density, tol);
    ABS_TEST(_fp_cp->v_from_p_T(p, T), 1 / density, tol);
    ABS_TEST(_fp_cp->e_from_p_T(p, T), e, tol);
    ABS_TEST(_fp_cp->e_from_p_rho(p, 1. / v), e, tol);
    ABS_TEST(_fp_cp->mu_from_p_T(p, T), visc, tol);
    ABS_TEST(_fp_cp->h_from_p_T(p, T), h, tol);
    ABS_TEST(_fp_cp->cp_from_v_e(v, e), cp, 1.5 * large_tol);
    ABS_TEST(_fp_cp->cv_from_v_e(v, e), cv, 1.5 * large_tol);
    ABS_TEST(_fp_cp->k_from_v_e(v, e), thermal_cond, large_tol);
    ABS_TEST(_fp_cp->T_from_v_e(v, e), T, 1.5 * large_tol);
    ABS_TEST(_fp_cp->T_from_p_h(p, _fp_cp->h_from_p_T(p, T)), T, large_tol);
    ABS_TEST(_fp_cp->T_from_p_rho(p, 1. / v), T, large_tol);
    ABS_TEST(_fp_cp->p_from_v_e(v, e), p, newton_tol);
    ABS_TEST(_fp_cp->mu_from_v_e(v, e), visc, large_tol);
  }
 
  // Testing the incompressible fluid
  {
    Real p = 1.06841E8;
    Real T = 300.0;
    Real cp = 3000. + 3 * T + 5e-4 * p;
    Real e = 17061150.0000000037;
    Real v = 1. / 2000.;
    Real h = e + p * v;
    Real cv = cp;
 
    // Some of these routines used to fail with cp_function input and incompressible fluid
    ABS_TEST(_fp_incompressible->cp_from_p_T(p, T), cp, tol);
    ABS_TEST(_fp_incompressible->cv_from_p_T(p, T), cv, tol);
    ABS_TEST(_fp_incompressible->e_from_p_T(p, T), e, tol);
    // ABS_TEST(_fp_incompressible->e_from_p_rho(p, 1. / v), e, tol); // T_from_p_rho cannot be done
    // with incompressible
    ABS_TEST(_fp_incompressible->T_from_p_h(p, _fp_incompressible->h_from_p_T(p, T)), T, large_tol);
    ABS_TEST(_fp_incompressible->h_from_p_T(p, T), h, tol);
 
    // Pressure derivative is not as accurate due to missing term in e(p,T) calculation
    DERIV_TEST(_fp_incompressible->e_from_p_T, p, T, 1e-5);
    DERIV_TEST(_fp_incompressible->h_from_p_T, p, T, 1e-5);
  }
}