TabulatedPVTFluidPropertiesTest.C File Reference

Go to the source code of this file.

Functions
	TEST_F (TabulatedPVTFluidPropertiesTest, fluidName)
	Test that the fluid name and molar mass are correctly returned.
	TEST_F (TabulatedPVTFluidPropertiesTest, molarMass)
	TEST_F (TabulatedPVTFluidPropertiesTest, tableValues)
	At the tabulated pressures, density and viscosity must reproduce the table exactly.
	TEST_F (TabulatedPVTFluidPropertiesTest, interpolation)
	Between tabulated pressures, 1/B and 1/(B mu) are interpolated linearly.
	TEST_F (TabulatedPVTFluidPropertiesTest, extrapolation)
	With out_of_bounds_behavior = ignore, 1/B and 1/(B mu) are extrapolated linearly from the end intervals outside the table.
	TEST_F (TabulatedPVTFluidPropertiesTest, setToClosestBound)
	With out_of_bounds_behavior = set_to_closest_bound, the end values are used outside the table.
	TEST_F (TabulatedPVTFluidPropertiesTest, derivatives)
	Verify the derivatives against finite differences, inside the table and in the extrapolated regions.
	TEST_F (TabulatedPVTFluidPropertiesTest, ad)
	The AD versions (used by PorousFlow AD kernels) must carry the same derivatives.
	TEST_F (TabulatedPVTFluidPropertiesTest, throwException)
	With out_of_bounds_behavior = throw (the default), an exception is thrown outside the table so that the time step can be cut, but not at the end points of the table.
	TEST_F (TabulatedPVTFluidPropertiesTest, invalidSolution)
	With out_of_bounds_behavior = declare_invalid or warn_invalid, the end values are used outside the table, as with set_to_closest_bound, and no exception is thrown.
	TEST_F (TabulatedPVTFluidPropertiesTest, extrapolationError)
	Extrapolating far enough below the table makes 1/B non-positive: an exception must be thrown so that the time step can be cut.
	TEST_F (TabulatedPVTFluidPropertiesTest, nonIncreasingPressure)
	Invalid tables must be rejected with an informative error.
	TEST_F (TabulatedPVTFluidPropertiesTest, tooFewPressures)
	TEST_F (TabulatedPVTFluidPropertiesTest, wrongSizeFormationVolumeFactor)
	TEST_F (TabulatedPVTFluidPropertiesTest, wrongSizeViscosity)
	TEST_F (TabulatedPVTFluidPropertiesTest, nonPositiveFormationVolumeFactor)
	TEST_F (TabulatedPVTFluidPropertiesTest, nonPositiveViscosity)

Function Documentation

TEST_F() [1/17]

TEST_F	(	TabulatedPVTFluidPropertiesTest	,
		ad
	)

The AD versions (used by PorousFlow AD kernels) must carry the same derivatives.

Definition at line 157 of file TabulatedPVTFluidPropertiesTest.C.

{
  const Real p = 1.5e7;
  const Real T = 350.0;
 
  ADReal ad_p = p;
  Moose::derivInsert(ad_p.derivatives(), 0, 1.0);
  ADReal ad_T = T;
  Moose::derivInsert(ad_T.derivatives(), 1, 1.0);
 
  Real rho, drho_dp, drho_dT, mu, dmu_dp, dmu_dT;
  _fp->rho_from_p_T(p, T, rho, drho_dp, drho_dT);
  _fp->mu_from_p_T(p, T, mu, dmu_dp, dmu_dT);
 
  const ADReal ad_rho = _fp->rho_from_p_T(ad_p, ad_T);
  const ADReal ad_mu = _fp->mu_from_p_T(ad_p, ad_T);
 
  REL_TEST(ad_rho.value(), rho, REL_TOL_CONSISTENCY);
  REL_TEST(ad_rho.derivatives()[0], drho_dp, REL_TOL_CONSISTENCY);
  ABS_TEST(ad_rho.derivatives()[1], drho_dT, REL_TOL_CONSISTENCY);
  REL_TEST(ad_mu.value(), mu, REL_TOL_CONSISTENCY);
  REL_TEST(ad_mu.derivatives()[0], dmu_dp, REL_TOL_CONSISTENCY);
  ABS_TEST(ad_mu.derivatives()[1], dmu_dT, REL_TOL_CONSISTENCY);
}

