Kernel = (desorped_mass - desorped_mass_old)/dt It is NOT lumped to the nodes. More...

#include <PorousFlowDesorpedMassTimeDerivative.h>

Inheritance diagram for PorousFlowDesorpedMassTimeDerivative:

Public Member Functions
	PorousFlowDesorpedMassTimeDerivative (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpResidual () override

virtual Real	computeQpJacobian () override

virtual Real	computeQpOffDiagJacobian (unsigned int jvar) override

Real	computeQpJac (unsigned int jvar) const
	Derivative of residual with respect to variable number jvar This is used by both computeQpJacobian and computeQpOffDiagJacobian. More...

Protected Attributes
const PorousFlowDictator &	_dictator
	PorousFlowDictator UserObject. More...

const unsigned int	_conc_var_number
	The MOOSE variable number of the concentration variable. More...

const VariableValue &	_conc
	The concentration variable. More...

const VariableValue &	_conc_old
	Old value of the concentration variable. More...

const MaterialProperty< Real > &	_porosity
	Porosity at the qps. More...

const MaterialProperty< Real > &	_porosity_old
	Old value of porosity. More...

const MaterialProperty< std::vector< Real > > &	_dporosity_dvar
	d(porosity)/d(PorousFlow variable) - these derivatives will be wrt variables at the qps More...

const MaterialProperty< std::vector< RealGradient > > &	_dporosity_dgradvar
	d(porosity)/d(grad PorousFlow variable) - these derivatives will be wrt grad(vars) at qps More...

Detailed Description

Kernel = (desorped_mass - desorped_mass_old)/dt It is NOT lumped to the nodes.

Definition at line 24 of file PorousFlowDesorpedMassTimeDerivative.h.

Constructor & Destructor Documentation

◆ PorousFlowDesorpedMassTimeDerivative()

PorousFlowDesorpedMassTimeDerivative::PorousFlowDesorpedMassTimeDerivative ( const InputParameters & parameters )

Definition at line 31 of file PorousFlowDesorpedMassTimeDerivative.C.

   : TimeKernel(parameters),
     _dictator(getUserObject<PorousFlowDictator>("PorousFlowDictator")),
     _conc_var_number(coupled("conc_var")),
     _conc(coupledValue("conc_var")),
     _conc_old(coupledValueOld("conc_var")),
     _porosity(getMaterialProperty<Real>("PorousFlow_porosity_qp")),
     _porosity_old(getMaterialPropertyOld<Real>("PorousFlow_porosity_qp")),
     _dporosity_dvar(getMaterialProperty<std::vector<Real>>("dPorousFlow_porosity_qp_dvar")),
     _dporosity_dgradvar(
         getMaterialProperty<std::vector<RealGradient>>("dPorousFlow_porosity_qp_dgradvar"))
 {
 }

Member Function Documentation

◆ computeQpJac()

Real PorousFlowDesorpedMassTimeDerivative::computeQpJac ( unsigned int jvar ) const

protected

Derivative of residual with respect to variable number jvar This is used by both computeQpJacobian and computeQpOffDiagJacobian.

Parameters

jvar	take the derivative of the residual wrt this Moose variable

Definition at line 67 of file PorousFlowDesorpedMassTimeDerivative.C.

 {
   Real deriv = 0.0;
  
   if (jvar == _conc_var_number)
     deriv = (1.0 - _porosity[_qp]) * _phi[_j][_qp];
  
   if (_dictator.notPorousFlowVariable(jvar))
     return _test[_i][_qp] * deriv / _dt;
   const unsigned int pvar = _dictator.porousFlowVariableNum(jvar);
  
   deriv -= _dporosity_dgradvar[_qp][pvar] * _grad_phi[_j][_qp] * _conc[_qp];
   deriv -= _dporosity_dvar[_qp][pvar] * _phi[_j][_qp] * _conc[_qp];
  
   return _test[_i][_qp] * deriv / _dt;
 }

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

◆ computeQpJacobian()

Real PorousFlowDesorpedMassTimeDerivative::computeQpJacobian ( )

overrideprotectedvirtual

Definition at line 55 of file PorousFlowDesorpedMassTimeDerivative.C.

 {
   return computeQpJac(_var.number());
 }

◆ computeQpOffDiagJacobian()

Real PorousFlowDesorpedMassTimeDerivative::computeQpOffDiagJacobian ( unsigned int jvar )

overrideprotectedvirtual

Definition at line 61 of file PorousFlowDesorpedMassTimeDerivative.C.

 {
   return computeQpJac(jvar);
 }

◆ computeQpResidual()

Real PorousFlowDesorpedMassTimeDerivative::computeQpResidual ( )

overrideprotectedvirtual

Definition at line 47 of file PorousFlowDesorpedMassTimeDerivative.C.

 {
   Real c = (1.0 - _porosity[_qp]) * _conc[_qp];
   Real c_old = (1.0 - _porosity_old[_qp]) * _conc_old[_qp];
   return _test[_i][_qp] * (c - c_old) / _dt;
 }