Go to the documentation of this file.
41 velSUPG(RealTensorValue perm, RealVectorValue gradp, Real
density, RealVectorValue gravity)
const;
56 dvelSUPG_dp(RealTensorValue perm, Real density_prime, RealVectorValue gravity)
const;
66 RealVectorValue
bb(RealVectorValue vel,
68 RealVectorValue xi_prime,
69 RealVectorValue eta_prime,
70 RealVectorValue zeta_prime)
const;
81 RealTensorValue dvel_dgradp,
82 RealVectorValue xi_prime,
83 RealVectorValue eta_prime,
84 RealVectorValue zeta_prime)
const;
94 Real
dbb2_dp(RealVectorValue vel,
95 RealVectorValue dvel_dp,
96 RealVectorValue xi_prime,
97 RealVectorValue eta_prime,
98 RealVectorValue zeta_prime)
const;
107 Real
tauSUPG(RealVectorValue vel, Real traceperm, RealVectorValue b)
const;
118 RealTensorValue dvel_dgradp,
121 RealVectorValue db2_dgradp)
const;
132 RealVectorValue dvel_dp,
RealVectorValue velSUPG(RealTensorValue perm, RealVectorValue gradp, Real density, RealVectorValue gravity) const
SUPG velocity = -perm*(gradp - density*gravity) This points in direction of information propagation.
Real cosh_relation_prime(Real alpha) const
derivative of cosh_relation wrt alpha
standard SUPG relationships valid for the Richards equation.
bool SUPG_trivial() const
returns false in this case since everything is trivial
Real dtauSUPG_dp(RealVectorValue vel, RealVectorValue dvel_dp, Real traceperm, RealVectorValue b, Real db2_dp) const
derivative of tau wrt porepressure (keeping gradp fixed)
RealVectorValue dvelSUPG_dp(RealTensorValue perm, Real density_prime, RealVectorValue gravity) const
derivative of SUPG velocity wrt porepressure (keeping gradp fixed)
RichardsSUPGstandard(const InputParameters ¶meters)
RealVectorValue dtauSUPG_dgradp(RealVectorValue vel, RealTensorValue dvel_dgradp, Real traceperm, RealVectorValue b, RealVectorValue db2_dgradp) const
derivative of tau wrt gradient of porepressure
RealTensorValue dvelSUPG_dgradp(RealTensorValue perm) const
derivative of SUPG velocity wrt gradient of porepressure
const std::string density
base class for SUPG of the Richards equation You must override all the functions below with your spec...
Real tauSUPG(RealVectorValue vel, Real traceperm, RealVectorValue b) const
The SUPG tau parameter.
RealVectorValue bb(RealVectorValue vel, int dimen, RealVectorValue xi_prime, RealVectorValue eta_prime, RealVectorValue zeta_prime) const
|bb| ~ 2*velocity/element_length
Real dbb2_dp(RealVectorValue vel, RealVectorValue dvel_dp, RealVectorValue xi_prime, RealVectorValue eta_prime, RealVectorValue zeta_prime) const
derivative of bb*bb wrt porepressure
RealVectorValue dbb2_dgradp(RealVectorValue vel, RealTensorValue dvel_dgradp, RealVectorValue xi_prime, RealVectorValue eta_prime, RealVectorValue zeta_prime) const
derivative of bb*bb wrt gradient of porepressure
Real _p_SUPG
the SUPG pressure parameter This dictates how strong the SUPG is _p_SUPG large means only a little SU...
InputParameters validParams< RichardsSUPGstandard >()
Real cosh_relation(Real alpha) const
cosh(alpha)/sinh(alpha) - 1/alpha, modified at extreme values of alpha to prevent overflows