TEST_F() [2/17]

TEST_F	(	TabulatedPVTFluidPropertiesTest	,
		derivatives
	)

Verify the derivatives against finite differences, inside the table and in the extrapolated regions.

Pressures are chosen away from the tabulated points, where the piecewise-linear interpolant has kinks.

Definition at line 128 of file TabulatedPVTFluidPropertiesTest.C.

{
  const Real T = 350.0;
  const Real tol = REL_TOL_DERIVATIVE;
 
  for (const Real p : {1.0e6, 5.0e6, 1.5e7, 2.5e7, 3.5e7})
  {
    DERIV_TEST(_fp->rho_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_bounds->rho_from_p_T, p, T, tol);
    DERIV_TEST(_fp_bounds->mu_from_p_T, p, T, tol);
    DERIV_TEST(_fp_bounds->h_from_p_T, p, T, tol);
  }
 
  // Inside the table, the default (throw) object behaves like the others
  for (const Real p : {5.0e6, 1.5e7, 2.5e7})
  {
    DERIV_TEST(_fp_throw->rho_from_p_T, p, T, tol);
    DERIV_TEST(_fp_throw->mu_from_p_T, p, T, tol);
    DERIV_TEST(_fp_throw->h_from_p_T, p, T, tol);
  }
}

TEST_F() [3/17]

TEST_F	(	TabulatedPVTFluidPropertiesTest	,
		extrapolation
	)

With out_of_bounds_behavior = ignore, 1/B and 1/(B mu) are extrapolated linearly from the end intervals outside the table.

Definition at line 75 of file TabulatedPVTFluidPropertiesTest.C.

{
  const Real T = 350.0;
 
  // Above the table: half an interval beyond the last tabulated pressure
  {
    const Real p = 3.5e7;
    const Real inv_B = 1.0 / _B[3] + 0.5 * (1.0 / _B[3] - 1.0 / _B[2]);
    const Real inv_B_mu =
        1.0 / (_B[3] * _mu[3]) + 0.5 * (1.0 / (_B[3] * _mu[3]) - 1.0 / (_B[2] * _mu[2]));
    REL_TEST(_fp->rho_from_p_T(p, T), _rho_sc * inv_B, REL_TOL_CONSISTENCY);
    REL_TEST(_fp->mu_from_p_T(p, T), inv_B / inv_B_mu, REL_TOL_CONSISTENCY);
  }
 
  // Below the table: an eighth of the first interval below the first tabulated pressure
  {
    const Real p = 1.0e6;
    const Real inv_B = 1.0 / _B[0] - 0.125 * (1.0 / _B[1] - 1.0 / _B[0]);
    const Real inv_B_mu =
        1.0 / (_B[0] * _mu[0]) - 0.125 * (1.0 / (_B[1] * _mu[1]) - 1.0 / (_B[0] * _mu[0]));
    REL_TEST(_fp->rho_from_p_T(p, T), _rho_sc * inv_B, REL_TOL_CONSISTENCY);
    REL_TEST(_fp->mu_from_p_T(p, T), inv_B / inv_B_mu, REL_TOL_CONSISTENCY);
  }
}

TEST_F() [4/17]

TEST_F	(	TabulatedPVTFluidPropertiesTest	,
		extrapolationError
	)

Extrapolating far enough below the table makes 1/B non-positive: an exception must be thrown so that the time step can be cut.

Definition at line 242 of file TabulatedPVTFluidPropertiesTest.C.

{
  // 1/B reaches zero at about -7E8 Pa for this table
  EXPECT_THROW(_fp->rho_from_p_T(-1.0e9, 350.0), MooseException);
  EXPECT_THROW(_fp->mu_from_p_T(-1.0e9, 350.0), MooseException);
}

TEST_F() [5/17]

TEST_F	(	TabulatedPVTFluidPropertiesTest	,
		fluidName
	)

