www.mooseframework.org
Functions
PorousFlow1PhaseP_BW.C File Reference

Go to the source code of this file.

Functions

template<>
InputParameters validParams< PorousFlow1PhaseP_BW > ()
 

Function Documentation

template<>
InputParameters validParams< PorousFlow1PhaseP_BW > ( )

Definition at line 13 of file PorousFlow1PhaseP_BW.C.

14 {
15  InputParameters params = validParams<PorousFlow1PhaseP>();
16  params.addRequiredRangeCheckedParam<Real>(
17  "Sn",
18  "Sn >= 0",
19  "Low saturation. This must be < Ss, and non-negative. This is BW's "
20  "initial effective saturation, below which effective saturation never goes "
21  "in their simulations/models. If Kn=0 then Sn is the immobile saturation. "
22  "This form of effective saturation is only correct for Kn small.");
23  params.addRangeCheckedParam<Real>(
24  "Ss",
25  1.0,
26  "Ss <= 1",
27  "High saturation. This must be > Sn and <= 1. Effective saturation "
28  "where porepressure = 0. Effective saturation never exceeds this "
29  "value in BW's simulations/models.");
30  params.addRequiredRangeCheckedParam<Real>(
31  "C", "C > 1", "BW's C parameter. Must be > 1. Typical value would be 1.05.");
32  params.addRequiredRangeCheckedParam<Real>(
33  "las",
34  "las > 0",
35  "BW's lambda_s parameter multiplied by (fluid_density * gravity). Must be "
36  "> 0. Typical value would be 1E5");
37  params.addClassDescription("Broadbridge-white form of effective saturation for negligable Kn. "
38  "Then porepressure = -las * ((1 - th) / th - (1 / c) * Ln((C - "
39  "th)/((C - 1) * th))), for th = (Seff - Sn) / (Ss - Sn). A Lambert-W "
40  "function must be evaluated to express Seff in terms of porepressure, "
41  "which can be expensive");
42  return params;
43 }
InputParameters validParams< PorousFlow1PhaseP >()