Base class for capillary pressure for multiphase flow in porous media. More...

#include <PorousFlowCapillaryPressure.h>

Inheritance diagram for PorousFlowCapillaryPressure:

Public Member Functions
	PorousFlowCapillaryPressure (const InputParameters &parameters)

virtual void	initialize () final

virtual void	initialSetup () override

virtual Real	capillaryPressure (Real saturation, unsigned qp=0) const
	Capillary pressure is calculated as a function of true saturation. More...

virtual Real	dCapillaryPressure (Real saturation, unsigned qp=0) const
	Derivative of capillary pressure wrt true saturation. More...

virtual Real	d2CapillaryPressure (Real saturation, unsigned qp=0) const
	Second derivative of capillary pressure wrt true saturation. More...

virtual DualReal	capillaryPressure (DualReal saturation, unsigned qp=0) const

virtual Real	effectiveSaturation (Real pc, unsigned qp=0) const =0
	Effective saturation as a function of capillary pressure. More...

virtual Real	dEffectiveSaturation (Real pc, unsigned qp=0) const =0
	Derivative of effective saturation wrt capillary pressure. More...

virtual Real	d2EffectiveSaturation (Real pc, unsigned qp=0) const =0
	Second derivative of effective saturation wrt capillary pressure. More...

Real	saturation (Real pc, unsigned qp=0) const
	Saturation as a function of capillary pressure. More...

Real	dSaturation (Real pc, unsigned qp=0) const
	Derivative of saturation wrt capillary pressure. More...

Real	d2Saturation (Real pc, unsigned qp=0) const
	Second derivative of saturation wrt capillary pressure. More...

Protected Member Functions
Real	effectiveSaturationFromSaturation (Real saturation) const
	Effective saturation of liquid phase given liquid saturation and residual liquid saturation. More...

Real	extensionSaturation () const
	Calculates the saturation where the logarithmic extension to capillary pressure meets the raw curve using Newton's method. More...

Real	interceptFunction (Real s) const
	Calculates the saturation where the logarithmic extension to capillary pressure at low saturation. More...

Real	interceptFunctionDeriv (Real s) const
	Calculates the saturation where the logarithmic extension to capillary pressure at low saturation. More...

Real	capillaryPressureLogExt (Real s) const
	The capillary pressure in the logarithmic extension This implementation assumes capillary-pressure is a function of saturation only, and not any other quad-point dependent quantities. More...

Real	dCapillaryPressureLogExt (Real s) const
	The derivative of capillary pressure in the logarithmic extension This implementation assumes capillary-pressure is a function of saturation only, and not any other quad-point dependent quantities. More...

Real	d2CapillaryPressureLogExt (Real s) const
	The second derivative of capillary pressure in the logarithmic extension This implementation assumes capillary-pressure is a function of saturation only, and not any other quad-point dependent quantities. More...

virtual Real	capillaryPressureCurve (Real saturation, unsigned qp=0) const =0
	Raw capillary pressure curve (does not include logarithmic extension) More...

virtual Real	dCapillaryPressureCurve (Real saturation, unsigned qp=0) const =0
	Derivative of raw capillary pressure wrt true saturation. More...

virtual Real	d2CapillaryPressureCurve (Real saturation, unsigned qp=0) const =0
	Second derivative of raw capillary pressure wrt true saturation. More...

Protected Attributes
const Real	_sat_lr
	Liquid residual saturation. More...

const Real	_dseff_ds
	Derivative of effective saturation with respect to saturation. More...

bool	_log_ext
	Flag to use a logarithmic extension for low saturation. More...

const Real	_pc_max
	Maximum capillary pressure (Pa). Note: must be <= 0. More...

Real	_sat_ext
	Saturation where the logarithmic extension meets the raw curve This is computed only for qp=0. More...

Real	_pc_ext
	Capillary pressure where the extension meets the raw curve. More...

Real	_slope_ext
	Gradient of the logarithmic extension This is computed only for qp=0. More...

