PorousFlowBroadbridgeWhiteTest.C File Reference

Go to the source code of this file.

Classes
class	PorousFlowBroadbridgeWhiteTest

Functions
	TEST_F (PorousFlowBroadbridgeWhiteTest, sat)
	TEST_F (PorousFlowBroadbridgeWhiteTest, dsat)
	TEST_F (PorousFlowBroadbridgeWhiteTest, d2sat)
	TEST_F (PorousFlowBroadbridgeWhiteTest, relperm)
	TEST_F (PorousFlowBroadbridgeWhiteTest, drelperm)
	TEST_F (PorousFlowBroadbridgeWhiteTest, d2relperm)
	TEST_F (PorousFlowBroadbridgeWhiteTest, adrelperm)

Function Documentation

TEST_F() [1/7]

TEST_F	(	PorousFlowBroadbridgeWhiteTest	,
		adrelperm
	)

Definition at line 105 of file PorousFlowBroadbridgeWhiteTest.C.

{
  ADReal sat = 0.3;
  sat.derivatives() = {};
  Moose::derivInsert(sat.derivatives(), 0, 1.0);
 
  const auto adrelperm =
      PorousFlowBroadbridgeWhite::relativePermeability(sat, _c, _sn, _ss, _kn, _ks);
 
  const auto relperm =
      PorousFlowBroadbridgeWhite::relativePermeability(sat.value(), _c, _sn, _ss, _kn, _ks);
  const auto drelperm =
      PorousFlowBroadbridgeWhite::dRelativePermeability(sat.value(), _c, _sn, _ss, _kn, _ks);
 
  EXPECT_NEAR(adrelperm.value(), relperm, 1.0E-5);
  EXPECT_NEAR(adrelperm.derivatives()[0], drelperm, 1.0E-5);
}

TEST_F() [2/7]

TEST_F	(	PorousFlowBroadbridgeWhiteTest	,
		d2relperm
	)

Definition at line 88 of file PorousFlowBroadbridgeWhiteTest.C.

{
  EXPECT_NEAR(0.0,
              PorousFlowBroadbridgeWhite::d2RelativePermeability(0.01, _c, _sn, _ss, _kn, _ks),
              1.0E-5);
  const Real sat = 0.6;
  const Real fd =
      (PorousFlowBroadbridgeWhite::dRelativePermeability(sat + _ep, _c, _sn, _ss, _kn, _ks) -
       PorousFlowBroadbridgeWhite::dRelativePermeability(sat, _c, _sn, _ss, _kn, _ks)) /
      _ep;
  EXPECT_NEAR(
      fd, PorousFlowBroadbridgeWhite::d2RelativePermeability(sat, _c, _sn, _ss, _kn, _ks), 1.0E-5);
  EXPECT_NEAR(0.0,
              PorousFlowBroadbridgeWhite::d2RelativePermeability(0.99, _c, _sn, _ss, _kn, _ks),
              1.0E-5);
}

TEST_F() [3/7]

TEST_F	(	PorousFlowBroadbridgeWhiteTest	,
		d2sat
	)

Definition at line 47 of file PorousFlowBroadbridgeWhiteTest.C.

{
  EXPECT_NEAR(
      0.0, PorousFlowBroadbridgeWhite::d2EffectiveSaturation(50, 0.7, 0.5, 0.6, 1.0), 1.0E-5);
  const Real fd =
      (PorousFlowBroadbridgeWhite::dEffectiveSaturation(-1.0 + _ep, _c, _sn, _ss, _las) -
       PorousFlowBroadbridgeWhite::dEffectiveSaturation(-1.0, _c, _sn, _ss, _las)) /
      _ep;
  EXPECT_NEAR(
      fd, PorousFlowBroadbridgeWhite::d2EffectiveSaturation(-1.0, _c, _sn, _ss, _las), 1.0E-5);
}

TEST_F() [4/7]

TEST_F	(	PorousFlowBroadbridgeWhiteTest	,
		drelperm
	)

Definition at line 73 of file PorousFlowBroadbridgeWhiteTest.C.

{
  EXPECT_NEAR(
      0.0, PorousFlowBroadbridgeWhite::dRelativePermeability(0.01, _c, _sn, _ss, _kn, _ks), 1.0E-5);
  const Real sat = 0.4;
  const Real fd =
      (PorousFlowBroadbridgeWhite::relativePermeability(sat + _ep, _c, _sn, _ss, _kn, _ks) -
       PorousFlowBroadbridgeWhite::relativePermeability(sat, _c, _sn, _ss, _kn, _ks)) /
      _ep;
  EXPECT_NEAR(
      fd, PorousFlowBroadbridgeWhite::dRelativePermeability(sat, _c, _sn, _ss, _kn, _ks), 1.0E-5);
  EXPECT_NEAR(
      0.0, PorousFlowBroadbridgeWhite::dRelativePermeability(0.99, _c, _sn, _ss, _kn, _ks), 1.0E-5);
}

TEST_F() [5/7]

TEST_F	(	PorousFlowBroadbridgeWhiteTest	,
		dsat
	)

Definition at line 36 of file PorousFlowBroadbridgeWhiteTest.C.

{
  EXPECT_NEAR(
      0.0, PorousFlowBroadbridgeWhite::dEffectiveSaturation(50, 0.7, 0.5, 0.6, 1.0), 1.0E-5);
  const Real fd = (PorousFlowBroadbridgeWhite::effectiveSaturation(-1.0 + _ep, _c, _sn, _ss, _las) -
                   PorousFlowBroadbridgeWhite::effectiveSaturation(-1.0, _c, _sn, _ss, _las)) /
                  _ep;
  EXPECT_NEAR(
      fd, PorousFlowBroadbridgeWhite::dEffectiveSaturation(-1.0, _c, _sn, _ss, _las), 1.0E-5);
}

TEST_F() [6/7]

TEST_F	(	PorousFlowBroadbridgeWhiteTest	,
		relperm
	)

Definition at line 59 of file PorousFlowBroadbridgeWhiteTest.C.

{
  EXPECT_NEAR(
      _kn, PorousFlowBroadbridgeWhite::relativePermeability(0.01, _c, _sn, _ss, _kn, _ks), 1.0E-5);
  const Real sat = 0.5;
  const Real th = (sat - _sn) / (_ss - _sn);
  const Real expect = _kn + (_ks - _kn) * th * th * (_c - 1.0) / (_c - th);
  EXPECT_NEAR(expect,
              PorousFlowBroadbridgeWhite::relativePermeability(sat, _c, _sn, _ss, _kn, _ks),
              1.0E-5);
  EXPECT_NEAR(
      _ks, PorousFlowBroadbridgeWhite::relativePermeability(0.99, _c, _sn, _ss, _kn, _ks), 1.0E-5);
}

TEST_F() [7/7]

TEST_F	(	PorousFlowBroadbridgeWhiteTest	,
		sat
	)

Definition at line 26 of file PorousFlowBroadbridgeWhiteTest.C.

{
  EXPECT_NEAR(1.0, PorousFlowBroadbridgeWhite::effectiveSaturation(50, _c, _sn, _ss, _las), 1.0E-5);
  const Real eff = PorousFlowBroadbridgeWhite::effectiveSaturation(-1.0, _c, _sn, _ss, _las);
  const Real th = (eff - _sn) / (_ss - _sn);
  const Real t1 = (1.0 / _c) * std::log((_c - th) / (_c - 1.0) / th);
  const Real t2 = (th - 1.0) / th;
  EXPECT_NEAR(-1.0, _las * (t2 - t1), 1.0E-5);
}