Porosity (typically the aperture of a lower-dimensional fracture) given by the PorousFlowPorosityLinear relationship plus a thermal-opening term that depends on how much heat the surrounding matrix has gained or lost: More...

#include <PorousFlowPorosityVinsomeWesterveld.h>

Inheritance diagram for PorousFlowPorosityVinsomeWesterveld:

[legend]

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	PorousFlowPorosityVinsomeWesterveld (const InputParameters &parameters)

virtual void	initialSetup () override

const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialProperty (const std::string &name)

const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialPropertyByName (const std::string &name)

Moose::Kokkos::MaterialProperty< U, dimension >	declareKokkosPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

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="")

Moose::Kokkos::MaterialProperty< U, dimension >	getKokkosMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

Moose::Kokkos::MaterialProperty< U, dimension >	getKokkosMaterialPropertyDerivative (const std::string &base, const SymbolName &c1)

Moose::Kokkos::MaterialProperty< U, dimension >	getKokkosMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

Moose::Kokkos::MaterialProperty< U, dimension >	getKokkosMaterialPropertyDerivative (const std::string &base, const SymbolName &c1)

Moose::Kokkos::MaterialProperty< U, dimension >	getKokkosMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)

Moose::Kokkos::MaterialProperty< U, dimension >	getKokkosMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)

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	initQpStatefulProperties () override

virtual void	computeUnboundedPorosityQp () override
	Computes porosity and its derivatives at the current node or quadrature point, without applying the porosity_min lower bound.

virtual void	computeQpProperties () override

void	applyPorosityMin ()
	If _porosity[_qp] < porosity_min, set _porosity[_qp] = porosity_min and multiply the derivatives of porosity by _zero_modifier (the derivative material properties for non-AD objects, and the AD derivatives for AD objects).

Protected Attributes
const Real	_thermal_opening_coeff
	Change of porosity per unit matrix heat content (units 1/K when porosity is an aperture)

const Real	_kappa
	Thermal diffusivity of the matrix in the direction normal to the fracture (m^2/s), assumed constant (independent of temperature and saturation)

const VariableValue &	_T_matrix_ref
	Undisturbed (reference) matrix temperature.

const Real	_start_time
	Time at which the matrix temperature profile is undisturbed.

MaterialProperty< Real > &	_theta
	Temperature change at the fracture face, theta = T - T_matrix_ref.

const MaterialProperty< Real > &	_theta_old
	Old value of theta.

MaterialProperty< Real > &	_heat_content
	Heat content of the adjacent matrix, I = int (T - T_matrix_ref) dz.

const MaterialProperty< Real > &	_heat_content_old
	Old value of the heat content.

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

const VariableValue &	_phi_ref
	Porosity at reference porepressure, temperature and volumetric strain.

const VariableValue &	_P_ref
	Reference effective porepressure.

const VariableValue &	_T_ref
	Reference temperature.

const VariableValue &	_epv_ref
	Reference volumetric strain.

const Real	_P_coeff
	coefficient of effective porepressure

const Real	_T_coeff
	coefficient of temperature

const Real	_epv_coeff
	coefficient of volumetric strain

const bool	_uses_volstrain
	whether epv_coeff has been set

const OptionalMaterialProperty< Real > &	_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 OptionalMaterialProperty< std::vector< RealGradient > > &	_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 bool	_uses_pf
	whether P_coeff has been set

const OptionalMaterialProperty< Real > &	_pf_nodal
	Effective porepressure at the quadpoints or nodes.

const OptionalMaterialProperty< Real > &	_pf_qp

const MaterialProperty< Real > *	_pf

const OptionalMaterialProperty< std::vector< Real > > &	_dpf_dvar_nodal
	d(effective porepressure)/(d porflow variable)

const OptionalMaterialProperty< std::vector< Real > > &	_dpf_dvar_qp

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

const bool	_uses_T
	whether T_coeff has been set

const OptionalMaterialProperty< Real > &	_temperature_nodal
	Temperature at the quadpoints or nodes.

const OptionalMaterialProperty< Real > &	_temperature_qp

const MaterialProperty< Real > *	_temperature

const OptionalMaterialProperty< std::vector< Real > > &	_dtemperature_dvar_nodal
	d(temperature)/(d porflow variable)

