Base material designed to calculate and store quantities relevant for hysteretic capillary pressure calculations. More...

#include <PorousFlowHystereticCapillaryPressure.h>

Inheritance diagram for PorousFlowHystereticCapillaryPressure:

[legend]

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	PorousFlowHystereticCapillaryPressure (const InputParameters &parameters)

const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialProperty (const std::string &name)

const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialPropertyByName (const std::string &name)

void	validateDerivativeMaterialPropertyBase (const std::string &base)

const MaterialPropertyName	derivativePropertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	derivativePropertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	derivativePropertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	derivativePropertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)

void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)

void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual void	initQpStatefulProperties () override

virtual void	computeQpProperties () override

Real	landSat (Real slDel) const

void	computeTurningPointInfo (unsigned tp_num, Real tp_sat, Real tp_pc)
	Compute all relevant quantities at the given turning point. More...

Real	capillaryPressureQp (Real sat) const

Real	dcapillaryPressureQp (Real sat) const

Real	d2capillaryPressureQp (Real sat) const

Real	liquidSaturationQp (Real pc) const

Real	dliquidSaturationQp (Real pc) const

Real	d2liquidSaturationQp (Real pc) const

Protected Attributes
const Real	_alpha_d
	van Genuchten alpha parameter for the primary drying curve More...

const Real	_alpha_w
	van Genuchten alpha parameter for the primary wetting curve More...

const Real	_n_d
	van Genuchten n parameter for the primary drying curve More...

const Real	_n_w
	van Genuchten n parameter for the primary wetting curve More...

const Real	_s_l_min
	Minimum liquid saturation for which the van Genuchten expression is valid (Pc(_s_l_min) = infinity) More...

const Real	_s_lr
	Liquid saturation below which the liquid relative permeability is zero. More...

const Real	_s_gr_max
	Residual gas saturation: 1 - _s_gr_max is the maximum saturation for which the van Genuchten expression is valid for the wetting curve. More...

const Real	_pc_max
	Maximum capillary pressure: for Pc above this value, a "lower" extension will be used. More...

const Real	_high_ratio
	The high-saturation extension to the wetting will commence at _high_ratio * (1 - _s_gr_del) More...

const PorousFlowVanGenuchten::LowCapillaryPressureExtension::ExtensionStrategy	_low_ext_type
	Type of low-saturation extension. More...

const Real	_s_low_d
	Saturation on the primary drying curve where low-saturation extension commences. More...

const Real	_dpc_low_d
	d(Pc)/dS on the primary drying curve at S = _s_low_d More...

const PorousFlowVanGenuchten::LowCapillaryPressureExtension	_low_ext_d
	Parameters involved in the low-saturation extension of the primary drying curve. More...

const Real	_s_low_w
	Saturation on the primary wetting curve where low-saturation extension commences. More...

const Real	_dpc_low_w
	d(Pc)/dS on the primary wetting curve at S = _s_low_w More...

const PorousFlowVanGenuchten::LowCapillaryPressureExtension	_low_ext_w
	Parameters involved in the low-saturation extension of the primary wetting curve. More...

const PorousFlowVanGenuchten::HighCapillaryPressureExtension::ExtensionStrategy	_high_ext_type
	Type of high-saturation extension of the wetting curves. More...

const Real	_s_high
	Saturation at the point of high-saturation extension. More...

const Real	_pc_high
	Pc at the point of high-saturation extension. More...

const Real	_dpc_high
	d(Pc)/d(S) at the point of high-saturation extension More...

const PorousFlowVanGenuchten::HighCapillaryPressureExtension	_high_ext
	Parameters involved in the high-saturation extension of the primary wetting curve. More...

const MaterialProperty< unsigned > &	_hys_order
	Hysteresis order, as computed by PorousFlowHysteresisOrder. More...

const MaterialProperty< unsigned > &	_hys_order_old
	Old value of hysteresis order, as computed by PorousFlowHysteresisOrder. More...

const MaterialProperty< std::array< Real, PorousFlowConstants::MAX_HYSTERESIS_ORDER > > &	_hys_sat_tps
	Saturation values at the turning points, as computed by PorousFlowHysteresisOrder. More...

const MaterialProperty< Real > &	_pc_older
	Older value of capillary pressure. More...

MaterialProperty< std::array< Real, PorousFlowConstants::MAX_HYSTERESIS_ORDER > > &	_pc_tps
	Nodal or quadpoint values of Pc at the turning points. More...

MaterialProperty< std::array< Real, PorousFlowConstants::MAX_HYSTERESIS_ORDER > > &	_s_d_tps
	Computed nodal or quadpoint values of saturation on the drying curve at _pc_tps. More...

MaterialProperty< std::array< Real, PorousFlowConstants::MAX_HYSTERESIS_ORDER > > &	_s_gr_tps
	Computed nodal or quadpoint values of S_gr_Del, ie, the Land expression, at the turning points. More...

MaterialProperty< std::array< PorousFlowVanGenuchten::LowCapillaryPressureExtension, PorousFlowConstants::MAX_HYSTERESIS_ORDER > > &	_w_low_ext_tps
	Nodal or quadpoint values of the low extension of the wetting curve defined by _s_gr_tps. More...

MaterialProperty< std::array< PorousFlowVanGenuchten::HighCapillaryPressureExtension, PorousFlowConstants::MAX_HYSTERESIS_ORDER > > &	_w_high_ext_tps
	Nodal or quadpoint values of the high extension of the wetting curve defined by _s_gr_tps. More...

MaterialProperty< std::array< Real, PorousFlowConstants::MAX_HYSTERESIS_ORDER > > &	_s_w_tps
	Computed nodal or quadpoint values of liquid saturation on the wetting curve defined by _s_gr_del, at pc = _pc_d_tps. More...

const unsigned int	_num_phases
	Number of phases. More...

const unsigned int	_num_components
	Number of components. More...

const unsigned int	_num_pf_vars
	Number of PorousFlow variables. More...

GenericMaterialProperty< std::vector< Real >, is_ad > &	_porepressure
	Computed nodal or quadpoint values of porepressure of the phases. More...

MaterialProperty< std::vector< std::vector< Real > > > *const	_dporepressure_dvar
	d(porepressure)/d(PorousFlow variable) More...

GenericMaterialProperty< std::vector< RealGradient >, is_ad > *const	_gradp_qp
	Grad(p) at the quadpoints. More...

MaterialProperty< std::vector< std::vector< Real > > > *const	_dgradp_qp_dgradv
	d(grad porepressure)/d(grad PorousFlow variable) at the quadpoints More...

MaterialProperty< std::vector< std::vector< RealGradient > > > *const	_dgradp_qp_dv
	d(grad porepressure)/d(PorousFlow variable) at the quadpoints More...

GenericMaterialProperty< std::vector< Real >, is_ad > &	_saturation
	Computed nodal or qp saturation of the phases. More...

MaterialProperty< std::vector< std::vector< Real > > > *const	_dsaturation_dvar
	d(saturation)/d(PorousFlow variable) More...

GenericMaterialProperty< std::vector< RealGradient >, is_ad > *const	_grads_qp
	Grad(s) at the quadpoints. More...

MaterialProperty< std::vector< std::vector< Real > > > *const	_dgrads_qp_dgradv
	d(grad saturation)/d(grad PorousFlow variable) at the quadpoints More...

MaterialProperty< std::vector< std::vector< RealGradient > > > *const	_dgrads_qp_dv
	d(grad saturation)/d(PorousFlow variable) at the quadpoints More...

Private Member Functions
Real	firstOrderWettingPc (Real sat) const

Real	dfirstOrderWettingPc (Real sat) const

Real	d2firstOrderWettingPc (Real sat) const

Real	secondOrderDryingPc (Real sat) const

Real	dsecondOrderDryingPc (Real sat) const

Real	d2secondOrderDryingPc (Real sat) const

Real	firstOrderWettingSat (Real pc) const

Real	dfirstOrderWettingSat (Real pc) const

Real	d2firstOrderWettingSat (Real pc) const

Real	secondOrderDryingSat (Real pc) const

Real	dsecondOrderDryingSat (Real pc) const

Real	d2secondOrderDryingSat (Real pc) const

Detailed Description

Base material designed to calculate and store quantities relevant for hysteretic capillary pressure calculations.

Definition at line 20 of file PorousFlowHystereticCapillaryPressure.h.

Constructor & Destructor Documentation

◆ PorousFlowHystereticCapillaryPressure()

PorousFlowHystereticCapillaryPressure::PorousFlowHystereticCapillaryPressure ( const InputParameters & parameters )

