www.mooseframework.org
Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
PorousFlowPorosity Class Reference

Material designed to provide the porosity in PorousFlow simulations chemistry + biot + (phi0 - reference_chemistry - biot) * exp(-vol_strain. More...

#include <PorousFlowPorosity.h>

Inheritance diagram for PorousFlowPorosity:
[legend]

Public Types

typedef DerivativeMaterialPropertyNameInterface::SymbolName SymbolName
 

Public Member Functions

 PorousFlowPorosity (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 atNegInfinityQp () const override
 Returns "a" at the quadpoint (porosity = a + (b - a) * exp(decay)) More...
 
virtual Real datNegInfinityQp (unsigned pvar) const override
 d(a)/d(PorousFlow variable pvar) More...
 
virtual Real atZeroQp () const override
 Returns "b" at the quadpoint (porosity = a + (b - a) * exp(decay)) More...
 
virtual Real datZeroQp (unsigned pvar) const override
 d(a)/d(PorousFlow variable pvar) More...
 
virtual Real decayQp () const override
 Returns "decay" at the quadpoint (porosity = a + (b - a) * exp(decay)) More...
 
virtual Real ddecayQp_dvar (unsigned pvar) const override
 d(decay)/d(PorousFlow variable pvar) More...
 
virtual RealGradient ddecayQp_dgradvar (unsigned pvar) const override
 d(decay)/d(grad(PorousFlow variable pvar)) More...
 
virtual void initQpStatefulProperties () override
 
virtual void computeQpProperties () override
 

Protected Attributes

const bool _mechanical
 Porosity is a function of volumetric strain. More...
 
const bool _fluid
 Porosity is a function of effective porepressure. More...
 
const bool _thermal
 Porosity is a function of temperature. More...
 
const bool _chemical
 Porosity is a function of chemistry. More...
 
const VariableValue_phi0
 Porosity at zero strain and zero porepressure and zero temperature. More...
 
const Real _biot
 Biot coefficient. More...
 
const Real _exp_coeff
 Thermal expansion coefficient of the solid porous skeleton. More...
 
const Real _solid_bulk
 Drained bulk modulus of the porous skeleton. More...
 
const Real _coeff
 Short-hand number (biot-1)/solid_bulk. More...
 
const VariableValue_t_reference
 Reference temperature. More...
 
const VariableValue_p_reference
 Reference porepressure. More...
 
const unsigned _num_c_ref
 Number of reference mineral concentrations provided by user. More...
 
std::vector< const VariableValue * > _c_reference
 Reference mineral concentrations. More...
 
const unsigned _num_initial_c
 Number of reference mineral concentrations provided by user. More...
 
std::vector< const VariableValue * > _initial_c
 Reference mineral concentrations. More...
 
std::vector< Real_c_weights
 Weights for the mineral concentrations. More...
 
const MaterialProperty< Real > *const _porosity_old
 Old value of porosity. More...
 
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)) More...
 
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)) More...
 
const MaterialProperty< Real > *const _pf
 Effective porepressure at the quadpoints or nodes. More...
 
const MaterialProperty< std::vector< Real > > *const _dpf_dvar
 d(effective porepressure)/(d porflow variable) More...
 
const MaterialProperty< Real > *const _temperature
 Temperature at the quadpoints or nodes. More...
 
const MaterialProperty< std::vector< Real > > *const _dtemperature_dvar
 d(temperature)/(d porflow variable) More...
 
const MaterialProperty< std::vector< Real > > *const _mineral_conc_old
 Old value of mineral concentration at the quadpoints or nodes. More...
 
const MaterialProperty< std::vector< Real > > *const _reaction_rate
 Reaction rate of mineralisation. More...
 
const MaterialProperty< std::vector< std::vector< Real > > > *const _dreaction_rate_dvar
 d(reaction_rate_conc)/d(porflow variable) More...
 
const unsigned int _aq_ph
 Aqueous phase number. More...
 
const MaterialProperty< std::vector< Real > > *const _saturation
 Saturation. More...
 
const MaterialProperty< std::vector< std::vector< Real > > > *const _dsaturation_dvar
 d(saturation)/d(PorousFlow var) More...
 
const bool _strain_at_nearest_qp
 When calculating nodal porosity, use the strain at the nearest quadpoint to the node. More...
 
const bool _ensure_positive
 for decayQp() > 0, porosity can be negative when using porosity = a + (b - a) * exp(decay). More...
 
GenericMaterialProperty< Real, is_ad > & _porosity
 Computed porosity at the nodes or quadpoints. More...
 
MaterialProperty< std::vector< Real > > *const _dporosity_dvar
 d(porosity)/d(PorousFlow variable) More...
 
MaterialProperty< std::vector< RealGradient > > *const _dporosity_dgradvar
 d(porosity)/d(grad PorousFlow variable) 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

