PoroFullSatTimeDerivative Class Reference

Kernel = biot_coefficient*d(volumetric_strain)/dt + (1/biot_modulus)*d(porepressure)/dt this is the time-derivative for poromechanics for a single-phase, fully-saturated fluid with constant bulk modulus. More...

#include <PoroFullSatTimeDerivative.h>

Inheritance diagram for PoroFullSatTimeDerivative:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	PoroFullSatTimeDerivative (const InputParameters &parameters)
const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialProperty (const std::string &name)
const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialPropertyByName (const std::string &name)
void	validateDerivativeMaterialPropertyBase (const std::string &base)
const MaterialPropertyName	derivativePropertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const
const MaterialPropertyName	derivativePropertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const
const MaterialPropertyName	derivativePropertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const
const MaterialPropertyName	derivativePropertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const
GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)
const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)
void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)
void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const
const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const
const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const
const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const
const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const
const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const
const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const
const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

Static Public Member Functions
static InputParameters	validParams ()

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

Private Attributes
const VariableValue &	_u_old
	old value of porepressure More...
const MaterialProperty< Real > &	_volstrain
	volumetric strain More...
const MaterialProperty< Real > &	_volstrain_old
	old value of volumetric strain More...
unsigned int	_ndisp
	number of displacement variables More...
std::vector< unsigned int >	_disp_var_num
	variable number of the displacements variables More...
const MaterialProperty< Real > &	_alpha
	Biot coefficient. More...
const MaterialProperty< Real > &	_one_over_biot_modulus
	1/M, where M is the Biot modulus More...
const MaterialProperty< Real > &	_done_over_biot_modulus_dP
	d(1/M)/d(porepressure) More...
const MaterialProperty< Real > &	_done_over_biot_modulus_dep
	d(1/M)/d(volumetric strain) More...

Detailed Description

Kernel = biot_coefficient*d(volumetric_strain)/dt + (1/biot_modulus)*d(porepressure)/dt this is the time-derivative for poromechanics for a single-phase, fully-saturated fluid with constant bulk modulus.

Definition at line 22 of file PoroFullSatTimeDerivative.h.

Constructor & Destructor Documentation

PoroFullSatTimeDerivative()

PoroFullSatTimeDerivative::PoroFullSatTimeDerivative

(

const InputParameters &

parameters

)

Definition at line 28 of file PoroFullSatTimeDerivative.C.

  : DerivativeMaterialInterface<TimeDerivative>(parameters),
    _u_old(valueOld()),
    _volstrain(getMaterialProperty<Real>("volumetric_strain")),
    _volstrain_old(getMaterialPropertyOld<Real>("volumetric_strain")),

    _ndisp(coupledComponents("displacements")),
    _disp_var_num(_ndisp),

    _alpha(getMaterialProperty<Real>("biot_coefficient")),

    _one_over_biot_modulus(getMaterialProperty<Real>("one_over_biot_modulus")),
    _done_over_biot_modulus_dP(
        getMaterialPropertyDerivative<Real>("one_over_biot_modulus", _var.name())),
    _done_over_biot_modulus_dep(
        getMaterialPropertyDerivative<Real>("one_over_biot_modulus", "volumetric_strain"))
{
  for (unsigned i = 0; i < _ndisp; ++i)
    _disp_var_num[i] = coupled("displacements", i);
}

Member Function Documentation

computeQpJacobian()

Real PoroFullSatTimeDerivative::computeQpJacobian ( )

protectedvirtual

Definition at line 59 of file PoroFullSatTimeDerivative.C.

{
  Real jac = _one_over_biot_modulus[_qp] * _phi[_j][_qp];
  jac += _done_over_biot_modulus_dP[_qp] * _phi[_j][_qp] * (_u[_qp] - _u_old[_qp]);
  return _test[_i][_qp] * jac / _dt;
}

computeQpOffDiagJacobian()

Real PoroFullSatTimeDerivative::computeQpOffDiagJacobian ( unsigned int  jvar )

protectedvirtual

