This Matrial evolves porosity so that the PorousFlow equations match the standard equations of poroelasticity theory with a constant BiotModulus. More...

#include <PorousFlowPorosityHMBiotModulus.h>

Inheritance diagram for PorousFlowPorosityHMBiotModulus:

[legend]

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	PorousFlowPorosityHMBiotModulus (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 void	computeQpProperties () override

virtual Real	atNegInfinityQp () const override
	Returns "a" at the quadpoint (porosity = a + (b - a) * exp(decay))

virtual Real	datNegInfinityQp (unsigned pvar) const override
	d(a)/d(PorousFlow variable pvar)

virtual Real	atZeroQp () const override
	Returns "b" at the quadpoint (porosity = a + (b - a) * exp(decay))

virtual Real	datZeroQp (unsigned pvar) const override
	d(a)/d(PorousFlow variable pvar)

virtual Real	decayQp () const override
	Returns "decay" at the quadpoint (porosity = a + (b - a) * exp(decay))

virtual Real	ddecayQp_dvar (unsigned pvar) const override
	d(decay)/d(PorousFlow variable pvar)

virtual RealGradient	ddecayQp_dgradvar (unsigned pvar) const override
	d(decay)/d(grad(PorousFlow variable pvar))

virtual void	initQpStatefulProperties () override

Protected Attributes
const MaterialProperty< Real > &	_porosity_old
	Old value of porosity.

const Real	_biot_modulus
	Constant biot modulus.

const Real	_fluid_bulk_modulus
	Constant fluid bulk modulus.

const MaterialProperty< Real > &	_pf_old
	Old value of effective fluid pressure.

const MaterialProperty< Real > &	_vol_strain_qp_old
	Old value of total volumetric strain.

const MaterialProperty< Real > &	_vol_strain_rate_qp
	Volumetric strain rate.

const MaterialProperty< std::vector< RealGradient > > &	_dvol_strain_rate_qp_dvar
	d(volumetric strain rate)/d(PorousFlow variable)

const bool	_mechanical
	Porosity is a function of volumetric strain.

const bool	_fluid
	Porosity is a function of effective porepressure.

const bool	_thermal
	Porosity is a function of temperature.

const bool	_chemical
	Porosity is a function of chemistry.

const VariableValue &	_phi0
	Porosity at zero strain and zero porepressure and zero temperature.

const Real	_biot
	Biot coefficient.

const Real	_exp_coeff
	Thermal expansion coefficient of the solid porous skeleton.

const Moose::Functor< Real > *	_solid_bulk
	Drained bulk modulus of the porous skeleton (as a functor; can be variable, function, etc.) Only used if _fluid is True.

const Real	_coeff
	Short-hand number (biot-1). Please note: this variable is not yet divided by solid_bulk.

const VariableValue &	_t_reference
	Reference temperature.

const VariableValue &	_p_reference
	Reference porepressure.

const unsigned	_num_c_ref
	Number of reference mineral concentrations provided by user.

std::vector< const VariableValue * >	_c_reference
	Reference mineral concentrations.

const unsigned	_num_initial_c
	Number of reference mineral concentrations provided by user.

std::vector< const VariableValue * >	_initial_c
	Reference mineral concentrations.

std::vector< Real >	_c_weights
	Weights for the mineral concentrations.

const MaterialProperty< Real > *const	_vol_strain_qp
	Strain (first const means we never want to dereference and change the value, second means we'll always be pointing to the same address after initialization (like a reference))

const MaterialProperty< std::vector< RealGradient > > *const	_dvol_strain_qp_dvar
	d(strain)/(dvar) (first const means we never want to dereference and change the value, second means we'll always be pointing to the same address after initialization (like a reference))

const MaterialProperty< Real > *const	_pf
	Effective porepressure at the quadpoints or nodes.

const MaterialProperty< std::vector< Real > > *const	_dpf_dvar
	d(effective porepressure)/(d porflow variable)

const MaterialProperty< Real > *const	_temperature
	Temperature at the quadpoints or nodes.

const MaterialProperty< std::vector< Real > > *const	_dtemperature_dvar
	d(temperature)/(d porflow variable)

const MaterialProperty< std::vector< Real > > *const	_mineral_conc_old
	Old value of mineral concentration at the quadpoints or nodes.

const MaterialProperty< std::vector< Real > > *const	_reaction_rate
	Reaction rate of mineralisation.

const MaterialProperty< std::vector< std::vector< Real > > > *const	_dreaction_rate_dvar
	d(reaction_rate_conc)/d(porflow variable)

const unsigned int	_aq_ph
	Aqueous phase number.

const MaterialProperty< std::vector< Real > > *const	_saturation
	Saturation.

const MaterialProperty< std::vector< std::vector< Real > > > *const	_dsaturation_dvar
	d(saturation)/d(PorousFlow var)

const bool	_strain_at_nearest_qp
	When calculating nodal porosity, use the strain at the nearest quadpoint to the node.

const bool	_ensure_positive
	for decayQp() > 0, porosity can be negative when using porosity = a + (b - a) * exp(decay).

const Real	_porosity_min
	Minimum allowed porosity.

const Real	_zero_modifier
	If the porosity_min floor is active the porosity derivatives are set to _zero_modifier times their unfloored values (rather than exactly zero) to hint to the Newton process that porosity is not strictly constant, which aids convergence.

GenericMaterialProperty< Real, is_ad > &	_porosity
	Computed porosity at the nodes or quadpoints.

MaterialProperty< std::vector< Real > > *const	_dporosity_dvar
	d(porosity)/d(PorousFlow variable)

MaterialProperty< std::vector< RealGradient > > *const	_dporosity_dgradvar
	d(porosity)/d(grad PorousFlow variable)

const unsigned int	_num_phases
	Number of phases.

const unsigned int	_num_components
	Number of fluid components.

const unsigned int	_num_var
	Number of PorousFlow variables.

Private Member Functions
bool	haveMaterialProperty (const std::string &prop_name)

std::vector< VariableName >	buildVariableVector (const VariableName &c1, const VariableName &c2, const VariableName &c3)

void	validateCouplingHelper (const MaterialPropertyName &base, const std::vector< VariableName > &c, const System &system, std::vector< VariableName > &missing)

bool	isNotObjectVariable (const VariableName &name)

Private Attributes
FEProblemBase &	_dmi_fe_problem

Detailed Description

This Matrial evolves porosity so that the PorousFlow equations match the standard equations of poroelasticity theory with a constant BiotModulus.

Even though a constant BiotModulus is not strictly correct, many analytical solutions of poroelasticity theory assume this, so PorousFlowPorosityHMBiotModulus is useful for comparing with those solutions. Otherwise it shouldn't generally be used because PorousFlowPorosity is physically more correct.

Definition at line 23 of file PorousFlowPorosityHMBiotModulus.h.

Constructor & Destructor Documentation

◆ PorousFlowPorosityHMBiotModulus()

PorousFlowPorosityHMBiotModulus::PorousFlowPorosityHMBiotModulus ( const InputParameters & parameters )

Definition at line 35 of file PorousFlowPorosityHMBiotModulus.C.

  : PorousFlowPorosity(parameters),
    _porosity_old(_nodal_material ? getMaterialPropertyOld<Real>("PorousFlow_porosity_nodal")
                                  : getMaterialPropertyOld<Real>("PorousFlow_porosity_qp")),
    _biot_modulus(getParam<Real>("constant_biot_modulus")),
    _fluid_bulk_modulus(getParam<Real>("constant_fluid_bulk_modulus")),
    _pf_old(_nodal_material
                ? getMaterialPropertyOld<Real>("PorousFlow_effective_fluid_pressure_nodal")
                : getMaterialPropertyOld<Real>("PorousFlow_effective_fluid_pressure_qp")),
    _vol_strain_qp_old(getMaterialPropertyOld<Real>("PorousFlow_total_volumetric_strain_qp")),
    _vol_strain_rate_qp(getMaterialProperty<Real>("PorousFlow_volumetric_strain_rate_qp")),
    _dvol_strain_rate_qp_dvar(getMaterialProperty<std::vector<RealGradient>>(
        "dPorousFlow_volumetric_strain_rate_qp_dvar"))
{
}

Member Function Documentation

◆ atNegInfinityQp()

Real PorousFlowPorosity::atNegInfinityQp ( ) const

overrideprotectedvirtualinherited

Returns "a" at the quadpoint (porosity = a + (b - a) * exp(decay))

Implements PorousFlowPorosityExponentialBase.

Definition at line 175 of file PorousFlowPorosity.C.

{
  /*
   *
   * Note the use of the OLD value of porosity here.
   * This strategy, which breaks the cyclic dependency between porosity
   * and mineral concentration, is used in
   * Kernel: PorousFlowPreDis
   * Material: PorousFlowPorosity
   * Material: PorousFlowAqueousPreDisChemistry
   * Material: PorousFlowAqueousPreDisMineral
   *
   */
  Real result = _biot;
  if (_chemical)
  {
    if (_t_step == 0 && !_app.isRestarting())
      for (unsigned i = 0; i < _num_c_ref; ++i)
        result -= _c_weights[i] * (*_initial_c[i])[_qp];
    else
      for (unsigned i = 0; i < _num_c_ref; ++i)
        result -= _c_weights[i] * ((*_mineral_conc_old)[_qp][i] + _dt * (*_porosity_old)[_qp] *
                                                                      (*_saturation)[_qp][_aq_ph] *
                                                                      (*_reaction_rate)[_qp][i]);
  }
  return result;
}

◆ atZeroQp()

Real PorousFlowPorosity::atZeroQp ( ) const

overrideprotectedvirtualinherited

Returns "b" at the quadpoint (porosity = a + (b - a) * exp(decay))

Implements PorousFlowPorosityExponentialBase.

Definition at line 216 of file PorousFlowPorosity.C.

{
  // note the [0] below: _phi0 is a constant monomial and we use [0] regardless of _nodal_material
  Real result = _phi0[0];
  if (_chemical)
  {
    if (_t_step == 0 && !_app.isRestarting())
      for (unsigned i = 0; i < _num_c_ref; ++i)
        result -= _c_weights[i] * ((*_initial_c[i])[_qp] - (*_c_reference[i])[_qp]);
    else
      for (unsigned i = 0; i < _num_c_ref; ++i)
        result -= _c_weights[i] * ((*_mineral_conc_old)[_qp][i] +
                                   _dt * (*_porosity_old)[_qp] * (*_saturation)[_qp][_aq_ph] *
                                       (*_reaction_rate)[_qp][i] -
                                   (*_c_reference[i])[_qp]);
  }
  return result;
}

◆ computeQpProperties()

void PorousFlowPorosityHMBiotModulus::computeQpProperties ( )

overrideprotectedvirtual

Reimplemented from PorousFlowPorosityExponentialBase.

Definition at line 52 of file PorousFlowPorosityHMBiotModulus.C.

{
  // Note that in the following _strain[_qp] is evaluated at q quadpoint
  // So _porosity_nodal[_qp], which should be the nodal value of porosity
  // actually uses the strain at a quadpoint.  This
  // is OK for LINEAR elements, as strain is constant over the element anyway.
 
  const unsigned qp_to_use =
      (_nodal_material && (_bnd || _strain_at_nearest_qp) ? nearestQP(_qp) : _qp);
 
  const Real denom = (1.0 + _vol_strain_rate_qp[qp_to_use] * _dt + _vol_strain_qp_old[qp_to_use]);
  const Real s_p = _porosity_old[_qp] * (1 + _vol_strain_qp_old[qp_to_use]);
  _porosity[_qp] = (s_p * std::exp(-((*_pf)[_qp] - _pf_old[_qp]) / _fluid_bulk_modulus) +
                    ((*_pf)[_qp] - _pf_old[_qp]) / _biot_modulus +
                    _biot * ((*_vol_strain_qp)[qp_to_use] - _vol_strain_qp_old[qp_to_use])) /
                   denom;
 
  (*_dporosity_dvar)[_qp].resize(_num_var);
  for (unsigned int v = 0; v < _num_var; ++v)
    (*_dporosity_dvar)[_qp][v] =
        (*_dpf_dvar)[_qp][v] *
        (-s_p * std::exp(-((*_pf)[_qp] - _pf_old[_qp]) / _fluid_bulk_modulus) /
             _fluid_bulk_modulus +
         1.0 / _biot_modulus) /
        denom;
 
  (*_dporosity_dgradvar)[_qp].resize(_num_var);
  for (unsigned int v = 0; v < _num_var; ++v)
    (*_dporosity_dgradvar)[_qp][v] =
        _biot * (*_dvol_strain_qp_dvar)[qp_to_use][v] / denom -
        _porosity[_qp] / denom * _dvol_strain_rate_qp_dvar[qp_to_use][v] * _dt;
}

◆ datNegInfinityQp()

Real PorousFlowPorosity::datNegInfinityQp ( unsigned pvar ) const

overrideprotectedvirtualinherited

d(a)/d(PorousFlow variable pvar)

Implements PorousFlowPorosityExponentialBase.

Definition at line 204 of file PorousFlowPorosity.C.

{
  Real result = 0.0;
  if (_chemical && (_t_step >= 1 || _app.isRestarting()))
    for (unsigned i = 0; i < _num_c_ref; ++i)
      result -= _c_weights[i] * _dt * (*_porosity_old)[_qp] *
                ((*_saturation)[_qp][_aq_ph] * (*_dreaction_rate_dvar)[_qp][i][pvar] +
                 (*_dsaturation_dvar)[_qp][_aq_ph][pvar] * (*_reaction_rate)[_qp][i]);
  return result;
}

◆ datZeroQp()

Real PorousFlowPorosity::datZeroQp ( unsigned pvar ) const

overrideprotectedvirtualinherited

d(a)/d(PorousFlow variable pvar)

Implements PorousFlowPorosityExponentialBase.

Definition at line 236 of file PorousFlowPorosity.C.

{
  Real result = 0.0;
  if (_chemical && (_t_step >= 1 || _app.isRestarting()))
    for (unsigned i = 0; i < _num_c_ref; ++i)
      result -= _c_weights[i] * _dt * (*_porosity_old)[_qp] *
                ((*_saturation)[_qp][_aq_ph] * (*_dreaction_rate_dvar)[_qp][i][pvar] +
                 (*_dsaturation_dvar)[_qp][_aq_ph][pvar] * (*_reaction_rate)[_qp][i]);
  return result;
}

◆ ddecayQp_dgradvar()

RealGradient PorousFlowPorosity::ddecayQp_dgradvar ( unsigned pvar ) const

overrideprotectedvirtualinherited

d(decay)/d(grad(PorousFlow variable pvar))

Implements PorousFlowPorosityExponentialBase.

Definition at line 337 of file PorousFlowPorosity.C.

{
  RealGradient result(0.0, 0.0, 0.0);
  if (_mechanical)
  {
    const unsigned qp_to_use =
        (_nodal_material && (_bnd || _strain_at_nearest_qp) ? nearestQP(_qp) : _qp);
    result += -(*_dvol_strain_qp_dvar)[qp_to_use][pvar];
  }
  return result;
}

◆ ddecayQp_dvar()

Real PorousFlowPorosity::ddecayQp_dvar ( unsigned pvar ) const

overrideprotectedvirtualinherited

d(decay)/d(PorousFlow variable pvar)

Implements PorousFlowPorosityExponentialBase.

Definition at line 298 of file PorousFlowPorosity.C.

{
  Real result = 0.0;
 
  if (_thermal)
    result += _exp_coeff * (*_dtemperature_dvar)[_qp][pvar];
 
  if (_fluid)
  {
    Real solid_bulk;
    // Using Qp 0 can leverage the functor caching
    // TODO: Find a way to effectively use subdomain-constant-ness
    unsigned int qp_used = (_constant_option == ConstantTypeEnum::NONE) ? _qp : 0;
    if (_nodal_material)
    {
      const std::set<SubdomainID> subdomain_set = {_current_elem->subdomain_id()};
      const Moose::NodeArg space_arg = {_current_elem->node_ptr(qp_used), &subdomain_set};
      solid_bulk = (*_solid_bulk)(space_arg, Moose::currentState());
    }
    else if (_bnd)
    {
      const Moose::ElemSideQpArg space_arg = {
          _current_elem, _current_side, qp_used, _qrule, _q_point[qp_used]};
      solid_bulk = (*_solid_bulk)(space_arg, Moose::currentState());
    }
    else
    {
      const Moose::ElemQpArg space_arg = {_current_elem, qp_used, _qrule, _q_point[qp_used]};
      solid_bulk = (*_solid_bulk)(space_arg, Moose::currentState());
    }
    if (solid_bulk <= 0)
      mooseError("PorousFlowPorosity: solid_bulk must be larger than Zero.");
    result += _coeff / solid_bulk * (*_dpf_dvar)[_qp][pvar];
  }
 
  return result;
}

◆ decayQp()

Real PorousFlowPorosity::decayQp ( ) const

overrideprotectedvirtualinherited

Returns "decay" at the quadpoint (porosity = a + (b - a) * exp(decay))

Implements PorousFlowPorosityExponentialBase.

Definition at line 248 of file PorousFlowPorosity.C.

{
  Real result = 0.0;
 
  if (_thermal)
    result += _exp_coeff * ((*_temperature)[_qp] - _t_reference[_qp]);
 
  if (_fluid)
  {
    Real solid_bulk;
    // Using Qp 0 can leverage the functor caching
    // TODO: Find a way to effectively use subdomain-constant-ness
    unsigned int qp_used = (_constant_option == ConstantTypeEnum::NONE) ? _qp : 0;
    if (_nodal_material)
    {
      const std::set<SubdomainID> subdomain_set = {_current_elem->subdomain_id()};
      const Moose::NodeArg space_arg = {_current_elem->node_ptr(qp_used), &subdomain_set};
      solid_bulk = (*_solid_bulk)(space_arg, Moose::currentState());
    }
    else if (_bnd)
    {
      const Moose::ElemSideQpArg space_arg = {
          _current_elem, _current_side, qp_used, _qrule, _q_point[qp_used]};
      solid_bulk = (*_solid_bulk)(space_arg, Moose::currentState());
    }
    else
    {
      const Moose::ElemQpArg space_arg = {_current_elem, qp_used, _qrule, _q_point[qp_used]};
      solid_bulk = (*_solid_bulk)(space_arg, Moose::currentState());
    }
    if (solid_bulk <= 0)
      mooseError("PorousFlowPorosity: solid_bulk must be larger than Zero");
    result += _coeff / solid_bulk * ((*_pf)[_qp] - _p_reference[_qp]);
  }
 
  if (_mechanical)
  {
    // Note that in the following _strain[_qp] is evaluated at q quadpoint
    // So _porosity_nodal[_qp], which should be the nodal value of porosity
    // actually uses the strain at a quadpoint.  This
    // is OK for LINEAR elements, as strain is constant over the element anyway.
    const unsigned qp_to_use =
        (_nodal_material && (_bnd || _strain_at_nearest_qp) ? nearestQP(_qp) : _qp);
    result += -(*_vol_strain_qp)[qp_to_use];
  }
 
  return result;
}

◆ initQpStatefulProperties()

void PorousFlowPorosityExponentialBase::initQpStatefulProperties ( )

overrideprotectedvirtualinherited

Definition at line 64 of file PorousFlowPorosityExponentialBase.C.

{
  const Real a = atNegInfinityQp();
  const Real b = atZeroQp();
  mooseAssert(a > b, "PorousFlowPorosityExponentialBase a must be larger than b");
  const Real decay = decayQp();
 
  if (decay <= 0.0 || !_ensure_positive)
    _porosity[_qp] = a + (b - a) * std::exp(decay);
  else
  {
    const Real c = std::log(a / (a - b));
    const Real expx = std::exp(-decay / c);
    _porosity[_qp] = a + (b - a) * std::exp(c * (1.0 - expx));
  }
 
  if (_porosity[_qp] < _porosity_min)
    _porosity[_qp] = _porosity_min;
}

◆ validParams()

InputParameters PorousFlowPorosityHMBiotModulus::validParams ( )

static

Definition at line 16 of file PorousFlowPorosityHMBiotModulus.C.

{
  InputParameters params = PorousFlowPorosity::validParams();
  params.set<bool>("mechanical") = true;
  params.set<bool>("fluid") = true;
  params.addRequiredRangeCheckedParam<Real>("constant_biot_modulus",
                                            "constant_biot_modulus>0",
                                            "Biot modulus, which is constant for this Material");
  params.addRequiredRangeCheckedParam<Real>(
      "constant_fluid_bulk_modulus",
      "constant_fluid_bulk_modulus>0",
      "Fluid bulk modulus, which is constant for this Material");
  params.addClassDescription(
      "This Material calculates the porosity for hydro-mechanical simulations, assuming that the "
      "Biot modulus and the fluid bulk modulus are both constant.  This is useful for comparing "
      "with solutions from poroelasticity theory, but is less accurate than PorousFlowPorosity");
  return params;
}

Member Data Documentation

◆ _aq_ph

const unsigned int PorousFlowPorosity::_aq_ph

protectedinherited

Aqueous phase number.

Definition at line 116 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosity::atNegInfinityQp(), PorousFlowPorosity::atZeroQp(), PorousFlowPorosity::datNegInfinityQp(), and PorousFlowPorosity::datZeroQp().

◆ _biot

const Real PorousFlowPorosity::_biot

protectedinherited

Biot coefficient.

Definition at line 52 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosity::atNegInfinityQp(), and computeQpProperties().

◆ _biot_modulus

const Real PorousFlowPorosityHMBiotModulus::_biot_modulus

protected

Constant biot modulus.

Definition at line 37 of file PorousFlowPorosityHMBiotModulus.h.

Referenced by computeQpProperties().

◆ _c_reference

std::vector<const VariableValue *> PorousFlowPorosity::_c_reference

protectedinherited

Reference mineral concentrations.

Definition at line 74 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosity::atZeroQp(), and PorousFlowPorosity::PorousFlowPorosity().

◆ _c_weights

std::vector<Real> PorousFlowPorosity::_c_weights

protectedinherited

Weights for the mineral concentrations.

Definition at line 83 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosity::atNegInfinityQp(), PorousFlowPorosity::atZeroQp(), PorousFlowPorosity::datNegInfinityQp(), PorousFlowPorosity::datZeroQp(), and PorousFlowPorosity::PorousFlowPorosity().

◆ _chemical

const bool PorousFlowPorosity::_chemical

protectedinherited

Porosity is a function of chemistry.

Definition at line 46 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosity::atNegInfinityQp(), PorousFlowPorosity::atZeroQp(), PorousFlowPorosity::datNegInfinityQp(), PorousFlowPorosity::datZeroQp(), and PorousFlowPorosity::PorousFlowPorosity().

◆ _coeff

const Real PorousFlowPorosity::_coeff

protectedinherited

Short-hand number (biot-1). Please note: this variable is not yet divided by solid_bulk.

Definition at line 62 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosity::ddecayQp_dvar(), and PorousFlowPorosity::decayQp().

◆ _dpf_dvar

const MaterialProperty<std::vector<Real> >* const PorousFlowPorosity::_dpf_dvar

protectedinherited

d(effective porepressure)/(d porflow variable)

Definition at line 98 of file PorousFlowPorosity.h.

◆ _dporosity_dgradvar

template<bool is_ad>

MaterialProperty<std::vector<RealGradient> >* const PorousFlowPorosityBaseTempl< is_ad >::_dporosity_dgradvar

protectedinherited

d(porosity)/d(grad PorousFlow variable)

Definition at line 33 of file PorousFlowPorosityBase.h.

Referenced by computeQpProperties(), and PorousFlowPorosityLinear::computeQpProperties().

◆ _dporosity_dvar

template<bool is_ad>

MaterialProperty<std::vector<Real> >* const PorousFlowPorosityBaseTempl< is_ad >::_dporosity_dvar

protectedinherited

d(porosity)/d(PorousFlow variable)

Definition at line 30 of file PorousFlowPorosityBase.h.

Referenced by computeQpProperties(), and PorousFlowPorosityLinear::computeQpProperties().

◆ _dreaction_rate_dvar

const MaterialProperty<std::vector<std::vector<Real> > >* const PorousFlowPorosity::_dreaction_rate_dvar

protectedinherited

d(reaction_rate_conc)/d(porflow variable)

Definition at line 113 of file PorousFlowPorosity.h.

◆ _dsaturation_dvar

const MaterialProperty<std::vector<std::vector<Real> > >* const PorousFlowPorosity::_dsaturation_dvar

protectedinherited

d(saturation)/d(PorousFlow var)

Definition at line 122 of file PorousFlowPorosity.h.

◆ _dtemperature_dvar

const MaterialProperty<std::vector<Real> >* const PorousFlowPorosity::_dtemperature_dvar

protectedinherited

d(temperature)/(d porflow variable)

Definition at line 104 of file PorousFlowPorosity.h.

◆ _dvol_strain_qp_dvar

const MaterialProperty<std::vector<RealGradient> >* const PorousFlowPorosity::_dvol_strain_qp_dvar

protectedinherited

d(strain)/(dvar) (first const means we never want to dereference and change the value, second means we'll always be pointing to the same address after initialization (like a reference))

Definition at line 92 of file PorousFlowPorosity.h.

◆ _dvol_strain_rate_qp_dvar

const MaterialProperty<std::vector<RealGradient> >& PorousFlowPorosityHMBiotModulus::_dvol_strain_rate_qp_dvar

protected

d(volumetric strain rate)/d(PorousFlow variable)

Definition at line 52 of file PorousFlowPorosityHMBiotModulus.h.

Referenced by computeQpProperties().

◆ _ensure_positive

const bool PorousFlowPorosityExponentialBase::_ensure_positive

protectedinherited

for decayQp() > 0, porosity can be negative when using porosity = a + (b - a) * exp(decay).

This expression is modified if ensure_positive = true to read porosity = a + (b - a) * exp(c * (1 - Exp(- decay / c))) where c = log(a/(b-a))

Definition at line 73 of file PorousFlowPorosityExponentialBase.h.

Referenced by PorousFlowPorosityExponentialBase::computeQpProperties(), and PorousFlowPorosityExponentialBase::initQpStatefulProperties().

◆ _exp_coeff

const Real PorousFlowPorosity::_exp_coeff

protectedinherited

Thermal expansion coefficient of the solid porous skeleton.

Definition at line 55 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosity::ddecayQp_dvar(), and PorousFlowPorosity::decayQp().

◆ _fluid

const bool PorousFlowPorosity::_fluid

protectedinherited

Porosity is a function of effective porepressure.

Definition at line 40 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosity::ddecayQp_dvar(), PorousFlowPorosity::decayQp(), and PorousFlowPorosity::PorousFlowPorosity().

◆ _fluid_bulk_modulus

const Real PorousFlowPorosityHMBiotModulus::_fluid_bulk_modulus

protected

Constant fluid bulk modulus.

Definition at line 40 of file PorousFlowPorosityHMBiotModulus.h.

Referenced by computeQpProperties().

◆ _initial_c

std::vector<const VariableValue *> PorousFlowPorosity::_initial_c

protectedinherited

Reference mineral concentrations.

Definition at line 80 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosity::atNegInfinityQp(), PorousFlowPorosity::atZeroQp(), and PorousFlowPorosity::PorousFlowPorosity().

◆ _mechanical

const bool PorousFlowPorosity::_mechanical

protectedinherited

Porosity is a function of volumetric strain.

Definition at line 37 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosity::ddecayQp_dgradvar(), and PorousFlowPorosity::decayQp().

◆ _mineral_conc_old

const MaterialProperty<std::vector<Real> >* const PorousFlowPorosity::_mineral_conc_old

protectedinherited

Old value of mineral concentration at the quadpoints or nodes.

Definition at line 107 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosity::atNegInfinityQp(), and PorousFlowPorosity::atZeroQp().

◆ _num_c_ref

const unsigned PorousFlowPorosity::_num_c_ref

protectedinherited

Number of reference mineral concentrations provided by user.

Definition at line 71 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosity::atNegInfinityQp(), PorousFlowPorosity::atZeroQp(), PorousFlowPorosity::datNegInfinityQp(), PorousFlowPorosity::datZeroQp(), and PorousFlowPorosity::PorousFlowPorosity().

◆ _num_components

const unsigned int PorousFlowMaterialVectorBase::_num_components

protectedinherited

Number of fluid components.

Definition at line 30 of file PorousFlowMaterialVectorBase.h.

Referenced by PorousFlowMassFractionAqueousEquilibriumChemistry::computeQpProperties(), PorousFlowAqueousPreDisChemistry::PorousFlowAqueousPreDisChemistry(), PorousFlowDiffusivityBaseTempl< is_ad >::PorousFlowDiffusivityBaseTempl(), and PorousFlowMassFractionTempl< is_ad >::PorousFlowMassFractionTempl().

◆ _num_initial_c

const unsigned PorousFlowPorosity::_num_initial_c

protectedinherited

Number of reference mineral concentrations provided by user.

Definition at line 77 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosity::PorousFlowPorosity().

◆ _num_phases

const unsigned int PorousFlowMaterialVectorBase::_num_phases

protectedinherited

Number of phases.

Definition at line 27 of file PorousFlowMaterialVectorBase.h.

Referenced by PorousFlowDiffusivityBaseTempl< is_ad >::PorousFlowDiffusivityBaseTempl(), PorousFlowDiffusivityConstTempl< is_ad >::PorousFlowDiffusivityConstTempl(), PorousFlowJoinerTempl< is_ad >::PorousFlowJoinerTempl(), PorousFlowMassFractionTempl< is_ad >::PorousFlowMassFractionTempl(), PorousFlowThermalConductivityFromPorosityTempl< is_ad >::PorousFlowThermalConductivityFromPorosityTempl(), PorousFlowThermalConductivityIdealTempl< is_ad >::PorousFlowThermalConductivityIdealTempl(), and PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity::PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity().

◆ _num_var

const unsigned int PorousFlowMaterialVectorBase::_num_var

protectedinherited

Number of PorousFlow variables.

Definition at line 33 of file PorousFlowMaterialVectorBase.h.

Referenced by PorousFlowAqueousPreDisChemistry::computeQpProperties(), PorousFlowMassFractionAqueousEquilibriumChemistry::computeQpProperties(), PorousFlowPorosityExponentialBase::computeQpProperties(), computeQpProperties(), PorousFlowPorosityLinear::computeQpProperties(), PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity::computeQpProperties(), PorousFlowVolumetricStrain::computeQpProperties(), PorousFlowAqueousPreDisChemistry::initQpStatefulProperties(), and PorousFlowPorosityLinear::initQpStatefulProperties().

◆ _p_reference

const VariableValue& PorousFlowPorosity::_p_reference

protectedinherited

Reference porepressure.

Definition at line 68 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosity::decayQp().

◆ _pf

const MaterialProperty<Real>* const PorousFlowPorosity::_pf

protectedinherited

Effective porepressure at the quadpoints or nodes.

Definition at line 95 of file PorousFlowPorosity.h.

Referenced by computeQpProperties().

◆ _pf_old

const MaterialProperty<Real>& PorousFlowPorosityHMBiotModulus::_pf_old

protected

Old value of effective fluid pressure.

Definition at line 43 of file PorousFlowPorosityHMBiotModulus.h.

Referenced by computeQpProperties().

◆ _phi0

const VariableValue& PorousFlowPorosity::_phi0

protectedinherited

Porosity at zero strain and zero porepressure and zero temperature.

Definition at line 49 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosity::atZeroQp().

◆ _porosity

template<bool is_ad>

GenericMaterialProperty<Real, is_ad>& PorousFlowPorosityBaseTempl< is_ad >::_porosity

protectedinherited

Computed porosity at the nodes or quadpoints.

Definition at line 27 of file PorousFlowPorosityBase.h.

Referenced by PorousFlowPorosityExponentialBase::computeQpProperties(), computeQpProperties(), PorousFlowPorosityLinear::computeQpProperties(), PorousFlowPorosityExponentialBase::initQpStatefulProperties(), and PorousFlowPorosityLinear::initQpStatefulProperties().

◆ _porosity_min

const Real PorousFlowPorosityExponentialBase::_porosity_min

protectedinherited

Minimum allowed porosity.

The ensure_positive transform above only acts for decay > 0, so chemistry-driven porosity (which enters through "a" and "b", not "decay") is otherwise unbounded below and can go negative once pore space is filled by precipitated mineral. If the computed porosity is less than this value it is set to this value instead. Defaults to no floor.

Definition at line 82 of file PorousFlowPorosityExponentialBase.h.

Referenced by PorousFlowPorosityExponentialBase::computeQpProperties(), and PorousFlowPorosityExponentialBase::initQpStatefulProperties().

◆ _porosity_old

const MaterialProperty<Real>& PorousFlowPorosityHMBiotModulus::_porosity_old

protected

Old value of porosity.

Definition at line 34 of file PorousFlowPorosityHMBiotModulus.h.

Referenced by computeQpProperties().

◆ _reaction_rate

const MaterialProperty<std::vector<Real> >* const PorousFlowPorosity::_reaction_rate

protectedinherited

Reaction rate of mineralisation.

Definition at line 110 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosity::atNegInfinityQp(), and PorousFlowPorosity::atZeroQp().

◆ _saturation

const MaterialProperty<std::vector<Real> >* const PorousFlowPorosity::_saturation

protectedinherited

Saturation.

Definition at line 119 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosity::atNegInfinityQp(), and PorousFlowPorosity::atZeroQp().

◆ _solid_bulk

const Moose::Functor<Real>* PorousFlowPorosity::_solid_bulk

protectedinherited

Drained bulk modulus of the porous skeleton (as a functor; can be variable, function, etc.) Only used if _fluid is True.

Definition at line 59 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosity::PorousFlowPorosity().

◆ _strain_at_nearest_qp

const bool PorousFlowPorosityExponentialBase::_strain_at_nearest_qp

protectedinherited

When calculating nodal porosity, use the strain at the nearest quadpoint to the node.

Definition at line 64 of file PorousFlowPorosityExponentialBase.h.

Referenced by computeQpProperties(), PorousFlowPorosity::ddecayQp_dgradvar(), and PorousFlowPorosity::decayQp().

◆ _t_reference

const VariableValue& PorousFlowPorosity::_t_reference

protectedinherited

Reference temperature.

Definition at line 65 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosity::decayQp().

◆ _temperature

const MaterialProperty<Real>* const PorousFlowPorosity::_temperature

protectedinherited

Temperature at the quadpoints or nodes.

Definition at line 101 of file PorousFlowPorosity.h.

◆ _thermal

const bool PorousFlowPorosity::_thermal

protectedinherited

Porosity is a function of temperature.

Definition at line 43 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosity::ddecayQp_dvar(), PorousFlowPorosity::decayQp(), and PorousFlowPorosity::PorousFlowPorosity().

◆ _vol_strain_qp

const MaterialProperty<Real>* const PorousFlowPorosity::_vol_strain_qp

protectedinherited

Strain (first const means we never want to dereference and change the value, second means we'll always be pointing to the same address after initialization (like a reference))

Definition at line 89 of file PorousFlowPorosity.h.

Referenced by computeQpProperties().

◆ _vol_strain_qp_old

const MaterialProperty<Real>& PorousFlowPorosityHMBiotModulus::_vol_strain_qp_old

protected

Old value of total volumetric strain.

Definition at line 46 of file PorousFlowPorosityHMBiotModulus.h.

Referenced by computeQpProperties().

◆ _vol_strain_rate_qp

const MaterialProperty<Real>& PorousFlowPorosityHMBiotModulus::_vol_strain_rate_qp

protected

Volumetric strain rate.

Definition at line 49 of file PorousFlowPorosityHMBiotModulus.h.

Referenced by computeQpProperties().

◆ _zero_modifier

const Real PorousFlowPorosityExponentialBase::_zero_modifier

protectedinherited

If the porosity_min floor is active the porosity derivatives are set to _zero_modifier times their unfloored values (rather than exactly zero) to hint to the Newton process that porosity is not strictly constant, which aids convergence.

Definition at line 90 of file PorousFlowPorosityExponentialBase.h.

Referenced by PorousFlowPorosityExponentialBase::computeQpProperties().

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

porous_flow/include/materials/PorousFlowPorosityHMBiotModulus.h
porous_flow/src/materials/PorousFlowPorosityHMBiotModulus.C

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ PorousFlowPorosityHMBiotModulus()

Member Function Documentation

◆ atNegInfinityQp()

◆ atZeroQp()

◆ computeQpProperties()

◆ datNegInfinityQp()

◆ datZeroQp()

◆ ddecayQp_dgradvar()

◆ ddecayQp_dvar()

◆ decayQp()

◆ initQpStatefulProperties()

◆ validParams()

Member Data Documentation

◆ _aq_ph

◆ _biot

◆ _biot_modulus

◆ _c_reference

◆ _c_weights

◆ _chemical

◆ _coeff

◆ _dpf_dvar

◆ _dporosity_dgradvar

◆ _dporosity_dvar

◆ _dreaction_rate_dvar

◆ _dsaturation_dvar

◆ _dtemperature_dvar

◆ _dvol_strain_qp_dvar

◆ _dvol_strain_rate_qp_dvar

◆ _ensure_positive

◆ _exp_coeff

◆ _fluid

◆ _fluid_bulk_modulus

◆ _initial_c

◆ _mechanical

◆ _mineral_conc_old

◆ _num_c_ref

◆ _num_components

◆ _num_initial_c

◆ _num_phases

◆ _num_var

◆ _p_reference

◆ _pf

◆ _pf_old

◆ _phi0

◆ _porosity

◆ _porosity_min

◆ _porosity_old

◆ _reaction_rate

◆ _saturation

◆ _solid_bulk

◆ _strain_at_nearest_qp

◆ _t_reference

◆ _temperature

◆ _thermal

◆ _vol_strain_qp

◆ _vol_strain_qp_old

◆ _vol_strain_rate_qp

◆ _zero_modifier