Material designed to calculate and store all the quantities needed for the fluid-flow part of poromechanics, assuming a fully-saturated, single-phase fluid with constant bulk modulus. More...

#include <PoroFullSatMaterial.h>

Inheritance diagram for PoroFullSatMaterial:

[legend]

Public Member Functions
	PoroFullSatMaterial (const InputParameters &parameters)

Protected Member Functions
virtual void	initQpStatefulProperties ()

virtual void	computeQpProperties ()

Protected Attributes
Real	_phi0
	porosity at zero porepressure and volumetric strain More...

Real	_alpha
	Biot coefficient. More...

Real	_one_over_K
	1/K, where K is the solid bulk modulus. Usually 1/K = C_iijj, where C is the compliance matrix: strain_ij = C_ijkl stress_kl More...

Real	_one_over_Kf
	1/Kf, where Kf is the fluid bulk modulus. More...

bool	_constant_porosity
	whether to use constant porosity (set _porosity = _phi0 always) More...

const VariableValue &	_porepressure
	porepressure variable More...

std::string	_porepressure_name
	name given by user to the porepressure variable More...

unsigned int	_ndisp
	number of displacement variables supplied More...

std::vector< const VariableGradient * >	_grad_disp
	grad(displacement) More...

MaterialProperty< Real > &	_vol_strain
	volumetric strain = strain_ii More...

MaterialProperty< Real > &	_biot_coefficient
	Biot coefficient. More...

MaterialProperty< Real > &	_porosity
	porosity More...

MaterialProperty< Real > &	_dporosity_dP
	d(porosity)/d(porepressure) More...

MaterialProperty< Real > &	_dporosity_dep
	d(porosity)/d(volumetric_strain) More...

MaterialProperty< Real > &	_one_over_biot_modulus
	1/M, where M is the Biot modulus More...

MaterialProperty< Real > &	_done_over_biot_modulus_dP
	d(1/M)/d(porepressure) More...

MaterialProperty< Real > &	_done_over_biot_modulus_dep
	d(1/M)/d(volumetric_strain) More...

Detailed Description

Material designed to calculate and store all the quantities needed for the fluid-flow part of poromechanics, assuming a fully-saturated, single-phase fluid with constant bulk modulus.

Definition at line 27 of file PoroFullSatMaterial.h.

Constructor & Destructor Documentation

◆ PoroFullSatMaterial()

PoroFullSatMaterial::PoroFullSatMaterial ( const InputParameters & parameters )

Definition at line 45 of file PoroFullSatMaterial.C.

   : DerivativeMaterialInterface<Material>(parameters),
  
     _phi0(getParam<Real>("porosity0")),
     _alpha(getParam<Real>("biot_coefficient")),
     _one_over_K(getParam<Real>("solid_bulk_compliance")),
     _one_over_Kf(getParam<Real>("fluid_bulk_compliance")),
     _constant_porosity(getParam<bool>("constant_porosity")),
  
     _porepressure(coupledValue("porepressure")),
     _porepressure_name(getVar("porepressure", 0)->name()),
  
     _ndisp(coupledComponents("displacements")),
     _grad_disp(_ndisp),
  
     _vol_strain(declareProperty<Real>("volumetric_strain")),
  
     _biot_coefficient(declareProperty<Real>("biot_coefficient")),
  
     _porosity(declareProperty<Real>("porosity")),
     _dporosity_dP(declarePropertyDerivative<Real>("porosity", _porepressure_name)),
     _dporosity_dep(declarePropertyDerivative<Real>("porosity", "volumetric_strain")),
  
     _one_over_biot_modulus(declareProperty<Real>("one_over_biot_modulus")),
     _done_over_biot_modulus_dP(
         declarePropertyDerivative<Real>("one_over_biot_modulus", _porepressure_name)),
     _done_over_biot_modulus_dep(
         declarePropertyDerivative<Real>("one_over_biot_modulus", "volumetric_strain"))
 {
   for (unsigned int i = 0; i < _ndisp; ++i)
     _grad_disp[i] = &coupledGradient("displacements", i);
 }

