Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : #include "SCMHTCKazimiCarelli.h"
      11             : 
      12             : registerMooseObject("SubChannelApp", SCMHTCKazimiCarelli);
      13             : 
      14             : InputParameters
      15          79 : SCMHTCKazimiCarelli::validParams()
      16             : {
      17          79 :   InputParameters params = SCMHTCClosureBase::validParams();
      18          79 :   params.addClassDescription(
      19             :       "Class that computes the convective heat transfer coefficient using the "
      20             :       "Kazimi-Carelli correlation. Only use for fuel-pins.");
      21          79 :   return params;
      22           0 : }
      23             : 
      24          44 : SCMHTCKazimiCarelli::SCMHTCKazimiCarelli(const InputParameters & parameters)
      25          44 :   : SCMHTCClosureBase(parameters)
      26             : {
      27             :   // Check that Kazimi-Carelli is not used for the duct (not supported yet)
      28          44 :   if (const auto * duct_uo = _scm_problem.getDuctHTCClosure(); duct_uo && duct_uo == this)
      29           0 :     mooseError("'Kazimi-Carelli' is not yet supported for the 'duct_htc_correlation'.");
      30          44 : }
      31             : 
      32             : Real
      33       57456 : SCMHTCKazimiCarelli::computeNusseltNumber(const FrictionStruct & /*friction_args*/,
      34             :                                           const NusseltStruct & nusselt_args) const
      35             : {
      36       57456 :   const auto pre = computeNusseltNumberPreInfo(nusselt_args);
      37       57456 :   const auto Pe = pre.Re * pre.Pr;
      38             : 
      39       57456 :   if (Pe < 10 || Pe > 5000)
      40           0 :     flagSolutionWarning("Peclet number (Pe) out of range for the Kazimi-Carelli correlation.");
      41             : 
      42       57456 :   if (pre.poD < 1.1 || pre.poD > 1.4)
      43           0 :     flagSolutionWarning(
      44             :         "Pitch over pin diameter ratio out of range for the Kazimi-Carelli correlation.");
      45             : 
      46       57456 :   const auto NuT = 4.0 + 0.33 * std::pow(pre.poD, 3.8) * std::pow((Pe / 1e2), 0.86) +
      47       57456 :                    0.16 * std::pow(pre.poD, 5);
      48       57456 :   return blendTurbulentNusseltNumber(pre, NuT);
      49             : }