Definition at line 99 of file PorousFlowHystereticCapillaryPressure.C.

   : PorousFlowVariableBase(parameters),
     _alpha_d(getParam<Real>("alpha_d")),
     _alpha_w(getParam<Real>("alpha_w")),
     _n_d(getParam<Real>("n_d")),
     _n_w(getParam<Real>("n_w")),
     _s_l_min(getParam<Real>("S_l_min")),
     _s_lr(getParam<Real>("S_lr")),
     _s_gr_max(getParam<Real>("S_gr_max")),
     _pc_max(getParam<Real>("Pc_max")),
     _high_ratio(getParam<Real>("high_ratio")),
     _low_ext_type(
         getParam<MooseEnum>("low_extension_type")
             .getEnum<PorousFlowVanGenuchten::LowCapillaryPressureExtension::ExtensionStrategy>()),
     _s_low_d(PorousFlowVanGenuchten::saturationHys(_pc_max, _s_l_min, 0.0, _alpha_d, _n_d)),
     _dpc_low_d(
         PorousFlowVanGenuchten::dcapillaryPressureHys(_s_low_d, _s_l_min, 0.0, _alpha_d, _n_d)),
     _low_ext_d(_low_ext_type, _s_low_d, _pc_max, _dpc_low_d),
     _s_low_w(PorousFlowVanGenuchten::saturationHys(_pc_max, _s_l_min, _s_gr_max, _alpha_w, _n_w)),
     _dpc_low_w(PorousFlowVanGenuchten::dcapillaryPressureHys(
         _s_low_w, _s_l_min, _s_gr_max, _alpha_w, _n_w)),
     _low_ext_w(_low_ext_type, _s_low_w, _pc_max, _dpc_low_w),
     _high_ext_type(
         getParam<MooseEnum>("high_extension_type")
             .getEnum<PorousFlowVanGenuchten::HighCapillaryPressureExtension::ExtensionStrategy>()),
     _s_high(_high_ratio * (1 - _s_gr_max)),
     _pc_high(
         PorousFlowVanGenuchten::capillaryPressureHys(_s_high, _s_l_min, _s_gr_max, _alpha_w, _n_w)),
     _dpc_high(PorousFlowVanGenuchten::dcapillaryPressureHys(
         _s_high, _s_l_min, _s_gr_max, _alpha_w, _n_w)),
     _high_ext(_high_ext_type, _s_high, _pc_high, _dpc_high),
     _hys_order(_nodal_material ? getMaterialProperty<unsigned>("PorousFlow_hysteresis_order_nodal")
                                : getMaterialProperty<unsigned>("PorousFlow_hysteresis_order_qp")),
     _hys_order_old(_nodal_material
                        ? getMaterialPropertyOld<unsigned>("PorousFlow_hysteresis_order_nodal")
                        : getMaterialPropertyOld<unsigned>("PorousFlow_hysteresis_order_qp")),
     _hys_sat_tps(
         _nodal_material
             ? getMaterialProperty<std::array<Real, PorousFlowConstants::MAX_HYSTERESIS_ORDER>>(
                   "PorousFlow_hysteresis_saturation_tps_nodal")
             : getMaterialProperty<std::array<Real, PorousFlowConstants::MAX_HYSTERESIS_ORDER>>(
                   "PorousFlow_hysteresis_saturation_tps_qp")),
     _pc_older(_nodal_material
                   ? getMaterialPropertyOlder<Real>("PorousFlow_hysteretic_capillary_pressure_nodal")
                   : getMaterialPropertyOlder<Real>("PorousFlow_hysteretic_capillary_pressure_qp")),
     _pc_tps(_nodal_material
                 ? declareProperty<std::array<Real, PorousFlowConstants::MAX_HYSTERESIS_ORDER>>(
                       "PorousFlow_hysteresis_Pc_tps_nodal")
                 : declareProperty<std::array<Real, PorousFlowConstants::MAX_HYSTERESIS_ORDER>>(
                       "PorousFlow_hysteresis_Pc_tps_qp")),
     _s_d_tps(_nodal_material
                  ? declareProperty<std::array<Real, PorousFlowConstants::MAX_HYSTERESIS_ORDER>>(
                        "PorousFlow_hysteresis_s_d_tps_nodal")
                  : declareProperty<std::array<Real, PorousFlowConstants::MAX_HYSTERESIS_ORDER>>(
                        "PorousFlow_hysteresis_s_d_tps_qp")),
     _s_gr_tps(_nodal_material
                   ? declareProperty<std::array<Real, PorousFlowConstants::MAX_HYSTERESIS_ORDER>>(
                         "PorousFlow_hysteresis_s_gr_tps_nodal")
                   : declareProperty<std::array<Real, PorousFlowConstants::MAX_HYSTERESIS_ORDER>>(
                         "PorousFlow_hysteresis_s_gr_tps_qp")),
     _w_low_ext_tps(
         _nodal_material
             ? declareProperty<std::array<PorousFlowVanGenuchten::LowCapillaryPressureExtension,
                                          PorousFlowConstants::MAX_HYSTERESIS_ORDER>>(
                   "PorousFlow_hysteresis_w_low_ext_tps_nodal")
             : declareProperty<std::array<PorousFlowVanGenuchten::LowCapillaryPressureExtension,
                                          PorousFlowConstants::MAX_HYSTERESIS_ORDER>>(
                   "PorousFlow_hysteresis_w_low_ext_tps_qp")),
     _w_high_ext_tps(
         _nodal_material
             ? declareProperty<std::array<PorousFlowVanGenuchten::HighCapillaryPressureExtension,
                                          PorousFlowConstants::MAX_HYSTERESIS_ORDER>>(
                   "PorousFlow_hysteresis_w_high_ext_tps_nodal")
             : declareProperty<std::array<PorousFlowVanGenuchten::HighCapillaryPressureExtension,
                                          PorousFlowConstants::MAX_HYSTERESIS_ORDER>>(
                   "PorousFlow_hysteresis_w_high_ext_tps_qp")),
     _s_w_tps(_nodal_material
                  ? declareProperty<std::array<Real, PorousFlowConstants::MAX_HYSTERESIS_ORDER>>(
                        "PorousFlow_hysteresis_s_w_tps_nodal")
                  : declareProperty<std::array<Real, PorousFlowConstants::MAX_HYSTERESIS_ORDER>>(
                        "PorousFlow_hysteresis_s_w_tps_qp"))
 {
   if (_s_gr_max >= 1 - _s_l_min)
     paramError("S_gr_max", "Must be less than 1 - S_l_min");
   if (_s_high < _s_low_w)
     paramError("high_ratio",
                "Should be chosen sufficiently close to 1 so that the high extension does not "
                "interfere with the low extension.  Instead, you may have chosen Pc_max too low");
   if (_s_lr <= _s_l_min)
     mooseWarning("S_lr should usually be greater than S_l_min");
 }

Member Function Documentation

◆ capillaryPressureQp()

Real PorousFlowHystereticCapillaryPressure::capillaryPressureQp ( Real sat ) const

protected

Returns: the value of capillary pressure, given the liquid saturation. This uses _hys_order[_qp]

Parameters

sat	liquid saturation

Definition at line 300 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by PorousFlow2PhaseHysPS::buildQpPPSS(), PorousFlowHystereticInfo::computeQpInfo(), PorousFlowHystereticInfo::computeQpProperties(), and PorousFlowHystereticInfo::initQpStatefulProperties().

 {
   Real pc = 0.0;
   if (_hys_order[_qp] == 0) // on primary drying curve
     pc = PorousFlowVanGenuchten::capillaryPressureHys(
         sat, _s_l_min, 0.0, _alpha_d, _n_d, _low_ext_d);
   else if (_hys_order[_qp] == 1) // first-order wetting
     pc = firstOrderWettingPc(sat);
   else if (_hys_order[_qp] == 2) // second-order drying
     pc = secondOrderDryingPc(sat);
   else // third order drying and wetting
   {
     const Real tp1 = _hys_sat_tps[_qp].at(1);
     const Real tp2 = _hys_sat_tps[_qp].at(2);
     const Real pc1 = firstOrderWettingPc(sat);
     const Real pc2 = secondOrderDryingPc(sat);
     // handle cases that occur just at the transition from 3rd to 2nd order, or 3rd to 1st order
     if (sat < tp2)
       pc = pc2;
     else if (sat > tp1)
       pc = pc1;
     else if (pc1 >= pc2)
       pc = pc1;
     else if (pc1 <= 0.0 || pc2 <= 0.0)
       pc = 0.0;
     else
       pc = std::exp(std::log(pc1) + (sat - tp1) * (std::log(pc2) - std::log(pc1)) / (tp2 - tp1));
   }
   return pc;
 }

◆ computeQpProperties()

void PorousFlowHystereticCapillaryPressure::computeQpProperties ( )

overrideprotectedvirtual

Reimplemented from PorousFlowVariableBaseTempl< is_ad >.

Reimplemented in PorousFlow1PhaseHysP, PorousFlow2PhaseHysPP, PorousFlow2PhaseHysPS, and PorousFlowHystereticInfo.

Definition at line 220 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by PorousFlowHystereticInfo::computeQpProperties(), PorousFlow2PhaseHysPP::computeQpProperties(), PorousFlow2PhaseHysPS::computeQpProperties(), and PorousFlow1PhaseHysP::computeQpProperties().

 {
   // size stuff correctly and prepare the derivative matrices with zeroes
   PorousFlowVariableBase::computeQpProperties();
 
   if (_hys_order[_qp] != _hys_order_old[_qp] && _hys_order[_qp] > 0)
   {
     const unsigned tp_num = _hys_order[_qp] - 1;
     computeTurningPointInfo(tp_num, _hys_sat_tps[_qp].at(tp_num), _pc_older[_qp]);
   }
 }

◆ computeTurningPointInfo()

void PorousFlowHystereticCapillaryPressure::computeTurningPointInfo	(	unsigned	tp_num,
		Real	tp_sat,
		Real	tp_pc
	)

protected

Compute all relevant quantities at the given turning point.

Parameters

tp_num	The turning point number. Eg, tp_num=0 upon transition from zeroth-order drying to first-order wetting
tp_sat	Liquid saturation at the turning point
tp_pc	Capillary pressure at the turning point