const OptionalMaterialProperty< std::vector< Real > > &	_dtemperature_dvar_qp

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

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 Real	_porosity_min
	Minimum allowed porosity (default zero, so porosity is never negative)

const Real	_zero_modifier
	When the porosity_min floor is active, the porosity derivatives are multiplied by this.

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

Porosity (typically the aperture of a lower-dimensional fracture) given by the PorousFlowPorosityLinear relationship plus a thermal-opening term that depends on how much heat the surrounding matrix has gained or lost:

porosity = linear terms - thermal_opening_coeff * I

where I = int_0^infinity (T(z) - T_ref) dz is the heat content (in K.m) of the matrix adjacent to the fracture. I is computed using the semi-analytic method of Vinsome and Westerveld (1980), which assumes the matrix temperature profile normal to the fracture is

T(z, t) - T_ref = (theta + p z + q z^2) exp(-z / d), d = sqrt(kappa t) / 2

with theta = T - T_ref the temperature change at the fracture face. Each time step, p is found from the energy balance and the heat-conduction equation at z = 0, giving

p = (a theta / d + I_old - d^3 (theta - theta_old) / a) / (3 d^2 + a), a = kappa dt I = I_old + a (theta / d - p)

so I is linear in the current temperature, and dI/dT is a scalar.

Definition at line 36 of file PorousFlowPorosityVinsomeWesterveld.h.

Constructor & Destructor Documentation

◆ PorousFlowPorosityVinsomeWesterveld()

PorousFlowPorosityVinsomeWesterveld::PorousFlowPorosityVinsomeWesterveld ( const InputParameters & parameters )

Definition at line 56 of file PorousFlowPorosityVinsomeWesterveld.C.

  : PorousFlowPorosityLinear(parameters),
    _thermal_opening_coeff(getParam<Real>("thermal_opening_coeff")),
    _kappa(getParam<Real>("matrix_thermal_diffusivity")),
    _T_matrix_ref(nodalOrQpValue("matrix_reference_temperature")),
    _start_time(getParam<Real>("vw_start_time")),
    _theta(_nodal_material ? declareProperty<Real>("PorousFlow_vw_theta_nodal")
                           : declareProperty<Real>("PorousFlow_vw_theta_qp")),
    _theta_old(_nodal_material ? getMaterialPropertyOld<Real>("PorousFlow_vw_theta_nodal")
                               : getMaterialPropertyOld<Real>("PorousFlow_vw_theta_qp")),
    _heat_content(_nodal_material ? declareProperty<Real>("PorousFlow_vw_heat_content_nodal")
                                  : declareProperty<Real>("PorousFlow_vw_heat_content_qp")),
    _heat_content_old(_nodal_material
                          ? getMaterialPropertyOld<Real>("PorousFlow_vw_heat_content_nodal")
                          : getMaterialPropertyOld<Real>("PorousFlow_vw_heat_content_qp"))
{
}

Member Function Documentation

◆ applyPorosityMin()

template<bool is_ad>

void PorousFlowPorosityBaseTempl< is_ad >::applyPorosityMin ( )

protectedinherited

If _porosity[_qp] < porosity_min, set _porosity[_qp] = porosity_min and multiply the derivatives of porosity by _zero_modifier (the derivative material properties for non-AD objects, and the AD derivatives for AD objects).

Otherwise do nothing. Derived classes call this at the end of computeQpProperties(), after the unfloored porosity and its derivatives have been computed.

Definition at line 58 of file PorousFlowPorosityBase.C.

{
  if (MetaPhysicL::raw_value(_porosity[_qp]) >= _porosity_min)
    return;
 
  if constexpr (!is_ad)
  {
    _porosity[_qp] = _porosity_min;
    for (unsigned int v = 0; v < _num_var; ++v)
    {
      (*_dporosity_dvar)[_qp][v] *= _zero_modifier;
      (*_dporosity_dgradvar)[_qp][v] *= _zero_modifier;
    }
  }
  else
  {
    // The AD path carries its derivatives inside the value and has no derivative material
    // properties, so the value and the derivatives are softened separately.  Assigning
    // _porosity_min to the ADReal would instead discard the derivatives entirely.
    _porosity[_qp].value() = _porosity_min;
    _porosity[_qp].derivatives() *= _zero_modifier;
  }
}

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

