HeatTransferUtils.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 "MooseTypes.h"

namespace HeatTransferUtils
{

template <typename T1, typename T2, typename T3, typename T4>
auto
reynolds(const T1 & rho, const T2 & vel, const T3 & L, const T4 & mu)
{
  return rho * std::fabs(vel) * L / mu;
}

template <typename T1, typename T2, typename T3>
auto
prandtl(const T1 & cp, const T2 & mu, const T3 & k)
{
  return cp * mu / k;
}

template <typename T1, typename T2, typename T3, typename T4, typename T5, typename T6>
auto
grashof(const T1 & beta,
        const T2 & T_s,
        const T3 & T_bulk,
        const T4 & L,
        const T5 & rho,
        const T6 & mu,
        Real gravity_magnitude)
{
  return gravity_magnitude * beta * std::abs(T_s - T_bulk) * std::pow(L, 3) * (rho * rho) /
         (mu * mu);
}

template <typename T1, typename T2, typename T3, typename T4>
auto
laplace(const T1 & sigma, const T2 & rho, const T3 & L, const T4 & mu)
{
  return sigma * rho * L / (mu * mu);
}

template <typename T1, typename T2, typename T3>
auto
thermalDiffusivity(const T1 & k, const T2 & rho, const T3 & cp)
{
  return k / (rho * cp);
}

template <typename T1, typename T2, typename T3>
auto
peclet(const T1 & vel, const T2 & L, const T3 & diffusivity)
{
  return std::fabs(vel) * L / diffusivity;
}

} // namespace HeatTransferUtils