Definition at line 240 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by computeQpProperties(), and initQpStatefulProperties().

 {
   _pc_tps[_qp].at(tp_num) = tp_pc;
 
   // Quantities on the drying curve:
   // pc on the drying curve at the turning point saturation
   const Real pc_d_tps = PorousFlowVanGenuchten::capillaryPressureHys(
       tp_sat, _s_l_min, 0.0, _alpha_d, _n_d, _low_ext_d);
   // saturation on the drying curve, at tp_pc
   _s_d_tps[_qp].at(tp_num) =
       PorousFlowVanGenuchten::saturationHys(tp_pc, _s_l_min, 0.0, _alpha_d, _n_d, _low_ext_d);
 
   // Quantities relevant to the wetting curve defined by the Land expression.
   // s_gr_tps is the Land expression as a function of the turning point saturation
   _s_gr_tps[_qp].at(tp_num) = landSat(tp_sat);
   // the low extension of the wetting curve defined by _s_gr_tps
   const Real s_w_low_ext = PorousFlowVanGenuchten::saturationHys(
       _pc_max, _s_l_min, _s_gr_tps[_qp].at(tp_num), _alpha_w, _n_w);
   const Real dpc_w_low_ext = PorousFlowVanGenuchten::dcapillaryPressureHys(
       s_w_low_ext, _s_l_min, _s_gr_tps[_qp].at(tp_num), _alpha_w, _n_w);
   _w_low_ext_tps[_qp].at(tp_num) = PorousFlowVanGenuchten::LowCapillaryPressureExtension(
       _low_ext_type, s_w_low_ext, _pc_max, dpc_w_low_ext);
   // the high extension of the wetting curve defined by _s_gr_tps
   const Real s_w_high_ext =
       (_high_ext_type == PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE)
           ? 1.0 - _s_gr_tps[_qp].at(tp_num)
           : _high_ratio *
                 (1.0 -
                  _s_gr_tps[_qp].at(
                      tp_num)); // if NONE then use the vanGenuchten all the way to 1 - _s_gr_tps
   const Real pc_w_high_ext =
       PorousFlowVanGenuchten::capillaryPressureHys(s_w_high_ext,
                                                    _s_l_min,
                                                    _s_gr_tps[_qp].at(tp_num),
                                                    _alpha_w,
                                                    _n_w,
                                                    _w_low_ext_tps[_qp].at(tp_num));
   const Real dpc_w_high_ext =
       PorousFlowVanGenuchten::dcapillaryPressureHys(s_w_high_ext,
                                                     _s_l_min,
                                                     _s_gr_tps[_qp].at(tp_num),
                                                     _alpha_w,
                                                     _n_w,
                                                     _w_low_ext_tps[_qp].at(tp_num));
   _w_high_ext_tps[_qp].at(tp_num) = PorousFlowVanGenuchten::HighCapillaryPressureExtension(
       _high_ext_type, s_w_high_ext, pc_w_high_ext, dpc_w_high_ext);
 
   // saturation on the wetting curve defined by _s_gr_tps, at pc = pc_d_tps
   _s_w_tps[_qp].at(tp_num) = PorousFlowVanGenuchten::saturationHys(pc_d_tps,
                                                                    _s_l_min,
                                                                    _s_gr_tps[_qp].at(tp_num),
                                                                    _alpha_w,
                                                                    _n_w,
                                                                    _w_low_ext_tps[_qp].at(tp_num),
                                                                    _w_high_ext_tps[_qp].at(tp_num));
 }

◆ d2capillaryPressureQp()

Real PorousFlowHystereticCapillaryPressure::d2capillaryPressureQp ( Real sat ) const

protected

Returns: d^2(capillary pressure)/d(sat)^2. This uses _hys_order[_qp]

Parameters

sat	liquid saturation

Definition at line 371 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by PorousFlowHystereticInfo::computeQpInfo(), and PorousFlow2PhaseHysPS::computeQpProperties().

 {
   Real d2pc = 0.0;
   if (_hys_order[_qp] == 0) // on primary drying curve
     d2pc = PorousFlowVanGenuchten::d2capillaryPressureHys(
         sat, _s_l_min, 0.0, _alpha_d, _n_d, _low_ext_d);
   else if (_hys_order[_qp] == 1) // first-order wetting
     d2pc = d2firstOrderWettingPc(sat);
   else if (_hys_order[_qp] == 2) // second-order drying
     d2pc = d2secondOrderDryingPc(sat);
   else // third order drying and wetting
   {
     const Real tp1 = _hys_sat_tps[_qp].at(1);
     const Real tp2 = _hys_sat_tps[_qp].at(2);
     const Real pc1 = firstOrderWettingPc(sat);
     const Real pc2 = secondOrderDryingPc(sat);
     // handle cases that occur just at the transition from 3rd to 2nd order, or 3rd to 1st order
     if (sat < tp2)
       d2pc = d2secondOrderDryingPc(sat);
     else if (sat > tp1)
       d2pc = d2firstOrderWettingPc(sat);
     else if (pc1 >= pc2)
       d2pc = d2firstOrderWettingPc(sat);
     else if (pc1 <= 0.0 || pc2 <= 0.0)
       d2pc = 0.0;
     else
     {
       const Real pc =
           std::exp(std::log(pc1) + (sat - tp1) * (std::log(pc2) - std::log(pc1)) / (tp2 - tp1));
       const Real dpc1 = dfirstOrderWettingPc(sat);
       const Real dpc2 = dsecondOrderDryingPc(sat);
       const Real dpc = pc * (dpc1 / pc1 + (std::log(pc2) - std::log(pc1)) / (tp2 - tp1) +
                              (sat - tp1) * (dpc2 / pc2 - dpc1 / pc1) / (tp2 - tp1));
       const Real d2pc1 = d2firstOrderWettingPc(sat);
       const Real d2pc2 = d2secondOrderDryingPc(sat);
       d2pc =
           dpc * (dpc1 / pc1 + (std::log(pc2) - std::log(pc1)) / (tp2 - tp1) +
                  (sat - tp1) * (dpc2 / pc2 - dpc1 / pc1) / (tp2 - tp1)) +
           pc *
               (d2pc1 / pc1 - dpc1 * dpc1 / pc1 / pc1 + (dpc2 / pc2 - dpc1 / pc1) / (tp2 - tp1) +
                (dpc2 / pc2 - dpc1 / pc1) / (tp2 - tp1) +
                (sat - tp1) *
                    (d2pc2 / pc2 - dpc2 * dpc2 / pc2 / pc2 - d2pc1 / pc1 + dpc1 * dpc1 / pc1 / pc1) /
                    (tp2 - tp1));
     }
   }
   return d2pc;
 }

◆ d2firstOrderWettingPc()

Real PorousFlowHystereticCapillaryPressure::d2firstOrderWettingPc ( Real sat ) const

private

Returns: d^2(capillary pressure on the first-order wetting curve)/d(sat)^2

Parameters

sat	liquid saturation

Definition at line 467 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by d2capillaryPressureQp().

 {
   const Real max_s = (_w_high_ext_tps[_qp].at(0).strategy ==
                       PorousFlowVanGenuchten::HighCapillaryPressureExtension::POWER)
                          ? 1.0
                          : 1.0 - _s_gr_tps[_qp].at(0);
   const Real sat_to_use = _s_w_tps[_qp].at(0) + (max_s - _s_w_tps[_qp].at(0)) *
                                                     (sat - _hys_sat_tps[_qp].at(0)) /
                                                     (max_s - _hys_sat_tps[_qp].at(0));
   const Real dsat_to_use = (max_s - _s_w_tps[_qp].at(0)) / (max_s - _hys_sat_tps[_qp].at(0));
   return PorousFlowVanGenuchten::d2capillaryPressureHys(sat_to_use,
                                                         _s_l_min,
                                                         _s_gr_tps[_qp].at(0),
                                                         _alpha_w,
                                                         _n_w,
                                                         _w_low_ext_tps[_qp].at(0),
                                                         _w_high_ext_tps[_qp].at(0)) *
          dsat_to_use * dsat_to_use;
 }

◆ d2firstOrderWettingSat()

Real PorousFlowHystereticCapillaryPressure::d2firstOrderWettingSat ( Real pc ) const

private

Returns: d^2(saturation on the first-order wetting curve)/d(pc)^2

Parameters

pc	capillary pressure

Definition at line 572 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by d2liquidSaturationQp().

 {
   const Real d2sat_to_use = PorousFlowVanGenuchten::d2saturationHys(pc,
                                                                     _s_l_min,
                                                                     _s_gr_tps[_qp].at(0),
                                                                     _alpha_w,
                                                                     _n_w,
                                                                     _w_low_ext_tps[_qp].at(0),
                                                                     _w_high_ext_tps[_qp].at(0));
   const Real max_s = (_w_high_ext_tps[_qp].at(0).strategy ==
                       PorousFlowVanGenuchten::HighCapillaryPressureExtension::POWER)
                          ? 1.0
                          : 1.0 - _s_gr_tps[_qp].at(0);
   if (pc <= 0.0 && _w_high_ext_tps[_qp].at(0).strategy ==
                        PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE)
     return 0.0;
   return d2sat_to_use * (max_s - _hys_sat_tps[_qp].at(0)) / (max_s - _s_w_tps[_qp].at(0));
 }

◆ d2liquidSaturationQp()

Real PorousFlowHystereticCapillaryPressure::d2liquidSaturationQp ( Real pc ) const

protected

Returns: d^2(liquid saturation)/d(pc)^2. This uses _hys_order[_qp]

Parameters

pc	capillary pressure

