19 "sat_lr >= 0 & sat_lr < 1",
20 "Liquid residual saturation. Must be between 0 and 1. Default is 0");
24 "Maximum capillary pressure (Pa). Must be > 0. Default is 1e9");
25 params.
addParam<
bool>(
"log_extension",
27 "Use a logarithmic extension for low saturation to avoid capillary " 28 "pressure going to infinity. Default is true. Set to false if your " 29 "capillary pressure depends on spatially-dependent variables other than " 30 "saturation, as the log-extension C++ code for this case has yet to be " 38 _sat_lr(getParam<
Real>(
"sat_lr")),
39 _dseff_ds(1.0 / (1.0 - _sat_lr)),
40 _log_ext(getParam<bool>(
"log_extension")),
41 _pc_max(getParam<
Real>(
"pc_max")),
78 result.derivatives() =
saturation.derivatives() * dPc_ds;
99 result.derivatives() =
saturation.derivatives() * d2Pc_ds2;
132 result.derivatives() = pc.derivatives() * ds_dpc;
150 result.derivatives() = pc.derivatives() * d2s_dpc2;
188 const unsigned int max_its = 20;
189 const Real nr_tol = 1.0e-8;
190 unsigned int iter = 0;
virtual Real dCapillaryPressureCurve(Real saturation, unsigned qp=0) const =0
Derivative of raw capillary pressure wrt true saturation.
virtual Real effectiveSaturation(Real pc, unsigned qp=0) const =0
Effective saturation as a function of capillary pressure.
virtual Real d2EffectiveSaturation(Real pc, unsigned qp=0) const =0
Second derivative of effective saturation wrt capillary pressure.
virtual void initialSetup() override
static InputParameters validParams()
Real effectiveSaturationFromSaturation(Real saturation) const
Effective saturation of liquid phase given liquid saturation and residual liquid saturation.
Real saturation(Real pc, unsigned qp=0) const
Saturation as a function of capillary pressure.
virtual Real capillaryPressureCurve(Real saturation, unsigned qp=0) const =0
Raw capillary pressure curve (does not include logarithmic extension)
Real _slope_ext
Gradient of the logarithmic extension This is computed only for qp=0.
ADRealEigenVector< T, D, asd > abs(const ADRealEigenVector< T, D, asd > &)
const Real _pc_max
Maximum capillary pressure (Pa). Note: must be <= 0.
Real interceptFunctionDeriv(Real s) const
Calculates the saturation where the logarithmic extension to capillary pressure at low saturation...
Real extensionSaturation() const
Calculates the saturation where the logarithmic extension to capillary pressure meets the raw curve u...
virtual Real capillaryPressure(Real saturation, unsigned qp=0) const
Capillary pressure is calculated as a function of true saturation.
Real dSaturation(Real pc, unsigned qp=0) const
Derivative of saturation wrt capillary pressure.
Real d2CapillaryPressureLogExt(Real s) const
The second derivative of capillary pressure in the logarithmic extension This implementation assumes ...
PorousFlowCapillaryPressure(const InputParameters ¶meters)
Real _pc_ext
Capillary pressure where the extension meets the raw curve.
virtual Real d2CapillaryPressure(Real saturation, unsigned qp=0) const
Second derivative of capillary pressure wrt true saturation.
virtual Real d2CapillaryPressureCurve(Real saturation, unsigned qp=0) const =0
Second derivative of raw capillary pressure wrt true saturation.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
Real interceptFunction(Real s) const
Calculates the saturation where the logarithmic extension to capillary pressure at low saturation...
Real capillaryPressureLogExt(Real s) const
The capillary pressure in the logarithmic extension This implementation assumes capillary-pressure is...
virtual Real dCapillaryPressure(Real saturation, unsigned qp=0) const
Derivative of capillary pressure wrt true saturation.
Real _sat_ext
Saturation where the logarithmic extension meets the raw curve This is computed only for qp=0...
const Real _sat_lr
Liquid residual saturation.
bool _log_ext
Flag to use a logarithmic extension for low saturation.
Real d2Saturation(Real pc, unsigned qp=0) const
Second derivative of saturation wrt capillary pressure.
static InputParameters validParams()
MooseUnits pow(const MooseUnits &, int)
Real dCapillaryPressureLogExt(Real s) const
The derivative of capillary pressure in the logarithmic extension This implementation assumes capilla...
virtual Real dEffectiveSaturation(Real pc, unsigned qp=0) const =0
Derivative of effective saturation wrt capillary pressure.