-fluid_mass_old/dt with the fluid mass being lumped to the nodes. More...

#include <Q2PNegativeNodalMassOld.h>

Inheritance diagram for Q2PNegativeNodalMassOld:

Public Member Functions
	Q2PNegativeNodalMassOld (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpResidual ()

Protected Attributes
const RichardsDensity &	_density

const VariableValue &	_other_var_nodal_old
	old value of the other variable (this is porepressure if the Variable is saturation) More...

bool	_var_is_pp
	whether the "other variable" is actually porepressure More...

const MaterialProperty< Real > &	_porosity_old
	value of the porosity at the start of the timestep More...

Detailed Description

-fluid_mass_old/dt with the fluid mass being lumped to the nodes.

Definition at line 26 of file Q2PNegativeNodalMassOld.h.

Constructor & Destructor Documentation

◆ Q2PNegativeNodalMassOld()

Q2PNegativeNodalMassOld::Q2PNegativeNodalMassOld ( const InputParameters & parameters )

Definition at line 43 of file Q2PNegativeNodalMassOld.C.

   : TimeKernel(parameters),
     _density(getUserObject<RichardsDensity>("fluid_density")),
     _other_var_nodal_old(coupledDofValuesOld("other_var")),
     _var_is_pp(getParam<bool>("var_is_porepressure")),
     _porosity_old(getMaterialProperty<Real>("porosity_old"))
 {
 }

Member Function Documentation

◆ computeQpResidual()

Real Q2PNegativeNodalMassOld::computeQpResidual ( )

protectedvirtual

Definition at line 53 of file Q2PNegativeNodalMassOld.C.

 {
   Real density_old;
   Real mass_old;
  
   if (_var_is_pp)
   {
     density_old = _density.density(_var.dofValuesOld()[_i]);
     mass_old = _porosity_old[_qp] * density_old * (1 - _other_var_nodal_old[_i]);
   }
   else
   {
     density_old = _density.density(_other_var_nodal_old[_i]);
     mass_old = _porosity_old[_qp] * density_old * _var.dofValuesOld()[_i];
   }
  
   return _test[_i][_qp] * (-mass_old) / _dt;
 }