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:

[legend]

Public Member Functions
	PoroFullSatTimeDerivative (const InputParameters &parameters)

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 26 of file PoroFullSatTimeDerivative.h.

Constructor & Destructor Documentation

◆ PoroFullSatTimeDerivative()

PoroFullSatTimeDerivative::PoroFullSatTimeDerivative ( const InputParameters & parameters )

Definition at line 29 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 60 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 68 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 51 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;
 }