PorousFlowVolumetricStrain Class Reference

PorousFlowVolumetricStrain computes volumetric strains, and derivatives thereof. More...

#include <PorousFlowVolumetricStrain.h>

Inheritance diagram for PorousFlowVolumetricStrain:

Public Member Functions
	PorousFlowVolumetricStrain (const InputParameters &parameters)

Protected Member Functions
virtual void	initQpStatefulProperties () override
virtual void	computeQpProperties () override

Protected Attributes
const bool	_consistent
	If true then the strain rate will include terms that ensure mass is conserved when doing integrals over the displaced mesh. More...
const unsigned int	_ndisp
	Number of displacements supplied (1 in 1D, 2 in 2D, 3 in 3D) More...
std::vector< const VariableValue * >	_disp
	Displacement variable values at the quad point. More...
std::vector< unsigned int >	_disp_var_num
	MOOSE variable number of the displacements variables provided. More...
std::vector< const VariableGradient * >	_grad_disp
	Gradient of the displacements. More...
std::vector< const VariableGradient * >	_grad_disp_old
	Old value of gradient of the displacements. More...
MaterialProperty< Real > &	_vol_strain_rate_qp
	The volumetric strain rate at the quadpoints. More...
MaterialProperty< std::vector< RealGradient > > &	_dvol_strain_rate_qp_dvar
	The derivative of the volumetric strain rate with respect to the porous flow variables. More...
MaterialProperty< Real > &	_vol_total_strain_qp
	The total volumetric strain at the quadpoints. More...
MaterialProperty< std::vector< RealGradient > > &	_dvol_total_strain_qp_dvar
	The derivative of the total volumetric strain with respect to the porous flow variables. More...
const unsigned int	_num_phases
	Number of phases. More...
const unsigned int	_num_components
	Number of fluid components. More...
const unsigned int	_num_var
	Number of PorousFlow variables. More...

Detailed Description

PorousFlowVolumetricStrain computes volumetric strains, and derivatives thereof.

Definition at line 24 of file PorousFlowVolumetricStrain.h.

Constructor & Destructor Documentation

PorousFlowVolumetricStrain()

PorousFlowVolumetricStrain::PorousFlowVolumetricStrain

(

const InputParameters &

parameters

)

Definition at line 39 of file PorousFlowVolumetricStrain.C.

  : PorousFlowMaterialVectorBase(parameters),
    _consistent(getParam<bool>("consistent_with_displaced_mesh")),
    _ndisp(coupledComponents("displacements")),
    _disp(3),
    _disp_var_num(3),
    _grad_disp(3),
    _grad_disp_old(3),

    _vol_strain_rate_qp(declareProperty<Real>("PorousFlow_volumetric_strain_rate_qp")),
    _dvol_strain_rate_qp_dvar(
        declareProperty<std::vector<RealGradient>>("dPorousFlow_volumetric_strain_rate_qp_dvar")),
    _vol_total_strain_qp(declareProperty<Real>("PorousFlow_total_volumetric_strain_qp")),
    _dvol_total_strain_qp_dvar(
        declareProperty<std::vector<RealGradient>>("dPorousFlow_total_volumetric_strain_qp_dvar"))
{
  if (_ndisp != _mesh.dimension())
    paramError("displacements", "The number of variables supplied must match the mesh dimension.");

  // fetch coupled variables and gradients (as stateful properties if necessary)
  for (unsigned int i = 0; i < _ndisp; ++i)
  {
    _disp[i] = &coupledValue("displacements", i);
    _disp_var_num[i] = coupled("displacements", i);
    _grad_disp[i] = &coupledGradient("displacements", i);
    _grad_disp_old[i] = &coupledGradientOld("displacements", i);
  }

  // set unused dimensions to zero
  for (unsigned i = _ndisp; i < 3; ++i)
  {
    _disp[i] = &_zero;
    _disp_var_num[i] = 0;
    while (_dictator.isPorousFlowVariable(_disp_var_num[i]))
      _disp_var_num[i]++; // increment until disp_var_num[i] is not a porflow var
    _grad_disp[i] = &_grad_zero;
    _grad_disp_old[i] = &_grad_zero;
  }
  if (_nodal_material == true)
    mooseError("PorousFlowVolumetricStrain classes are only defined for at_nodes = false");
}

Member Function Documentation

computeQpProperties()

void PorousFlowVolumetricStrain::computeQpProperties ( )

