www.mooseframework.org
Functions
PorousFlowFLACrelperm Namespace Reference

FLAC inspired relative permeability relationship. More...

Functions

Real relativePermeability (Real seff, Real m)
 Relative permeability as a function of effective saturation. More...
 
Real dRelativePermeability (Real seff, Real m)
 Derivative of relative permeability with respect to effective saturation. More...
 
Real d2RelativePermeability (Real seff, Real m)
 Second derivative of relative permeability with respect to effective saturation. More...
 

Detailed Description

FLAC inspired relative permeability relationship.

Function Documentation

Real PorousFlowFLACrelperm::d2RelativePermeability ( Real  seff,
Real  m 
)

Second derivative of relative permeability with respect to effective saturation.

Parameters
seffeffective saturation
mvan Genuchten exponent
Returns
second derivative of relative permeability wrt effective saturation

Definition at line 31 of file PorousFlowFLACrelperm.C.

32 {
33  if (seff <= 0.0 || seff >= 1.0)
34  return 0.0;
35  return (1.0 + m) * m * (m - 1.0) * std::pow(seff, m - 2.0) -
36  m * (m + 1.0) * m * std::pow(seff, m - 1.0);
37 }
ExpressionBuilder::EBTerm pow(const ExpressionBuilder::EBTerm &left, T exponent)
Real PorousFlowFLACrelperm::dRelativePermeability ( Real  seff,
Real  m 
)

Derivative of relative permeability with respect to effective saturation.

Parameters
seffeffective saturation
mvan Genuchten exponent
Returns
derivative of relative permeability wrt effective saturation

Definition at line 23 of file PorousFlowFLACrelperm.C.

Referenced by PorousFlowRelativePermeabilityFLAC::dRelativePermeability().

24 {
25  if (seff <= 0.0 || seff >= 1.0)
26  return 0.0;
27  return (1.0 + m) * m * std::pow(seff, m - 1.0) - m * (m + 1.0) * std::pow(seff, m);
28 }
ExpressionBuilder::EBTerm pow(const ExpressionBuilder::EBTerm &left, T exponent)
Real PorousFlowFLACrelperm::relativePermeability ( Real  seff,
Real  m 
)

Relative permeability as a function of effective saturation.

Parameters
seffeffective saturation
mvan Genuchten exponent
Returns
relative permeability

Definition at line 13 of file PorousFlowFLACrelperm.C.

Referenced by PorousFlowRelativePermeabilityFLAC::relativePermeability().

14 {
15  if (seff <= 0.0)
16  return 0.0;
17  else if (seff >= 1.0)
18  return 1.0;
19  return (1.0 + m) * std::pow(seff, m) - m * std::pow(seff, m + 1.0);
20 }
ExpressionBuilder::EBTerm pow(const ExpressionBuilder::EBTerm &left, T exponent)