Kernel = (mass - mass_old)/dt where mass = porosity*density*saturation This is used for the time derivative in Richards simulations Note that it is not lumped, so usually you want to use RichardsLumpedMassChange instead. More...

#include <RichardsMassChange.h>

Inheritance diagram for RichardsMassChange:

[legend]

Public Member Functions
	RichardsMassChange (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpResidual ()

virtual Real	computeQpJacobian ()

virtual Real	computeQpOffDiagJacobian (unsigned int jvar)

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

Protected Attributes
const RichardsVarNames &	_richards_name_UO
	holds info on the Richards variables More...

unsigned int	_pvar
	the Richards variable number eg, if richards name = 'pwater pgas poil', and this kernel is for pgas, then _pvar = 1 More...

bool	_use_supg
	whether to use SUPG for this kernel (not recommended) More...

const MaterialProperty< std::vector< Real > > &	_mass
	fluid mass (or fluid masses in multiphase) at quadpoints More...

const MaterialProperty< std::vector< std::vector< Real > > > &	_dmass
	d(fluid mass_i)/d(var_j) More...

const MaterialProperty< std::vector< Real > > &	_mass_old
	old value of fluid mass (or fluid masses in multiphase) at quadpoints More...

const MaterialProperty< std::vector< RealVectorValue > > &	_tauvel_SUPG
	tau_SUPG More...

const MaterialProperty< std::vector< std::vector< RealTensorValue > > > &	_dtauvel_SUPG_dgradv
	derivative of tau_SUPG wrt grad(variable) More...

const MaterialProperty< std::vector< std::vector< RealVectorValue > > > &	_dtauvel_SUPG_dv
	deriv of tau_SUPG wrt variable More...

Detailed Description

Kernel = (mass - mass_old)/dt where mass = porosity*density*saturation This is used for the time derivative in Richards simulations Note that it is not lumped, so usually you want to use RichardsLumpedMassChange instead.

Definition at line 27 of file RichardsMassChange.h.

Constructor & Destructor Documentation

◆ RichardsMassChange()

RichardsMassChange::RichardsMassChange ( const InputParameters & parameters )

Definition at line 35 of file RichardsMassChange.C.

   : TimeDerivative(parameters),
     _richards_name_UO(getUserObject<RichardsVarNames>("richardsVarNames_UO")),
     _pvar(_richards_name_UO.richards_var_num(_var.number())),
  
     _use_supg(getParam<bool>("use_supg")),
  
     _mass(getMaterialProperty<std::vector<Real>>("mass")),
     _dmass(getMaterialProperty<std::vector<std::vector<Real>>>("dmass")),
     _mass_old(getMaterialProperty<std::vector<Real>>("mass_old")),
  
     _tauvel_SUPG(getMaterialProperty<std::vector<RealVectorValue>>("tauvel_SUPG")),
     _dtauvel_SUPG_dgradv(
         getMaterialProperty<std::vector<std::vector<RealTensorValue>>>("dtauvel_SUPG_dgradv")),
     _dtauvel_SUPG_dv(
         getMaterialProperty<std::vector<std::vector<RealVectorValue>>>("dtauvel_SUPG_dv"))
 {
 }

Member Function Documentation

◆ computeQpJac()

Real RichardsMassChange::computeQpJac ( unsigned int wrt_num )

protected

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

Parameters

wrt_num take the derivative of the residual wrt this Richards variable

Definition at line 64 of file RichardsMassChange.C.

 {
   Real mass = _mass[_qp][_pvar];
   Real mass_old = _mass_old[_qp][_pvar];
   Real mass_prime = _phi[_j][_qp] * _dmass[_qp][_pvar][wrt_num];
  
   Real test_fcn = _test[_i][_qp];
   Real test_fcn_prime = 0;
  
   if (_use_supg)
   {
     test_fcn += _tauvel_SUPG[_qp][_pvar] * _grad_test[_i][_qp];
     test_fcn_prime +=
         _grad_phi[_j][_qp] * (_dtauvel_SUPG_dgradv[_qp][_pvar][wrt_num] * _grad_test[_i][_qp]) +
         _phi[_j][_qp] * _dtauvel_SUPG_dv[_qp][_pvar][wrt_num] * _grad_test[_i][_qp];
   }
   return (test_fcn * mass_prime + test_fcn_prime * (mass - mass_old)) / _dt;
 }

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

◆ computeQpJacobian()

Real RichardsMassChange::computeQpJacobian ( )

protectedvirtual

Definition at line 84 of file RichardsMassChange.C.

 {
   return computeQpJac(_pvar);
 }

◆ computeQpOffDiagJacobian()

Real RichardsMassChange::computeQpOffDiagJacobian ( unsigned int jvar )

protectedvirtual

Definition at line 90 of file RichardsMassChange.C.

 {
   if (_richards_name_UO.not_richards_var(jvar))
     return 0.0;
   unsigned int dvar = _richards_name_UO.richards_var_num(jvar);
   return computeQpJac(dvar);
 }

◆ computeQpResidual()

Real RichardsMassChange::computeQpResidual ( )

protectedvirtual

Definition at line 55 of file RichardsMassChange.C.

 {
   Real test_fcn = _test[_i][_qp];
   if (_use_supg)
     test_fcn += _tauvel_SUPG[_qp][_pvar] * _grad_test[_i][_qp];
   return test_fcn * (_mass[_qp][_pvar] - _mass_old[_qp][_pvar]) / _dt;
 }

Member Data Documentation

◆ _dmass

const MaterialProperty<std::vector<std::vector<Real> > >& RichardsMassChange::_dmass

protected

d(fluid mass_i)/d(var_j)

Definition at line 56 of file RichardsMassChange.h.

Referenced by computeQpJac().

◆ _dtauvel_SUPG_dgradv

const MaterialProperty<std::vector<std::vector<RealTensorValue> > >& RichardsMassChange::_dtauvel_SUPG_dgradv

protected

derivative of tau_SUPG wrt grad(variable)

Definition at line 65 of file RichardsMassChange.h.

Referenced by computeQpJac().

◆ _dtauvel_SUPG_dv

const MaterialProperty<std::vector<std::vector<RealVectorValue> > >& RichardsMassChange::_dtauvel_SUPG_dv

protected

deriv of tau_SUPG wrt variable

Definition at line 68 of file RichardsMassChange.h.

Referenced by computeQpJac().

◆ _mass

const MaterialProperty<std::vector<Real> >& RichardsMassChange::_mass

protected

fluid mass (or fluid masses in multiphase) at quadpoints

Definition at line 53 of file RichardsMassChange.h.

Referenced by computeQpJac(), and computeQpResidual().

◆ _mass_old

const MaterialProperty<std::vector<Real> >& RichardsMassChange::_mass_old

protected

old value of fluid mass (or fluid masses in multiphase) at quadpoints

Definition at line 59 of file RichardsMassChange.h.

Referenced by computeQpJac(), and computeQpResidual().

◆ _pvar

unsigned int RichardsMassChange::_pvar

protected

the Richards variable number eg, if richards name = 'pwater pgas poil', and this kernel is for pgas, then _pvar = 1

Definition at line 47 of file RichardsMassChange.h.

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

◆ _richards_name_UO

const RichardsVarNames& RichardsMassChange::_richards_name_UO

protected

holds info on the Richards variables

Definition at line 40 of file RichardsMassChange.h.

Referenced by computeQpOffDiagJacobian().

◆ _tauvel_SUPG

const MaterialProperty<std::vector<RealVectorValue> >& RichardsMassChange::_tauvel_SUPG

protected

tau_SUPG

Definition at line 62 of file RichardsMassChange.h.

Referenced by computeQpJac(), and computeQpResidual().

◆ _use_supg

bool RichardsMassChange::_use_supg

protected

whether to use SUPG for this kernel (not recommended)

Definition at line 50 of file RichardsMassChange.h.

Referenced by computeQpJac(), and computeQpResidual().

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

richards/include/kernels/RichardsMassChange.h
richards/src/kernels/RichardsMassChange.C

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ RichardsMassChange()

Member Function Documentation

◆ computeQpJac()

◆ computeQpJacobian()

◆ computeQpOffDiagJacobian()

◆ computeQpResidual()

Member Data Documentation

◆ _dmass

◆ _dtauvel_SUPG_dgradv

◆ _dtauvel_SUPG_dv

◆ _mass

◆ _mass_old

◆ _pvar

◆ _richards_name_UO

◆ _tauvel_SUPG

◆ _use_supg