PorousFlowBrooksCoreyTest.C File Reference

Go to the source code of this file.

Functions
	TEST (PorousFlowBrooksCoreyTest, sat)
	TEST (PorousFlowBrooksCoreyTest, dsat)
	TEST (PorousFlowBrooksCoreyTest, d2sat)
	TEST (PorousFlowBrooksCoreyTest, cap)
	TEST (PorousFlowBrooksCoreyTest, dcap)
	TEST (PorousFlowBrooksCoreyTest, d2cap)
	TEST (PorousFlowBrooksCoreyTest, relpermw)
	TEST (PorousFlowBrooksCoreyTest, drelpermw)
	TEST (PorousFlowBrooksCoreyTest, relpermnw)
	TEST (PorousFlowBrooksCoreyTest, drelpermnw)
	TEST (PorousFlowBrooksCoreyTest, adrelpermw)
	TEST (PorousFlowBrooksCoreyTest, adrelpermnw)

Variables
const Real	eps = 1.0e-8
const Real	tol = 1.0e-6

Function Documentation

TEST() [1/12]

TEST	(	PorousFlowBrooksCoreyTest	,
		adrelpermnw
	)

Definition at line 136 of file PorousFlowBrooksCoreyTest.C.

{
  ADReal sat = 0.3;
  Moose::derivInsert(sat.derivatives(), 0, 1.0);
 
  const auto adrelperm = PorousFlowBrooksCorey::relativePermeabilityNW(sat, 2.5);
 
  const auto relperm = PorousFlowBrooksCorey::relativePermeabilityNW(sat.value(), 2.5);
  const auto drelperm = PorousFlowBrooksCorey::dRelativePermeabilityNW(sat.value(), 2.5);
 
  EXPECT_NEAR(adrelperm.value(), relperm, tol);
  EXPECT_NEAR(adrelperm.derivatives()[0], drelperm, tol);
}

TEST() [2/12]

TEST	(	PorousFlowBrooksCoreyTest	,
		adrelpermw
	)

Definition at line 122 of file PorousFlowBrooksCoreyTest.C.

{
  ADReal sat = 0.3;
  Moose::derivInsert(sat.derivatives(), 0, 1.0);
 
  const auto adrelperm = PorousFlowBrooksCorey::relativePermeabilityW(sat, 2.5);
 
  const auto relperm = PorousFlowBrooksCorey::relativePermeabilityW(sat.value(), 2.5);
  const auto drelperm = PorousFlowBrooksCorey::dRelativePermeabilityW(sat.value(), 2.5);
 
  EXPECT_NEAR(adrelperm.value(), relperm, tol);
  EXPECT_NEAR(adrelperm.derivatives()[0], drelperm, tol);
}

TEST() [3/12]

TEST	(	PorousFlowBrooksCoreyTest	,
		cap
	)

Definition at line 53 of file PorousFlowBrooksCoreyTest.C.

{
  EXPECT_NEAR(1.2, PorousFlowBrooksCorey::capillaryPressure(1.3, 1.2, 2.0, 1.0e9), tol);
  const Real pc = 1.2 * std::pow(0.3, -1.0 / 2.0);
  EXPECT_NEAR(pc, PorousFlowBrooksCorey::capillaryPressure(0.3, 1.2, 2.0, 1.0e9), tol);
  EXPECT_NEAR(1000.0, PorousFlowBrooksCorey::capillaryPressure(0.0, 1.2, 2.0, 1000.0), tol);
  EXPECT_NEAR(1.2, PorousFlowBrooksCorey::capillaryPressure(1.0, 1.2, 2.0, 1000.0), tol);
}

TEST() [4/12]

TEST	(	PorousFlowBrooksCoreyTest	,
		d2cap
	)

Definition at line 73 of file PorousFlowBrooksCoreyTest.C.

