Kernel = (heat_energy - heat_energy_old)/dt It is lumped to the nodes. More...

#include <PorousFlowEnergyTimeDerivative.h>

Inheritance diagram for PorousFlowEnergyTimeDerivative:

Public Member Functions
	PorousFlowEnergyTimeDerivative (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpResidual () override

virtual Real	computeQpJacobian () override

virtual Real	computeQpOffDiagJacobian (unsigned int jvar) override

Real	computeQpJac (unsigned int pvar) const
	Derivative of residual with respect to PorousFlow variable number pvar This is used by both computeQpJacobian and computeQpOffDiagJacobian. 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 _num_phases > 0 (ie. there is a fluid 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 at the nodes, but it can depend on grad(variables) which are actually evaluated at the qps. More...

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

const MaterialProperty< std::vector< Real > > &	_dporosity_dvar
	d(porosity)/d(PorousFlow variable) - these derivatives will be wrt variables at the nodes More...

const MaterialProperty< std::vector< RealGradient > > &	_dporosity_dgradvar
	d(porosity)/d(grad PorousFlow variable) - remember these derivatives will be wrt grad(vars) at qps 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< Real > &	_rock_energy_nodal_old
	Old value of 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< Real > > *const	_fluid_density_old
	Old value of 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< Real > > *const	_fluid_saturation_nodal_old
	Old value of 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< Real > > *const	_energy_nodal_old
	Old value of 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...

Detailed Description

Kernel = (heat_energy - heat_energy_old)/dt It is lumped to the nodes.

Definition at line 24 of file PorousFlowEnergyTimeDerivative.h.

Constructor & Destructor Documentation

◆ PorousFlowEnergyTimeDerivative()

PorousFlowEnergyTimeDerivative::PorousFlowEnergyTimeDerivative ( const InputParameters & parameters )

Definition at line 34 of file PorousFlowEnergyTimeDerivative.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")),
     _porosity_old(getMaterialPropertyOld<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")),
     _rock_energy_nodal_old(getMaterialPropertyOld<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),
     _fluid_density_old(_fluid_present ? &getMaterialPropertyOld<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),
     _fluid_saturation_nodal_old(
         _fluid_present ? &getMaterialPropertyOld<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),
     _energy_nodal_old(_fluid_present ? &getMaterialPropertyOld<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)
 {
 }

Member Function Documentation

◆ computeQpJac()

Real PorousFlowEnergyTimeDerivative::computeQpJac ( unsigned int pvar ) const

protected

Derivative of residual with respect to PorousFlow variable number pvar This is used by both computeQpJacobian and computeQpOffDiagJacobian.

Parameters

pvar	take the derivative of the residual wrt this PorousFlow variable

Definition at line 119 of file PorousFlowEnergyTimeDerivative.C.

 {
   const unsigned nearest_qp = (_strain_at_nearest_qp ? (*_nearest_qp)[_i] : _i);
  
   // porosity is dependent on variables that are lumped to the nodes,
   // but it can depend on the gradient
   // of variables, which are NOT lumped to the nodes, hence:
   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 / _dt;
  
   // As the fluid energy is lumped to the nodes, only non-zero terms are for _i==_j
   denergy += -_dporosity_dvar[_i][pvar] * _rock_energy_nodal[_i];
   denergy += (1.0 - _porosity[_i]) * _drock_energy_nodal_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 / _dt;
 }

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

◆ computeQpJacobian()

Real PorousFlowEnergyTimeDerivative::computeQpJacobian ( )

overrideprotectedvirtual

Definition at line 101 of file PorousFlowEnergyTimeDerivative.C.

 {
   // If the variable is not a PorousFlow variable (very unusual), the diag Jacobian terms are 0
   if (!_var_is_porflow_var)
     return 0.0;
   return computeQpJac(_dictator.porousFlowVariableNum(_var.number()));
 }

◆ computeQpOffDiagJacobian()

Real PorousFlowEnergyTimeDerivative::computeQpOffDiagJacobian ( unsigned int jvar )

overrideprotectedvirtual

Definition at line 110 of file PorousFlowEnergyTimeDerivative.C.

 {
   // If the variable is not a PorousFlow variable, the OffDiag Jacobian terms are 0
   if (_dictator.notPorousFlowVariable(jvar))
     return 0.0;
   return computeQpJac(_dictator.porousFlowVariableNum(jvar));
 }

◆ computeQpResidual()

Real PorousFlowEnergyTimeDerivative::computeQpResidual ( )

overrideprotectedvirtual

Definition at line 85 of file PorousFlowEnergyTimeDerivative.C.

 {
   Real energy = (1.0 - _porosity[_i]) * _rock_energy_nodal[_i];
   Real energy_old = (1.0 - _porosity_old[_i]) * _rock_energy_nodal_old[_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];
     energy_old += (*_fluid_density_old)[_i][ph] * (*_fluid_saturation_nodal_old)[_i][ph] *
                   (*_energy_nodal_old)[_i][ph] * _porosity_old[_i];
   }
  
   return _test[_i][_qp] * (energy - energy_old) / _dt;
 }

Member Data Documentation

◆ _denergy_nodal_dvar

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

protected

d(internal energy)/d(PorousFlow variable)

Definition at line 98 of file PorousFlowEnergyTimeDerivative.h.

◆ _dfluid_density_dvar

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

protected

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

Definition at line 80 of file PorousFlowEnergyTimeDerivative.h.

◆ _dfluid_saturation_nodal_dvar

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

protected

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

Definition at line 89 of file PorousFlowEnergyTimeDerivative.h.

◆ _dictator

const PorousFlowDictator& PorousFlowEnergyTimeDerivative::_dictator

protected

PorousFlowDictator UserObject.

Definition at line 35 of file PorousFlowEnergyTimeDerivative.h.

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

◆ _dporosity_dgradvar

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

protected

d(porosity)/d(grad PorousFlow variable) - remember these derivatives will be wrt grad(vars) at qps

Definition at line 59 of file PorousFlowEnergyTimeDerivative.h.

Referenced by computeQpJac().

◆ _dporosity_dvar

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

protected

d(porosity)/d(PorousFlow variable) - these derivatives will be wrt variables at the nodes

Definition at line 56 of file PorousFlowEnergyTimeDerivative.h.

Referenced by computeQpJac().

◆ _drock_energy_nodal_dvar

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

protected

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

Definition at line 71 of file PorousFlowEnergyTimeDerivative.h.

Referenced by computeQpJac().

◆ _energy_nodal

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

protected

Internal energy of the phases, evaluated at the nodes.

Definition at line 92 of file PorousFlowEnergyTimeDerivative.h.

◆ _energy_nodal_old

const MaterialProperty<std::vector<Real> >* const PorousFlowEnergyTimeDerivative::_energy_nodal_old

protected

Old value of internal energy of the phases, evaluated at the nodes.

Definition at line 95 of file PorousFlowEnergyTimeDerivative.h.

◆ _fluid_density

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

protected

Nodal fluid density.

Definition at line 74 of file PorousFlowEnergyTimeDerivative.h.

Referenced by computeQpJac().

◆ _fluid_density_old

const MaterialProperty<std::vector<Real> >* const PorousFlowEnergyTimeDerivative::_fluid_density_old

protected

Old value of nodal fluid density.

Definition at line 77 of file PorousFlowEnergyTimeDerivative.h.

◆ _fluid_present

const bool PorousFlowEnergyTimeDerivative::_fluid_present

protected

Whether _num_phases > 0 (ie. there is a fluid present)

Definition at line 44 of file PorousFlowEnergyTimeDerivative.h.

◆ _fluid_saturation_nodal

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

protected

Nodal fluid saturation.

Definition at line 83 of file PorousFlowEnergyTimeDerivative.h.

◆ _fluid_saturation_nodal_old

const MaterialProperty<std::vector<Real> >* const PorousFlowEnergyTimeDerivative::_fluid_saturation_nodal_old

protected

Old value of fluid saturation.

Definition at line 86 of file PorousFlowEnergyTimeDerivative.h.

◆ _nearest_qp

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

protected

The nearest qp to the node.

Definition at line 62 of file PorousFlowEnergyTimeDerivative.h.

◆ _num_phases

const unsigned int PorousFlowEnergyTimeDerivative::_num_phases

protected

Number of fluid phases.

Definition at line 41 of file PorousFlowEnergyTimeDerivative.h.

Referenced by computeQpJac(), and computeQpResidual().

◆ _porosity

const MaterialProperty<Real>& PorousFlowEnergyTimeDerivative::_porosity

protected

Porosity at the nodes, but it can depend on grad(variables) which are actually evaluated at the qps.

Definition at line 50 of file PorousFlowEnergyTimeDerivative.h.

Referenced by computeQpJac(), and computeQpResidual().

◆ _porosity_old

const MaterialProperty<Real>& PorousFlowEnergyTimeDerivative::_porosity_old

protected

Old value of porosity.

Definition at line 53 of file PorousFlowEnergyTimeDerivative.h.

Referenced by computeQpResidual().

◆ _rock_energy_nodal

const MaterialProperty<Real>& PorousFlowEnergyTimeDerivative::_rock_energy_nodal

protected

Nodal rock energy density.

Definition at line 65 of file PorousFlowEnergyTimeDerivative.h.

Referenced by computeQpJac(), and computeQpResidual().

◆ _rock_energy_nodal_old

const MaterialProperty<Real>& PorousFlowEnergyTimeDerivative::_rock_energy_nodal_old

protected

Old value of nodal rock energy density.

Definition at line 68 of file PorousFlowEnergyTimeDerivative.h.

Referenced by computeQpResidual().

◆ _strain_at_nearest_qp

const bool PorousFlowEnergyTimeDerivative::_strain_at_nearest_qp

protected

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

Definition at line 47 of file PorousFlowEnergyTimeDerivative.h.

Referenced by computeQpJac().

◆ _var_is_porflow_var

const bool PorousFlowEnergyTimeDerivative::_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 PorousFlowEnergyTimeDerivative.h.

Referenced by computeQpJacobian().

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

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

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ PorousFlowEnergyTimeDerivative()

Member Function Documentation

◆ computeQpJac()

◆ 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

◆ _energy_nodal

◆ _energy_nodal_old

◆ _fluid_density

◆ _fluid_density_old

◆ _fluid_present

◆ _fluid_saturation_nodal

◆ _fluid_saturation_nodal_old

◆ _nearest_qp

◆ _num_phases

◆ _porosity

◆ _porosity_old

◆ _rock_energy_nodal

◆ _rock_energy_nodal_old

◆ _strain_at_nearest_qp

◆ _var_is_porflow_var