Member Function Documentation

◆ computeQpProperties()

void PoroFullSatMaterial::computeQpProperties ( )

protectedvirtual

Definition at line 85 of file PoroFullSatMaterial.C.

 {
   _biot_coefficient[_qp] = _alpha;
  
   _vol_strain[_qp] = 0;
   for (unsigned i = 0; i < _ndisp; ++i)
     _vol_strain[_qp] += (*_grad_disp[i])[_qp](i); // cartesian coordinates?
  
   if (_constant_porosity)
   {
     _porosity[_qp] = _phi0;
     _dporosity_dP[_qp] = 0;
     _dporosity_dep[_qp] = 0;
  
     _one_over_biot_modulus[_qp] =
         (1 - _alpha) * (_alpha - _porosity[_qp]) * _one_over_K + _porosity[_qp] * _one_over_Kf;
     _done_over_biot_modulus_dP[_qp] = 0;
     _done_over_biot_modulus_dep[_qp] = 0;
   }
   else
   {
     _porosity[_qp] = _alpha +
                      (_phi0 - _alpha) * std::exp(-(1 - _alpha) * _one_over_K * _porepressure[_qp] -
                                                  _vol_strain[_qp]);
     _dporosity_dP[_qp] =
         (_phi0 - _alpha) * (_alpha - 1) * _one_over_K *
         std::exp(-(1 - _alpha) * _one_over_K * _porepressure[_qp] - _vol_strain[_qp]);
     _dporosity_dep[_qp] =
         -(_phi0 - _alpha) *
         std::exp(-(1 - _alpha) * _one_over_K * _porepressure[_qp] - _vol_strain[_qp]);
  
     _one_over_biot_modulus[_qp] =
         (1 - _alpha) * (_alpha - _porosity[_qp]) * _one_over_K + _porosity[_qp] * _one_over_Kf;
     _done_over_biot_modulus_dP[_qp] =
         -(1 - _alpha) * _dporosity_dP[_qp] * _one_over_K + _dporosity_dP[_qp] * _one_over_Kf;
     _done_over_biot_modulus_dep[_qp] =
         -(1 - _alpha) * _dporosity_dep[_qp] * _one_over_K + _dporosity_dep[_qp] * _one_over_Kf;
   }
 }

◆ initQpStatefulProperties()

void PoroFullSatMaterial::initQpStatefulProperties ( )

protectedvirtual

Definition at line 79 of file PoroFullSatMaterial.C.

 {
   _vol_strain[_qp] = 0.0;
 }

Member Data Documentation

◆ _alpha

Real PoroFullSatMaterial::_alpha

protected

Biot coefficient.

Definition at line 37 of file PoroFullSatMaterial.h.

Referenced by computeQpProperties().

◆ _biot_coefficient

MaterialProperty<Real>& PoroFullSatMaterial::_biot_coefficient

protected

Biot coefficient.

Definition at line 64 of file PoroFullSatMaterial.h.

Referenced by computeQpProperties().

◆ _constant_porosity

bool PoroFullSatMaterial::_constant_porosity

protected

whether to use constant porosity (set _porosity = _phi0 always)

Definition at line 46 of file PoroFullSatMaterial.h.

Referenced by computeQpProperties().

◆ _done_over_biot_modulus_dep

MaterialProperty<Real>& PoroFullSatMaterial::_done_over_biot_modulus_dep

protected

d(1/M)/d(volumetric_strain)

Definition at line 82 of file PoroFullSatMaterial.h.

Referenced by computeQpProperties().

◆ _done_over_biot_modulus_dP

MaterialProperty<Real>& PoroFullSatMaterial::_done_over_biot_modulus_dP

protected

