PorousFlowFluidStateFlashTest.C File Reference

Go to the source code of this file.

Functions
	TEST_F (PorousFlowFluidStateFlashTest, rachfordRice)
	Verify construction of the Rachford-Rice equation. More...
	TEST_F (PorousFlowFluidStateFlashTest, vaporMassFraction)
	Verify iterative solution of the Rachford-Rice equation to calculate the vapor mass fraction. More...
	TEST_F (PorousFlowFluidStateFlashTest, rachfordRiceDeriv)
	Verify calculation of the derivative of the Rachford-Rice equation. More...

Function Documentation

TEST_F() [1/3]

TEST_F	(	PorousFlowFluidStateFlashTest	,
		rachfordRice
	)

Verify construction of the Rachford-Rice equation.

Definition at line 16 of file PorousFlowFluidStateFlashTest.C.

{
  // Two fluid components
  Real vmf = 0.2316623869599;
  std::vector<Real> zi = {0.6};
  std::vector<Real> Ki = {1.338, 0.576};
  ABS_TEST(_fp->rachfordRice(vmf, zi, Ki), 0.0, 1.0e-8);

  // Four fluid components
  vmf = 0.20329862165314910428;
  zi = {0.6, 0.01, 0.01};
  Ki = {1.338, 0.613, 0.222, 0.576};
  ABS_TEST(_fp->rachfordRice(vmf, zi, Ki), 0.0, 1.0e-8);
}

TEST_F() [2/3]

TEST_F	(	PorousFlowFluidStateFlashTest	,
		vaporMassFraction
	)

Verify iterative solution of the Rachford-Rice equation to calculate the vapor mass fraction.

Definition at line 35 of file PorousFlowFluidStateFlashTest.C.

{
  // Test calculation using example data
  std::vector<Real> zi = {0.6};
  std::vector<Real> Ki = {1.338, 0.576};
  ABS_TEST(_fp->vaporMassFraction(zi, Ki), 0.2316623869599, 1.0e-8);

  // Four fluid components
  zi = {0.6, 0.01, 0.01};
  Ki = {1.338, 0.613, 0.222, 0.576};
  ABS_TEST(_fp->vaporMassFraction(zi, Ki), 0.20329862165314910428, 1.0e-8);
}

TEST_F() [3/3]

TEST_F	(	PorousFlowFluidStateFlashTest	,
		rachfordRiceDeriv
	)

Verify calculation of the derivative of the Rachford-Rice equation.

Definition at line 51 of file PorousFlowFluidStateFlashTest.C.

{
  // Two fluid components
  Real vmf = 0.25;
  const Real dvmf = 1.0e-8;
  std::vector<Real> zi = {0.6};
  std::vector<Real> Ki = {1.338, 0.576};

  Real rr1 = _fp->rachfordRice(vmf - dvmf, zi, Ki);
  Real rr2 = _fp->rachfordRice(vmf + dvmf, zi, Ki);
  Real fd = (rr2 - rr1) / (2.0 * dvmf);

  ABS_TEST(_fp->rachfordRiceDeriv(vmf, zi, Ki), fd, 1.0e-8);

  // Four fluid components
  vmf = 0.6;
  zi = {0.6, 0.01, 0.01};
  Ki = {1.338, 0.613, 0.222, 0.576};

  rr1 = _fp->rachfordRice(vmf - dvmf, zi, Ki);
  rr2 = _fp->rachfordRice(vmf + dvmf, zi, Ki);
  fd = (rr2 - rr1) / (2.0 * dvmf);

  ABS_TEST(_fp->rachfordRiceDeriv(vmf, zi, Ki), fd, 1.0e-8);
}