◆ computeQpProperties()

void PorousFlowPorosityLinear::computeQpProperties ( )

overrideprotectedvirtualinherited

Definition at line 121 of file PorousFlowPorosityLinear.C.

{
  computeUnboundedPorosityQp();
  applyPorosityMin();
}

◆ computeUnboundedPorosityQp()

void PorousFlowPorosityVinsomeWesterveld::computeUnboundedPorosityQp ( )

overrideprotectedvirtual

Computes porosity and its derivatives at the current node or quadrature point, without applying the porosity_min lower bound.

Derived classes may override this to add further terms to _porosity and its derivatives; computeQpProperties() then applies porosity_min

Reimplemented from PorousFlowPorosityLinear.

Definition at line 93 of file PorousFlowPorosityVinsomeWesterveld.C.

{
  PorousFlowPorosityLinear::computeUnboundedPorosityQp();
 
  const Real theta = (*_temperature)[_qp] - _T_matrix_ref[_qp];
  _theta[_qp] = theta;
 
  // dI/dT, used for the Jacobian (equal to dI/dtheta, because theta = T - T_matrix_ref)
  Real dI_dT = 0.0;
 
  const Real time = _t - _start_time;
  if (time <= 0.0)
  {
    // Before vw_start_time the matrix is undisturbed, so any difference between the fracture-face
    // temperature and matrix_reference_temperature at vw_start_time is treated as a step change
    _theta[_qp] = 0.0;
    _heat_content[_qp] = 0.0;
  }
  else
  {
    const Real a = _kappa * std::min(_dt, time);
    // Vinsome and Westerveld (1980): the profile is
    // (theta + p z + q z^2) exp(-z / d), with d = sqrt(kappa t) / 2.
    // Imposing the heat-conduction equation at z = 0 and the energy balance
    // dI/dt = kappa * (theta / d - p) over the time step gives p and I below.
    const Real d = 0.5 * std::sqrt(_kappa * time);
    const Real denom = 3.0 * d * d + a;
    const Real p = (a * theta / d + _heat_content_old[_qp] -
                    Utility::pow<3>(d) * (theta - _theta_old[_qp]) / a) /
                   denom;
    _heat_content[_qp] = _heat_content_old[_qp] + a * (theta / d - p);
 
    const Real dp_dT = (a / d - Utility::pow<3>(d) / a) / denom;
    dI_dT = a * (1.0 / d - dp_dT);
  }
 
  _porosity[_qp] -= _thermal_opening_coeff * _heat_content[_qp];
  for (unsigned int pvar = 0; pvar < _num_var; ++pvar)
    (*_dporosity_dvar)[_qp][pvar] -=
        _thermal_opening_coeff * dI_dT * (*_dtemperature_dvar)[_qp][pvar];
}

◆ initialSetup()

void PorousFlowPorosityVinsomeWesterveld::initialSetup ( )

overridevirtual

Reimplemented from PorousFlowPorosityLinear.

Definition at line 76 of file PorousFlowPorosityVinsomeWesterveld.C.

{
  PorousFlowPorosityLinear::initialSetup();
  if (!(_temperature && _dtemperature_dvar))
    mooseError("PorousFlowPorosityVinsomeWesterveld: a temperature Material is required");
}

◆ initQpStatefulProperties()

void PorousFlowPorosityVinsomeWesterveld::initQpStatefulProperties ( )

overrideprotectedvirtual

Reimplemented from PorousFlowPorosityLinear.

Definition at line 84 of file PorousFlowPorosityVinsomeWesterveld.C.

{
  PorousFlowPorosityLinear::initQpStatefulProperties();
  // The matrix starts at its reference temperature
  _theta[_qp] = 0.0;
  _heat_content[_qp] = 0.0;
}

◆ validParams()

InputParameters PorousFlowPorosityVinsomeWesterveld::validParams ( )

static

Definition at line 19 of file PorousFlowPorosityVinsomeWesterveld.C.

