Kernel = energy_density * d(volumetric_strain)/dt which is lumped to the nodes. More...

#include <PorousFlowHeatVolumetricExpansion.h>

Inheritance diagram for PorousFlowHeatVolumetricExpansion:

Public Member Functions
	PorousFlowHeatVolumetricExpansion (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpResidual () override

virtual Real	computeQpJacobian () override

virtual Real	computeQpOffDiagJacobian (unsigned int jvar) override

Real	computedEnergyQpJac (unsigned int jvar)
	Derivative of energy part of the residual with respect to the Variable with variable number jvar. More...

Real	computedVolQpJac (unsigned int jvar)
	Derivative of volumetric-strain part of the residual with respect to the Variable with variable number jvar. More...

Protected Attributes
const PorousFlowDictator &	_dictator
	PorousFlowDictator UserObject. More...

const bool	_var_is_porflow_var
	Whether the Variable for this Kernel is a PorousFlow variable according to the Dictator. More...

const unsigned int	_num_phases
	Number of fluid phases. More...

const bool	_fluid_present
	Whether fluid is present. More...

const bool	_strain_at_nearest_qp
	Whether the porosity uses the volumetric strain at the closest quadpoint. More...

const MaterialProperty< Real > &	_porosity
	Porosity. More...

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

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

const MaterialProperty< unsigned int > *const	_nearest_qp
	The nearest qp to the node. More...

const MaterialProperty< Real > &	_rock_energy_nodal
	Nodal rock energy density. More...

const MaterialProperty< std::vector< Real > > &	_drock_energy_nodal_dvar
	d(nodal rock energy density)/d(PorousFlow variable) More...

const MaterialProperty< std::vector< Real > > *const	_fluid_density
	Nodal fluid density. More...

const MaterialProperty< std::vector< std::vector< Real > > > *const	_dfluid_density_dvar
	d(nodal fluid density)/d(PorousFlow variable) More...

const MaterialProperty< std::vector< Real > > *const	_fluid_saturation_nodal
	Nodal fluid saturation. More...

const MaterialProperty< std::vector< std::vector< Real > > > *const	_dfluid_saturation_nodal_dvar
	d(nodal fluid saturation)/d(PorousFlow variable) More...

const MaterialProperty< std::vector< Real > > *const	_energy_nodal
	Internal energy of the phases, evaluated at the nodes. More...

const MaterialProperty< std::vector< std::vector< Real > > > *const	_denergy_nodal_dvar
	d(internal energy)/d(PorousFlow variable) More...

const MaterialProperty< Real > &	_strain_rate_qp
	Strain rate. More...

const MaterialProperty< std::vector< RealGradient > > &	_dstrain_rate_qp_dvar
	d(strain rate)/d(PorousFlow variable) More...

Detailed Description

Kernel = energy_density * d(volumetric_strain)/dt which is lumped to the nodes.

Definition at line 24 of file PorousFlowHeatVolumetricExpansion.h.

Constructor & Destructor Documentation

◆ PorousFlowHeatVolumetricExpansion()

PorousFlowHeatVolumetricExpansion::PorousFlowHeatVolumetricExpansion ( const InputParameters & parameters )

Definition at line 34 of file PorousFlowHeatVolumetricExpansion.C.

   : TimeKernel(parameters),
     _dictator(getUserObject<PorousFlowDictator>("PorousFlowDictator")),
     _var_is_porflow_var(_dictator.isPorousFlowVariable(_var.number())),
     _num_phases(_dictator.numPhases()),
     _fluid_present(_num_phases > 0),
     _strain_at_nearest_qp(getParam<bool>("strain_at_nearest_qp")),
     _porosity(getMaterialProperty<Real>("PorousFlow_porosity_nodal")),
     _dporosity_dvar(getMaterialProperty<std::vector<Real>>("dPorousFlow_porosity_nodal_dvar")),
     _dporosity_dgradvar(
         getMaterialProperty<std::vector<RealGradient>>("dPorousFlow_porosity_nodal_dgradvar")),
     _nearest_qp(_strain_at_nearest_qp
                     ? &getMaterialProperty<unsigned int>("PorousFlow_nearestqp_nodal")
                     : nullptr),
     _rock_energy_nodal(getMaterialProperty<Real>("PorousFlow_matrix_internal_energy_nodal")),
     _drock_energy_nodal_dvar(
         getMaterialProperty<std::vector<Real>>("dPorousFlow_matrix_internal_energy_nodal_dvar")),
     _fluid_density(_fluid_present ? &getMaterialProperty<std::vector<Real>>(
                                         "PorousFlow_fluid_phase_density_nodal")
                                   : nullptr),
     _dfluid_density_dvar(_fluid_present ? &getMaterialProperty<std::vector<std::vector<Real>>>(
                                               "dPorousFlow_fluid_phase_density_nodal_dvar")
                                         : nullptr),
     _fluid_saturation_nodal(
         _fluid_present ? &getMaterialProperty<std::vector<Real>>("PorousFlow_saturation_nodal")
                        : nullptr),
     _dfluid_saturation_nodal_dvar(_fluid_present
                                       ? &getMaterialProperty<std::vector<std::vector<Real>>>(
                                             "dPorousFlow_saturation_nodal_dvar")
                                       : nullptr),
     _energy_nodal(_fluid_present ? &getMaterialProperty<std::vector<Real>>(
                                        "PorousFlow_fluid_phase_internal_energy_nodal")
                                  : nullptr),
     _denergy_nodal_dvar(_fluid_present ? &getMaterialProperty<std::vector<std::vector<Real>>>(
                                              "dPorousFlow_fluid_phase_internal_energy_nodal_dvar")
                                        : nullptr),
     _strain_rate_qp(getMaterialProperty<Real>("PorousFlow_volumetric_strain_rate_qp")),
     _dstrain_rate_qp_dvar(getMaterialProperty<std::vector<RealGradient>>(
         "dPorousFlow_volumetric_strain_rate_qp_dvar"))
 {
 }

Member Function Documentation

◆ computedEnergyQpJac()

Real PorousFlowHeatVolumetricExpansion::computedEnergyQpJac ( unsigned int jvar )

protected

Derivative of energy part of the residual with respect to the Variable with variable number jvar.

This is used by both computeQpJacobian and computeQpOffDiagJacobian

Parameters

jvar	take the derivative of the energy part of the residual wrt this variable number

Definition at line 117 of file PorousFlowHeatVolumetricExpansion.C.

 {
   if (_dictator.notPorousFlowVariable(jvar))
     return 0.0;
  
   const unsigned int pvar = _dictator.porousFlowVariableNum(jvar);
   const unsigned nearest_qp = (_strain_at_nearest_qp ? (*_nearest_qp)[_i] : _i);
  
   Real denergy = -_dporosity_dgradvar[_i][pvar] * _grad_phi[_j][_i] * _rock_energy_nodal[_i];
   for (unsigned ph = 0; ph < _num_phases; ++ph)
     denergy += (*_fluid_density)[_i][ph] * (*_fluid_saturation_nodal)[_i][ph] *
                (*_energy_nodal)[_i][ph] * _dporosity_dgradvar[_i][pvar] * _grad_phi[_j][nearest_qp];
  
   if (_i != _j)
     return _test[_i][_qp] * denergy * _strain_rate_qp[_qp];
  
   denergy += _drock_energy_nodal_dvar[_i][pvar] * (1.0 - _porosity[_i]);
   denergy -= _rock_energy_nodal[_i] * _dporosity_dvar[_i][pvar];
   for (unsigned ph = 0; ph < _num_phases; ++ph)
   {
     denergy += (*_dfluid_density_dvar)[_i][ph][pvar] * (*_fluid_saturation_nodal)[_i][ph] *
                (*_energy_nodal)[_i][ph] * _porosity[_i];
     denergy += (*_fluid_density)[_i][ph] * (*_dfluid_saturation_nodal_dvar)[_i][ph][pvar] *
                (*_energy_nodal)[_i][ph] * _porosity[_i];
     denergy += (*_fluid_density)[_i][ph] * (*_fluid_saturation_nodal)[_i][ph] *
                (*_denergy_nodal_dvar)[_i][ph][pvar] * _porosity[_i];
     denergy += (*_fluid_density)[_i][ph] * (*_fluid_saturation_nodal)[_i][ph] *
                (*_energy_nodal)[_i][ph] * _dporosity_dvar[_i][pvar];
   }
  
   return _test[_i][_qp] * denergy * _strain_rate_qp[_qp];
 }

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

◆ computedVolQpJac()

Real PorousFlowHeatVolumetricExpansion::computedVolQpJac ( unsigned int jvar )

protected

Derivative of volumetric-strain part of the residual with respect to the Variable with variable number jvar.

This is used by both computeQpJacobian and computeQpOffDiagJacobian

Parameters

jvar	take the derivative of the volumetric-strain part of the residual wrt this variable number

Definition at line 101 of file PorousFlowHeatVolumetricExpansion.C.

 {
   if (_dictator.notPorousFlowVariable(jvar))
     return 0.0;
  
   Real energy = (1.0 - _porosity[_i]) * _rock_energy_nodal[_i];
   for (unsigned ph = 0; ph < _num_phases; ++ph)
     energy += (*_fluid_density)[_i][ph] * (*_fluid_saturation_nodal)[_i][ph] *
               (*_energy_nodal)[_i][ph] * _porosity[_i];
  
   const unsigned int pvar = _dictator.porousFlowVariableNum(jvar);
   Real dvol = _dstrain_rate_qp_dvar[_qp][pvar] * _grad_phi[_j][_qp];
  
   return _test[_i][_qp] * energy * dvol;
 }

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

◆ computeQpJacobian()

Real PorousFlowHeatVolumetricExpansion::computeQpJacobian ( )

overrideprotectedvirtual

Definition at line 89 of file PorousFlowHeatVolumetricExpansion.C.

 {
   return computedEnergyQpJac(_var.number()) + computedVolQpJac(_var.number());
 }

◆ computeQpOffDiagJacobian()

Real PorousFlowHeatVolumetricExpansion::computeQpOffDiagJacobian ( unsigned int jvar )

overrideprotectedvirtual

Definition at line 95 of file PorousFlowHeatVolumetricExpansion.C.

 {
   return computedEnergyQpJac(jvar) + computedVolQpJac(jvar);
 }

◆ computeQpResidual()

Real PorousFlowHeatVolumetricExpansion::computeQpResidual ( )

overrideprotectedvirtual

Definition at line 78 of file PorousFlowHeatVolumetricExpansion.C.

 {
   Real energy = (1.0 - _porosity[_i]) * _rock_energy_nodal[_i];
   for (unsigned ph = 0; ph < _num_phases; ++ph)
     energy += (*_fluid_density)[_i][ph] * (*_fluid_saturation_nodal)[_i][ph] *
               (*_energy_nodal)[_i][ph] * _porosity[_i];
  
   return _test[_i][_qp] * energy * _strain_rate_qp[_qp];
 }

Member Data Documentation

◆ _denergy_nodal_dvar

const MaterialProperty<std::vector<std::vector<Real> > >* const PorousFlowHeatVolumetricExpansion::_denergy_nodal_dvar

protected

d(internal energy)/d(PorousFlow variable)

Definition at line 83 of file PorousFlowHeatVolumetricExpansion.h.

◆ _dfluid_density_dvar

const MaterialProperty<std::vector<std::vector<Real> > >* const PorousFlowHeatVolumetricExpansion::_dfluid_density_dvar

protected

d(nodal fluid density)/d(PorousFlow variable)

Definition at line 71 of file PorousFlowHeatVolumetricExpansion.h.

◆ _dfluid_saturation_nodal_dvar

const MaterialProperty<std::vector<std::vector<Real> > >* const PorousFlowHeatVolumetricExpansion::_dfluid_saturation_nodal_dvar

protected

d(nodal fluid saturation)/d(PorousFlow variable)

Definition at line 77 of file PorousFlowHeatVolumetricExpansion.h.

◆ _dictator

const PorousFlowDictator& PorousFlowHeatVolumetricExpansion::_dictator

protected

PorousFlowDictator UserObject.

Definition at line 35 of file PorousFlowHeatVolumetricExpansion.h.

Referenced by computedEnergyQpJac(), and computedVolQpJac().

◆ _dporosity_dgradvar

const MaterialProperty<std::vector<RealGradient> >& PorousFlowHeatVolumetricExpansion::_dporosity_dgradvar

protected

d(porosity)/d(grad PorousFlow variable)

Definition at line 56 of file PorousFlowHeatVolumetricExpansion.h.

Referenced by computedEnergyQpJac().

◆ _dporosity_dvar

const MaterialProperty<std::vector<Real> >& PorousFlowHeatVolumetricExpansion::_dporosity_dvar

protected

d(porosity)/d(PorousFlow variable)

Definition at line 53 of file PorousFlowHeatVolumetricExpansion.h.

Referenced by computedEnergyQpJac().

◆ _drock_energy_nodal_dvar

const MaterialProperty<std::vector<Real> >& PorousFlowHeatVolumetricExpansion::_drock_energy_nodal_dvar

protected

d(nodal rock energy density)/d(PorousFlow variable)

Definition at line 65 of file PorousFlowHeatVolumetricExpansion.h.

Referenced by computedEnergyQpJac().

◆ _dstrain_rate_qp_dvar

const MaterialProperty<std::vector<RealGradient> >& PorousFlowHeatVolumetricExpansion::_dstrain_rate_qp_dvar

protected

d(strain rate)/d(PorousFlow variable)

Definition at line 89 of file PorousFlowHeatVolumetricExpansion.h.

Referenced by computedVolQpJac().

◆ _energy_nodal

const MaterialProperty<std::vector<Real> >* const PorousFlowHeatVolumetricExpansion::_energy_nodal

protected

Internal energy of the phases, evaluated at the nodes.

Definition at line 80 of file PorousFlowHeatVolumetricExpansion.h.

◆ _fluid_density

const MaterialProperty<std::vector<Real> >* const PorousFlowHeatVolumetricExpansion::_fluid_density

protected

Nodal fluid density.

Definition at line 68 of file PorousFlowHeatVolumetricExpansion.h.

Referenced by computedEnergyQpJac(), computedVolQpJac(), and computeQpResidual().

◆ _fluid_present

const bool PorousFlowHeatVolumetricExpansion::_fluid_present

protected

Whether fluid is present.

Definition at line 44 of file PorousFlowHeatVolumetricExpansion.h.

◆ _fluid_saturation_nodal

const MaterialProperty<std::vector<Real> >* const PorousFlowHeatVolumetricExpansion::_fluid_saturation_nodal

protected

Nodal fluid saturation.

Definition at line 74 of file PorousFlowHeatVolumetricExpansion.h.

◆ _nearest_qp

const MaterialProperty<unsigned int>* const PorousFlowHeatVolumetricExpansion::_nearest_qp

protected

The nearest qp to the node.

Definition at line 59 of file PorousFlowHeatVolumetricExpansion.h.

◆ _num_phases

const unsigned int PorousFlowHeatVolumetricExpansion::_num_phases

protected

Number of fluid phases.

Definition at line 41 of file PorousFlowHeatVolumetricExpansion.h.

Referenced by computedEnergyQpJac(), computedVolQpJac(), and computeQpResidual().

◆ _porosity

const MaterialProperty<Real>& PorousFlowHeatVolumetricExpansion::_porosity

protected

Porosity.

Definition at line 50 of file PorousFlowHeatVolumetricExpansion.h.

Referenced by computedEnergyQpJac(), computedVolQpJac(), and computeQpResidual().

◆ _rock_energy_nodal

const MaterialProperty<Real>& PorousFlowHeatVolumetricExpansion::_rock_energy_nodal

protected

Nodal rock energy density.

Definition at line 62 of file PorousFlowHeatVolumetricExpansion.h.

Referenced by computedEnergyQpJac(), computedVolQpJac(), and computeQpResidual().

◆ _strain_at_nearest_qp

const bool PorousFlowHeatVolumetricExpansion::_strain_at_nearest_qp

protected

Whether the porosity uses the volumetric strain at the closest quadpoint.

Definition at line 47 of file PorousFlowHeatVolumetricExpansion.h.

Referenced by computedEnergyQpJac().

◆ _strain_rate_qp

const MaterialProperty<Real>& PorousFlowHeatVolumetricExpansion::_strain_rate_qp

protected

Strain rate.

Definition at line 86 of file PorousFlowHeatVolumetricExpansion.h.

Referenced by computedEnergyQpJac(), and computeQpResidual().

◆ _var_is_porflow_var

const bool PorousFlowHeatVolumetricExpansion::_var_is_porflow_var

protected

Whether the Variable for this Kernel is a PorousFlow variable according to the Dictator.

Definition at line 38 of file PorousFlowHeatVolumetricExpansion.h.

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

porous_flow/include/kernels/PorousFlowHeatVolumetricExpansion.h
porous_flow/src/kernels/PorousFlowHeatVolumetricExpansion.C

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ PorousFlowHeatVolumetricExpansion()

Member Function Documentation

◆ computedEnergyQpJac()

◆ computedVolQpJac()

◆ computeQpJacobian()

◆ computeQpOffDiagJacobian()

◆ computeQpResidual()

Member Data Documentation

◆ _denergy_nodal_dvar

◆ _dfluid_density_dvar

◆ _dfluid_saturation_nodal_dvar

◆ _dictator

◆ _dporosity_dgradvar

◆ _dporosity_dvar

◆ _drock_energy_nodal_dvar

◆ _dstrain_rate_qp_dvar

◆ _energy_nodal

◆ _fluid_density

◆ _fluid_present

◆ _fluid_saturation_nodal

◆ _nearest_qp

◆ _num_phases

◆ _porosity

◆ _rock_energy_nodal

◆ _strain_at_nearest_qp

◆ _strain_rate_qp

◆ _var_is_porflow_var