https://mooseframework.inl.gov
Loading...
Searching...
No Matches
PorousFlowVanGenuchtenTest.C
Go to the documentation of this file.
1//* This file is part of the MOOSE framework
2//* https://mooseframework.inl.gov
3//*
4//* All rights reserved, see COPYRIGHT for full restrictions
5//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
6//*
7//* Licensed under LGPL 2.1, please see LICENSE for details
8//* https://www.gnu.org/licenses/lgpl-2.1.html
9
10#include "gtest/gtest.h"
11
12#include "PorousFlowVanGenuchten.h"
13
14const double eps = 1.0E-8;
15
16const PorousFlowVanGenuchten::LowCapillaryPressureExtension
17 no_low_ext(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, -1.0, 123.0, 0.0);
18const Real low_ext_Pc = 1.5;
19const Real low_ext_S = PorousFlowVanGenuchten::saturationHys(low_ext_Pc, 0.1, 0.3, 10.0, 1.9);
20const Real low_ext_dPc =
21 PorousFlowVanGenuchten::dcapillaryPressureHys(low_ext_S, 0.1, 0.3, 10.0, 1.9);
22const PorousFlowVanGenuchten::LowCapillaryPressureExtension
23 low_ext_exp(PorousFlowVanGenuchten::LowCapillaryPressureExtension::EXPONENTIAL,
24 low_ext_S,
25 low_ext_Pc,
26 low_ext_dPc);
27const PorousFlowVanGenuchten::LowCapillaryPressureExtension
28 low_ext_quad(PorousFlowVanGenuchten::LowCapillaryPressureExtension::QUADRATIC,
29 low_ext_S,
30 low_ext_Pc,
31 low_ext_dPc);
32const PorousFlowVanGenuchten::LowCapillaryPressureExtension
33 low_ext_none(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE,
34 low_ext_S,
35 low_ext_Pc,
36 low_ext_dPc);
37const PorousFlowVanGenuchten::HighCapillaryPressureExtension
38 no_high_ext(PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, 2.0, 0.0, 0.0);
39const Real high_ext_S = 0.9 * (1.0 - 0.3); // here 0.9 is the so-called ratio
40const Real high_ext_Pc =
41 PorousFlowVanGenuchten::capillaryPressureHys(high_ext_S, 0.1, 0.3, 10.0, 1.9);
42const Real high_ext_dPc =
43 PorousFlowVanGenuchten::dcapillaryPressureHys(high_ext_S, 0.1, 0.3, 10.0, 1.9);
44const PorousFlowVanGenuchten::HighCapillaryPressureExtension
45 high_ext_power(PorousFlowVanGenuchten::HighCapillaryPressureExtension::POWER,
46 high_ext_S,
47 high_ext_Pc,
48 high_ext_dPc);
49const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_none(
50 PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, 0.8, high_ext_Pc, high_ext_dPc);
51const PorousFlowVanGenuchten::HighCapillaryPressureExtension
52 high_ext_none2(PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE,
53 high_ext_S,
54 high_ext_Pc,
55 high_ext_dPc);
56
57TEST(PorousFlowVanGenuchtenTest, sat)
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}
67
68TEST(PorousFlowVanGenuchtenTest, dsat)
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}
85
86TEST(PorousFlowVanGenuchtenTest, d2sat)
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}
103
104TEST(PorousFlowVanGenuchtenTest, cap)
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}
113
114TEST(PorousFlowVanGenuchtenTest, dcap)
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}
124
125TEST(PorousFlowVanGenuchtenTest, d2cap)
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}
135
136TEST(PorousFlowVanGenuchtenTest, relperm)
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}
142
143TEST(PorousFlowVanGenuchtenTest, drelperm)
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}
157
158TEST(PorousFlowVanGenuchtenTest, d2relperm)
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}
172
173TEST(PorousFlowVanGenuchtenTest, adrelperm)
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}
186
187TEST(PorousFlowVanGenuchtenTest, relpermNW)
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}
194
195TEST(PorousFlowVanGenuchtenTest, drelpermNW)
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}
209
210TEST(PorousFlowVanGenuchtenTest, d2relpermNW)
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}
224
225TEST(PorousFlowVanGenuchtenTest, adrelpermnw)
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}
238
240TEST(PorousFlowVanGenuchtenTest, caphys)
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}
338
340TEST(PorousFlowVanGenuchtenTest, sathys)
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}
425
427TEST(PorousFlowVanGenuchtenTest, dcaphys)
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}
522
524TEST(PorousFlowVanGenuchtenTest, d2caphys)
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}
644
646TEST(PorousFlowVanGenuchtenTest, dsathys)
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}
754
756TEST(PorousFlowVanGenuchtenTest, d2sathys)
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}
864
866TEST(PorousFlowVanGenuchtenTest, relativePermeabilityHys)
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}
936
938TEST(PorousFlowVanGenuchtenTest, drelativePermeabilityHys)
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}
1040
1042TEST(PorousFlowVanGenuchtenTest, relativePermeabilityNWHys)
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}
1116
1118TEST(PorousFlowVanGenuchtenTest, drelativePermeabilityNWHys)
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}
ADReal
DualNumber< Real, DNDerivativeType, true > ADReal
TEST
TEST(PorousFlowVanGenuchtenTest, sat)
Definition PorousFlowVanGenuchtenTest.C:57
low_ext_Pc
const Real low_ext_Pc
Definition PorousFlowVanGenuchtenTest.C:18
low_ext_exp
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_exp(PorousFlowVanGenuchten::LowCapillaryPressureExtension::EXPONENTIAL, low_ext_S, low_ext_Pc, low_ext_dPc)
high_ext_dPc
const Real high_ext_dPc
Definition PorousFlowVanGenuchtenTest.C:42
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)
eps
const double eps
Definition PorousFlowVanGenuchtenTest.C:14
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)
high_ext_none2
const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_none2(PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, high_ext_S, high_ext_Pc, high_ext_dPc)
no_low_ext
const PorousFlowVanGenuchten::LowCapillaryPressureExtension no_low_ext(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, -1.0, 123.0, 0.0)
high_ext_S
const Real high_ext_S
Definition PorousFlowVanGenuchtenTest.C:39
high_ext_Pc
const Real high_ext_Pc
Definition PorousFlowVanGenuchtenTest.C:40
low_ext_none
const PorousFlowVanGenuchten::LowCapillaryPressureExtension low_ext_none(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, low_ext_S, low_ext_Pc, low_ext_dPc)
low_ext_S
const Real low_ext_S
Definition PorousFlowVanGenuchtenTest.C:19
low_ext_dPc
const Real low_ext_dPc
Definition PorousFlowVanGenuchtenTest.C:20
Moose::derivInsert
void derivInsert(SemiDynamicSparseNumberArray< Real, libMesh::dof_id_type, NWrapper< N > > &derivs, libMesh::dof_id_type index, Real value)
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::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::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::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::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
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
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
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
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::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
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::capillaryPressure
Real capillaryPressure(Real seff, Real alpha, Real m, Real pc_max)
Capillary pressure as a function of effective saturation.
Definition PorousFlowVanGenuchten.C:63
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::d2CapillaryPressure
Real d2CapillaryPressure(Real seff, Real alpha, Real m, Real pc_max)
Second derivative of capillary pressure wrt effective saturation.
Definition PorousFlowVanGenuchten.C:93
PorousFlowVanGenuchten::effectiveSaturation
Real effectiveSaturation(Real p, Real alpha, Real m)
Effective saturation as a function of porepressure.
Definition PorousFlowVanGenuchten.C:16
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
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
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::d2RelativePermeability
Real d2RelativePermeability(Real seff, Real m)
Second derivative of relative permeability with respect to effective saturation.
Definition PorousFlowVanGenuchten.C:129
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
PorousFlowVanGenuchten::HighCapillaryPressureExtension
Parameters associated with the extension of the hysteretic wetting capillary pressure function to hig...
Definition PorousFlowVanGenuchten.h:216
PorousFlowVanGenuchten::LowCapillaryPressureExtension
Parameters associated with the extension of the hysteretic capillary pressure function to low saturat...
Definition PorousFlowVanGenuchten.h:181
Generated on Fri Aug 21 2026 13:41:31 for https://mooseframework.inl.gov by doxygen 1.9.8