Postprocessor that records the mass flux from porespace to a half-gaussian sink. More...

#include <RichardsHalfGaussianSinkFlux.h>

Inheritance diagram for RichardsHalfGaussianSinkFlux:

Public Member Functions
	RichardsHalfGaussianSinkFlux (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpIntegral ()

Protected Attributes
FEProblemBase &	_feproblem

Real	_maximum
	flux out = maxexp((-0.5(p - centre)/sd)^2) for p<centre, and flux out = max otherwise More...

Real	_sd
	flux out = maxexp((-0.5(p - centre)/sd)^2) for p<centre, and flux out = max otherwise More...

Real	_centre
	flux out = maxexp((-0.5(p - centre)/sd)^2) for p<centre, and flux out = max otherwise More...

const RichardsVarNames &	_richards_name_UO
	holds info regarding the names of the Richards variables and methods for extracting values of these variables More...

unsigned int	_pvar
	the index of this variable in the list of Richards variables held by _richards_name_UO. More...

const Function &	_m_func
	the multiplier function More...

const MaterialProperty< std::vector< Real > > &	_pp
	porepressure (or porepressure vector for multiphase problems) More...

Detailed Description

Postprocessor that records the mass flux from porespace to a half-gaussian sink.

(Positive if fluid is being removed from porespace.) flux out = max*exp((-0.5*(p - centre)/sd)^2) for p<centre, and flux out = max otherwise If a function, _m_func, is used then the flux is multiplied by _m_func. The result is the flux integrated over the specified sideset.

Definition at line 30 of file RichardsHalfGaussianSinkFlux.h.

Constructor & Destructor Documentation

◆ RichardsHalfGaussianSinkFlux()

RichardsHalfGaussianSinkFlux::RichardsHalfGaussianSinkFlux ( const InputParameters & parameters )

Definition at line 45 of file RichardsHalfGaussianSinkFlux.C.

   : SideIntegralVariablePostprocessor(parameters),
     _feproblem(dynamic_cast<FEProblemBase &>(_subproblem)),
     _maximum(getParam<Real>("max")),
     _sd(getParam<Real>("sd")),
     _centre(getParam<Real>("centre")),
     _richards_name_UO(getUserObject<RichardsVarNames>("richardsVarNames_UO")),
     _pvar(_richards_name_UO.richards_var_num(coupled("variable"))),
     _m_func(getFunction("multiplying_fcn")),
     _pp(getMaterialProperty<std::vector<Real>>("porepressure"))
 {
 }

Member Function Documentation

◆ computeQpIntegral()

Real RichardsHalfGaussianSinkFlux::computeQpIntegral ( )

protectedvirtual

Definition at line 59 of file RichardsHalfGaussianSinkFlux.C.

 {
   if (_pp[_qp][_pvar] >= _centre)
     return _maximum * _dt * _m_func.value(_t, _q_point[_qp]);
   else
     return _maximum * exp(-0.5 * std::pow((_pp[_qp][_pvar] - _centre) / _sd, 2)) * _dt *
            _m_func.value(_t, _q_point[_qp]);
 }

Member Data Documentation

◆ _centre

Real RichardsHalfGaussianSinkFlux::_centre

protected

flux out = max*exp((-0.5*(p - centre)/sd)^2) for p<centre, and flux out = max otherwise

Definition at line 47 of file RichardsHalfGaussianSinkFlux.h.

Referenced by computeQpIntegral().

◆ _feproblem

FEProblemBase& RichardsHalfGaussianSinkFlux::_feproblem

protected

Definition at line 38 of file RichardsHalfGaussianSinkFlux.h.

◆ _m_func

const Function& RichardsHalfGaussianSinkFlux::_m_func

protected

the multiplier function

Definition at line 65 of file RichardsHalfGaussianSinkFlux.h.

Referenced by computeQpIntegral().

◆ _maximum

Real RichardsHalfGaussianSinkFlux::_maximum

protected

flux out = max*exp((-0.5*(p - centre)/sd)^2) for p<centre, and flux out = max otherwise

Definition at line 41 of file RichardsHalfGaussianSinkFlux.h.

Referenced by computeQpIntegral().

◆ _pp

const MaterialProperty<std::vector<Real> >& RichardsHalfGaussianSinkFlux::_pp

protected

porepressure (or porepressure vector for multiphase problems)

Definition at line 68 of file RichardsHalfGaussianSinkFlux.h.

Referenced by computeQpIntegral().

◆ _pvar

unsigned int RichardsHalfGaussianSinkFlux::_pvar

protected

the index of this variable in the list of Richards variables held by _richards_name_UO.

Eg if richards_vars = 'pwater pgas poil' in the _richards_name_UO and this kernel has variable = pgas, then _pvar = 1 This is used to index correctly into _viscosity, _seff, etc

Definition at line 62 of file RichardsHalfGaussianSinkFlux.h.

Referenced by computeQpIntegral().