const Real	_log10
	log(10) More...

Detailed Description

Base class for capillary pressure for multiphase flow in porous media.

To implement an effective saturation formulation, override effectiveSaturation() and derivatives. To implement a capillary pressure curve that will include the logarithmic extension, override capillaryPressureCurve() and derivatives.

Note: Capillary pressure is calculated as a function of true saturation, not effective saturation (saturation minus residual). Derivatives are returned wrt true saturation, so no scaling of them is required in objects using these methods.

Includes an optional logarithmic extension in the low saturation region where capillary pressure can go to infinity as saturation tends to 0. Calculation of logarithmic extension from Webb, A simple extension of two-phase characteristic curves to include the dry region, Water Resources Research 36, 1425-1430 (2000) The logarithmic extension is only computed for qp=0. If your overriding class has no qp dependence (capillary pressure is a function only of saturation, for example) then no changes are required. Otherwise, if you use log_extension, you will have to modify this base class to correctly calculate the log-extension parameters at each quadpoint.

Definition at line 39 of file PorousFlowCapillaryPressure.h.

Constructor & Destructor Documentation

◆ PorousFlowCapillaryPressure()

PorousFlowCapillaryPressure::PorousFlowCapillaryPressure ( const InputParameters & parameters )

Definition at line 37 of file PorousFlowCapillaryPressure.C.

   : DiscreteElementUserObject(parameters),
     _sat_lr(getParam<Real>("sat_lr")),
     _dseff_ds(1.0 / (1.0 - _sat_lr)),
     _log_ext(getParam<bool>("log_extension")),
     _pc_max(getParam<Real>("pc_max")),
     _sat_ext(0.0),
     _pc_ext(0.0),
     _slope_ext(0.0),
     _log10(std::log(10.0))
 {
 }

Member Function Documentation

◆ capillaryPressure() [1/2]

DualReal PorousFlowCapillaryPressure::capillaryPressure	(	DualReal	saturation,
		unsigned	qp = `0`
	)		const

virtual

Definition at line 177 of file PorousFlowCapillaryPressure.C.

 {
   const Real Pc = capillaryPressure(saturation.value(), qp);
   const Real dPc_ds = dCapillaryPressure(saturation.value(), qp);
  
   DualReal result = Pc;
   result.derivatives() = saturation.derivatives() * dPc_ds;
  
   return result;
 }

◆ capillaryPressure() [2/2]

Real PorousFlowCapillaryPressure::capillaryPressure	(	Real	saturation,
		unsigned	qp = `0`
	)		const

virtual

Capillary pressure is calculated as a function of true saturation.

Note that this method includes the ability to use a logarithmic extension at low saturation.

Parameters

saturation	true saturation
qp	quadpoint to use (when capillary-pressure depends on coupled variables, not just saturation)

Returns: capillary pressure (Pa)

Definition at line 64 of file PorousFlowCapillaryPressure.C.

 {
   if (_log_ext && saturation < _sat_ext)
     return capillaryPressureLogExt(saturation);
   else
     return capillaryPressureCurve(saturation, qp);
 }

Referenced by PorousFlow2PhasePS::buildQpPPSS(), capillaryPressure(), initialSetup(), PorousFlowWaterVapor::thermophysicalProperties(), PorousFlowWaterNCG::thermophysicalProperties(), PorousFlowBrineCO2::thermophysicalProperties(), PorousFlowWaterNCG::totalMassFraction(), PorousFlowBrineCO2::totalMassFraction(), PorousFlowWaterNCG::twoPhaseProperties(), and PorousFlowBrineCO2::twoPhaseProperties().

◆ capillaryPressureCurve()

virtual Real PorousFlowCapillaryPressure::capillaryPressureCurve	(	Real	saturation,
		unsigned	qp = `0`
	)		const

protectedpure virtual

Raw capillary pressure curve (does not include logarithmic extension)

Parameters