Definition at line 698 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by PorousFlowHystereticInfo::computeQpInfo(), PorousFlow2PhaseHysPP::computeQpProperties(), and PorousFlow1PhaseHysP::computeQpProperties().

 {
   Real d2sat = 0.0;
   if (_hys_order[_qp] == 0) // on primary drying curve
     d2sat = PorousFlowVanGenuchten::d2saturationHys(pc, _s_l_min, 0.0, _alpha_d, _n_d, _low_ext_d);
   else if (_hys_order[_qp] == 1) // first-order wetting
     d2sat = d2firstOrderWettingSat(pc);
   else if (_hys_order[_qp] == 2) // second-order drying
     d2sat = d2secondOrderDryingSat(pc);
   else // third order drying and wetting
   {
     const Real pc_tp1 = _pc_tps[_qp].at(1); // pc on first-order wetting at TP_1
     const Real pc_tp2 = _pc_tps[_qp].at(2); // pc on second-order drying at TP_2
     const Real sat1 = firstOrderWettingSat(pc);
     const Real sat2 = secondOrderDryingSat(pc);
     const Real dsat1 = dfirstOrderWettingSat(pc);
     const Real dsat2 = dsecondOrderDryingSat(pc);
     const Real d2sat1 = d2firstOrderWettingSat(pc);
     const Real d2sat2 = d2secondOrderDryingSat(pc);
     if (pc_tp1 <= 0.0 || pc <= 0.0 || pc_tp2 <= 0.0)
       d2sat = d2sat2;
     else if (pc > pc_tp2)
       d2sat = d2sat2;
     else if (pc < pc_tp1)
       d2sat = d2sat1;
     else
       d2sat = d2sat1 - 1.0 / pc / pc * (sat2 - sat1) / std::log(pc_tp2 / pc_tp1) +
               2.0 / pc * (dsat2 - dsat1) / std::log(pc_tp2 / pc_tp1) +
               (std::log(pc) - std::log(pc_tp1)) * (d2sat2 - d2sat1) / std::log(pc_tp2 / pc_tp1);
   }
   return d2sat;
 }

◆ d2secondOrderDryingPc()

Real PorousFlowHystereticCapillaryPressure::d2secondOrderDryingPc ( Real sat ) const

private

Returns: d^2(capillary pressure on the second-order drying curve)/d(sat)^2

Parameters

sat	liquid saturation

Definition at line 517 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by d2capillaryPressureQp().

 {
   const Real tp0 = _hys_sat_tps[_qp].at(0);
   const Real tp1 = _hys_sat_tps[_qp].at(1);
   const Real s1 = _s_d_tps[_qp].at(1);
   const Real sat_to_use = (sat >= tp0) ? tp0 + (sat - tp0) * (s1 - tp0) / (tp1 - tp0) : sat;
   const Real dsat_to_use = (sat >= tp0) ? (s1 - tp0) / (tp1 - tp0) : 1.0;
   return PorousFlowVanGenuchten::d2capillaryPressureHys(
              sat_to_use, _s_l_min, 0.0, _alpha_d, _n_d, _low_ext_d) *
          dsat_to_use * dsat_to_use;
 }

◆ d2secondOrderDryingSat()

Real PorousFlowHystereticCapillaryPressure::d2secondOrderDryingSat ( Real pc ) const

private

Returns: d^2(saturation on the second-order drying curve)/d(pc)^2

Parameters

pc	capillary pressure

Definition at line 617 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by d2liquidSaturationQp().

 {
   const Real sat_to_use =
       PorousFlowVanGenuchten::saturationHys(pc, _s_l_min, 0.0, _alpha_d, _n_d, _low_ext_d);
   const Real d2sat_to_use =
       PorousFlowVanGenuchten::d2saturationHys(pc, _s_l_min, 0.0, _alpha_d, _n_d, _low_ext_d);
   const Real tp0 = _hys_sat_tps[_qp].at(0);
   const Real tp1 = _hys_sat_tps[_qp].at(1);
   const Real s1 = _s_d_tps[_qp].at(1);
   return (sat_to_use >= tp0) ? d2sat_to_use * (tp1 - tp0) / (s1 - tp0) : d2sat_to_use;
 }

◆ dcapillaryPressureQp()

Real PorousFlowHystereticCapillaryPressure::dcapillaryPressureQp ( Real sat ) const

protected

Returns: d(capillary pressure)/d(sat). This uses _hys_order[_qp]

Parameters

sat	liquid saturation

Definition at line 332 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by PorousFlowHystereticInfo::computeQpInfo(), and PorousFlow2PhaseHysPS::computeQpProperties().

 {
   Real dpc = 0.0;
   if (_hys_order[_qp] == 0) // on primary drying curve
     dpc = PorousFlowVanGenuchten::dcapillaryPressureHys(
         sat, _s_l_min, 0.0, _alpha_d, _n_d, _low_ext_d);
   else if (_hys_order[_qp] == 1) // first-order wetting
     dpc = dfirstOrderWettingPc(sat);
   else if (_hys_order[_qp] == 2) // second-order drying
     dpc = dsecondOrderDryingPc(sat);
   else // third order drying and wetting
   {
     const Real tp1 = _hys_sat_tps[_qp].at(1);
     const Real tp2 = _hys_sat_tps[_qp].at(2);
     const Real pc1 = firstOrderWettingPc(sat);
     const Real pc2 = secondOrderDryingPc(sat);
     // handle cases that occur just at the transition from 3rd to 2nd order, or 3rd to 1st order
     if (sat < tp2)
       dpc = dsecondOrderDryingPc(sat);
     else if (sat > tp1)
       dpc = dfirstOrderWettingPc(sat);
     else if (pc1 >= pc2)
       dpc = dfirstOrderWettingPc(sat);
     else if (pc1 <= 0.0 || pc2 <= 0.0)
       dpc = 0.0;
     else
     {
       const Real pc =
           std::exp(std::log(pc1) + (sat - tp1) * (std::log(pc2) - std::log(pc1)) / (tp2 - tp1));
       const Real dpc1 = dfirstOrderWettingPc(sat);
       const Real dpc2 = dsecondOrderDryingPc(sat);
       dpc = pc * (dpc1 / pc1 + (std::log(pc2) - std::log(pc1)) / (tp2 - tp1) +
                   (sat - tp1) * (dpc2 / pc2 - dpc1 / pc1) / (tp2 - tp1));
     }
   }
   return dpc;
 }

◆ dfirstOrderWettingPc()

Real PorousFlowHystereticCapillaryPressure::dfirstOrderWettingPc ( Real sat ) const

private

Returns: d(capillary pressure on the first-order wetting curve)/d(sat)

Parameters

sat	liquid saturation

Definition at line 446 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by d2capillaryPressureQp(), and dcapillaryPressureQp().

 {
   const Real max_s = (_w_high_ext_tps[_qp].at(0).strategy ==
                       PorousFlowVanGenuchten::HighCapillaryPressureExtension::POWER)
                          ? 1.0
                          : 1.0 - _s_gr_tps[_qp].at(0);
   const Real sat_to_use = _s_w_tps[_qp].at(0) + (max_s - _s_w_tps[_qp].at(0)) *
                                                     (sat - _hys_sat_tps[_qp].at(0)) /
                                                     (max_s - _hys_sat_tps[_qp].at(0));
   const Real dsat_to_use = (max_s - _s_w_tps[_qp].at(0)) / (max_s - _hys_sat_tps[_qp].at(0));
   return PorousFlowVanGenuchten::dcapillaryPressureHys(sat_to_use,
                                                        _s_l_min,
                                                        _s_gr_tps[_qp].at(0),
                                                        _alpha_w,
                                                        _n_w,
                                                        _w_low_ext_tps[_qp].at(0),
                                                        _w_high_ext_tps[_qp].at(0)) *
          dsat_to_use;
 }

◆ dfirstOrderWettingSat()

Real PorousFlowHystereticCapillaryPressure::dfirstOrderWettingSat ( Real pc ) const

private

Returns: d(saturation on the first-order wetting curve)/d(pc)

Parameters

pc	capillary pressure

Definition at line 552 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by d2liquidSaturationQp(), and dliquidSaturationQp().

 {
   const Real dsat_to_use = PorousFlowVanGenuchten::dsaturationHys(pc,
                                                                   _s_l_min,
                                                                   _s_gr_tps[_qp].at(0),
                                                                   _alpha_w,
                                                                   _n_w,
                                                                   _w_low_ext_tps[_qp].at(0),
                                                                   _w_high_ext_tps[_qp].at(0));
   const Real max_s = (_w_high_ext_tps[_qp].at(0).strategy ==
                       PorousFlowVanGenuchten::HighCapillaryPressureExtension::POWER)
                          ? 1.0
                          : 1.0 - _s_gr_tps[_qp].at(0);
   if (pc <= 0.0 && _w_high_ext_tps[_qp].at(0).strategy ==
                        PorousFlowVanGenuchten::HighCapillaryPressureExtension::NONE)
     return 0.0;
   return dsat_to_use * (max_s - _hys_sat_tps[_qp].at(0)) / (max_s - _s_w_tps[_qp].at(0));
 }

◆ dliquidSaturationQp()

Real PorousFlowHystereticCapillaryPressure::dliquidSaturationQp ( Real pc ) const

protected

Returns: d(liquid saturation)/d(pc). This uses _hys_order[_qp]

Parameters

pc	capillary pressure

Definition at line 667 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by PorousFlowHystereticInfo::computeQpInfo(), PorousFlow2PhaseHysPP::computeQpProperties(), and PorousFlow1PhaseHysP::computeQpProperties().

 {
   Real dsat = 0.0;
   if (_hys_order[_qp] == 0) // on primary drying curve
     dsat = PorousFlowVanGenuchten::dsaturationHys(pc, _s_l_min, 0.0, _alpha_d, _n_d, _low_ext_d);
   else if (_hys_order[_qp] == 1) // first-order wetting
     dsat = dfirstOrderWettingSat(pc);
   else if (_hys_order[_qp] == 2) // second-order drying
     dsat = dsecondOrderDryingSat(pc);
   else // third order drying and wetting
   {
     const Real pc_tp1 = _pc_tps[_qp].at(1); // pc on first-order wetting at TP_1
     const Real pc_tp2 = _pc_tps[_qp].at(2); // pc on second-order drying at TP_2
     const Real sat1 = firstOrderWettingSat(pc);
     const Real sat2 = secondOrderDryingSat(pc);
     const Real dsat1 = dfirstOrderWettingSat(pc);
     const Real dsat2 = dsecondOrderDryingSat(pc);
     if (pc_tp1 <= 0.0 || pc <= 0.0 || pc_tp2 <= 0.0)
       dsat = dsat2;
     else if (pc > pc_tp2)
       dsat = dsat2;
     else if (pc < pc_tp1)
       dsat = dsat1;
     else
       dsat = dsat1 + 1.0 / pc * (sat2 - sat1) / std::log(pc_tp2 / pc_tp1) +
              (std::log(pc) - std::log(pc_tp1)) * (dsat2 - dsat1) / std::log(pc_tp2 / pc_tp1);
   }
   return dsat;
 }

