Rogers-Stallybrass-Clements form of capillary pressure. More...

#include <PorousFlowCapillaryPressureRSC.h>

Inheritance diagram for PorousFlowCapillaryPressureRSC:

Public Member Functions
	PorousFlowCapillaryPressureRSC (const InputParameters &parameters)

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

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

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

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

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

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

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 DualReal	capillaryPressure (DualReal saturation, unsigned qp=0) const

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...

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...

Protected Attributes
const Real	_oil_viscosity
	Oil viscosity (which must be twice the water viscocity in this formulation) More...

const Real	_scale_ratio
	Scale ratio: porosity/permeability/beta^2, where beta is chosen by the user. More...

const Real	_shift
	Shift. seff_water = 1/Sqrt(1 + Exp((Pc - shift)/scale)), where scale = 0.25 * scale_ratio * oil_viscosity. More...

const Real	_scale
	Scale = 0.25 * scale_ratio * oil_viscosity. More...

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

Rogers-Stallybrass-Clements form of capillary pressure.

Definition at line 22 of file PorousFlowCapillaryPressureRSC.h.

Constructor & Destructor Documentation

◆ PorousFlowCapillaryPressureRSC()

PorousFlowCapillaryPressureRSC::PorousFlowCapillaryPressureRSC ( const InputParameters & parameters )

Definition at line 36 of file PorousFlowCapillaryPressureRSC.C.

   : PorousFlowCapillaryPressure(parameters),
     _oil_viscosity(getParam<Real>("oil_viscosity")),
     _scale_ratio(getParam<Real>("scale_ratio")),
     _shift(getParam<Real>("shift")),
     _scale(0.25 * _scale_ratio * _oil_viscosity)
 {
   // Set _log_ext to false as no capillary pressure curves are implmented in this class
   _log_ext = false;
 }

Member Function Documentation

◆ capillaryPressure() [1/2]

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

virtualinherited

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

virtualinherited

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(), PorousFlowCapillaryPressure::capillaryPressure(), PorousFlowCapillaryPressure::initialSetup(), PorousFlowWaterVapor::thermophysicalProperties(), PorousFlowWaterNCG::thermophysicalProperties(), PorousFlowBrineCO2::thermophysicalProperties(), PorousFlowWaterNCG::totalMassFraction(), PorousFlowBrineCO2::totalMassFraction(), PorousFlowWaterNCG::twoPhaseProperties(), and PorousFlowBrineCO2::twoPhaseProperties().

◆ capillaryPressureCurve()

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

overridevirtual

Raw capillary pressure curve (does not include logarithmic extension)

Parameters

saturation true saturation

Returns: capillary pressure (Pa)

Implements PorousFlowCapillaryPressure.

Definition at line 48 of file PorousFlowCapillaryPressureRSC.C.

 {
   mooseError("PorousFlowCapillaryPressureRSC::capillaryPressure not implemented");
   return 0.0;
 }

◆ capillaryPressureLogExt()

Real PorousFlowCapillaryPressure::capillaryPressureLogExt ( Real s ) const

protectedinherited

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 PorousFlowCapillaryPressure::capillaryPressure().

◆ d2CapillaryPressure()

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

virtualinherited

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()

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

overridevirtual

Second derivative of raw capillary pressure wrt true saturation.

Parameters

saturation true saturation

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

Implements PorousFlowCapillaryPressure.

Definition at line 62 of file PorousFlowCapillaryPressureRSC.C.

 {
   mooseError("PorousFlowCapillaryPressureRSC::d2CapillaryPressure not implemented");
   return 0.0;
 }

◆ d2CapillaryPressureLogExt()

Real PorousFlowCapillaryPressure::d2CapillaryPressureLogExt ( Real s ) const

protectedinherited

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 PorousFlowCapillaryPressure::d2CapillaryPressure().

◆ d2EffectiveSaturation()

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

overridevirtual

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

Implements PorousFlowCapillaryPressure.

Definition at line 81 of file PorousFlowCapillaryPressureRSC.C.

 {
   return PorousFlowRogersStallybrassClements::d2EffectiveSaturation(-pc, _shift, _scale);
 }

◆ d2Saturation()

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

inherited

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

virtualinherited

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 PorousFlowCapillaryPressure::capillaryPressure(), PorousFlow2PhasePS::computeQpProperties(), PorousFlowFluidStateSingleComponent::computeQpProperties(), and PorousFlowFluidState::computeQpProperties().

◆ dCapillaryPressureCurve()

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

overridevirtual

Derivative of raw capillary pressure wrt true saturation.

Parameters

saturation true saturation

Returns: derivative of capillary pressure with respect to true saturation

Implements PorousFlowCapillaryPressure.

Definition at line 55 of file PorousFlowCapillaryPressureRSC.C.

 {
   mooseError("PorousFlowCapillaryPressureRSC::dCapillaryPressure not implemented");
   return 0.0;
 }

◆ dCapillaryPressureLogExt()

Real PorousFlowCapillaryPressure::dCapillaryPressureLogExt ( Real s ) const

protectedinherited

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 PorousFlowCapillaryPressure::dCapillaryPressure().

◆ dEffectiveSaturation()

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

