Q2PNodalMass Class Reference

fluid_mass/dt lumped to the nodes More...

#include <Q2PNodalMass.h>

Inheritance diagram for Q2PNodalMass:

Public Member Functions
	Q2PNodalMass (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpResidual ()
virtual Real	computeQpJacobian ()
virtual Real	computeQpOffDiagJacobian (unsigned int jvar)

Protected Attributes
const RichardsDensity &	_density
const VariableValue &	_other_var_nodal
	the other variable (this is porepressure if the Variable is saturation) More...
unsigned int	_other_var_num
	variable number of the other variable More...
bool	_var_is_pp
	whether the "other variable" is actually porepressure More...
const MaterialProperty< Real > &	_porosity
	current value of the porosity More...

Detailed Description

fluid_mass/dt lumped to the nodes

Definition at line 25 of file Q2PNodalMass.h.

Constructor & Destructor Documentation

Q2PNodalMass()

Q2PNodalMass::Q2PNodalMass

(

const InputParameters &

parameters

)

Definition at line 43 of file Q2PNodalMass.C.

  : TimeKernel(parameters),
    _density(getUserObject<RichardsDensity>("fluid_density")),
    _other_var_nodal(coupledDofValues("other_var")),
    _other_var_num(coupled("other_var")),
    _var_is_pp(getParam<bool>("var_is_porepressure")),
    _porosity(getMaterialProperty<Real>("porosity"))
{
}

Member Function Documentation

computeQpJacobian()

Real Q2PNodalMass::computeQpJacobian ( )

protectedvirtual

Definition at line 74 of file Q2PNodalMass.C.

{
  if (_i != _j)
    return 0.0;
 
  Real mass_prime;
 
  if (_var_is_pp)
  {
    // we're calculating the derivative wrt porepressure
    Real ddensity = _density.ddensity(_var.dofValues()[_i]);
    mass_prime = _porosity[_qp] * ddensity * (1 - _other_var_nodal[_i]);
  }
  else
  {
    // we're calculating the deriv wrt saturation
    Real density = _density.density(_other_var_nodal[_i]);
    mass_prime = _porosity[_qp] * density;
  }
 
  return _test[_i][_qp] * mass_prime / _dt;
}

computeQpOffDiagJacobian()

Real Q2PNodalMass::computeQpOffDiagJacobian ( unsigned int  jvar )

protectedvirtual

Definition at line 98 of file Q2PNodalMass.C.

{
  if (jvar != _other_var_num)
    return 0.0;
  if (_i != _j)
    return 0.0;
 
  Real mass_prime;
 
  if (_var_is_pp)
  {
    // we're calculating the deriv wrt saturation
    Real density = _density.density(_var.dofValues()[_i]);
    mass_prime = -_porosity[_qp] * density;
  }
  else
  {
    // we're calculating the deriv wrt porepressure
    Real ddensity = _density.ddensity(_other_var_nodal[_i]);
    mass_prime = _porosity[_qp] * ddensity * _var.dofValues()[_i];
  }
 
  return _test[_i][_qp] * mass_prime / _dt;
}

computeQpResidual()

Real Q2PNodalMass::computeQpResidual ( )

protectedvirtual

Definition at line 54 of file Q2PNodalMass.C.

{
  Real density;
  Real mass;
 
  if (_var_is_pp)
  {
    density = _density.density(_var.dofValues()[_i]);
    mass = _porosity[_qp] * density * (1 - _other_var_nodal[_i]);
  }
  else
  {
    density = _density.density(_other_var_nodal[_i]);
    mass = _porosity[_qp] * density * _var.dofValues()[_i];
  }
 
  return _test[_i][_qp] * mass / _dt;
}

Member Data Documentation

_density

const RichardsDensity& Q2PNodalMass::_density

protected

Definition at line 37 of file Q2PNodalMass.h.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), and computeQpResidual().

_other_var_nodal

const VariableValue& Q2PNodalMass::_other_var_nodal

protected

the other variable (this is porepressure if the Variable is saturation)

Definition at line 40 of file Q2PNodalMass.h.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), and computeQpResidual().

_other_var_num

unsigned int Q2PNodalMass::_other_var_num

protected

variable number of the other variable

Definition at line 43 of file Q2PNodalMass.h.

Referenced by computeQpOffDiagJacobian().

_porosity

const MaterialProperty<Real>& Q2PNodalMass::_porosity

protected

current value of the porosity

Definition at line 49 of file Q2PNodalMass.h.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), and computeQpResidual().

_var_is_pp

bool Q2PNodalMass::_var_is_pp

protected

whether the "other variable" is actually porepressure

Definition at line 46 of file Q2PNodalMass.h.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), and computeQpResidual().

---

The documentation for this class was generated from the following files:

• richards/include/kernels/Q2PNodalMass.h
• richards/src/kernels/Q2PNodalMass.C