The following wall functions are those appropriate to a near-wall layer in local equilibrium:

U^{+} = ln(E*y^{+})/κ

k = U_{τ}^{2}/√ℂ_{μ}

ε = ℂ_{μ}^{3/4}*k^{3/2}/(κ*δ)

ω=√k/(κ*ℂ_{μ}^{1/4}*δ)

(8.2.1)

(8.2.2)

(8.2.3)

(8.2.4)

where:

- U
^{+}= U_{r}/U_{Τ}; - U
_{r}is the absolute value of the resultant velocity parallel to the wall at the first grid node; - U
_{τ}is the resultant friction velocity ( = √(τ_{w}/ρ) ); - δ is the normal distance of the first grid point from the wall;
- y
^{+}is the dimensionless wall distance ( = U_{τ}*δ/ν ); - ν
_{l}is the laminar kinematic viscosity; - ℂ
_{μ}=C_{μ}C_{d}=0.09; - κ=0.41 is von Karman's constant; and
- E is a roughness parameter, which is equal to 8.6 for smooth walls.

Equation (8.2.1) is the well-known logarithmic law of the wall, and this law
should be applied to a point whose y^{+} value is in the range
30 ≤ y^{+} <130.

The boundary condition (8.2.2) for the turbulent kinetic energy k assumes that the turbulence is in local equilibrium, and consequently, this set of wall functions are not really suitable under separated conditions because turbulent energy diffusion towards the wall is significant, leading to appreciable departures from local equilibrium.

The wall function defined by eqn (8.2.1) is implemented in the momentum equations by
way of source terms per unit wall surface area, S_{m}, which take the form:

S_{mo} = ρ*s*|U_{r}|*(V_{w}-V)

(8.2.5)

where ρ denotes the fluid density, s the friction factor ( = τ_{w}/(ρ*U_{r}^{2}) ), V
denotes the in-cell value of velocity, and V_{w} denotes the value of the
velocity at the wall. The friction factor is s determined from the laminar and turbulent friction factors,
s_{l} and s_{t}, as follows:

s = max(s_{l}, s_{t})

s_{l} = 1./Re

s_{t} = [κ/ln(E*Re*√s_{t})]^{2}

Re = U_{r}*δ/ν_{l}

(8.2.6)

(8.2.7)

(8.2.8)

(8.2.9)

and the expression for s_{t} follows from equation (8.2.1), and its
value is determined by iteration from equation (8.2.8).

For heat and mass transfer at the wall, the flux of the variable φ between the fluid and the wall is given by:

S_{φ} = ρ*St*|U_{r}|*(φ_{w}-φ)

(8.2.10)

in which φ denotes the in-cell value of the heat- or mass-transfer
variable, and φ_{w} denotes the value of φ at the wall. The Stanton
number St is defined by:

St ≡ -q_{w}/(ρ*C_{p}*|U_{r}|*(φ-φ_{w}) = h/(ρ*C_{p}*|U_{r}|)

(8.2.11)

where h is the local heat transfer coefficient; q_{w} is the wall heat flux, and C_{p} is the specific heat.

The Stanton number in equation (8.2.10) is computed from:

St = max(St_{l}, St_{t})

St_{l} = 1./(σ_{t}*Re)

St_{t} = s_{t}/[σ_{t}*(1.+P_{m}*√s_{t})]

P_{m} = 9.*(σ_{l}/σ_{t}-1)*(σ_{t}/σ_{l})^{1/4}

(8.2.12)

(8.2.13)

(8.2.14)

(8.2.15)

In the foregoing equations, σ_{t} and σ_{l} denote
the turbulent and laminar Prandtl numbers respectively, and equation (8.2.15)
defines the smooth-wall sublayer-resistance function P_{m} of Jayatilleke
[1969], which is valid for moderate to high Prandtl number fluids, i.e.
for σ_{l} ≥0.5.

Back to wall-functions contents list