RichardsRelPermVG Class Reference

Van-Genuchten form of relative permeability as a function of effective saturation. More...

#include <RichardsRelPermVG.h>

Inheritance diagram for RichardsRelPermVG:

Public Member Functions
	RichardsRelPermVG (const InputParameters &parameters)
Real	relperm (Real seff) const
	relative permeability as a function of effective saturation More...
Real	drelperm (Real seff) const
	derivative of relative permeability wrt effective saturation More...
Real	d2relperm (Real seff) const
	second derivative of relative permeability wrt effective saturation More...
void	initialize ()
void	execute ()
void	finalize ()

Protected Attributes
Real	_simm
	immobile saturation More...
Real	_m
	van Genuchten m parameter More...

Detailed Description

Van-Genuchten form of relative permeability as a function of effective saturation.

Definition at line 23 of file RichardsRelPermVG.h.

Constructor & Destructor Documentation

RichardsRelPermVG()

RichardsRelPermVG::RichardsRelPermVG

(

const InputParameters &

parameters

)

Definition at line 38 of file RichardsRelPermVG.C.

  : RichardsRelPerm(parameters), _simm(getParam<Real>("simm")), _m(getParam<Real>("m"))
{
}

Member Function Documentation

d2relperm()

Real RichardsRelPermVG::d2relperm ( Real  seff ) const

virtual

second derivative of relative permeability wrt effective saturation

Parameters

seff	effective saturation

Implements RichardsRelPerm.

Reimplemented in RichardsRelPermVG1.

Definition at line 86 of file RichardsRelPermVG.C.

{
  if (seff >= 1.0)
    return 0.0;
 
  if (seff <= _simm)
    return 0.0;
 
  Real s_internal = (seff - _simm) / (1.0 - _simm);
  Real tmp = 1.0 - std::pow(s_internal, 1.0 / _m);
  Real tmpp = -1.0 / _m * std::pow(s_internal, 1.0 / _m - 1.0);
  Real tmppp = -1.0 / _m * (1.0 / _m - 1.0) * std::pow(s_internal, 1.0 / _m - 2);
  Real tmp2 = 1.0 - std::pow(tmp, _m);
  Real tmp2p = -_m * std::pow(tmp, _m - 1.0) * tmpp;
  Real tmp2pp = -_m * (_m - 1.0) * std::pow(tmp, _m - 2.0) * tmpp * tmpp -
                _m * std::pow(tmp, _m - 1.0) * tmppp;
  // Real krel = std::sqrt(s_internal)*std::pow(tmp2, 2);
  // Real krelp = 0.5 * std::pow(s_internal, -0.5)*std::pow(tmp2, 2) +
  // 2*std::sqrt(s_internal)*tmp2*tmp2p;
  Real krelpp = -0.25 * std::pow(s_internal, -1.5) * tmp2 * tmp2 +
                2.0 * 0.5 * std::pow(s_internal, -0.5) * 2.0 * tmp2 * tmp2p +
                2.0 * std::sqrt(s_internal) * (tmp2p * tmp2p + tmp2 * tmp2pp);
 
  return krelpp / Utility::pow<2>(1.0 - _simm);
}

Referenced by RichardsRelPermVG1::d2relperm(), and RichardsRelPermVG1::RichardsRelPermVG1().

drelperm()

Real RichardsRelPermVG::drelperm ( Real  seff ) const

virtual

derivative of relative permeability wrt effective saturation

Parameters

seff	effective saturation

Implements RichardsRelPerm.

Reimplemented in RichardsRelPermVG1.

Definition at line 66 of file RichardsRelPermVG.C.

{
  if (seff >= 1.0)
    return 0.0;
 
  if (seff <= _simm)
    return 0.0;
 
  Real s_internal = (seff - _simm) / (1.0 - _simm);
  Real tmp = 1.0 - std::pow(s_internal, 1.0 / _m);
  Real tmpp = -1.0 / _m * std::pow(s_internal, 1.0 / _m - 1.0);
  Real tmp2 = 1.0 - std::pow(tmp, _m);
  Real tmp2p = -_m * std::pow(tmp, _m - 1.0) * tmpp;
  // Real krel = std::sqrt(s_internal)*std::pow(tmp2, 2);
  Real krelp =
      0.5 * std::pow(s_internal, -0.5) * tmp2 * tmp2 + 2.0 * std::sqrt(s_internal) * tmp2 * tmp2p;
  return krelp / (1.0 - _simm);
}

Referenced by RichardsRelPermVG1::drelperm(), and RichardsRelPermVG1::RichardsRelPermVG1().

execute()

void RichardsRelPerm::execute ( )

inherited

Definition at line 34 of file RichardsRelPerm.C.

{
}

finalize()

void RichardsRelPerm::finalize ( )

inherited

Definition at line 39 of file RichardsRelPerm.C.

{
}

initialize()

void RichardsRelPerm::initialize ( )

inherited

Definition at line 29 of file RichardsRelPerm.C.

{
}

relperm()

Real RichardsRelPermVG::relperm ( Real  seff ) const

virtual

relative permeability as a function of effective saturation

Parameters

seff	effective saturation

Implements RichardsRelPerm.

Reimplemented in RichardsRelPermVG1.

Definition at line 44 of file RichardsRelPermVG.C.

{
  if (seff >= 1.0)
    return 1.0;
 
  if (seff <= _simm)
    return 0.0;
 
  Real s_internal = (seff - _simm) / (1.0 - _simm);
  Real krel = std::sqrt(s_internal) *
              Utility::pow<2>(1.0 - std::pow(1.0 - std::pow(s_internal, 1.0 / _m), _m));
 
  // bound, just in case
  if (krel < 0.0)
    krel = 0.0;
  if (krel > 1.0)
    krel = 1.0;
 
  return krel;
}

Referenced by RichardsRelPermVG1::relperm(), and RichardsRelPermVG1::RichardsRelPermVG1().

Member Data Documentation

_m

Real RichardsRelPermVG::_m

protected

van Genuchten m parameter

Definition at line 51 of file RichardsRelPermVG.h.

Referenced by d2relperm(), drelperm(), and relperm().

_simm

Real RichardsRelPermVG::_simm

protected

immobile saturation

Definition at line 48 of file RichardsRelPermVG.h.

Referenced by d2relperm(), drelperm(), and relperm().

---

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

• richards/include/userobjects/RichardsRelPermVG.h
• richards/src/userobjects/RichardsRelPermVG.C