◆ dsecondOrderDryingPc()

Real PorousFlowHystereticCapillaryPressure::dsecondOrderDryingPc ( Real sat ) const

private

Returns: d(capillary pressure on the second-order drying curve)/d(sat)

Parameters

sat	liquid saturation

Definition at line 504 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by d2capillaryPressureQp(), and dcapillaryPressureQp().

 {
   const Real tp0 = _hys_sat_tps[_qp].at(0);
   const Real tp1 = _hys_sat_tps[_qp].at(1);
   const Real s1 = _s_d_tps[_qp].at(1);
   const Real sat_to_use = (sat >= tp0) ? tp0 + (sat - tp0) * (s1 - tp0) / (tp1 - tp0) : sat;
   const Real dsat_to_use = (sat >= tp0) ? (s1 - tp0) / (tp1 - tp0) : 1.0;
   return PorousFlowVanGenuchten::dcapillaryPressureHys(
              sat_to_use, _s_l_min, 0.0, _alpha_d, _n_d, _low_ext_d) *
          dsat_to_use;
 }

◆ dsecondOrderDryingSat()

Real PorousFlowHystereticCapillaryPressure::dsecondOrderDryingSat ( Real pc ) const

private

Returns: d(saturation on the second-order drying curve)/d(pc)

Parameters

pc	capillary pressure

Definition at line 604 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by d2liquidSaturationQp(), and dliquidSaturationQp().

 {
   const Real sat_to_use =
       PorousFlowVanGenuchten::saturationHys(pc, _s_l_min, 0.0, _alpha_d, _n_d, _low_ext_d);
   const Real dsat_to_use =
       PorousFlowVanGenuchten::dsaturationHys(pc, _s_l_min, 0.0, _alpha_d, _n_d, _low_ext_d);
   const Real tp0 = _hys_sat_tps[_qp].at(0);
   const Real tp1 = _hys_sat_tps[_qp].at(1);
   const Real s1 = _s_d_tps[_qp].at(1);
   return (sat_to_use >= tp0) ? dsat_to_use * (tp1 - tp0) / (s1 - tp0) : dsat_to_use;
 }

◆ firstOrderWettingPc()

Real PorousFlowHystereticCapillaryPressure::firstOrderWettingPc ( Real sat ) const

private

Returns: capillary pressure on the first-order wetting curve

Parameters

sat	liquid saturation

Definition at line 421 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by capillaryPressureQp(), d2capillaryPressureQp(), dcapillaryPressureQp(), and initQpStatefulProperties().

 {
   // Simplest version is to just use the wetting curve defined by _s_gr_tps[0], but want
   // continuity at sat = _hys_sat_tps[0] (the turning point), so use the following process. The
   // wetting curve is defined for S <= max_s, where
   const Real max_s = (_w_high_ext_tps[_qp].at(0).strategy ==
                       PorousFlowVanGenuchten::HighCapillaryPressureExtension::POWER)
                          ? 1.0
                          : 1.0 - _s_gr_tps[_qp].at(0);
   // define an interpolation: s_to_use smoothly transitions from _s_w_tps[0] (the value of liquid
   // saturation on the wetting curve defined by _s_gr_tps[0]) when sat = _hys_sat_tps[0], to max_s
   // when sat = max_s
   const Real sat_to_use = _s_w_tps[_qp].at(0) + (max_s - _s_w_tps[_qp].at(0)) *
                                                     (sat - _hys_sat_tps[_qp].at(0)) /
                                                     (max_s - _hys_sat_tps[_qp].at(0));
   return PorousFlowVanGenuchten::capillaryPressureHys(sat_to_use,
                                                       _s_l_min,
                                                       _s_gr_tps[_qp].at(0),
                                                       _alpha_w,
                                                       _n_w,
                                                       _w_low_ext_tps[_qp].at(0),
                                                       _w_high_ext_tps[_qp].at(0));
 }

◆ firstOrderWettingSat()

Real PorousFlowHystereticCapillaryPressure::firstOrderWettingSat ( Real pc ) const

private

Returns: saturation on the first-order wetting curve

Parameters

pc	capillary pressure

Definition at line 530 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by d2liquidSaturationQp(), dliquidSaturationQp(), and liquidSaturationQp().

 {
   // this is inverse of firstOrderWettingPc: see that method for comments
   const Real sat_to_use = PorousFlowVanGenuchten::saturationHys(pc,
                                                                 _s_l_min,
                                                                 _s_gr_tps[_qp].at(0),
                                                                 _alpha_w,
                                                                 _n_w,
                                                                 _w_low_ext_tps[_qp].at(0),
                                                                 _w_high_ext_tps[_qp].at(0));
   const Real max_s = (_w_high_ext_tps[_qp].at(0).strategy ==
                       PorousFlowVanGenuchten::HighCapillaryPressureExtension::POWER)
                          ? 1.0
                          : 1.0 - _s_gr_tps[_qp].at(0);
   if (sat_to_use > max_s) // this occurs when using no high extension and pc = 0
     return max_s;
   return (sat_to_use - _s_w_tps[_qp].at(0)) * (max_s - _hys_sat_tps[_qp].at(0)) /
              (max_s - _s_w_tps[_qp].at(0)) +
          _hys_sat_tps[_qp].at(0);
 }

◆ initQpStatefulProperties()

void PorousFlowHystereticCapillaryPressure::initQpStatefulProperties ( )

overrideprotectedvirtual

Reimplemented from PorousFlowVariableBaseTempl< is_ad >.

Reimplemented in PorousFlow1PhaseHysP, PorousFlow2PhaseHysPP, PorousFlow2PhaseHysPS, and PorousFlowHystereticInfo.

Definition at line 193 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by PorousFlowHystereticInfo::initQpStatefulProperties(), PorousFlow2PhaseHysPP::initQpStatefulProperties(), PorousFlow2PhaseHysPS::initQpStatefulProperties(), and PorousFlow1PhaseHysP::initQpStatefulProperties().

 {
   PorousFlowVariableBase::initQpStatefulProperties();
   Real tp_sat;
   Real tp_pc;
   if (_hys_order[_qp] >= 1)
   {
     tp_sat = _hys_sat_tps[_qp].at(0);
     tp_pc = PorousFlowVanGenuchten::capillaryPressureHys(
         tp_sat, _s_l_min, 0.0, _alpha_d, _n_d, _low_ext_d);
     computeTurningPointInfo(0, tp_sat, tp_pc);
   }
   if (_hys_order[_qp] >= 2)
   {
     tp_sat = _hys_sat_tps[_qp].at(1);
     tp_pc = firstOrderWettingPc(tp_sat);
     computeTurningPointInfo(1, tp_sat, tp_pc);
   }
   if (_hys_order[_qp] >= 3)
   {
     tp_sat = _hys_sat_tps[_qp].at(2);
     tp_pc = secondOrderDryingPc(tp_sat);
     computeTurningPointInfo(2, tp_sat, tp_pc);
   }
 }

◆ landSat()

Real PorousFlowHystereticCapillaryPressure::landSat ( Real slDel ) const

protected

Returns: the value of gas saturation (called S_gr^Delta in the markdown documentation) using the Land expression. (This is a function of the liquid saturation at the turning point)

Parameters

slDel the value of the liquid saturation at the turning point

Definition at line 233 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by computeTurningPointInfo().

 {
   const Real a = 1.0 / _s_gr_max - 1.0 / (1.0 - _s_lr);
   return (1.0 - slDel) / (1.0 + a * (1.0 - slDel));
 }

◆ liquidSaturationQp()

Real PorousFlowHystereticCapillaryPressure::liquidSaturationQp ( Real pc ) const

protected

Returns: the value of liquid saturation, given the capillary pressure. This uses _hys_order[_qp]

Parameters

pc	capillary pressure

Definition at line 630 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by PorousFlow1PhaseHysP::buildQpPPSS(), PorousFlow2PhaseHysPP::buildQpPPSS(), and PorousFlowHystereticInfo::computeQpInfo().

 {
   Real sat = 0.0;
   if (_hys_order[_qp] == 0) // on primary drying curve
     sat = PorousFlowVanGenuchten::saturationHys(pc, _s_l_min, 0.0, _alpha_d, _n_d, _low_ext_d);
   else if (_hys_order[_qp] == 1) // first-order wetting
     sat = firstOrderWettingSat(pc);
   else if (_hys_order[_qp] == 2) // second-order drying
     sat = secondOrderDryingSat(pc);
   else // third order drying and wetting
   {
     // NOTE: this is not the exact inverse of the third-order capillary-pressure formula, but only
     // the approximate inverse.  In any simulation, only liquidSaturationQp or capillaryPressureQp
     // are used (not both) so having a slightly different formulation for these two functions is OK
     // Rationale: when pc is close to pc_tp1 then use the first-order wetting curve; when pc is
     // close to pc_tp2 then use the second-order drying curve
     const Real pc_tp1 = _pc_tps[_qp].at(1); // pc on first-order wetting at TP_1
     const Real pc_tp2 = _pc_tps[_qp].at(2); // pc on second-order drying at TP_2
     const Real sat1 = firstOrderWettingSat(pc);
     const Real sat2 = secondOrderDryingSat(pc);
     // handle cases that occur just at the transition from 3rd to 2nd order, or 3rd to 1st order
     if (pc_tp1 <= 0.0 || pc <= 0.0 ||
         pc_tp2 <= 0.0) // only the first condition is strictly necessary as cannot get pc==0 or
                        // pc_tp2=0 without pc_tp1=0 in reality.  The other conditions are added in
                        // case of numerical strangenesses
       sat = sat2;
     else if (pc > pc_tp2)
       sat = sat2;
     else if (pc < pc_tp1)
       sat = sat1;
     else
       sat = sat1 + (std::log(pc) - std::log(pc_tp1)) * (sat2 - sat1) / std::log(pc_tp2 / pc_tp1);
   }
   return sat;
 }

