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

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

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

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.

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

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

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.

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

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

Member Data Documentation

_m

Real RichardsRelPermVG::_m

protected

van Genuchten m parameter

Definition at line 52 of file RichardsRelPermVG.h.

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

_simm

Real RichardsRelPermVG::_simm

protected

immobile saturation

Definition at line 49 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