22 "simm >= 0 & simm < 1",
23 "Immobile saturation. Must be between 0 and 1. Define s = " 24 "(seff - simm)/(1 - simm). Then relperm = (n+1)s^n - " 28 "Exponent. Must be >= 2. Define s = " 29 "(seff - simm)/(1 - simm). Then " 30 "relperm = (n+1)s^n - ns^(n+1)");
31 params.
addClassDescription(
"Power form of relative permeability. Define s = (seff - simm)/(1 - " 32 "simm). Then relperm = (n+1)s^n - ns^(n+1) if s<simm, otherwise " 78 return krelp / (1.0 -
_simm);
static InputParameters validParams()
Real d2relperm(Real seff) const
second derivative of relative permeability wrt effective saturation
Base class for Richards relative permeability classes that provide relative permeability as a functio...
Real drelperm(Real seff) const
derivative of relative permeability wrt effective saturation
Real relperm(Real seff) const
relative permeability as a function of effective saturation
Power form of relative permeability, usually used for water.
static InputParameters validParams()
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
Real _n
exponent used in the power relationship
registerMooseObject("RichardsApp", RichardsRelPermPower)
MooseUnits pow(const MooseUnits &, int)
Real _simm
immobile saturation
RichardsRelPermPower(const InputParameters ¶meters)