◆ secondOrderDryingPc()

Real PorousFlowHystereticCapillaryPressure::secondOrderDryingPc ( Real sat ) const

private

Returns: capillary pressure on the second-order drying curve

Parameters

sat	liquid saturation

Definition at line 488 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by capillaryPressureQp(), d2capillaryPressureQp(), dcapillaryPressureQp(), and initQpStatefulProperties().

 {
   // Simplest version is to just use the primary drying curve, but want
   // continuity at sat = _hys_sat_tps[0] (the dry-to-wet turning point) and sat = _hys_sat_tps[1]
   // (the wet-to-dry turning point), so use the following process.
   const Real tp0 = _hys_sat_tps[_qp].at(0);
   const Real tp1 = _hys_sat_tps[_qp].at(1);
   const Real s1 = _s_d_tps[_qp].at(1);
   const Real sat_to_use =
       (sat >= tp0) ? tp0 + (sat - tp0) * (s1 - tp0) / (tp1 - tp0)
                    : sat; // final case can occur just at transition from 2nd to 0th order
   return PorousFlowVanGenuchten::capillaryPressureHys(
       sat_to_use, _s_l_min, 0.0, _alpha_d, _n_d, _low_ext_d);
 }

◆ secondOrderDryingSat()

Real PorousFlowHystereticCapillaryPressure::secondOrderDryingSat ( Real pc ) const

private

Returns: saturation on the second-order drying curve

Parameters

pc	capillary pressure

Definition at line 592 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by d2liquidSaturationQp(), dliquidSaturationQp(), and liquidSaturationQp().

 {
   // This is the inverse of secondOrderDryingPc: see that method for comments
   const Real sat_to_use =
       PorousFlowVanGenuchten::saturationHys(pc, _s_l_min, 0.0, _alpha_d, _n_d, _low_ext_d);
   const Real tp0 = _hys_sat_tps[_qp].at(0);
   const Real tp1 = _hys_sat_tps[_qp].at(1);
   const Real s1 = _s_d_tps[_qp].at(1);
   return (sat_to_use >= tp0) ? (sat_to_use - tp0) * (tp1 - tp0) / (s1 - tp0) + tp0 : sat_to_use;
 }

◆ validParams()

InputParameters PorousFlowHystereticCapillaryPressure::validParams ( )

static

Definition at line 16 of file PorousFlowHystereticCapillaryPressure.C.

Referenced by PorousFlowHystereticInfo::validParams(), PorousFlow2PhaseHysPP::validParams(), PorousFlow2PhaseHysPS::validParams(), and PorousFlow1PhaseHysP::validParams().

 {
   InputParameters params = PorousFlowVariableBase::validParams();
   params.addRequiredRangeCheckedParam<Real>(
       "alpha_d",
       "alpha_d > 0",
       "van Genuchten alpha parameter for the primary drying curve.  If using standard units, this "
       "is measured in Pa^-1.  Suggested value is around 1E-5");
   params.addRequiredRangeCheckedParam<Real>(
       "alpha_w",
       "alpha_w > 0",
       "van Genuchten alpha parameter for the primary wetting curve.  If using standard units, this "
       "is measured in Pa^-1.  Suggested value is around 1E-5");
   params.addRequiredRangeCheckedParam<Real>("n_d",
                                             "n_d > 1",
                                             "van Genuchten n parameter for the primary drying "
                                             "curve.  Dimensionless.  Suggested value is around 2");
   params.addRequiredRangeCheckedParam<Real>("n_w",
                                             "n_w > 1",
                                             "van Genuchten n parameter for the primary wetting "
                                             "curve.  Dimensionless.  Suggested value is around 2");
   params.addRequiredRangeCheckedParam<Real>(
       "S_l_min",
       "S_l_min >= 0 & S_l_min < 1",
       "Minimum liquid saturation for which the van Genuchten expression is valid.  If no lower "
       "extension is used then Pc = Pc_max for liquid saturation <= S_l_min");
   params.addRangeCheckedParam<Real>(
       "S_lr",
       0.0,
       "S_lr >= 0 & S_lr < 1",
       "Liquid residual saturation where the liquid relative permeability is zero.  This is used in "
       "the Land expression to find S_gr_del.  Almost definitely you need to set S_lr equal to the "
       "quantity used for your relative-permeability curves.  Almost definitely you should set S_lr "
       "> S_l_min");
   params.addRequiredRangeCheckedParam<Real>(
       "S_gr_max",
       "S_gr_max >= 0",
       "Residual gas saturation.  1 - S_gr_max is the maximum saturation for which the "
       "van Genuchten expression is valid for the wetting curve.  You must ensure S_gr_max < 1 - "
       "S_l_min.  Often S_gr_max = -0.3136 * ln(porosity) - 0.1334 is used");
   params.addRangeCheckedParam<Real>(
       "Pc_max",
       std::numeric_limits<Real>::max(),
       "Pc_max > 0",
       "Value of capillary pressure at which the lower extension commences.  The default value "
       "means capillary pressure uses the van Genuchten expression for S > S_l_min and is "
       "'infinity' for S <= S_l_min.  This will result in poor convergence around S = S_l_min");
   MooseEnum low_ext_enum("none quadratic exponential", "exponential");
   params.addParam<MooseEnum>(
       "low_extension_type",
       low_ext_enum,
       "Type of extension to use for small liquid saturation values.  The extensions modify the "
       "capillary pressure for all saturation values less than S(Pc_max).  That is, if the "
       "van Genuchten "
       "expression would produce Pc > Pc_max, then the extension is used instead.  NONE: Simply "
       "cut-off the capillary-pressure at Pc_max, so that Pc <= Pc_max for all S.  QUADRATIC: Pc is "
       "a quadratic in S that is continuous and differentiable at S(Pc_max) and has zero derivative "
       "at S = 0 (hence, its value at S = 0 will be greater than Pc_max).  EXPONENTIAL: Pc is an "
       "exponential in S that is continuous and differentiable at S(Pc_max) (hence, its value at S "
       "= 0 will be much greater than Pc_max");
   params.addRangeCheckedParam<Real>(
       "high_ratio",
       0.9,
       "high_ratio > 0 & high_ratio < 1",
       "The extension to the wetting curves commences at high_ratio * (1 - S_gr_del), where 1 - "
       "S_gr_del is the value of the liquid saturation when Pc = 0 (on the wetting curve)");
   MooseEnum high_ext_enum("none power", "power");
   params.addParam<MooseEnum>(
       "high_extension_type",
       high_ext_enum,
       "Type of extension to use for the wetting curves when the liquid saturation is around 1.  "
       "The extensions modify the wetting capillary pressure for all saturation values greater than "
       "high_ratio * (1 - S_gr_del), where 1 - S_gr_del is the value of liquid saturation when the "
       "van Genuchten expression gives Pc = 0.  NONE: use the van Genuchten expression and when S > "
       "1 - S_gr_del, set Pc = 0.  POWER: Pc is proportional to (1 - S)^power, where the "
       "coefficient of proportionality and the power are chosen so the resulting curve is "
       "continuous and differentiable");
   params.addClassDescription(
       "This Material computes information that is required for the computation of hysteretic "
       "capillary pressures in single and multi phase situations");
   return params;
 }

Member Data Documentation

◆ _alpha_d

const Real PorousFlowHystereticCapillaryPressure::_alpha_d

protected

van Genuchten alpha parameter for the primary drying curve

Definition at line 86 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by capillaryPressureQp(), computeTurningPointInfo(), d2capillaryPressureQp(), d2liquidSaturationQp(), d2secondOrderDryingPc(), d2secondOrderDryingSat(), dcapillaryPressureQp(), dliquidSaturationQp(), dsecondOrderDryingPc(), dsecondOrderDryingSat(), initQpStatefulProperties(), liquidSaturationQp(), secondOrderDryingPc(), and secondOrderDryingSat().

◆ _alpha_w

const Real PorousFlowHystereticCapillaryPressure::_alpha_w

protected

van Genuchten alpha parameter for the primary wetting curve

Definition at line 88 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by computeTurningPointInfo(), d2firstOrderWettingPc(), d2firstOrderWettingSat(), dfirstOrderWettingPc(), dfirstOrderWettingSat(), firstOrderWettingPc(), and firstOrderWettingSat().

◆ _dgradp_qp_dgradv

template<bool is_ad>

MaterialProperty<std::vector<std::vector<Real> > >* const PorousFlowVariableBaseTempl< is_ad >::_dgradp_qp_dgradv

protectedinherited

d(grad porepressure)/d(grad PorousFlow variable) at the quadpoints

Definition at line 51 of file PorousFlowVariableBase.h.

◆ _dgradp_qp_dv

template<bool is_ad>

MaterialProperty<std::vector<std::vector<RealGradient> > >* const PorousFlowVariableBaseTempl< is_ad >::_dgradp_qp_dv

protectedinherited

d(grad porepressure)/d(PorousFlow variable) at the quadpoints

Definition at line 54 of file PorousFlowVariableBase.h.

