RichardsRelPermPowerGas Class Reference

PowerGas form of relative permeability Define s = (seff - simm)/(1 - simm). More...

#include <RichardsRelPermPowerGas.h>

Inheritance diagram for RichardsRelPermPowerGas:

Public Member Functions
	RichardsRelPermPowerGas (const InputParameters &parameters)
Real	relperm (Real seff) const
	Relative permeability. More...
Real	drelperm (Real seff) const
	Derivative of relative permeability wrt seff. More...
Real	d2relperm (Real seff) const
	Second derivative of relative permeability wrt seff. More...
void	initialize ()
void	execute ()
void	finalize ()

Protected Attributes
Real	_simm
	immobile saturation More...
Real	_n
	exponent More...

Detailed Description

PowerGas form of relative permeability Define s = (seff - simm)/(1 - simm).

Then relperm = 1 - (n+1)(1-s)^n + n(1-s)^(n+1) if s<simm, otherwise relperm=1

Definition at line 24 of file RichardsRelPermPowerGas.h.

Constructor & Destructor Documentation

RichardsRelPermPowerGas()

RichardsRelPermPowerGas::RichardsRelPermPowerGas

(

const InputParameters &

parameters

)

Definition at line 38 of file RichardsRelPermPowerGas.C.

  : RichardsRelPerm(parameters), _simm(getParam<Real>("simm")), _n(getParam<Real>("n"))
{
}

Member Function Documentation

d2relperm()

Real RichardsRelPermPowerGas::d2relperm ( Real  seff ) const

virtual

Second derivative of relative permeability wrt seff.

Parameters

seff	effective saturation

Implements RichardsRelPerm.

Definition at line 83 of file RichardsRelPermPowerGas.C.

{
  if (seff >= 1.0)
    return 0.0;
 
  if (seff <= _simm)
    return 0.0;
 
  Real s_internal = (seff - _simm) / (1.0 - _simm);
  Real krelpp = -(_n + 1) * _n * (_n - 1) * std::pow(1 - s_internal, _n - 2) +
                _n * (_n + 1) * _n * std::pow(1 - s_internal, _n - 1);
  return krelpp / std::pow(1.0 - _simm, 2);
}

drelperm()

Real RichardsRelPermPowerGas::drelperm ( Real  seff ) const

virtual

Derivative of relative permeability wrt seff.

Parameters

seff	effective saturation

Implements RichardsRelPerm.

Definition at line 68 of file RichardsRelPermPowerGas.C.

{
  if (seff >= 1.0)
    return 0.0;
 
  if (seff <= _simm)
    return 0.0;
 
  Real s_internal = (seff - _simm) / (1.0 - _simm);
  Real krelp = (_n + 1) * _n * std::pow(1 - s_internal, _n - 1) -
               _n * (_n + 1) * std::pow(1 - s_internal, _n);
  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 RichardsRelPermPowerGas::relperm ( Real  seff ) const

virtual

Relative permeability.

Parameters

seff	effective saturation

Implements RichardsRelPerm.

Definition at line 44 of file RichardsRelPermPowerGas.C.

{
  if (seff >= 1.0)
    return 1.0;
 
  if (seff <= _simm)
    return 0.0;
 
  Real s_internal = (seff - _simm) / (1.0 - _simm);
  Real krel = 1 - (_n + 1) * std::pow(1 - s_internal, _n) + _n * std::pow(1 - s_internal, _n + 1);
 
  // bound, just in case
  if (krel < 0)
  {
    krel = 0;
  }
  if (krel > 1)
  {
    krel = 1;
  }
  return krel;
}

Member Data Documentation

_n

Real RichardsRelPermPowerGas::_n

protected

exponent

Definition at line 52 of file RichardsRelPermPowerGas.h.

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

_simm

Real RichardsRelPermPowerGas::_simm

protected

immobile saturation

Definition at line 49 of file RichardsRelPermPowerGas.h.

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

---

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

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