This postprocessor computes the fluid flux to a RichardsPiecewiseLinearSink. More...

#include <RichardsPiecewiseLinearSinkFlux.h>

Inheritance diagram for RichardsPiecewiseLinearSinkFlux:

Public Member Functions
	RichardsPiecewiseLinearSinkFlux (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpIntegral ()

Protected Attributes
LinearInterpolation	_sink_func
	the sink function, which is a piecewise linear function of porepressure values More...

bool	_use_mobility
	whether to include density*permeability_nn/viscosity in the flux More...

bool	_use_relperm
	whether to include relative permeability in the flux More...

const Function &	_m_func
	the multiplier function More...

const RichardsVarNames &	_richards_name_UO
	holds info regarding the Richards variable names, and their values in the simulation More...

unsigned int	_pvar
	the index into _richards_name_UO corresponding to this Postprocessor's variable eg, if the richards names are 'pwater pgas poil pplasma' and the variable of this Postprocessor is pgas, then _pvar=1 More...

const MaterialProperty< std::vector< Real > > &	_pp
	porepressure values (only the _pvar component is used) More...

const MaterialProperty< std::vector< Real > > &	_viscosity
	fluid viscosity More...

const MaterialProperty< RealTensorValue > &	_permeability
	medium permeability More...

const MaterialProperty< std::vector< Real > > &	_rel_perm
	fluid relative permeability More...

const MaterialProperty< std::vector< Real > > &	_density
	fluid density More...

Detailed Description

This postprocessor computes the fluid flux to a RichardsPiecewiseLinearSink.

The flux is integral_over_boundary of _sink_func*_dt (here _sink_func is a function of porepressure) and if _use_relperm = true, this integrand is multiplied by _rel_perm and if _m_func is entered, this integrand is multiplied by _m_func at the quad point and if _use_mobility = true, this integrand is multiplied by density*knn/viscosity, where knn is n.permeability.n where n is the normal to the boundary

Definition at line 33 of file RichardsPiecewiseLinearSinkFlux.h.

Constructor & Destructor Documentation

◆ RichardsPiecewiseLinearSinkFlux()

RichardsPiecewiseLinearSinkFlux::RichardsPiecewiseLinearSinkFlux ( const InputParameters & parameters )

Definition at line 55 of file RichardsPiecewiseLinearSinkFlux.C.

   : SideIntegralVariablePostprocessor(parameters),
     _sink_func(getParam<std::vector<Real>>("pressures"),
                getParam<std::vector<Real>>("bare_fluxes")),
  
     _use_mobility(getParam<bool>("use_mobility")),
     _use_relperm(getParam<bool>("use_relperm")),
  
     _m_func(getFunction("multiplying_fcn")),
  
     _richards_name_UO(getUserObject<RichardsVarNames>("richardsVarNames_UO")),
     _pvar(_richards_name_UO.richards_var_num(coupled("variable"))),
  
     _pp(getMaterialProperty<std::vector<Real>>("porepressure")),
  
     _viscosity(getMaterialProperty<std::vector<Real>>("viscosity")),
     _permeability(getMaterialProperty<RealTensorValue>("permeability")),
     _rel_perm(getMaterialProperty<std::vector<Real>>("rel_perm")),
     _density(getMaterialProperty<std::vector<Real>>("density"))
 {
 }

Member Function Documentation

◆ computeQpIntegral()

Real RichardsPiecewiseLinearSinkFlux::computeQpIntegral ( )

protectedvirtual

Definition at line 78 of file RichardsPiecewiseLinearSinkFlux.C.

 {
   Real flux = _sink_func.sample(_pp[_qp][_pvar]);
  
   flux *= _m_func.value(_t, _q_point[_qp]);
  
   if (_use_mobility)
   {
     Real k = (_permeability[_qp] * _normals[_qp]) * _normals[_qp];
     flux *= _density[_qp][_pvar] * k / _viscosity[_qp][_pvar];
   }
   if (_use_relperm)
     flux *= _rel_perm[_qp][_pvar];
  
   return flux * _dt;
 }

Member Data Documentation

◆ _density

const MaterialProperty<std::vector<Real> >& RichardsPiecewiseLinearSinkFlux::_density

protected

fluid density

Definition at line 76 of file RichardsPiecewiseLinearSinkFlux.h.

Referenced by computeQpIntegral().

◆ _m_func

const Function& RichardsPiecewiseLinearSinkFlux::_m_func

protected

the multiplier function

Definition at line 51 of file RichardsPiecewiseLinearSinkFlux.h.

Referenced by computeQpIntegral().

◆ _permeability

const MaterialProperty<RealTensorValue>& RichardsPiecewiseLinearSinkFlux::_permeability

protected

medium permeability

Definition at line 70 of file RichardsPiecewiseLinearSinkFlux.h.

Referenced by computeQpIntegral().

◆ _pp

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

protected

porepressure values (only the _pvar component is used)

Definition at line 64 of file RichardsPiecewiseLinearSinkFlux.h.

Referenced by computeQpIntegral().

◆ _pvar

unsigned int RichardsPiecewiseLinearSinkFlux::_pvar

protected

the index into _richards_name_UO corresponding to this Postprocessor's variable eg, if the richards names are 'pwater pgas poil pplasma' and the variable of this Postprocessor is pgas, then _pvar=1

Definition at line 61 of file RichardsPiecewiseLinearSinkFlux.h.

Referenced by computeQpIntegral().