{
  InputParameters params = PorousFlowPorosityLinear::validParams();
  params.addRequiredParam<Real>(
      "thermal_opening_coeff",
      "Change in porosity per unit heat content of the adjacent matrix.  The heat content, "
      "I = int (T - matrix_reference_temperature) dz, has units K.m, so when porosity represents "
      "the aperture of a lower-dimensional fracture this coefficient has units 1/K.  Porosity "
      "changes by -thermal_opening_coeff * I, so cooling (I < 0) opens the fracture when this "
      "coefficient is positive");
  params.addRequiredRangeCheckedParam<Real>(
      "matrix_thermal_diffusivity",
      "matrix_thermal_diffusivity > 0",
      "Thermal diffusivity of the matrix adjacent to the "
      "fracture (thermal conductivity divided by "
      "volumetric heat capacity), m^2/s.  This is the diffusivity in the "
      "direction normal to the fracture (or the matrix is "
      "assumed isotropic), and it is assumed constant");
  params.addRequiredCoupledVar(
      "matrix_reference_temperature",
      "Undisturbed temperature of the matrix adjacent to the fracture.  This should usually be a "
      "linear-lagrange variable or a constant");
  params.addParam<Real>(
      "vw_start_time",
      0.0,
      "Time at which the temperature at the fracture face first starts to change "
      "(for instance, when injection begins).  The matrix adjacent to the fracture "
      "is at matrix_reference_temperature at this time, and the thermal-opening "
      "term is zero before it");
  params.addClassDescription(
      "This Material calculates porosity (usually the aperture of a lower-dimensional fracture) "
      "using the PorousFlowPorosityLinear relationship plus a thermal-opening term proportional to "
      "the heat content of the adjacent matrix, which is computed using the semi-analytic method "
      "of Vinsome and Westerveld (1980)");
  return params;
}

Member Data Documentation

◆ _dpf_dvar

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

protectedinherited

Definition at line 83 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::initialSetup().

◆ _dpf_dvar_nodal

const OptionalMaterialProperty<std::vector<Real> >& PorousFlowPorosityLinear::_dpf_dvar_nodal

protectedinherited

d(effective porepressure)/(d porflow variable)

Definition at line 81 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::initialSetup().

◆ _dpf_dvar_qp

const OptionalMaterialProperty<std::vector<Real> >& PorousFlowPorosityLinear::_dpf_dvar_qp

protectedinherited

Definition at line 82 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::initialSetup().

◆ _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 42 of file PorousFlowPorosityBase.h.

Referenced by 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 39 of file PorousFlowPorosityBase.h.

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

◆ _dtemperature_dvar

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

protectedinherited

Definition at line 96 of file PorousFlowPorosityLinear.h.

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

◆ _dtemperature_dvar_nodal

const OptionalMaterialProperty<std::vector<Real> >& PorousFlowPorosityLinear::_dtemperature_dvar_nodal

protectedinherited

d(temperature)/(d porflow variable)

Definition at line 94 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::initialSetup().

◆ _dtemperature_dvar_qp

const OptionalMaterialProperty<std::vector<Real> >& PorousFlowPorosityLinear::_dtemperature_dvar_qp

protectedinherited

Definition at line 95 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::initialSetup().

◆ _dvol_strain_qp_dvar

const OptionalMaterialProperty<std::vector<RealGradient> >& PorousFlowPorosityLinear::_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 70 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::computeUnboundedPorosityQp(), and PorousFlowPorosityLinear::initialSetup().

◆ _epv_coeff

const Real PorousFlowPorosityLinear::_epv_coeff

protectedinherited

coefficient of volumetric strain

Definition at line 61 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::computeUnboundedPorosityQp().

◆ _epv_ref

const VariableValue& PorousFlowPorosityLinear::_epv_ref

protectedinherited

Reference volumetric strain.

Definition at line 52 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::computeUnboundedPorosityQp().

◆ _heat_content

MaterialProperty<Real>& PorousFlowPorosityVinsomeWesterveld::_heat_content

protected

Heat content of the adjacent matrix, I = int (T - T_matrix_ref) dz.

Definition at line 69 of file PorousFlowPorosityVinsomeWesterveld.h.

Referenced by computeUnboundedPorosityQp(), and initQpStatefulProperties().

◆ _heat_content_old

const MaterialProperty<Real>& PorousFlowPorosityVinsomeWesterveld::_heat_content_old

protected

Old value of the heat content.

Definition at line 72 of file PorousFlowPorosityVinsomeWesterveld.h.

Referenced by computeUnboundedPorosityQp().