d(1/M)/d(porepressure)

Definition at line 79 of file PoroFullSatMaterial.h.

Referenced by computeQpProperties().

◆ _dporosity_dep

MaterialProperty<Real>& PoroFullSatMaterial::_dporosity_dep

protected

d(porosity)/d(volumetric_strain)

Definition at line 73 of file PoroFullSatMaterial.h.

Referenced by computeQpProperties().

◆ _dporosity_dP

MaterialProperty<Real>& PoroFullSatMaterial::_dporosity_dP

protected

d(porosity)/d(porepressure)

Definition at line 70 of file PoroFullSatMaterial.h.

Referenced by computeQpProperties().

◆ _grad_disp

std::vector<const VariableGradient *> PoroFullSatMaterial::_grad_disp

protected

grad(displacement)

Definition at line 58 of file PoroFullSatMaterial.h.

Referenced by computeQpProperties(), and PoroFullSatMaterial().

◆ _ndisp

unsigned int PoroFullSatMaterial::_ndisp

protected

number of displacement variables supplied

Definition at line 55 of file PoroFullSatMaterial.h.

Referenced by computeQpProperties(), and PoroFullSatMaterial().

◆ _one_over_biot_modulus

MaterialProperty<Real>& PoroFullSatMaterial::_one_over_biot_modulus

protected

1/M, where M is the Biot modulus

Definition at line 76 of file PoroFullSatMaterial.h.

Referenced by computeQpProperties().

◆ _one_over_K

Real PoroFullSatMaterial::_one_over_K

protected

1/K, where K is the solid bulk modulus. Usually 1/K = C_iijj, where C is the compliance matrix: strain_ij = C_ijkl stress_kl

Definition at line 40 of file PoroFullSatMaterial.h.

Referenced by computeQpProperties().

◆ _one_over_Kf

Real PoroFullSatMaterial::_one_over_Kf

protected

1/Kf, where Kf is the fluid bulk modulus.

Definition at line 43 of file PoroFullSatMaterial.h.

Referenced by computeQpProperties().

◆ _phi0

Real PoroFullSatMaterial::_phi0

protected

porosity at zero porepressure and volumetric strain

Definition at line 34 of file PoroFullSatMaterial.h.

Referenced by computeQpProperties().

◆ _porepressure

const VariableValue& PoroFullSatMaterial::_porepressure

protected

porepressure variable

Definition at line 49 of file PoroFullSatMaterial.h.

Referenced by computeQpProperties().

◆ _porepressure_name

std::string PoroFullSatMaterial::_porepressure_name

protected

name given by user to the porepressure variable

Definition at line 52 of file PoroFullSatMaterial.h.

◆ _porosity

MaterialProperty<Real>& PoroFullSatMaterial::_porosity

protected

porosity

Definition at line 67 of file PoroFullSatMaterial.h.

Referenced by computeQpProperties().

◆ _vol_strain

MaterialProperty<Real>& PoroFullSatMaterial::_vol_strain

protected

volumetric strain = strain_ii

Definition at line 61 of file PoroFullSatMaterial.h.

Referenced by computeQpProperties(), and initQpStatefulProperties().

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

richards/include/materials/PoroFullSatMaterial.h
richards/src/materials/PoroFullSatMaterial.C

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ PoroFullSatMaterial()

Member Function Documentation

◆ computeQpProperties()

◆ initQpStatefulProperties()

Member Data Documentation

◆ _alpha

◆ _biot_coefficient

◆ _constant_porosity

◆ _done_over_biot_modulus_dep

◆ _done_over_biot_modulus_dP

◆ _dporosity_dep

◆ _dporosity_dP

◆ _grad_disp

◆ _ndisp

◆ _one_over_biot_modulus

◆ _one_over_K

◆ _one_over_Kf

◆ _phi0

◆ _porepressure

◆ _porepressure_name

◆ _porosity

◆ _vol_strain