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 "DataIO.h"
13#include "PorousFlowVanGenuchten.h"
14
15const double eps = 1.0E-8;
16
17const PorousFlowVanGenuchten::LowCapillaryPressureExtension
18 no_low_ext(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE, -1.0, 123.0, 0.0);
19const Real low_ext_Pc = 1.5;
20const Real low_ext_S = PorousFlowVanGenuchten::saturationHys(low_ext_Pc, 0.1, 0.3, 10.0, 1.9);
21const Real low_ext_dPc =
22 PorousFlowVanGenuchten::dcapillaryPressureHys(low_ext_S, 0.1, 0.3, 10.0, 1.9);
23const PorousFlowVanGenuchten::LowCapillaryPressureExtension
24 low_ext_exp(PorousFlowVanGenuchten::LowCapillaryPressureExtension::EXPONENTIAL,
25 low_ext_S,
26 low_ext_Pc,
27 low_ext_dPc);
28const PorousFlowVanGenuchten::LowCapillaryPressureExtension
29 low_ext_quad(PorousFlowVanGenuchten::LowCapillaryPressureExtension::QUADRATIC,
30 low_ext_S,
31 low_ext_Pc,
32 low_ext_dPc);
33const PorousFlowVanGenuchten::LowCapillaryPressureExtension
34 low_ext_none(PorousFlowVanGenuchten::LowCapillaryPressureExtension::NONE,
35 low_ext_S,
36 low_ext_Pc,
37 low_ext_dPc);
38const PorousFlowVanGenuchten::HighCapillaryPressureExtension
39 no_high_ext(PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, 2.0, 0.0, 0.0);
40const Real high_ext_S = 0.9 * (1.0 - 0.3); // here 0.9 is the so-called ratio
41const Real high_ext_Pc =
42 PorousFlowVanGenuchten::capillaryPressureHys(high_ext_S, 0.1, 0.3, 10.0, 1.9);
43const Real high_ext_dPc =
44 PorousFlowVanGenuchten::dcapillaryPressureHys(high_ext_S, 0.1, 0.3, 10.0, 1.9);
45const PorousFlowVanGenuchten::HighCapillaryPressureExtension
46 high_ext_power(PorousFlowVanGenuchten::HighCapillaryPressureExtension::POWER,
47 high_ext_S,
48 high_ext_Pc,
49 high_ext_dPc);
50const PorousFlowVanGenuchten::HighCapillaryPressureExtension high_ext_none(
51 PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE, 0.8, high_ext_Pc, high_ext_dPc);
52const PorousFlowVanGenuchten::HighCapillaryPressureExtension
53 high_ext_none2(PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE,
54 high_ext_S,
55 high_ext_Pc,
56 high_ext_dPc);
57
58TEST(PorousFlowVanGenuchtenTest, extensionsDataIO)
59{
60 std::array<PorousFlowVanGenuchten::LowCapillaryPressureExtension, 1> stored_low = {low_ext_exp};
61 std::array<PorousFlowVanGenuchten::HighCapillaryPressureExtension, 1> stored_high = {
62 high_ext_power};
63 std::stringstream stream;
64 dataStore(stream, stored_low, nullptr);
65 dataStore(stream, stored_high, nullptr);
66
67 stream.seekg(0, std::ios::beg);
68 std::array<PorousFlowVanGenuchten::LowCapillaryPressureExtension, 1> loaded_low;
69 std::array<PorousFlowVanGenuchten::HighCapillaryPressureExtension, 1> loaded_high;
70 dataLoad(stream, loaded_low, nullptr);
71 dataLoad(stream, loaded_high, nullptr);
72
73 EXPECT_EQ(stored_low[0].strategy, loaded_low[0].strategy);
74 EXPECT_EQ(stored_low[0].S, loaded_low[0].S);
75 EXPECT_EQ(stored_low[0].Pc, loaded_low[0].Pc);
76 EXPECT_EQ(stored_low[0].dPc, loaded_low[0].dPc);
77 EXPECT_EQ(stored_high[0].strategy, loaded_high[0].strategy);
78 EXPECT_EQ(stored_high[0].S, loaded_high[0].S);
79 EXPECT_EQ(stored_high[0].Pc, loaded_high[0].Pc);
80 EXPECT_EQ(stored_high[0].dPc, loaded_high[0].dPc);
81}
82
83TEST(PorousFlowVanGenuchtenTest, sat)
84{
85 EXPECT_NEAR(1.0, PorousFlowVanGenuchten::effectiveSaturation(1.0E30, 0.7, 0.5), 1.0E-5);
86 EXPECT_NEAR(1.0, PorousFlowVanGenuchten::effectiveSaturation(1.0, 0.7, 0.5), 1.0E-5);
87 EXPECT_NEAR(1.0, PorousFlowVanGenuchten::effectiveSaturation(0.0, 0.7, 0.5), 1.0E-5);
88 EXPECT_NEAR(1.0, PorousFlowVanGenuchten::effectiveSaturation(1.0E-10, 0.7, 0.5), 1.0E-5);
89 EXPECT_NEAR(
90 0.486841442435055, PorousFlowVanGenuchten::effectiveSaturation(-2.0, 0.7, 0.6), 1.0E-5);
91 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::effectiveSaturation(-1.0E30, 0.7, 0.5), 1.0E-5);
92}
93
94TEST(PorousFlowVanGenuchtenTest, dsat)
95{
96 Real fd;
97 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dEffectiveSaturation(1.0E30, 0.7, 0.5), 1.0E-5);
98 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dEffectiveSaturation(1.0, 0.7, 0.5), 1.0E-5);
99 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dEffectiveSaturation(0.0, 0.7, 0.5), 1.0E-5);
100 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dEffectiveSaturation(1.0E-10, 0.7, 0.5), 1.0E-5);
101 fd = (PorousFlowVanGenuchten::effectiveSaturation(-2.0 + eps, 0.7, 0.6) -
102 PorousFlowVanGenuchten::effectiveSaturation(-2.0, 0.7, 0.6)) /
103 eps;
104 EXPECT_NEAR(fd, PorousFlowVanGenuchten::dEffectiveSaturation(-2.0, 0.7, 0.6), 1.0E-5);
105 fd = (PorousFlowVanGenuchten::effectiveSaturation(-1.1 + eps, 0.9, 0.66) -
106 PorousFlowVanGenuchten::effectiveSaturation(-1.1, 0.9, 0.66)) /
107 eps;
108 EXPECT_NEAR(fd, PorousFlowVanGenuchten::dEffectiveSaturation(-1.1, 0.9, 0.66), 1.0E-5);
109 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dEffectiveSaturation(-1.0E30, 0.7, 0.5), 1.0E-5);
110}
111
112TEST(PorousFlowVanGenuchtenTest, d2sat)
113{
114 Real fd;
115 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2EffectiveSaturation(1.0E30, 0.7, 0.5), 1.0E-5);
116 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2EffectiveSaturation(1.0, 0.7, 0.5), 1.0E-5);
117 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2EffectiveSaturation(0.0, 0.7, 0.5), 1.0E-5);
118 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2EffectiveSaturation(1.0E-10, 0.7, 0.5), 1.0E-5);
119 fd = (PorousFlowVanGenuchten::dEffectiveSaturation(-2.0 + eps, 0.7, 0.6) -
120 PorousFlowVanGenuchten::dEffectiveSaturation(-2.0, 0.7, 0.6)) /
121 eps;
122 EXPECT_NEAR(fd, PorousFlowVanGenuchten::d2EffectiveSaturation(-2.0, 0.7, 0.6), 1.0E-5);
123 fd = (PorousFlowVanGenuchten::dEffectiveSaturation(-1.1 + eps, 2.3, 0.67) -
124 PorousFlowVanGenuchten::dEffectiveSaturation(-1.1, 2.3, 0.67)) /
125 eps;
126 EXPECT_NEAR(fd, PorousFlowVanGenuchten::d2EffectiveSaturation(-1.1, 2.3, 0.67), 1.0E-5);
127 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2EffectiveSaturation(-1.0E30, 0.7, 0.5), 1.0E-5);
128}
129
130TEST(PorousFlowVanGenuchtenTest, cap)
131{
132 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::capillaryPressure(1.1, 1.0, 0.55, 1.0E30), 1.0E-5);
133 EXPECT_NEAR(4.06172392297447,
134 PorousFlowVanGenuchten::capillaryPressure(0.3, 0.625, 0.55, 1.0E30),
135 1.0E-5);
136 EXPECT_NEAR(2.9, PorousFlowVanGenuchten::capillaryPressure(0.001, 0.625, 0.55, 2.9), 1.0E-5);
137 EXPECT_NEAR(1000.0, PorousFlowVanGenuchten::capillaryPressure(0.0, 0.625, 0.55, 1000.0), 1.0E-5);
138}
139
140TEST(PorousFlowVanGenuchtenTest, dcap)
141{
142 Real fd;
143 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dCapillaryPressure(0.0, 1.0, 0.55, 1.0E30), 1.0E-5);
144 fd = (PorousFlowVanGenuchten::capillaryPressure(0.3 + eps, 0.625, 0.55, 1.0E30) -
145 PorousFlowVanGenuchten::capillaryPressure(0.3, 0.625, 0.55, 1.0E30)) /
146 eps;
147 EXPECT_NEAR(fd, PorousFlowVanGenuchten::dCapillaryPressure(0.3, 0.625, 0.55, 1.0E30), 1.0E-5);
148 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dCapillaryPressure(1.0, 1.0, 0.55, 1.0E30), 1.0E-5);
149}
150
151TEST(PorousFlowVanGenuchtenTest, d2cap)
152{
153 Real fd;
154 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2CapillaryPressure(0.0, 1.0, 0.55, 1.0E30), 1.0E-5);
155 fd = (PorousFlowVanGenuchten::dCapillaryPressure(0.3 + eps, 0.625, 0.55, 1.0E30) -
156 PorousFlowVanGenuchten::dCapillaryPressure(0.3, 0.625, 0.55, 1.0E30)) /
157 eps;
158 EXPECT_NEAR(fd, PorousFlowVanGenuchten::d2CapillaryPressure(0.3, 0.625, 0.55, 1.0E30), 1.0E-5);
159 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2CapillaryPressure(1.0, 0.625, 0.55, 1.0E30), 1.0E-5);
160}
161
162TEST(PorousFlowVanGenuchtenTest, relperm)
163{
164 EXPECT_NEAR(1.0, PorousFlowVanGenuchten::relativePermeability(1.0E30, 0.7), 1.0E-5);
165 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::relativePermeability(-1.0, 0.7), 1.0E-5);
166 EXPECT_NEAR(0.0091160727, PorousFlowVanGenuchten::relativePermeability(0.3, 0.7), 1.0E-5);
167}
168
169TEST(PorousFlowVanGenuchtenTest, drelperm)
170{
171 Real fd;
172 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dRelativePermeability(1.0E30, 0.7), 1.0E-5);
173 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dRelativePermeability(-1.0, 0.7), 1.0E-5);
174 fd = (PorousFlowVanGenuchten::relativePermeability(0.3 + eps, 0.7) -
175 PorousFlowVanGenuchten::relativePermeability(0.3, 0.7)) /
176 eps;
177 EXPECT_NEAR(fd, PorousFlowVanGenuchten::dRelativePermeability(0.3, 0.7), 1.0E-5);
178 fd = (PorousFlowVanGenuchten::relativePermeability(0.8 + eps, 0.65) -
179 PorousFlowVanGenuchten::relativePermeability(0.8, 0.65)) /
180 eps;
181 EXPECT_NEAR(fd, PorousFlowVanGenuchten::dRelativePermeability(0.8, 0.65), 1.0E-5);
182}
183
184TEST(PorousFlowVanGenuchtenTest, d2relperm)
185{
186 Real fd;
187 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2RelativePermeability(1.0E30, 0.7), 1.0E-5);
188 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2RelativePermeability(-1.0, 0.7), 1.0E-5);
189 fd = (PorousFlowVanGenuchten::dRelativePermeability(0.3 + eps, 0.7) -
190 PorousFlowVanGenuchten::dRelativePermeability(0.3, 0.7)) /
191 eps;
192 EXPECT_NEAR(fd, PorousFlowVanGenuchten::d2RelativePermeability(0.3, 0.7), 1.0E-5);
193 fd = (PorousFlowVanGenuchten::dRelativePermeability(0.8 + eps, 0.65) -
194 PorousFlowVanGenuchten::dRelativePermeability(0.8, 0.65)) /
195 eps;
196 EXPECT_NEAR(fd, PorousFlowVanGenuchten::d2RelativePermeability(0.8, 0.65), 1.0E-5);
197}
198
199TEST(PorousFlowVanGenuchtenTest, adrelperm)
200{
201 ADReal sat = 0.3;
202 Moose::derivInsert(sat.derivatives(), 0, 1.0);
203
204 const auto adrelperm = PorousFlowVanGenuchten::relativePermeability(sat, 0.7);
205
206 const auto relperm = PorousFlowVanGenuchten::relativePermeability(sat.value(), 0.7);
207 const auto drelperm = PorousFlowVanGenuchten::dRelativePermeability(sat.value(), 0.7);
208
209 EXPECT_NEAR(adrelperm.value(), relperm, 1.0E-5);
210 EXPECT_NEAR(adrelperm.derivatives()[0], drelperm, 1.0E-5);
211}
212
213TEST(PorousFlowVanGenuchtenTest, relpermNW)
214{
215 EXPECT_NEAR(1.0, PorousFlowVanGenuchten::relativePermeabilityNW(1.0E30, 0.7), 1.0E-5);
216 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::relativePermeabilityNW(-1.0, 0.7), 1.0E-5);
217 EXPECT_NEAR(0.664138097, PorousFlowVanGenuchten::relativePermeabilityNW(0.5, 0.3), 1.0E-5);
218 EXPECT_NEAR(0.559879012, PorousFlowVanGenuchten::relativePermeabilityNW(0.7, 0.8), 1.0E-5);
219}
220
221TEST(PorousFlowVanGenuchtenTest, drelpermNW)
222{
223 Real fd;
224 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dRelativePermeabilityNW(1.0E30, 0.7), 1.0E-5);
225 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dRelativePermeabilityNW(-1.0, 0.7), 1.0E-5);
226 fd = (PorousFlowVanGenuchten::relativePermeabilityNW(0.3 + eps, 0.7) -
227 PorousFlowVanGenuchten::relativePermeabilityNW(0.3, 0.7)) /
228 eps;
229 EXPECT_NEAR(fd, PorousFlowVanGenuchten::dRelativePermeabilityNW(0.3, 0.7), 1.0E-5);
230 fd = (PorousFlowVanGenuchten::relativePermeabilityNW(0.8 + eps, 0.65) -
231 PorousFlowVanGenuchten::relativePermeabilityNW(0.8, 0.65)) /
232 eps;
233 EXPECT_NEAR(fd, PorousFlowVanGenuchten::dRelativePermeabilityNW(0.8, 0.65), 1.0E-5);
234}
235
236TEST(PorousFlowVanGenuchtenTest, d2relpermNW)
237{
238 Real fd;
239 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2RelativePermeabilityNW(1.0E30, 0.7), 1.0E-5);
240 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2RelativePermeabilityNW(-1.0, 0.7), 1.0E-5);
241 fd = (PorousFlowVanGenuchten::dRelativePermeabilityNW(0.3 + eps, 0.7) -
242 PorousFlowVanGenuchten::dRelativePermeabilityNW(0.3, 0.7)) /
243 eps;
244 EXPECT_NEAR(fd, PorousFlowVanGenuchten::d2RelativePermeabilityNW(0.3, 0.7), 1.0E-5);
245 fd = (PorousFlowVanGenuchten::dRelativePermeabilityNW(0.8 + eps, 0.65) -
246 PorousFlowVanGenuchten::dRelativePermeabilityNW(0.8, 0.65)) /
247 eps;
248 EXPECT_NEAR(fd, PorousFlowVanGenuchten::d2RelativePermeabilityNW(0.8, 0.65), 1.0E-5);
249}
250
251TEST(PorousFlowVanGenuchtenTest, adrelpermnw)
252{
253 ADReal sat = 0.3;
254 Moose::derivInsert(sat.derivatives(), 0, 1.0);
255
256 const auto adrelperm = PorousFlowVanGenuchten::relativePermeabilityNW(sat, 0.7);
257
258 const auto relperm = PorousFlowVanGenuchten::relativePermeabilityNW(sat.value(), 0.7);
259 const auto drelperm = PorousFlowVanGenuchten::dRelativePermeabilityNW(sat.value(), 0.7);
260
261 EXPECT_NEAR(adrelperm.value(), relperm, 1.0E-5);
262 EXPECT_NEAR(adrelperm.derivatives()[0], drelperm, 1.0E-5);
263}
264
266TEST(PorousFlowVanGenuchtenTest, caphys)
267{
268 // first define some extensions that do not actually induce any extension, so the non-extended Pc
269 // can be checked
270 const std::vector<Real> no_ext_sats{1.1, 0.75, 0.2, 0.4, 0.6};
271 const std::vector<Real> expected_pcs{
272 0.0, 0.0, 0.7233512030263158, 0.18806139488299345, 0.06716785873826017};
273 for (unsigned i = 0; i < no_ext_sats.size(); ++i)
274 {
275 EXPECT_NEAR(expected_pcs[i],
276 PorousFlowVanGenuchten::capillaryPressureHys(
277 no_ext_sats[i], 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
278 1.0E-5);
279 EXPECT_NEAR(expected_pcs[i],
280 PorousFlowVanGenuchten::capillaryPressureHys(no_ext_sats[i], 0.1, 0.3, 10.0, 1.9),
281 1.0E-5);
282 }
283 EXPECT_NEAR(123.0,
284 PorousFlowVanGenuchten::capillaryPressureHys(
285 -1.0, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
286 1.0E-5);
287 EXPECT_NEAR(123.0,
288 PorousFlowVanGenuchten::capillaryPressureHys(
289 0.05, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
290 1.0E-5);
291 EXPECT_NEAR(std::numeric_limits<Real>::max(),
292 PorousFlowVanGenuchten::capillaryPressureHys(-1.0, 0.1, 0.3, 10.0, 1.9),
293 1.0E-5);
294 EXPECT_NEAR(std::numeric_limits<Real>::max(),
295 PorousFlowVanGenuchten::capillaryPressureHys(0.05, 0.1, 0.3, 10.0, 1.9),
296 1.0E-5);
297
298 // check the lower extension values
299 EXPECT_NEAR(0.15229665668299805, low_ext_exp.S, 1.0E-5);
300 EXPECT_NEAR(-32.05516428734052, low_ext_exp.dPc, 1.0E-5);
301
302 // check the high extension
303 EXPECT_NEAR(0.05300654102157442, high_ext_power.Pc, 1.0E-5);
304 EXPECT_NEAR(-0.48230528124844047, high_ext_power.dPc, 1.0E-5);
305
306 // using low_ext_exp and high_ext_power
307 const std::vector<Real> sats{0.01, 0.1, 0.15, 0.2, 0.5, 0.63, 0.8, 0.99, 1.0};
308 std::vector<Real> pc;
309 pc = {31.385947046636815,
310 4.5861935551774735,
311 1.5754562501536364,
312 0.7233512030263158,
313 0.11727884570711045,
314 0.05300654102157442,
315 0.006681337544884095,
316 2.7847615834816514e-07,
317 0.0};
318 for (unsigned i = 0; i < sats.size(); ++i)
319 EXPECT_NEAR(pc[i],
320 PorousFlowVanGenuchten::capillaryPressureHys(
321 sats[i], 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power),
322 1.0E-5);
323
324 // using low_ext_quad and high_ext_power
325 pc = {3.9304232526771696,
326 2.8885549236001244,
327 1.5730646091089062,
328 0.7233512030263158,
329 0.11727884570711045,
330 0.05300654102157442,
331 0.006681337544884095,
332 2.7847615834816514e-07,
333 0.0};
334 for (unsigned i = 0; i < sats.size(); ++i)
335 EXPECT_NEAR(pc[i],
336 PorousFlowVanGenuchten::capillaryPressureHys(
337 sats[i], 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power),
338 1.0E-5);
339
340 // using low_ext_none and high_ext_power
341 pc = {1.5,
342 1.5,
343 1.5,
344 0.7233512030263158,
345 0.11727884570711045,
346 0.05300654102157442,
347 0.006681337544884095,
348 2.7847615834816514e-07,
349 0.0};
350 for (unsigned i = 0; i < sats.size(); ++i)
351 EXPECT_NEAR(pc[i],
352 PorousFlowVanGenuchten::capillaryPressureHys(
353 sats[i], 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power),
354 1.0E-5);
355
356 // using low_ext_none and high_ext_none
357 pc = {1.5, 1.5, 1.5, 0.7233512030263158, 0.11727884570711045, 0.05300654102157442, 0.0, 0.0, 0.0};
358 for (unsigned i = 0; i < sats.size(); ++i)
359 EXPECT_NEAR(pc[i],
360 PorousFlowVanGenuchten::capillaryPressureHys(
361 sats[i], 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none),
362 1.0E-5);
363}
364
366TEST(PorousFlowVanGenuchtenTest, sathys)
367{
368 // first define some extensions that do not actually induce any extension, so the non-extended Pc
369 // can be checked
370 const std::vector<Real> expected_sats{1.0, 0.2, 0.4, 0.6};
371 const std::vector<Real> no_ext_pcs{
372 -1.1, 0.7233512030263158, 0.18806139488299345, 0.06716785873826017};
373 for (unsigned i = 0; i < no_ext_pcs.size(); ++i)
374 {
375 EXPECT_NEAR(expected_sats[i],
376 PorousFlowVanGenuchten::saturationHys(
377 no_ext_pcs[i], 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
378 1.0E-5);
379 EXPECT_NEAR(expected_sats[i],
380 PorousFlowVanGenuchten::saturationHys(no_ext_pcs[i], 0.1, 0.3, 10.0, 1.9),
381 1.0E-5);
382 }
383 EXPECT_NEAR(
384 -1.0,
385 PorousFlowVanGenuchten::saturationHys(1234.0, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
386 1.0E-5);
387 EXPECT_NEAR(
388 0.0,
389 PorousFlowVanGenuchten::saturationHys(std::numeric_limits<Real>::max(), 0.1, 0.3, 10.0, 1.9),
390 1.0E-5);
391
392 // now with extensions
393 std::vector<Real> sats{0.01, 0.1, 0.15, 0.2, 0.5, 0.63, 0.8, 0.99, 1.0};
394 std::vector<Real> pc;
395 pc = {31.385947046636815,
396 4.5861935551774735,
397 1.5754562501536364,
398 0.7233512030263158,
399 0.11727884570711045,
400 0.05300654102157442,
401 0.006681337544884095,
402 2.7847615834816514e-07,
403 0.0};
404 for (unsigned i = 0; i < sats.size(); ++i)
405 EXPECT_NEAR(sats[i],
406 PorousFlowVanGenuchten::saturationHys(
407 pc[i], 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power),
408 1.0E-5);
409
410 // different lower extension
411 pc = {3.9304232526771696,
412 2.8885549236001244,
413 1.5730646091089062,
414 0.7233512030263158,
415 0.11727884570711045,
416 0.05300654102157442,
417 0.006681337544884095,
418 2.7847615834816514e-07,
419 0.0};
420 for (unsigned i = 0; i < sats.size(); ++i)
421 EXPECT_NEAR(sats[i],
422 PorousFlowVanGenuchten::saturationHys(
423 pc[i], 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power),
424 1.0E-5);
425
426 // different low extension
427 sats = {low_ext_none.S, low_ext_none.S, 0.2, 0.5, 0.63, 0.8, 0.99, 1.0};
428 pc = {2.5,
429 1.5,
430 0.7233512030263158,
431 0.11727884570711045,
432 0.05300654102157442,
433 0.006681337544884095,
434 2.7847615834816514e-07,
435 0.0};
436 for (unsigned i = 0; i < sats.size(); ++i)
437 EXPECT_NEAR(sats[i],
438 PorousFlowVanGenuchten::saturationHys(
439 pc[i], 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power),
440 1.0E-5);
441
442 // different high extension
443 pc = {2.5, 1.5, 0.7233512030263158, 0.11727884570711045, 0.05300654102157442, 1E-10};
444 sats = {low_ext_none.S, low_ext_none.S, 0.2, 0.5, 0.63, 0.8};
445 for (unsigned i = 0; i < sats.size(); ++i)
446 EXPECT_NEAR(sats[i],
447 PorousFlowVanGenuchten::saturationHys(
448 pc[i], 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none),
449 1.0E-5);
450}
451
453TEST(PorousFlowVanGenuchtenTest, dcaphys)
454{
455 const Real eps = 1.0E-9;
456
457 // first define some extensions that do not actually induce any extension, so the non-extended
458 // Pc can be checked
459 std::vector<Real> sats{1.1, 0.75, -1.0, 0.05};
460 for (const auto & sat : sats)
461 {
462 EXPECT_NEAR(0.0,
463 PorousFlowVanGenuchten::dcapillaryPressureHys(
464 sat, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
465 1.0E-5);
466 EXPECT_NEAR(
467 0.0, PorousFlowVanGenuchten::dcapillaryPressureHys(sat, 0.1, 0.3, 10.0, 1.9), 1.0E-5);
468 }
469
470 sats = {0.2, 0.4, 0.6};
471 for (const auto & sat : sats)
472 {
473 const Real fd = 0.5 *
474 (PorousFlowVanGenuchten::capillaryPressureHys(sat + eps, 0.1, 0.3, 10.0, 1.9) -
475 PorousFlowVanGenuchten::capillaryPressureHys(sat - eps, 0.1, 0.3, 10.0, 1.9)) /
476 eps;
477 EXPECT_NEAR(fd,
478 PorousFlowVanGenuchten::dcapillaryPressureHys(
479 sat, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
480 1.0E-5);
481 EXPECT_NEAR(
482 fd, PorousFlowVanGenuchten::dcapillaryPressureHys(sat, 0.1, 0.3, 10.0, 1.9), 1.0E-5);
483 };
484
485 // now with extensions
486 sats = {0.01, 0.1, 0.15, 0.2, 0.5, 0.63, 0.8, 0.99, 1.0};
487
488 for (const auto & sat : sats)
489 {
490 const Real fd = 0.5 *
491 (PorousFlowVanGenuchten::capillaryPressureHys(
492 sat + eps, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power) -
493 PorousFlowVanGenuchten::capillaryPressureHys(
494 sat - eps, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power)) /
495 eps;
496 EXPECT_NEAR(fd,
497 PorousFlowVanGenuchten::dcapillaryPressureHys(
498 sat, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power),
499 1.0E-5);
500 };
501
502 // different low extension
503 for (const auto & sat : sats)
504 {
505 const Real fd = 0.5 *
506 (PorousFlowVanGenuchten::capillaryPressureHys(
507 sat + eps, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power) -
508 PorousFlowVanGenuchten::capillaryPressureHys(
509 sat - eps, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power)) /
510 eps;
511 EXPECT_NEAR(fd,
512 PorousFlowVanGenuchten::dcapillaryPressureHys(
513 sat, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power),
514 1.0E-5);
515 };
516
517 // different lower extension
518 for (const auto & sat : sats)
519 {
520 const Real fd = 0.5 *
521 (PorousFlowVanGenuchten::capillaryPressureHys(
522 sat + eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power) -
523 PorousFlowVanGenuchten::capillaryPressureHys(
524 sat - eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power)) /
525 eps;
526 EXPECT_NEAR(fd,
527 PorousFlowVanGenuchten::dcapillaryPressureHys(
528 sat, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power),
529 1.0E-5);
530 };
531
532 // different upper extension (cannot evaluate derivative at s = 0.63)
533 sats = {0.01, 0.1, 0.15, 0.2, 0.5, 0.8, 0.99, 1.0};
534 for (const auto & sat : sats)
535 {
536 const Real fd = 0.5 *
537 (PorousFlowVanGenuchten::capillaryPressureHys(
538 sat + eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none2) -
539 PorousFlowVanGenuchten::capillaryPressureHys(
540 sat - eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none2)) /
541 eps;
542 EXPECT_NEAR(fd,
543 PorousFlowVanGenuchten::dcapillaryPressureHys(
544 sat, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none2),
545 1.0E-5);
546 };
547}
548
550TEST(PorousFlowVanGenuchtenTest, d2caphys)
551{
552 const Real eps = 1.0E-9;
553
554 // first define some extensions that do not actually induce any extension, so the non-extended
555 // Pc can be checked
556 std::vector<Real> sats{1.1, 0.75, -1.0, 0.05};
557 for (const auto & sat : sats)
558 {
559 EXPECT_NEAR(0.0,
560 PorousFlowVanGenuchten::d2capillaryPressureHys(
561 sat, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
562 1.0E-5);
563 EXPECT_NEAR(
564 0.0, PorousFlowVanGenuchten::d2capillaryPressureHys(sat, 0.1, 0.3, 10.0, 1.9), 1.0E-5);
565 }
566
567 sats = {0.2, 0.4, 0.6};
568 for (const auto & sat : sats)
569 {
570 const Real fd =
571 0.5 *
572 (PorousFlowVanGenuchten::dcapillaryPressureHys(sat + eps, 0.1, 0.3, 10.0, 1.9) -
573 PorousFlowVanGenuchten::dcapillaryPressureHys(sat - eps, 0.1, 0.3, 10.0, 1.9)) /
574 eps;
575 EXPECT_NEAR(fd,
576 PorousFlowVanGenuchten::d2capillaryPressureHys(
577 sat, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
578 1.0E-5);
579 EXPECT_NEAR(
580 fd, PorousFlowVanGenuchten::d2capillaryPressureHys(sat, 0.1, 0.3, 10.0, 1.9), 1.0E-5);
581 };
582
583 // now with extensions
584 sats = {0.01, 0.1, 0.15, 0.2, 0.5, 0.6, 0.7, 0.8, 0.99, 1.0};
585
586 for (const auto & sat : sats)
587 {
588 const Real fd = 0.5 *
589 (PorousFlowVanGenuchten::dcapillaryPressureHys(
590 sat + eps, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power) -
591 PorousFlowVanGenuchten::dcapillaryPressureHys(
592 sat - eps, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power)) /
593 eps;
594 if (std::abs(fd) > 10)
595 EXPECT_NEAR(1.0,
596 fd / PorousFlowVanGenuchten::d2capillaryPressureHys(
597 sat, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power),
598 1.0E-5);
599 else
600 EXPECT_NEAR(fd,
601 PorousFlowVanGenuchten::d2capillaryPressureHys(
602 sat, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power),
603 1.0E-5);
604 };
605
606 // different low extension
607 for (const auto & sat : sats)
608 {
609 const Real fd = 0.5 *
610 (PorousFlowVanGenuchten::dcapillaryPressureHys(
611 sat + eps, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power) -
612 PorousFlowVanGenuchten::dcapillaryPressureHys(
613 sat - eps, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power)) /
614 eps;
615 if (std::abs(fd) > 10)
616 EXPECT_NEAR(1.0,
617 fd / PorousFlowVanGenuchten::d2capillaryPressureHys(
618 sat, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power),
619 1.0E-5);
620 else
621 EXPECT_NEAR(fd,
622 PorousFlowVanGenuchten::d2capillaryPressureHys(
623 sat, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power),
624 1.0E-5);
625 };
626
627 // different lower extension
628 for (const auto & sat : sats)
629 {
630 const Real fd = 0.5 *
631 (PorousFlowVanGenuchten::dcapillaryPressureHys(
632 sat + eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power) -
633 PorousFlowVanGenuchten::dcapillaryPressureHys(
634 sat - eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power)) /
635 eps;
636 if (std::abs(fd) > 10)
637 EXPECT_NEAR(1.0,
638 fd / PorousFlowVanGenuchten::d2capillaryPressureHys(
639 sat, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power),
640 1.0E-5);
641 else
642 EXPECT_NEAR(fd,
643 PorousFlowVanGenuchten::d2capillaryPressureHys(
644 sat, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power),
645 1.0E-5);
646 };
647
648 // different upper extension (cannot evaluate derivative at s = 0.63)
649 sats = {0.01, 0.1, 0.15, 0.2, 0.5, 0.8, 0.99, 1.0};
650 for (const auto & sat : sats)
651 {
652 const Real fd = 0.5 *
653 (PorousFlowVanGenuchten::dcapillaryPressureHys(
654 sat + eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none2) -
655 PorousFlowVanGenuchten::dcapillaryPressureHys(
656 sat - eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none2)) /
657 eps;
658 if (std::abs(fd) > 10)
659 EXPECT_NEAR(1.0,
660 fd / PorousFlowVanGenuchten::d2capillaryPressureHys(
661 sat, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none2),
662 1.0E-5);
663 else
664 EXPECT_NEAR(fd,
665 PorousFlowVanGenuchten::d2capillaryPressureHys(
666 sat, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none2),
667 1.0E-5);
668 };
669}
670
672TEST(PorousFlowVanGenuchtenTest, dsathys)
673{
674 const Real eps = 1E-9;
675
676 // first define some extensions that do not actually induce any extension, so the non-extended Pc
677 // can be checked
678 std::vector<Real> pcs{-1.1, 1234.0};
679 for (const auto & pc : pcs)
680 {
681 EXPECT_NEAR(
682 0.0,
683 PorousFlowVanGenuchten::dsaturationHys(pc, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
684 1.0E-5);
685 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::dsaturationHys(pc, 0.1, 0.3, 10.0, 1.9), 1.0E-5);
686 }
687
688 pcs = {0.7, 0.2, 0.07};
689 for (const auto & pc : pcs)
690 {
691 const Real fd = 0.5 *
692 (PorousFlowVanGenuchten::saturationHys(
693 pc + eps, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext) -
694 PorousFlowVanGenuchten::saturationHys(
695 pc - eps, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext)) /
696 eps;
697 EXPECT_NEAR(
698 fd,
699 PorousFlowVanGenuchten::dsaturationHys(pc, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
700 1.0E-5);
701 EXPECT_NEAR(fd, PorousFlowVanGenuchten::dsaturationHys(pc, 0.1, 0.3, 10.0, 1.9), 1.0E-5);
702 }
703
704 // now with low and high extensions
705 pcs = {31.385947046636815,
706 4.5861935551774735,
707 1.5754562501536364,
708 0.7233512030263158,
709 0.11727884570711045,
710 0.05300654102157442,
711 0.006681337544884095};
712 for (const auto & pc : pcs)
713 {
714 const Real fd = 0.5 *
715 (PorousFlowVanGenuchten::saturationHys(
716 pc + eps, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power) -
717 PorousFlowVanGenuchten::saturationHys(
718 pc - eps, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power)) /
719 eps;
720 EXPECT_NEAR(fd,
721 PorousFlowVanGenuchten::dsaturationHys(
722 pc, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power),
723 1.0E-5);
724 }
725
726 // different lower extension
727 pcs = {3.9304232526771696,
728 2.8885549236001244,
729 1.5730646091089062,
730 0.7233512030263158,
731 0.11727884570711045,
732 0.05300654102157442,
733 0.006681337544884095};
734 for (const auto & pc : pcs)
735 {
736 const Real fd = 0.5 *
737 (PorousFlowVanGenuchten::saturationHys(
738 pc + eps, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power) -
739 PorousFlowVanGenuchten::saturationHys(
740 pc - eps, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power)) /
741 eps;
742 EXPECT_NEAR(fd,
743 PorousFlowVanGenuchten::dsaturationHys(
744 pc, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power),
745 1.0E-5);
746 }
747
748 // different lower extension
749 pcs = {2.5, 0.7233512030263158, 0.11727884570711045, 0.05300654102157442, 0.006681337544884095};
750 for (const auto & pc : pcs)
751 {
752 const Real fd = 0.5 *
753 (PorousFlowVanGenuchten::saturationHys(
754 pc + eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power) -
755 PorousFlowVanGenuchten::saturationHys(
756 pc - eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power)) /
757 eps;
758 EXPECT_NEAR(fd,
759 PorousFlowVanGenuchten::dsaturationHys(
760 pc, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power),
761 1.0E-5);
762 }
763
764 // different upper extension
765 pcs = {2.5, 0.7233512030263158, 0.11727884570711045, 0.05, 0.04};
766 for (const auto & pc : pcs)
767 {
768 const Real fd = 0.5 *
769 (PorousFlowVanGenuchten::saturationHys(
770 pc + eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none) -
771 PorousFlowVanGenuchten::saturationHys(
772 pc - eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none)) /
773 eps;
774 EXPECT_NEAR(fd,
775 PorousFlowVanGenuchten::dsaturationHys(
776 pc, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none),
777 1.0E-5);
778 }
779}
780
782TEST(PorousFlowVanGenuchtenTest, d2sathys)
783{
784 const Real eps = 1E-9;
785
786 // first define some extensions that do not actually induce any extension, so the non-extended Pc
787 // can be checked
788 std::vector<Real> pcs{-1.1, 1234.0};
789 for (const auto & pc : pcs)
790 {
791 EXPECT_NEAR(
792 0.0,
793 PorousFlowVanGenuchten::d2saturationHys(pc, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
794 1.0E-5);
795 EXPECT_NEAR(0.0, PorousFlowVanGenuchten::d2saturationHys(pc, 0.1, 0.3, 10.0, 1.9), 1.0E-5);
796 }
797
798 pcs = {0.7, 0.2, 0.07};
799 for (const auto & pc : pcs)
800 {
801 const Real fd = 0.5 *
802 (PorousFlowVanGenuchten::dsaturationHys(
803 pc + eps, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext) -
804 PorousFlowVanGenuchten::dsaturationHys(
805 pc - eps, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext)) /
806 eps;
807 EXPECT_NEAR(
808 fd,
809 PorousFlowVanGenuchten::d2saturationHys(pc, 0.1, 0.3, 10.0, 1.9, no_low_ext, no_high_ext),
810 1.0E-5);
811 EXPECT_NEAR(fd, PorousFlowVanGenuchten::d2saturationHys(pc, 0.1, 0.3, 10.0, 1.9), 1.0E-5);
812 }
813
814 // now with low and high extensions
815 pcs = {31.385947046636815,
816 4.5861935551774735,
817 1.5754562501536364,
818 0.7233512030263158,
819 0.11727884570711045,
820 0.06300654102157442,
821 0.006681337544884095};
822 for (const auto & pc : pcs)
823 {
824 const Real fd = 0.5 *
825 (PorousFlowVanGenuchten::dsaturationHys(
826 pc + eps, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power) -
827 PorousFlowVanGenuchten::dsaturationHys(
828 pc - eps, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power)) /
829 eps;
830 EXPECT_NEAR(fd,
831 PorousFlowVanGenuchten::d2saturationHys(
832 pc, 0.1, 0.3, 10.0, 1.9, low_ext_exp, high_ext_power),
833 1.0E-5);
834 }
835
836 // different lower extension
837 pcs = {3.9304232526771696,
838 2.8885549236001244,
839 1.5730646091089062,
840 0.7233512030263158,
841 0.11727884570711045,
842 0.06300654102157442,
843 0.006681337544884095};
844 for (const auto & pc : pcs)
845 {
846 const Real fd = 0.5 *
847 (PorousFlowVanGenuchten::dsaturationHys(
848 pc + eps, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power) -
849 PorousFlowVanGenuchten::dsaturationHys(
850 pc - eps, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power)) /
851 eps;
852 EXPECT_NEAR(fd,
853 PorousFlowVanGenuchten::d2saturationHys(
854 pc, 0.1, 0.3, 10.0, 1.9, low_ext_quad, high_ext_power),
855 1.0E-5);
856 }
857
858 // different lower extension
859 pcs = {2.5, 0.7233512030263158, 0.11727884570711045, 0.06300654102157442, 0.006681337544884095};
860 for (const auto & pc : pcs)
861 {
862 const Real fd = 0.5 *
863 (PorousFlowVanGenuchten::dsaturationHys(
864 pc + eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power) -
865 PorousFlowVanGenuchten::dsaturationHys(
866 pc - eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power)) /
867 eps;
868 EXPECT_NEAR(fd,
869 PorousFlowVanGenuchten::d2saturationHys(
870 pc, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_power),
871 1.0E-5);
872 }
873
874 // different upper extension
875 pcs = {2.5, 0.7233512030263158, 0.11727884570711045, 0.05, 0.04};
876 for (const auto & pc : pcs)
877 {
878 const Real fd = 0.5 *
879 (PorousFlowVanGenuchten::dsaturationHys(
880 pc + eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none) -
881 PorousFlowVanGenuchten::dsaturationHys(
882 pc - eps, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none)) /
883 eps;
884 EXPECT_NEAR(fd,
885 PorousFlowVanGenuchten::d2saturationHys(
886 pc, 0.1, 0.3, 10.0, 1.9, low_ext_none, high_ext_none),
887 1.0E-5);
888 }
889}
890
892TEST(PorousFlowVanGenuchtenTest, relativePermeabilityHys)
893{
894 // Tests for sldel = 0.5 >= slr = 0.2
895 EXPECT_NEAR(0.0,
896 PorousFlowVanGenuchten::relativePermeabilityHys(
897 0.1, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
898 1.0E-12); // zero because sl < slr
899 EXPECT_NEAR(0.0,
900 PorousFlowVanGenuchten::relativePermeabilityHys(
901 0.2, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
902 1.0E-12); // zero because sl = slr
903 EXPECT_NEAR(0.002847974503642625,
904 PorousFlowVanGenuchten::relativePermeabilityHys(
905 0.3, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
906 1.0E-12); // drying because sl <= sldel
907 EXPECT_NEAR(0.002847974503642625,
908 PorousFlowVanGenuchten::relativePermeabilityHys(
909 0.3, 0.2, 0.0, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
910 1.0E-12); // drying because sgrdel = 0
911 EXPECT_NEAR(0.05828443549816974,
912 PorousFlowVanGenuchten::relativePermeabilityHys(
913 0.5, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
914 1.0E-12); // drying because sl <= sldel
915 EXPECT_NEAR(0.10056243969693253,
916 PorousFlowVanGenuchten::relativePermeabilityHys(
917 0.55, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
918 1.0E-12); // wetting
919 EXPECT_NEAR(0.08943153676247108,
920 PorousFlowVanGenuchten::relativePermeabilityHys(
921 0.55, 0.2, 0.0, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
922 1.0E-12); // drying because sgrdel = 0
923 EXPECT_NEAR(0.3159149169921876,
924 PorousFlowVanGenuchten::relativePermeabilityHys(
925 0.8, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
926 1.0E-12); // cubic
927 EXPECT_NEAR(0.8011290515690178,
928 PorousFlowVanGenuchten::relativePermeabilityHys(
929 0.95, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
930 1.0E-12); // drying, because sl > cubic modification region
931
932 // Tests for sldel = 0.15 < slr = 0.2
933 EXPECT_NEAR(0.0,
934 PorousFlowVanGenuchten::relativePermeabilityHys(
935 0.1, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
936 1.0E-12); // zero because sl < slr
937 EXPECT_NEAR(0.0,
938 PorousFlowVanGenuchten::relativePermeabilityHys(
939 0.2, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
940 1.0E-12); // zero because sl = slr
941 EXPECT_NEAR(0.005858393312913491,
942 PorousFlowVanGenuchten::relativePermeabilityHys(
943 0.3, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
944 1.0E-12); // wetting
945 EXPECT_NEAR(0.002847974503642625,
946 PorousFlowVanGenuchten::relativePermeabilityHys(
947 0.3, 0.2, 0.0, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
948 1.0E-12); // drying because sgrdel = 0
949 EXPECT_NEAR(0.1363562523627124,
950 PorousFlowVanGenuchten::relativePermeabilityHys(
951 0.55, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
952 1.0E-12); // wetting
953 EXPECT_NEAR(0.31,
954 PorousFlowVanGenuchten::relativePermeabilityHys(
955 0.7, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
956 1.0E-12); // cubic
957 EXPECT_NEAR(0.8011290515690178,
958 PorousFlowVanGenuchten::relativePermeabilityHys(
959 0.95, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
960 1.0E-12); // drying, because sl > cubic modification region
961}
962
964TEST(PorousFlowVanGenuchtenTest, drelativePermeabilityHys)
965{
966 // essentially follow the testing strategy in the relativePermeabilityHys Test, but do not test
967 // the derivative at the points where it is discontinuous
968 EXPECT_NEAR(0.0,
969 PorousFlowVanGenuchten::drelativePermeabilityHys(
970 0.1, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
971 1.0E-12);
972 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
973 0.3, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
974 (PorousFlowVanGenuchten::relativePermeabilityHys(
975 0.3 + eps, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
976 PorousFlowVanGenuchten::relativePermeabilityHys(
977 0.3, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
978 eps,
979 1.0E-5);
980 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
981 0.3, 0.2, 0.0, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
982 (PorousFlowVanGenuchten::relativePermeabilityHys(
983 0.3 + eps, 0.2, 0.0, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
984 PorousFlowVanGenuchten::relativePermeabilityHys(
985 0.3, 0.2, 0.0, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
986 eps,
987 1.0E-5);
988 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
989 0.55, 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.55 + eps, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
992 PorousFlowVanGenuchten::relativePermeabilityHys(
993 0.55, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
994 eps,
995 1.0E-5);
996 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
997 0.55, 0.2, 0.0, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
998 (PorousFlowVanGenuchten::relativePermeabilityHys(
999 0.55 + eps, 0.2, 0.0, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
1000 PorousFlowVanGenuchten::relativePermeabilityHys(
1001 0.55, 0.2, 0.0, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
1002 eps,
1003 1.0E-5);
1004 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
1005 0.8, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
1006 (PorousFlowVanGenuchten::relativePermeabilityHys(
1007 0.8 + eps, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
1008 PorousFlowVanGenuchten::relativePermeabilityHys(
1009 0.8, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
1010 eps,
1011 1.0E-5);
1012 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
1013 0.95, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
1014 (PorousFlowVanGenuchten::relativePermeabilityHys(
1015 0.95 + eps, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
1016 PorousFlowVanGenuchten::relativePermeabilityHys(
1017 0.95, 0.2, 0.15, 0.25, 0.5, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
1018 eps,
1019 1.0E-5);
1020
1021 EXPECT_NEAR(0.0,
1022 PorousFlowVanGenuchten::drelativePermeabilityHys(
1023 0.1, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
1024 1.0E-5);
1025 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
1026 0.3, 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.3 + eps, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
1029 PorousFlowVanGenuchten::relativePermeabilityHys(
1030 0.3, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
1031 eps,
1032 1.0E-5);
1033 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
1034 0.3, 0.2, 0.0, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
1035 (PorousFlowVanGenuchten::relativePermeabilityHys(
1036 0.3 + eps, 0.2, 0.0, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
1037 PorousFlowVanGenuchten::relativePermeabilityHys(
1038 0.3, 0.2, 0.0, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
1039 eps,
1040 1.0E-5);
1041 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
1042 0.55, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
1043 (PorousFlowVanGenuchten::relativePermeabilityHys(
1044 0.55 + eps, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
1045 PorousFlowVanGenuchten::relativePermeabilityHys(
1046 0.55, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
1047 eps,
1048 1.0E-5);
1049 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
1050 0.7, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
1051 (PorousFlowVanGenuchten::relativePermeabilityHys(
1052 0.7 + eps, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
1053 PorousFlowVanGenuchten::relativePermeabilityHys(
1054 0.7, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
1055 eps,
1056 1.0E-5);
1057 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityHys(
1058 0.95, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2),
1059 (PorousFlowVanGenuchten::relativePermeabilityHys(
1060 0.95 + eps, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2) -
1061 PorousFlowVanGenuchten::relativePermeabilityHys(
1062 0.95, 0.2, 0.15, 0.25, 0.15, 0.9, 0.9, 0.3, 0.5, 0.45, 2.2)) /
1063 eps,
1064 1.0E-5);
1065}
1066
1068TEST(PorousFlowVanGenuchtenTest, relativePermeabilityNWHys)
1069{
1070 // Tests for sldel = 0.5 >= slr = 0.2
1071 EXPECT_NEAR(0.91875,
1072 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1073 0.1, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1074 1.0E-12); // cubic extension region because sl < slr
1075 EXPECT_NEAR(1.0,
1076 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1077 0.1, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 1.0, -0.75),
1078 1.0E-12); // cubic extension region because sl < slr, but k_rg_max = 1
1079 EXPECT_NEAR(0.8,
1080 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1081 0.2, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1082 1.0E-12); // at boundary of extension since sl = slr
1083 EXPECT_NEAR(0.6368139633459233,
1084 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1085 0.3, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1086 1.0E-12); // drying because sl <= sldel
1087 EXPECT_NEAR(0.6368139633459233,
1088 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1089 0.3, 0.2, 0.0, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1090 1.0E-12); // drying because sgrdel = 0
1091 EXPECT_NEAR(0.33222054246699,
1092 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1093 0.5, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1094 1.0E-12); // drying because sl <= sldel
1095 EXPECT_NEAR(0.24397255389213568,
1096 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1097 0.55, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1098 1.0E-12); // wetting
1099 EXPECT_NEAR(0.2691316236913443,
1100 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1101 0.55, 0.2, 0.0, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1102 1.0E-12); // drying because sgrdel = 0
1103 EXPECT_NEAR(0.01509148277290631,
1104 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1105 0.85, 0.2, 0.05, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1106 1.0E-12); // wetting
1107 EXPECT_NEAR(0.0,
1108 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1109 0.85, 0.2, 0.18, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1110 1.0E-12); // sl > 1 - s_gr_del
1111
1112 // Tests for sldel = 0.15 < slr = 0.2
1113 EXPECT_NEAR(0.91875,
1114 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1115 0.1, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1116 1.0E-12); // cubic extension region because sl < slr
1117 EXPECT_NEAR(1.0,
1118 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1119 0.1, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 1.0, -0.75),
1120 1.0E-12); // cubic extension region because sl < slr, but k_rg_max = 1
1121 EXPECT_NEAR(0.8,
1122 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1123 0.2, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1124 1.0E-12); // at boundary of extension since sl = slr
1125 EXPECT_NEAR(0.5616820966427329,
1126 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1127 0.3, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1128 1.0E-12); // wetting
1129 EXPECT_NEAR(0.6368139633459233,
1130 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1131 0.3, 0.2, 0.0, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1132 1.0E-12); // wetting = drying because sgrdel = 0
1133 EXPECT_NEAR(0.11086138435665097,
1134 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1135 0.55, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1136 1.0E-12); // wetting
1137 EXPECT_NEAR(0.0,
1138 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1139 0.76, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1140 1.0E-12); // always zero for large saturation
1141}
1142
1144TEST(PorousFlowVanGenuchtenTest, drelativePermeabilityNWHys)
1145{
1146 // Testing strategy follows the relativePermeabilityNWHys strategy, except the points at which the
1147 // derivative is not defined are not tested
1148 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1149 0.1, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1150 (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1151 0.1 + eps, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75) -
1152 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1153 0.1, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75)) /
1154 eps,
1155 1.0E-5); // cubic extension region because sl < slr
1156 EXPECT_NEAR(0.0,
1157 PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1158 0.1, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 1.0, -0.75),
1159 1.0E-12); // cubic extension region because sl < slr, but k_rg_max = 1
1160 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1161 0.3, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1162 (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1163 0.3 + eps, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75) -
1164 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1165 0.3, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75)) /
1166 eps,
1167 1.0E-5); // drying because sl <= sldel
1168 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1169 0.3, 0.2, 0.0, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1170 (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1171 0.3 + eps, 0.2, 0.0, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75) -
1172 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1173 0.3, 0.2, 0.0, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75)) /
1174 eps,
1175 1.0E-5); // drying because sgrdel = 0
1176 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1177 0.45, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1178 (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1179 0.45 + eps, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75) -
1180 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1181 0.45, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75)) /
1182 eps,
1183 1.0E-5); // drying because sl <= sldel
1184 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1185 0.55, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1186 (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1187 0.55 + eps, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75) -
1188 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1189 0.55, 0.2, 0.15, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75)) /
1190 eps,
1191 1.0E-5); // wetting
1192 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1193 0.55, 0.2, 0.0, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1194 (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1195 0.55 + eps, 0.2, 0.0, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75) -
1196 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1197 0.55, 0.2, 0.0, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75)) /
1198 eps,
1199 1.0E-5); // drying because sgrdel = 0
1200 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1201 0.85, 0.2, 0.05, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1202 (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1203 0.85 + eps, 0.2, 0.05, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75) -
1204 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1205 0.85, 0.2, 0.05, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75)) /
1206 eps,
1207 1.0E-5); // wetting
1208 EXPECT_NEAR(0.0,
1209 PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1210 0.85, 0.2, 0.18, 0.25, 0.5, 0.9, 0.3, 0.8, -0.75),
1211 1.0E-12); // sl > 1 - s_gr_del
1212
1213 // Tests for sldel = 0.15 < slr = 0.2
1214 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1215 0.1, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1216 (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1217 0.1 + eps, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75) -
1218 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1219 0.1, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75)) /
1220 eps,
1221 1.0E-5); // cubic extension region because sl < slr
1222 EXPECT_NEAR(0.0,
1223 PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1224 0.1, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 1.0, -0.75),
1225 1.0E-12); // cubic extension region because sl < slr, but k_rg_max = 1
1226 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1227 0.3, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1228 (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1229 0.3 + eps, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75) -
1230 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1231 0.3, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75)) /
1232 eps,
1233 1.0E-5); // wetting
1234 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1235 0.3, 0.2, 0.0, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1236 (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1237 0.3 + eps, 0.2, 0.0, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75) -
1238 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1239 0.3, 0.2, 0.0, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75)) /
1240 eps,
1241 1.0E-5); // wetting = drying because sgrdel = 0
1242 EXPECT_NEAR(PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1243 0.55, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1244 (PorousFlowVanGenuchten::relativePermeabilityNWHys(
1245 0.55 + eps, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75) -
1246 PorousFlowVanGenuchten::relativePermeabilityNWHys(
1247 0.55, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75)) /
1248 eps,
1249 1.0E-5); // wetting
1250 EXPECT_NEAR(0.0,
1251 PorousFlowVanGenuchten::drelativePermeabilityNWHys(
1252 0.76, 0.2, 0.18, 0.25, 0.15, 0.9, 0.3, 0.8, -0.75),
1253 1.0E-12); // always zero for large saturation
1254}
ADReal
DualNumber< Real, DNDerivativeType, true > ADReal
dataLoad
void dataLoad(std::istream &stream, LineSegment &l, void *context)
dataStore
void dataStore(std::ostream &stream, LineSegment &l, void *context)
low_ext_Pc
const Real low_ext_Pc
Definition PorousFlowVanGenuchtenTest.C:19
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:43
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:15
TEST
TEST(PorousFlowVanGenuchtenTest, extensionsDataIO)
Definition PorousFlowVanGenuchtenTest.C:58
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:40
high_ext_Pc
const Real high_ext_Pc
Definition PorousFlowVanGenuchtenTest.C:41
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:20
low_ext_dPc
const Real low_ext_dPc
Definition PorousFlowVanGenuchtenTest.C:21
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:554
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:142
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:346
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:811
PorousFlowVanGenuchten::dRelativePermeability
Real dRelativePermeability(Real seff, Real m)
Derivative of relative permeability with respect to effective saturation.
Definition PorousFlowVanGenuchten.C:161
PorousFlowVanGenuchten::relativePermeability
T relativePermeability(const T &seff, Real m)
Relative permeability as a function of effective saturation.
Definition PorousFlowVanGenuchten.h:103
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:411
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:767
PorousFlowVanGenuchten::d2EffectiveSaturation
Real d2EffectiveSaturation(Real p, Real alpha, Real m)
Second derivative of effective saturation wrt porepressure.
Definition PorousFlowVanGenuchten.C:93
PorousFlowVanGenuchten::dEffectiveSaturation
Real dEffectiveSaturation(Real p, Real alpha, Real m)
Derivative of effective saturation wrt porepressure.
Definition PorousFlowVanGenuchten.C:78
PorousFlowVanGenuchten::dCapillaryPressure
Real dCapillaryPressure(Real seff, Real alpha, Real m, Real pc_max)
Derivative of capillary pressure wrt effective saturation.
Definition PorousFlowVanGenuchten.C:124
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:227
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:689
PorousFlowVanGenuchten::capillaryPressure
Real capillaryPressure(Real seff, Real alpha, Real m, Real pc_max)
Capillary pressure as a function of effective saturation.
Definition PorousFlowVanGenuchten.C:110
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:285
PorousFlowVanGenuchten::d2CapillaryPressure
Real d2CapillaryPressure(Real seff, Real alpha, Real m, Real pc_max)
Second derivative of capillary pressure wrt effective saturation.
Definition PorousFlowVanGenuchten.C:140
PorousFlowVanGenuchten::effectiveSaturation
Real effectiveSaturation(Real p, Real alpha, Real m)
Effective saturation as a function of porepressure.
Definition PorousFlowVanGenuchten.C:63
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:633
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:209
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:480
PorousFlowVanGenuchten::d2RelativePermeability
Real d2RelativePermeability(Real seff, Real m)
Second derivative of relative permeability with respect to effective saturation.
Definition PorousFlowVanGenuchten.C:176
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:194
PorousFlowVanGenuchten::HighCapillaryPressureExtension
Parameters associated with the extension of the hysteretic wetting capillary pressure function to hig...
Definition PorousFlowVanGenuchten.h:218
PorousFlowVanGenuchten::LowCapillaryPressureExtension
Parameters associated with the extension of the hysteretic capillary pressure function to low saturat...
Definition PorousFlowVanGenuchten.h:183
Generated on Fri Sep 11 2026 20:36:38 for https://mooseframework.inl.gov by doxygen 1.9.8