27 "Viscosity of oil (gas) phase. It is assumed this is double the water-phase viscosity");
29 "This is porosity/permeability/beta^2, where beta " 30 "may be chosen by the user (RSC define beta<0, but " 31 "MOOSE only uses beta^2, so its sign is " 32 "irrelevant). It has dimensions [time]");
34 params.
addClassDescription(
"Rogers-Stallybrass-Clements version of effective saturation for the " 35 "water phase, valid for residual saturations = 0, and viscosityOil = " 36 "2*viscosityWater. seff_water = 1/Sqrt(1 + Exp(Pc - shift)/scale)), " 37 "where scale = 0.25*scale_ratio*oil_viscosity");
43 _oil_viscosity(getParam<
Real>(
"oil_viscosity")),
44 _scale_ratio(getParam<
Real>(
"scale_ratio")),
45 _shift(getParam<
Real>(
"shift")),
46 _scale(0.25 * _scale_ratio * _oil_viscosity)
53 Real pc = (*p[1])[qp] - (*p[0])[qp];
60 std::vector<Real> & result)
const 62 Real pc = (*p[1])[qp] - (*p[0])[qp];
64 result[0] = -result[1];
70 std::vector<std::vector<Real>> & result)
const 72 Real pc = (*p[1])[qp] - (*p[0])[qp];
74 result[0][1] = -result[1][1];
75 result[1][0] = -result[1][1];
76 result[0][0] = result[1][1];
static Real d2seff(Real pc, Real shift, Real scale)
2nd derivative of effective saturation wrt capillary pressure
Base class for effective saturation as a function of porepressure(s) The functions seff...
static InputParameters validParams()
void dseff(std::vector< const VariableValue *> p, unsigned int qp, std::vector< Real > &result) const
derivative of effective saturation as a function of porepressure
registerMooseObject("RichardsApp", RichardsSeff2waterRSC)
static InputParameters validParams()
RichardsSeff2waterRSC(const InputParameters ¶meters)
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
static Real seff(Real pc, Real shift, Real scale)
effective saturation as a function of capillary pressure
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
Rogers-Stallybrass-Clements version of effective saturation of water phase as a function of (Pwater...
Real seff(std::vector< const VariableValue *> p, unsigned int qp) const
water effective saturation
static Real dseff(Real pc, Real shift, Real scale)
derivative of effective saturation wrt capillary pressure