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 turbulentReynoldsNumber(const NusseltPreInfo & nusselt_info) const;
 
  Real blendTurbulentNusseltNumber(const NusseltPreInfo & nusselt_info,
                                   const Real turbulent_nusselt) const;
 
  std::unique_ptr<SolutionHandle> _Dpin_soln;
};