[Fluent Inc. Logo] return to home search
next up previous contents index

12.10.2 Standard Wall Functions

The standard wall functions in FLUENT are based on the proposal of Launder and Spalding [ 197], and have been most widely used for industrial flows. They are provided as a default option in FLUENT.

Momentum

The law-of-the-wall for mean velocity yields

 U^* = \frac{1}{\kappa} \, \ln (E y^*) (12.10-1)

where
 U^* \equiv \frac{U_P C_{\mu}^{1/4} k_P^{1/2}}{\tau_w/\rho} (12.10-2)


 y^* \equiv \frac{\rho C_{\mu}^{1/4} k_P^{1/2} y_P}{\mu} (12.10-3)


and $\kappa$ = von Kármán constant (= 0.4187)
  $E$ = empirical constant (= 9.793)
  $U_P$ = mean velocity of the fluid at point $P$
  $k_P$ = turbulence kinetic energy at point $P$
  $y_P$ = distance from point $P$ to the wall
  $\mu$ = dynamic viscosity of the fluid

The logarithmic law for mean velocity is known to be valid for $30 < y^* < 300$. In FLUENT, the log-law is employed when $y^* > 11.225$.

When the mesh is such that $y^* < 11.225$ at the wall-adjacent cells, FLUENT applies the laminar stress-strain relationship that can be written as
 U^* = y^* (12.10-4)

It should be noted that, in FLUENT, the laws-of-the-wall for mean velocity and temperature are based on the wall unit, $y^*$, rather than $y^+$ ( $\equiv \rho u_{\tau} y/\mu$). These quantities are approximately equal in equilibrium turbulent boundary layers.

Energy

Reynolds' analogy between momentum and energy transport gives a similar logarithmic law for mean temperature. As in the law-of-the-wall for mean velocity, the law-of-the-wall for temperature employed in FLUENT comprises the following two different laws:

The thickness of the thermal conduction layer is, in general, different from the thickness of the (momentum) viscous sublayer, and changes from fluid to fluid. For example, the thickness of the thermal sublayer for a high-Prandtl-number fluid (e.g., oil) is much less than its momentum sublayer thickness. For fluids of low Prandtl numbers (e.g., liquid metal), on the contrary, it is much larger than the momentum sublayer thickness.

In highly compressible flows, the temperature distribution in the near-wall region can be significantly different from that of low subsonic flows, due to the heating by viscous dissipation. In FLUENT, the temperature wall functions include the contribution from the viscous heating [ 382].

The law-of-the-wall implemented in FLUENT has the following composite form:
 T^* \equiv \frac{\left(T_w - T_P\right) \rho c_p C_{\mu} ^... ...ht\} & \; \; \; (y^* > y^*_T) \\ \end{array} \right . (12.10-5)

