standard SUPG relationships valid for the Richards equation. More...
#include <RichardsSUPGstandard.h>
Public Member Functions | |
RichardsSUPGstandard (const InputParameters ¶meters) | |
RealVectorValue | velSUPG (RealTensorValue perm, RealVectorValue gradp, Real density, RealVectorValue gravity) const |
SUPG velocity = -perm*(gradp - density*gravity) This points in direction of information propagation. More... | |
RealTensorValue | dvelSUPG_dgradp (RealTensorValue perm) const |
derivative of SUPG velocity wrt gradient of porepressure More... | |
RealVectorValue | dvelSUPG_dp (RealTensorValue perm, Real density_prime, RealVectorValue gravity) const |
derivative of SUPG velocity wrt porepressure (keeping gradp fixed) More... | |
RealVectorValue | bb (RealVectorValue vel, int dimen, RealVectorValue xi_prime, RealVectorValue eta_prime, RealVectorValue zeta_prime) const |
|bb| ~ 2*velocity/element_length More... | |
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 More... | |
Real | dbb2_dp (RealVectorValue vel, RealVectorValue dvel_dp, RealVectorValue xi_prime, RealVectorValue eta_prime, RealVectorValue zeta_prime) const |
derivative of bb*bb wrt porepressure More... | |
Real | tauSUPG (RealVectorValue vel, Real traceperm, RealVectorValue b) const |
The SUPG tau parameter. More... | |
RealVectorValue | dtauSUPG_dgradp (RealVectorValue vel, RealTensorValue dvel_dgradp, Real traceperm, RealVectorValue b, RealVectorValue db2_dgradp) const |
derivative of tau wrt gradient of porepressure More... | |
Real | dtauSUPG_dp (RealVectorValue vel, RealVectorValue dvel_dp, Real traceperm, RealVectorValue b, Real db2_dp) const |
derivative of tau wrt porepressure (keeping gradp fixed) More... | |
bool | SUPG_trivial () const |
returns false in this case since everything is trivial More... | |
void | initialize () |
void | execute () |
void | finalize () |
Protected Attributes | |
Real | _p_SUPG |
the SUPG pressure parameter This dictates how strong the SUPG is _p_SUPG large means only a little SUPG _p_SUPG small means close to fully-upwind More... | |
Private Member Functions | |
Real | cosh_relation (Real alpha) const |
cosh(alpha)/sinh(alpha) - 1/alpha, modified at extreme values of alpha to prevent overflows More... | |
Real | cosh_relation_prime (Real alpha) const |
derivative of cosh_relation wrt alpha More... | |
standard SUPG relationships valid for the Richards equation.
here i use a formula for "tau" presented in Appendix A of TJR Hughes, M Mallet and A Mizukami `‘A new finite element formulation for computational fluid dynamics:: II. Behond SUPG’' Computer Methods in Applied Mechanics and Engineering 54 (1986) 341–355
Definition at line 28 of file RichardsSUPGstandard.h.
RichardsSUPGstandard::RichardsSUPGstandard | ( | const InputParameters & | parameters | ) |
Definition at line 37 of file RichardsSUPGstandard.C.
|
virtual |
|bb| ~ 2*velocity/element_length
vel | SUPG velocity |
dimen | dimension of problem |
xi_prime | spatial gradient of the isoparametric coordinate xi |
eta_prime | spatial gradient of the isoparametric coordinate eta |
zeta_prime | spatial gradient of the isoparametric coordinate zeta |
Implements RichardsSUPG.
Definition at line 87 of file RichardsSUPGstandard.C.
|
private |
cosh(alpha)/sinh(alpha) - 1/alpha, modified at extreme values of alpha to prevent overflows
Definition at line 66 of file RichardsSUPGstandard.C.
Referenced by dtauSUPG_dgradp(), dtauSUPG_dp(), and tauSUPG().
|
private |
derivative of cosh_relation wrt alpha
Definition at line 76 of file RichardsSUPGstandard.C.
Referenced by dtauSUPG_dgradp(), and dtauSUPG_dp().
|
virtual |
derivative of bb*bb wrt gradient of porepressure
vel | SUPG velocity |
dvel_dgradp | derivative of velocity wrt gradient of porepressure |
xi_prime | spatial gradient of the isoparametric coordinate xi |
eta_prime | spatial gradient of the isoparametric coordinate eta |
zeta_prime | spatial gradient of the isoparametric coordinate zeta |
Implements RichardsSUPG.
Definition at line 104 of file RichardsSUPGstandard.C.
|
virtual |
derivative of bb*bb wrt porepressure
vel | SUPG velocity |
dvel_dp | derivative of velocity wrt porepressure |
xi_prime | spatial gradient of the isoparametric coordinate xi |
eta_prime | spatial gradient of the isoparametric coordinate eta |
zeta_prime | spatial gradient of the isoparametric coordinate zeta |
Implements RichardsSUPG.
Definition at line 118 of file RichardsSUPGstandard.C.
|
virtual |
derivative of tau wrt gradient of porepressure
vel | SUPG velocity |
dvel_dgradp | derivative of the SUPG velocity wrt gradient of porepressure |
traceperm | trace of the permeability tensor for the material |
b | the b parameter: |b| ~ 2*SUPGvelocity/element_length |
db2_dgradp | derivative of b*b wrt gradient of porepressure |
Implements RichardsSUPG.
Definition at line 150 of file RichardsSUPGstandard.C.
|
virtual |
derivative of tau wrt porepressure (keeping gradp fixed)
vel | SUPG velocity |
dvel_dp | derivative of the SUPG velocity wrt porepressure |
traceperm | trace of the permeability tensor for the material |
b | the b parameter: |b| ~ 2*SUPGvelocity/element_length |
db2_dp | derivative of b*b wrt porepressure |
Implements RichardsSUPG.
Definition at line 186 of file RichardsSUPGstandard.C.
|
virtual |
derivative of SUPG velocity wrt gradient of porepressure
perm | permeability tensor |
Implements RichardsSUPG.
Definition at line 52 of file RichardsSUPGstandard.C.
|
virtual |
derivative of SUPG velocity wrt porepressure (keeping gradp fixed)
perm | permeability tensor |
density_prime | derivative of fluid density wrt porepressure |
gravity | gravitational acceleration vector |
Implements RichardsSUPG.
Definition at line 58 of file RichardsSUPGstandard.C.
|
inherited |
Definition at line 31 of file RichardsSUPG.C.
|
inherited |
Definition at line 36 of file RichardsSUPG.C.
|
inherited |
Definition at line 26 of file RichardsSUPG.C.
|
virtual |
returns false in this case since everything is trivial
Implements RichardsSUPG.
Definition at line 220 of file RichardsSUPGstandard.C.
|
virtual |
The SUPG tau parameter.
This dictates how strong the SUPG is
vel | SUPG velocity |
traceperm | trace of the permeability tensor for the material |
b | the b parameter: |b| ~ 2*SUPGvelocity/element_length |
Implements RichardsSUPG.
Definition at line 130 of file RichardsSUPGstandard.C.
|
virtual |
SUPG velocity = -perm*(gradp - density*gravity) This points in direction of information propagation.
perm | permeability tensor |
gradp | gradient of porepressure |
density | fluid density |
gravity | gravitational acceleration vector |
Implements RichardsSUPG.
Definition at line 43 of file RichardsSUPGstandard.C.
|
protected |
the SUPG pressure parameter This dictates how strong the SUPG is _p_SUPG large means only a little SUPG _p_SUPG small means close to fully-upwind
Definition at line 148 of file RichardsSUPGstandard.h.
Referenced by dtauSUPG_dgradp(), dtauSUPG_dp(), and tauSUPG().