van Genuchten effective saturation, capillary pressure and relative permeability functions. More...
Classes | |
struct | HighCapillaryPressureExtension |
Parameters associated with the extension of the hysteretic wetting capillary pressure function to high saturation values @ ExtensionStrategy the type of extension used @ S liquid saturation at the point of extension @ Pc capillary pressure at the point of extension @ dPc d(Pc)/dS at the point of extension. More... | |
struct | LowCapillaryPressureExtension |
Parameters associated with the extension of the hysteretic capillary pressure function to low saturation values @ ExtensionStrategy the type of extension used @ S liquid saturation at the point of extension @ Pc capillary pressure at the point of extension @ dPc d(Pc)/dS at the point of extension. More... | |
Functions | |
Real | effectiveSaturation (Real p, Real alpha, Real m) |
Effective saturation as a function of porepressure. More... | |
Real | dEffectiveSaturation (Real p, Real alpha, Real m) |
Derivative of effective saturation wrt porepressure. More... | |
Real | d2EffectiveSaturation (Real p, Real alpha, Real m) |
Second derivative of effective saturation wrt porepressure. More... | |
Real | capillaryPressure (Real seff, Real alpha, Real m, Real pc_max) |
Capillary pressure as a function of effective saturation. More... | |
Real | dCapillaryPressure (Real seff, Real alpha, Real m, Real pc_max) |
Derivative of capillary pressure wrt effective saturation. More... | |
Real | d2CapillaryPressure (Real seff, Real alpha, Real m, Real pc_max) |
Second derivative of capillary pressure wrt effective saturation. More... | |
template<typename T > | |
T | relativePermeability (const T &seff, Real m) |
Relative permeability as a function of effective saturation. More... | |
Real | dRelativePermeability (Real seff, Real m) |
Derivative of relative permeability with respect to effective saturation. More... | |
Real | d2RelativePermeability (Real seff, Real m) |
Second derivative of relative permeability with respect to effective saturation. More... | |
template<typename T > | |
T | relativePermeabilityNW (const T &seff, Real m) |
Relative permeability for a non-wetting phase as a function of effective saturation. More... | |
Real | dRelativePermeabilityNW (Real seff, Real m) |
Derivative of relative permeability for a non-wetting phase with respect to effective saturation. More... | |
Real | d2RelativePermeabilityNW (Real seff, Real m) |
Second derivative of relative permeability for a non-wetting phase with respect to effective saturation. More... | |
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 Doughty2008). More... | |
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. More... | |
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. More... | |
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), which is the inverse of capillaryPressureHys. More... | |
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. More... | |
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. More... | |
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. More... | |
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. More... | |
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. More... | |
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. More... | |
van Genuchten effective saturation, capillary pressure and relative permeability functions.
Note: effective saturation is provided as a function of porepressure, not capillary pressure. Note: capillary pressure and relative permeability are functions of effective saturation. The derivatives are therefore given wrt effective saturation. These derivatives must be multiplied by the derivative of effective saturation wrt the true saturation in objects using these relations.
Based on van Genuchten, M. Th., A closed for equation for predicting the hydraulic conductivity of unsaturated soils, Soil Sci. Soc., 44, 892-898 (1980).
Capillary pressure as a function of effective saturation.
seff | effective saturation |
alpha | van Genuchten parameter |
m | van Genuchten exponent |
pc_max | maximum capillary pressure (Pa) |
Definition at line 63 of file PorousFlowVanGenuchten.C.
Referenced by PorousFlowCapillaryPressureVG::capillaryPressureCurve(), and TEST().
Real PorousFlowVanGenuchten::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 Doughty2008).
NOTE: this function is undefined for sl < 0 and sl > 1, so you MUST ensure 0 <= sl <= 1 NOTE: this returns a non-negative quantity.
sl | liquid saturation. 0 <= sl <= 1 |
slmin | value of liquid sat where the van Genuchten expression -> infinity. 0 <= slmin < 1 |
sgrdel | value of gas saturation where van Genuchten expression -> 0. slmin < 1 - Sgrdel <= 1 |
alpha | van Genuchten alpha parameter, with dimensions 1/Pa. alpha > 0 |
n | van Genuchten n parameter. n > 1 |
low_ext | strategy and parameters to use for the extension in the small-saturation region (defaults to no extension: this default is not recommended for simulations of real phenomena) |
high_ext | strategy and parameters to use for the extension in the high-saturation region (defaults to no extension: this default is not recommended for simulations of real phenomena) |
Definition at line 180 of file PorousFlowVanGenuchten.C.
Referenced by PorousFlowHystereticCapillaryPressure::capillaryPressureQp(), PorousFlowHystereticCapillaryPressure::computeTurningPointInfo(), PorousFlowHystereticCapillaryPressure::firstOrderWettingPc(), PorousFlowHystereticCapillaryPressure::initQpStatefulProperties(), PorousFlowHystereticCapillaryPressure::secondOrderDryingPc(), and TEST().
Second derivative of capillary pressure wrt effective saturation.
seff | effective saturation |
alpha | van Genuchten parameter |
m | van Genuchten exponent |
pc_max | maximum capillary pressure (Pa) |
Definition at line 93 of file PorousFlowVanGenuchten.C.
Referenced by PorousFlowCapillaryPressureVG::d2CapillaryPressureCurve(), and TEST().
Real PorousFlowVanGenuchten::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.
NOTE: this function is undefined for sl < 0 and sl > 1, so you MUST ensure 0 <= sl <= 1
sl | liquid saturation. 0 <= sl <= 1 |
slmin | value of liquid sat where the van Genuchten expression -> infinity. 0 <= slmin < 1 |
sgrdel | value of gas saturation where van Genuchten expression -> 0. slmin < 1 - Sgrdel <= 1 |
alpha | van Genuchten alpha parameter, with dimensions 1/Pa. alpha > 0 |
n | van Genuchten n parameter. n > 1 |
low_ext | strategy and parameters to use for the extension in the small-saturation region (defaults to no extension: this default is not recommended for simulations of real phenomena) |
high_ext | strategy and parameters to use for the extension in the high-saturation region (defaults to no extension: this default is not recommended for simulations of real phenomena) |
Definition at line 299 of file PorousFlowVanGenuchten.C.
Referenced by PorousFlowHystereticCapillaryPressure::d2capillaryPressureQp(), PorousFlowHystereticCapillaryPressure::d2firstOrderWettingPc(), PorousFlowHystereticCapillaryPressure::d2secondOrderDryingPc(), and TEST().
Second derivative of effective saturation wrt porepressure.
p | porepressure |
alpha | van Genuchten parameter |
m | van Genuchten exponent |
Definition at line 46 of file PorousFlowVanGenuchten.C.
Referenced by PorousFlowCapillaryPressureVG::d2EffectiveSaturation(), and TEST().
Second derivative of relative permeability with respect to effective saturation.
seff | effective saturation |
m | van Genuchten exponent |
Definition at line 129 of file PorousFlowVanGenuchten.C.
Referenced by TEST().
Second derivative of relative permeability for a non-wetting phase with respect to effective saturation.
seff | effective saturation |
m | van Genuchten exponent |
Definition at line 162 of file PorousFlowVanGenuchten.C.
Referenced by TEST().
Real PorousFlowVanGenuchten::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.
pc | capillary pressure. 0 <= pc |
slmin | value of liquid sat where the van Genuchten expression -> infinity. 0 <= slmin < 1 |
sgrdel | value of gas saturation where van Genuchten expression -> 0. slmin < 1 - Sgrdel <= 1 |
alpha | van Genuchten alpha parameter, with dimensions 1/Pa. alpha > 0 |
n | van Genuchten n parameter. n > 1 |
low_ext | strategy and parameters to use for the extension in the small-saturation region (defaults to no extension: this default is not recommended for simulations of real phenomena) |
high_ext | strategy and parameters to use for the extension in the high-saturation region (defaults to no extension: this default is not recommended for simulations of real phenomena) |
Definition at line 507 of file PorousFlowVanGenuchten.C.
Referenced by PorousFlowHystereticCapillaryPressure::d2firstOrderWettingSat(), PorousFlowHystereticCapillaryPressure::d2liquidSaturationQp(), PorousFlowHystereticCapillaryPressure::d2secondOrderDryingSat(), and TEST().
Derivative of capillary pressure wrt effective saturation.
seff | effective saturation |
alpha | van Genuchten parameter |
m | van Genuchten exponent |
pc_max | maximum capillary pressure (Pa) |
Definition at line 77 of file PorousFlowVanGenuchten.C.
Referenced by PorousFlowCapillaryPressureVG::dCapillaryPressureCurve(), and TEST().
Real PorousFlowVanGenuchten::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.
NOTE: this function is undefined for sl < 0 and sl > 1, so you MUST ensure 0 <= sl <= 1 NOTE: this returns a negative quantity.
sl | liquid saturation. 0 <= sl <= 1 |
slmin | value of liquid sat where the van Genuchten expression -> infinity. 0 <= slmin < 1 |
sgrdel | value of gas saturation where van Genuchten expression -> 0. slmin < 1 - Sgrdel <= 1 |
alpha | van Genuchten alpha parameter, with dimensions 1/Pa. alpha > 0 |
n | van Genuchten n parameter. n > 1 |
low_ext | strategy and parameters to use for the extension in the small-saturation region (defaults to no extension: this default is not recommended for simulations of real phenomena) |
high_ext | strategy and parameters to use for the extension in the high-saturation region (defaults to no extension: this default is not recommended for simulations of real phenomena) |
Definition at line 238 of file PorousFlowVanGenuchten.C.
Referenced by PorousFlowHystereticCapillaryPressure::computeTurningPointInfo(), PorousFlowHystereticCapillaryPressure::dcapillaryPressureQp(), PorousFlowHystereticCapillaryPressure::dfirstOrderWettingPc(), PorousFlowHystereticCapillaryPressure::dsecondOrderDryingPc(), and TEST().
Derivative of effective saturation wrt porepressure.
p | porepressure |
alpha | van Genuchten parameter |
m | van Genuchten exponent |
Definition at line 31 of file PorousFlowVanGenuchten.C.
Referenced by PorousFlowCapillaryPressureVG::dEffectiveSaturation(), and TEST().
Derivative of relative permeability with respect to effective saturation.
seff | effective saturation |
m | van Genuchten exponent |
Definition at line 114 of file PorousFlowVanGenuchten.C.
Referenced by PorousFlowRelativePermeabilityVGTempl< is_ad >::dRelativePermeability(), and TEST().
Real PorousFlowVanGenuchten::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.
sl | liquid saturation |
slr | residual liquid saturation. For sl < slr, this function will always return 0 |
sgrdel | value of gas saturation where van Genuchten wetting capillary-pressure expression -> 0. This depends on the turning-point saturation when drying became wetting, using the Land equation |
sgrmax | maximum value possible for sgrdel. This will be equal to sgrdel if the turning-point saturation is small |
sldel | value of the turning-point saturation when drying became wetting |
m | van-Genuchten m parameter |
upper_liquid_param | cubic modification of the wetting relative permeability will occur between upper_liquid_param * (1 - sgrdel) and 1 - 0.5 * sgrdel. 0 < upper_liquid_param <= 1. Usually upper_liquid_param is close to 1 (eg 0.9) |
y0 | value of the unmodified wetting relative permeability at sl = upper_liquid_param * (1 - sgrdel) |
y0p | value of the derivtaive of the unmodified wetting relative permeability at sl = upper_liquid_param * (1 - sgrdel) |
y1 | value of the unmodified wetting relative permeability at sl = 1 - 0.5 * sgrdel |
y1p | value of the derivtaive of the unmodified wetting relative permeability at sl = 1 - 0.5 * sgrdel |
Definition at line 642 of file PorousFlowVanGenuchten.C.
Referenced by PorousFlowHystereticRelativePermeabilityLiquid::computeRelPermQp(), PorousFlowHystereticRelativePermeabilityLiquid::computeTurningPoint0Info(), and TEST().
Derivative of relative permeability for a non-wetting phase with respect to effective saturation.
seff | effective saturation |
m | van Genuchten exponent |
Definition at line 147 of file PorousFlowVanGenuchten.C.
Referenced by PorousFlowRelativePermeabilityVGTempl< is_ad >::dRelativePermeability(), and TEST().
Real PorousFlowVanGenuchten::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.
sl | liquid saturation |
slr | residual liquid saturation. For sl < slr, this function will always return 0 |
sgrdel | value of gas saturation where van Genuchten wetting capillary-pressure expression -> 0. This depends on the turning-point saturation when drying became wetting, using the Land equation |
sgrmax | maximum value possible for sgrdel. This will be equal to sgrdel if the turning-point saturation is small |
sldel | value of the turning-point saturation when drying became wetting |
m | van-Genuchten m parameter |
gamma | index satisfying gamma > 0. Usually gamma = 1/3. |
k_rg_max | Maximum value of unextended gas relative permeability. If no low-saturation extension is used then gas relative permeability = k_rg_max for sl <= slr. Satisfies 0 < k_rg_max <= 1. Frequently k_rg_max = 1 is used |
y0p | Value of the derivative of the low-saturation extension at sl = slr. If an extension is used then the gas relative permeability in the region sl <= slr is a cubic whose value is unity at sl = 0, and derivative is zero at sl = 0 |
Definition at line 764 of file PorousFlowVanGenuchten.C.
Referenced by PorousFlowHystereticRelativePermeabilityGas::computeRelPermQp(), and TEST().
Real PorousFlowVanGenuchten::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.
pc | capillary pressure. 0 <= pc |
slmin | value of liquid sat where the van Genuchten expression -> infinity. 0 <= slmin < 1 |
sgrdel | value of gas saturation where van Genuchten expression -> 0. slmin < 1 - Sgrdel <= 1 |
alpha | van Genuchten alpha parameter, with dimensions 1/Pa. alpha > 0 |
n | van Genuchten n parameter. n > 1 |
low_ext | strategy and parameters to use for the extension in the small-saturation region (defaults to no extension: this default is not recommended for simulations of real phenomena) |
high_ext | strategy and parameters to use for the extension in the high-saturation region (defaults to no extension: this default is not recommended for simulations of real phenomena) |
Definition at line 433 of file PorousFlowVanGenuchten.C.
Referenced by PorousFlowHystereticCapillaryPressure::dfirstOrderWettingSat(), PorousFlowHystereticCapillaryPressure::dliquidSaturationQp(), PorousFlowHystereticCapillaryPressure::dsecondOrderDryingSat(), and TEST().
Effective saturation as a function of porepressure.
Note: seff = 1 for p >= 0
p | porepressure |
alpha | van Genuchten parameter |
m | van Genuchten exponent |
Definition at line 16 of file PorousFlowVanGenuchten.C.
Referenced by PorousFlowCapillaryPressureVG::effectiveSaturation(), and TEST().
T PorousFlowVanGenuchten::relativePermeability | ( | const T & | seff, |
Real | m | ||
) |
Relative permeability as a function of effective saturation.
seff | effective saturation |
m | van Genuchten exponent |
Definition at line 101 of file PorousFlowVanGenuchten.h.
Referenced by PorousFlowRelativePermeabilityVGTempl< is_ad >::relativePermeability(), and TEST().
Real PorousFlowVanGenuchten::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.
sl | liquid saturation |
slr | residual liquid saturation. For sl < slr, this function will always return 0 |
sgrdel | value of gas saturation where van Genuchten wetting capillary-pressure expression -> 0. This depends on the turning-point saturation when drying became wetting, using the Land equation |
sgrmax | maximum value possible for sgrdel. This will be equal to sgrdel if the turning-point saturation is small |
sldel | value of the turning-point saturation when drying became wetting |
m | van-Genuchten m parameter |
upper_liquid_param | cubic modification of the wetting relative permeability will occur between upper_liquid_param * (1 - sgrdel) and 1 - 0.5 * sgrdel. 0 < upper_liquid_param <= 1. Usually upper_liquid_param is close to 1 (eg 0.9) |
y0 | value of the unmodified wetting relative permeability at sl = upper_liquid_param * (1 - sgrdel) |
y0p | value of the derivtaive of the unmodified wetting relative permeability at sl = upper_liquid_param * (1 - sgrdel) |
y1 | value of the unmodified wetting relative permeability at sl = 1 - 0.5 * sgrdel |
y1p | value of the derivtaive of the unmodified wetting relative permeability at sl = 1 - 0.5 * sgrdel |
Definition at line 586 of file PorousFlowVanGenuchten.C.
Referenced by PorousFlowHystereticRelativePermeabilityLiquid::computeRelPermQp(), PorousFlowHystereticRelativePermeabilityLiquid::computeTurningPoint0Info(), and TEST().
T PorousFlowVanGenuchten::relativePermeabilityNW | ( | const T & | seff, |
Real | m | ||
) |
Relative permeability for a non-wetting phase as a function of effective saturation.
seff | effective saturation |
m | van Genuchten exponent |
Definition at line 138 of file PorousFlowVanGenuchten.h.
Referenced by PorousFlowRelativePermeabilityVGTempl< is_ad >::relativePermeability(), and TEST().
Real PorousFlowVanGenuchten::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.
sl | liquid saturation |
slr | residual liquid saturation. For sl < slr, this function will always return 0 |
sgrdel | value of gas saturation where van Genuchten wetting capillary-pressure expression -> 0. This depends on the turning-point saturation when drying became wetting, using the Land equation |
sgrmax | maximum value possible for sgrdel. This will be equal to sgrdel if the turning-point saturation is small |
sldel | value of the turning-point saturation when drying became wetting |
m | van-Genuchten m parameter |
gamma | index satisfying gamma > 0. Usually gamma = 1/3. |
k_rg_max | Maximum value of unextended gas relative permeability. If no low-saturation extension is used then gas relative permeability = k_rg_max for sl <= slr. Satisfies 0 < k_rg_max <= 1. Frequently k_rg_max = 1 is used |
y0p | Value of the derivative of the low-saturation extension at sl = slr. If an extension is used then the gas relative permeability in the region sl <= slr is a cubic whose value is unity at sl = 0, and derivative is zero at sl = 0 |
Definition at line 720 of file PorousFlowVanGenuchten.C.
Referenced by PorousFlowHystereticRelativePermeabilityGas::computeRelPermQp(), and TEST().
Real PorousFlowVanGenuchten::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), which is the inverse of capillaryPressureHys.
pc | capillary pressure. 0 <= pc |
slmin | value of liquid sat where the van Genuchten expression -> infinity. 0 <= slmin < 1 |
sgrdel | value of gas saturation where van Genuchten expression -> 0. slmin < 1 - Sgrdel <= 1 |
alpha | van Genuchten alpha parameter, with dimensions 1/Pa. alpha > 0 |
n | van Genuchten n parameter. n > 1 |
low_ext | strategy and parameters to use for the extension in the small-saturation region (defaults to no extension: this default is not recommended for simulations of real phenomena) |
high_ext | strategy and parameters to use for the extension in the high-saturation region (defaults to no extension: this default is not recommended for simulations of real phenomena) |
Definition at line 364 of file PorousFlowVanGenuchten.C.
Referenced by PorousFlowHystereticCapillaryPressure::computeTurningPointInfo(), PorousFlowHystereticCapillaryPressure::d2secondOrderDryingSat(), PorousFlowHystereticCapillaryPressure::dsecondOrderDryingSat(), PorousFlowHystereticCapillaryPressure::firstOrderWettingSat(), PorousFlowHystereticCapillaryPressure::liquidSaturationQp(), PorousFlowHystereticCapillaryPressure::secondOrderDryingSat(), and TEST().