17 if (seff <= 0.0 || seff >= 1.0)
19 return (1.0 + m) * m *
std::pow(seff, m - 1.0) - m * (m + 1.0) *
std::pow(seff, m);
25 if (seff <= 0.0 || seff >= 1.0)
27 return (1.0 + m) * m * (m - 1.0) *
std::pow(seff, m - 2.0) -
28 m * (m + 1.0) * m *
std::pow(seff, m - 1.0);
Real d2RelativePermeability(Real seff, Real m)
Second derivative of relative permeability with respect to effective saturation.
Real dRelativePermeability(Real seff, Real m)
Derivative of relative permeability with respect to effective saturation.
FLAC inspired relative permeability relationship.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
MooseUnits pow(const MooseUnits &, int)