◆ _dgrads_qp_dgradv

template<bool is_ad>

MaterialProperty<std::vector<std::vector<Real> > >* const PorousFlowVariableBaseTempl< is_ad >::_dgrads_qp_dgradv

protectedinherited

d(grad saturation)/d(grad PorousFlow variable) at the quadpoints

Definition at line 66 of file PorousFlowVariableBase.h.

◆ _dgrads_qp_dv

template<bool is_ad>

MaterialProperty<std::vector<std::vector<RealGradient> > >* const PorousFlowVariableBaseTempl< is_ad >::_dgrads_qp_dv

protectedinherited

d(grad saturation)/d(PorousFlow variable) at the quadpoints

Definition at line 69 of file PorousFlowVariableBase.h.

◆ _dpc_high

const Real PorousFlowHystereticCapillaryPressure::_dpc_high

protected

d(Pc)/d(S) at the point of high-saturation extension

Definition at line 124 of file PorousFlowHystereticCapillaryPressure.h.

◆ _dpc_low_d

const Real PorousFlowHystereticCapillaryPressure::_dpc_low_d

protected

d(Pc)/dS on the primary drying curve at S = _s_low_d

Definition at line 108 of file PorousFlowHystereticCapillaryPressure.h.

◆ _dpc_low_w

const Real PorousFlowHystereticCapillaryPressure::_dpc_low_w

protected

d(Pc)/dS on the primary wetting curve at S = _s_low_w

Definition at line 114 of file PorousFlowHystereticCapillaryPressure.h.

◆ _dporepressure_dvar

template<bool is_ad>

MaterialProperty<std::vector<std::vector<Real> > >* const PorousFlowVariableBaseTempl< is_ad >::_dporepressure_dvar

protectedinherited

d(porepressure)/d(PorousFlow variable)

Definition at line 45 of file PorousFlowVariableBase.h.

◆ _dsaturation_dvar

template<bool is_ad>

MaterialProperty<std::vector<std::vector<Real> > >* const PorousFlowVariableBaseTempl< is_ad >::_dsaturation_dvar

protectedinherited

d(saturation)/d(PorousFlow variable)

Definition at line 60 of file PorousFlowVariableBase.h.

◆ _gradp_qp

template<bool is_ad>

GenericMaterialProperty<std::vector<RealGradient>, is_ad>* const PorousFlowVariableBaseTempl< is_ad >::_gradp_qp

protectedinherited

Grad(p) at the quadpoints.

Definition at line 48 of file PorousFlowVariableBase.h.

◆ _grads_qp

template<bool is_ad>

GenericMaterialProperty<std::vector<RealGradient>, is_ad>* const PorousFlowVariableBaseTempl< is_ad >::_grads_qp

protectedinherited

Grad(s) at the quadpoints.

Definition at line 63 of file PorousFlowVariableBase.h.

◆ _high_ext

const PorousFlowVanGenuchten::HighCapillaryPressureExtension PorousFlowHystereticCapillaryPressure::_high_ext

protected

Parameters involved in the high-saturation extension of the primary wetting curve.

Definition at line 126 of file PorousFlowHystereticCapillaryPressure.h.

◆ _high_ext_type

const PorousFlowVanGenuchten::HighCapillaryPressureExtension::ExtensionStrategy PorousFlowHystereticCapillaryPressure::_high_ext_type

protected

Type of high-saturation extension of the wetting curves.

Definition at line 118 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by computeTurningPointInfo().

◆ _high_ratio

const Real PorousFlowHystereticCapillaryPressure::_high_ratio

protected

The high-saturation extension to the wetting will commence at _high_ratio * (1 - _s_gr_del)

Definition at line 102 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by computeTurningPointInfo().

◆ _hys_order

const MaterialProperty<unsigned>& PorousFlowHystereticCapillaryPressure::_hys_order

protected

Hysteresis order, as computed by PorousFlowHysteresisOrder.

Definition at line 129 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by capillaryPressureQp(), computeQpProperties(), d2capillaryPressureQp(), d2liquidSaturationQp(), dcapillaryPressureQp(), dliquidSaturationQp(), initQpStatefulProperties(), and liquidSaturationQp().

◆ _hys_order_old

const MaterialProperty<unsigned>& PorousFlowHystereticCapillaryPressure::_hys_order_old

protected

Old value of hysteresis order, as computed by PorousFlowHysteresisOrder.

Definition at line 132 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by computeQpProperties().

◆ _hys_sat_tps

const MaterialProperty<std::array<Real, PorousFlowConstants::MAX_HYSTERESIS_ORDER> >& PorousFlowHystereticCapillaryPressure::_hys_sat_tps

protected

Saturation values at the turning points, as computed by PorousFlowHysteresisOrder.

Definition at line 136 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by capillaryPressureQp(), computeQpProperties(), d2capillaryPressureQp(), d2firstOrderWettingPc(), d2firstOrderWettingSat(), d2secondOrderDryingPc(), d2secondOrderDryingSat(), dcapillaryPressureQp(), dfirstOrderWettingPc(), dfirstOrderWettingSat(), dsecondOrderDryingPc(), dsecondOrderDryingSat(), firstOrderWettingPc(), firstOrderWettingSat(), initQpStatefulProperties(), secondOrderDryingPc(), and secondOrderDryingSat().

◆ _low_ext_d

const PorousFlowVanGenuchten::LowCapillaryPressureExtension PorousFlowHystereticCapillaryPressure::_low_ext_d

protected

Parameters involved in the low-saturation extension of the primary drying curve.

Definition at line 110 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by capillaryPressureQp(), computeTurningPointInfo(), d2capillaryPressureQp(), d2liquidSaturationQp(), d2secondOrderDryingPc(), d2secondOrderDryingSat(), dcapillaryPressureQp(), dliquidSaturationQp(), dsecondOrderDryingPc(), dsecondOrderDryingSat(), initQpStatefulProperties(), liquidSaturationQp(), secondOrderDryingPc(), and secondOrderDryingSat().

◆ _low_ext_type

const PorousFlowVanGenuchten::LowCapillaryPressureExtension::ExtensionStrategy PorousFlowHystereticCapillaryPressure::_low_ext_type

protected

Type of low-saturation extension.

Definition at line 104 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by computeTurningPointInfo().

◆ _low_ext_w

const PorousFlowVanGenuchten::LowCapillaryPressureExtension PorousFlowHystereticCapillaryPressure::_low_ext_w

protected

Parameters involved in the low-saturation extension of the primary wetting curve.

Definition at line 116 of file PorousFlowHystereticCapillaryPressure.h.

◆ _n_d

const Real PorousFlowHystereticCapillaryPressure::_n_d

protected

van Genuchten n parameter for the primary drying curve

Definition at line 90 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by capillaryPressureQp(), computeTurningPointInfo(), d2capillaryPressureQp(), d2liquidSaturationQp(), d2secondOrderDryingPc(), d2secondOrderDryingSat(), dcapillaryPressureQp(), dliquidSaturationQp(), dsecondOrderDryingPc(), dsecondOrderDryingSat(), initQpStatefulProperties(), liquidSaturationQp(), secondOrderDryingPc(), and secondOrderDryingSat().

◆ _n_w

const Real PorousFlowHystereticCapillaryPressure::_n_w

protected

van Genuchten n parameter for the primary wetting curve

Definition at line 92 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by computeTurningPointInfo(), d2firstOrderWettingPc(), d2firstOrderWettingSat(), dfirstOrderWettingPc(), dfirstOrderWettingSat(), firstOrderWettingPc(), and firstOrderWettingSat().

◆ _num_components

template<bool is_ad>

const unsigned int PorousFlowVariableBaseTempl< is_ad >::_num_components

protectedinherited

Number of components.

Definition at line 36 of file PorousFlowVariableBase.h.

Referenced by PorousFlowFluidStateBaseMaterialTempl< is_ad >::PorousFlowFluidStateBaseMaterialTempl().

◆ _num_pf_vars

template<bool is_ad>

const unsigned int PorousFlowVariableBaseTempl< is_ad >::_num_pf_vars

protectedinherited

Number of PorousFlow variables.

Definition at line 39 of file PorousFlowVariableBase.h.

◆ _num_phases

template<bool is_ad>

const unsigned int PorousFlowVariableBaseTempl< is_ad >::_num_phases

protectedinherited

Number of phases.

Definition at line 33 of file PorousFlowVariableBase.h.

Referenced by PorousFlow2PhaseHysPS::computeQpProperties(), PorousFlow1PhaseFullySaturatedTempl< is_ad >::PorousFlow1PhaseFullySaturatedTempl(), PorousFlow1PhaseHysP::PorousFlow1PhaseHysP(), PorousFlow1PhaseMD_Gaussian::PorousFlow1PhaseMD_Gaussian(), PorousFlow1PhasePTempl< is_ad >::PorousFlow1PhasePTempl(), PorousFlow2PhaseHysPP::PorousFlow2PhaseHysPP(), PorousFlow2PhaseHysPS::PorousFlow2PhaseHysPS(), PorousFlow2PhasePPTempl< is_ad >::PorousFlow2PhasePPTempl(), PorousFlowFluidStateBaseMaterialTempl< is_ad >::PorousFlowFluidStateBaseMaterialTempl(), PorousFlowFluidStateSingleComponentTempl< is_ad >::PorousFlowFluidStateSingleComponentTempl(), and PorousFlowFluidStateTempl< is_ad >::PorousFlowFluidStateTempl().

◆ _pc_high

const Real PorousFlowHystereticCapillaryPressure::_pc_high

protected

Pc at the point of high-saturation extension.

Definition at line 122 of file PorousFlowHystereticCapillaryPressure.h.

◆ _pc_max

const Real PorousFlowHystereticCapillaryPressure::_pc_max