Test that the fluid name and molar mass are correctly returned.

Definition at line 18 of file TabulatedPVTFluidPropertiesTest.C.

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

TEST_F() [6/17]

TEST_F	(	TabulatedPVTFluidPropertiesTest	,
		interpolation
	)

Between tabulated pressures, 1/B and 1/(B mu) are interpolated linearly.

Definition at line 45 of file TabulatedPVTFluidPropertiesTest.C.

{
  const Real T = 350.0;
  const Real p = 1.5e7; // midpoint of the second interval
 
  const Real inv_B = 0.5 * (1.0 / _B[1] + 1.0 / _B[2]);
  const Real inv_B_mu = 0.5 * (1.0 / (_B[1] * _mu[1]) + 1.0 / (_B[2] * _mu[2]));
  const Real rho = _rho_sc * inv_B;
  const Real mu = inv_B / inv_B_mu;
 
  REL_TEST(_fp->rho_from_p_T(p, T), rho, REL_TOL_CONSISTENCY);
  REL_TEST(_fp->mu_from_p_T(p, T), mu, REL_TOL_CONSISTENCY);
  REL_TEST(_fp->e_from_p_T(p, T), _cv * T, REL_TOL_CONSISTENCY);
  REL_TEST(_fp->h_from_p_T(p, T), _cv * T + p / rho, REL_TOL_CONSISTENCY);
 
  // Saved values, also used in the regression test fluids/tabulated_pvt.i
  REL_TEST(_fp->rho_from_p_T(p, T), 814.20313790256, REL_TOL_SAVED_VALUE);
  REL_TEST(_fp->mu_from_p_T(p, T), 1.1341216986149e-3, REL_TOL_SAVED_VALUE);
 
  // rho_mu_from_p_T (used by PorousFlow) must be consistent with the individual methods
  Real rho2, mu2;
  _fp->rho_mu_from_p_T(p, T, rho2, mu2);
  REL_TEST(rho2, rho, REL_TOL_CONSISTENCY);
  REL_TEST(mu2, mu, REL_TOL_CONSISTENCY);
}

TEST_F() [7/17]

TEST_F	(	TabulatedPVTFluidPropertiesTest	,
		invalidSolution
	)

With out_of_bounds_behavior = declare_invalid or warn_invalid, the end values are used outside the table, as with set_to_closest_bound, and no exception is thrown.

Definition at line 216 of file TabulatedPVTFluidPropertiesTest.C.

{
  const Real T = 350.0;
 
  for (const auto * fp : {_fp_declare_invalid, _fp_warn_invalid})
  {
    Real rho, drho_dp, drho_dT, mu, dmu_dp, dmu_dT;
 
    fp->rho_from_p_T(3.5e7, T, rho, drho_dp, drho_dT);
    fp->mu_from_p_T(3.5e7, T, mu, dmu_dp, dmu_dT);
    REL_TEST(rho, _rho_sc / _B[3], REL_TOL_CONSISTENCY);
    REL_TEST(mu, _mu[3], REL_TOL_CONSISTENCY);
    ABS_TEST(drho_dp, 0.0, REL_TOL_CONSISTENCY);
    ABS_TEST(dmu_dp, 0.0, REL_TOL_CONSISTENCY);
 
    fp->rho_from_p_T(1.0e6, T, rho, drho_dp, drho_dT);
    fp->mu_from_p_T(1.0e6, T, mu, dmu_dp, dmu_dT);
    REL_TEST(rho, _rho_sc / _B[0], REL_TOL_CONSISTENCY);
    REL_TEST(mu, _mu[0], REL_TOL_CONSISTENCY);
  }
}

TEST_F() [8/17]

TEST_F	(	TabulatedPVTFluidPropertiesTest	,
		molarMass
	)

Definition at line 20 of file TabulatedPVTFluidPropertiesTest.C.

{
  ABS_TEST(_fp->molarMass(), _molar_mass, REL_TOL_SAVED_VALUE);
}

TEST_F() [9/17]