saturation true saturation

Returns: capillary pressure (Pa)

Implemented in PorousFlowCapillaryPressureBC, PorousFlowCapillaryPressureVG, PorousFlowCapillaryPressureBW, PorousFlowCapillaryPressureConst, and PorousFlowCapillaryPressureRSC.

Referenced by capillaryPressure(), interceptFunction(), and interceptFunctionDeriv().

◆ capillaryPressureLogExt()

Real PorousFlowCapillaryPressure::capillaryPressureLogExt ( Real s ) const

protected

The capillary pressure in the logarithmic extension This implementation assumes capillary-pressure is a function of saturation only, and not any other quad-point dependent quantities.

Parameters

s	liquid saturation

Returns: capillary pressure in logarithmic extension

Definition at line 115 of file PorousFlowCapillaryPressure.C.

 {
   return _pc_ext * std::pow(10.0, _slope_ext * (saturation - _sat_ext));
 }

Referenced by capillaryPressure().

◆ d2CapillaryPressure()

Real PorousFlowCapillaryPressure::d2CapillaryPressure	(	Real	saturation,
		unsigned	qp = `0`
	)		const

virtual

Second derivative of capillary pressure wrt true saturation.

Parameters

saturation	true saturation
qp	quadpoint to use (when capillary-pressure depends on coupled variables, not just saturation)

Returns: second derivative of capillary pressure with respect to true saturation

Definition at line 82 of file PorousFlowCapillaryPressure.C.

 {
   if (_log_ext && saturation < _sat_ext)
     return d2CapillaryPressureLogExt(saturation);
   else
     return d2CapillaryPressureCurve(saturation, qp);
 }

Referenced by PorousFlow2PhasePS::computeQpProperties(), PorousFlowFluidStateSingleComponent::computeQpProperties(), and PorousFlowFluidState::computeQpProperties().

◆ d2CapillaryPressureCurve()

virtual Real PorousFlowCapillaryPressure::d2CapillaryPressureCurve	(	Real	saturation,
		unsigned	qp = `0`
	)		const

protectedpure virtual

Second derivative of raw capillary pressure wrt true saturation.

Parameters

saturation true saturation

Returns: second derivative of capillary pressure with respect to true saturation

Implemented in PorousFlowCapillaryPressureBC, PorousFlowCapillaryPressureVG, PorousFlowCapillaryPressureBW, PorousFlowCapillaryPressureConst, and PorousFlowCapillaryPressureRSC.

Referenced by d2CapillaryPressure(), and interceptFunctionDeriv().

◆ d2CapillaryPressureLogExt()

Real PorousFlowCapillaryPressure::d2CapillaryPressureLogExt ( Real s ) const

protected

The second derivative of capillary pressure in the logarithmic extension This implementation assumes capillary-pressure is a function of saturation only, and not any other quad-point dependent quantities.

Parameters

s	liquid saturation

Returns: second derivative of capillary pressure in logarithmic extension

Definition at line 127 of file PorousFlowCapillaryPressure.C.

 {
   return _pc_ext * _slope_ext * _slope_ext * _log10 * _log10 *
          std::pow(10.0, _slope_ext * (saturation - _sat_ext));
 }

Referenced by d2CapillaryPressure().

◆ d2EffectiveSaturation()

virtual Real PorousFlowCapillaryPressure::d2EffectiveSaturation	(	Real	pc,
		unsigned	qp = `0`
	)		const

pure virtual

Second derivative of effective saturation wrt capillary pressure.

Parameters

pc	capillary pressure
qp	quadpoint to use (when effective saturation depends on coupled variables, not just pc)

Returns: second derivative of effective saturation wrt capillary pressure

Implemented in PorousFlowCapillaryPressureBC, PorousFlowCapillaryPressureVG, PorousFlowCapillaryPressureBW, PorousFlowCapillaryPressureConst, and PorousFlowCapillaryPressureRSC.

Referenced by d2Saturation().