◆ _kappa

const Real PorousFlowPorosityVinsomeWesterveld::_kappa

protected

Thermal diffusivity of the matrix in the direction normal to the fracture (m^2/s), assumed constant (independent of temperature and saturation)

Definition at line 54 of file PorousFlowPorosityVinsomeWesterveld.h.

Referenced by computeUnboundedPorosityQp().

◆ _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_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(), PorousFlowPorosityHMBiotModulus::computeQpProperties(), PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity::computeQpProperties(), PorousFlowPorosityLinear::computeUnboundedPorosityQp(), computeUnboundedPorosityQp(), PorousFlowAqueousPreDisChemistry::initQpStatefulProperties(), and PorousFlowPorosityLinear::initQpStatefulProperties().

◆ _P_coeff

const Real PorousFlowPorosityLinear::_P_coeff

protectedinherited

coefficient of effective porepressure

Definition at line 55 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::computeUnboundedPorosityQp().

◆ _P_ref

const VariableValue& PorousFlowPorosityLinear::_P_ref

protectedinherited

Reference effective porepressure.

Definition at line 46 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::computeUnboundedPorosityQp().

◆ _pf

const MaterialProperty<Real>* PorousFlowPorosityLinear::_pf

protectedinherited

Definition at line 78 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::initialSetup().

◆ _pf_nodal

const OptionalMaterialProperty<Real>& PorousFlowPorosityLinear::_pf_nodal

protectedinherited

Effective porepressure at the quadpoints or nodes.

Definition at line 76 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::initialSetup().

◆ _pf_qp

const OptionalMaterialProperty<Real>& PorousFlowPorosityLinear::_pf_qp

protectedinherited

Definition at line 77 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::initialSetup().

◆ _phi_ref

const VariableValue& PorousFlowPorosityLinear::_phi_ref

protectedinherited

Porosity at reference porepressure, temperature and volumetric strain.

Definition at line 43 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::computeUnboundedPorosityQp(), and PorousFlowPorosityLinear::initQpStatefulProperties().

◆ _porosity

template<bool is_ad>

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

protectedinherited

Computed porosity at the nodes or quadpoints.

Definition at line 36 of file PorousFlowPorosityBase.h.

Referenced by PorousFlowPorosityHMBiotModulus::computeQpProperties(), PorousFlowPorosityLinear::computeUnboundedPorosityQp(), computeUnboundedPorosityQp(), and PorousFlowPorosityLinear::initQpStatefulProperties().

◆ _porosity_min

template<bool is_ad>

const Real PorousFlowPorosityBaseTempl< is_ad >::_porosity_min

protectedinherited

Minimum allowed porosity (default zero, so porosity is never negative)

Definition at line 45 of file PorousFlowPorosityBase.h.

◆ _start_time

const Real PorousFlowPorosityVinsomeWesterveld::_start_time

protected

Time at which the matrix temperature profile is undisturbed.

Definition at line 60 of file PorousFlowPorosityVinsomeWesterveld.h.

Referenced by computeUnboundedPorosityQp().

◆ _strain_at_nearest_qp

const bool PorousFlowPorosityLinear::_strain_at_nearest_qp

protectedinherited

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

Definition at line 40 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::computeUnboundedPorosityQp().

◆ _T_coeff

const Real PorousFlowPorosityLinear::_T_coeff

protectedinherited

coefficient of temperature

Definition at line 58 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::computeUnboundedPorosityQp().

◆ _T_matrix_ref

const VariableValue& PorousFlowPorosityVinsomeWesterveld::_T_matrix_ref

protected

Undisturbed (reference) matrix temperature.

Definition at line 57 of file PorousFlowPorosityVinsomeWesterveld.h.

Referenced by computeUnboundedPorosityQp().

◆ _T_ref

const VariableValue& PorousFlowPorosityLinear::_T_ref

protectedinherited

Reference temperature.

Definition at line 49 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::computeUnboundedPorosityQp().

◆ _temperature

const MaterialProperty<Real>* PorousFlowPorosityLinear::_temperature

protectedinherited

Definition at line 91 of file PorousFlowPorosityLinear.h.

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

◆ _temperature_nodal

const OptionalMaterialProperty<Real>& PorousFlowPorosityLinear::_temperature_nodal

protectedinherited

