Q2PSaturationDiffusion Class Reference

Diffusive Kernel that models nonzero capillary pressure in Q2P models The Variable of this Kernel should be the saturation. More...

#include <Q2PSaturationDiffusion.h>

Inheritance diagram for Q2PSaturationDiffusion:

Public Member Functions
	Q2PSaturationDiffusion (const InputParameters &parameters)

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

Protected Attributes
const RichardsDensity &	_density
	water density More...
const RichardsRelPerm &	_relperm
	water relative permeability More...
const VariableValue &	_pp
	porepressure at the quadpoints More...
unsigned int	_pp_var_num
	variable number of the porepressure variable More...
Real	_viscosity
	fluid viscosity More...
const MaterialProperty< RealTensorValue > &	_permeability
	permeability More...
Real	_diffusivity

Detailed Description

Diffusive Kernel that models nonzero capillary pressure in Q2P models The Variable of this Kernel should be the saturation.

Definition at line 28 of file Q2PSaturationDiffusion.h.

Constructor & Destructor Documentation

Q2PSaturationDiffusion()

Q2PSaturationDiffusion::Q2PSaturationDiffusion

(

const InputParameters &

parameters

)

Definition at line 35 of file Q2PSaturationDiffusion.C.

  : Kernel(parameters),
    _density(getUserObject<RichardsDensity>("fluid_density")),
    _relperm(getUserObject<RichardsRelPerm>("fluid_relperm")),
    _pp(coupledValue("porepressure_variable")),
    _pp_var_num(coupled("porepressure_variable")),
    _viscosity(getParam<Real>("fluid_viscosity")),
    _permeability(getMaterialProperty<RealTensorValue>("permeability")),
    _diffusivity(getParam<Real>("diffusivity"))
{
}

Member Function Documentation

computeQpJacobian()

Real Q2PSaturationDiffusion::computeQpJacobian ( )

protectedvirtual

Definition at line 55 of file Q2PSaturationDiffusion.C.

{
  Real coef = _diffusivity * _relperm.relperm(_u[_qp]) * _density.density(_pp[_qp]) / _viscosity;
  Real coefp = _diffusivity * _relperm.drelperm(_u[_qp]) * _density.density(_pp[_qp]) / _viscosity;
  return coefp * _phi[_j][_qp] * _grad_test[_i][_qp] * (_permeability[_qp] * _grad_u[_qp]) +
         coef * _grad_test[_i][_qp] * (_permeability[_qp] * _grad_phi[_j][_qp]);
}

computeQpOffDiagJacobian()

Real Q2PSaturationDiffusion::computeQpOffDiagJacobian ( unsigned int  jvar )

protectedvirtual

Definition at line 64 of file Q2PSaturationDiffusion.C.

{
  if (jvar != _pp_var_num)
    return 0.0;
  Real coefp = _diffusivity * _relperm.relperm(_u[_qp]) * _density.ddensity(_pp[_qp]) / _viscosity;
  return coefp * _phi[_j][_qp] * (_grad_test[_i][_qp] * (_permeability[_qp] * _grad_u[_qp]));
}

computeQpResidual()

Real Q2PSaturationDiffusion::computeQpResidual ( )

protectedvirtual

Definition at line 48 of file Q2PSaturationDiffusion.C.

{
  Real coef = _diffusivity * _relperm.relperm(_u[_qp]) * _density.density(_pp[_qp]) / _viscosity;
  return coef * _grad_test[_i][_qp] * (_permeability[_qp] * _grad_u[_qp]);
}

Member Data Documentation

_density

const RichardsDensity& Q2PSaturationDiffusion::_density

protected

water density

Definition at line 41 of file Q2PSaturationDiffusion.h.

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

_diffusivity

Real Q2PSaturationDiffusion::_diffusivity

protected

Definition at line 58 of file Q2PSaturationDiffusion.h.

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

_permeability

const MaterialProperty<RealTensorValue>& Q2PSaturationDiffusion::_permeability

protected

permeability

Definition at line 56 of file Q2PSaturationDiffusion.h.

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

_pp

const VariableValue& Q2PSaturationDiffusion::_pp

protected

porepressure at the quadpoints

Definition at line 47 of file Q2PSaturationDiffusion.h.

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

_pp_var_num

unsigned int Q2PSaturationDiffusion::_pp_var_num

protected

variable number of the porepressure variable

Definition at line 50 of file Q2PSaturationDiffusion.h.

Referenced by computeQpOffDiagJacobian().

_relperm

const RichardsRelPerm& Q2PSaturationDiffusion::_relperm

protected

water relative permeability

Definition at line 44 of file Q2PSaturationDiffusion.h.

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

_viscosity

Real Q2PSaturationDiffusion::_viscosity

protected

fluid viscosity

Definition at line 53 of file Q2PSaturationDiffusion.h.

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

---

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

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