23 "Maximum of the flux (measured in kg.m^-2.s^-1). Flux out " 24 "= max*exp((-0.5*(p - centre)/sd)^2) for p<centre, and Flux " 25 "out = max for p>centre. Note, to make this a source " 26 "rather than a sink, let max<0");
28 "Standard deviation of the Gaussian (measured in Pa). Flux " 29 "out = max*exp((-0.5*(p - centre)/sd)^2) for p<centre, and " 30 "Flux out = max for p>centre.");
32 "Centre of the Gaussian (measured in Pa). Flux out = " 33 "max*exp((-0.5*(p - centre)/sd)^2) for p<centre, and " 34 "Flux out = max for p>centre.");
38 "The flux will be multiplied by this spatially-and-temporally varying function.");
40 "richardsVarNames_UO",
"The UserObject that holds the list of Richards variable names.");
47 _maximum(getParam<
Real>(
"max")),
48 _sd(getParam<
Real>(
"sd")),
49 _centre(getParam<
Real>(
"centre")),
51 _pvar(_richards_name_UO.richards_var_num(coupled(
"variable"))),
52 _m_func(getFunction(
"multiplying_fcn")),
53 _pp(getMaterialProperty<
std::vector<
Real>>(
"porepressure"))
static InputParameters validParams()
const MooseArray< Point > & _q_point
This holds maps between pressure_var or pressure_var, sat_var used in RichardsMaterial and kernels...
Real _maximum
flux out = max*exp((-0.5*(p - centre)/sd)^2) for p<centre, and flux out = max otherwise ...
Postprocessor that records the mass flux from porespace to a half-gaussian sink.
RichardsHalfGaussianSinkFlux(const InputParameters ¶meters)
registerMooseObject("RichardsApp", RichardsHalfGaussianSinkFlux)
const Function & _m_func
the multiplier function
virtual Real computeQpIntegral()
static InputParameters validParams()
const MaterialProperty< std::vector< Real > > & _pp
porepressure (or porepressure vector for multiphase problems)
Real _centre
flux out = max*exp((-0.5*(p - centre)/sd)^2) for p<centre, and flux out = max otherwise ...
unsigned int _pvar
the index of this variable in the list of Richards variables held by _richards_name_UO.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
virtual Real value(Real t, const Point &p) const
MooseUnits pow(const MooseUnits &, int)
Real _sd
flux out = max*exp((-0.5*(p - centre)/sd)^2) for p<centre, and flux out = max otherwise ...