Loading [MathJax]/extensions/tex2jax.js
https://mooseframework.inl.gov
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends
Functions | Variables
PorousFlowVanGenuchtenTest.C File Reference

Go to the source code of this file.

Functions

 TEST (PorousFlowVanGenuchtenTest, sat)
 
 TEST (PorousFlowVanGenuchtenTest, dsat)
 
 TEST (PorousFlowVanGenuchtenTest, d2sat)
 
 TEST (PorousFlowVanGenuchtenTest, cap)
 
 TEST (PorousFlowVanGenuchtenTest, dcap)
 
 TEST (PorousFlowVanGenuchtenTest, d2cap)
 
 TEST (PorousFlowVanGenuchtenTest, relperm)
 
 TEST (PorousFlowVanGenuchtenTest, drelperm)
 
 TEST (PorousFlowVanGenuchtenTest, d2relperm)
 
 TEST (PorousFlowVanGenuchtenTest, adrelperm)
 
 TEST (PorousFlowVanGenuchtenTest, relpermNW)
 
 TEST (PorousFlowVanGenuchtenTest, drelpermNW)
 
 TEST (PorousFlowVanGenuchtenTest, d2relpermNW)
 
 TEST (PorousFlowVanGenuchtenTest, adrelpermnw)
 
 TEST (PorousFlowVanGenuchtenTest, caphys)
 Test the hysteretic capillary pressure. More...
 
 TEST (PorousFlowVanGenuchtenTest, sathys)
 Test the hysteretic saturation. More...
 
 TEST (PorousFlowVanGenuchtenTest, dcaphys)
 Test the first derivative of the hysteretic capillary pressure. More...
 
 TEST (PorousFlowVanGenuchtenTest, d2caphys)
 Test the second derivative of the hysteretic capillary pressure. More...
 
 TEST (PorousFlowVanGenuchtenTest, dsathys)
 Test the derivative of the hysteretic saturation. More...
 
 TEST (PorousFlowVanGenuchtenTest, d2sathys)
 Test the second derivative of the hysteretic saturation. More...
 
 TEST (PorousFlowVanGenuchtenTest, relativePermeabilityHys)
 Test relativePermeabilityHys. More...
 
 TEST (PorousFlowVanGenuchtenTest, drelativePermeabilityHys)
 Test drelativePermeabilityHys. More...
 
 TEST (PorousFlowVanGenuchtenTest, relativePermeabilityNWHys)
 Test relativePermeabilityNWHys. More...
 
 TEST (PorousFlowVanGenuchtenTest, drelativePermeabilityNWHys)
 Test drelativePermeabilityNWHys. More...
 

Variables

const double eps = 1.0E-8
 
const PorousFlowVanGenuchten::LowCapillaryPressureExtension no_low_ext (PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, -1.0, 123.0, 0.0)
 
const Real low_ext_Pc = 1.5
 
const Real low_ext_S = PorousFlowVanGenuchten::saturationHys(low_ext_Pc, 0.1, 0.3, 10.0, 1.9)
 
const Real low_ext_dPc
 
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_exp (PorousFlowVanGenuchten::LowCapillaryPressureExtension::EXPONENTIAL, low_ext_S, low_ext_Pc, low_ext_dPc)
 
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_quad (PorousFlowVanGenuchten::LowCapillaryPressureExtension::QUADRATIC, low_ext_S, low_ext_Pc, low_ext_dPc)
 
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_none (PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, low_ext_S, low_ext_Pc, low_ext_dPc)
 
const PorousFlowVanGenuchten::HighCapillaryPressureExtension no_high_ext (PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, 2.0, 0.0, 0.0)
 
const Real high_ext_S = 0.9 * (1.0 - 0.3)
 
const Real high_ext_Pc
 
const Real high_ext_dPc
 
const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_power (PorousFlowVanGenuchten::HighCapillaryPressureExtension::POWER, high_ext_S, high_ext_Pc, high_ext_dPc)
 
const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_none (PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, 0.8, high_ext_Pc, high_ext_dPc)
 
const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_none2 (PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, high_ext_S, high_ext_Pc, high_ext_dPc)
 

Function Documentation

◆ TEST() [1/24]

TEST ( PorousFlowVanGenuchtenTest  ,
sat   
)

Definition at line 57 of file PorousFlowVanGenuchtenTest.C.

58 {
59  EXPECT_NEAR(1.0, PorousFlowVanGenuchten::effectiveSaturation(1.0E30, 0.7, 0.5), 1.0E-5);
60  EXPECT_NEAR(1.0, PorousFlowVanGenuchten::effectiveSaturation(1.0, 0.7, 0.5), 1.0E-5);
61  EXPECT_NEAR(1.0, PorousFlowVanGenuchten::effectiveSaturation(0.0, 0.7, 0.5), 1.0E-5);
62  EXPECT_NEAR(1.0, PorousFlowVanGenuchten::effectiveSaturation(1.0E-10, 0.7, 0.5), 1.0E-5);
63  EXPECT_NEAR(
64  0.486841442435055, PorousFlowVanGenuchten::effectiveSaturation(-2.0, 0.7, 0.6), 1.0E-5);
65  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::effectiveSaturation(-1.0E30, 0.7, 0.5), 1.0E-5);
66 }
PorousFlowVanGenuchten::effectiveSaturation
Real effectiveSaturation(Real p, Real alpha, Real m)
Effective saturation as a function of porepressure.
Definition: PorousFlowVanGenuchten.C:16

◆ TEST() [2/24]

TEST ( PorousFlowVanGenuchtenTest  ,
dsat   
)

Definition at line 68 of file PorousFlowVanGenuchtenTest.C.

69 {
70  Real fd;
71  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dEffectiveSaturation(1.0E30, 0.7, 0.5), 1.0E-5);
72  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dEffectiveSaturation(1.0, 0.7, 0.5), 1.0E-5);
73  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dEffectiveSaturation(0.0, 0.7, 0.5), 1.0E-5);
74  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dEffectiveSaturation(1.0E-10, 0.7, 0.5), 1.0E-5);
75  fd = (PorousFlowVanGenuchten::effectiveSaturation(-2.0 + eps, 0.7, 0.6) -
76  PorousFlowVanGenuchten::effectiveSaturation(-2.0, 0.7, 0.6)) /
77  eps;
78  EXPECT_NEAR(fd, PorousFlowVanGenuchten::dEffectiveSaturation(-2.0, 0.7, 0.6), 1.0E-5);
79  fd = (PorousFlowVanGenuchten::effectiveSaturation(-1.1 + eps, 0.9, 0.66) -
80  PorousFlowVanGenuchten::effectiveSaturation(-1.1, 0.9, 0.66)) /
81  eps;
82  EXPECT_NEAR(fd, PorousFlowVanGenuchten::dEffectiveSaturation(-1.1, 0.9, 0.66), 1.0E-5);
83  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dEffectiveSaturation(-1.0E30, 0.7, 0.5), 1.0E-5);
84 }
eps
const double eps
Definition: PorousFlowVanGenuchtenTest.C:14
PorousFlowVanGenuchten::effectiveSaturation
Real effectiveSaturation(Real p, Real alpha, Real m)
Effective saturation as a function of porepressure.
Definition: PorousFlowVanGenuchten.C:16
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
PorousFlowVanGenuchten::dEffectiveSaturation
Real dEffectiveSaturation(Real p, Real alpha, Real m)
Derivative of effective saturation wrt porepressure.
Definition: PorousFlowVanGenuchten.C:31

◆ TEST() [3/24]

TEST ( PorousFlowVanGenuchtenTest  ,
d2sat   
)

Definition at line 86 of file PorousFlowVanGenuchtenTest.C.

87 {
88  Real fd;
89  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2EffectiveSaturation(1.0E30, 0.7, 0.5), 1.0E-5);
90  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2EffectiveSaturation(1.0, 0.7, 0.5), 1.0E-5);
91  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2EffectiveSaturation(0.0, 0.7, 0.5), 1.0E-5);
92  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2EffectiveSaturation(1.0E-10, 0.7, 0.5), 1.0E-5);
93  fd = (PorousFlowVanGenuchten::dEffectiveSaturation(-2.0 + eps, 0.7, 0.6) -
94  PorousFlowVanGenuchten::dEffectiveSaturation(-2.0, 0.7, 0.6)) /
95  eps;
96  EXPECT_NEAR(fd, PorousFlowVanGenuchten::d2EffectiveSaturation(-2.0, 0.7, 0.6), 1.0E-5);
97  fd = (PorousFlowVanGenuchten::dEffectiveSaturation(-1.1 + eps, 2.3, 0.67) -
98  PorousFlowVanGenuchten::dEffectiveSaturation(-1.1, 2.3, 0.67)) /
99  eps;
100  EXPECT_NEAR(fd, PorousFlowVanGenuchten::d2EffectiveSaturation(-1.1, 2.3, 0.67), 1.0E-5);
101  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2EffectiveSaturation(-1.0E30, 0.7, 0.5), 1.0E-5);
102 }
eps
const double eps
Definition: PorousFlowVanGenuchtenTest.C:14
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
PorousFlowVanGenuchten::d2EffectiveSaturation
Real d2EffectiveSaturation(Real p, Real alpha, Real m)
Second derivative of effective saturation wrt porepressure.
Definition: PorousFlowVanGenuchten.C:46
PorousFlowVanGenuchten::dEffectiveSaturation
Real dEffectiveSaturation(Real p, Real alpha, Real m)
Derivative of effective saturation wrt porepressure.
Definition: PorousFlowVanGenuchten.C:31

◆ TEST() [4/24]

TEST ( PorousFlowVanGenuchtenTest  ,
cap   
)

Definition at line 104 of file PorousFlowVanGenuchtenTest.C.

105 {
106  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::capillaryPressure(1.1, 1.0, 0.55, 1.0E30), 1.0E-5);
107  EXPECT_NEAR(4.06172392297447,
108  PorousFlowVanGenuchten::capillaryPressure(0.3, 0.625, 0.55, 1.0E30),
109  1.0E-5);
110  EXPECT_NEAR(2.9, PorousFlowVanGenuchten::capillaryPressure(0.001, 0.625, 0.55, 2.9), 1.0E-5);
111  EXPECT_NEAR(1000.0, PorousFlowVanGenuchten::capillaryPressure(0.0, 0.625, 0.55, 1000.0), 1.0E-5);
112 }
PorousFlowVanGenuchten::capillaryPressure
Real capillaryPressure(Real seff, Real alpha, Real m, Real pc_max)
Capillary pressure as a function of effective saturation.
Definition: PorousFlowVanGenuchten.C:63

