22 "simm >= 0 & simm < 1",
23 "Immobile saturation. Must be between 0 and 1. Define s = " 24 "(seff - simm)/(1 - simm). Then relperm = 1 - (n+1)(1-s)^n " 29 "Exponent. Must be >= 2. Define s = (seff - simm)/(1 - simm). Then " 30 "relperm = 1 - (n+1)(1-s)^n + n(1-s)^(n+1)");
31 params.
addClassDescription(
"Power form of relative permeability that might be useful for gases. " 32 "Define s = (seff - simm)/(1 - simm). Then relperm = 1 - " 33 "(n+1)(1-s)^n + n(1-s)^(n+1) if s<simm, otherwise relperm=1");
78 return krelp / (1.0 -
_simm);
RichardsRelPermPowerGas(const InputParameters ¶meters)
Real relperm(Real seff) const
Relative permeability.
static InputParameters validParams()
Base class for Richards relative permeability classes that provide relative permeability as a functio...
Real d2relperm(Real seff) const
Second derivative of relative permeability wrt seff.
static InputParameters validParams()
PowerGas form of relative permeability Define s = (seff - simm)/(1 - simm).
Real _simm
immobile saturation
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
registerMooseObject("RichardsApp", RichardsRelPermPowerGas)
MooseUnits pow(const MooseUnits &, int)
Real drelperm(Real seff) const
Derivative of relative permeability wrt seff.