Definition at line 67 of file PoroFullSatTimeDerivative.C.

{
  Real jac = 0;
  for (unsigned i = 0; i < _ndisp; ++i)
    if (jvar == _disp_var_num[i])
      jac = _grad_phi[_j][_qp](i);

  jac *= _done_over_biot_modulus_dep[_qp] * (_u[_qp] - _u_old[_qp]) + _alpha[_qp];

  return _test[_i][_qp] * jac / _dt;
}

computeQpResidual()

Real PoroFullSatTimeDerivative::computeQpResidual ( )

protectedvirtual

Definition at line 50 of file PoroFullSatTimeDerivative.C.

{
  // here, "_u" is the porepressure
  Real res = _one_over_biot_modulus[_qp] * (_u[_qp] - _u_old[_qp]);
  res += _alpha[_qp] * (_volstrain[_qp] - _volstrain_old[_qp]);
  return _test[_i][_qp] * res / _dt;
}

validParams()

InputParameters PoroFullSatTimeDerivative::validParams ( )

static

Definition at line 15 of file PoroFullSatTimeDerivative.C.

{
  InputParameters params = TimeDerivative::validParams();
  params.addRequiredCoupledVar(
      "displacements",
      "The displacements appropriate for the simulation geometry and coordinate system");
  params.addClassDescription("Kernel = biot_coefficient*d(volumetric_strain)/dt + "
                             "(1/biot_modulus)*d(porepressure)/dt.  This is the time-derivative "
                             "for poromechanics for a single-phase, fully-saturated fluid with "
                             "constant bulk modulus");
  return params;
}

Member Data Documentation

_alpha

const MaterialProperty<Real>& PoroFullSatTimeDerivative::_alpha

private

Biot coefficient.

Definition at line 53 of file PoroFullSatTimeDerivative.h.

Referenced by computeQpOffDiagJacobian(), and computeQpResidual().

_disp_var_num

std::vector<unsigned int> PoroFullSatTimeDerivative::_disp_var_num

private

variable number of the displacements variables

Definition at line 50 of file PoroFullSatTimeDerivative.h.

Referenced by computeQpOffDiagJacobian(), and PoroFullSatTimeDerivative().

_done_over_biot_modulus_dep

const MaterialProperty<Real>& PoroFullSatTimeDerivative::_done_over_biot_modulus_dep

private

d(1/M)/d(volumetric strain)

Definition at line 62 of file PoroFullSatTimeDerivative.h.

Referenced by computeQpOffDiagJacobian().

_done_over_biot_modulus_dP

const MaterialProperty<Real>& PoroFullSatTimeDerivative::_done_over_biot_modulus_dP

private

d(1/M)/d(porepressure)

Definition at line 59 of file PoroFullSatTimeDerivative.h.

Referenced by computeQpJacobian().

_ndisp

unsigned int PoroFullSatTimeDerivative::_ndisp

private

number of displacement variables

Definition at line 47 of file PoroFullSatTimeDerivative.h.

Referenced by computeQpOffDiagJacobian(), and PoroFullSatTimeDerivative().

_one_over_biot_modulus

const MaterialProperty<Real>& PoroFullSatTimeDerivative::_one_over_biot_modulus

private

1/M, where M is the Biot modulus

Definition at line 56 of file PoroFullSatTimeDerivative.h.

Referenced by computeQpJacobian(), and computeQpResidual().

_u_old

const VariableValue& PoroFullSatTimeDerivative::_u_old

private

old value of porepressure

Definition at line 38 of file PoroFullSatTimeDerivative.h.

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

_volstrain

const MaterialProperty<Real>& PoroFullSatTimeDerivative::_volstrain

private

volumetric strain

Definition at line 41 of file PoroFullSatTimeDerivative.h.

Referenced by computeQpResidual().

_volstrain_old

const MaterialProperty<Real>& PoroFullSatTimeDerivative::_volstrain_old

private

old value of volumetric strain

Definition at line 44 of file PoroFullSatTimeDerivative.h.

Referenced by computeQpResidual().

---

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

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