overrideprotectedvirtual

Definition at line 88 of file PorousFlowVolumetricStrain.C.

{
  RankTwoTensor A(
      (*_grad_disp[0])[_qp], (*_grad_disp[1])[_qp], (*_grad_disp[2])[_qp]); // Deformation gradient
  RankTwoTensor Fbar((*_grad_disp_old[0])[_qp],
                     (*_grad_disp_old[1])[_qp],
                     (*_grad_disp_old[2])[_qp]); // Old Deformation gradient

  _vol_total_strain_qp[_qp] = A.trace();

  A -= Fbar; // A = grad_disp - grad_disp_old

  RankTwoTensor total_strain_increment = 0.5 * (A + A.transpose());
  const Real andy = (_consistent ? 1.0 + (*_grad_disp_old[0])[_qp](0) +
                                       (*_grad_disp_old[1])[_qp](1) + (*_grad_disp_old[2])[_qp](2)
                                 : 1.0);
  _vol_strain_rate_qp[_qp] = total_strain_increment.trace() / _dt / andy;

  // prepare the derivatives with zeroes
  _dvol_strain_rate_qp_dvar[_qp].resize(_num_var, RealGradient());
  _dvol_total_strain_qp_dvar[_qp].resize(_num_var, RealGradient());
  for (unsigned i = 0; i < _ndisp; ++i)
    if (_dictator.isPorousFlowVariable(_disp_var_num[i]))
    {
      // the i_th displacement is a PorousFlow variable
      const unsigned int pvar = _dictator.porousFlowVariableNum(_disp_var_num[i]);
      _dvol_strain_rate_qp_dvar[_qp][pvar](i) = 1.0 / _dt / andy;
      _dvol_total_strain_qp_dvar[_qp][pvar](i) = 1.0;
    }
}

initQpStatefulProperties()

void PorousFlowVolumetricStrain::initQpStatefulProperties ( )

overrideprotectedvirtual

Definition at line 82 of file PorousFlowVolumetricStrain.C.

{
  _vol_total_strain_qp[_qp] = 0.0;
}

Member Data Documentation

_consistent

const bool PorousFlowVolumetricStrain::_consistent

protected

If true then the strain rate will include terms that ensure mass is conserved when doing integrals over the displaced mesh.

Definition at line 34 of file PorousFlowVolumetricStrain.h.

Referenced by computeQpProperties().

_disp

std::vector<const VariableValue *> PorousFlowVolumetricStrain::_disp

protected

Displacement variable values at the quad point.

Definition at line 40 of file PorousFlowVolumetricStrain.h.

Referenced by PorousFlowVolumetricStrain().

_disp_var_num

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

protected

MOOSE variable number of the displacements variables provided.

Definition at line 43 of file PorousFlowVolumetricStrain.h.

Referenced by computeQpProperties(), and PorousFlowVolumetricStrain().

_dvol_strain_rate_qp_dvar

MaterialProperty<std::vector<RealGradient> >& PorousFlowVolumetricStrain::_dvol_strain_rate_qp_dvar

protected

The derivative of the volumetric strain rate with respect to the porous flow variables.

Since the volumetric strain rate depends on derivatives of the displacement variables, this should be multiplied by _grad_phi in kernels

Definition at line 59 of file PorousFlowVolumetricStrain.h.

Referenced by computeQpProperties().

_dvol_total_strain_qp_dvar

MaterialProperty<std::vector<RealGradient> >& PorousFlowVolumetricStrain::_dvol_total_strain_qp_dvar

protected

The derivative of the total volumetric strain with respect to the porous flow variables.

Since the total volumetric strain depends on derivatives of the displacement variables, this should be multiplied by _grad_phi in kernels

Definition at line 69 of file PorousFlowVolumetricStrain.h.

Referenced by computeQpProperties().

_grad_disp

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

protected

Gradient of the displacements.

Definition at line 46 of file PorousFlowVolumetricStrain.h.

Referenced by computeQpProperties(), and PorousFlowVolumetricStrain().

_grad_disp_old

std::vector<const VariableGradient *> PorousFlowVolumetricStrain::_grad_disp_old

protected

Old value of gradient of the displacements.

Definition at line 49 of file PorousFlowVolumetricStrain.h.

Referenced by computeQpProperties(), and PorousFlowVolumetricStrain().

_ndisp

const unsigned int PorousFlowVolumetricStrain::_ndisp