Temperature at the quadpoints or nodes.

Definition at line 89 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::initialSetup().

◆ _temperature_qp

const OptionalMaterialProperty<Real>& PorousFlowPorosityLinear::_temperature_qp

protectedinherited

Definition at line 90 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::initialSetup().

◆ _thermal_opening_coeff

const Real PorousFlowPorosityVinsomeWesterveld::_thermal_opening_coeff

protected

Change of porosity per unit matrix heat content (units 1/K when porosity is an aperture)

Definition at line 50 of file PorousFlowPorosityVinsomeWesterveld.h.

Referenced by computeUnboundedPorosityQp().

◆ _theta

MaterialProperty<Real>& PorousFlowPorosityVinsomeWesterveld::_theta

protected

Temperature change at the fracture face, theta = T - T_matrix_ref.

Definition at line 63 of file PorousFlowPorosityVinsomeWesterveld.h.

Referenced by computeUnboundedPorosityQp(), and initQpStatefulProperties().

◆ _theta_old

const MaterialProperty<Real>& PorousFlowPorosityVinsomeWesterveld::_theta_old

protected

Old value of theta.

Definition at line 66 of file PorousFlowPorosityVinsomeWesterveld.h.

Referenced by computeUnboundedPorosityQp().

◆ _uses_pf

const bool PorousFlowPorosityLinear::_uses_pf

protectedinherited

whether P_coeff has been set

Definition at line 73 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::computeUnboundedPorosityQp(), and PorousFlowPorosityLinear::initialSetup().

◆ _uses_T

const bool PorousFlowPorosityLinear::_uses_T

protectedinherited

whether T_coeff has been set

Definition at line 86 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::computeUnboundedPorosityQp(), and PorousFlowPorosityLinear::initialSetup().

◆ _uses_volstrain

const bool PorousFlowPorosityLinear::_uses_volstrain

protectedinherited

whether epv_coeff has been set

Definition at line 64 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::computeUnboundedPorosityQp(), and PorousFlowPorosityLinear::initialSetup().

◆ _vol_strain_qp

const OptionalMaterialProperty<Real>& PorousFlowPorosityLinear::_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 67 of file PorousFlowPorosityLinear.h.

Referenced by PorousFlowPorosityLinear::computeUnboundedPorosityQp(), and PorousFlowPorosityLinear::initialSetup().

◆ _zero_modifier

template<bool is_ad>

const Real PorousFlowPorosityBaseTempl< is_ad >::_zero_modifier

protectedinherited

When the porosity_min floor is active, the porosity derivatives are multiplied by this.

Definition at line 48 of file PorousFlowPorosityBase.h.

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

porous_flow/include/materials/PorousFlowPorosityVinsomeWesterveld.h
porous_flow/src/materials/PorousFlowPorosityVinsomeWesterveld.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

◆ PorousFlowPorosityVinsomeWesterveld()

Member Function Documentation

◆ applyPorosityMin()

◆ computeQpProperties()

◆ computeUnboundedPorosityQp()

◆ initialSetup()

◆ initQpStatefulProperties()

◆ validParams()

Member Data Documentation

◆ _dpf_dvar

◆ _dpf_dvar_nodal

◆ _dpf_dvar_qp

◆ _dporosity_dgradvar

◆ _dporosity_dvar

◆ _dtemperature_dvar

◆ _dtemperature_dvar_nodal

◆ _dtemperature_dvar_qp

◆ _dvol_strain_qp_dvar

◆ _epv_coeff

◆ _epv_ref

◆ _heat_content

◆ _heat_content_old

◆ _kappa

◆ _num_components

◆ _num_phases

◆ _num_var

◆ _P_coeff

◆ _P_ref

◆ _pf

◆ _pf_nodal

◆ _pf_qp

◆ _phi_ref

◆ _porosity

◆ _porosity_min

◆ _start_time

◆ _strain_at_nearest_qp

◆ _T_coeff

◆ _T_matrix_ref

◆ _T_ref

◆ _temperature

◆ _temperature_nodal

◆ _temperature_qp

◆ _thermal_opening_coeff

◆ _theta

◆ _theta_old

◆ _uses_pf

◆ _uses_T

◆ _uses_volstrain

◆ _vol_strain_qp

◆ _zero_modifier