22 "simm >=0 & simm < 1",
23 "Immobile saturation. Must be between 0 and 1. Define s = " 24 "(seff - simm)/(1 - simm). Then relperm = s^(1/2) * (1 - (1 " 29 "van-Genuchten m parameter. Must be between 0 and 1, and optimally " 30 "should be set >0.5. Define s = (seff - simm)/(1 - simm). Then " 31 "relperm = s^(1/2) * (1 - (1 - s^(1/m))^m)^2");
33 "scut",
"scut > 0 & scut < 1",
"cutoff in effective saturation.");
34 params.
addClassDescription(
"VG1 form of relative permeability. Define s = (seff - simm)/(1 - " 35 "simm). Then relperm = s^(1/2) * (1 - (1 - s^(1/m))^m)^2, if s>0, " 36 "and relperm=0 otherwise");
42 _simm(getParam<
Real>(
"simm")),
43 _m(getParam<
Real>(
"m")),
44 _scut(getParam<
Real>(
"scut")),
60 _console <<
"Relative permeability of VG1 type has cubic coefficients " <<
_vg1_const <<
" " 75 if (s_internal <
_scut)
105 if (s_internal <
_scut)
110 return krelp / (1.0 -
_simm);
124 if (s_internal <
_scut)
Real d2relperm(Real seff) const
second derivative of relative permeability wrt effective saturation
static InputParameters validParams()
registerMooseObject("RichardsApp", RichardsRelPermVG1)
Van-Genuchten form of relative permeability as a function of effective saturation.
Real drelperm(Real seff) const
derivative of relative permeability wrt effective saturation
Real drelperm(Real seff) const
derivative of relative permeability wrt effective saturation
Real _vg1_quad
coefficient of quadratic term in cubic relperm relation
static InputParameters validParams()
Van-Genuchten form of relative permeability when seff <= _scut cubic relative permeability for seff >...
Real _vg1_cub
coefficient of cubic term in cubic relperm relation
RichardsRelPermVG1(const InputParameters ¶meters)
Real _simm
immobile saturation
Real d2relperm(Real seff) const
second derivative of relative permeability wrt effective saturation
Real relperm(Real seff) const
relative permeability as a function of effective saturation
void initialSetup()
just prints some (maybe) useful info to the console
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
Real _vg1_linear
coefficient of linear term in cubic relperm relation
Real _scut
for seff > _scut use cubic relative permeability, otherwise use van Genuchten
const ConsoleStream _console
Real _vg1_const
constant in cubic relperm relation
Real relperm(Real seff) const
relative permeability as a function of effective saturation
MooseUnits pow(const MooseUnits &, int)