PorousFlowBrineCO2Test.C File Reference

Go to the source code of this file.

Functions
	TEST_F (PorousFlowBrineCO2Test, name)
	Verify that the correct name is supplied. More...
	TEST_F (PorousFlowBrineCO2Test, indices)
	Verify that the correct values for the fluid phase and component indices are supplied. More...
	TEST_F (PorousFlowBrineCO2Test, equilibriumConstants)
	TEST_F (PorousFlowBrineCO2Test, fugacityCoefficients)
	TEST_F (PorousFlowBrineCO2Test, activityCoefficients)
	TEST_F (PorousFlowBrineCO2Test, activityCoefficientCO2Brine)
	TEST_F (PorousFlowBrineCO2Test, partialDensity)
	TEST_F (PorousFlowBrineCO2Test, equilibriumMassFraction)
	TEST_F (PorousFlowBrineCO2Test, MassFraction)
	TEST_F (PorousFlowBrineCO2Test, gasProperties)
	TEST_F (PorousFlowBrineCO2Test, liquidProperties)
	TEST_F (PorousFlowBrineCO2Test, saturation)
	TEST_F (PorousFlowBrineCO2Test, totalMassFraction)
	TEST_F (PorousFlowBrineCO2Test, henryConstant)
	TEST_F (PorousFlowBrineCO2Test, enthalpyOfDissolutionGas)
	TEST_F (PorousFlowBrineCO2Test, enthalpyOfDissolution)
	TEST_F (PorousFlowBrineCO2Test, solveEquilibriumMoleFractionHighTemp)
	Test iterative calculation of equilibrium mass fractions in the elevated temperature regime. More...
	TEST_F (PorousFlowBrineCO2Test, equilibriumMoleFractions)
	Test calculation of equilibrium mole fractions over entire temperature range. More...

Function Documentation

TEST_F() [1/18]

TEST_F	(	PorousFlowBrineCO2Test	,
		name
	)

Verify that the correct name is supplied.

Definition at line 16 of file PorousFlowBrineCO2Test.C.

{ EXPECT_EQ("brine-co2", _fs->fluidStateName()); }

TEST_F() [2/18]

TEST_F	(	PorousFlowBrineCO2Test	,
		indices
	)

Verify that the correct values for the fluid phase and component indices are supplied.

Definition at line 21 of file PorousFlowBrineCO2Test.C.

{
  EXPECT_EQ((unsigned int)2, _fs->numPhases());
  EXPECT_EQ((unsigned int)3, _fs->numComponents());
  EXPECT_EQ((unsigned int)0, _fs->aqueousPhaseIndex());
  EXPECT_EQ((unsigned int)1, _fs->gasPhaseIndex());
  EXPECT_EQ((unsigned int)0, _fs->aqueousComponentIndex());
  EXPECT_EQ((unsigned int)1, _fs->gasComponentIndex());
  EXPECT_EQ((unsigned int)2, _fs->saltComponentIndex());
}

TEST_F() [3/18]

TEST_F	(	PorousFlowBrineCO2Test	,
		equilibriumConstants
	)

Definition at line 35 of file PorousFlowBrineCO2Test.C.

{
  ADReal T = 350.0;
  Moose::derivInsert(T.derivatives(), _Tidx, 1.0);

  const Real dT = 1.0e-6;

  ADReal K0H2O = _fs->equilibriumConstantH2O(T);
  ADReal K0CO2 = _fs->equilibriumConstantCO2(T);

  ABS_TEST(K0H2O.value(), 0.412597711705, 1.0e-10);
  ABS_TEST(K0CO2.value(), 74.0435888596, 1.0e-10);

  Real K0H2O_2 = _fs->equilibriumConstantH2O(T + dT).value();
  Real K0CO2_2 = _fs->equilibriumConstantCO2(T + dT).value();

  Real dK0H2O_dT_fd = (K0H2O_2 - K0H2O.value()) / dT;
  Real dK0CO2_dT_fd = (K0CO2_2 - K0CO2.value()) / dT;
  REL_TEST(K0H2O.derivatives()[_Tidx], dK0H2O_dT_fd, 1.0e-6);
  REL_TEST(K0CO2.derivatives()[_Tidx], dK0CO2_dT_fd, 1.0e-6);

  // Test the high temperature formulation
  T = 450.0;
  Moose::derivInsert(T.derivatives(), _Tidx, 1.0);

  K0H2O = _fs->equilibriumConstantH2O(T);
  K0CO2 = _fs->equilibriumConstantCO2(T);

  ABS_TEST(K0H2O.value(), 8.75944517916, 1.0e-10);
  ABS_TEST(K0CO2.value(), 105.013867434, 1.0e-10);

  K0H2O_2 = _fs->equilibriumConstantH2O(T + dT).value();
  K0CO2_2 = _fs->equilibriumConstantCO2(T + dT).value();

  dK0H2O_dT_fd = (K0H2O_2 - K0H2O.value()) / dT;
  dK0CO2_dT_fd = (K0CO2_2 - K0CO2.value()) / dT;
  REL_TEST(K0H2O.derivatives()[_Tidx], dK0H2O_dT_fd, 1.0e-6);
  REL_TEST(K0CO2.derivatives()[_Tidx], dK0CO2_dT_fd, 1.0e-6);
}

TEST_F() [4/18]

TEST_F	(	PorousFlowBrineCO2Test	,
		fugacityCoefficients
	)

Definition at line 78 of file PorousFlowBrineCO2Test.C.

{
  // Test the low temperature formulation
  ADReal p = 40.0e6;
  Moose::derivInsert(p.derivatives(), _pidx, 1.0);

  ADReal T = 350.0;
  Moose::derivInsert(T.derivatives(), _Tidx, 1.0);

  ADReal Xnacl = 0.0;
  Moose::derivInsert(Xnacl.derivatives(), _Xidx, 1.0);

  ADReal co2_density = _co2_fp->rho_from_p_T(p, T);

  ADReal phiH2O, phiCO2;
  _fs->fugacityCoefficientsLowTemp(p, T, co2_density, phiCO2, phiH2O);
  ABS_TEST(phiCO2.value(), 0.401939386415, 1.0e-8);
  ABS_TEST(phiH2O.value(), 0.0898968578757, 1.0e-8);

  // Test the high temperature formulation
  T = 423.15;

  ADReal xco2, yh2o;
  co2_density = _co2_fp->rho_from_p_T(p, T);

  _fs->solveEquilibriumMoleFractionHighTemp(
      p.value(), T.value(), Xnacl.value(), co2_density.value(), xco2.value(), yh2o.value());
  phiH2O = _fs->fugacityCoefficientH2OHighTemp(p, T, co2_density, xco2, yh2o);
  phiCO2 = _fs->fugacityCoefficientCO2HighTemp(p, T, co2_density, xco2, yh2o);

  ABS_TEST(phiH2O.value(), 0.156303437579, 1.0e-8);
  ABS_TEST(phiCO2.value(), 0.641936639599, 1.0e-8);

  // Test that the same results are returned in fugacityCoefficientsHighTemp()
  ADReal phiCO2_2, phiH2O_2;
  _fs->fugacityCoefficientsHighTemp(p, T, co2_density, xco2, yh2o, phiCO2_2, phiH2O_2);

  ABS_TEST(phiH2O, phiH2O_2, 1.0e-12);
  ABS_TEST(phiCO2, phiCO2_2, 1.0e-12);
}

