PorousFlowMatrixInternalEnergyTempl< is_ad > Class Template Reference

This material computes internal energy (J/m^3) for a rock matrix assuming constant grain density, specific heat capacity, and a linear relationship with temperature. More...

#include <PorousFlowMatrixInternalEnergy.h>

Inheritance diagram for PorousFlowMatrixInternalEnergyTempl< is_ad >:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

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

Protected Attributes
const Real	_cp
	Specific heat capacity of rock grains. More...
const Real	_density
	Density of rock grains (equals the density of the matrix if porosity=0) More...
const Real	_heat_cap
	Heat capacity = _cp * _density. More...
const GenericMaterialProperty< Real, is_ad > &	_temperature
	Temperature at the nodes. More...
const MaterialProperty< std::vector< Real > > *const	_dtemperature_dvar
	d(temperature at the nodes)/d(PorousFlow variable) More...
GenericMaterialProperty< Real, is_ad > &	_energy
	Matrix internal_energy at the nodes. More...
MaterialProperty< std::vector< Real > > *const	_denergy_dvar
	d(matrix internal energy)/d(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

template<bool is_ad>
class PorousFlowMatrixInternalEnergyTempl< is_ad >

This material computes internal energy (J/m^3) for a rock matrix assuming constant grain density, specific heat capacity, and a linear relationship with temperature.

To get the volumetric heat capacity of the rock in a rock-fluid system, the result must be multiplied by (1 - porosity).

Definition at line 22 of file PorousFlowMatrixInternalEnergy.h.

Constructor & Destructor Documentation

PorousFlowMatrixInternalEnergyTempl()

template<bool is_ad>

PorousFlowMatrixInternalEnergyTempl< is_ad >::PorousFlowMatrixInternalEnergyTempl

(

const InputParameters &

parameters

)

Definition at line 33 of file PorousFlowMatrixInternalEnergy.C.

  : PorousFlowMaterialVectorBase(parameters),
    _cp(getParam<Real>("specific_heat_capacity")),
    _density(getParam<Real>("density")),
    _heat_cap(_cp * _density),
    _temperature(is_ad ? getGenericMaterialProperty<Real, is_ad>("PorousFlow_temperature_qp")
                       : getGenericMaterialProperty<Real, is_ad>("PorousFlow_temperature_nodal")),
    _dtemperature_dvar(
        is_ad ? nullptr
              : &getMaterialProperty<std::vector<Real>>("dPorousFlow_temperature_nodal_dvar")),
    _energy(declareGenericProperty<Real, is_ad>("PorousFlow_matrix_internal_energy_nodal")),
    _denergy_dvar(is_ad ? nullptr
                        : &declareProperty<std::vector<Real>>(
                              "dPorousFlow_matrix_internal_energy_nodal_dvar"))
{
  if (!is_ad && _nodal_material != true)
    mooseError("PorousFlowMatrixInternalEnergy classes are only defined for at_nodes = true");
}

Member Function Documentation

computeQpProperties()

template<bool is_ad>

void PorousFlowMatrixInternalEnergyTempl< is_ad >::computeQpProperties ( )

overrideprotectedvirtual

Definition at line 62 of file PorousFlowMatrixInternalEnergy.C.

{
  _energy[_qp] = _heat_cap * _temperature[_qp];

  if constexpr (!is_ad)
  {
    (*_denergy_dvar)[_qp].assign(_num_var, 0.0);
    for (unsigned v = 0; v < _num_var; ++v)
      (*_denergy_dvar)[_qp][v] = _heat_cap * (*_dtemperature_dvar)[_qp][v];
  }
}

initQpStatefulProperties()

template<bool is_ad>

void PorousFlowMatrixInternalEnergyTempl< is_ad >::initQpStatefulProperties ( )

overrideprotectedvirtual

Definition at line 55 of file PorousFlowMatrixInternalEnergy.C.

{
  _energy[_qp] = _heat_cap * _temperature[_qp];
}

validParams()

template<bool is_ad>

InputParameters PorousFlowMatrixInternalEnergyTempl< is_ad >::validParams ( )

static

Definition at line 17 of file PorousFlowMatrixInternalEnergy.C.

{
  InputParameters params = PorousFlowMaterialVectorBase::validParams();
  params.addRequiredParam<Real>("specific_heat_capacity",
                                "Specific heat capacity of the rock grains (J/kg/K).");
  params.addRequiredParam<Real>("density", "Density of the rock grains");
  params.set<bool>("at_nodes") = is_ad ? false : true;
  params.addPrivateParam<std::string>("pf_material_type", "matrix_internal_energy");
  params.addClassDescription("This Material calculates the internal energy of solid rock grains, "
                             "which is specific_heat_capacity * density * temperature.  Kernels "
                             "multiply this by (1 - porosity) to find the energy density of the "
                             "porous rock in a rock-fluid system");
  return params;
}

Member Data Documentation

_cp

template<bool is_ad>

const Real PorousFlowMatrixInternalEnergyTempl< is_ad >::_cp

protected

Specific heat capacity of rock grains.

Definition at line 34 of file PorousFlowMatrixInternalEnergy.h.

_denergy_dvar

template<bool is_ad>

MaterialProperty<std::vector<Real> >* const PorousFlowMatrixInternalEnergyTempl< is_ad >::_denergy_dvar

protected

d(matrix internal energy)/d(PorousFlow variable)

Definition at line 52 of file PorousFlowMatrixInternalEnergy.h.

_density

template<bool is_ad>

const Real PorousFlowMatrixInternalEnergyTempl< is_ad >::_density

protected

Density of rock grains (equals the density of the matrix if porosity=0)

Definition at line 37 of file PorousFlowMatrixInternalEnergy.h.

_dtemperature_dvar

template<bool is_ad>

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

protected

d(temperature at the nodes)/d(PorousFlow variable)

Definition at line 46 of file PorousFlowMatrixInternalEnergy.h.

_energy

template<bool is_ad>

GenericMaterialProperty<Real, is_ad>& PorousFlowMatrixInternalEnergyTempl< is_ad >::_energy

protected

Matrix internal_energy at the nodes.

Definition at line 49 of file PorousFlowMatrixInternalEnergy.h.

_heat_cap

template<bool is_ad>

const Real PorousFlowMatrixInternalEnergyTempl< is_ad >::_heat_cap

protected

Heat capacity = _cp * _density.

Definition at line 40 of file PorousFlowMatrixInternalEnergy.h.

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

_temperature

template<bool is_ad>

const GenericMaterialProperty<Real, is_ad>& PorousFlowMatrixInternalEnergyTempl< is_ad >::_temperature

protected

Temperature at the nodes.

Definition at line 43 of file PorousFlowMatrixInternalEnergy.h.

---

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

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