overridevirtual

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

Implements PorousFlowCapillaryPressure.

Definition at line 75 of file PorousFlowCapillaryPressureRSC.C.

 {
   return -PorousFlowRogersStallybrassClements::dEffectiveSaturation(-pc, _shift, _scale);
 }

◆ dSaturation()

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

inherited

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()

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

overridevirtual

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

Implements PorousFlowCapillaryPressure.

Definition at line 69 of file PorousFlowCapillaryPressureRSC.C.

 {
   return PorousFlowRogersStallybrassClements::effectiveSaturation(-pc, _shift, _scale);
 }

◆ effectiveSaturationFromSaturation()

Real PorousFlowCapillaryPressure::effectiveSaturationFromSaturation ( Real saturation ) const

protectedinherited

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

protectedinherited

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 PorousFlowCapillaryPressure::initialSetup().

◆ initialize()

virtual void PorousFlowCapillaryPressure::initialize ( )

inlinefinalvirtualinherited

Definition at line 44 of file PorousFlowCapillaryPressure.h.

44 {};

◆ initialSetup()

void PorousFlowCapillaryPressure::initialSetup ( )

overridevirtualinherited

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

protectedinherited

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 PorousFlowCapillaryPressure::extensionSaturation().

◆ interceptFunctionDeriv()

Real PorousFlowCapillaryPressure::interceptFunctionDeriv ( Real s ) const

protectedinherited

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 PorousFlowCapillaryPressure::extensionSaturation().

◆ saturation()

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

inherited

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;
 }

Member Data Documentation

◆ _dseff_ds

const Real PorousFlowCapillaryPressure::_dseff_ds

protectedinherited

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

protectedinherited

log(10)

Definition at line 249 of file PorousFlowCapillaryPressure.h.

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

◆ _log_ext

bool PorousFlowCapillaryPressure::_log_ext

protectedinherited

Flag to use a logarithmic extension for low saturation.

Definition at line 229 of file PorousFlowCapillaryPressure.h.

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

◆ _oil_viscosity

const Real PorousFlowCapillaryPressureRSC::_oil_viscosity

protected

Oil viscosity (which must be twice the water viscocity in this formulation)

Definition at line 37 of file PorousFlowCapillaryPressureRSC.h.

◆ _pc_ext

Real PorousFlowCapillaryPressure::_pc_ext

protectedinherited

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 PorousFlowCapillaryPressure::capillaryPressureLogExt(), PorousFlowCapillaryPressure::d2CapillaryPressureLogExt(), PorousFlowCapillaryPressure::dCapillaryPressureLogExt(), and PorousFlowCapillaryPressure::initialSetup().

◆ _pc_max

const Real PorousFlowCapillaryPressure::_pc_max

protectedinherited

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(), PorousFlowCapillaryPressure::initialSetup(), and PorousFlowCapillaryPressure::interceptFunction().

◆ _sat_ext

Real PorousFlowCapillaryPressure::_sat_ext

protectedinherited

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 PorousFlowCapillaryPressure::capillaryPressure(), PorousFlowCapillaryPressure::capillaryPressureLogExt(), PorousFlowCapillaryPressure::d2CapillaryPressure(), PorousFlowCapillaryPressure::d2CapillaryPressureLogExt(), PorousFlowCapillaryPressure::dCapillaryPressure(), PorousFlowCapillaryPressure::dCapillaryPressureLogExt(), and PorousFlowCapillaryPressure::initialSetup().

◆ _sat_lr

const Real PorousFlowCapillaryPressure::_sat_lr

protectedinherited

Liquid residual saturation.

Definition at line 225 of file PorousFlowCapillaryPressure.h.

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

◆ _scale

const Real PorousFlowCapillaryPressureRSC::_scale

protected

Scale = 0.25 * scale_ratio * oil_viscosity.

Definition at line 43 of file PorousFlowCapillaryPressureRSC.h.

Referenced by d2EffectiveSaturation(), dEffectiveSaturation(), and effectiveSaturation().

◆ _scale_ratio

const Real PorousFlowCapillaryPressureRSC::_scale_ratio

protected

Scale ratio: porosity/permeability/beta^2, where beta is chosen by the user.

Definition at line 39 of file PorousFlowCapillaryPressureRSC.h.

◆ _shift

const Real PorousFlowCapillaryPressureRSC::_shift

protected

Shift. seff_water = 1/Sqrt(1 + Exp((Pc - shift)/scale)), where scale = 0.25 * scale_ratio * oil_viscosity.

Definition at line 41 of file PorousFlowCapillaryPressureRSC.h.

Referenced by d2EffectiveSaturation(), dEffectiveSaturation(), and effectiveSaturation().

◆ _slope_ext

Real PorousFlowCapillaryPressure::_slope_ext

protectedinherited

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

Definition at line 247 of file PorousFlowCapillaryPressure.h.

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

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

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

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ PorousFlowCapillaryPressureRSC()

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

◆ _oil_viscosity

◆ _pc_ext

◆ _pc_max

◆ _sat_ext

◆ _sat_lr

◆ _scale

◆ _scale_ratio

◆ _shift

◆ _slope_ext