TEST_F() [5/18]

TEST_F	(	PorousFlowBrineCO2Test	,
		activityCoefficients
	)

Definition at line 122 of file PorousFlowBrineCO2Test.C.

{
  ADReal xco2 = 0.01;
  ADReal T = 350.0;

  ADReal gammaH2O = _fs->activityCoefficientH2O(T, xco2);
  ADReal gammaCO2 = _fs->activityCoefficientCO2(T, xco2);

  ABS_TEST(gammaH2O.value(), 1.0, 1.0e-10);
  ABS_TEST(gammaCO2.value(), 1.0, 1.0e-10);

  T = 450.0;

  gammaH2O = _fs->activityCoefficientH2O(T, xco2);
  gammaCO2 = _fs->activityCoefficientCO2(T, xco2);

  ABS_TEST(gammaH2O.value(), 1.00022113664, 1.0e-10);
  ABS_TEST(gammaCO2.value(), 0.956736800255, 1.0e-10);
}

TEST_F() [6/18]

TEST_F	(	PorousFlowBrineCO2Test	,
		activityCoefficientCO2Brine
	)

Definition at line 146 of file PorousFlowBrineCO2Test.C.

{
  ADReal p = 10.0e6;
  Moose::derivInsert(p.derivatives(), _pidx, 1.0);

  ADReal T = 350.0;
  Moose::derivInsert(T.derivatives(), _Tidx, 1.0);

  ADReal Xnacl = 0.1;
  Moose::derivInsert(Xnacl.derivatives(), _Xidx, 1.0);

  const Real dp = 1.0e-1;
  const Real dT = 1.0e-6;
  const Real dx = 1.0e-8;

  // Low temperature regime
  ADReal gamma = _fs->activityCoefficient(p, T, Xnacl);
  ABS_TEST(gamma.value(), 1.43276649338, 1.0e-8);

  ADReal gamma_2 = _fs->activityCoefficient(p + dp, T, Xnacl);

  Real dgamma_dp_fd = (gamma_2.value() - gamma.value()) / dp;
  REL_TEST(gamma.derivatives()[_pidx], dgamma_dp_fd, 1.0e-6);

  gamma_2 = _fs->activityCoefficient(p, T + dT, Xnacl);

  Real dgamma_dT_fd = (gamma_2.value() - gamma.value()) / dT;
  REL_TEST(gamma.derivatives()[_Tidx], dgamma_dT_fd, 1.0e-6);

  gamma_2 = _fs->activityCoefficient(p, T, Xnacl + dx);

  Real dgamma_dX_fd = (gamma_2.value() - gamma.value()) / dx;
  REL_TEST(gamma.derivatives()[_Xidx], dgamma_dX_fd, 1.0e-6);

  // High temperature regime
  T = 523.15;
  Moose::derivInsert(T.derivatives(), _Tidx, 1.0);

  gamma = _fs->activityCoefficientHighTemp(T, Xnacl);
  ABS_TEST(gamma.value(), 1.50047006243, 1.0e-8);

  gamma_2 = _fs->activityCoefficientHighTemp(T + dT, Xnacl);
  dgamma_dT_fd = (gamma_2.value() - gamma.value()) / dT;
  REL_TEST(gamma.derivatives()[_Tidx], dgamma_dT_fd, 1.0e-6);

  gamma_2 = _fs->activityCoefficientHighTemp(T, Xnacl + dx);
  dgamma_dX_fd = (gamma_2.value() - gamma.value()) / dx;
  REL_TEST(gamma.derivatives()[_Xidx], dgamma_dX_fd, 1.0e-6);

  // Check that both formulations return gamma = 1 for Xnacl = 0
  Xnacl = 0.0;
  T = 350.0;
  gamma = _fs->activityCoefficient(p, T, Xnacl);
  ABS_TEST(gamma.value(), 1.0, 1.0e-12);

  T = 523.15;
  gamma = _fs->activityCoefficientHighTemp(T, Xnacl);
  ABS_TEST(gamma.value(), 1.0, 1.0e-12);
}

TEST_F() [7/18]

TEST_F	(	PorousFlowBrineCO2Test	,
		partialDensity
	)

Definition at line 209 of file PorousFlowBrineCO2Test.C.

{
  ADReal T = 473.15;
  Moose::derivInsert(T.derivatives(), _Tidx, 1.0);

  const Real dT = 1.0e-6;

  ADReal partial_density = _fs->partialDensityCO2(T);
  ABS_TEST(partial_density.value(), 893.332, 1.0e-3);

  ADReal partial_density_2 = _fs->partialDensityCO2(T + dT);

  Real dpartial_density_dT_fd = (partial_density_2.value() - partial_density.value()) / dT;
  REL_TEST(partial_density.derivatives()[_Tidx], dpartial_density_dT_fd, 1.0e-6);
}

TEST_F() [8/18]

TEST_F	(	PorousFlowBrineCO2Test	,
		equilibriumMassFraction
	)

Definition at line 228 of file PorousFlowBrineCO2Test.C.