◆ TEST() [5/24]

TEST ( PorousFlowVanGenuchtenTest  ,
dcap   
)

Definition at line 114 of file PorousFlowVanGenuchtenTest.C.

115 {
116  Real fd;
117  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dCapillaryPressure(0.0, 1.0, 0.55, 1.0E30), 1.0E-5);
118  fd = (PorousFlowVanGenuchten::capillaryPressure(0.3 + eps, 0.625, 0.55, 1.0E30) -
119  PorousFlowVanGenuchten::capillaryPressure(0.3, 0.625, 0.55, 1.0E30)) /
120  eps;
121  EXPECT_NEAR(fd, PorousFlowVanGenuchten::dCapillaryPressure(0.3, 0.625, 0.55, 1.0E30), 1.0E-5);
122  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dCapillaryPressure(1.0, 1.0, 0.55, 1.0E30), 1.0E-5);
123 }
eps
const double eps
Definition: PorousFlowVanGenuchtenTest.C:14
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
PorousFlowVanGenuchten::dCapillaryPressure
Real dCapillaryPressure(Real seff, Real alpha, Real m, Real pc_max)
Derivative of capillary pressure wrt effective saturation.
Definition: PorousFlowVanGenuchten.C:77
PorousFlowVanGenuchten::capillaryPressure
Real capillaryPressure(Real seff, Real alpha, Real m, Real pc_max)
Capillary pressure as a function of effective saturation.
Definition: PorousFlowVanGenuchten.C:63

◆ TEST() [6/24]

TEST ( PorousFlowVanGenuchtenTest  ,
d2cap   
)

Definition at line 125 of file PorousFlowVanGenuchtenTest.C.

126 {
127  Real fd;
128  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2CapillaryPressure(0.0, 1.0, 0.55, 1.0E30), 1.0E-5);
129  fd = (PorousFlowVanGenuchten::dCapillaryPressure(0.3 + eps, 0.625, 0.55, 1.0E30) -
130  PorousFlowVanGenuchten::dCapillaryPressure(0.3, 0.625, 0.55, 1.0E30)) /
131  eps;
132  EXPECT_NEAR(fd, PorousFlowVanGenuchten::d2CapillaryPressure(0.3, 0.625, 0.55, 1.0E30), 1.0E-5);
133  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2CapillaryPressure(1.0, 0.625, 0.55, 1.0E30), 1.0E-5);
134 }
eps
const double eps
Definition: PorousFlowVanGenuchtenTest.C:14
PorousFlowVanGenuchten::d2CapillaryPressure
Real d2CapillaryPressure(Real seff, Real alpha, Real m, Real pc_max)
Second derivative of capillary pressure wrt effective saturation.
Definition: PorousFlowVanGenuchten.C:93
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
PorousFlowVanGenuchten::dCapillaryPressure
Real dCapillaryPressure(Real seff, Real alpha, Real m, Real pc_max)
Derivative of capillary pressure wrt effective saturation.
Definition: PorousFlowVanGenuchten.C:77

◆ TEST() [7/24]

TEST ( PorousFlowVanGenuchtenTest  ,
relperm   
)

Definition at line 136 of file PorousFlowVanGenuchtenTest.C.

137 {
138  EXPECT_NEAR(1.0, PorousFlowVanGenuchten::relativePermeability(1.0E30, 0.7), 1.0E-5);
139  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::relativePermeability(-1.0, 0.7), 1.0E-5);
140  EXPECT_NEAR(0.0091160727, PorousFlowVanGenuchten::relativePermeability(0.3, 0.7), 1.0E-5);
141 }
PorousFlowVanGenuchten::relativePermeability
T relativePermeability(const T &seff, Real m)
Relative permeability as a function of effective saturation.
Definition: PorousFlowVanGenuchten.h:101

◆ TEST() [8/24]

TEST ( PorousFlowVanGenuchtenTest  ,
drelperm   
)

Definition at line 143 of file PorousFlowVanGenuchtenTest.C.

144 {
145  Real fd;
146  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dRelativePermeability(1.0E30, 0.7), 1.0E-5);
147  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dRelativePermeability(-1.0, 0.7), 1.0E-5);
148  fd = (PorousFlowVanGenuchten::relativePermeability(0.3 + eps, 0.7) -
149  PorousFlowVanGenuchten::relativePermeability(0.3, 0.7)) /
150  eps;
151  EXPECT_NEAR(fd, PorousFlowVanGenuchten::dRelativePermeability(0.3, 0.7), 1.0E-5);
152  fd = (PorousFlowVanGenuchten::relativePermeability(0.8 + eps, 0.65) -
153  PorousFlowVanGenuchten::relativePermeability(0.8, 0.65)) /
154  eps;
155  EXPECT_NEAR(fd, PorousFlowVanGenuchten::dRelativePermeability(0.8, 0.65), 1.0E-5);
156 }
eps
const double eps
Definition: PorousFlowVanGenuchtenTest.C:14
PorousFlowVanGenuchten::dRelativePermeability
Real dRelativePermeability(Real seff, Real m)
Derivative of relative permeability with respect to effective saturation.
Definition: PorousFlowVanGenuchten.C:114
PorousFlowVanGenuchten::relativePermeability
T relativePermeability(const T &seff, Real m)
Relative permeability as a function of effective saturation.
Definition: PorousFlowVanGenuchten.h:101
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ TEST() [9/24]

TEST ( PorousFlowVanGenuchtenTest  ,
d2relperm   
)

Definition at line 158 of file PorousFlowVanGenuchtenTest.C.

159 {
160  Real fd;
161  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2RelativePermeability(1.0E30, 0.7), 1.0E-5);
162  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2RelativePermeability(-1.0, 0.7), 1.0E-5);
163  fd = (PorousFlowVanGenuchten::dRelativePermeability(0.3 + eps, 0.7) -
164  PorousFlowVanGenuchten::dRelativePermeability(0.3, 0.7)) /
165  eps;
166  EXPECT_NEAR(fd, PorousFlowVanGenuchten::d2RelativePermeability(0.3, 0.7), 1.0E-5);
167  fd = (PorousFlowVanGenuchten::dRelativePermeability(0.8 + eps, 0.65) -
168  PorousFlowVanGenuchten::dRelativePermeability(0.8, 0.65)) /
169  eps;
170  EXPECT_NEAR(fd, PorousFlowVanGenuchten::d2RelativePermeability(0.8, 0.65), 1.0E-5);
171 }
eps
const double eps
Definition: PorousFlowVanGenuchtenTest.C:14
PorousFlowVanGenuchten::dRelativePermeability
Real dRelativePermeability(Real seff, Real m)
Derivative of relative permeability with respect to effective saturation.
Definition: PorousFlowVanGenuchten.C:114
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
PorousFlowVanGenuchten::d2RelativePermeability
Real d2RelativePermeability(Real seff, Real m)
Second derivative of relative permeability with respect to effective saturation.
Definition: PorousFlowVanGenuchten.C:129

◆ TEST() [10/24]

TEST ( PorousFlowVanGenuchtenTest  ,
adrelperm   
)

Definition at line 173 of file PorousFlowVanGenuchtenTest.C.

174 {
175  ADReal sat = 0.3;
176  Moose::derivInsert(sat.derivatives(), 0, 1.0);
177 
178  const auto adrelperm = PorousFlowVanGenuchten::relativePermeability(sat, 0.7);
179 
180  const auto relperm = PorousFlowVanGenuchten::relativePermeability(sat.value(), 0.7);
181  const auto drelperm = PorousFlowVanGenuchten::dRelativePermeability(sat.value(), 0.7);
182 
183  EXPECT_NEAR(adrelperm.value(), relperm, 1.0E-5);
184  EXPECT_NEAR(adrelperm.derivatives()[0], drelperm, 1.0E-5);
185 }
ADReal
DualNumber< Real, DNDerivativeType, true > ADReal
PorousFlowVanGenuchten::dRelativePermeability
Real dRelativePermeability(Real seff, Real m)
Derivative of relative permeability with respect to effective saturation.
Definition: PorousFlowVanGenuchten.C:114
PorousFlowVanGenuchten::relativePermeability
T relativePermeability(const T &seff, Real m)
Relative permeability as a function of effective saturation.
Definition: PorousFlowVanGenuchten.h:101
Moose::derivInsert
void derivInsert(SemiDynamicSparseNumberArray< Real, libMesh::dof_id_type, NWrapper< N >> &derivs, libMesh::dof_id_type index, Real value)

◆ TEST() [11/24]

TEST ( PorousFlowVanGenuchtenTest  ,
relpermNW   
)

Definition at line 187 of file PorousFlowVanGenuchtenTest.C.

188 {
189  EXPECT_NEAR(1.0, PorousFlowVanGenuchten::relativePermeabilityNW(1.0E30, 0.7), 1.0E-5);
190  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::relativePermeabilityNW(-1.0, 0.7), 1.0E-5);
191  EXPECT_NEAR(0.664138097, PorousFlowVanGenuchten::relativePermeabilityNW(0.5, 0.3), 1.0E-5);
192  EXPECT_NEAR(0.559879012, PorousFlowVanGenuchten::relativePermeabilityNW(0.7, 0.8), 1.0E-5);
193 }
PorousFlowVanGenuchten::relativePermeabilityNW
T relativePermeabilityNW(const T &seff, Real m)
Relative permeability for a non-wetting phase as a function of effective saturation.
Definition: PorousFlowVanGenuchten.h:138

◆ TEST() [12/24]

TEST ( PorousFlowVanGenuchtenTest  ,
drelpermNW   
)

Definition at line 195 of file PorousFlowVanGenuchtenTest.C.