protected

Number of displacements supplied (1 in 1D, 2 in 2D, 3 in 3D)

Definition at line 37 of file PorousFlowVolumetricStrain.h.

Referenced by computeQpProperties(), and PorousFlowVolumetricStrain().

_num_components

const unsigned int PorousFlowMaterialVectorBase::_num_components

protectedinherited

Number of fluid components.

Definition at line 34 of file PorousFlowMaterialVectorBase.h.

Referenced by PorousFlowDiffusivityBase::computeQpProperties(), PorousFlowMassFractionAqueousEquilibriumChemistry::computeQpProperties(), PorousFlowMassFraction::computeQpProperties(), PorousFlowAqueousPreDisChemistry::PorousFlowAqueousPreDisChemistry(), PorousFlowDiffusivityBase::PorousFlowDiffusivityBase(), and PorousFlowMassFraction::PorousFlowMassFraction().

_num_phases

const unsigned int PorousFlowMaterialVectorBase::_num_phases

protectedinherited

Number of phases.

Definition at line 31 of file PorousFlowMaterialVectorBase.h.

Referenced by PorousFlowDiffusivityBase::computeQpProperties(), PorousFlowDiffusivityMillingtonQuirk::computeQpProperties(), PorousFlowEffectiveFluidPressure::computeQpProperties(), PorousFlowMassFraction::computeQpProperties(), PorousFlowJoiner::computeQpProperties(), PorousFlowEffectiveFluidPressure::initQpStatefulProperties(), PorousFlowJoiner::initQpStatefulProperties(), PorousFlowDiffusivityBase::PorousFlowDiffusivityBase(), PorousFlowDiffusivityConst::PorousFlowDiffusivityConst(), PorousFlowJoiner::PorousFlowJoiner(), PorousFlowMassFraction::PorousFlowMassFraction(), PorousFlowThermalConductivityFromPorosity::PorousFlowThermalConductivityFromPorosity(), PorousFlowThermalConductivityIdeal::PorousFlowThermalConductivityIdeal(), and PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity::PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity().

_num_var

const unsigned int PorousFlowMaterialVectorBase::_num_var

protectedinherited

Number of PorousFlow variables.

Definition at line 37 of file PorousFlowMaterialVectorBase.h.

Referenced by PorousFlowDiffusivityBase::computeQpProperties(), PorousFlowPermeabilityConst::computeQpProperties(), computeQpProperties(), PorousFlowDiffusivityMillingtonQuirk::computeQpProperties(), PorousFlowAqueousPreDisChemistry::computeQpProperties(), PorousFlowPermeabilityTensorFromVar::computeQpProperties(), PorousFlowPorosityConst::computeQpProperties(), PorousFlowMassFractionAqueousEquilibriumChemistry::computeQpProperties(), PorousFlowThermalConductivityIdeal::computeQpProperties(), PorousFlowEffectiveFluidPressure::computeQpProperties(), PorousFlowThermalConductivityFromPorosity::computeQpProperties(), PorousFlowMatrixInternalEnergy::computeQpProperties(), PorousFlowPorosityHMBiotModulus::computeQpProperties(), PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity::computeQpProperties(), PorousFlowPermeabilityConstFromVar::computeQpProperties(), PorousFlowPermeabilityExponential::computeQpProperties(), PorousFlowMassFraction::computeQpProperties(), PorousFlowPermeabilityKozenyCarman::computeQpProperties(), PorousFlowJoiner::computeQpProperties(), PorousFlowPorosityExponentialBase::computeQpProperties(), and PorousFlowAqueousPreDisChemistry::initQpStatefulProperties().

_vol_strain_rate_qp

MaterialProperty<Real>& PorousFlowVolumetricStrain::_vol_strain_rate_qp

protected

The volumetric strain rate at the quadpoints.

Definition at line 52 of file PorousFlowVolumetricStrain.h.

Referenced by computeQpProperties().

_vol_total_strain_qp

MaterialProperty<Real>& PorousFlowVolumetricStrain::_vol_total_strain_qp

protected

The total volumetric strain at the quadpoints.

Definition at line 62 of file PorousFlowVolumetricStrain.h.

Referenced by computeQpProperties(), and initQpStatefulProperties().

---

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

• porous_flow/include/materials/PorousFlowVolumetricStrain.h
• porous_flow/src/materials/PorousFlowVolumetricStrain.C