{
  // Low temperature regime
  ADReal p = 1.0e6;
  Moose::derivInsert(p.derivatives(), _pidx, 1.0);

  ADReal T = 350.0;
  Moose::derivInsert(T.derivatives(), _Tidx, 1.0);

  ADReal Xnacl = 0.1;
  Moose::derivInsert(Xnacl.derivatives(), _Xidx, 1.0);

  const Real dp = 1.0e-2;
  const Real dT = 1.0e-6;
  const Real dx = 1.0e-8;

  ADReal X, Y;
  _fs->equilibriumMassFractions(p, T, Xnacl, X, Y);

  ABS_TEST(X.value(), 0.0035573020148, 1.0e-10);
  ABS_TEST(Y.value(), 0.0171977397214, 1.0e-10);

  // Derivative wrt pressure
  ADReal X1, Y1, X2, Y2;
  _fs->equilibriumMassFractions(p - dp, T, Xnacl, X1, Y1);
  _fs->equilibriumMassFractions(p + dp, T, Xnacl, X2, Y2);

  Real dX_dp_fd = (X2.value() - X1.value()) / (2.0 * dp);
  Real dY_dp_fd = (Y2.value() - Y1.value()) / (2.0 * dp);

  REL_TEST(X.derivatives()[_pidx], dX_dp_fd, 1.0e-6);
  REL_TEST(Y.derivatives()[_pidx], dY_dp_fd, 1.0e-6);

  // Derivative wrt temperature
  _fs->equilibriumMassFractions(p, T - dT, Xnacl, X1, Y1);
  _fs->equilibriumMassFractions(p, T + dT, Xnacl, X2, Y2);

  Real dX_dT_fd = (X2.value() - X1.value()) / (2.0 * dT);
  Real dY_dT_fd = (Y2.value() - Y1.value()) / (2.0 * dT);

  REL_TEST(X.derivatives()[_Tidx], dX_dT_fd, 1.0e-6);
  REL_TEST(Y.derivatives()[_Tidx], dY_dT_fd, 1.0e-6);

  // Derivative wrt salt mass fraction
  _fs->equilibriumMassFractions(p, T, Xnacl - dx, X1, Y1);
  _fs->equilibriumMassFractions(p, T, Xnacl + dx, X2, Y2);

  Real dX_dX_fd = (X2.value() - X1.value()) / (2.0 * dx);
  Real dY_dX_fd = (Y2.value() - Y1).value() / (2.0 * dx);

  REL_TEST(X.derivatives()[_Xidx], dX_dX_fd, 1.0e-6);
  REL_TEST(Y.derivatives()[_Xidx], dY_dX_fd, 1.0e-6);

  // High temperature regime
  p = 10.0e6;
  Moose::derivInsert(p.derivatives(), _pidx, 1.0);

  T = 525.15;
  Moose::derivInsert(T.derivatives(), _Tidx, 1.0);

  _fs->equilibriumMassFractions(p, T, Xnacl, X, Y);

  ABS_TEST(X.value(), 0.016299479086, 1.0e-10);
  ABS_TEST(Y.value(), 0.249471400766, 1.0e-10);
}

TEST_F() [9/18]

TEST_F	(	PorousFlowBrineCO2Test	,
		MassFraction
	)

Definition at line 298 of file PorousFlowBrineCO2Test.C.

{
  ADReal p = 1.0e6;
  Moose::derivInsert(p.derivatives(), _pidx, 1.0);

  ADReal T = 350.0;
  Moose::derivInsert(T.derivatives(), _Tidx, 1.0);

  ADReal Xnacl = 0.1;
  Moose::derivInsert(Xnacl.derivatives(), _Xidx, 1.0);

  FluidStatePhaseEnum phase_state;
  const unsigned int np = _fs->numPhases();
  const unsigned int nc = _fs->numComponents();
  std::vector<FluidStateProperties> fsp(np, FluidStateProperties(nc));

  // Liquid region
  ADReal Z = 0.0001;
  Moose::derivInsert(Z.derivatives(), _Zidx, 1.0);

  _fs->massFractions(p, T, Xnacl, Z, phase_state, fsp);
  EXPECT_EQ(phase_state, FluidStatePhaseEnum::LIQUID);

  // Verfify mass fraction values
  ADReal Xco2 = fsp[0].mass_fraction[1];
  ADReal Yco2 = fsp[1].mass_fraction[1];
  ADReal Xh2o = fsp[0].mass_fraction[0];
  ADReal Yh2o = fsp[1].mass_fraction[0];
  ADReal Xnacl2 = fsp[0].mass_fraction[2];
  ABS_TEST(Xco2.value(), Z.value(), 1.0e-8);
  ABS_TEST(Yco2.value(), 0.0, 1.0e-8);
  ABS_TEST(Xh2o.value(), 1.0 - Z.value(), 1.0e-8);
  ABS_TEST(Yh2o.value(), 0.0, 1.0e-8);
  ABS_TEST(Xnacl2.value(), Xnacl.value(), 1.0e-8);

  // Verify derivatives
  ABS_TEST(Xco2.derivatives()[_pidx], 0.0, 1.0e-8);
  ABS_TEST(Xco2.derivatives()[_Tidx], 0.0, 1.0e-8);
  ABS_TEST(Xco2.derivatives()[_Xidx], 0.0, 1.0e-8);
  ABS_TEST(Xco2.derivatives()[_Zidx], 1.0, 1.0e-8);
  ABS_TEST(Yco2.derivatives()[_pidx], 0.0, 1.0e-8);
  ABS_TEST(Yco2.derivatives()[_Tidx], 0.0, 1.0e-8);
  ABS_TEST(Yco2.derivatives()[_Xidx], 0.0, 1.0e-8);
  ABS_TEST(Yco2.derivatives()[_Zidx], 0.0, 1.0e-8);
  ABS_TEST(Xnacl.derivatives()[_Xidx], 1.0, 1.0e-8);

  // Gas region
  Z = 0.995;
  Moose::derivInsert(Z.derivatives(), _Zidx, 1.0);

  _fs->massFractions(p, T, Xnacl, Z, phase_state, fsp);
  EXPECT_EQ(phase_state, FluidStatePhaseEnum::GAS);

  // Verfify mass fraction values
  Xco2 = fsp[0].mass_fraction[1];
  Yco2 = fsp[1].mass_fraction[1];
  Xh2o = fsp[0].mass_fraction[0];
  Yh2o = fsp[1].mass_fraction[0];
  ADReal Ynacl = fsp[1].mass_fraction[2];
  ABS_TEST(Xco2.value(), 0.0, 1.0e-8);
  ABS_TEST(Yco2.value(), Z.value(), 1.0e-8);
  ABS_TEST(Xh2o.value(), 0.0, 1.0e-8);
  ABS_TEST(Yh2o.value(), 1.0 - Z.value(), 1.0e-8);
  ABS_TEST(Ynacl.value(), 0.0, 1.0e-8);

  // Verify derivatives
  ABS_TEST(Xco2.derivatives()[_pidx], 0.0, 1.0e-8);
  ABS_TEST(Xco2.derivatives()[_Tidx], 0.0, 1.0e-8);
  ABS_TEST(Xco2.derivatives()[_Xidx], 0.0, 1.0e-8);
  ABS_TEST(Xco2.derivatives()[_Zidx], 0.0, 1.0e-8);
  ABS_TEST(Yco2.derivatives()[_pidx], 0.0, 1.0e-8);
  ABS_TEST(Yco2.derivatives()[_Tidx], 0.0, 1.0e-8);
  ABS_TEST(Yco2.derivatives()[_Xidx], 0.0, 1.0e-8);
  ABS_TEST(Yco2.derivatives()[_Zidx], 1.0, 1.0e-8);
  ABS_TEST(Ynacl.derivatives()[_Xidx], 0.0, 1.0e-8);

  // Two phase region. In this region, the mass fractions and derivatives can
  // be verified using the equilibrium mass fraction derivatives that have
  // been verified above
  Z = 0.45;
  Moose::derivInsert(Z.derivatives(), _Zidx, 1.0);

  _fs->massFractions(p, T, Xnacl, Z, phase_state, fsp);
  EXPECT_EQ(phase_state, FluidStatePhaseEnum::TWOPHASE);

  // Equilibrium mass fractions and derivatives
  ADReal Xco2_eq, Yh2o_eq;
  _fs->equilibriumMassFractions(p, T, Xnacl, Xco2_eq, Yh2o_eq);

  // Verfify mass fraction values - comparing ADReals verifies that derivatives
  // are also identical
  Xco2 = fsp[0].mass_fraction[1];
  Yco2 = fsp[1].mass_fraction[1];
  Xh2o = fsp[0].mass_fraction[0];
  Yh2o = fsp[1].mass_fraction[0];
  ABS_TEST(Xco2, Xco2_eq, 1.0e-8);
  ABS_TEST(Yco2, 1.0 - Yh2o_eq, 1.0e-8);
  ABS_TEST(Xh2o, 1.0 - Xco2_eq, 1.0e-8);
  ABS_TEST(Yh2o, Yh2o_eq, 1.0e-8);

  // Use finite differences to verify derivative wrt Z is unaffected by Z
  const Real dZ = 1.0e-8;
  _fs->massFractions(p, T, Xnacl, Z + dZ, phase_state, fsp);
  ADReal Xco21 = fsp[0].mass_fraction[1];
  ADReal Yco21 = fsp[1].mass_fraction[1];
  _fs->massFractions(p, T, Xnacl, Z - dZ, phase_state, fsp);
  ADReal Xco22 = fsp[0].mass_fraction[1];
  ADReal Yco22 = fsp[1].mass_fraction[1];

  ABS_TEST(Xco2.derivatives()[_Zidx], (Xco21.value() - Xco22.value()) / (2.0 * dZ), 1.0e-8);
  ABS_TEST(Yco2.derivatives()[_Zidx], (Yco21.value() - Yco22.value()) / (2.0 * dZ), 1.0e-8);
}