where $P$ is computed by using the formula given by Jayatilleke [ 165]:
 P = 9.24 \left[ \left(\frac{\rm Pr}{{\rm Pr}_t}\right)^{3/4... ...ight] \left[1 + 0.28 e^{-0.007 {\rm Pr}/{\rm Pr}_t} \right] (12.10-6)

and

  $k_P$ = turbulent kinetic energy at point P
  $\rho$ = density of fluid
  $c_p$ = specific heat of fluid
  $\dot{q}$ = wall heat flux
  $T_P$ = temperature at the cell adjacent to wall
  $T_w$ = temperature at the wall
  ${\rm Pr}$ = molecular Prandtl number ( $\mu c_p/k_f$)
  ${\rm Pr}_t$ = turbulent Prandtl number (0.85 at the wall)
  $A$ = Van Driest constant (= 26)
  $U_c$ = mean velocity magnitude at $y^{*}=y_T^{*}$

Note that, for the pressure-based solver, the terms

\frac{1}{2} \rho {\rm Pr} \frac{C_\mu^{1/4}k_P^{1/2}}{\dot{q}} U_P^2

and

\frac{1}{2} \rho \frac{C_\mu^{1/4}k_P^{1/2}}{\dot{q}} \left\{{\rm Pr}_t U_P^2 + ({\rm Pr} - {\rm Pr}_t)U_c^2 \right\}

will be included in Equation  12.10-5 only for compressible flow calculations.

The non-dimensional thermal sublayer thickness, $y^*_T$, in Equation  12.10-5 is computed as the $y^*$ value at which the linear law and the logarithmic law intersect, given the molecular Prandtl number of the fluid being modeled.

The procedure of applying the law-of-the-wall for temperature is as follows. Once the physical properties of the fluid being modeled are specified, its molecular Prandtl number is computed. Then, given the molecular Prandtl number, the thermal sublayer thickness, $y^*_T$, is computed from the intersection of the linear and logarithmic profiles, and stored.

During the iteration, depending on the $y^*$ value at the near-wall cell, either the linear or the logarithmic profile in Equation  12.10-5 is applied to compute the wall temperature $T_w$ or heat flux $\dot{q}$ (depending on the type of the thermal boundary conditions).

The function for $P$ given by equation Equation  12.10-6 is relevant for the smooth walls. For the rough walls, however, this function is modified as follows:
 P_{\rm rough} = 3.15 {\rm Pr}^{0.695} \left( \frac{1}{E'} -... ...}{E} \right)^{0.359} + \left( \frac{E'}{E} \right)^{0.6} P (12.10-7)

where $E'$ is the wall function constant modified for the rough walls, defined by $E' = E/f_r$. To find a description of the roughness function $f_r$, you may refer to Equation  7.13-3 in Section  7.13.1.

Species

When using wall functions for species transport, FLUENT assumes that species transport behaves analogously to heat transfer. Similarly to Equation  12.10-5, the law-of-the-wall for species can be expressed for constant property flow with no viscous dissipation as

 Y^* \equiv \frac{\left(Y_{i,w} - Y_{i}\right) \rho C_{\mu}^... ...P_c \right] & \; \; \; (y^* > y^*_c) \end{array} \right . (12.10-8)

where $Y_{i}$ is the local species mass fraction, ${\rm Sc}$ and ${\rm Sc}_t$ are molecular and turbulent Schmidt numbers, and $J_{i,w}$ is the diffusion flux of species $i$ at the wall. Note that $P_c$ and $y^*_c$ are calculated in a similar way as $P$ and $y^*_T$, with the difference being that the Prandtl numbers are always replaced by the corresponding Schmidt numbers.

Turbulence

In the $k$- $\epsilon$ models and in the RSM (if the option to obtain wall boundary conditions from the $k$ equation is enabled), the $k$ equation is solved in the whole domain including the wall-adjacent cells. The boundary condition for $k$ imposed at the wall is

 \frac{\partial k}{\partial n} = 0 (12.10-9)

where $n$ is the local coordinate normal to the wall.

The production of kinetic energy, $G_k$, and its dissipation rate, $\epsilon$, at the wall-adjacent cells, which are the source terms in the $k$ equation, are computed on the basis of the local equilibrium hypothesis. Under this assumption, the production of $k$ and its dissipation rate are assumed to be equal in the wall-adjacent control volume.

Thus, the production of $k$ is computed from
 G_k \approx \tau_w \frac{\partial U}{\partial y} = \tau_w \frac{\tau_w}{\kappa \rho C_{\mu}^{1/4} k_P^{1/2} y_P} (12.10-10)

and $\epsilon$ is computed from
 \epsilon_P = \frac{C_{\mu}^{3/4} k_P^{3/2}}{\kappa y_P} (12.10-11)

The $\epsilon$ equation is not solved at the wall-adjacent cells, but instead is computed using Equation  12.10-11. $\omega$ and Reynolds stress equations are solved as detailed in Sections  12.5.3 and 12.7.9, respectively.

Note that, as shown here, the wall boundary conditions for the solution variables, including mean velocity, temperature, species concentration, $k$, and $\epsilon$, are all taken care of by the wall functions. Therefore, you do not need to be concerned about the boundary conditions at the walls.

The standard wall functions described so far are provided as a default option in FLUENT. The standard wall functions work reasonably well for a broad range of wall-bounded flows. However, they tend to become less reliable when the flow situations depart too much from the ideal conditions that are assumed in their derivation. Among others, the constant-shear and local equilibrium hypotheses are the ones that most restrict the universality of the standard wall functions. Accordingly, when the near-wall flows are subjected to severe pressure gradients, and when the flows are in strong non-equilibrium, the quality of the predictions is likely to be compromised.

The non-equilibrium wall functions offered as an additional option can improve the results in such situations.

figure   

Standard wall functions are available with the following viscous models:
  • K-epsilon

  • Reynolds Stress


next up previous contents index Previous: 12.10.1 Overview
Up: 12.10 Near-Wall Treatments for
Next: 12.10.3 Non-Equilibrium Wall Functions
© Fluent Inc. 2006-09-20