◆ d2Saturation()

Real PorousFlowCapillaryPressure::d2Saturation	(	Real	pc,
		unsigned	qp = `0`
	)		const

Second derivative of saturation wrt capillary pressure.

Parameters

pc	capillary pressure
qp	quadpoint to use (when saturation depends on coupled variables, not just pc)

Returns: second derivative of saturation wrt capillary pressure

Definition at line 109 of file PorousFlowCapillaryPressure.C.

 {
   return d2EffectiveSaturation(pc, qp) * (1.0 - _sat_lr);
 }

Referenced by PorousFlow2PhasePP::computeQpProperties(), and PorousFlow1PhaseP::computeQpProperties().

◆ dCapillaryPressure()

Real PorousFlowCapillaryPressure::dCapillaryPressure	(	Real	saturation,
		unsigned	qp = `0`
	)		const

virtual

Derivative of capillary pressure wrt true saturation.

Parameters

saturation	true saturation
qp	quadpoint to use (when capillary-pressure depends on coupled variables, not just saturation)

Returns: derivative of capillary pressure with respect to true saturation

Definition at line 73 of file PorousFlowCapillaryPressure.C.

 {
   if (_log_ext && saturation < _sat_ext)
     return dCapillaryPressureLogExt(saturation);
   else
     return dCapillaryPressureCurve(saturation, qp);
 }

Referenced by capillaryPressure(), PorousFlow2PhasePS::computeQpProperties(), PorousFlowFluidStateSingleComponent::computeQpProperties(), and PorousFlowFluidState::computeQpProperties().

◆ dCapillaryPressureCurve()

virtual Real PorousFlowCapillaryPressure::dCapillaryPressureCurve	(	Real	saturation,
		unsigned	qp = `0`
	)		const

protectedpure virtual

Derivative of raw capillary pressure wrt true saturation.

Parameters

saturation true saturation

Returns: derivative of capillary pressure with respect to true saturation

Implemented in PorousFlowCapillaryPressureBC, PorousFlowCapillaryPressureVG, PorousFlowCapillaryPressureBW, PorousFlowCapillaryPressureConst, and PorousFlowCapillaryPressureRSC.

Referenced by dCapillaryPressure(), interceptFunction(), and interceptFunctionDeriv().

◆ dCapillaryPressureLogExt()

Real PorousFlowCapillaryPressure::dCapillaryPressureLogExt ( Real s ) const

protected

The derivative of capillary pressure in the logarithmic extension This implementation assumes capillary-pressure is a function of saturation only, and not any other quad-point dependent quantities.

Parameters

s	liquid saturation

Returns: derivative of capillary pressure in logarithmic extension

Definition at line 121 of file PorousFlowCapillaryPressure.C.

 {
   return _pc_ext * _slope_ext * _log10 * std::pow(10.0, _slope_ext * (saturation - _sat_ext));
 }

Referenced by dCapillaryPressure().

◆ dEffectiveSaturation()

virtual Real PorousFlowCapillaryPressure::dEffectiveSaturation	(	Real	pc,
		unsigned	qp = `0`
	)		const

pure virtual

Derivative of effective saturation wrt capillary pressure.

Parameters

pc	capillary pressure (Pa)
qp	quadpoint to use (when effective saturation depends on coupled variables, not just pc)

Returns: derivative of effective saturation wrt capillary pressure

Implemented in PorousFlowCapillaryPressureBC, PorousFlowCapillaryPressureVG, PorousFlowCapillaryPressureBW, PorousFlowCapillaryPressureConst, and PorousFlowCapillaryPressureRSC.

Referenced by dSaturation().

◆ dSaturation()

Real PorousFlowCapillaryPressure::dSaturation	(	Real	pc,
		unsigned	qp = `0`
	)		const

Derivative of saturation wrt capillary pressure.

Parameters

pc	capillary pressure (Pa)
qp	quadpoint to use (when saturation depends on coupled variables, not just pc)

