www.mooseframework.org
Classes | Functions
RichardsRelPermVG1.h File Reference

Go to the source code of this file.

Classes

class  RichardsRelPermVG1
 Van-Genuchten form of relative permeability when seff <= _scut cubic relative permeability for seff >= _scut These two match in value and derivative at seff = _scut and relperm = 1 for seff = 1. More...
 

Functions

template<>
InputParameters validParams< RichardsRelPermVG1 > ()
 

Function Documentation

template<>
InputParameters validParams< RichardsRelPermVG1 > ( )

Definition at line 14 of file RichardsRelPermVG1.C.

15 {
16  InputParameters params = validParams<RichardsRelPermVG>();
17  params.addRequiredRangeCheckedParam<Real>(
18  "simm",
19  "simm >=0 & simm < 1",
20  "Immobile saturation. Must be between 0 and 1. Define s = "
21  "(seff - simm)/(1 - simm). Then relperm = s^(1/2) * (1 - (1 "
22  "- s^(1/m))^m)^2");
23  params.addRequiredRangeCheckedParam<Real>(
24  "m",
25  "m > 0 & m < 1",
26  "van-Genuchten m parameter. Must be between 0 and 1, and optimally "
27  "should be set >0.5. Define s = (seff - simm)/(1 - simm). Then "
28  "relperm = s^(1/2) * (1 - (1 - s^(1/m))^m)^2");
29  params.addRequiredRangeCheckedParam<Real>(
30  "scut", "scut > 0 & scut < 1", "cutoff in effective saturation.");
31  params.addClassDescription("VG1 form of relative permeability. Define s = (seff - simm)/(1 - "
32  "simm). Then relperm = s^(1/2) * (1 - (1 - s^(1/m))^m)^2, if s>0, "
33  "and relperm=0 otherwise");
34  return params;
35 }
InputParameters validParams< RichardsRelPermVG >()