TEST_F	(	TabulatedPVTFluidPropertiesTest	,
		nonIncreasingPressure
	)

Invalid tables must be rejected with an informative error.

Each test constructs only one invalid object: an object whose constructor throws is not removed from the application's interface registry, so constructing a second object at the same address in the same test would trigger the "Interface already registered" assertion in debug builds

Definition at line 255 of file TabulatedPVTFluidPropertiesTest.C.

{
  InputParameters params = validFPParams();
  params.set<std::vector<Real>>("pressure") = {2.0e6, 1.0e7, 1.0e7, 3.0e7};
  EXPECT_MOOSEERROR_MSG_CONTAINS(
      _fe_problem->addUserObject("TabulatedPVTFluidProperties", "invalid_fp", params),
      "Tabulated pressures must be strictly increasing");
}

TEST_F() [10/17]

TEST_F	(	TabulatedPVTFluidPropertiesTest	,
		nonPositiveFormationVolumeFactor
	)

Definition at line 293 of file TabulatedPVTFluidPropertiesTest.C.

{
  InputParameters params = validFPParams();
  params.set<std::vector<Real>>("formation_volume_factor") = {1.062, 0.0, 1.038, 1.028};
  EXPECT_MOOSEERROR_MSG_CONTAINS(
      _fe_problem->addUserObject("TabulatedPVTFluidProperties", "invalid_fp", params),
      "All entries must be positive");
}

TEST_F() [11/17]

TEST_F	(	TabulatedPVTFluidPropertiesTest	,
		nonPositiveViscosity
	)

Definition at line 302 of file TabulatedPVTFluidPropertiesTest.C.

{
  InputParameters params = validFPParams();
  params.set<std::vector<Real>>("viscosity") = {1.04e-3, -1.10e-3, 1.17e-3, 1.24e-3};
  EXPECT_MOOSEERROR_MSG_CONTAINS(
      _fe_problem->addUserObject("TabulatedPVTFluidProperties", "invalid_fp", params),
      "All entries must be positive");
}

TEST_F() [12/17]

TEST_F	(	TabulatedPVTFluidPropertiesTest	,
		setToClosestBound
	)

With out_of_bounds_behavior = set_to_closest_bound, the end values are used outside the table.

Definition at line 103 of file TabulatedPVTFluidPropertiesTest.C.

{
  const Real T = 350.0;
  Real rho, drho_dp, drho_dT, mu, dmu_dp, dmu_dT;
 
  _fp_bounds->rho_from_p_T(3.5e7, T, rho, drho_dp, drho_dT);
  _fp_bounds->mu_from_p_T(3.5e7, T, mu, dmu_dp, dmu_dT);
  REL_TEST(rho, _rho_sc / _B[3], REL_TOL_CONSISTENCY);
  REL_TEST(mu, _mu[3], REL_TOL_CONSISTENCY);
  ABS_TEST(drho_dp, 0.0, REL_TOL_CONSISTENCY);
  ABS_TEST(dmu_dp, 0.0, REL_TOL_CONSISTENCY);
 
  _fp_bounds->rho_from_p_T(1.0e6, T, rho, drho_dp, drho_dT);
  _fp_bounds->mu_from_p_T(1.0e6, T, mu, dmu_dp, dmu_dT);
  REL_TEST(rho, _rho_sc / _B[0], REL_TOL_CONSISTENCY);
  REL_TEST(mu, _mu[0], REL_TOL_CONSISTENCY);
  ABS_TEST(drho_dp, 0.0, REL_TOL_CONSISTENCY);
  ABS_TEST(dmu_dp, 0.0, REL_TOL_CONSISTENCY);
}

TEST_F() [13/17]

TEST_F	(	TabulatedPVTFluidPropertiesTest	,
		tableValues
	)

At the tabulated pressures, density and viscosity must reproduce the table exactly.

Definition at line 28 of file TabulatedPVTFluidPropertiesTest.C.

