"Broadbridge-White" form of relative permeability as a function of effective saturation P Broadbridge and I White `‘Constant rate rainfall infiltration: A versatile nonlinear model 1. More...

#include <RichardsRelPermBW.h>

Inheritance diagram for RichardsRelPermBW:

[legend]

Public Member Functions
	RichardsRelPermBW (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	_sn

Real	_ss

Real	_kn

Real	_ks

Real	_c

Real	_coef

Detailed Description

"Broadbridge-White" form of relative permeability as a function of effective saturation P Broadbridge and I White `‘Constant rate rainfall infiltration: A versatile nonlinear model 1.

Analytic Solution’', Water Resources Research 24 (1988) 145-154)

Definition at line 24 of file RichardsRelPermBW.h.

Constructor & Destructor Documentation

◆ RichardsRelPermBW()

RichardsRelPermBW::RichardsRelPermBW ( const InputParameters & parameters )

Definition at line 53 of file RichardsRelPermBW.C.

   : RichardsRelPerm(parameters),
     _sn(getParam<Real>("Sn")),
     _ss(getParam<Real>("Ss")),
     _kn(getParam<Real>("Kn")),
     _ks(getParam<Real>("Ks")),
     _c(getParam<Real>("C"))
 {
   if (_ss <= _sn)
     mooseError("In BW relative permeability Sn set to ",
                _sn,
                " and Ss set to ",
                _ss,
                " but these must obey Ss > Sn");
   if (_ks <= _kn)
     mooseError("In BW relative permeability Kn set to ",
                _kn,
                " and Ks set to ",
                _ks,
                " but these must obey Ks > Kn");
   _coef = (_ks - _kn) * (_c - 1); // shorthand coefficient
 }

Member Function Documentation

◆ d2relperm()

Real RichardsRelPermBW::d2relperm ( Real seff ) const

virtual

second derivative of relative permeability wrt effective saturation

Parameters

seff	effective saturation

Implements RichardsRelPerm.

Definition at line 106 of file RichardsRelPermBW.C.

 {
   if (seff <= _sn)
     return 0.0;
  
   if (seff >= _ss)
     return 0.0;
  
   const Real s_internal = (seff - _sn) / (_ss - _sn);
   const Real krelpp =
       _coef * (2.0 / (_c - s_internal) + 4.0 * s_internal / Utility::pow<2>(_c - s_internal) +
                2.0 * Utility::pow<2>(s_internal) / Utility::pow<3>(_c - s_internal));
   return krelpp / Utility::pow<2>(_ss - _sn);
 }

◆ drelperm()

Real RichardsRelPermBW::drelperm ( Real seff ) const

virtual

derivative of relative permeability wrt effective saturation

Parameters

seff	effective saturation

Implements RichardsRelPerm.

Definition at line 91 of file RichardsRelPermBW.C.

 {
   if (seff <= _sn)
     return 0.0;
  
   if (seff >= _ss)
     return 0.0;
  
   const Real s_internal = (seff - _sn) / (_ss - _sn);
   const Real krelp = _coef * (2.0 * s_internal / (_c - s_internal) +
                               Utility::pow<2>(s_internal) / Utility::pow<2>(_c - s_internal));
   return krelp / (_ss - _sn);
 }

◆ execute()

void RichardsRelPerm::execute ( )

inherited

Definition at line 34 of file RichardsRelPerm.C.

35 {

36 }

◆ finalize()

void RichardsRelPerm::finalize ( )

inherited

Definition at line 39 of file RichardsRelPerm.C.

40 {

41 }

◆ initialize()

void RichardsRelPerm::initialize ( )

inherited

Definition at line 29 of file RichardsRelPerm.C.

30 {

31 }

◆ relperm()

Real RichardsRelPermBW::relperm ( Real seff ) const

virtual

relative permeability as a function of effective saturation

Parameters

seff	effective saturation

Implements RichardsRelPerm.

Definition at line 77 of file RichardsRelPermBW.C.

 {
   if (seff <= _sn)
     return _kn;
  
   if (seff >= _ss)
     return _ks;
  
   const Real s_internal = (seff - _sn) / (_ss - _sn);
   const Real krel = _kn + _coef * Utility::pow<2>(s_internal) / (_c - s_internal);
   return krel;
 }