TEST_F() [10/18]

TEST_F	(	PorousFlowBrineCO2Test	,
		gasProperties
	)

Definition at line 416 of file PorousFlowBrineCO2Test.C.

{
  ADReal p = 1.0e6;
  Moose::derivInsert(p.derivatives(), _pidx, 1.0);

  ADReal T = 350.0;
  Moose::derivInsert(T.derivatives(), _Tidx, 1.0);

  ADReal Xnacl = 0.1;
  Moose::derivInsert(Xnacl.derivatives(), _Xidx, 1.0);

  FluidStatePhaseEnum phase_state;
  const unsigned int np = _fs->numPhases();
  const unsigned int nc = _fs->numComponents();
  std::vector<FluidStateProperties> fsp(np, FluidStateProperties(nc));

  // Gas region
  ADReal Z = 0.995;
  Moose::derivInsert(Z.derivatives(), _Zidx, 1.0);

  _fs->massFractions(p, T, Xnacl, Z, phase_state, fsp);
  EXPECT_EQ(phase_state, FluidStatePhaseEnum::GAS);

  // Verify fluid density, viscosity and enthalpy
  _fs->gasProperties(p, T, fsp);
  ADReal gas_density = fsp[1].density;
  ADReal gas_viscosity = fsp[1].viscosity;
  ADReal gas_enthalpy = fsp[1].enthalpy;

  Real density = _co2_fp->rho_from_p_T(p.value(), T.value());
  Real viscosity = _co2_fp->mu_from_p_T(p.value(), T.value());
  Real enthalpy = _co2_fp->h_from_p_T(p.value(), T.value());

  ABS_TEST(gas_density.value(), density, 1.0e-8);
  ABS_TEST(gas_viscosity.value(), viscosity, 1.0e-8);
  ABS_TEST(gas_enthalpy.value(), enthalpy, 1.0e-8);

  // Verify derivatives
  const Real dp = 1.0e-1;
  _fs->gasProperties(p + dp, T, fsp);
  Real rho1 = fsp[1].density.value();
  Real mu1 = fsp[1].viscosity.value();
  Real h1 = fsp[1].enthalpy.value();

  _fs->gasProperties(p - dp, T, fsp);
  Real rho2 = fsp[1].density.value();
  Real mu2 = fsp[1].viscosity.value();
  Real h2 = fsp[1].enthalpy.value();

  REL_TEST(gas_density.derivatives()[_pidx], (rho1 - rho2) / (2.0 * dp), 1.0e-6);
  REL_TEST(gas_viscosity.derivatives()[_pidx], (mu1 - mu2) / (2.0 * dp), 1.0e-6);
  REL_TEST(gas_enthalpy.derivatives()[_pidx], (h1 - h2) / (2.0 * dp), 1.0e-6);

  const Real dT = 1.0e-3;
  _fs->gasProperties(p, T + dT, fsp);
  rho1 = fsp[1].density.value();
  mu1 = fsp[1].viscosity.value();
  h1 = fsp[1].enthalpy.value();

  _fs->gasProperties(p, T - dT, fsp);
  rho2 = fsp[1].density.value();
  mu2 = fsp[1].viscosity.value();
  h2 = fsp[1].enthalpy.value();

  REL_TEST(gas_density.derivatives()[_Tidx], (rho1 - rho2) / (2.0 * dT), 1.0e-6);
  REL_TEST(gas_viscosity.derivatives()[_Tidx], (mu1 - mu2) / (2.0 * dT), 1.0e-6);
  REL_TEST(gas_enthalpy.derivatives()[_Tidx], (h1 - h2) / (2.0 * dT), 1.0e-6);

  // Note: mass fraction changes with Z
  const Real dZ = 1.0e-8;
  _fs->massFractions(p, T, Xnacl, Z + dZ, phase_state, fsp);
  _fs->gasProperties(p, T, fsp);
  rho1 = fsp[1].density.value();
  mu1 = fsp[1].viscosity.value();
  h1 = fsp[1].enthalpy.value();

  _fs->massFractions(p, T, Xnacl, Z - dZ, phase_state, fsp);
  _fs->gasProperties(p, T, fsp);
  rho2 = fsp[1].density.value();
  mu2 = fsp[1].viscosity.value();
  h2 = fsp[1].enthalpy.value();

  ABS_TEST(gas_density.derivatives()[_Zidx], (rho1 - rho2) / (2.0 * dZ), 1.0e-8);
  ABS_TEST(gas_viscosity.derivatives()[_Zidx], (mu1 - mu2) / (2.0 * dZ), 1.0e-8);
  ABS_TEST(gas_enthalpy.derivatives()[_Zidx], (h1 - h2) / (2.0 * dZ), 1.0e-6);
}