{
  Real fd;
  EXPECT_NEAR(0.0, PorousFlowBrooksCorey::d2CapillaryPressure(1.3, 1.2, 2.0, 1.0e9), tol);
  fd = (PorousFlowBrooksCorey::dCapillaryPressure(0.3 + eps, 1.2, 2.0, 1.0e9) -
        PorousFlowBrooksCorey::dCapillaryPressure(0.3, 1.2, 2.0, 1.0e9)) /
       eps;
  EXPECT_NEAR(fd, PorousFlowBrooksCorey::d2CapillaryPressure(0.3, 1.2, 2.0, 1.0e9), tol);
  EXPECT_NEAR(0.0, PorousFlowBrooksCorey::d2CapillaryPressure(1.01, 1.2, 2.0, 1.0e9), tol);
}

TEST() [5/12]

TEST	(	PorousFlowBrooksCoreyTest	,
		d2sat
	)

Definition at line 40 of file PorousFlowBrooksCoreyTest.C.

{
  Real fd;
  EXPECT_NEAR(0.0, PorousFlowBrooksCorey::d2EffectiveSaturation(0.99, 1.0, 2.0), tol);
 
  fd = (PorousFlowBrooksCorey::dEffectiveSaturation(2.5 + eps, 1.2, 2.0) -
        PorousFlowBrooksCorey::dEffectiveSaturation(2.5, 1.2, 2.0)) /
       eps;
  EXPECT_NEAR(fd, PorousFlowBrooksCorey::d2EffectiveSaturation(2.5, 1.2, 2.0), tol);
 
  EXPECT_NEAR(0.0, PorousFlowBrooksCorey::d2EffectiveSaturation(1.0e9, 1.2, 2.0), tol);
}

TEST() [6/12]

TEST	(	PorousFlowBrooksCoreyTest	,
		dcap
	)

Definition at line 62 of file PorousFlowBrooksCoreyTest.C.

{
  Real fd;
  EXPECT_NEAR(0.0, PorousFlowBrooksCorey::dCapillaryPressure(1.3, 1.2, 2.0, 1.0e9), tol);
  fd = (PorousFlowBrooksCorey::capillaryPressure(0.3 + eps, 1.2, 2.0, 1.0e9) -
        PorousFlowBrooksCorey::capillaryPressure(0.3, 1.2, 2.0, 1.0e9)) /
       eps;
  EXPECT_NEAR(fd, PorousFlowBrooksCorey::dCapillaryPressure(0.3, 1.2, 2.0, 1.0e9), tol);
  EXPECT_NEAR(0.0, PorousFlowBrooksCorey::dCapillaryPressure(1.01, 1.2, 2.0, 1.0e9), tol);
}

TEST() [7/12]

TEST	(	PorousFlowBrooksCoreyTest	,
		drelpermnw
	)

Definition at line 111 of file PorousFlowBrooksCoreyTest.C.

{
  Real fd;
  EXPECT_NEAR(0.0, PorousFlowBrooksCorey::dRelativePermeabilityNW(1.1, 2.5), tol);
  EXPECT_NEAR(0.0, PorousFlowBrooksCorey::dRelativePermeabilityNW(-1.0, 2.5), tol);
  fd = (PorousFlowBrooksCorey::relativePermeabilityNW(0.3 + eps, 2.5) -
        PorousFlowBrooksCorey::relativePermeabilityNW(0.3, 2.5)) /
       eps;
  EXPECT_NEAR(fd, PorousFlowBrooksCorey::dRelativePermeabilityNW(0.3, 2.5), tol);
}

TEST() [8/12]

TEST	(	PorousFlowBrooksCoreyTest	,
		drelpermw
	)

Definition at line 92 of file PorousFlowBrooksCoreyTest.C.