196 {
197  Real fd;
198  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dRelativePermeabilityNW(1.0E30, 0.7), 1.0E-5);
199  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dRelativePermeabilityNW(-1.0, 0.7), 1.0E-5);
200  fd = (PorousFlowVanGenuchten::relativePermeabilityNW(0.3 + eps, 0.7) -
201  PorousFlowVanGenuchten::relativePermeabilityNW(0.3, 0.7)) /
202  eps;
203  EXPECT_NEAR(fd, PorousFlowVanGenuchten::dRelativePermeabilityNW(0.3, 0.7), 1.0E-5);
204  fd = (PorousFlowVanGenuchten::relativePermeabilityNW(0.8 + eps, 0.65) -
205  PorousFlowVanGenuchten::relativePermeabilityNW(0.8, 0.65)) /
206  eps;
207  EXPECT_NEAR(fd, PorousFlowVanGenuchten::dRelativePermeabilityNW(0.8, 0.65), 1.0E-5);
208 }
PorousFlowVanGenuchten::relativePermeabilityNW
T relativePermeabilityNW(const T &seff, Real m)
Relative permeability for a non-wetting phase as a function of effective saturation.
Definition: PorousFlowVanGenuchten.h:138
eps
const double eps
Definition: PorousFlowVanGenuchtenTest.C:14
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
PorousFlowVanGenuchten::dRelativePermeabilityNW
Real dRelativePermeabilityNW(Real seff, Real m)
Derivative of relative permeability for a non-wetting phase with respect to effective saturation...
Definition: PorousFlowVanGenuchten.C:147

◆ TEST() [13/24]

TEST ( PorousFlowVanGenuchtenTest  ,
d2relpermNW   
)

Definition at line 210 of file PorousFlowVanGenuchtenTest.C.

211 {
212  Real fd;
213  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2RelativePermeabilityNW(1.0E30, 0.7), 1.0E-5);
214  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2RelativePermeabilityNW(-1.0, 0.7), 1.0E-5);
215  fd = (PorousFlowVanGenuchten::dRelativePermeabilityNW(0.3 + eps, 0.7) -
216  PorousFlowVanGenuchten::dRelativePermeabilityNW(0.3, 0.7)) /
217  eps;
218  EXPECT_NEAR(fd, PorousFlowVanGenuchten::d2RelativePermeabilityNW(0.3, 0.7), 1.0E-5);
219  fd = (PorousFlowVanGenuchten::dRelativePermeabilityNW(0.8 + eps, 0.65) -
220  PorousFlowVanGenuchten::dRelativePermeabilityNW(0.8, 0.65)) /
221  eps;
222  EXPECT_NEAR(fd, PorousFlowVanGenuchten::d2RelativePermeabilityNW(0.8, 0.65), 1.0E-5);
223 }
eps
const double eps
Definition: PorousFlowVanGenuchtenTest.C:14
PorousFlowVanGenuchten::d2RelativePermeabilityNW
Real d2RelativePermeabilityNW(Real seff, Real m)
Second derivative of relative permeability for a non-wetting phase with respect to effective saturati...
Definition: PorousFlowVanGenuchten.C:162
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
PorousFlowVanGenuchten::dRelativePermeabilityNW
Real dRelativePermeabilityNW(Real seff, Real m)
Derivative of relative permeability for a non-wetting phase with respect to effective saturation...
Definition: PorousFlowVanGenuchten.C:147

◆ TEST() [14/24]

TEST ( PorousFlowVanGenuchtenTest  ,
adrelpermnw   
)

Definition at line 225 of file PorousFlowVanGenuchtenTest.C.

226 {
227  ADReal sat = 0.3;
228  Moose::derivInsert(sat.derivatives(), 0, 1.0);
229 
230  const auto adrelperm = PorousFlowVanGenuchten::relativePermeabilityNW(sat, 0.7);
231 
232  const auto relperm = PorousFlowVanGenuchten::relativePermeabilityNW(sat.value(), 0.7);
233  const auto drelperm = PorousFlowVanGenuchten::dRelativePermeabilityNW(sat.value(), 0.7);
234 
235  EXPECT_NEAR(adrelperm.value(), relperm, 1.0E-5);
236  EXPECT_NEAR(adrelperm.derivatives()[0], drelperm, 1.0E-5);
237 }
PorousFlowVanGenuchten::relativePermeabilityNW
T relativePermeabilityNW(const T &seff, Real m)
Relative permeability for a non-wetting phase as a function of effective saturation.
Definition: PorousFlowVanGenuchten.h:138
ADReal
DualNumber< Real, DNDerivativeType, true > ADReal
Moose::derivInsert
void derivInsert(SemiDynamicSparseNumberArray< Real, libMesh::dof_id_type, NWrapper< N >> &derivs, libMesh::dof_id_type index, Real value)
PorousFlowVanGenuchten::dRelativePermeabilityNW
Real dRelativePermeabilityNW(Real seff, Real m)
Derivative of relative permeability for a non-wetting phase with respect to effective saturation...
Definition: PorousFlowVanGenuchten.C:147

◆ TEST() [15/24]

TEST ( PorousFlowVanGenuchtenTest  ,
caphys   
)

Test the hysteretic capillary pressure.

Definition at line 240 of file PorousFlowVanGenuchtenTest.C.

