Loading [MathJax]/extensions/tex2jax.js
OGS
Toggle main menu visibility
Main Page
Related Pages
Namespaces
Namespace List
Namespace Members
All
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
u
v
w
x
y
z
Functions
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
u
v
w
x
y
z
Variables
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
u
v
w
Typedefs
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
v
w
z
Enumerations
b
c
d
e
f
g
i
l
m
n
o
p
s
t
u
v
Enumerator
a
b
c
d
e
f
g
h
l
m
n
p
r
s
t
v
y
Classes
Class List
Class Index
Class Hierarchy
Files
File List
File Members
All
a
b
c
d
e
f
g
i
k
m
n
o
p
r
s
t
v
w
x
Functions
a
c
d
e
f
g
i
m
o
p
r
s
t
v
w
Variables
Typedefs
Enumerations
Macros
b
c
g
m
n
o
p
r
s
t
▼
OGS
►
OpenGeoSys 6.5.4-399-g4e01e7328 source code documentation
►
OGS CTests—Project Files
►
OGS Input File Parameters—Quality Assurance
►
OGS Input File Parameters
BulkMappingDocuPage
Todo List
OGS Input File Parameters—List of incomplete documentation pages
Bibliography
►
Namespaces
►
Classes
▼
Files
▼
File List
►
Applications
►
BaseLib
►
build
►
ChemistryLib
Documentation
►
GeoLib
►
InfoLib
▼
MaterialLib
►
Adsorption
►
Fluid
►
FractureModels
▼
MPL
►
Components
▼
Properties
►
CapillaryPressureSaturation
►
Density
►
Enthalpy
▼
RelativePermeability
►
CreateRelPermBrooksCorey.cpp
►
CreateRelPermBrooksCorey.h
►
CreateRelPermBrooksCoreyNonwettingPhase.cpp
►
CreateRelPermBrooksCoreyNonwettingPhase.h
►
CreateRelPermGeneralizedPower.cpp
►
CreateRelPermGeneralizedPower.h
►
CreateRelPermGeneralizedPowerNonwettingPhase.cpp
►
CreateRelPermGeneralizedPowerNonwettingPhase.h
►
CreateRelPermLiakopoulos.cpp
►
CreateRelPermLiakopoulos.h
►
CreateRelPermNonWettingPhaseVanGenuchtenMualem.cpp
►
CreateRelPermNonWettingPhaseVanGenuchtenMualem.h
►
CreateRelPermUdell.cpp
►
CreateRelPermUdell.h
►
CreateRelPermUdellNonwettingPhase.cpp
►
CreateRelPermUdellNonwettingPhase.h
►
CreateRelPermVanGenuchten.cpp
►
CreateRelPermVanGenuchten.h
RelPermBrooksCorey.cpp
►
RelPermBrooksCorey.h
RelPermBrooksCoreyNonwettingPhase.cpp
►
RelPermBrooksCoreyNonwettingPhase.h
RelPermGeneralizedPower.cpp
►
RelPermGeneralizedPower.h
RelPermGeneralizedPowerNonwettingPhase.cpp
►
RelPermGeneralizedPowerNonwettingPhase.h
RelPermLiakopoulos.cpp
►
RelPermLiakopoulos.h
►
RelPermNonWettingPhaseVanGenuchtenMualem.cpp
►
RelPermNonWettingPhaseVanGenuchtenMualem.h
RelPermUdell.cpp
►
RelPermUdell.h
RelPermUdellNonwettingPhase.cpp
►
RelPermUdellNonwettingPhase.h
RelPermVanGenuchten.cpp
►
RelPermVanGenuchten.h
►
SwellingStress
►
ThermalConductivity
►
VapourDiffusion
►
Viscosity
AverageMolarMass.cpp
►
AverageMolarMass.h
BishopsPowerLaw.cpp
►
BishopsPowerLaw.h
BishopsSaturationCutoff.cpp
►
BishopsSaturationCutoff.h
ClausiusClapeyron.cpp
►
ClausiusClapeyron.h
►
Constant.cpp
►
Constant.h
►
CreateAverageMolarMass.cpp
►
CreateAverageMolarMass.h
►
CreateBishopsPowerLaw.cpp
►
CreateBishopsPowerLaw.h
►
CreateBishopsSaturationCutoff.cpp
►
CreateBishopsSaturationCutoff.h
►
CreateClausiusClapeyron.cpp
►
CreateClausiusClapeyron.h
►
CreateConstant.cpp
►
CreateConstant.h
►
CreateCubicLawPermeability.cpp
►
CreateCubicLawPermeability.h
►
CreateCurve.cpp
►
CreateCurve.h
►
CreateDupuitPermeability.cpp
►
CreateDupuitPermeability.h
►
CreateEffectiveThermalConductivityPorosityMixing.cpp
►
CreateEffectiveThermalConductivityPorosityMixing.h
►
CreateEmbeddedFracturePermeability.cpp
►
CreateEmbeddedFracturePermeability.h
►
CreateExponential.cpp
►
CreateExponential.h
►
CreateFunction.cpp
►
CreateFunction.h
►
CreateGasPressureDependentPermeability.cpp
►
CreateGasPressureDependentPermeability.h
►
CreateIdealGasLaw.cpp
►
CreateIdealGasLaw.h
►
CreateIdealGasLawBinaryMixture.cpp
►
CreateIdealGasLawBinaryMixture.h
►
CreateKozenyCarmanModel.cpp
►
CreateKozenyCarmanModel.h
►
CreateLinear.cpp
►
CreateLinear.h
►
CreateOrthotropicEmbeddedFracturePermeability.cpp
►
CreateOrthotropicEmbeddedFracturePermeability.h
►
CreateParameter.cpp
►
CreateParameter.h
►
CreatePengRobinson.cpp
►
CreatePengRobinson.h
►
CreatePermeabilityMohrCoulombFailureIndexModel.cpp
►
CreatePermeabilityMohrCoulombFailureIndexModel.h
►
CreatePermeabilityOrthotropicPowerLaw.cpp
►
CreatePermeabilityOrthotropicPowerLaw.h
►
CreatePorosityFromMassBalance.cpp
►
CreatePorosityFromMassBalance.h
CreateProperties.h
►
CreateSaturationDependentSwelling.cpp
►
CreateSaturationDependentSwelling.h
►
CreateSaturationDependentThermalConductivity.cpp
►
CreateSaturationDependentThermalConductivity.h
►
CreateSpecificHeatCapacityWithLatentHeat.cpp
►
CreateSpecificHeatCapacityWithLatentHeat.h
►
CreateStrainDependentPermeability.cpp
►
CreateStrainDependentPermeability.h
►
CreateTemperatureDependentDiffusion.cpp
►
CreateTemperatureDependentDiffusion.h
►
CreateTemperatureDependentFraction.cpp
►
CreateTemperatureDependentFraction.h
►
CreateTransportPorosityFromMassBalance.cpp
►
CreateTransportPorosityFromMassBalance.h
►
CreateVermaPruessModel.cpp
►
CreateVermaPruessModel.h
►
CreateVolumeFractionAverage.cpp
►
CreateVolumeFractionAverage.h
►
CreateWaterSaturationTemperatureIAPWSIF97Region4.cpp
►
CreateWaterSaturationTemperatureIAPWSIF97Region4.h
►
CreateWaterTemperatureIAPWSIF97Region1.cpp
►
CreateWaterTemperatureIAPWSIF97Region1.h
CubicLawPermeability.cpp
►
CubicLawPermeability.h
Curve.cpp
►
Curve.h
DupuitPermeability.cpp
►
DupuitPermeability.h
EffectiveThermalConductivityPorosityMixing.cpp
►
EffectiveThermalConductivityPorosityMixing.h
EmbeddedFracturePermeability.cpp
►
EmbeddedFracturePermeability.h
Exponential.cpp
►
Exponential.h
►
Function.cpp
►
Function.h
GasPressureDependentPermeability.cpp
►
GasPressureDependentPermeability.h
IdealGasLaw.cpp
►
IdealGasLaw.h
IdealGasLawBinaryMixture.cpp
►
IdealGasLawBinaryMixture.h
KozenyCarmanModel.cpp
►
KozenyCarmanModel.h
Linear.cpp
►
Linear.h
OrthotropicEmbeddedFracturePermeability.cpp
►
OrthotropicEmbeddedFracturePermeability.h
Parameter.cpp
►
Parameter.h
PengRobinson.cpp
►
PengRobinson.h
PermeabilityMohrCoulombFailureIndexModel.cpp
►
PermeabilityMohrCoulombFailureIndexModel.h
PermeabilityOrthotropicPowerLaw.cpp
►
PermeabilityOrthotropicPowerLaw.h
PorosityFromMassBalance.cpp
►
PorosityFromMassBalance.h
Properties.h
SaturationDependentSwelling.cpp
►
SaturationDependentSwelling.h
SpecificHeatCapacityWithLatentHeat.cpp
►
SpecificHeatCapacityWithLatentHeat.h
StrainDependentPermeability.cpp
►
StrainDependentPermeability.h
TemperatureDependentDiffusion.cpp
►
TemperatureDependentDiffusion.h
TemperatureDependentFraction.cpp
►
TemperatureDependentFraction.h
TransportPorosityFromMassBalance.cpp
►
TransportPorosityFromMassBalance.h
VermaPruessModel.cpp
►
VermaPruessModel.h
VolumeFractionAverage.cpp
►
VolumeFractionAverage.h
►
WaterSaturationCurveIAPWSIF97Region4.h
WaterSaturationTemperatureIAPWSIF97Region4.cpp
►
WaterSaturationTemperatureIAPWSIF97Region4.h
►
WaterTemperatureIAPWSIF97Region1.cpp
►
WaterTemperatureIAPWSIF97Region1.h
►
Utils
►
CheckMaterialSpatialDistributionMap.h
►
Component.cpp
►
Component.h
►
CreateComponent.cpp
►
CreateComponent.h
►
CreateMaterialSpatialDistributionMap.cpp
►
CreateMaterialSpatialDistributionMap.h
►
CreateMedium.cpp
►
CreateMedium.h
►
CreatePhase.cpp
►
CreatePhase.h
►
CreateProperty.cpp
►
CreateProperty.h
MaterialSpatialDistributionMap.cpp
►
MaterialSpatialDistributionMap.h
►
Medium.cpp
►
Medium.h
►
Phase.cpp
►
Phase.h
►
Property.cpp
►
Property.h
►
PropertyType.cpp
►
PropertyType.h
►
VariableType.cpp
►
VariableType.h
►
PorousMedium
►
SolidModels
►
Utils
►
PhysicalConstant.h
►
MathLib
►
MeshGeoToolsLib
►
MeshLib
►
MeshToolsLib
►
NumLib
►
ParameterLib
►
ProcessLib
►
File Members
RelPermNonWettingPhaseVanGenuchtenMualem.cpp
Go to the documentation of this file.
1
12
#include "
RelPermNonWettingPhaseVanGenuchtenMualem.h
"
13
14
#include <algorithm>
15
#include <cmath>
16
#include <limits>
17
18
#include "
BaseLib/Error.h
"
19
#include "
MaterialLib/MPL/Medium.h
"
20
#include "
MaterialLib/MPL/Utils/CheckVanGenuchtenExponentRange.h
"
21
#include "
MathLib/Nonlinear/Root1D.h
"
22
23
namespace
MaterialPropertyLib
24
{
25
double
computeVanGenuchtenMualemValue
(
const
double
S_L,
const
double
S_L_r,
26
const
double
S_L_max,
const
double
m)
27
{
28
const
double
Se = (S_L - S_L_r) / (S_L_max - S_L_r);
29
return
std::sqrt(1.0 - Se) * std::pow(1.0 - std::pow(Se, 1.0 / m), 2.0 * m);
30
}
25
double
computeVanGenuchtenMualemValue
(
const
double
S_L,
const
double
S_L_r, {
…
}
31
32
RelPermNonWettingPhaseVanGenuchtenMualem::
33
RelPermNonWettingPhaseVanGenuchtenMualem
(std::string
name
,
34
const
double
S_L_r,
35
const
double
S_n_r,
36
const
double
m,
37
const
double
krel_min,
38
const
double
a)
39
: S_L_r_(S_L_r),
40
S_L_max_(1. - S_n_r),
41
m_(m),
42
krel_min_(krel_min),
43
S_L_for_krel_min_(computeSaturationForMinimumRelativePermeability()),
44
a_(a)
45
{
46
name_
= std::move(
name
);
47
checkVanGenuchtenExponentRange
(
m_
);
48
}
32
RelPermNonWettingPhaseVanGenuchtenMualem::
{
…
}
49
50
PropertyDataType
RelPermNonWettingPhaseVanGenuchtenMualem::value
(
51
VariableArray
const
& variable_array,
52
ParameterLib::SpatialPosition
const
&
/*pos*/
,
double
const
/*t*/
,
53
double
const
/*dt*/
)
const
54
{
55
const
double
S_L =
56
std::clamp(variable_array.
liquid_saturation
,
S_L_r_
,
S_L_max_
);
57
58
const
double
krel =
59
computeVanGenuchtenMualemValue
(S_L,
S_L_r_
,
S_L_max_
,
m_
);
60
return
std::max(
krel_min_
,
a_
* krel);
61
}
50
PropertyDataType
RelPermNonWettingPhaseVanGenuchtenMualem::value
( {
…
}
62
63
PropertyDataType
RelPermNonWettingPhaseVanGenuchtenMualem::dValue
(
64
VariableArray
const
& variable_array,
Variable
const
variable,
65
ParameterLib::SpatialPosition
const
&
/*pos*/
,
double
const
/*t*/
,
66
double
const
/*dt*/
)
const
67
{
68
if
(variable !=
Variable::liquid_saturation
)
69
{
70
OGS_FATAL
(
71
"RelPermNonWettingPhaseVanGenuchtenMualem::dValue is implemented "
72
"for the derivative with respect to liquid saturation only."
);
73
}
74
75
const
double
S_L = variable_array.
liquid_saturation
;
76
if
(S_L < S_L_r_ || S_L >
S_L_for_krel_min_
)
77
{
78
return
0.0;
79
}
80
81
if
(std::fabs(S_L -
S_L_max_
) < std::numeric_limits<double>::epsilon())
82
{
83
return
0.0;
84
}
85
86
const
double
Se = (S_L -
S_L_r_
) / (
S_L_max_
-
S_L_r_
);
87
88
const
double
val1 = std::sqrt(1.0 - Se);
89
const
double
val2 = 1.0 - std::pow(Se, 1.0 /
m_
);
90
91
return
a_
*
92
(-0.5 * std::pow(val2, 2.0 *
m_
) / val1 -
93
2.0 * std::pow(Se, -1.0 + 1.0 /
m_
) * val1 *
94
std::pow(val2, 2.0 *
m_
- 1.0)) /
95
(
S_L_max_
-
S_L_r_
);
96
}
63
PropertyDataType
RelPermNonWettingPhaseVanGenuchtenMualem::dValue
( {
…
}
97
98
double
RelPermNonWettingPhaseVanGenuchtenMualem::
99
computeSaturationForMinimumRelativePermeability
()
const
100
{
101
auto
f = [
this
](
double
S_L) ->
double
102
{
103
return
computeVanGenuchtenMualemValue
(S_L, this->
S_L_r_
, this->
S_L_max_
,
104
this->
m_
) -
105
this->
krel_min_
;
106
};
107
108
auto
root_finder =
109
MathLib::Nonlinear::makeRegulaFalsi<MathLib::Nonlinear::Pegasus>
(
110
f,
S_L_r_
,
S_L_max_
);
111
112
root_finder.step(1000);
113
114
return
root_finder.getResult();
115
}
98
double
RelPermNonWettingPhaseVanGenuchtenMualem::
{
…
}
116
}
// namespace MaterialPropertyLib
CheckVanGenuchtenExponentRange.h
Error.h
OGS_FATAL
#define OGS_FATAL(...)
Definition
Error.h:26
Medium.h
RelPermNonWettingPhaseVanGenuchtenMualem.h
Root1D.h
MaterialPropertyLib::Property::value
virtual PropertyDataType value() const
Definition
Property.cpp:76
MaterialPropertyLib::Property::name_
std::string name_
Definition
Property.h:290
MaterialPropertyLib::RelPermNonWettingPhaseVanGenuchtenMualem::S_L_for_krel_min_
const double S_L_for_krel_min_
Liquid saturation that gives krel_min_.
Definition
RelPermNonWettingPhaseVanGenuchtenMualem.h:102
MaterialPropertyLib::RelPermNonWettingPhaseVanGenuchtenMualem::m_
const double m_
Exponent .
Definition
RelPermNonWettingPhaseVanGenuchtenMualem.h:99
MaterialPropertyLib::RelPermNonWettingPhaseVanGenuchtenMualem::S_L_r_
const double S_L_r_
Residual saturation of wetting phase.
Definition
RelPermNonWettingPhaseVanGenuchtenMualem.h:97
MaterialPropertyLib::RelPermNonWettingPhaseVanGenuchtenMualem::computeSaturationForMinimumRelativePermeability
double computeSaturationForMinimumRelativePermeability() const
Definition
RelPermNonWettingPhaseVanGenuchtenMualem.cpp:99
MaterialPropertyLib::RelPermNonWettingPhaseVanGenuchtenMualem::RelPermNonWettingPhaseVanGenuchtenMualem
RelPermNonWettingPhaseVanGenuchtenMualem(std::string name, const double S_L_r, const double S_n_r, const double m, const double krel_min, const double a)
Definition
RelPermNonWettingPhaseVanGenuchtenMualem.cpp:33
MaterialPropertyLib::RelPermNonWettingPhaseVanGenuchtenMualem::a_
const double a_
Definition
RelPermNonWettingPhaseVanGenuchtenMualem.h:103
MaterialPropertyLib::RelPermNonWettingPhaseVanGenuchtenMualem::S_L_max_
const double S_L_max_
Maximum saturation of wetting phase.
Definition
RelPermNonWettingPhaseVanGenuchtenMualem.h:98
MaterialPropertyLib::RelPermNonWettingPhaseVanGenuchtenMualem::dValue
PropertyDataType dValue(VariableArray const &variable_array, Variable const variable, ParameterLib::SpatialPosition const &pos, double const t, double const dt) const override
Definition
RelPermNonWettingPhaseVanGenuchtenMualem.cpp:63
MaterialPropertyLib::RelPermNonWettingPhaseVanGenuchtenMualem::krel_min_
const double krel_min_
Minimum relative permeability.
Definition
RelPermNonWettingPhaseVanGenuchtenMualem.h:100
MaterialPropertyLib::VariableArray
Definition
VariableType.h:101
MaterialPropertyLib::VariableArray::liquid_saturation
double liquid_saturation
Definition
VariableType.h:184
ParameterLib::SpatialPosition
Definition
SpatialPosition.h:28
MaterialPropertyLib
Definition
ChemicalSolverInterface.h:21
MaterialPropertyLib::computeVanGenuchtenMualemValue
double computeVanGenuchtenMualemValue(const double S_L, const double S_L_r, const double S_L_max, const double m)
Definition
RelPermNonWettingPhaseVanGenuchtenMualem.cpp:25
MaterialPropertyLib::Variable
Variable
Definition
VariableType.h:30
MaterialPropertyLib::Variable::liquid_saturation
@ liquid_saturation
MaterialPropertyLib::name
@ name
Definition
PropertyType.h:66
MaterialPropertyLib::checkVanGenuchtenExponentRange
void checkVanGenuchtenExponentRange(const double m)
Definition
CheckVanGenuchtenExponentRange.cpp:18
MaterialPropertyLib::PropertyDataType
std::variant< double, Eigen::Matrix< double, 2, 1 >, Eigen::Matrix< double, 3, 1 >, Eigen::Matrix< double, 2, 2 >, Eigen::Matrix< double, 3, 3 >, Eigen::Matrix< double, 4, 1 >, Eigen::Matrix< double, 6, 1 >, Eigen::MatrixXd > PropertyDataType
Definition
Property.h:31
MathLib::Nonlinear::makeRegulaFalsi
RegulaFalsi< SubType, Function > makeRegulaFalsi(Function &&f, double const a, double const b)
Definition
Root1D.h:137
MaterialLib
MPL
Properties
RelativePermeability
RelPermNonWettingPhaseVanGenuchtenMualem.cpp
Generated by
1.12.0