Q2PRelPermPowerGas Class Reference

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

#include <Q2PRelPermPowerGas.h>

Inheritance diagram for Q2PRelPermPowerGas:

Public Member Functions
	Q2PRelPermPowerGas (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/(1 - simm).

Then relperm = 1 - (n+1)s^n + ns^(n+1) if s<simm, otherwise relperm=1 or 0 This is suitable for gas relative permeability as a function of water saturation in Q2P models

Definition at line 25 of file Q2PRelPermPowerGas.h.

Constructor & Destructor Documentation

Q2PRelPermPowerGas()

Q2PRelPermPowerGas::Q2PRelPermPowerGas

(

const InputParameters &

parameters

)

Definition at line 38 of file Q2PRelPermPowerGas.C.

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

Member Function Documentation

d2relperm()

Real Q2PRelPermPowerGas::d2relperm ( Real  seff ) const

virtual

Second derivative of relative permeability wrt seff.

Parameters

seff	effective saturation

Implements RichardsRelPerm.

Definition at line 83 of file Q2PRelPermPowerGas.C.

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

drelperm()

Real Q2PRelPermPowerGas::drelperm ( Real  seff ) const

virtual

Derivative of relative permeability wrt seff.

Parameters

seff	effective saturation

Implements RichardsRelPerm.

Definition at line 68 of file Q2PRelPermPowerGas.C.

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

virtual

Relative permeability.

Parameters

seff	effective saturation

Implements RichardsRelPerm.

Definition at line 44 of file Q2PRelPermPowerGas.C.

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

protected

exponent

Definition at line 53 of file Q2PRelPermPowerGas.h.

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

_simm

Real Q2PRelPermPowerGas::_simm

protected

immobile saturation

Definition at line 50 of file Q2PRelPermPowerGas.h.

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

---

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

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