protected

Maximum capillary pressure: for Pc above this value, a "lower" extension will be used.

Definition at line 100 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by computeTurningPointInfo().

◆ _pc_older

const MaterialProperty<Real>& PorousFlowHystereticCapillaryPressure::_pc_older

protected

Older value of capillary pressure.

Definition at line 139 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by computeQpProperties().

◆ _pc_tps

MaterialProperty<std::array<Real, PorousFlowConstants::MAX_HYSTERESIS_ORDER> >& PorousFlowHystereticCapillaryPressure::_pc_tps

protected

Nodal or quadpoint values of Pc at the turning points.

Definition at line 142 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by computeTurningPointInfo(), d2liquidSaturationQp(), dliquidSaturationQp(), and liquidSaturationQp().

◆ _porepressure

template<bool is_ad>

GenericMaterialProperty<std::vector<Real>, is_ad>& PorousFlowVariableBaseTempl< is_ad >::_porepressure

protectedinherited

Computed nodal or quadpoint values of porepressure of the phases.

Definition at line 42 of file PorousFlowVariableBase.h.

Referenced by PorousFlow1PhaseMD_Gaussian::buildPS(), PorousFlow1PhaseHysP::buildQpPPSS(), PorousFlow2PhaseHysPS::buildQpPPSS(), PorousFlow2PhaseHysPP::buildQpPPSS(), and PorousFlow1PhaseMD_Gaussian::computeQpProperties().

◆ _s_d_tps

MaterialProperty<std::array<Real, PorousFlowConstants::MAX_HYSTERESIS_ORDER> >& PorousFlowHystereticCapillaryPressure::_s_d_tps

protected

Computed nodal or quadpoint values of saturation on the drying curve at _pc_tps.

Definition at line 145 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by computeTurningPointInfo(), d2secondOrderDryingPc(), d2secondOrderDryingSat(), dsecondOrderDryingPc(), dsecondOrderDryingSat(), secondOrderDryingPc(), and secondOrderDryingSat().

◆ _s_gr_max

const Real PorousFlowHystereticCapillaryPressure::_s_gr_max

protected

Residual gas saturation: 1 - _s_gr_max is the maximum saturation for which the van Genuchten expression is valid for the wetting curve.

Definition at line 98 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by landSat(), and PorousFlowHystereticCapillaryPressure().

◆ _s_gr_tps

MaterialProperty<std::array<Real, PorousFlowConstants::MAX_HYSTERESIS_ORDER> >& PorousFlowHystereticCapillaryPressure::_s_gr_tps

protected

Computed nodal or quadpoint values of S_gr_Del, ie, the Land expression, at the turning points.

Definition at line 148 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by computeTurningPointInfo(), d2firstOrderWettingPc(), d2firstOrderWettingSat(), dfirstOrderWettingPc(), dfirstOrderWettingSat(), firstOrderWettingPc(), and firstOrderWettingSat().

◆ _s_high

const Real PorousFlowHystereticCapillaryPressure::_s_high

protected

Saturation at the point of high-saturation extension.

Definition at line 120 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by PorousFlowHystereticCapillaryPressure().

◆ _s_l_min

const Real PorousFlowHystereticCapillaryPressure::_s_l_min

protected

Minimum liquid saturation for which the van Genuchten expression is valid (Pc(_s_l_min) = infinity)

Definition at line 94 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by capillaryPressureQp(), computeTurningPointInfo(), d2capillaryPressureQp(), d2firstOrderWettingPc(), d2firstOrderWettingSat(), d2liquidSaturationQp(), d2secondOrderDryingPc(), d2secondOrderDryingSat(), dcapillaryPressureQp(), dfirstOrderWettingPc(), dfirstOrderWettingSat(), dliquidSaturationQp(), dsecondOrderDryingPc(), dsecondOrderDryingSat(), firstOrderWettingPc(), firstOrderWettingSat(), initQpStatefulProperties(), liquidSaturationQp(), PorousFlowHystereticCapillaryPressure(), secondOrderDryingPc(), and secondOrderDryingSat().

◆ _s_low_d

const Real PorousFlowHystereticCapillaryPressure::_s_low_d

protected

Saturation on the primary drying curve where low-saturation extension commences.

Definition at line 106 of file PorousFlowHystereticCapillaryPressure.h.

◆ _s_low_w

const Real PorousFlowHystereticCapillaryPressure::_s_low_w

protected

Saturation on the primary wetting curve where low-saturation extension commences.

Definition at line 112 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by PorousFlowHystereticCapillaryPressure().

◆ _s_lr

const Real PorousFlowHystereticCapillaryPressure::_s_lr

protected

Liquid saturation below which the liquid relative permeability is zero.

Definition at line 96 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by landSat(), and PorousFlowHystereticCapillaryPressure().

◆ _s_w_tps

MaterialProperty<std::array<Real, PorousFlowConstants::MAX_HYSTERESIS_ORDER> >& PorousFlowHystereticCapillaryPressure::_s_w_tps

protected

Computed nodal or quadpoint values of liquid saturation on the wetting curve defined by _s_gr_del, at pc = _pc_d_tps.

Definition at line 159 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by computeTurningPointInfo(), d2firstOrderWettingPc(), d2firstOrderWettingSat(), dfirstOrderWettingPc(), dfirstOrderWettingSat(), firstOrderWettingPc(), and firstOrderWettingSat().

◆ _saturation

template<bool is_ad>

GenericMaterialProperty<std::vector<Real>, is_ad>& PorousFlowVariableBaseTempl< is_ad >::_saturation

protectedinherited

Computed nodal or qp saturation of the phases.

Definition at line 57 of file PorousFlowVariableBase.h.

Referenced by PorousFlow1PhaseMD_Gaussian::buildPS(), PorousFlow1PhaseHysP::buildQpPPSS(), PorousFlow2PhaseHysPS::buildQpPPSS(), PorousFlow2PhaseHysPP::buildQpPPSS(), PorousFlowHystereticInfo::computeQpProperties(), PorousFlow1PhaseMD_Gaussian::computeQpProperties(), and PorousFlowHystereticInfo::initQpStatefulProperties().

◆ _w_high_ext_tps

MaterialProperty<std::array<PorousFlowVanGenuchten::HighCapillaryPressureExtension, PorousFlowConstants::MAX_HYSTERESIS_ORDER> >& PorousFlowHystereticCapillaryPressure::_w_high_ext_tps

protected

Nodal or quadpoint values of the high extension of the wetting curve defined by _s_gr_tps.

Definition at line 156 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by computeTurningPointInfo(), d2firstOrderWettingPc(), d2firstOrderWettingSat(), dfirstOrderWettingPc(), dfirstOrderWettingSat(), firstOrderWettingPc(), and firstOrderWettingSat().

◆ _w_low_ext_tps

MaterialProperty<std::array<PorousFlowVanGenuchten::LowCapillaryPressureExtension, PorousFlowConstants::MAX_HYSTERESIS_ORDER> >& PorousFlowHystereticCapillaryPressure::_w_low_ext_tps

protected

Nodal or quadpoint values of the low extension of the wetting curve defined by _s_gr_tps.

Definition at line 152 of file PorousFlowHystereticCapillaryPressure.h.

Referenced by computeTurningPointInfo(), d2firstOrderWettingPc(), d2firstOrderWettingSat(), dfirstOrderWettingPc(), dfirstOrderWettingSat(), firstOrderWettingPc(), and firstOrderWettingSat().

The documentation for this class was generated from the following files:

porous_flow/include/materials/PorousFlowHystereticCapillaryPressure.h
porous_flow/src/materials/PorousFlowHystereticCapillaryPressure.C

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ PorousFlowHystereticCapillaryPressure()

Member Function Documentation

◆ capillaryPressureQp()

◆ computeQpProperties()

◆ computeTurningPointInfo()

◆ d2capillaryPressureQp()

◆ d2firstOrderWettingPc()

◆ d2firstOrderWettingSat()

◆ d2liquidSaturationQp()

◆ d2secondOrderDryingPc()

◆ d2secondOrderDryingSat()

◆ dcapillaryPressureQp()

◆ dfirstOrderWettingPc()

◆ dfirstOrderWettingSat()

◆ dliquidSaturationQp()

◆ dsecondOrderDryingPc()

◆ dsecondOrderDryingSat()

◆ firstOrderWettingPc()

◆ firstOrderWettingSat()

◆ initQpStatefulProperties()

◆ landSat()

◆ liquidSaturationQp()

◆ secondOrderDryingPc()

◆ secondOrderDryingSat()

◆ validParams()

Member Data Documentation

◆ _alpha_d

◆ _alpha_w

◆ _dgradp_qp_dgradv

◆ _dgradp_qp_dv

◆ _dgrads_qp_dgradv

◆ _dgrads_qp_dv

◆ _dpc_high

◆ _dpc_low_d

◆ _dpc_low_w

◆ _dporepressure_dvar

◆ _dsaturation_dvar

◆ _gradp_qp

◆ _grads_qp

◆ _high_ext

◆ _high_ext_type

◆ _high_ratio

◆ _hys_order

◆ _hys_order_old

◆ _hys_sat_tps

◆ _low_ext_d

◆ _low_ext_type

◆ _low_ext_w

◆ _n_d

◆ _n_w

◆ _num_components

◆ _num_pf_vars

◆ _num_phases

◆ _pc_high

◆ _pc_max

◆ _pc_older

◆ _pc_tps

◆ _porepressure

◆ _s_d_tps

◆ _s_gr_max

◆ _s_gr_tps

◆ _s_high

◆ _s_l_min

◆ _s_low_d

◆ _s_low_w

◆ _s_lr

◆ _s_w_tps

◆ _saturation

◆ _w_high_ext_tps

◆ _w_low_ext_tps