TEST_F() [11/18]

TEST_F	(	PorousFlowBrineCO2Test	,
		liquidProperties
	)

Definition at line 506 of file PorousFlowBrineCO2Test.C.

{
  ADReal p = 1.0e6;
  Moose::derivInsert(p.derivatives(), _pidx, 1.0);

  ADReal T = 350.0;
  Moose::derivInsert(T.derivatives(), _Tidx, 1.0);

  ADReal Xnacl = 0.1;
  Moose::derivInsert(Xnacl.derivatives(), _Xidx, 1.0);

  FluidStatePhaseEnum phase_state;
  const unsigned int np = _fs->numPhases();
  const unsigned int nc = _fs->numComponents();
  std::vector<FluidStateProperties> fsp(np, FluidStateProperties(nc));

  // Liquid region
  ADReal Z = 0.0001;
  Moose::derivInsert(Z.derivatives(), _Zidx, 1.0);

  _fs->massFractions(p, T, Xnacl, Z, phase_state, fsp);
  EXPECT_EQ(phase_state, FluidStatePhaseEnum::LIQUID);

  // Verify fluid density and viscosity
  _fs->liquidProperties(p, T, Xnacl, fsp);

  ADReal liquid_density = fsp[0].density;
  ADReal liquid_viscosity = fsp[0].viscosity;
  ADReal liquid_enthalpy = fsp[0].enthalpy;

  Real co2_partial_density = _fs->partialDensityCO2(T).value();
  Real brine_density = _brine_fp->rho_from_p_T_X(p.value(), T.value(), Xnacl.value());

  Real density = 1.0 / (Z.value() / co2_partial_density + (1.0 - Z.value()) / brine_density);

  Real viscosity = _brine_fp->mu_from_p_T_X(p.value(), T.value(), Xnacl.value());

  Real brine_enthalpy = _brine_fp->h_from_p_T_X(p.value(), T.value(), Xnacl.value());
  Real hdis = _fs->enthalpyOfDissolution(T).value();
  Real co2_enthalpy = _co2_fp->h_from_p_T(p.value(), T.value());
  Real enthalpy = (1.0 - Z.value()) * brine_enthalpy + Z.value() * (co2_enthalpy + hdis);

  ABS_TEST(liquid_density.value(), density, 1.0e-12);
  ABS_TEST(liquid_viscosity.value(), viscosity, 1.0e-12);
  // Tolerance loosened to pass on Apple M4
  ABS_TEST(liquid_enthalpy.value(), enthalpy, 3.e-10);

  // Verify derivatives
  // Derivatives wrt pressure
  const Real dp = 1.0;
  _fs->liquidProperties(p + dp, T, Xnacl, fsp);
  Real rho1 = fsp[0].density.value();
  Real mu1 = fsp[0].viscosity.value();
  Real h1 = fsp[0].enthalpy.value();

  _fs->liquidProperties(p - dp, T, Xnacl, fsp);
  Real rho2 = fsp[0].density.value();
  Real mu2 = fsp[0].viscosity.value();
  Real h2 = fsp[0].enthalpy.value();

  REL_TEST(liquid_density.derivatives()[_pidx], (rho1 - rho2) / (2.0 * dp), 1.0e-6);
  REL_TEST(liquid_viscosity.derivatives()[_pidx], (mu1 - mu2) / (2.0 * dp), 2.0e-6);
  REL_TEST(liquid_enthalpy.derivatives()[_pidx], (h1 - h2) / (2.0 * dp), 1.0e-6);

  // Derivatives wrt temperature
  const Real dT = 1.0e-4;
  _fs->liquidProperties(p, T + dT, Xnacl, fsp);
  rho1 = fsp[0].density.value();
  mu1 = fsp[0].viscosity.value();
  h1 = fsp[0].enthalpy.value();

  _fs->liquidProperties(p, T - dT, Xnacl, fsp);
  rho2 = fsp[0].density.value();
  mu2 = fsp[0].viscosity.value();
  h2 = fsp[0].enthalpy.value();

  REL_TEST(liquid_density.derivatives()[_Tidx], (rho1 - rho2) / (2.0 * dT), 1.0e-6);
  REL_TEST(liquid_viscosity.derivatives()[_Tidx], (mu1 - mu2) / (2.0 * dT), 1.0e-6);
  REL_TEST(liquid_enthalpy.derivatives()[_Tidx], (h1 - h2) / (2.0 * dT), 1.0e-6);

  // Derivatives wrt Xnacl
  const Real dx = 1.0e-8;
  _fs->liquidProperties(p, T, Xnacl + dx, fsp);
  rho1 = fsp[0].density.value();
  mu1 = fsp[0].viscosity.value();
  h1 = fsp[0].enthalpy.value();

  _fs->liquidProperties(p, T, Xnacl - dx, fsp);
  rho2 = fsp[0].density.value();
  mu2 = fsp[0].viscosity.value();
  h2 = fsp[0].enthalpy.value();

  REL_TEST(liquid_density.derivatives()[_Xidx], (rho1 - rho2) / (2.0 * dx), 1.0e-6);
  REL_TEST(liquid_viscosity.derivatives()[_Xidx], (mu1 - mu2) / (2.0 * dx), 1.0e-6);
  REL_TEST(liquid_enthalpy.derivatives()[_Xidx], (h1 - h2) / (2.0 * dx), 1.0e-6);

  // Derivatives wrt Z
  const Real dZ = 1.0e-8;
  _fs->massFractions(p, T, Xnacl, Z + dZ, phase_state, fsp);
  _fs->liquidProperties(p, T, Xnacl, fsp);
  rho1 = fsp[0].density.value();
  mu1 = fsp[0].viscosity.value();
  h1 = fsp[0].enthalpy.value();

  _fs->massFractions(p, T, Xnacl, Z - dZ, phase_state, fsp);
  _fs->liquidProperties(p, T, Xnacl, fsp);
  rho2 = fsp[0].density.value();
  mu2 = fsp[0].viscosity.value();
  h2 = fsp[0].enthalpy.value();

  REL_TEST(liquid_density.derivatives()[_Zidx], (rho1 - rho2) / (2.0 * dZ), 1.0e-6);
  ABS_TEST(liquid_viscosity.derivatives()[_Zidx], (mu1 - mu2) / (2.0 * dZ), 1.0e-6);
  REL_TEST(liquid_enthalpy.derivatives()[_Zidx], (h1 - h2) / (2.0 * dZ), 1.0e-6);

  // Two-phase region
  Z = 0.045;
  Moose::derivInsert(Z.derivatives(), _Zidx, 1.0);

  _fs->massFractions(p, T, Xnacl, Z, phase_state, fsp);
  EXPECT_EQ(phase_state, FluidStatePhaseEnum::TWOPHASE);

  // Verify fluid density and viscosity derivatives
  _fs->liquidProperties(p, T, Xnacl, fsp);

  liquid_density = fsp[0].density;
  liquid_viscosity = fsp[0].viscosity;
  liquid_enthalpy = fsp[0].enthalpy;

  // Derivatives wrt pressure
  _fs->massFractions(p + dp, T, Xnacl, Z, phase_state, fsp);
  _fs->liquidProperties(p + dp, T, Xnacl, fsp);
  rho1 = fsp[0].density.value();
  mu1 = fsp[0].viscosity.value();
  h1 = fsp[0].enthalpy.value();

  _fs->massFractions(p - dp, T, Xnacl, Z, phase_state, fsp);
  _fs->liquidProperties(p - dp, T, Xnacl, fsp);
  rho2 = fsp[0].density.value();
  mu2 = fsp[0].viscosity.value();
  h2 = fsp[0].enthalpy.value();

  REL_TEST(liquid_density.derivatives()[_pidx], (rho1 - rho2) / (2.0 * dp), 1.0e-6);
  REL_TEST(liquid_viscosity.derivatives()[_pidx], (mu1 - mu2) / (2.0 * dp), 2.0e-6);
  REL_TEST(liquid_enthalpy.derivatives()[_pidx], (h1 - h2) / (2.0 * dp), 1.0e-6);

  // Derivatives wrt temperature
  _fs->massFractions(p, T + dT, Xnacl, Z, phase_state, fsp);
  _fs->liquidProperties(p, T + dT, Xnacl, fsp);
  rho1 = fsp[0].density.value();
  mu1 = fsp[0].viscosity.value();
  h1 = fsp[0].enthalpy.value();

  _fs->massFractions(p, T - dT, Xnacl, Z, phase_state, fsp);
  _fs->liquidProperties(p, T - dT, Xnacl, fsp);
  rho2 = fsp[0].density.value();
  mu2 = fsp[0].viscosity.value();
  h2 = fsp[0].enthalpy.value();

  REL_TEST(liquid_density.derivatives()[_Tidx], (rho1 - rho2) / (2.0 * dT), 1.0e-6);
  REL_TEST(liquid_viscosity.derivatives()[_Tidx], (mu1 - mu2) / (2.0 * dT), 1.0e-6);
  REL_TEST(liquid_enthalpy.derivatives()[_Tidx], (h1 - h2) / (2.0 * dT), 1.0e-6);

  // Derivatives wrt Xnacl
  _fs->massFractions(p, T, Xnacl + dx, Z, phase_state, fsp);
  _fs->liquidProperties(p, T, Xnacl + dx, fsp);
  rho1 = fsp[0].density.value();
  mu1 = fsp[0].viscosity.value();
  h1 = fsp[0].enthalpy.value();

  _fs->massFractions(p, T, Xnacl - dx, Z, phase_state, fsp);
  _fs->liquidProperties(p, T, Xnacl - dx, fsp);
  rho2 = fsp[0].density.value();
  mu2 = fsp[0].viscosity.value();
  h2 = fsp[0].enthalpy.value();

  REL_TEST(liquid_density.derivatives()[_Xidx], (rho1 - rho2) / (2.0 * dx), 1.0e-6);
  REL_TEST(liquid_viscosity.derivatives()[_Xidx], (mu1 - mu2) / (2.0 * dx), 1.0e-6);
  REL_TEST(liquid_enthalpy.derivatives()[_Xidx], (h1 - h2) / (2.0 * dx), 1.0e-6);

  // Derivatives wrt Z
  _fs->massFractions(p, T, Xnacl, Z + dZ, phase_state, fsp);
  _fs->liquidProperties(p, T, Xnacl, fsp);
  rho1 = fsp[0].density.value();
  mu1 = fsp[0].viscosity.value();
  h1 = fsp[0].enthalpy.value();

  _fs->massFractions(p, T, Xnacl, Z - dZ, phase_state, fsp);
  _fs->liquidProperties(p, T, Xnacl, fsp);
  rho2 = fsp[0].density.value();
  mu2 = fsp[0].viscosity.value();
  h2 = fsp[0].enthalpy.value();

  ABS_TEST(liquid_density.derivatives()[_Zidx], (rho1 - rho2) / (2.0 * dZ), 1.0e-6);
  ABS_TEST(liquid_viscosity.derivatives()[_Zidx], (mu1 - mu2) / (2.0 * dZ), 1.0e-6);
  ABS_TEST(liquid_enthalpy.derivatives()[_Zidx], (h1 - h2) / (2.0 * dZ), 1.0e-6);
}

