Shifted van-Genuchten water effective saturation as a function of (Pwater, Pgas), and its derivs wrt to those pressures. More...

#include <RichardsSeff2gasVGshifted.h>

Inheritance diagram for RichardsSeff2gasVGshifted:

[legend]

Public Member Functions
	RichardsSeff2gasVGshifted (const InputParameters &parameters)

Real	seff (std::vector< const VariableValue * > p, unsigned int qp) const
	gas effective saturation More...

void	dseff (std::vector< const VariableValue * > p, unsigned int qp, std::vector< Real > &result) const
	derivative of effective saturation as a function of porepressure More...

void	d2seff (std::vector< const VariableValue * > p, unsigned int qp, std::vector< std::vector< Real >> &result) const
	second derivative of effective saturation as a function of porepressure More...

void	initialize ()

void	execute ()

void	finalize ()

Protected Attributes
Real	_al
	van Genuchten alpha parameter More...

Real	_m
	van Genuchten m parameter More...

Real	_shift
	shift More...

Real	_scale
	scale More...

Detailed Description

Shifted van-Genuchten water effective saturation as a function of (Pwater, Pgas), and its derivs wrt to those pressures.

Note that the water pressure appears first in the tuple (Pwater, Pgas) This takes the original van-Genuchten Seff = Seff(Pwater-Pgas), and shifts it to the right by "shift", and scales the result so 0<=Seff<=1. The purpose of this is so dSeff/dP>0 at P=0.

Definition at line 28 of file RichardsSeff2gasVGshifted.h.

Constructor & Destructor Documentation

◆ RichardsSeff2gasVGshifted()

RichardsSeff2gasVGshifted::RichardsSeff2gasVGshifted ( const InputParameters & parameters )

Definition at line 47 of file RichardsSeff2gasVGshifted.C.

   : RichardsSeff(parameters),
     _al(getParam<Real>("al")),
     _m(getParam<Real>("m")),
     _shift(getParam<Real>("shift"))
 {
   _scale = RichardsSeffVG::seff(-_shift, _al, _m);
 }

Member Function Documentation

◆ d2seff()

void RichardsSeff2gasVGshifted::d2seff	(	std::vector< const VariableValue * >	p,
		unsigned int	qp,
		std::vector< std::vector< Real >> &	result
	)		const

virtual

second derivative of effective saturation as a function of porepressure

Parameters

p	porepressure in the element. Note that (*p[0])[qp] is the porepressure at quadpoint qp
qp	the quad point to evaluate effective saturation at
result	the derivtives will be placed in this array

Implements RichardsSeff.

Definition at line 76 of file RichardsSeff2gasVGshifted.C.

 {
   Real negpc = (*p[0])[qp] - (*p[1])[qp];
   negpc = negpc - _shift;
   result[0][0] = -RichardsSeffVG::d2seff(negpc, _al, _m) / _scale;
   result[0][1] = -result[0][0];
   result[1][0] = -result[0][0];
   result[1][1] = result[0][0];
 }

◆ dseff()

void RichardsSeff2gasVGshifted::dseff	(	std::vector< const VariableValue * >	p,
		unsigned int	qp,
		std::vector< Real > &	result
	)		const

virtual

derivative of effective saturation as a function of porepressure

Parameters

p	porepressure in the element. Note that (*p[0])[qp] is the porepressure at quadpoint qp
qp	the quad point to evaluate effective saturation at
result	the derivtives will be placed in this array

Implements RichardsSeff.

Definition at line 65 of file RichardsSeff2gasVGshifted.C.

 {
   Real negpc = (*p[0])[qp] - (*p[1])[qp];
   negpc = negpc - _shift;
   result[0] = -RichardsSeffVG::dseff(negpc, _al, _m) / _scale;
   result[1] = -result[0];
 }

◆ execute()

void RichardsSeff::execute ( )

inherited

Definition at line 32 of file RichardsSeff.C.

33 {

34 }

◆ finalize()

void RichardsSeff::finalize ( )

inherited

Definition at line 37 of file RichardsSeff.C.

38 {

39 }

◆ initialize()

void RichardsSeff::initialize ( )

inherited

Definition at line 27 of file RichardsSeff.C.

28 {

29 }

◆ seff()

Real RichardsSeff2gasVGshifted::seff	(	std::vector< const VariableValue * >	p,
		unsigned int	qp
	)		const

virtual

gas effective saturation

Parameters

p	porepressures. Here (p[0])[qp] is the water pressure at quadpoint qp, and (p[1])[qp] is the gas porepressure
qp	the quadpoint to evaluate effective saturation at

Implements RichardsSeff.

Definition at line 57 of file RichardsSeff2gasVGshifted.C.

 {
   Real negpc = (*p[0])[qp] - (*p[1])[qp];
   negpc = negpc - _shift;
   return std::max(1 - RichardsSeffVG::seff(negpc, _al, _m) / _scale, 0.0);
 }