Returns: derivative of saturation wrt capillary pressure

Definition at line 103 of file PorousFlowCapillaryPressure.C.

 {
   return dEffectiveSaturation(pc, qp) * (1.0 - _sat_lr);
 }

Referenced by PorousFlow2PhasePP::computeQpProperties(), and PorousFlow1PhaseP::computeQpProperties().

◆ effectiveSaturation()

virtual Real PorousFlowCapillaryPressure::effectiveSaturation	(	Real	pc,
		unsigned	qp = `0`
	)		const

pure virtual

Effective saturation as a function of capillary pressure.

Parameters

pc	capillary pressure (Pa)
qp	quadpoint to use (when effective saturation depends on coupled variables, not just pc)

Returns: effective saturation

Implemented in PorousFlowCapillaryPressureBC, PorousFlowCapillaryPressureVG, PorousFlowCapillaryPressureBW, PorousFlowCapillaryPressureConst, and PorousFlowCapillaryPressureRSC.

Referenced by saturation().

◆ effectiveSaturationFromSaturation()

Real PorousFlowCapillaryPressure::effectiveSaturationFromSaturation ( Real saturation ) const

protected

Effective saturation of liquid phase given liquid saturation and residual liquid saturation.

Note: not to be mistaken with effectiveSaturation(pc) which is a function of capillary pressure.

Parameters

saturation true saturation

Returns: effective saturation

Definition at line 91 of file PorousFlowCapillaryPressure.C.

 {
   return (saturation - _sat_lr) / (1.0 - _sat_lr);
 }

Referenced by PorousFlowCapillaryPressureVG::capillaryPressureCurve(), PorousFlowCapillaryPressureBC::capillaryPressureCurve(), PorousFlowCapillaryPressureVG::d2CapillaryPressureCurve(), PorousFlowCapillaryPressureBC::d2CapillaryPressureCurve(), PorousFlowCapillaryPressureVG::dCapillaryPressureCurve(), and PorousFlowCapillaryPressureBC::dCapillaryPressureCurve().

◆ extensionSaturation()

Real PorousFlowCapillaryPressure::extensionSaturation ( ) const

protected

Calculates the saturation where the logarithmic extension to capillary pressure meets the raw curve using Newton's method.

This implementation assumes capillary-pressure is a function of saturation only, and not any other quad-point dependent quantities

Returns: saturation where logarithmic extension begins

Definition at line 134 of file PorousFlowCapillaryPressure.C.

 {
   // Initial guess for saturation where extension matches curve
   Real sat = _sat_lr + 0.01;
  
   // Calculate the saturation where the extension matches the derivative of the
   // raw capillary pressure curve
   const unsigned int max_its = 20;
   const Real nr_tol = 1.0e-8;
   unsigned int iter = 0;
  
   while (std::abs(interceptFunction(sat)) > nr_tol)
   {
     sat = sat - interceptFunction(sat) / interceptFunctionDeriv(sat);
  
     iter++;
     if (iter > max_its)
       break;
   }
  
   return sat;
 }

Referenced by initialSetup().

◆ initialize()

virtual void PorousFlowCapillaryPressure::initialize ( )

inlinefinalvirtual

Definition at line 44 of file PorousFlowCapillaryPressure.h.

44 {};

◆ initialSetup()

void PorousFlowCapillaryPressure::initialSetup ( )

overridevirtual

Definition at line 51 of file PorousFlowCapillaryPressure.C.

 {
   // If _log_ext = true, calculate the saturation where the the logarithmic
   // extension meets the raw capillary pressure curve.
   if (_log_ext)
   {
     _sat_ext = extensionSaturation();
     _pc_ext = capillaryPressure(_sat_ext);
     _slope_ext = (std::log10(_pc_ext) - std::log10(_pc_max)) / _sat_ext;
   }
 }

◆ interceptFunction()

Real PorousFlowCapillaryPressure::interceptFunction ( Real s ) const

