LinearFVEnthalpyVolumetricHeatTransfer.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html
 
#include "LinearFVEnthalpyVolumetricHeatTransfer.h"
#include "NS.h"
 
registerMooseObject("NavierStokesApp", LinearFVEnthalpyVolumetricHeatTransfer);
 
InputParameters
LinearFVEnthalpyVolumetricHeatTransfer::validParams()
{
  InputParameters params = LinearFVVolumetricHeatTransfer::validParams();
  params.addClassDescription(
      "Represents a volumetric heat transfer term between two porous-medium energy equations when "
      "the solved variable is specific enthalpy and the current-side temperature is linearized "
      "with dT/dh ~= 1/cp.");
  params.addRequiredParam<MooseFunctorName>(
      NS::cp,
      "Specific heat used to linearize the current-side temperature from the enthalpy variable.");
  return params;
}
 
LinearFVEnthalpyVolumetricHeatTransfer::LinearFVEnthalpyVolumetricHeatTransfer(
    const InputParameters & params)
  : LinearFVVolumetricHeatTransfer(params), _cp(getFunctor<Real>(NS::cp))
{
}
 
Real
LinearFVEnthalpyVolumetricHeatTransfer::positiveCp(const Moose::ElemArg & elem,
                                                   const Moose::StateArg & state) const
{
  const auto cp = _cp(elem, state);
  if (cp <= 0.0)
    mooseError("cp must stay positive when solving for enthalpy, but cp = ",
               cp,
               " in element ",
               elem.elem->id(),
               " (subdomain ",
               elem.elem->subdomain_id(),
               ", vertex average ",
               elem.elem->vertex_average(),
               ").");
  return cp;
}
 
Real
LinearFVEnthalpyVolumetricHeatTransfer::computeMatrixContribution()
{
  const auto elem = makeElemArg(_current_elem_info->elem());
  const auto state = determineState();
  const auto cp = positiveCp(elem, state);
 
  return _h_solid_fluid(elem, state) * _current_elem_volume / cp;
}
 
Real
LinearFVEnthalpyVolumetricHeatTransfer::computeRightHandSideContribution()
{
  const auto elem = makeElemArg(_current_elem_info->elem());
  const auto state = determineState();
  const auto cp = positiveCp(elem, state);
 
  const auto & current_temperature = _is_solid ? _temp_solid : _temp_fluid;
  const auto & other_temperature = _is_solid ? _temp_fluid : _temp_solid;
  const auto h_current = _var.getElemValue(*_current_elem_info, state);
  const auto linearized_current_temperature = h_current / cp;
 
  // linearized temperature is added to the matrix and removed here
  return _h_solid_fluid(elem, state) *
         (other_temperature.getElemValue(*_current_elem_info, state) -
          current_temperature.getElemValue(*_current_elem_info, state) +
          linearized_current_temperature) *
         _current_elem_volume;
}