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

d2RelativePermeability()

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

Second derivative of relative permeability with respect to effective saturation.

Parameters

seff	effective saturation
m	van Genuchten exponent

Returns

second derivative of relative permeability wrt effective saturation

Definition at line 33 of file PorousFlowFLACrelperm.C.

{
  if (seff <= 0.0 || seff >= 1.0)
    return 0.0;
  return (1.0 + m) * m * (m - 1.0) * std::pow(seff, m - 2.0) -
         m * (m + 1.0) * m * std::pow(seff, m - 1.0);
}

dRelativePermeability()

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

Derivative of relative permeability with respect to effective saturation.

Parameters

seff	effective saturation
m	van Genuchten exponent

Returns

derivative of relative permeability wrt effective saturation

Definition at line 25 of file PorousFlowFLACrelperm.C.

{
  if (seff <= 0.0 || seff >= 1.0)
    return 0.0;
  return (1.0 + m) * m * std::pow(seff, m - 1.0) - m * (m + 1.0) * std::pow(seff, m);
}

Referenced by PorousFlowRelativePermeabilityFLAC::dRelativePermeability().

relativePermeability()

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

Relative permeability as a function of effective saturation.

Parameters

seff	effective saturation
m	van Genuchten exponent

Returns

relative permeability

Definition at line 15 of file PorousFlowFLACrelperm.C.

{
  if (seff <= 0.0)
    return 0.0;
  else if (seff >= 1.0)
    return 1.0;
  return (1.0 + m) * std::pow(seff, m) - m * std::pow(seff, m + 1.0);
}

Referenced by PorousFlowRelativePermeabilityFLAC::relativePermeability().