PorousFlowHeatEnergy Class Reference

Postprocessor produces the sum of heat energy of the porous skeleton and/or fluid components in a region. More...

#include <PorousFlowHeatEnergy.h>

Inheritance diagram for PorousFlowHeatEnergy:

Public Member Functions
	PorousFlowHeatEnergy (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeIntegral () override
virtual Real	computeQpIntegral () override

Protected Attributes
const PorousFlowDictator &	_dictator
	PorousFlowDictator UserObject. More...
const unsigned int	_num_phases
	Number of fluid phases. More...
const bool	_fluid_present
	Whether fluid is present. More...
const bool	_include_porous_skeleton
	Whether to include the heat energy of the porous skeleton in the calculations. More...
std::vector< unsigned int >	_phase_index
	The phase indices that this Postprocessor is restricted to. More...
const MaterialProperty< Real > &	_porosity
	Porosity. More...
const MaterialProperty< Real > &	_rock_energy_nodal
	Nodal rock energy density. More...
const MaterialProperty< std::vector< Real > > *const	_fluid_density
	Nodal fluid density. More...
const MaterialProperty< std::vector< Real > > *const	_fluid_saturation_nodal
	Nodal fluid saturation. More...
const MaterialProperty< std::vector< Real > > *const	_energy_nodal
	Internal energy of the phases, evaluated at the nodes. More...
MooseVariable *const	_var
	The variable for the corresponding PorousFlowEnergyTimeDerivative Kernel: this provides test functions. More...

Detailed Description

Postprocessor produces the sum of heat energy of the porous skeleton and/or fluid components in a region.

Definition at line 24 of file PorousFlowHeatEnergy.h.

Constructor & Destructor Documentation

PorousFlowHeatEnergy()

PorousFlowHeatEnergy::PorousFlowHeatEnergy

(

const InputParameters &

parameters

)

Definition at line 43 of file PorousFlowHeatEnergy.C.

  : ElementIntegralPostprocessor(parameters),
    _dictator(getUserObject<PorousFlowDictator>("PorousFlowDictator")),
    _num_phases(_dictator.numPhases()),
    _fluid_present(_num_phases > 0),
    _include_porous_skeleton(getParam<bool>("include_porous_skeleton")),
    _phase_index(getParam<std::vector<unsigned int>>("phase")),
    _porosity(getMaterialProperty<Real>("PorousFlow_porosity_nodal")),
    _rock_energy_nodal(getMaterialProperty<Real>("PorousFlow_matrix_internal_energy_nodal")),
    _fluid_density(_fluid_present ? &getMaterialProperty<std::vector<Real>>(
                                        "PorousFlow_fluid_phase_density_nodal")
                                  : nullptr),
    _fluid_saturation_nodal(
        _fluid_present ? &getMaterialProperty<std::vector<Real>>("PorousFlow_saturation_nodal")
                       : nullptr),
    _energy_nodal(_fluid_present ? &getMaterialProperty<std::vector<Real>>(
                                       "PorousFlow_fluid_phase_internal_energy_nodal")
                                 : nullptr),
    _var(getParam<unsigned>("kernel_variable_number") < _dictator.numVariables()
             ? &_fe_problem.getStandardVariable(
                   _tid,
                   _dictator
                       .getCoupledStandardMooseVars()[getParam<unsigned>("kernel_variable_number")]
                       ->name())
             : nullptr)
{
  if (!_phase_index.empty())
  {
    // Check that the phase indices entered are not greater than the number of phases
    const unsigned int max_phase_num = *std::max_element(_phase_index.begin(), _phase_index.end());
    if (max_phase_num > _num_phases - 1)
      paramError("phase",
                 "The Dictator proclaims that the phase index ",
                 max_phase_num,
                 " is greater than the largest phase index possible, which is ",
                 _num_phases - 1);
  }
 
  // Check that kernel_variable_number is OK
  if (getParam<unsigned>("kernel_variable_number") >= _dictator.numVariables())
    paramError("kernel_variable_number",
               "The Dictator pronounces that the number of PorousFlow variables is ",
               _dictator.numVariables(),
               ", however you have used ",
               getParam<unsigned>("kernel_variable_number"),
               ". This is an error");
 
  // Now that we know kernel_variable_number is OK, _var must be OK,
  // so ensure that reinit is called on _var:
  addMooseVariableDependency(_var);
}

Member Function Documentation

computeIntegral()

Real PorousFlowHeatEnergy::computeIntegral ( )

overrideprotectedvirtual

Definition at line 96 of file PorousFlowHeatEnergy.C.

{
  Real sum = 0;
 
  // The use of _test in the loops below mean that the
  // integral is exactly the same as the one computed
  // by the PorousFlowMassTimeDerivative Kernel.  Because that
  // Kernel is lumped, this Postprocessor also needs to
  // be lumped.  Hence the use of the "nodal" Material
  // Properties
  const VariableTestValue & test = _var->phi();
 
  for (unsigned node = 0; node < test.size(); ++node)
  {
    Real nodal_volume = 0.0;
    for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
      nodal_volume += _JxW[_qp] * _coord[_qp] * test[node][_qp];
 
    Real energy = 0.0;
    if (_include_porous_skeleton)
      energy += (1.0 - _porosity[node]) * _rock_energy_nodal[node];
 
    for (auto ph : _phase_index)
      energy += (*_fluid_density)[node][ph] * (*_fluid_saturation_nodal)[node][ph] *
                (*_energy_nodal)[node][ph] * _porosity[node];
 
    sum += nodal_volume * energy;
  }
 
  return sum;
}

computeQpIntegral()

Real PorousFlowHeatEnergy::computeQpIntegral ( )

overrideprotectedvirtual

Definition at line 129 of file PorousFlowHeatEnergy.C.

{
  return 0.0;
}

Member Data Documentation

_dictator

const PorousFlowDictator& PorousFlowHeatEnergy::_dictator

protected

PorousFlowDictator UserObject.

Definition at line 34 of file PorousFlowHeatEnergy.h.

Referenced by PorousFlowHeatEnergy().

_energy_nodal

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

protected

Internal energy of the phases, evaluated at the nodes.

Definition at line 61 of file PorousFlowHeatEnergy.h.

_fluid_density

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

protected

Nodal fluid density.

Definition at line 55 of file PorousFlowHeatEnergy.h.

_fluid_present

const bool PorousFlowHeatEnergy::_fluid_present

protected

Whether fluid is present.

Definition at line 40 of file PorousFlowHeatEnergy.h.

_fluid_saturation_nodal

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

protected

Nodal fluid saturation.

Definition at line 58 of file PorousFlowHeatEnergy.h.

_include_porous_skeleton

const bool PorousFlowHeatEnergy::_include_porous_skeleton

protected

Whether to include the heat energy of the porous skeleton in the calculations.

Definition at line 43 of file PorousFlowHeatEnergy.h.

Referenced by computeIntegral().

_num_phases

const unsigned int PorousFlowHeatEnergy::_num_phases

protected

Number of fluid phases.

Definition at line 37 of file PorousFlowHeatEnergy.h.

Referenced by PorousFlowHeatEnergy().

_phase_index

std::vector<unsigned int> PorousFlowHeatEnergy::_phase_index

protected

The phase indices that this Postprocessor is restricted to.

Definition at line 46 of file PorousFlowHeatEnergy.h.

Referenced by computeIntegral(), and PorousFlowHeatEnergy().

_porosity

const MaterialProperty<Real>& PorousFlowHeatEnergy::_porosity

protected

Porosity.

Definition at line 49 of file PorousFlowHeatEnergy.h.

Referenced by computeIntegral().

_rock_energy_nodal

const MaterialProperty<Real>& PorousFlowHeatEnergy::_rock_energy_nodal

protected

Nodal rock energy density.

Definition at line 52 of file PorousFlowHeatEnergy.h.

Referenced by computeIntegral().

_var

MooseVariable* const PorousFlowHeatEnergy::_var

protected

The variable for the corresponding PorousFlowEnergyTimeDerivative Kernel: this provides test functions.

Definition at line 64 of file PorousFlowHeatEnergy.h.

Referenced by computeIntegral(), and PorousFlowHeatEnergy().

---

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

• porous_flow/include/postprocessors/PorousFlowHeatEnergy.h
• porous_flow/src/postprocessors/PorousFlowHeatEnergy.C