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

#include <Q2PPiecewiseLinearSinkFlux.h>

Inheritance diagram for Q2PPiecewiseLinearSinkFlux:

Public Member Functions
	Q2PPiecewiseLinearSinkFlux (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...

const Function &	_m_func
	the multiplier function More...

const VariableValue &	_pp
	the porepressure variable 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 RichardsDensity *	_density
	fluid density, optional More...

Real	_viscosity
	fluid viscosity, optional More...

const RichardsRelPerm *	_relperm
	fluid relative permeaility, optional More...

const VariableValue &	_sat
	saturation variable, optional More...

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

Detailed Description

This postprocessor computes the fluid flux to a Q2PPiecewiseLinearSink.

The flux is integral_over_boundary of _sink_func*_dt (here _sink_func is a function of porepressure) and if relative permeaility and saturation are given, 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 density and viscosity are given, 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 34 of file Q2PPiecewiseLinearSinkFlux.h.

Constructor & Destructor Documentation

◆ Q2PPiecewiseLinearSinkFlux()

Q2PPiecewiseLinearSinkFlux::Q2PPiecewiseLinearSinkFlux ( const InputParameters & parameters )

Definition at line 59 of file Q2PPiecewiseLinearSinkFlux.C.

   : SideIntegralPostprocessor(parameters),
     _sink_func(getParam<std::vector<Real>>("pressures"),
                getParam<std::vector<Real>>("bare_fluxes")),
     _m_func(getFunction("multiplying_fcn")),
     _pp(coupledValue("porepressure")),
     _use_mobility(isParamValid("fluid_density") && isParamValid("fluid_viscosity")),
     _use_relperm(isParamValid("fluid_relperm") && isCoupled("saturation")),
     _density(isParamValid("fluid_density") ? &getUserObject<RichardsDensity>("fluid_density")
                                            : NULL),
     _viscosity(isParamValid("fluid_viscosity") ? getParam<Real>("fluid_viscosity") : 1),
     _relperm(isParamValid("fluid_relperm") ? &getUserObject<RichardsRelPerm>("fluid_relperm")
                                            : NULL),
     _sat(isCoupled("saturation") ? coupledValue("saturation") : _zero),
     _permeability(getMaterialProperty<RealTensorValue>("permeability"))
 {
   if ((isParamValid("fluid_density") && !isParamValid("fluid_viscosity")) ||
       (!isParamValid("fluid_density") && isParamValid("fluid_viscosity")))
     mooseError("Q2PPiecewiseLinearSink: you must supply both of fluid_density and fluid_viscosity "
                "if you wish to multiply by the mobility");
   if ((isParamValid("fluid_relperm") && !isCoupled("saturation")) ||
       (!isParamValid("fluid_relperm") && isCoupled("saturation")))
     mooseError("Q2PPiecewiseLinearSink: you must supply both of fluid_relperm and saturation if "
                "you wish to multiply by the relative permeaility");
 }

Member Function Documentation

◆ computeQpIntegral()

Real Q2PPiecewiseLinearSinkFlux::computeQpIntegral ( )

protectedvirtual

Definition at line 86 of file Q2PPiecewiseLinearSinkFlux.C.

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