TEST_F() [12/18]

TEST_F	(	PorousFlowBrineCO2Test	,
		saturation
	)

Definition at line 706 of file PorousFlowBrineCO2Test.C.

{
  ADReal p = 1.0e6;
  Moose::derivInsert(p.derivatives(), _pidx, 1.0);

  ADReal T = 350.0;
  Moose::derivInsert(T.derivatives(), _Tidx, 1.0);

  ADReal Xnacl = 0.1;
  Moose::derivInsert(Xnacl.derivatives(), _Xidx, 1.0);

  FluidStatePhaseEnum phase_state;
  const unsigned int np = _fs->numPhases();
  const unsigned int nc = _fs->numComponents();
  std::vector<FluidStateProperties> fsp(np, FluidStateProperties(nc));

  // In the two-phase region, the mass fractions are the equilibrium values, so
  // a temporary value of Z can be used (as long as it corresponds to the two-phase
  // region)
  ADReal Z = 0.45;
  Moose::derivInsert(Z.derivatives(), _Zidx, 1.0);

  _fs->massFractions(p, T, Xnacl, Z, phase_state, fsp);
  EXPECT_EQ(phase_state, FluidStatePhaseEnum::TWOPHASE);

  // Calculate Z that gives a saturation of 0.25
  ADReal gas_saturation = 0.25;
  ADReal liquid_pressure = p - _pc->capillaryPressure(1.0 - gas_saturation);

  // Calculate gas density and liquid density
  _fs->gasProperties(p, T, fsp);
  _fs->liquidProperties(liquid_pressure, T, Xnacl, fsp);

  // The mass fraction that corresponds to a gas_saturation = 0.25
  ADReal Zc = (gas_saturation * fsp[1].density * fsp[1].mass_fraction[1] +
               (1.0 - gas_saturation) * fsp[0].density * fsp[0].mass_fraction[1]) /
              (gas_saturation * fsp[1].density + (1.0 - gas_saturation) * fsp[0].density);

  // Calculate the gas saturation and derivatives
  ADReal saturation = _fs->saturation(p, T, Xnacl, Zc, fsp);
  ABS_TEST(saturation.value(), gas_saturation.value(), 1.0e-6);

  // Test the derivatives of gas saturation
  gas_saturation = _fs->saturation(p, T, Xnacl, Z, fsp);

  const Real dp = 1.0e-1;

  // Derivative wrt pressure
  _fs->massFractions(p + dp, T, Xnacl, Z, phase_state, fsp);
  Real gsat1 = _fs->saturation(p + dp, T, Xnacl, Z, fsp).value();

  _fs->massFractions(p - dp, T, Xnacl, Z, phase_state, fsp);
  Real gsat2 = _fs->saturation(p - dp, T, Xnacl, Z, fsp).value();

  REL_TEST(gas_saturation.derivatives()[_pidx], (gsat1 - gsat2) / (2.0 * dp), 1.0e-6);

  // Derivative wrt temperature
  const Real dT = 1.0e-4;
  _fs->massFractions(p, T + dT, Xnacl, Z, phase_state, fsp);
  gsat1 = _fs->saturation(p, T + dT, Xnacl, Z, fsp).value();

  _fs->massFractions(p, T - dT, Xnacl, Z, phase_state, fsp);
  gsat2 = _fs->saturation(p, T - dT, Xnacl, Z, fsp).value();

  REL_TEST(gas_saturation.derivatives()[_Tidx], (gsat1 - gsat2) / (2.0 * dT), 1.0e-6);

  // Derivative wrt Xnacl
  const Real dx = 1.0e-8;
  _fs->massFractions(p, T, Xnacl + dx, Z, phase_state, fsp);
  gsat1 = _fs->saturation(p, T, Xnacl + dx, Z, fsp).value();

  _fs->massFractions(p, T, Xnacl - dx, Z, phase_state, fsp);
  gsat2 = _fs->saturation(p, T, Xnacl - dx, Z, fsp).value();

  REL_TEST(gas_saturation.derivatives()[_Xidx], (gsat1 - gsat2) / (2.0 * dx), 1.0e-6);

  // Derivative wrt Z
  const Real dZ = 1.0e-8;

  _fs->massFractions(p, T, Xnacl, Z, phase_state, fsp);
  gsat1 = _fs->saturation(p, T, Xnacl, Z + dZ, fsp).value();

  gsat2 = _fs->saturation(p, T, Xnacl, Z - dZ, fsp).value();

  REL_TEST(gas_saturation.derivatives()[_Zidx], (gsat1 - gsat2) / (2.0 * dZ), 1.0e-6);
}