Material designed to provide the porosity in PorousFlow simulations chemistry + biot + (phi0 - reference_chemistry - biot) * exp(-vol_strain.

Definition at line 20 of file PorousFlowPorosity.h.

Constructor & Destructor Documentation

◆ PorousFlowPorosity()

PorousFlowPorosity::PorousFlowPorosity ( const InputParameters parameters)

Definition at line 63 of file PorousFlowPorosity.C.

65 
66  _mechanical(getParam<bool>("mechanical")),
67  _fluid(getParam<bool>("fluid")),
68  _thermal(getParam<bool>("thermal")),
69  _chemical(getParam<bool>("chemical")),
70  _phi0(coupledValue("porosity_zero")),
71  _biot(getParam<Real>("biot_coefficient")),
72  _exp_coeff(isParamValid("thermal_expansion_coeff") ? getParam<Real>("thermal_expansion_coeff")
73  : 0.0),
74  _solid_bulk(isParamValid("solid_bulk") ? getParam<Real>("solid_bulk")
75  : std::numeric_limits<Real>::max()),
76  _coeff(isParamValid("biot_coefficient_prime")
77  ? (getParam<Real>("biot_coefficient_prime") - 1.0) / _solid_bulk
78  : (_biot - 1.0) / _solid_bulk),
79 
80  _t_reference(_nodal_material ? coupledDofValues("reference_temperature")
81  : coupledValue("reference_temperature")),
82  _p_reference(_nodal_material ? coupledDofValues("reference_porepressure")
83  : coupledValue("reference_porepressure")),
84  _num_c_ref(coupledComponents("reference_chemistry")),
86  _num_initial_c(coupledComponents("initial_mineral_concentrations")),
88  _c_weights(isParamValid("chemical_weights") ? getParam<std::vector<Real>>("chemical_weights")
89  : std::vector<Real>(_num_c_ref, 1.0)),
90 
91  _porosity_old(_chemical ? (_nodal_material
92  ? &getMaterialPropertyOld<Real>("PorousFlow_porosity_nodal")
93  : &getMaterialPropertyOld<Real>("PorousFlow_porosity_qp"))
94  : nullptr),
95  _vol_strain_qp(_mechanical ? &getMaterialProperty<Real>("PorousFlow_total_volumetric_strain_qp")
96  : nullptr),
97  _dvol_strain_qp_dvar(_mechanical ? &getMaterialProperty<std::vector<RealGradient>>(
98  "dPorousFlow_total_volumetric_strain_qp_dvar")
99  : nullptr),
100 
101  _pf(_fluid ? (_nodal_material
102  ? &getMaterialProperty<Real>("PorousFlow_effective_fluid_pressure_nodal")
103  : &getMaterialProperty<Real>("PorousFlow_effective_fluid_pressure_qp"))
104  : nullptr),
105  _dpf_dvar(_fluid ? (_nodal_material ? &getMaterialProperty<std::vector<Real>>(
106  "dPorousFlow_effective_fluid_pressure_nodal_dvar")
107  : &getMaterialProperty<std::vector<Real>>(
108  "dPorousFlow_effective_fluid_pressure_qp_dvar"))
109  : nullptr),
110 
112  ? (_nodal_material ? &getMaterialProperty<Real>("PorousFlow_temperature_nodal")
113  : &getMaterialProperty<Real>("PorousFlow_temperature_qp"))
114  : nullptr),
116  _thermal
117  ? (_nodal_material
118  ? &getMaterialProperty<std::vector<Real>>("dPorousFlow_temperature_nodal_dvar")
119  : &getMaterialProperty<std::vector<Real>>("dPorousFlow_temperature_qp_dvar"))
120  : nullptr),
121 
122  _mineral_conc_old(_chemical ? (_nodal_material ? &getMaterialPropertyOld<std::vector<Real>>(
123  "PorousFlow_mineral_concentration_nodal")
124  : &getMaterialPropertyOld<std::vector<Real>>(
125  "PorousFlow_mineral_concentration_qp"))
126  : nullptr),
127  _reaction_rate(_chemical ? (_nodal_material ? &getMaterialProperty<std::vector<Real>>(
128  "PorousFlow_mineral_reaction_rate_nodal")
129  : &getMaterialProperty<std::vector<Real>>(
130  "PorousFlow_mineral_reaction_rate_qp"))
131  : nullptr),
132  _dreaction_rate_dvar(_chemical ? (_nodal_material
133  ? &getMaterialProperty<std::vector<std::vector<Real>>>(
134  "dPorousFlow_mineral_reaction_rate_nodal_dvar")
135  : &getMaterialProperty<std::vector<std::vector<Real>>>(
136  "dPorousFlow_mineral_reaction_rate_qp_dvar"))
137  : nullptr),
138  _aq_ph(_dictator.aqueousPhaseNumber()),
140  ? (_nodal_material
141  ? &getMaterialProperty<std::vector<Real>>("PorousFlow_saturation_nodal")
142  : &getMaterialProperty<std::vector<Real>>("PorousFlow_saturation_qp"))
143  : nullptr),
145  ? (_nodal_material ? &getMaterialProperty<std::vector<std::vector<Real>>>(
146  "dPorousFlow_saturation_nodal_dvar")
147  : &getMaterialProperty<std::vector<std::vector<Real>>>(
148  "dPorousFlow_saturation_qp_dvar"))
149  : nullptr)
150 {
151  if (_thermal && !isParamValid("thermal_expansion_coeff"))
152  mooseError("PorousFlowPorosity: When thermal=true you must provide a thermal_expansion_coeff");
153  if (_fluid && !isParamValid("solid_bulk"))
154  mooseError("PorousFlowPorosity: When fluid=true you must provide a solid_bulk");
155  if (_chemical && _num_c_ref != _dictator.numAqueousKinetic())
156  mooseError("PorousFlowPorosity: When chemical=true you must provide the reference_chemistry "
157  "values. The Dictator proclaims there should be ",
158  _dictator.numAqueousKinetic(),
159  " of these");
160  if (_chemical && _num_initial_c != _dictator.numAqueousKinetic())
161  mooseError("PorousFlowPorosity: When chemical=true you must provide the "
162  "initial_mineral_concentrations. "
163  "The Dictator proclaims there should be ",
164  _dictator.numAqueousKinetic(),
165  " of these");
166  if (_chemical && _c_weights.size() != _dictator.numAqueousKinetic())
167  mooseError(
168  "PorousFlowPorosity: When chemical=true you must provde the correct number of "
169  "chemical_weights (which the Dictator knows is ",
170  _dictator.numAqueousKinetic(),
171  ") or do not provide any chemical_weights and use the default value of 1 for each mineral");
172 
173  for (unsigned i = 0; i < _num_c_ref; ++i)
174  {
175  _c_reference[i] = (_nodal_material ? &coupledDofValues("reference_chemistry", i)
176  : &coupledValue("reference_chemistry", i));
177  _initial_c[i] = (_nodal_material ? &coupledDofValues("initial_mineral_concentrations", i)
178  : &coupledValue("initial_mineral_concentrations", i));
179  }
180 }
const Real _coeff
Short-hand number (biot-1)/solid_bulk.
const MaterialProperty< std::vector< Real > > *const _dtemperature_dvar
d(temperature)/(d porflow variable)
const Real _biot
Biot coefficient.
const VariableValue & _t_reference
Reference temperature.
const VariableValue & _phi0
Porosity at zero strain and zero porepressure and zero temperature.
const MaterialProperty< std::vector< Real > > *const _reaction_rate
Reaction rate of mineralisation.
const Real _solid_bulk
Drained bulk modulus of the porous skeleton.
void mooseError(Args &&... args)
const MaterialProperty< std::vector< Real > > *const _saturation
Saturation.
const MaterialProperty< std::vector< Real > > *const _mineral_conc_old
Old value of mineral concentration at the quadpoints or nodes.
const bool _thermal
Porosity is a function of temperature.
std::vector< Real > _c_weights
Weights for the mineral concentrations.
const bool _chemical
Porosity is a function of chemistry.
const MaterialProperty< std::vector< std::vector< Real > > > *const _dreaction_rate_dvar
d(reaction_rate_conc)/d(porflow variable)
const unsigned _num_c_ref
Number of reference mineral concentrations provided by user.
const MaterialProperty< Real > *const _porosity_old
Old value of porosity.
PorousFlowPorosityExponentialBase(const InputParameters &parameters)
const unsigned int _aq_ph
Aqueous phase number.
const MaterialProperty< Real > *const _pf
Effective porepressure at the quadpoints or nodes.
std::vector< const VariableValue * > _initial_c
Reference mineral concentrations.
const MaterialProperty< std::vector< std::vector< Real > > > *const _dsaturation_dvar
d(saturation)/d(PorousFlow var)
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 Real _exp_coeff
Thermal expansion coefficient of the solid porous skeleton.
const VariableValue & _p_reference
Reference porepressure.
const unsigned _num_initial_c
Number of reference mineral concentrations provided by user.
const bool _fluid
Porosity is a function of effective porepressure.
std::vector< const VariableValue * > _c_reference
Reference mineral concentrations.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
const bool _mechanical
Porosity is a function of volumetric strain.
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&#39;ll always be pointing to the same address after initialization (like a reference))
const MaterialProperty< Real > *const _vol_strain_qp
Strain (first const means we never want to dereference and change the value, second means we&#39;ll alway...

Member Function Documentation

◆ atNegInfinityQp()

Real PorousFlowPorosity::atNegInfinityQp ( ) const
overrideprotectedvirtual

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

Implements PorousFlowPorosityExponentialBase.

Definition at line 183 of file PorousFlowPorosity.C.

184 {
185  /*
186  *
187  * Note the use of the OLD value of porosity here.
188  * This strategy, which breaks the cyclic dependency between porosity
189  * and mineral concentration, is used in
190  * Kernel: PorousFlowPreDis
191  * Material: PorousFlowPorosity
192  * Material: PorousFlowAqueousPreDisChemistry
193  * Material: PorousFlowAqueousPreDisMineral
194  *
195  */
196  Real result = _biot;
197  if (_chemical)
198  {
199  if (_t_step == 0 && !_app.isRestarting())
200  for (unsigned i = 0; i < _num_c_ref; ++i)
201  result -= _c_weights[i] * (*_initial_c[i])[_qp];
202  else
203  for (unsigned i = 0; i < _num_c_ref; ++i)
204  result -= _c_weights[i] * ((*_mineral_conc_old)[_qp][i] + _dt * (*_porosity_old)[_qp] *
205  (*_saturation)[_qp][_aq_ph] *
206  (*_reaction_rate)[_qp][i]);
207  }
208  return result;
209 }
const Real _biot
Biot coefficient.
const MaterialProperty< std::vector< Real > > *const _mineral_conc_old
Old value of mineral concentration at the quadpoints or nodes.
std::vector< Real > _c_weights
Weights for the mineral concentrations.
const bool _chemical
Porosity is a function of chemistry.
const unsigned _num_c_ref
Number of reference mineral concentrations provided by user.
const unsigned int _aq_ph
Aqueous phase number.
std::vector< const VariableValue * > _initial_c
Reference mineral concentrations.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ atZeroQp()

Real PorousFlowPorosity::atZeroQp ( ) const
overrideprotectedvirtual

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

Implements PorousFlowPorosityExponentialBase.

Definition at line 224 of file PorousFlowPorosity.C.

225 {
226  // note the [0] below: _phi0 is a constant monomial and we use [0] regardless of _nodal_material
227  Real result = _phi0[0];
228  if (_chemical)
229  {
230  if (_t_step == 0 && !_app.isRestarting())
231  for (unsigned i = 0; i < _num_c_ref; ++i)
232  result -= _c_weights[i] * ((*_initial_c[i])[_qp] - (*_c_reference[i])[_qp]);
233  else
234  for (unsigned i = 0; i < _num_c_ref; ++i)
235  result -= _c_weights[i] * ((*_mineral_conc_old)[_qp][i] +
236  _dt * (*_porosity_old)[_qp] * (*_saturation)[_qp][_aq_ph] *
237  (*_reaction_rate)[_qp][i] -
238  (*_c_reference[i])[_qp]);
239  }
240  return result;
241 }
const VariableValue & _phi0
Porosity at zero strain and zero porepressure and zero temperature.
const MaterialProperty< std::vector< Real > > *const _mineral_conc_old
Old value of mineral concentration at the quadpoints or nodes.
std::vector< Real > _c_weights
Weights for the mineral concentrations.
const bool _chemical
Porosity is a function of chemistry.
const unsigned _num_c_ref
Number of reference mineral concentrations provided by user.
const unsigned int _aq_ph
Aqueous phase number.
std::vector< const VariableValue * > _initial_c
Reference mineral concentrations.
std::vector< const VariableValue * > _c_reference
Reference mineral concentrations.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ computeQpProperties()

void PorousFlowPorosityExponentialBase::computeQpProperties ( )
overrideprotectedvirtualinherited

Reimplemented in PorousFlowPorosityHMBiotModulus.

Definition at line 61 of file PorousFlowPorosityExponentialBase.C.

62 {
63  const Real a = atNegInfinityQp();
64  const Real b = atZeroQp();
65  const Real decay = decayQp();
66  Real exp_term = 1.0; // set appropriately below
67 
68  Real deriv = 0.0; // = d(porosity)/d(decay)
69  if (decay <= 0.0 || !_ensure_positive)
70  {
71  exp_term = std::exp(decay);
72  _porosity[_qp] = a + (b - a) * exp_term;
73  deriv = _porosity[_qp] - a;
74  }
75  else
76  {
77  const Real c = std::log(a / (a - b));
78  const Real expx = std::exp(-decay / c);
79  // note that at decay = 0, we have expx = 1, so porosity = a + b - a = b
80  // and at decay = infinity, expx = 0, so porosity = a + (b - a) * a / (a - b) = 0
81  exp_term = std::exp(c * (1.0 - expx));
82  _porosity[_qp] = a + (b - a) * exp_term;
83  deriv = (_porosity[_qp] - a) * expx;
84  }
85 
86  (*_dporosity_dvar)[_qp].resize(_num_var);
87  (*_dporosity_dgradvar)[_qp].resize(_num_var);
88  for (unsigned int v = 0; v < _num_var; ++v)
89  {
90  (*_dporosity_dvar)[_qp][v] = ddecayQp_dvar(v) * deriv;
91  (*_dporosity_dgradvar)[_qp][v] = ddecayQp_dgradvar(v) * deriv;
92 
93  const Real da = datNegInfinityQp(v);
94  const Real db = datZeroQp(v);
95  (*_dporosity_dvar)[_qp][v] += da * (1 - exp_term) + db * exp_term;
96 
97  if (!(decay <= 0.0 || !_ensure_positive))
98  {
99  const Real c = std::log(a / (a - b));
100  const Real expx = std::exp(-decay / c);
101  const Real dc = (a - b) * (da * b / a - db) / std::pow(a, 2);
102  (*_dporosity_dvar)[_qp][v] += (b - a) * exp_term * dc * (1 - expx - expx / c);
103  }
104  }
105 }
const bool _ensure_positive
for decayQp() > 0, porosity can be negative when using porosity = a + (b - a) * exp(decay).
virtual Real atZeroQp() const =0
Returns "b" at the quadpoint (porosity = a + (b - a) * exp(decay))
const GeochemicalDatabaseReader db("database/moose_testdb.json", true, true, false)
Real deriv(unsigned n, unsigned alpha, unsigned beta, Real x)
virtual Real datZeroQp(unsigned pvar) const =0
d(a)/d(PorousFlow variable pvar)
virtual Real decayQp() const =0
Returns "decay" at the quadpoint (porosity = a + (b - a) * exp(decay))
const unsigned int _num_var
Number of PorousFlow variables.
virtual Real atNegInfinityQp() const =0
Returns "a" at the quadpoint (porosity = a + (b - a) * exp(decay))
virtual RealGradient ddecayQp_dgradvar(unsigned pvar) const =0
d(decay)/d(grad(PorousFlow variable pvar))
GenericMaterialProperty< Real, is_ad > & _porosity
Computed porosity at the nodes or quadpoints.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
static const std::string v
Definition: NS.h:82
virtual Real datNegInfinityQp(unsigned pvar) const =0
d(a)/d(PorousFlow variable pvar)
virtual Real ddecayQp_dvar(unsigned pvar) const =0
d(decay)/d(PorousFlow variable pvar)
MooseUnits pow(const MooseUnits &, int)

◆ datNegInfinityQp()

Real PorousFlowPorosity::datNegInfinityQp ( unsigned  pvar) const
overrideprotectedvirtual

d(a)/d(PorousFlow variable pvar)

Implements PorousFlowPorosityExponentialBase.

Definition at line 212 of file PorousFlowPorosity.C.

213 {
214  Real result = 0.0;
215  if (_chemical && (_t_step >= 1 || _app.isRestarting()))
216  for (unsigned i = 0; i < _num_c_ref; ++i)
217  result -= _c_weights[i] * _dt * (*_porosity_old)[_qp] *
218  ((*_saturation)[_qp][_aq_ph] * (*_dreaction_rate_dvar)[_qp][i][pvar] +
219  (*_dsaturation_dvar)[_qp][_aq_ph][pvar] * (*_reaction_rate)[_qp][i]);
220  return result;
221 }
std::vector< Real > _c_weights
Weights for the mineral concentrations.
const bool _chemical
Porosity is a function of chemistry.
const unsigned _num_c_ref
Number of reference mineral concentrations provided by user.
const MaterialProperty< Real > *const _porosity_old
Old value of porosity.
const unsigned int _aq_ph
Aqueous phase number.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ datZeroQp()

Real PorousFlowPorosity::datZeroQp ( unsigned  pvar) const
overrideprotectedvirtual

d(a)/d(PorousFlow variable pvar)

Implements PorousFlowPorosityExponentialBase.

Definition at line 244 of file PorousFlowPorosity.C.

245 {
246  Real result = 0.0;
247  if (_chemical && (_t_step >= 1 || _app.isRestarting()))
248  for (unsigned i = 0; i < _num_c_ref; ++i)
249  result -= _c_weights[i] * _dt * (*_porosity_old)[_qp] *
250  ((*_saturation)[_qp][_aq_ph] * (*_dreaction_rate_dvar)[_qp][i][pvar] +
251  (*_dsaturation_dvar)[_qp][_aq_ph][pvar] * (*_reaction_rate)[_qp][i]);
252  return result;
253 }
std::vector< Real > _c_weights
Weights for the mineral concentrations.
const bool _chemical
Porosity is a function of chemistry.
const unsigned _num_c_ref
Number of reference mineral concentrations provided by user.
const MaterialProperty< Real > *const _porosity_old
Old value of porosity.
const unsigned int _aq_ph
Aqueous phase number.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ ddecayQp_dgradvar()

RealGradient PorousFlowPorosity::ddecayQp_dgradvar ( unsigned  pvar) const
overrideprotectedvirtual

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

Implements PorousFlowPorosityExponentialBase.

Definition at line 295 of file PorousFlowPorosity.C.

296 {
297  RealGradient result(0.0, 0.0, 0.0);
298  if (_mechanical)
299  {
300  const unsigned qp_to_use =
301  (_nodal_material && (_bnd || _strain_at_nearest_qp) ? nearestQP(_qp) : _qp);
302  result += -(*_dvol_strain_qp_dvar)[qp_to_use][pvar];
303  }
304  return result;
305 }
const bool _strain_at_nearest_qp
When calculating nodal porosity, use the strain at the nearest quadpoint to the node.
const bool _mechanical
Porosity is a function of volumetric strain.

◆ ddecayQp_dvar()

Real PorousFlowPorosity::ddecayQp_dvar ( unsigned  pvar) const
overrideprotectedvirtual

d(decay)/d(PorousFlow variable pvar)

Implements PorousFlowPorosityExponentialBase.

Definition at line 281 of file PorousFlowPorosity.C.

282 {
283  Real result = 0.0;
284 
285  if (_thermal)
286  result += _exp_coeff * (*_dtemperature_dvar)[_qp][pvar];
287 
288  if (_fluid)
289  result += _coeff * (*_dpf_dvar)[_qp][pvar];
290 
291  return result;
292 }
const Real _coeff
Short-hand number (biot-1)/solid_bulk.
const bool _thermal
Porosity is a function of temperature.
const Real _exp_coeff
Thermal expansion coefficient of the solid porous skeleton.
const bool _fluid
Porosity is a function of effective porepressure.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ decayQp()

Real PorousFlowPorosity::decayQp ( ) const
overrideprotectedvirtual

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

Implements PorousFlowPorosityExponentialBase.

Definition at line 256 of file PorousFlowPorosity.C.

257 {
258  Real result = 0.0;
259 
260  if (_thermal)
261  result += _exp_coeff * ((*_temperature)[_qp] - _t_reference[_qp]);
262 
263  if (_fluid)
264  result += _coeff * ((*_pf)[_qp] - _p_reference[_qp]);
265 
266  if (_mechanical)
267  {
268  // Note that in the following _strain[_qp] is evaluated at q quadpoint
269  // So _porosity_nodal[_qp], which should be the nodal value of porosity
270  // actually uses the strain at a quadpoint. This
271  // is OK for LINEAR elements, as strain is constant over the element anyway.
272  const unsigned qp_to_use =
273  (_nodal_material && (_bnd || _strain_at_nearest_qp) ? nearestQP(_qp) : _qp);
274  result += -(*_vol_strain_qp)[qp_to_use];
275  }
276 
277  return result;
278 }
const Real _coeff
Short-hand number (biot-1)/solid_bulk.
const VariableValue & _t_reference
Reference temperature.
const bool _strain_at_nearest_qp
When calculating nodal porosity, use the strain at the nearest quadpoint to the node.
const bool _thermal
Porosity is a function of temperature.
const Real _exp_coeff
Thermal expansion coefficient of the solid porous skeleton.
const VariableValue & _p_reference
Reference porepressure.
const bool _fluid
Porosity is a function of effective porepressure.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
const bool _mechanical
Porosity is a function of volumetric strain.

◆ initQpStatefulProperties()

void PorousFlowPorosityExponentialBase::initQpStatefulProperties ( )
overrideprotectedvirtualinherited

Definition at line 43 of file PorousFlowPorosityExponentialBase.C.

44 {
45  const Real a = atNegInfinityQp();
46  const Real b = atZeroQp();
47  mooseAssert(a > b, "PorousFlowPorosityExponentialBase a must be larger than b");
48  const Real decay = decayQp();
49 
50  if (decay <= 0.0 || !_ensure_positive)
51  _porosity[_qp] = a + (b - a) * std::exp(decay);
52  else
53  {
54  const Real c = std::log(a / (a - b));
55  const Real expx = std::exp(-decay / c);
56  _porosity[_qp] = a + (b - a) * std::exp(c * (1.0 - expx));
57  }
58 }
const bool _ensure_positive
for decayQp() > 0, porosity can be negative when using porosity = a + (b - a) * exp(decay).
virtual Real atZeroQp() const =0
Returns "b" at the quadpoint (porosity = a + (b - a) * exp(decay))
virtual Real decayQp() const =0
Returns "decay" at the quadpoint (porosity = a + (b - a) * exp(decay))
virtual Real atNegInfinityQp() const =0
Returns "a" at the quadpoint (porosity = a + (b - a) * exp(decay))
GenericMaterialProperty< Real, is_ad > & _porosity
Computed porosity at the nodes or quadpoints.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ validParams()

InputParameters PorousFlowPorosity::validParams ( )
static

Definition at line 15 of file PorousFlowPorosity.C.

Referenced by PorousFlowPorosityHMBiotModulus::validParams().

16 {
18  params.addParam<bool>(
19  "mechanical", false, "If true, porosity will be a function of total volumetric strain");
20  params.addParam<bool>(
21  "fluid", false, "If true, porosity will be a function of effective porepressure");
22  params.addParam<bool>("thermal", false, "If true, porosity will be a function of temperature");
23  params.addParam<bool>("chemical", false, "If true, porosity will be a function of precipitate");
24  params.addRequiredCoupledVar("porosity_zero",
25  "The porosity at zero volumetric strain and "
26  "reference temperature and reference effective "
27  "porepressure and reference chemistry. This must be a real number "
28  "or a constant monomial variable (not a linear lagrange or other "
29  "type of variable)");
30  params.addParam<Real>("thermal_expansion_coeff",
31  "Volumetric thermal expansion coefficient of the drained porous solid "
32  "skeleton (only used if thermal=true)");
33  params.addRangeCheckedParam<Real>(
34  "biot_coefficient", 1, "biot_coefficient>=0 & biot_coefficient<=1", "Biot coefficient");
35  params.addParam<Real>("biot_coefficient_prime",
36  "Biot coefficient that appears in the term (biot_coefficient_prime - 1) * "
37  "(P - reference_porepressure) / solid_bulk. If not provided, this "
38  "defaults to the standard biot_coefficient");
39  params.addRangeCheckedParam<Real>(
40  "solid_bulk",
41  "solid_bulk>0",
42  "Bulk modulus of the drained porous solid skeleton (only used if fluid=true)");
43  params.addCoupledVar(
44  "reference_temperature", 0.0, "Reference temperature (only used if thermal=true)");
45  params.addCoupledVar(
46  "reference_porepressure", 0.0, "Reference porepressure (only used if fluid=true)");
47  params.addCoupledVar("reference_chemistry",
48  "Reference values of the solid mineral concentrations "
49  "(m^3(precipitate)/m^3(porous material)), entered as "
50  "a vector (one value per mineral). (Only used if chemical=true)");
51  params.addCoupledVar(
52  "initial_mineral_concentrations",
53  "Initial mineral concentrations (m^3(precipitate)/m^3(porous material)), entered as "
54  "a vector (one value per mineral). (Only used if chemical=true)");
55  params.addParam<std::vector<Real>>("chemical_weights",
56  "When chemical=true, porosity is a linear combination of the "
57  "solid mineral concentrations multiplied by these weights. "
58  "Default=1 for all minerals.");
59  params.addClassDescription("This Material calculates the porosity PorousFlow simulations");
60  return params;
61 }
void addParam(const std::string &name, const std::initializer_list< typename T::value_type > &value, const std::string &doc_string)
void addCoupledVar(const std::string &name, const std::string &doc_string)
void addRequiredCoupledVar(const std::string &name, const std::string &doc_string)
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
void addClassDescription(const std::string &doc_string)
void addRangeCheckedParam(const std::string &name, const T &value, const std::string &parsed_function, const std::string &doc_string)

Member Data Documentation

◆ _aq_ph

const unsigned int PorousFlowPorosity::_aq_ph
protected

Aqueous phase number.

Definition at line 115 of file PorousFlowPorosity.h.

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

◆ _biot

const Real PorousFlowPorosity::_biot
protected

Biot coefficient.

Definition at line 52 of file PorousFlowPorosity.h.

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

◆ _c_reference

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

Reference mineral concentrations.

Definition at line 73 of file PorousFlowPorosity.h.

Referenced by atZeroQp(), and PorousFlowPorosity().

◆ _c_weights

std::vector<Real> PorousFlowPorosity::_c_weights
protected

Weights for the mineral concentrations.

Definition at line 82 of file PorousFlowPorosity.h.

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

◆ _chemical

const bool PorousFlowPorosity::_chemical
protected

Porosity is a function of chemistry.

Definition at line 46 of file PorousFlowPorosity.h.

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

◆ _coeff

const Real PorousFlowPorosity::_coeff
protected

Short-hand number (biot-1)/solid_bulk.

Definition at line 61 of file PorousFlowPorosity.h.

Referenced by ddecayQp_dvar(), and decayQp().

◆ _dpf_dvar

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

d(effective porepressure)/(d porflow variable)

Definition at line 97 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 PorousFlowPorosityLinear::computeQpProperties(), and PorousFlowPorosityHMBiotModulus::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 PorousFlowPorosityLinear::computeQpProperties(), and PorousFlowPorosityHMBiotModulus::computeQpProperties().

◆ _dreaction_rate_dvar

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

d(reaction_rate_conc)/d(porflow variable)

Definition at line 112 of file PorousFlowPorosity.h.

◆ _dsaturation_dvar

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

d(saturation)/d(PorousFlow var)

Definition at line 121 of file PorousFlowPorosity.h.

◆ _dtemperature_dvar

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

d(temperature)/(d porflow variable)

Definition at line 103 of file PorousFlowPorosity.h.

◆ _dvol_strain_qp_dvar

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

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 91 of file PorousFlowPorosity.h.

◆ _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
protected

Thermal expansion coefficient of the solid porous skeleton.

Definition at line 55 of file PorousFlowPorosity.h.

Referenced by ddecayQp_dvar(), and decayQp().

◆ _fluid

const bool PorousFlowPorosity::_fluid
protected

Porosity is a function of effective porepressure.

Definition at line 40 of file PorousFlowPorosity.h.

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

◆ _initial_c

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

Reference mineral concentrations.

Definition at line 79 of file PorousFlowPorosity.h.

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

◆ _mechanical

const bool PorousFlowPorosity::_mechanical
protected

Porosity is a function of volumetric strain.

Definition at line 37 of file PorousFlowPorosity.h.

Referenced by ddecayQp_dgradvar(), and decayQp().

◆ _mineral_conc_old

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

Old value of mineral concentration at the quadpoints or nodes.

Definition at line 106 of file PorousFlowPorosity.h.

Referenced by atNegInfinityQp(), and atZeroQp().

◆ _num_c_ref

const unsigned PorousFlowPorosity::_num_c_ref
protected

Number of reference mineral concentrations provided by user.

Definition at line 70 of file PorousFlowPorosity.h.

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

◆ _num_components

const unsigned int PorousFlowMaterialVectorBase::_num_components
protectedinherited

◆ _num_initial_c

const unsigned PorousFlowPorosity::_num_initial_c
protected

Number of reference mineral concentrations provided by user.

Definition at line 76 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosity().

◆ _num_phases

const unsigned int PorousFlowMaterialVectorBase::_num_phases
protectedinherited

◆ _num_var

const unsigned int PorousFlowMaterialVectorBase::_num_var
protectedinherited

◆ _p_reference

const VariableValue& PorousFlowPorosity::_p_reference
protected

Reference porepressure.

Definition at line 67 of file PorousFlowPorosity.h.

Referenced by decayQp().

◆ _pf

const MaterialProperty<Real>* const PorousFlowPorosity::_pf
protected

Effective porepressure at the quadpoints or nodes.

Definition at line 94 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosityHMBiotModulus::computeQpProperties().

◆ _phi0

const VariableValue& PorousFlowPorosity::_phi0
protected

Porosity at zero strain and zero porepressure and zero temperature.

Definition at line 49 of file PorousFlowPorosity.h.

Referenced by atZeroQp().

◆ _porosity

template<bool is_ad>
GenericMaterialProperty<Real, is_ad>& PorousFlowPorosityBaseTempl< is_ad >::_porosity
protectedinherited

◆ _porosity_old

const MaterialProperty<Real>* const PorousFlowPorosity::_porosity_old
protected

Old value of porosity.

Definition at line 85 of file PorousFlowPorosity.h.

Referenced by datNegInfinityQp(), and datZeroQp().

◆ _reaction_rate

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

Reaction rate of mineralisation.

Definition at line 109 of file PorousFlowPorosity.h.

◆ _saturation

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

Saturation.

Definition at line 118 of file PorousFlowPorosity.h.

◆ _solid_bulk

const Real PorousFlowPorosity::_solid_bulk
protected

Drained bulk modulus of the porous skeleton.

Definition at line 58 of file PorousFlowPorosity.h.

◆ _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 PorousFlowPorosityHMBiotModulus::computeQpProperties(), ddecayQp_dgradvar(), and decayQp().

◆ _t_reference

const VariableValue& PorousFlowPorosity::_t_reference
protected

Reference temperature.

Definition at line 64 of file PorousFlowPorosity.h.

Referenced by decayQp().

◆ _temperature

const MaterialProperty<Real>* const PorousFlowPorosity::_temperature
protected

Temperature at the quadpoints or nodes.

Definition at line 100 of file PorousFlowPorosity.h.

◆ _thermal

const bool PorousFlowPorosity::_thermal
protected

Porosity is a function of temperature.

Definition at line 43 of file PorousFlowPorosity.h.

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

◆ _vol_strain_qp

const MaterialProperty<Real>* const PorousFlowPorosity::_vol_strain_qp
protected

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 88 of file PorousFlowPorosity.h.

Referenced by PorousFlowPorosityHMBiotModulus::computeQpProperties().


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