89 std::vector<Real>
_pp;
110 Real jac(
unsigned int wrt_num);
virtual void computeOffDiagJacobian(unsigned int jvar) override
void prepareNodalValues()
calculates the nodal values of pressure, mobility, and derivatives thereof
virtual void computeJacobian() override
virtual void computeResidual() override
std::vector< Real > _nodal_density
nodal values of fluid density
LinearInterpolation _sink_func
piecewise-linear function of porepressure (this defines the strength of the sink) ...
bool _use_mobility
whether to multiply the sink flux by permeability*density/viscosity
Q2PPiecewiseLinearSink(const InputParameters ¶meters)
Base class for Richards relative permeability classes that provide relative permeability as a functio...
Real jac(unsigned int wrt_num)
derivative of residual wrt the wrt_num variable
unsigned int _other_var_num
the variable number of the other variable
static InputParameters validParams()
Applies a fully-upwinded flux sink to a boundary The sink is a piecewise linear function of porepress...
virtual Real computeQpResidual() override
bool _var_is_pp
whether the Variable for this BC is porepressure or not
virtual Real computeQpOffDiagJacobian(unsigned int jvar) override
std::vector< Real > _dnodal_density_dp
d(_nodal_density)/d(porepressure)
std::vector< Real > _pp
nodal values of porepressure
const RichardsRelPerm & _relperm
fluid relative permeability
OutputTools< Real >::VariableValue VariableValue
virtual Real computeQpJacobian() override
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
const MaterialProperty< RealTensorValue > & _permeability
permeability
bool _use_relperm
whether to multiply the sink flux by relative permeability
const InputParameters & parameters() const
std::vector< Real > _dnodal_relperm_ds
d(_nodal_relperm)/d(saturation)
Base class for fluid density as a function of porepressure The functions density, ddensity and d2dens...
unsigned int _num_nodes
number of nodes in this element.
std::vector< Real > _nodal_relperm
nodal values of relative permeability
std::vector< Real > _sat
nodal values of saturation
const RichardsDensity & _density
fluid density
const Function & _m_func
sink flux gets multiplied by this function
const VariableValue & _other_var_nodal
the other variable in the 2-phase system (this is saturation if Variable=porepressure, and viceversa)