PorousFlowEffectiveStressCoupling computes -coefficient*effective_porepressure*grad_component(test) where component is the spatial component (not a fluid component!) More...

#include <PorousFlowEffectiveStressCoupling.h>

Inheritance diagram for PorousFlowEffectiveStressCoupling:

[legend]

Public Member Functions
	PorousFlowEffectiveStressCoupling (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpResidual () override

virtual Real	computeQpJacobian () override

virtual Real	computeQpOffDiagJacobian (unsigned int jvar) override

Protected Attributes
const PorousFlowDictator &	_dictator
	The PorousFlow dictator that holds global info about the simulation. More...

const Real	_coefficient
	Biot coefficient. More...

const unsigned int	_component
	The spatial component. More...

const MaterialProperty< Real > &	_pf
	Effective porepressure. More...

const MaterialProperty< std::vector< Real > > &	_dpf_dvar
	d(effective porepressure)/(d porflow variable) More...

const bool	_rz
	Whether an RZ coordinate system is being used. More...

Detailed Description

PorousFlowEffectiveStressCoupling computes -coefficient*effective_porepressure*grad_component(test) where component is the spatial component (not a fluid component!)

Definition at line 26 of file PorousFlowEffectiveStressCoupling.h.

Constructor & Destructor Documentation

◆ PorousFlowEffectiveStressCoupling()

PorousFlowEffectiveStressCoupling::PorousFlowEffectiveStressCoupling ( const InputParameters & parameters )

Definition at line 34 of file PorousFlowEffectiveStressCoupling.C.

   : Kernel(parameters),
     _dictator(getUserObject<PorousFlowDictator>("PorousFlowDictator")),
     _coefficient(getParam<Real>("biot_coefficient")),
     _component(getParam<unsigned int>("component")),
     _pf(getMaterialProperty<Real>("PorousFlow_effective_fluid_pressure_qp")),
     _dpf_dvar(
         getMaterialProperty<std::vector<Real>>("dPorousFlow_effective_fluid_pressure_qp_dvar")),
     _rz(getBlockCoordSystem() == Moose::COORD_RZ)
 {
   if (_component >= _mesh.dimension())
     paramError("component", "The component cannot be greater than the mesh dimension");
 }

Member Function Documentation

◆ computeQpJacobian()

Real PorousFlowEffectiveStressCoupling::computeQpJacobian ( )

overrideprotectedvirtual

Definition at line 58 of file PorousFlowEffectiveStressCoupling.C.

 {
   if (_dictator.notPorousFlowVariable(_var.number()))
     return 0.0;
   const unsigned int pvar = _dictator.porousFlowVariableNum(_var.number());
   if (_rz && _component == 0)
     return -_coefficient * _phi[_j][_qp] * _dpf_dvar[_qp][pvar] *
            (_grad_test[_i][_qp](0) + _test[_i][_qp] / _q_point[_qp](0));
   return -_coefficient * _phi[_j][_qp] * _dpf_dvar[_qp][pvar] * _grad_test[_i][_qp](_component);
 }

◆ computeQpOffDiagJacobian()

Real PorousFlowEffectiveStressCoupling::computeQpOffDiagJacobian ( unsigned int jvar )

overrideprotectedvirtual

Definition at line 70 of file PorousFlowEffectiveStressCoupling.C.

 {
   if (_dictator.notPorousFlowVariable(jvar))
     return 0.0;
   const unsigned int pvar = _dictator.porousFlowVariableNum(jvar);
   if (_rz && _component == 0)
     return -_coefficient * _phi[_j][_qp] * _dpf_dvar[_qp][pvar] *
            (_grad_test[_i][_qp](0) + _test[_i][_qp] / _q_point[_qp](0));
   return -_coefficient * _phi[_j][_qp] * _dpf_dvar[_qp][pvar] * _grad_test[_i][_qp](_component);
 }

◆ computeQpResidual()

Real PorousFlowEffectiveStressCoupling::computeQpResidual ( )

overrideprotectedvirtual

Definition at line 50 of file PorousFlowEffectiveStressCoupling.C.

 {
   if (_rz && _component == 0)
     return -_coefficient * _pf[_qp] * (_grad_test[_i][_qp](0) + _test[_i][_qp] / _q_point[_qp](0));
   return -_coefficient * _pf[_qp] * _grad_test[_i][_qp](_component);
 }