protected

Calculates the saturation where the logarithmic extension to capillary pressure at low saturation.

This is computed for qp=0 only. This implementation assumes capillary-pressure is a function of saturation only, and not any other quad-point dependent quantities

Parameters

s	effective saturation

Returns: capillary pressure function in the logarithmic extension

Definition at line 158 of file PorousFlowCapillaryPressure.C.

 {
   Real pc = capillaryPressureCurve(saturation);
   Real dpc = dCapillaryPressureCurve(saturation);
  
   return std::log10(pc) - saturation * dpc / (_log10 * pc) - std::log10(_pc_max);
 }

Referenced by extensionSaturation().

◆ interceptFunctionDeriv()

Real PorousFlowCapillaryPressure::interceptFunctionDeriv ( Real s ) const

protected

Calculates the saturation where the logarithmic extension to capillary pressure at low saturation.

This implementation assumes capillary-pressure is a function of saturation only, and not any other quad-point dependent quantities

Parameters

s	effective saturation

Returns: derivative of logarithmic extension function

Definition at line 167 of file PorousFlowCapillaryPressure.C.

 {
   Real pc = capillaryPressureCurve(saturation);
   Real dpc = dCapillaryPressureCurve(saturation);
   Real d2pc = d2CapillaryPressureCurve(saturation);
  
   return saturation * (dpc * dpc / pc - d2pc) / (_log10 * pc);
 }

Referenced by extensionSaturation().

◆ saturation()

Real PorousFlowCapillaryPressure::saturation	(	Real	pc,
		unsigned	qp = `0`
	)		const

Saturation as a function of capillary pressure.

Parameters

pc	capillary pressure (Pa)
qp	quadpoint to use (when saturation depends on coupled variables, not just pc)

Returns: saturation

Definition at line 97 of file PorousFlowCapillaryPressure.C.

 {
   return effectiveSaturation(pc, qp) * (1.0 - _sat_lr) + _sat_lr;
 }

Referenced by PorousFlow1PhaseP::buildQpPPSS(), PorousFlow2PhasePP::buildQpPPSS(), capillaryPressure(), PorousFlowCapillaryPressureVG::capillaryPressureCurve(), PorousFlowCapillaryPressureBC::capillaryPressureCurve(), capillaryPressureLogExt(), d2CapillaryPressure(), PorousFlowCapillaryPressureVG::d2CapillaryPressureCurve(), PorousFlowCapillaryPressureBC::d2CapillaryPressureCurve(), d2CapillaryPressureLogExt(), dCapillaryPressure(), PorousFlowCapillaryPressureVG::dCapillaryPressureCurve(), PorousFlowCapillaryPressureBC::dCapillaryPressureCurve(), dCapillaryPressureLogExt(), effectiveSaturationFromSaturation(), interceptFunction(), and interceptFunctionDeriv().

Member Data Documentation

◆ _dseff_ds

const Real PorousFlowCapillaryPressure::_dseff_ds

protected

Derivative of effective saturation with respect to saturation.

Definition at line 227 of file PorousFlowCapillaryPressure.h.

Referenced by PorousFlowCapillaryPressureVG::d2CapillaryPressureCurve(), PorousFlowCapillaryPressureBC::d2CapillaryPressureCurve(), PorousFlowCapillaryPressureVG::dCapillaryPressureCurve(), and PorousFlowCapillaryPressureBC::dCapillaryPressureCurve().

◆ _log10

const Real PorousFlowCapillaryPressure::_log10

protected

log(10)

Definition at line 249 of file PorousFlowCapillaryPressure.h.

Referenced by d2CapillaryPressureLogExt(), dCapillaryPressureLogExt(), interceptFunction(), and interceptFunctionDeriv().

◆ _log_ext

bool PorousFlowCapillaryPressure::_log_ext

protected

Flag to use a logarithmic extension for low saturation.

Definition at line 229 of file PorousFlowCapillaryPressure.h.

