programmer's documentation
 All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Macros Groups Pages
Functions/Subroutines
hturbp.f90 File Reference

Functions/Subroutines

subroutine hturbp (prl, prt, ckarm, yplus, dplus, htur, yplim)
 Compute the correction of the exchange coefficient between the fluid and the wall for a turbulent flow. More...
 

Function/Subroutine Documentation

subroutine hturbp ( double precision  prl,
double precision  prt,
double precision  ckarm,
double precision  yplus,
double precision  dplus,
double precision  htur,
double precision  yplim 
)

Compute the correction of the exchange coefficient between the fluid and the wall for a turbulent flow.

This is function of the dimensionless distance to the wall $ y^+ = \dfrac{\centip \centf u_\star}{\nu}$.

Then the return coefficient reads:

\[ h_{tur} = Pr \dfrac{y^+}{T^+} \]

This coefficient is computed thanks to a similarity model between dynamic viscous sub-layer and themal sub-layer.

$ T^+ $ is computed as follows:

  • For a laminar Prandtl number smaller than 0.1 (such as liquid metals), the standard model with two sub-layers (Prandtl-Taylor) is used.
  • For a laminar Prandtl number larger than 0.1 (such as liquids and gaz), a model with three sub-layers (Arpaci-Larsen) is used.

The final exchange coefficient is:

\[ h = \dfrac{K}{\centip \centf} h_{tur} \]

Parameters
[in]prllaminar Prandtl number
[in]prtturbulent Prandtl number
[in]ckarmVon Karman constant
[in]yplusdimensionless distance to the wall
[in]dplusdimensionless distance for scalable wall functions
[out]hturcorrected exchange coefficient
[out]yplimvalue of the limit for $ y^+ $