CHTHandler.h

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

#pragma once

#include "MooseObject.h"
#include "FaceCenteredMapFunctor.h"
#include "SystemBase.h"
#include "NS.h"

class LinearFVFluxKernel;
class LinearFVBoundaryCondition;

namespace NS
{
namespace FV
{

class CHTHandler : public MooseObject
{
public:
  CHTHandler(const InputParameters & parameters);

  static InputParameters validParams();

  void linkEnergySystems(SystemBase * solid_energy_system,
                         SystemBase * fluid_energy_system,
                         std::vector<SystemBase *> pm_radiation_systems);

  void setupConjugateHeatTransferContainers();

  void deduceCHTBoundaryCoupling();

  void updateCHTBoundaryCouplingFields(const NS::CHTSide side);

  void initializeCHTCouplingFields();

  bool converged() const;

  void resetCHTConvergence();

  void incrementCHTIterators();

  void sumIntegratedFluxes();

  void printIntegratedFluxes() const;

  void resetIntegratedFluxes();

  virtual bool enabled() const override final;

protected:
  FEProblemBase & _problem;

  MooseMesh & _mesh;

  SystemBase * _energy_system;

  SystemBase * _solid_energy_system;

  std::vector<SystemBase *> _pm_radiation_systems;

  std::vector<BoundaryName> _cht_boundary_names;

  std::vector<BoundaryID> _cht_boundary_ids;

  const unsigned int _max_cht_fpi;

  const Real _cht_heat_flux_tolerance;

  std::vector<std::vector<Real>> _cht_flux_relaxation_factor;

  std::vector<std::vector<Real>> _cht_temperature_relaxation_factor;

  std::vector<unsigned int> _cht_system_numbers;

  std::vector<unsigned int> _cht_pm_radiation_system_numbers;

  std::vector<std::vector<const FaceInfo *>> _cht_face_info;

  std::vector<LinearFVFluxKernel *> _cht_conduction_kernels;

  std::vector<LinearFVFluxKernel *> _cht_pm_radiation_kernels;

  std::vector<std::vector<LinearFVBoundaryCondition *>> _cht_boundary_conditions;

  std::vector<std::vector<LinearFVBoundaryCondition *>> _cht_pm_radiation_boundary_conditions;

  std::vector<std::vector<FaceCenteredMapFunctor<Real, std::unordered_map<dof_id_type, Real>>>>
      _boundary_heat_flux;

  std::vector<std::vector<Real>> _integrated_boundary_heat_flux;

  std::vector<std::vector<FaceCenteredMapFunctor<Real, std::unordered_map<dof_id_type, Real>>>>
      _boundary_temperature;

private:
  unsigned int _fpi_it;
};

inline bool
CHTHandler::enabled() const
{
  return !_cht_boundary_names.empty();
}

inline void
CHTHandler::resetCHTConvergence()
{
  _fpi_it = 0;
}

inline void
CHTHandler::incrementCHTIterators()
{
  _fpi_it++;
}

} // End FV namespace
} // End Moose namespace