Referenced by capillaryPressure(), d2CapillaryPressure(), dCapillaryPressure(), initialSetup(), PorousFlowCapillaryPressureBW::PorousFlowCapillaryPressureBW(), PorousFlowCapillaryPressureConst::PorousFlowCapillaryPressureConst(), and PorousFlowCapillaryPressureRSC::PorousFlowCapillaryPressureRSC().

◆ _pc_ext

Real PorousFlowCapillaryPressure::_pc_ext

protected

Capillary pressure where the extension meets the raw curve.

This is computed only for qp=0.

Definition at line 242 of file PorousFlowCapillaryPressure.h.

Referenced by capillaryPressureLogExt(), d2CapillaryPressureLogExt(), dCapillaryPressureLogExt(), and initialSetup().

◆ _pc_max

const Real PorousFlowCapillaryPressure::_pc_max

protected

Maximum capillary pressure (Pa). Note: must be <= 0.

Definition at line 231 of file PorousFlowCapillaryPressure.h.

Referenced by PorousFlowCapillaryPressureVG::capillaryPressureCurve(), PorousFlowCapillaryPressureBC::capillaryPressureCurve(), PorousFlowCapillaryPressureVG::d2CapillaryPressureCurve(), PorousFlowCapillaryPressureBC::d2CapillaryPressureCurve(), PorousFlowCapillaryPressureVG::dCapillaryPressureCurve(), PorousFlowCapillaryPressureBC::dCapillaryPressureCurve(), initialSetup(), and interceptFunction().

◆ _sat_ext

Real PorousFlowCapillaryPressure::_sat_ext

protected

Saturation where the logarithmic extension meets the raw curve This is computed only for qp=0.

Definition at line 237 of file PorousFlowCapillaryPressure.h.

Referenced by capillaryPressure(), capillaryPressureLogExt(), d2CapillaryPressure(), d2CapillaryPressureLogExt(), dCapillaryPressure(), dCapillaryPressureLogExt(), and initialSetup().

◆ _sat_lr

const Real PorousFlowCapillaryPressure::_sat_lr

protected

Liquid residual saturation.

Definition at line 225 of file PorousFlowCapillaryPressure.h.

Referenced by d2Saturation(), dSaturation(), effectiveSaturationFromSaturation(), extensionSaturation(), and saturation().

◆ _slope_ext

Real PorousFlowCapillaryPressure::_slope_ext

protected

Gradient of the logarithmic extension This is computed only for qp=0.

Definition at line 247 of file PorousFlowCapillaryPressure.h.

Referenced by capillaryPressureLogExt(), d2CapillaryPressureLogExt(), dCapillaryPressureLogExt(), and initialSetup().

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

porous_flow/include/userobjects/PorousFlowCapillaryPressure.h
porous_flow/src/userobjects/PorousFlowCapillaryPressure.C

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ PorousFlowCapillaryPressure()

Member Function Documentation

◆ capillaryPressure() [1/2]

◆ capillaryPressure() [2/2]

◆ capillaryPressureCurve()

◆ capillaryPressureLogExt()

◆ d2CapillaryPressure()

◆ d2CapillaryPressureCurve()

◆ d2CapillaryPressureLogExt()

◆ d2EffectiveSaturation()

◆ d2Saturation()

◆ dCapillaryPressure()

◆ dCapillaryPressureCurve()

◆ dCapillaryPressureLogExt()

◆ dEffectiveSaturation()

◆ dSaturation()

◆ effectiveSaturation()

◆ effectiveSaturationFromSaturation()

◆ extensionSaturation()

◆ initialize()

◆ initialSetup()

◆ interceptFunction()

◆ interceptFunctionDeriv()

◆ saturation()

Member Data Documentation

◆ _dseff_ds

◆ _log10

◆ _log_ext

◆ _pc_ext

◆ _pc_max

◆ _sat_ext

◆ _sat_lr

◆ _slope_ext