TEST_F() [13/18]

TEST_F	(	PorousFlowBrineCO2Test	,
		totalMassFraction
	)

Definition at line 796 of file PorousFlowBrineCO2Test.C.

{
  const Real p = 1.0e6;
  const Real T = 350.0;
  const Real Xnacl = 0.1;
  const Real s = 0.2;
  const unsigned qp = 0;

  Real Z = _fs->totalMassFraction(p, T, Xnacl, s, qp);

  // Test that the saturation calculated in this fluid state using Z is equal to s
  FluidStatePhaseEnum phase_state;
  const unsigned int np = _fs->numPhases();
  const unsigned int nc = _fs->numComponents();
  std::vector<FluidStateProperties> fsp(np, FluidStateProperties(nc));

  ADReal pressure = p;
  ADReal temperature = T;
  ADReal z = Z;
  ADReal x = Xnacl;
  _fs->massFractions(pressure, temperature, x, z, phase_state, fsp);
  EXPECT_EQ(phase_state, FluidStatePhaseEnum::TWOPHASE);

  ADReal gas_saturation = _fs->saturation(pressure, temperature, x, z, fsp);
  ABS_TEST(gas_saturation, s, 1.0e-6);
}

TEST_F() [14/18]

TEST_F	(	PorousFlowBrineCO2Test	,
		henryConstant
	)

Definition at line 829 of file PorousFlowBrineCO2Test.C.

{
  ADReal T = 373.15;
  Moose::derivInsert(T.derivatives(), _Tidx, 1.0);

  ADReal Xnacl = 0.1;
  Moose::derivInsert(Xnacl.derivatives(), _Xidx, 1.0);

  ADReal Kh = _fs->henryConstant(T, Xnacl);
  REL_TEST(Kh.value(), 7.46559e+08, 1.0e-3);

  T = 473.15;
  Moose::derivInsert(T.derivatives(), _Tidx, 1.0);

  Xnacl = 0.2;
  Moose::derivInsert(Xnacl.derivatives(), _Xidx, 1.0);

  Kh = _fs->henryConstant(T, Xnacl);
  REL_TEST(Kh.value(), 1.66069e+09, 1.0e-3);

  // Test the derivative wrt temperature
  const Real dT = 1.0e-4;
  ADReal Kh1 = _fs->henryConstant(T + dT, Xnacl);
  ADReal Kh2 = _fs->henryConstant(T - dT, Xnacl);

  REL_TEST(Kh.derivatives()[_Tidx], (Kh1.value() - Kh2.value()) / (2.0 * dT), 1.0e-6);

  // Test the derivative wrt Xnacl
  const Real dx = 1.0e-8;
  Kh1 = _fs->henryConstant(T, Xnacl + dx);
  Kh2 = _fs->henryConstant(T, Xnacl - dx);

  REL_TEST(Kh.derivatives()[_Xidx], (Kh1.value() - Kh2.value()) / (2.0 * dx), 1.0e-6);
}

TEST_F() [15/18]

TEST_F	(	PorousFlowBrineCO2Test	,
		enthalpyOfDissolutionGas
	)

Definition at line 869 of file PorousFlowBrineCO2Test.C.

