Real _cut_limit
where the cubic starts
Real d2density(Real p) const
second derivative of fluid density wrt porepressure
static InputParameters validParams()
Real _dens0
density = _dens0*exp(p/_bulk), modified by cubic
Real _bulk
density = _dens0*exp(p/_bulk), modified by cubic
Real density(Real p) const
fluid density as a function of porepressure
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
Real _c3
(cut_limit-zero_point)^3
const InputParameters & parameters() const
Base class for fluid density as a function of porepressure The functions density, ddensity and d2dens...
Fluid density assuming constant bulk modulus, for p>cut_limit Then following a cubic for zero_point <...
Real ddensity(Real p) const
derivative of fluid density wrt porepressure
Real _zero_point
where the density is zero
RichardsDensityConstBulkCut(const InputParameters ¶meters)