40 Real seff(std::vector<const VariableValue *> p,
unsigned int qp)
const;
50 dseff(std::vector<const VariableValue *> p,
unsigned int qp, std::vector<Real> & result)
const;
59 void d2seff(std::vector<const VariableValue *> p,
61 std::vector<std::vector<Real>> & result)
const;
RichardsSeff1BWsmall(const InputParameters ¶meters)
void dseff(std::vector< const VariableValue *> p, unsigned int qp, std::vector< Real > &result) const
derivative of effective saturation as a function of porepressure
Base class for effective saturation as a function of porepressure(s) The functions seff...
"Broadbridge-White" form of effective saturation for Kn small as a function of porepressure (not capi...
void d2seff(std::vector< const VariableValue *> p, unsigned int qp, std::vector< std::vector< Real >> &result) const
second derivative of effective saturation as a function of porepressure
Real _las
BW's lambda_s parameter multiplied by (fluiddensity*gravity)
static InputParameters validParams()
Real _sn
BW's initial effective saturation.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
Real _c
BW's C parameter.
Real LambertW(const double z) const
LambertW function, returned value satisfies W(z)*exp(W(z))=z.
const InputParameters & parameters() const
Real seff(std::vector< const VariableValue *> p, unsigned int qp) const
effective saturation as a function of porepressure
Real _ss
Effective saturation where porepressure = 0.