241 {
242  // first define some extensions that do not actually induce any extension, so the non-extended Pc
243  // can be checked
244  const std::vector<Real> no_ext_sats{1.1, 0.75, 0.2, 0.4, 0.6};
245  const std::vector<Real> expected_pcs{
246  0.0, 0.0, 0.7233512030263158, 0.18806139488299345, 0.06716785873826017};
247  for (unsigned i = 0; i < no_ext_sats.size(); ++i)
248  {
249  EXPECT_NEAR(expected_pcs[i],
250  PorousFlowVanGenuchten::capillaryPressureHys(
251  no_ext_sats[i], 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
252  1.0E-5);
253  EXPECT_NEAR(expected_pcs[i],
254  PorousFlowVanGenuchten::capillaryPressureHys(no_ext_sats[i], 0.1, 0.3, 10.0, 1.9),
255  1.0E-5);
256  }
257  EXPECT_NEAR(123.0,
258  PorousFlowVanGenuchten::capillaryPressureHys(
259  -1.0, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
260  1.0E-5);
261  EXPECT_NEAR(123.0,
262  PorousFlowVanGenuchten::capillaryPressureHys(
263  0.05, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
264  1.0E-5);
265  EXPECT_NEAR(std::numeric_limits<Real>::max(),
266  PorousFlowVanGenuchten::capillaryPressureHys(-1.0, 0.1, 0.3, 10.0, 1.9),
267  1.0E-5);
268  EXPECT_NEAR(std::numeric_limits<Real>::max(),
269  PorousFlowVanGenuchten::capillaryPressureHys(0.05, 0.1, 0.3, 10.0, 1.9),
270  1.0E-5);
271 
272  // check the lower extension values
273  EXPECT_NEAR(0.15229665668299805, low_ext_exp.S, 1.0E-5);
274  EXPECT_NEAR(-32.05516428734052, low_ext_exp.dPc, 1.0E-5);
275 
276  // check the high extension
277  EXPECT_NEAR(0.05300654102157442, high_ext_power.Pc, 1.0E-5);
278  EXPECT_NEAR(-0.48230528124844047, high_ext_power.dPc, 1.0E-5);
279 
280  // using low_ext_exp and high_ext_power
281  const std::vector<Real> sats{0.01, 0.1, 0.15, 0.2, 0.5, 0.63, 0.8, 0.99, 1.0};
282  std::vector<Real> pc;
283  pc = {31.385947046636815,
284  4.5861935551774735,
285  1.5754562501536364,
286  0.7233512030263158,
287  0.11727884570711045,
288  0.05300654102157442,
289  0.006681337544884095,
290  2.7847615834816514e-07,
291  0.0};
292  for (unsigned i = 0; i < sats.size(); ++i)
293  EXPECT_NEAR(pc[i],
294  PorousFlowVanGenuchten::capillaryPressureHys(
295  sats[i], 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power),
296  1.0E-5);
297 
298  // using low_ext_quad and high_ext_power
299  pc = {3.9304232526771696,
300  2.8885549236001244,
301  1.5730646091089062,
302  0.7233512030263158,
303  0.11727884570711045,
304  0.05300654102157442,
305  0.006681337544884095,
306  2.7847615834816514e-07,
307  0.0};
308  for (unsigned i = 0; i < sats.size(); ++i)
309  EXPECT_NEAR(pc[i],
310  PorousFlowVanGenuchten::capillaryPressureHys(
311  sats[i], 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power),
312  1.0E-5);
313 
314  // using low_ext_none and high_ext_power
315  pc = {1.5,
316  1.5,
317  1.5,
318  0.7233512030263158,
319  0.11727884570711045,
320  0.05300654102157442,
321  0.006681337544884095,
322  2.7847615834816514e-07,
323  0.0};
324  for (unsigned i = 0; i < sats.size(); ++i)
325  EXPECT_NEAR(pc[i],
326  PorousFlowVanGenuchten::capillaryPressureHys(
327  sats[i], 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power),
328  1.0E-5);
329 
330  // using low_ext_none and high_ext_none
331  pc = {1.5, 1.5, 1.5, 0.7233512030263158, 0.11727884570711045, 0.05300654102157442, 0.0, 0.0, 0.0};
332  for (unsigned i = 0; i < sats.size(); ++i)
333  EXPECT_NEAR(pc[i],
334  PorousFlowVanGenuchten::capillaryPressureHys(
335  sats[i], 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none),
336  1.0E-5);
337 }
low_ext_exp
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_exp(PorousFlowVanGenuchten::LowCapillaryPressureExtension::EXPONENTIAL, low_ext_S, low_ext_Pc, low_ext_dPc)
no_high_ext
const PorousFlowVanGenuchten::HighCapillaryPressureExtension no_high_ext(PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, 2.0, 0.0, 0.0)
no_low_ext
const PorousFlowVanGenuchten::LowCapillaryPressureExtension no_low_ext(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, -1.0, 123.0, 0.0)
high_ext_power
const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_power(PorousFlowVanGenuchten::HighCapillaryPressureExtension::POWER, high_ext_S, high_ext_Pc, high_ext_dPc)
PorousFlowVanGenuchten::capillaryPressureHys
Real capillaryPressureHys(Real sl, Real slmin, Real sgrdel, Real alpha, Real n, const LowCapillaryPressureExtension &low_ext=LowCapillaryPressureExtension(), const HighCapillaryPressureExtension &high_ext=HighCapillaryPressureExtension())
Hysteretic capillary pressure function (Eqn(1) of Doughty2007) with extensions (page5 and Fig1 of Dou...
Definition: PorousFlowVanGenuchten.C:180
high_ext_none
const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_none(PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, 0.8, high_ext_Pc, high_ext_dPc)
low_ext_none
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_none(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, low_ext_S, low_ext_Pc, low_ext_dPc)
low_ext_quad
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_quad(PorousFlowVanGenuchten::LowCapillaryPressureExtension::QUADRATIC, low_ext_S, low_ext_Pc, low_ext_dPc)

◆ TEST() [16/24]

TEST ( PorousFlowVanGenuchtenTest  ,
sathys   
)

Test the hysteretic saturation.

Definition at line 340 of file PorousFlowVanGenuchtenTest.C.

341 {
342  // first define some extensions that do not actually induce any extension, so the non-extended Pc
343  // can be checked
344  const std::vector<Real> expected_sats{1.0, 0.2, 0.4, 0.6};
345  const std::vector<Real> no_ext_pcs{
346  -1.1, 0.7233512030263158, 0.18806139488299345, 0.06716785873826017};
347  for (unsigned i = 0; i < no_ext_pcs.size(); ++i)
348  {
349  EXPECT_NEAR(expected_sats[i],
350  PorousFlowVanGenuchten::saturationHys(
351  no_ext_pcs[i], 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
352  1.0E-5);
353  EXPECT_NEAR(expected_sats[i],
354  PorousFlowVanGenuchten::saturationHys(no_ext_pcs[i], 0.1, 0.3, 10.0, 1.9),
355  1.0E-5);
356  }
357  EXPECT_NEAR(
358  -1.0,
359  PorousFlowVanGenuchten::saturationHys(1234.0, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
360  1.0E-5);
361  EXPECT_NEAR(
362  0.0,
363  PorousFlowVanGenuchten::saturationHys(std::numeric_limits<Real>::max(), 0.1, 0.3, 10.0, 1.9),
364  1.0E-5);
365 
366  // now with extensions
367  std::vector<Real> sats{0.01, 0.1, 0.15, 0.2, 0.5, 0.63, 0.8, 0.99, 1.0};
368  std::vector<Real> pc;
369  pc = {31.385947046636815,
370  4.5861935551774735,
371  1.5754562501536364,
372  0.7233512030263158,
373  0.11727884570711045,
374  0.05300654102157442,
375  0.006681337544884095,
376  2.7847615834816514e-07,
377  0.0};
378  for (unsigned i = 0; i < sats.size(); ++i)
379  EXPECT_NEAR(sats[i],
380  PorousFlowVanGenuchten::saturationHys(
381  pc[i], 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power),
382  1.0E-5);
383 
384  // different lower extension
385  pc = {3.9304232526771696,
386  2.8885549236001244,
387  1.5730646091089062,
388  0.7233512030263158,
389  0.11727884570711045,
390  0.05300654102157442,
391  0.006681337544884095,
392  2.7847615834816514e-07,
393  0.0};
394  for (unsigned i = 0; i < sats.size(); ++i)
395  EXPECT_NEAR(sats[i],
396  PorousFlowVanGenuchten::saturationHys(
397  pc[i], 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power),
398  1.0E-5);
399 
400  // different low extension
401  sats = {low_ext_none.S, low_ext_none.S, 0.2, 0.5, 0.63, 0.8, 0.99, 1.0};
402  pc = {2.5,
403  1.5,
404  0.7233512030263158,
405  0.11727884570711045,
406  0.05300654102157442,
407  0.006681337544884095,
408  2.7847615834816514e-07,
409  0.0};
410  for (unsigned i = 0; i < sats.size(); ++i)
411  EXPECT_NEAR(sats[i],
412  PorousFlowVanGenuchten::saturationHys(
413  pc[i], 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power),
414  1.0E-5);
415 
416  // different high extension
417  pc = {2.5, 1.5, 0.7233512030263158, 0.11727884570711045, 0.05300654102157442, 1E-10};
418  sats = {low_ext_none.S, low_ext_none.S, 0.2, 0.5, 0.63, 0.8};
419  for (unsigned i = 0; i < sats.size(); ++i)
420  EXPECT_NEAR(sats[i],
421  PorousFlowVanGenuchten::saturationHys(
422  pc[i], 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none),
423  1.0E-5);
424 }
low_ext_exp
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_exp(PorousFlowVanGenuchten::LowCapillaryPressureExtension::EXPONENTIAL, low_ext_S, low_ext_Pc, low_ext_dPc)
no_high_ext
const PorousFlowVanGenuchten::HighCapillaryPressureExtension no_high_ext(PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, 2.0, 0.0, 0.0)
no_low_ext
const PorousFlowVanGenuchten::LowCapillaryPressureExtension no_low_ext(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, -1.0, 123.0, 0.0)
high_ext_power
const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_power(PorousFlowVanGenuchten::HighCapillaryPressureExtension::POWER, high_ext_S, high_ext_Pc, high_ext_dPc)
high_ext_none
const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_none(PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, 0.8, high_ext_Pc, high_ext_dPc)
low_ext_none
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_none(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, low_ext_S, low_ext_Pc, low_ext_dPc)
PorousFlowVanGenuchten::saturationHys
Real saturationHys(Real pc, Real slmin, Real sgrdel, Real alpha, Real n, const LowCapillaryPressureExtension &low_ext=LowCapillaryPressureExtension(), const HighCapillaryPressureExtension &high_ext=HighCapillaryPressureExtension())
Hysteretic saturation function (Eqn(1) of Doughty2007) with extensions (page5 and Fig1 of Doughty2008...
Definition: PorousFlowVanGenuchten.C:364
low_ext_quad
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_quad(PorousFlowVanGenuchten::LowCapillaryPressureExtension::QUADRATIC, low_ext_S, low_ext_Pc, low_ext_dPc)

◆ TEST() [17/24]

TEST ( PorousFlowVanGenuchtenTest  ,
dcaphys   
)

Test the first derivative of the hysteretic capillary pressure.

Definition at line 427 of file PorousFlowVanGenuchtenTest.C.

428 {
429  const Real eps = 1.0E-9;
430 
431  // first define some extensions that do not actually induce any extension, so the non-extended
432  // Pc can be checked
433  std::vector<Real> sats{1.1, 0.75, -1.0, 0.05};
434  for (const auto & sat : sats)
435  {
436  EXPECT_NEAR(0.0,
437  PorousFlowVanGenuchten::dcapillaryPressureHys(
438  sat, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
439  1.0E-5);
440  EXPECT_NEAR(
441  0.0, PorousFlowVanGenuchten::dcapillaryPressureHys(sat, 0.1, 0.3, 10.0, 1.9), 1.0E-5);
442  }
443 
444  sats = {0.2, 0.4, 0.6};
445  for (const auto & sat : sats)
446  {
447  const Real fd = 0.5 *
448  (PorousFlowVanGenuchten::capillaryPressureHys(sat + eps, 0.1, 0.3, 10.0, 1.9) -
449  PorousFlowVanGenuchten::capillaryPressureHys(sat - eps, 0.1, 0.3, 10.0, 1.9)) /
450  eps;
451  EXPECT_NEAR(fd,
452  PorousFlowVanGenuchten::dcapillaryPressureHys(
453  sat, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
454  1.0E-5);
455  EXPECT_NEAR(
456  fd, PorousFlowVanGenuchten::dcapillaryPressureHys(sat, 0.1, 0.3, 10.0, 1.9), 1.0E-5);
457  };
458 
459  // now with extensions
460  sats = {0.01, 0.1, 0.15, 0.2, 0.5, 0.63, 0.8, 0.99, 1.0};
461 
462  for (const auto & sat : sats)
463  {
464  const Real fd = 0.5 *
465  (PorousFlowVanGenuchten::capillaryPressureHys(
466  sat + eps, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power) -
467  PorousFlowVanGenuchten::capillaryPressureHys(
468  sat - eps, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power)) /
469  eps;
470  EXPECT_NEAR(fd,
471  PorousFlowVanGenuchten::dcapillaryPressureHys(
472  sat, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power),
473  1.0E-5);
474  };
475 
476  // different low extension
477  for (const auto & sat : sats)
478  {
479  const Real fd = 0.5 *
480  (PorousFlowVanGenuchten::capillaryPressureHys(
481  sat + eps, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power) -
482  PorousFlowVanGenuchten::capillaryPressureHys(
483  sat - eps, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power)) /
484  eps;
485  EXPECT_NEAR(fd,
486  PorousFlowVanGenuchten::dcapillaryPressureHys(
487  sat, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power),
488  1.0E-5);
489  };
490 
491  // different lower extension
492  for (const auto & sat : sats)
493  {
494  const Real fd = 0.5 *
495  (PorousFlowVanGenuchten::capillaryPressureHys(
496  sat + eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power) -
497  PorousFlowVanGenuchten::capillaryPressureHys(
498  sat - eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power)) /
499  eps;
500  EXPECT_NEAR(fd,
501  PorousFlowVanGenuchten::dcapillaryPressureHys(
502  sat, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power),
503  1.0E-5);
504  };
505 
506  // different upper extension (cannot evaluate derivative at s = 0.63)
507  sats = {0.01, 0.1, 0.15, 0.2, 0.5, 0.8, 0.99, 1.0};
508  for (const auto & sat : sats)
509  {
510  const Real fd = 0.5 *
511  (PorousFlowVanGenuchten::capillaryPressureHys(
512  sat + eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none2) -
513  PorousFlowVanGenuchten::capillaryPressureHys(
514  sat - eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none2)) /
515  eps;
516  EXPECT_NEAR(fd,
517  PorousFlowVanGenuchten::dcapillaryPressureHys(
518  sat, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none2),
519  1.0E-5);
520  };
521 }
low_ext_exp
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_exp(PorousFlowVanGenuchten::LowCapillaryPressureExtension::EXPONENTIAL, low_ext_S, low_ext_Pc, low_ext_dPc)
PorousFlowVanGenuchten::dcapillaryPressureHys
Real dcapillaryPressureHys(Real sl, Real slmin, Real sgrdel, Real alpha, Real n, const LowCapillaryPressureExtension &low_ext=LowCapillaryPressureExtension(), const HighCapillaryPressureExtension &high_ext=HighCapillaryPressureExtension())
Derivative of capillaryPressureHys with respect to sl.
Definition: PorousFlowVanGenuchten.C:238
no_high_ext
const PorousFlowVanGenuchten::HighCapillaryPressureExtension no_high_ext(PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, 2.0, 0.0, 0.0)
no_low_ext
const PorousFlowVanGenuchten::LowCapillaryPressureExtension no_low_ext(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, -1.0, 123.0, 0.0)
eps
const double eps
Definition: PorousFlowVanGenuchtenTest.C:14
high_ext_power
const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_power(PorousFlowVanGenuchten::HighCapillaryPressureExtension::POWER, high_ext_S, high_ext_Pc, high_ext_dPc)
PorousFlowVanGenuchten::capillaryPressureHys
Real capillaryPressureHys(Real sl, Real slmin, Real sgrdel, Real alpha, Real n, const LowCapillaryPressureExtension &low_ext=LowCapillaryPressureExtension(), const HighCapillaryPressureExtension &high_ext=HighCapillaryPressureExtension())
Hysteretic capillary pressure function (Eqn(1) of Doughty2007) with extensions (page5 and Fig1 of Dou...
Definition: PorousFlowVanGenuchten.C:180
high_ext_none2
const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_none2(PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, high_ext_S, high_ext_Pc, high_ext_dPc)
low_ext_none
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_none(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, low_ext_S, low_ext_Pc, low_ext_dPc)
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
low_ext_quad
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_quad(PorousFlowVanGenuchten::LowCapillaryPressureExtension::QUADRATIC, low_ext_S, low_ext_Pc, low_ext_dPc)

◆ TEST() [18/24]

TEST ( PorousFlowVanGenuchtenTest  ,
d2caphys   
)

Test the second derivative of the hysteretic capillary pressure.

Definition at line 524 of file PorousFlowVanGenuchtenTest.C.

525 {
526  const Real eps = 1.0E-9;
527 
528  // first define some extensions that do not actually induce any extension, so the non-extended
529  // Pc can be checked
530  std::vector<Real> sats{1.1, 0.75, -1.0, 0.05};
531  for (const auto & sat : sats)
532  {
533  EXPECT_NEAR(0.0,
534  PorousFlowVanGenuchten::d2capillaryPressureHys(
535  sat, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
536  1.0E-5);
537  EXPECT_NEAR(
538  0.0, PorousFlowVanGenuchten::d2capillaryPressureHys(sat, 0.1, 0.3, 10.0, 1.9), 1.0E-5);
539  }
540 
541  sats = {0.2, 0.4, 0.6};
542  for (const auto & sat : sats)
543  {
544  const Real fd =
545  0.5 *
546  (PorousFlowVanGenuchten::dcapillaryPressureHys(sat + eps, 0.1, 0.3, 10.0, 1.9) -
547  PorousFlowVanGenuchten::dcapillaryPressureHys(sat - eps, 0.1, 0.3, 10.0, 1.9)) /
548  eps;
549  EXPECT_NEAR(fd,
550  PorousFlowVanGenuchten::d2capillaryPressureHys(
551  sat, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
552  1.0E-5);
553  EXPECT_NEAR(
554  fd, PorousFlowVanGenuchten::d2capillaryPressureHys(sat, 0.1, 0.3, 10.0, 1.9), 1.0E-5);
555  };
556 
557  // now with extensions
558  sats = {0.01, 0.1, 0.15, 0.2, 0.5, 0.6, 0.7, 0.8, 0.99, 1.0};
559 
560  for (const auto & sat : sats)
561  {
562  const Real fd = 0.5 *
563  (PorousFlowVanGenuchten::dcapillaryPressureHys(
564  sat + eps, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power) -
565  PorousFlowVanGenuchten::dcapillaryPressureHys(
566  sat - eps, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power)) /
567  eps;
568  if (std::abs(fd) > 10)
569  EXPECT_NEAR(1.0,
570  fd / PorousFlowVanGenuchten::d2capillaryPressureHys(
571  sat, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power),
572  1.0E-5);
573  else
574  EXPECT_NEAR(fd,
575  PorousFlowVanGenuchten::d2capillaryPressureHys(
576  sat, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power),
577  1.0E-5);
578  };
579 
580  // different low extension
581  for (const auto & sat : sats)
582  {
583  const Real fd = 0.5 *
584  (PorousFlowVanGenuchten::dcapillaryPressureHys(
585  sat + eps, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power) -
586  PorousFlowVanGenuchten::dcapillaryPressureHys(
587  sat - eps, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power)) /
588  eps;
589  if (std::abs(fd) > 10)
590  EXPECT_NEAR(1.0,
591  fd / PorousFlowVanGenuchten::d2capillaryPressureHys(
592  sat, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power),
593  1.0E-5);
594  else
595  EXPECT_NEAR(fd,
596  PorousFlowVanGenuchten::d2capillaryPressureHys(
597  sat, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power),
598  1.0E-5);
599  };
600 
601  // different lower extension
602  for (const auto & sat : sats)
603  {
604  const Real fd = 0.5 *
605  (PorousFlowVanGenuchten::dcapillaryPressureHys(
606  sat + eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power) -
607  PorousFlowVanGenuchten::dcapillaryPressureHys(
608  sat - eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power)) /
609  eps;
610  if (std::abs(fd) > 10)
611  EXPECT_NEAR(1.0,
612  fd / PorousFlowVanGenuchten::d2capillaryPressureHys(
613  sat, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power),
614  1.0E-5);
615  else
616  EXPECT_NEAR(fd,
617  PorousFlowVanGenuchten::d2capillaryPressureHys(
618  sat, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power),
619  1.0E-5);
620  };
621 
622  // different upper extension (cannot evaluate derivative at s = 0.63)
623  sats = {0.01, 0.1, 0.15, 0.2, 0.5, 0.8, 0.99, 1.0};
624  for (const auto & sat : sats)
625  {
626  const Real fd = 0.5 *
627  (PorousFlowVanGenuchten::dcapillaryPressureHys(
628  sat + eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none2) -
629  PorousFlowVanGenuchten::dcapillaryPressureHys(
630  sat - eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none2)) /
631  eps;
632  if (std::abs(fd) > 10)
633  EXPECT_NEAR(1.0,
634  fd / PorousFlowVanGenuchten::d2capillaryPressureHys(
635  sat, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none2),
636  1.0E-5);
637  else
638  EXPECT_NEAR(fd,
639  PorousFlowVanGenuchten::d2capillaryPressureHys(
640  sat, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none2),
641  1.0E-5);
642  };
643 }
low_ext_exp
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_exp(PorousFlowVanGenuchten::LowCapillaryPressureExtension::EXPONENTIAL, low_ext_S, low_ext_Pc, low_ext_dPc)
PorousFlowVanGenuchten::dcapillaryPressureHys
Real dcapillaryPressureHys(Real sl, Real slmin, Real sgrdel, Real alpha, Real n, const LowCapillaryPressureExtension &low_ext=LowCapillaryPressureExtension(), const HighCapillaryPressureExtension &high_ext=HighCapillaryPressureExtension())
Derivative of capillaryPressureHys with respect to sl.
Definition: PorousFlowVanGenuchten.C:238
PorousFlowVanGenuchten::d2capillaryPressureHys
Real d2capillaryPressureHys(Real sl, Real slmin, Real sgrdel, Real alpha, Real n, const LowCapillaryPressureExtension &low_ext=LowCapillaryPressureExtension(), const HighCapillaryPressureExtension &high_ext=HighCapillaryPressureExtension())
Second derivative of capillaryPressureHys with respect to sl.
Definition: PorousFlowVanGenuchten.C:299
no_high_ext
const PorousFlowVanGenuchten::HighCapillaryPressureExtension no_high_ext(PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, 2.0, 0.0, 0.0)
no_low_ext
const PorousFlowVanGenuchten::LowCapillaryPressureExtension no_low_ext(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, -1.0, 123.0, 0.0)
eps
const double eps
Definition: PorousFlowVanGenuchtenTest.C:14
high_ext_power
const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_power(PorousFlowVanGenuchten::HighCapillaryPressureExtension::POWER, high_ext_S, high_ext_Pc, high_ext_dPc)
high_ext_none2
const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_none2(PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, high_ext_S, high_ext_Pc, high_ext_dPc)
low_ext_none
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_none(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, low_ext_S, low_ext_Pc, low_ext_dPc)
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
low_ext_quad
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_quad(PorousFlowVanGenuchten::LowCapillaryPressureExtension::QUADRATIC, low_ext_S, low_ext_Pc, low_ext_dPc)

◆ TEST() [19/24]

TEST ( PorousFlowVanGenuchtenTest  ,
dsathys   
)

Test the derivative of the hysteretic saturation.

Definition at line 646 of file PorousFlowVanGenuchtenTest.C.

647 {
648  const Real eps = 1E-9;
649 
650  // first define some extensions that do not actually induce any extension, so the non-extended Pc
651  // can be checked
652  std::vector<Real> pcs{-1.1, 1234.0};
653  for (const auto & pc : pcs)
654  {
655  EXPECT_NEAR(
656  0.0,
657  PorousFlowVanGenuchten::dsaturationHys(pc, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
658  1.0E-5);
659  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dsaturationHys(pc, 0.1, 0.3, 10.0, 1.9), 1.0E-5);
660  }
661 
662  pcs = {0.7, 0.2, 0.07};
663  for (const auto & pc : pcs)
664  {
665  const Real fd = 0.5 *
666  (PorousFlowVanGenuchten::saturationHys(
667  pc + eps, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext) -
668  PorousFlowVanGenuchten::saturationHys(
669  pc - eps, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext)) /
670  eps;
671  EXPECT_NEAR(
672  fd,
673  PorousFlowVanGenuchten::dsaturationHys(pc, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
674  1.0E-5);
675  EXPECT_NEAR(fd, PorousFlowVanGenuchten::dsaturationHys(pc, 0.1, 0.3, 10.0, 1.9), 1.0E-5);
676  }
677 
678  // now with low and high extensions
679  pcs = {31.385947046636815,
680  4.5861935551774735,
681  1.5754562501536364,
682  0.7233512030263158,
683  0.11727884570711045,
684  0.05300654102157442,
685  0.006681337544884095};
686  for (const auto & pc : pcs)
687  {
688  const Real fd = 0.5 *
689  (PorousFlowVanGenuchten::saturationHys(
690  pc + eps, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power) -
691  PorousFlowVanGenuchten::saturationHys(
692  pc - eps, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power)) /
693  eps;
694  EXPECT_NEAR(fd,
695  PorousFlowVanGenuchten::dsaturationHys(
696  pc, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power),
697  1.0E-5);
698  }
699 
700  // different lower extension
701  pcs = {3.9304232526771696,
702  2.8885549236001244,
703  1.5730646091089062,
704  0.7233512030263158,
705  0.11727884570711045,
706  0.05300654102157442,
707  0.006681337544884095};
708  for (const auto & pc : pcs)
709  {
710  const Real fd = 0.5 *
711  (PorousFlowVanGenuchten::saturationHys(
712  pc + eps, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power) -
713  PorousFlowVanGenuchten::saturationHys(
714  pc - eps, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power)) /
715  eps;
716  EXPECT_NEAR(fd,
717  PorousFlowVanGenuchten::dsaturationHys(
718  pc, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power),
719  1.0E-5);
720  }
721 
722  // different lower extension
723  pcs = {2.5, 0.7233512030263158, 0.11727884570711045, 0.05300654102157442, 0.006681337544884095};
724  for (const auto & pc : pcs)
725  {
726  const Real fd = 0.5 *
727  (PorousFlowVanGenuchten::saturationHys(
728  pc + eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power) -
729  PorousFlowVanGenuchten::saturationHys(
730  pc - eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power)) /
731  eps;
732  EXPECT_NEAR(fd,
733  PorousFlowVanGenuchten::dsaturationHys(
734  pc, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power),
735  1.0E-5);
736  }
737 
738  // different upper extension
739  pcs = {2.5, 0.7233512030263158, 0.11727884570711045, 0.05, 0.04};
740  for (const auto & pc : pcs)
741  {
742  const Real fd = 0.5 *
743  (PorousFlowVanGenuchten::saturationHys(
744  pc + eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none) -
745  PorousFlowVanGenuchten::saturationHys(
746  pc - eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none)) /
747  eps;
748  EXPECT_NEAR(fd,
749  PorousFlowVanGenuchten::dsaturationHys(
750  pc, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none),
751  1.0E-5);
752  }
753 }
low_ext_exp
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_exp(PorousFlowVanGenuchten::LowCapillaryPressureExtension::EXPONENTIAL, low_ext_S, low_ext_Pc, low_ext_dPc)
no_high_ext
const PorousFlowVanGenuchten::HighCapillaryPressureExtension no_high_ext(PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, 2.0, 0.0, 0.0)
no_low_ext
const PorousFlowVanGenuchten::LowCapillaryPressureExtension no_low_ext(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, -1.0, 123.0, 0.0)
eps
const double eps
Definition: PorousFlowVanGenuchtenTest.C:14
high_ext_power
const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_power(PorousFlowVanGenuchten::HighCapillaryPressureExtension::POWER, high_ext_S, high_ext_Pc, high_ext_dPc)
high_ext_none
const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_none(PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, 0.8, high_ext_Pc, high_ext_dPc)
low_ext_none
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_none(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, low_ext_S, low_ext_Pc, low_ext_dPc)
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
PorousFlowVanGenuchten::dsaturationHys
Real dsaturationHys(Real pc, Real slmin, Real sgrdel, Real alpha, Real n, const LowCapillaryPressureExtension &low_ext=LowCapillaryPressureExtension(), const HighCapillaryPressureExtension &high_ext=HighCapillaryPressureExtension())
Derivative of Hysteretic saturation function with respect to pc.
Definition: PorousFlowVanGenuchten.C:433
PorousFlowVanGenuchten::saturationHys
Real saturationHys(Real pc, Real slmin, Real sgrdel, Real alpha, Real n, const LowCapillaryPressureExtension &low_ext=LowCapillaryPressureExtension(), const HighCapillaryPressureExtension &high_ext=HighCapillaryPressureExtension())
Hysteretic saturation function (Eqn(1) of Doughty2007) with extensions (page5 and Fig1 of Doughty2008...
Definition: PorousFlowVanGenuchten.C:364
low_ext_quad
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_quad(PorousFlowVanGenuchten::LowCapillaryPressureExtension::QUADRATIC, low_ext_S, low_ext_Pc, low_ext_dPc)

◆ TEST() [20/24]

TEST ( PorousFlowVanGenuchtenTest  ,
d2sathys   
)

Test the second derivative of the hysteretic saturation.

Definition at line 756 of file PorousFlowVanGenuchtenTest.C.

757 {
758  const Real eps = 1E-9;
759 
760  // first define some extensions that do not actually induce any extension, so the non-extended Pc
761  // can be checked
762  std::vector<Real> pcs{-1.1, 1234.0};
763  for (const auto & pc : pcs)
764  {
765  EXPECT_NEAR(
766  0.0,
767  PorousFlowVanGenuchten::d2saturationHys(pc, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
768  1.0E-5);
769  EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2saturationHys(pc, 0.1, 0.3, 10.0, 1.9), 1.0E-5);
770  }
771 
772  pcs = {0.7, 0.2, 0.07};
773  for (const auto & pc : pcs)
774  {
775  const Real fd = 0.5 *
776  (PorousFlowVanGenuchten::dsaturationHys(
777  pc + eps, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext) -
778  PorousFlowVanGenuchten::dsaturationHys(
779  pc - eps, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext)) /
780  eps;
781  EXPECT_NEAR(
782  fd,
783  PorousFlowVanGenuchten::d2saturationHys(pc, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
784  1.0E-5);
785  EXPECT_NEAR(fd, PorousFlowVanGenuchten::d2saturationHys(pc, 0.1, 0.3, 10.0, 1.9), 1.0E-5);
786  }
787 
788  // now with low and high extensions
789  pcs = {31.385947046636815,
790  4.5861935551774735,
791  1.5754562501536364,
792  0.7233512030263158,
793  0.11727884570711045,
794  0.06300654102157442,
795  0.006681337544884095};
796  for (const auto & pc : pcs)
797  {
798  const Real fd = 0.5 *
799  (PorousFlowVanGenuchten::dsaturationHys(
800  pc + eps, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power) -
801  PorousFlowVanGenuchten::dsaturationHys(
802  pc - eps, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power)) /
803  eps;
804  EXPECT_NEAR(fd,
805  PorousFlowVanGenuchten::d2saturationHys(
806  pc, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power),
807  1.0E-5);
808  }
809 
810  // different lower extension
811  pcs = {3.9304232526771696,
812  2.8885549236001244,
813  1.5730646091089062,
814  0.7233512030263158,
815  0.11727884570711045,
816  0.06300654102157442,
817  0.006681337544884095};
818  for (const auto & pc : pcs)
819  {
820  const Real fd = 0.5 *
821  (PorousFlowVanGenuchten::dsaturationHys(
822  pc + eps, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power) -
823  PorousFlowVanGenuchten::dsaturationHys(
824  pc - eps, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power)) /
825  eps;
826  EXPECT_NEAR(fd,
827  PorousFlowVanGenuchten::d2saturationHys(
828  pc, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power),
829  1.0E-5);
830  }
831 
832  // different lower extension
833  pcs = {2.5, 0.7233512030263158, 0.11727884570711045, 0.06300654102157442, 0.006681337544884095};
834  for (const auto & pc : pcs)
835  {
836  const Real fd = 0.5 *
837  (PorousFlowVanGenuchten::dsaturationHys(
838  pc + eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power) -
839  PorousFlowVanGenuchten::dsaturationHys(
840  pc - eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power)) /
841  eps;
842  EXPECT_NEAR(fd,
843  PorousFlowVanGenuchten::d2saturationHys(
844  pc, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power),
845  1.0E-5);
846  }
847 
848  // different upper extension
849  pcs = {2.5, 0.7233512030263158, 0.11727884570711045, 0.05, 0.04};
850  for (const auto & pc : pcs)
851  {
852  const Real fd = 0.5 *
853  (PorousFlowVanGenuchten::dsaturationHys(
854  pc + eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none) -
855  PorousFlowVanGenuchten::dsaturationHys(
856  pc - eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none)) /
857  eps;
858  EXPECT_NEAR(fd,
859  PorousFlowVanGenuchten::d2saturationHys(
860  pc, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none),
861  1.0E-5);
862  }
863 }
low_ext_exp
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_exp(PorousFlowVanGenuchten::LowCapillaryPressureExtension::EXPONENTIAL, low_ext_S, low_ext_Pc, low_ext_dPc)
no_high_ext
const PorousFlowVanGenuchten::HighCapillaryPressureExtension no_high_ext(PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, 2.0, 0.0, 0.0)
no_low_ext
const PorousFlowVanGenuchten::LowCapillaryPressureExtension no_low_ext(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, -1.0, 123.0, 0.0)
eps
const double eps
Definition: PorousFlowVanGenuchtenTest.C:14
high_ext_power
const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_power(PorousFlowVanGenuchten::HighCapillaryPressureExtension::POWER, high_ext_S, high_ext_Pc, high_ext_dPc)
high_ext_none
const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_none(PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, 0.8, high_ext_Pc, high_ext_dPc)
low_ext_none
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_none(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, low_ext_S, low_ext_Pc, low_ext_dPc)
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
PorousFlowVanGenuchten::dsaturationHys
Real dsaturationHys(Real pc, Real slmin, Real sgrdel, Real alpha, Real n, const LowCapillaryPressureExtension &low_ext=LowCapillaryPressureExtension(), const HighCapillaryPressureExtension &high_ext=HighCapillaryPressureExtension())
Derivative of Hysteretic saturation function with respect to pc.
Definition: PorousFlowVanGenuchten.C:433
PorousFlowVanGenuchten::d2saturationHys
Real d2saturationHys(Real pc, Real slmin, Real sgrdel, Real alpha, Real n, const LowCapillaryPressureExtension &low_ext=LowCapillaryPressureExtension(), const HighCapillaryPressureExtension &high_ext=HighCapillaryPressureExtension())
Second derivative of Hysteretic saturation function with respect to pc.
Definition: PorousFlowVanGenuchten.C:507
low_ext_quad
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_quad(PorousFlowVanGenuchten::LowCapillaryPressureExtension::QUADRATIC, low_ext_S, low_ext_Pc, low_ext_dPc)

◆ TEST() [21/24]

TEST ( PorousFlowVanGenuchtenTest  ,
relativePermeabilityHys   
)

Test relativePermeabilityHys.

Definition at line 866 of file PorousFlowVanGenuchtenTest.C.

867 {
868  // Tests for sldel = 0.5 >= slr = 0.2
869  EXPECT_NEAR(0.0,
870  PorousFlowVanGenuchten::relativePermeabilityHys(
871  0.1, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
872  1.0E-12); // zero because sl < slr
873  EXPECT_NEAR(0.0,
874  PorousFlowVanGenuchten::relativePermeabilityHys(
875  0.2, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
876  1.0E-12); // zero because sl = slr
877  EXPECT_NEAR(0.002847974503642625,
878  PorousFlowVanGenuchten::relativePermeabilityHys(
879  0.3, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
880  1.0E-12); // drying because sl <= sldel
881  EXPECT_NEAR(0.002847974503642625,
882  PorousFlowVanGenuchten::relativePermeabilityHys(
883  0.3, 0.2, 0.0, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
884  1.0E-12); // drying because sgrdel = 0
885  EXPECT_NEAR(0.05828443549816974,
886  PorousFlowVanGenuchten::relativePermeabilityHys(
887  0.5, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
888  1.0E-12); // drying because sl <= sldel
889  EXPECT_NEAR(0.10056243969693253,
890  PorousFlowVanGenuchten::relativePermeabilityHys(
891  0.55, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
892  1.0E-12); // wetting
893  EXPECT_NEAR(0.08943153676247108,
894  PorousFlowVanGenuchten::relativePermeabilityHys(
895  0.55, 0.2, 0.0, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
896  1.0E-12); // drying because sgrdel = 0
897  EXPECT_NEAR(0.3159149169921876,
898  PorousFlowVanGenuchten::relativePermeabilityHys(
899  0.8, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
900  1.0E-12); // cubic
901  EXPECT_NEAR(0.8011290515690178,
902  PorousFlowVanGenuchten::relativePermeabilityHys(
903  0.95, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
904  1.0E-12); // drying, because sl > cubic modification region
905 
906  // Tests for sldel = 0.15 < slr = 0.2
907  EXPECT_NEAR(0.0,
908  PorousFlowVanGenuchten::relativePermeabilityHys(
909  0.1, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
910  1.0E-12); // zero because sl < slr
911  EXPECT_NEAR(0.0,
912  PorousFlowVanGenuchten::relativePermeabilityHys(
913  0.2, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
914  1.0E-12); // zero because sl = slr
915  EXPECT_NEAR(0.005858393312913491,
916  PorousFlowVanGenuchten::relativePermeabilityHys(
917  0.3, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
918  1.0E-12); // wetting
919  EXPECT_NEAR(0.002847974503642625,
920  PorousFlowVanGenuchten::relativePermeabilityHys(
921  0.3, 0.2, 0.0, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
922  1.0E-12); // drying because sgrdel = 0
923  EXPECT_NEAR(0.1363562523627124,
924  PorousFlowVanGenuchten::relativePermeabilityHys(
925  0.55, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
926  1.0E-12); // wetting
927  EXPECT_NEAR(0.31,
928  PorousFlowVanGenuchten::relativePermeabilityHys(
929  0.7, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
930  1.0E-12); // cubic
931  EXPECT_NEAR(0.8011290515690178,
932  PorousFlowVanGenuchten::relativePermeabilityHys(
933  0.95, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
934  1.0E-12); // drying, because sl > cubic modification region
935 }
PorousFlowVanGenuchten::relativePermeabilityHys
Real relativePermeabilityHys(Real sl, Real slr, Real sgrdel, Real sgrmax, Real sldel, Real m, Real upper_liquid_param, Real y0, Real y0p, Real y1, Real y1p)
Hysteretic relative permeability for liquid.
Definition: PorousFlowVanGenuchten.C:586

◆ TEST() [22/24]

TEST ( PorousFlowVanGenuchtenTest  ,
drelativePermeabilityHys   
)

Test drelativePermeabilityHys.

Definition at line 938 of file PorousFlowVanGenuchtenTest.C.

939 {
940  // essentially follow the testing strategy in the relativePermeabilityHys Test, but do not test
941  // the derivative at the points where it is discontinuous
942  EXPECT_NEAR(0.0,
943  PorousFlowVanGenuchten::drelativePermeabilityHys(
944  0.1, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
945  1.0E-12);
946  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
947  0.3, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
948  (PorousFlowVanGenuchten::relativePermeabilityHys(
949  0.3 + eps, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
950  PorousFlowVanGenuchten::relativePermeabilityHys(
951  0.3, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
952  eps,
953  1.0E-5);
954  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
955  0.3, 0.2, 0.0, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
956  (PorousFlowVanGenuchten::relativePermeabilityHys(
957  0.3 + eps, 0.2, 0.0, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
958  PorousFlowVanGenuchten::relativePermeabilityHys(
959  0.3, 0.2, 0.0, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
960  eps,
961  1.0E-5);
962  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
963  0.55, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
964  (PorousFlowVanGenuchten::relativePermeabilityHys(
965  0.55 + eps, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
966  PorousFlowVanGenuchten::relativePermeabilityHys(
967  0.55, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
968  eps,
969  1.0E-5);
970  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
971  0.55, 0.2, 0.0, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
972  (PorousFlowVanGenuchten::relativePermeabilityHys(
973  0.55 + eps, 0.2, 0.0, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
974  PorousFlowVanGenuchten::relativePermeabilityHys(
975  0.55, 0.2, 0.0, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
976  eps,
977  1.0E-5);
978  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
979  0.8, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
980  (PorousFlowVanGenuchten::relativePermeabilityHys(
981  0.8 + eps, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
982  PorousFlowVanGenuchten::relativePermeabilityHys(
983  0.8, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
984  eps,
985  1.0E-5);
986  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
987  0.95, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
988  (PorousFlowVanGenuchten::relativePermeabilityHys(
989  0.95 + eps, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
990  PorousFlowVanGenuchten::relativePermeabilityHys(
991  0.95, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
992  eps,
993  1.0E-5);
994 
995  EXPECT_NEAR(0.0,
996  PorousFlowVanGenuchten::drelativePermeabilityHys(
997  0.1, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
998  1.0E-5);
999  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
1000  0.3, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
1001  (PorousFlowVanGenuchten::relativePermeabilityHys(
1002  0.3 + eps, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
1003  PorousFlowVanGenuchten::relativePermeabilityHys(
1004  0.3, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
1005  eps,
1006  1.0E-5);
1007  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
1008  0.3, 0.2, 0.0, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
1009  (PorousFlowVanGenuchten::relativePermeabilityHys(
1010  0.3 + eps, 0.2, 0.0, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
1011  PorousFlowVanGenuchten::relativePermeabilityHys(
1012  0.3, 0.2, 0.0, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
1013  eps,
1014  1.0E-5);
1015  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
1016  0.55, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
1017  (PorousFlowVanGenuchten::relativePermeabilityHys(
1018  0.55 + eps, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
1019  PorousFlowVanGenuchten::relativePermeabilityHys(
1020  0.55, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
1021  eps,
1022  1.0E-5);
1023  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
1024  0.7, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
1025  (PorousFlowVanGenuchten::relativePermeabilityHys(
1026  0.7 + eps, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
1027  PorousFlowVanGenuchten::relativePermeabilityHys(
1028  0.7, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
1029  eps,
1030  1.0E-5);
1031  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
1032  0.95, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
1033  (PorousFlowVanGenuchten::relativePermeabilityHys(
1034  0.95 + eps, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
1035  PorousFlowVanGenuchten::relativePermeabilityHys(
1036  0.95, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
1037  eps,
1038  1.0E-5);
1039 }
PorousFlowVanGenuchten::drelativePermeabilityHys
Real drelativePermeabilityHys(Real sl, Real slr, Real sgrdel, Real sgrmax, Real sldel, Real m, Real upper_liquid_param, Real y0, Real y0p, Real y1, Real y1p)
Derivative of Hysteretic relative permeability for liquid, with respect to liquid saturation...
Definition: PorousFlowVanGenuchten.C:642
eps
const double eps
Definition: PorousFlowVanGenuchtenTest.C:14
PorousFlowVanGenuchten::relativePermeabilityHys
Real relativePermeabilityHys(Real sl, Real slr, Real sgrdel, Real sgrmax, Real sldel, Real m, Real upper_liquid_param, Real y0, Real y0p, Real y1, Real y1p)
Hysteretic relative permeability for liquid.
Definition: PorousFlowVanGenuchten.C:586

◆ TEST() [23/24]

TEST ( PorousFlowVanGenuchtenTest  ,
relativePermeabilityNWHys   
)

Test relativePermeabilityNWHys.

Definition at line 1042 of file PorousFlowVanGenuchtenTest.C.

1043 {
1044  // Tests for sldel = 0.5 >= slr = 0.2
1045  EXPECT_NEAR(0.91875,
1046  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1047  0.1, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1048  1.0E-12); // cubic extension region because sl < slr
1049  EXPECT_NEAR(1.0,
1050  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1051  0.1, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 1.0, -0.75),
1052  1.0E-12); // cubic extension region because sl < slr, but k_rg_max = 1
1053  EXPECT_NEAR(0.8,
1054  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1055  0.2, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1056  1.0E-12); // at boundary of extension since sl = slr
1057  EXPECT_NEAR(0.6368139633459233,
1058  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1059  0.3, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1060  1.0E-12); // drying because sl <= sldel
1061  EXPECT_NEAR(0.6368139633459233,
1062  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1063  0.3, 0.2, 0.0, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1064  1.0E-12); // drying because sgrdel = 0
1065  EXPECT_NEAR(0.33222054246699,
1066  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1067  0.5, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1068  1.0E-12); // drying because sl <= sldel
1069  EXPECT_NEAR(0.24397255389213568,
1070  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1071  0.55, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1072  1.0E-12); // wetting
1073  EXPECT_NEAR(0.2691316236913443,
1074  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1075  0.55, 0.2, 0.0, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1076  1.0E-12); // drying because sgrdel = 0
1077  EXPECT_NEAR(0.01509148277290631,
1078  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1079  0.85, 0.2, 0.05, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1080  1.0E-12); // wetting
1081  EXPECT_NEAR(0.0,
1082  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1083  0.85, 0.2, 0.18, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1084  1.0E-12); // sl > 1 - s_gr_del
1085 
1086  // Tests for sldel = 0.15 < slr = 0.2
1087  EXPECT_NEAR(0.91875,
1088  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1089  0.1, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1090  1.0E-12); // cubic extension region because sl < slr
1091  EXPECT_NEAR(1.0,
1092  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1093  0.1, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 1.0, -0.75),
1094  1.0E-12); // cubic extension region because sl < slr, but k_rg_max = 1
1095  EXPECT_NEAR(0.8,
1096  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1097  0.2, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1098  1.0E-12); // at boundary of extension since sl = slr
1099  EXPECT_NEAR(0.5616820966427329,
1100  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1101  0.3, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1102  1.0E-12); // wetting
1103  EXPECT_NEAR(0.6368139633459233,
1104  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1105  0.3, 0.2, 0.0, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1106  1.0E-12); // wetting = drying because sgrdel = 0
1107  EXPECT_NEAR(0.11086138435665097,
1108  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1109  0.55, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1110  1.0E-12); // wetting
1111  EXPECT_NEAR(0.0,
1112  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1113  0.76, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1114  1.0E-12); // always zero for large saturation
1115 }
PorousFlowVanGenuchten::relativePermeabilityNWHys
Real relativePermeabilityNWHys(Real sl, Real slr, Real sgrdel, Real sgrmax, Real sldel, Real m, Real gamma, Real k_rg_max, Real y0p)
Hysteretic relative permeability for gas.
Definition: PorousFlowVanGenuchten.C:720

◆ TEST() [24/24]

TEST ( PorousFlowVanGenuchtenTest  ,
drelativePermeabilityNWHys   
)

Test drelativePermeabilityNWHys.

Definition at line 1118 of file PorousFlowVanGenuchtenTest.C.

1119 {
1120  // Testing strategy follows the relativePermeabilityNWHys strategy, except the points at which the
1121  // derivative is not defined are not tested
1122  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1123  0.1, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1124  (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1125  0.1 + eps, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75) -
1126  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1127  0.1, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75)) /
1128  eps,
1129  1.0E-5); // cubic extension region because sl < slr
1130  EXPECT_NEAR(0.0,
1131  PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1132  0.1, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 1.0, -0.75),
1133  1.0E-12); // cubic extension region because sl < slr, but k_rg_max = 1
1134  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1135  0.3, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1136  (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1137  0.3 + eps, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75) -
1138  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1139  0.3, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75)) /
1140  eps,
1141  1.0E-5); // drying because sl <= sldel
1142  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1143  0.3, 0.2, 0.0, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1144  (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1145  0.3 + eps, 0.2, 0.0, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75) -
1146  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1147  0.3, 0.2, 0.0, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75)) /
1148  eps,
1149  1.0E-5); // drying because sgrdel = 0
1150  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1151  0.45, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1152  (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1153  0.45 + eps, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75) -
1154  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1155  0.45, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75)) /
1156  eps,
1157  1.0E-5); // drying because sl <= sldel
1158  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1159  0.55, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1160  (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1161  0.55 + eps, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75) -
1162  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1163  0.55, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75)) /
1164  eps,
1165  1.0E-5); // wetting
1166  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1167  0.55, 0.2, 0.0, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1168  (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1169  0.55 + eps, 0.2, 0.0, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75) -
1170  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1171  0.55, 0.2, 0.0, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75)) /
1172  eps,
1173  1.0E-5); // drying because sgrdel = 0
1174  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1175  0.85, 0.2, 0.05, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1176  (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1177  0.85 + eps, 0.2, 0.05, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75) -
1178  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1179  0.85, 0.2, 0.05, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75)) /
1180  eps,
1181  1.0E-5); // wetting
1182  EXPECT_NEAR(0.0,
1183  PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1184  0.85, 0.2, 0.18, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1185  1.0E-12); // sl > 1 - s_gr_del
1186 
1187  // Tests for sldel = 0.15 < slr = 0.2
1188  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1189  0.1, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1190  (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1191  0.1 + eps, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75) -
1192  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1193  0.1, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75)) /
1194  eps,
1195  1.0E-5); // cubic extension region because sl < slr
1196  EXPECT_NEAR(0.0,
1197  PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1198  0.1, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 1.0, -0.75),
1199  1.0E-12); // cubic extension region because sl < slr, but k_rg_max = 1
1200  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1201  0.3, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1202  (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1203  0.3 + eps, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75) -
1204  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1205  0.3, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75)) /
1206  eps,
1207  1.0E-5); // wetting
1208  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1209  0.3, 0.2, 0.0, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1210  (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1211  0.3 + eps, 0.2, 0.0, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75) -
1212  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1213  0.3, 0.2, 0.0, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75)) /
1214  eps,
1215  1.0E-5); // wetting = drying because sgrdel = 0
1216  EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1217  0.55, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1218  (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1219  0.55 + eps, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75) -
1220  PorousFlowVanGenuchten::relativePermeabilityNWHys(
1221  0.55, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75)) /
1222  eps,
1223  1.0E-5); // wetting
1224  EXPECT_NEAR(0.0,
1225  PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1226  0.76, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1227  1.0E-12); // always zero for large saturation
1228 }
PorousFlowVanGenuchten::relativePermeabilityNWHys
Real relativePermeabilityNWHys(Real sl, Real slr, Real sgrdel, Real sgrmax, Real sldel, Real m, Real gamma, Real k_rg_max, Real y0p)
Hysteretic relative permeability for gas.
Definition: PorousFlowVanGenuchten.C:720
eps
const double eps
Definition: PorousFlowVanGenuchtenTest.C:14
PorousFlowVanGenuchten::drelativePermeabilityNWHys
Real drelativePermeabilityNWHys(Real sl, Real slr, Real sgrdel, Real sgrmax, Real sldel, Real m, Real gamma, Real k_rg_max, Real y0p)
Derivative of hysteretic relative permeability for gas with respect to the liquid saturation...
Definition: PorousFlowVanGenuchten.C:764

Variable Documentation

◆ eps

const double eps = 1.0E-8

Definition at line 14 of file PorousFlowVanGenuchtenTest.C.

Referenced by TEST().

◆ high_ext_dPc

const Real high_ext_dPc
Initial value:
=
PorousFlowVanGenuchten::dcapillaryPressureHys(high_ext_S, 0.1, 0.3, 10.0, 1.9)
PorousFlowVanGenuchten::dcapillaryPressureHys
Real dcapillaryPressureHys(Real sl, Real slmin, Real sgrdel, Real alpha, Real n, const LowCapillaryPressureExtension &low_ext=LowCapillaryPressureExtension(), const HighCapillaryPressureExtension &high_ext=HighCapillaryPressureExtension())
Derivative of capillaryPressureHys with respect to sl.
Definition: PorousFlowVanGenuchten.C:238
high_ext_S
const Real high_ext_S
Definition: PorousFlowVanGenuchtenTest.C:39

Definition at line 42 of file PorousFlowVanGenuchtenTest.C.

◆ high_ext_none

const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_none(PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, 0.8, high_ext_Pc, high_ext_dPc)

Referenced by TEST().

◆ high_ext_none2

const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_none2(PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, high_ext_S, high_ext_Pc, high_ext_dPc)

Referenced by TEST().

◆ high_ext_Pc

const Real high_ext_Pc
Initial value:
=
PorousFlowVanGenuchten::capillaryPressureHys(high_ext_S, 0.1, 0.3, 10.0, 1.9)
high_ext_S
const Real high_ext_S
Definition: PorousFlowVanGenuchtenTest.C:39
PorousFlowVanGenuchten::capillaryPressureHys
Real capillaryPressureHys(Real sl, Real slmin, Real sgrdel, Real alpha, Real n, const LowCapillaryPressureExtension &low_ext=LowCapillaryPressureExtension(), const HighCapillaryPressureExtension &high_ext=HighCapillaryPressureExtension())
Hysteretic capillary pressure function (Eqn(1) of Doughty2007) with extensions (page5 and Fig1 of Dou...
Definition: PorousFlowVanGenuchten.C:180

Definition at line 40 of file PorousFlowVanGenuchtenTest.C.

◆ high_ext_power

const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_power(PorousFlowVanGenuchten::HighCapillaryPressureExtension::POWER, high_ext_S, high_ext_Pc, high_ext_dPc)

Referenced by TEST().

◆ high_ext_S

const Real high_ext_S = 0.9 * (1.0 - 0.3)

Definition at line 39 of file PorousFlowVanGenuchtenTest.C.

◆ low_ext_dPc

const Real low_ext_dPc
Initial value:
=
PorousFlowVanGenuchten::dcapillaryPressureHys(low_ext_S, 0.1, 0.3, 10.0, 1.9)
PorousFlowVanGenuchten::dcapillaryPressureHys
Real dcapillaryPressureHys(Real sl, Real slmin, Real sgrdel, Real alpha, Real n, const LowCapillaryPressureExtension &low_ext=LowCapillaryPressureExtension(), const HighCapillaryPressureExtension &high_ext=HighCapillaryPressureExtension())
Derivative of capillaryPressureHys with respect to sl.
Definition: PorousFlowVanGenuchten.C:238
low_ext_S
const Real low_ext_S
Definition: PorousFlowVanGenuchtenTest.C:19

Definition at line 20 of file PorousFlowVanGenuchtenTest.C.

◆ low_ext_exp

const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_exp(PorousFlowVanGenuchten::LowCapillaryPressureExtension::EXPONENTIAL, low_ext_S, low_ext_Pc, low_ext_dPc)

Referenced by TEST().

◆ low_ext_none

const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_none(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, low_ext_S, low_ext_Pc, low_ext_dPc)

Referenced by TEST().

◆ low_ext_Pc

const Real low_ext_Pc = 1.5

Definition at line 18 of file PorousFlowVanGenuchtenTest.C.

◆ low_ext_quad

const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_quad(PorousFlowVanGenuchten::LowCapillaryPressureExtension::QUADRATIC, low_ext_S, low_ext_Pc, low_ext_dPc)

Referenced by TEST().

◆ low_ext_S

const Real low_ext_S = PorousFlowVanGenuchten::saturationHys(low_ext_Pc, 0.1, 0.3, 10.0, 1.9)

Definition at line 19 of file PorousFlowVanGenuchtenTest.C.

◆ no_high_ext

const PorousFlowVanGenuchten::HighCapillaryPressureExtension no_high_ext(PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, 2.0, 0.0, 0.0)

Referenced by TEST().

◆ no_low_ext

const PorousFlowVanGenuchten::LowCapillaryPressureExtension no_low_ext(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, -1.0, 123.0, 0.0)

Referenced by TEST().

Generated on Fri Jun 27 2025 10:08:00 for https://mooseframework.inl.gov by   doxygen 1.8.14