{
  const Real T = 350.0;
  for (std::size_t i = 0; i < _p.size(); ++i)
  {
    REL_TEST(_fp->rho_from_p_T(_p[i], T), _rho_sc / _B[i], REL_TOL_CONSISTENCY);
    REL_TEST(_fp->mu_from_p_T(_p[i], T), _mu[i], REL_TOL_CONSISTENCY);
    REL_TEST(_fp_bounds->rho_from_p_T(_p[i], T), _rho_sc / _B[i], REL_TOL_CONSISTENCY);
    REL_TEST(_fp_bounds->mu_from_p_T(_p[i], T), _mu[i], REL_TOL_CONSISTENCY);
    REL_TEST(_fp_throw->rho_from_p_T(_p[i], T), _rho_sc / _B[i], REL_TOL_CONSISTENCY);
    REL_TEST(_fp_throw->mu_from_p_T(_p[i], T), _mu[i], REL_TOL_CONSISTENCY);
  }
}

TEST_F() [14/17]

TEST_F	(	TabulatedPVTFluidPropertiesTest	,
		throwException
	)

With out_of_bounds_behavior = throw (the default), an exception is thrown outside the table so that the time step can be cut, but not at the end points of the table.

Definition at line 186 of file TabulatedPVTFluidPropertiesTest.C.

{
  const Real T = 350.0;
 
  EXPECT_THROW(_fp_throw->rho_from_p_T(3.5e7, T), MooseException);
  EXPECT_THROW(_fp_throw->mu_from_p_T(3.5e7, T), MooseException);
  EXPECT_THROW(_fp_throw->h_from_p_T(3.5e7, T), MooseException);
  EXPECT_THROW(_fp_throw->rho_from_p_T(1.0e6, T), MooseException);
  EXPECT_THROW(_fp_throw->mu_from_p_T(1.0e6, T), MooseException);
 
  EXPECT_NO_THROW(_fp_throw->rho_from_p_T(_p.front(), T));
  EXPECT_NO_THROW(_fp_throw->rho_from_p_T(_p.back(), T));
 
  try
  {
    _fp_throw->rho_from_p_T(3.5e7, T);
    FAIL() << "missing expected exception";
  }
  catch (const MooseException & err)
  {
    EXPECT_NE(std::string(err.what()).find("is outside the range of tabulated pressure"),
              std::string::npos)
        << err.what();
  }
}

TEST_F() [15/17]

TEST_F	(	TabulatedPVTFluidPropertiesTest	,
		tooFewPressures
	)

Definition at line 264 of file TabulatedPVTFluidPropertiesTest.C.

{
  InputParameters params = validFPParams();
  params.set<std::vector<Real>>("pressure") = {2.0e6};
  params.set<std::vector<Real>>("formation_volume_factor") = {1.062};
  params.set<std::vector<Real>>("viscosity") = {1.04e-3};
  EXPECT_MOOSEERROR_MSG_CONTAINS(
      _fe_problem->addUserObject("TabulatedPVTFluidProperties", "invalid_fp", params),
      "At least two tabulated pressures are required");
}

TEST_F() [16/17]

TEST_F	(	TabulatedPVTFluidPropertiesTest	,
		wrongSizeFormationVolumeFactor
	)

Definition at line 275 of file TabulatedPVTFluidPropertiesTest.C.

{
  InputParameters params = validFPParams();
  params.set<std::vector<Real>>("formation_volume_factor") = {1.062, 1.050, 1.038};
  EXPECT_MOOSEERROR_MSG_CONTAINS(
      _fe_problem->addUserObject("TabulatedPVTFluidProperties", "invalid_fp", params),
      "Must have the same number of entries as 'pressure'");
}

TEST_F() [17/17]

TEST_F	(	TabulatedPVTFluidPropertiesTest	,
		wrongSizeViscosity
	)

Definition at line 284 of file TabulatedPVTFluidPropertiesTest.C.

{
  InputParameters params = validFPParams();
  params.set<std::vector<Real>>("viscosity") = {1.04e-3, 1.10e-3};
  EXPECT_MOOSEERROR_MSG_CONTAINS(
      _fe_problem->addUserObject("TabulatedPVTFluidProperties", "invalid_fp", params),
      "Must have the same number of entries as 'pressure'");
}