https://mooseframework.inl.gov
Loading...
Searching...
No Matches
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.
 
 TEST (PorousFlowVanGenuchtenTest, sathys)
 Test the hysteretic saturation.
 
 TEST (PorousFlowVanGenuchtenTest, dcaphys)
 Test the first derivative of the hysteretic capillary pressure.
 
 TEST (PorousFlowVanGenuchtenTest, d2caphys)
 Test the second derivative of the hysteretic capillary pressure.
 
 TEST (PorousFlowVanGenuchtenTest, dsathys)
 Test the derivative of the hysteretic saturation.
 
 TEST (PorousFlowVanGenuchtenTest, d2sathys)
 Test the second derivative of the hysteretic saturation.
 
 TEST (PorousFlowVanGenuchtenTest, relativePermeabilityHys)
 Test relativePermeabilityHys.
 
 TEST (PorousFlowVanGenuchtenTest, drelativePermeabilityHys)
 Test drelativePermeabilityHys.
 
 TEST (PorousFlowVanGenuchtenTest, relativePermeabilityNWHys)
 Test relativePermeabilityNWHys.
 
 TEST (PorousFlowVanGenuchtenTest, drelativePermeabilityNWHys)
 Test drelativePermeabilityNWHys.
 

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  ,
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
Moose::derivInsert
void derivInsert(SemiDynamicSparseNumberArray< Real, libMesh::dof_id_type, NWrapper< N > > &derivs, libMesh::dof_id_type index, Real value)
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

◆ TEST() [2/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:140
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() [3/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() [4/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)
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_quad
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_quad(PorousFlowVanGenuchten::LowCapillaryPressureExtension::QUADRATIC, low_ext_S, low_ext_Pc, low_ext_dPc)
no_low_ext
const PorousFlowVanGenuchten::LowCapillaryPressureExtension no_low_ext(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, -1.0, 123.0, 0.0)
low_ext_none
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_none(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, low_ext_S, low_ext_Pc, low_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

◆ TEST() [5/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::dCapillaryPressure
Real dCapillaryPressure(Real seff, Real alpha, Real m, Real pc_max)
Derivative of capillary pressure wrt effective saturation.
Definition PorousFlowVanGenuchten.C:77
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
libMesh::Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ TEST() [6/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}
high_ext_none2
const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_none2(PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, high_ext_S, high_ext_Pc, high_ext_dPc)
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
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

◆ TEST() [7/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}
PorousFlowVanGenuchten::d2RelativePermeability
Real d2RelativePermeability(Real seff, Real m)
Second derivative of relative permeability with respect to effective saturation.
Definition PorousFlowVanGenuchten.C:129

◆ TEST() [8/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}
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

◆ TEST() [9/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}
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() [10/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}
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
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

◆ TEST() [11/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}

◆ TEST() [12/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}

◆ TEST() [13/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
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() [14/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::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
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() [15/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}

◆ TEST() [16/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}

◆ TEST() [17/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}
PorousFlowVanGenuchten::effectiveSaturation
Real effectiveSaturation(Real p, Real alpha, Real m)
Effective saturation as a function of porepressure.
Definition PorousFlowVanGenuchten.C:16

◆ TEST() [18/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}
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

◆ TEST() [19/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}

◆ TEST() [20/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}

◆ TEST() [21/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}

◆ TEST() [22/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}

◆ TEST() [23/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}

◆ TEST() [24/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}

Variable Documentation

◆ eps

const double eps = 1.0E-8

Definition at line 14 of file PorousFlowVanGenuchtenTest.C.

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

◆ high_ext_dPc

const Real high_ext_dPc
Initial value:
=
PorousFlowVanGenuchten::dcapillaryPressureHys(high_ext_S, 0.1, 0.3, 10.0, 1.9)
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) ( PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE  ,
0.  8,
high_ext_Pc  ,
high_ext_dPc   
)

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

◆ high_ext_none2

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

Referenced by TEST(), and TEST().

◆ high_ext_Pc

const Real high_ext_Pc
Initial value:
=
PorousFlowVanGenuchten::capillaryPressureHys(high_ext_S, 0.1, 0.3, 10.0, 1.9)

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) ( PorousFlowVanGenuchten::HighCapillaryPressureExtension::POWER  ,
high_ext_S  ,
high_ext_Pc  ,
high_ext_dPc   
)

Referenced by TEST(), TEST(), TEST(), TEST(), TEST(), and 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)
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) ( PorousFlowVanGenuchten::LowCapillaryPressureExtension::EXPONENTIAL  ,
low_ext_S  ,
low_ext_Pc  ,
low_ext_dPc   
)

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

◆ low_ext_none

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

Referenced by TEST(), TEST(), TEST(), TEST(), TEST(), and 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) ( PorousFlowVanGenuchten::LowCapillaryPressureExtension::QUADRATIC  ,
low_ext_S  ,
low_ext_Pc  ,
low_ext_dPc   
)

Referenced by TEST(), TEST(), TEST(), TEST(), TEST(), and 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) ( PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE  ,
2.  0,
0.  0,
0.  0 
)

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

◆ no_low_ext

const PorousFlowVanGenuchten::LowCapillaryPressureExtension no_low_ext(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, -1.0, 123.0, 0.0) ( PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE  ,
-1.  0,
123.  0,
0.  0 
)

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

Generated on Fri Aug 21 2026 13:42:05 for https://mooseframework.inl.gov by doxygen 1.9.8