{
  Real fd;
  EXPECT_NEAR(0.0, PorousFlowBrooksCorey::dRelativePermeabilityW(1.1, 2.5), tol);
  EXPECT_NEAR(0.0, PorousFlowBrooksCorey::dRelativePermeabilityW(-1.0, 2.5), tol);
  fd = (PorousFlowBrooksCorey::relativePermeabilityW(0.3 + eps, 2.5) -
        PorousFlowBrooksCorey::relativePermeabilityW(0.3, 2.5)) /
       eps;
  EXPECT_NEAR(fd, PorousFlowBrooksCorey::dRelativePermeabilityW(0.3, 2.5), tol);
}

TEST() [9/12]

TEST	(	PorousFlowBrooksCoreyTest	,
		dsat
	)

Definition at line 27 of file PorousFlowBrooksCoreyTest.C.

{
  Real fd;
  EXPECT_NEAR(0.0, PorousFlowBrooksCorey::dEffectiveSaturation(0.99, 1.0, 2.0), tol);
 
  fd = (PorousFlowBrooksCorey::effectiveSaturation(2.5 + eps, 1.2, 2.0) -
        PorousFlowBrooksCorey::effectiveSaturation(2.5, 1.2, 2.0)) /
       eps;
  EXPECT_NEAR(fd, PorousFlowBrooksCorey::dEffectiveSaturation(2.5, 1.2, 2.0), tol);
 
  EXPECT_NEAR(0.0, PorousFlowBrooksCorey::dEffectiveSaturation(1.0e9, 1.2, 2.0), tol);
}

TEST() [10/12]

TEST	(	PorousFlowBrooksCoreyTest	,
		relpermnw
	)

Definition at line 103 of file PorousFlowBrooksCoreyTest.C.

{
  EXPECT_NEAR(1.0, PorousFlowBrooksCorey::relativePermeabilityNW(1.1, 2.5), tol);
  EXPECT_NEAR(0.0, PorousFlowBrooksCorey::relativePermeabilityNW(-1.0, 2.5), tol);
  const Real relperm = 0.3 * 0.3 * (1.0 - std::pow(1.0 - 0.3, (2.0 + 2.5) / 2.5));
  EXPECT_NEAR(relperm, PorousFlowBrooksCorey::relativePermeabilityNW(0.3, 2.5), tol);
}

TEST() [11/12]

TEST	(	PorousFlowBrooksCoreyTest	,
		relpermw
	)

Definition at line 84 of file PorousFlowBrooksCoreyTest.C.

{
  EXPECT_NEAR(1.0, PorousFlowBrooksCorey::relativePermeabilityW(1.1, 2.5), tol);
  EXPECT_NEAR(0.0, PorousFlowBrooksCorey::relativePermeabilityW(-1.0, 2.5), tol);
  const Real relperm = std::pow(0.3, (2.0 + 3.0 * 2.5) / 2.5);
  EXPECT_NEAR(relperm, PorousFlowBrooksCorey::relativePermeabilityW(0.3, 2.5), tol);
}

TEST() [12/12]

TEST	(	PorousFlowBrooksCoreyTest	,
		sat
	)

Definition at line 17 of file PorousFlowBrooksCoreyTest.C.

{
  EXPECT_NEAR(1.0, PorousFlowBrooksCorey::effectiveSaturation(0.99, 1.0, 2.0), tol);
  EXPECT_NEAR(1.0, PorousFlowBrooksCorey::effectiveSaturation(1.0, 1.0, 2.0), tol);
 
  const Real seff = std::pow(2.5 / 1.2, -2.0);
  EXPECT_NEAR(seff, PorousFlowBrooksCorey::effectiveSaturation(2.5, 1.2, 2.0), tol);
  EXPECT_NEAR(0.0, PorousFlowBrooksCorey::effectiveSaturation(1.0e9, 1.2, 2.0), tol);
}

Variable Documentation

eps

const Real eps = 1.0e-8

Definition at line 14 of file PorousFlowBrooksCoreyTest.C.

Referenced by TEST(), TEST(), TEST(), TEST(), TEST(), and TEST().

tol

const Real tol = 1.0e-6

Definition at line 15 of file PorousFlowBrooksCoreyTest.C.

Referenced by TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), TEST(), and TEST().