{
  // T = 50C
  ADReal T = 323.15;
  Moose::derivInsert(T.derivatives(), _Tidx, 1.0);

  ADReal Xnacl = 0.1;
  Moose::derivInsert(Xnacl.derivatives(), _Xidx, 1.0);

  // Enthalpy of dissolution of CO2 in brine
  ADReal hdis = _fs->enthalpyOfDissolutionGas(T, Xnacl);
  REL_TEST(hdis, -3.20130e5, 1.0e-3);

  // T = 350C
  T = 623.15;
  Moose::derivInsert(T.derivatives(), _Tidx, 1.0);

  // Enthalpy of dissolution of CO2 in brine
  hdis = _fs->enthalpyOfDissolutionGas(T, Xnacl);
  REL_TEST(hdis, 9.83813e+05, 1.0e-3);

  // Test the derivative wrt temperature
  const Real dT = 1.0e-4;
  Real hdis1 = _fs->enthalpyOfDissolutionGas(T + dT, Xnacl).value();
  Real hdis2 = _fs->enthalpyOfDissolutionGas(T - dT, Xnacl).value();

  REL_TEST(hdis.derivatives()[_Tidx], (hdis1 - hdis2) / (2.0 * dT), 1.0e-5);

  // Test the derivative wrt salt mass fraction
  const Real dx = 1.0e-8;
  hdis1 = _fs->enthalpyOfDissolutionGas(T, Xnacl + dx).value();
  hdis2 = _fs->enthalpyOfDissolutionGas(T, Xnacl - dx).value();

  REL_TEST(hdis.derivatives()[_Xidx], (hdis1 - hdis2) / (2.0 * dx), 1.0e-5);
}

TEST_F() [16/18]

TEST_F	(	PorousFlowBrineCO2Test	,
		enthalpyOfDissolution
	)

Definition at line 911 of file PorousFlowBrineCO2Test.C.

{
  // T = 50C
  ADReal T = 323.15;
  Moose::derivInsert(T.derivatives(), _Tidx, 1.0);

  // Enthalpy of dissolution of CO2 in water
  ADReal hdis = _fs->enthalpyOfDissolution(T);
  REL_TEST(hdis.value(), -3.38185e5, 1.0e-3);

  // T = 350C
  T = 623.15;
  Moose::derivInsert(T.derivatives(), _Tidx, 1.0);

  // Enthalpy of dissolution of CO2 in water
  hdis = _fs->enthalpyOfDissolution(T);
  REL_TEST(hdis.value(), 5.78787e5, 1.0e-3);

  // Test the derivative wrt temperature
  const Real dT = 1.0e-4;
  Real hdis1 = _fs->enthalpyOfDissolution(T + dT).value();
  Real hdis2 = _fs->enthalpyOfDissolution(T - dT).value();

  REL_TEST(hdis.derivatives()[_Tidx], (hdis1 - hdis2) / (2.0 * dT), 1.0e-5);
}

TEST_F() [17/18]

TEST_F	(	PorousFlowBrineCO2Test	,
		solveEquilibriumMoleFractionHighTemp
	)

Test iterative calculation of equilibrium mass fractions in the elevated temperature regime.

Definition at line 941 of file PorousFlowBrineCO2Test.C.

{
  Real p = 20.0e6;
  Real T = 473.15;
  Real Xnacl = 0.0;

  Real yh2o, xco2;
  Real co2_density = _co2_fp->rho_from_p_T(p, T);
  _fs->solveEquilibriumMoleFractionHighTemp(p, T, Xnacl, co2_density, xco2, yh2o);
  ABS_TEST(yh2o, 0.161429743509, 1.0e-10);
  ABS_TEST(xco2, 0.0236966821858, 1.0e-10);

  // Mass fraction equivalent to molality of 2
  p = 30.0e6;
  T = 423.15;
  Xnacl = 0.105;

  _fs->solveEquilibriumMoleFractionHighTemp(p, T, Xnacl, co2_density, xco2, yh2o);
  ABS_TEST(yh2o, 0.0468197608955, 1.0e-10);
  ABS_TEST(xco2, 0.0236644599437, 1.0e-10);

  // Mass fraction equivalent to molality of 4
  T = 523.15;
  Xnacl = 0.189;

  co2_density = _co2_fp->rho_from_p_T(p, T);
  _fs->solveEquilibriumMoleFractionHighTemp(p, T, Xnacl, co2_density, xco2, yh2o);
  ABS_TEST(yh2o, 0.253292728991, 1.0e-10);
  ABS_TEST(xco2, 0.0168344513321, 1.0e-10);
}

TEST_F() [18/18]

TEST_F	(	PorousFlowBrineCO2Test	,
		equilibriumMoleFractions
	)

Test calculation of equilibrium mole fractions over entire temperature range.

Definition at line 975 of file PorousFlowBrineCO2Test.C.

{
  // Test pure water (Xnacl = 0)
  // Low temperature regime
  ADReal p = 20.0e6;
  ADReal T = 323.15;
  ADReal Xnacl = 0.0;

  ADReal x, y;
  _fs->equilibriumMoleFractions(p, T, Xnacl, x, y);
  ABS_TEST(y.value(), 0.00696393845155, 1.0e-8);
  ABS_TEST(x.value(), 0.0236554537395, 1.0e-8);

  // Intermediate temperature regime
  T = 373.15;

  _fs->equilibriumMoleFractions(p, T, Xnacl, x, y);
  ABS_TEST(y.value(), 0.0194394631944, 1.0e-8);
  ABS_TEST(x.value(), 0.020195776139, 1.0e-8);

  // High temperature regime
  p = 30.0e6;
  T = 523.15;

  _fs->equilibriumMoleFractions(p, T, Xnacl, x, y);
  ABS_TEST(y.value(), 0.286117195565, 1.0e-8);
  ABS_TEST(x.value(), 0.0409622051253, 1.0e-8);

  // High temperature regime with low pressure
  p = 1.0e6;
  T = 523.15;

  _fs->equilibriumMoleFractions(p, T, Xnacl, x, y);
  ABS_TEST(y.value(), 1.0, 1.0e-8);
  ABS_TEST(x.value(), 0.0, 1.0e-8);

  // Test brine (Xnacl = 0.1)
  // Low temperature regime
  p = 20.0e6;
  T = 323.15;
  Xnacl = 0.1;

  _fs->equilibriumMoleFractions(p, T, Xnacl, x, y);
  ABS_TEST(y.value(), 0.00657335846826, 1.0e-8);
  ABS_TEST(x.value(), 0.0152863933134, 1.0e-8);

  // Intermediate temperature regime
  T = 373.15;

  _fs->equilibriumMoleFractions(p, T, Xnacl, x, y);
  ABS_TEST(y.value(), 0.01831360857, 1.0e-8);
  ABS_TEST(x.value(), 0.0132653916293, 1.0e-8);

  // High temperature regime
  p = 30.0e6;
  T = 523.15;

  _fs->equilibriumMoleFractions(p, T, Xnacl, x, y);
  ABS_TEST(y.value(), 0.270258924237, 1.0e-8);
  ABS_TEST(x.value(), 0.0246589135776, 1.0e-8);
}