SCMHTCClosureBase.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 "SCMClosureBase.h"
#include "SubChannel1PhaseProblem.h"

#include <memory>

struct NusseltPreInfo
{
  Real Re;                 // Reynolds number
  Real Pr;                 // Prandtl number
  Real poD;                // Pitch over diameter ratio
  Real ReL;                // Laminar Reynolds number limit
  Real ReT;                // Turbulent Reynolds number limit
  Real laminar_Nu;         // Laminar Nusselt number
  EChannelType subch_type; // Subchannel type (corner, edge, center)
};

class SCMHTCClosureBase : public SCMClosureBase
{
public:
  static InputParameters validParams();

  SCMHTCClosureBase(const InputParameters & parameters);

  typedef SubChannel1PhaseProblem::FrictionStruct FrictionStruct;
  typedef SubChannel1PhaseProblem::NusseltStruct NusseltStruct;

  virtual Real computeNusseltNumber(const FrictionStruct & friction_info,
                                    const NusseltStruct & nusselt_info) const = 0;

  NusseltPreInfo computeNusseltNumberPreInfo(const NusseltStruct & nusselt_info) const;

  Real computeHTC(const FrictionStruct & friction_info,
                  const NusseltStruct & nusselt_info,
                  const Real conduction_k) const;

protected:
  Real blendTurbulentNusseltNumber(const NusseltPreInfo & nusselt_info,
                                   const Real turbulent_nusselt) const;

  std::unique_ptr<SolutionHandle> _Dpin_soln;
};