Law of the wall explained

In fluid dynamics, the law of the wall (also known as the logarithmic law of the wall) states that the average velocity of a turbulent flow at a certain point is proportional to the logarithm of the distance from that point to the "wall", or the boundary of the fluid region. This law of the wall was first published in 1930 by Hungarian-American mathematician, aerospace engineer, and physicist Theodore von Kármán.[1] It is only technically applicable to parts of the flow that are close to the wall (<20% of the height of the flow), though it is a good approximation for the entire velocity profile of natural streams.[2]

General logarithmic formulation

The logarithmic law of the wall is a self similar solution for the mean velocity parallel to the wall, and is valid for flows at high Reynolds numbers — in an overlap region with approximately constant shear stress and far enough from the wall for (direct) viscous effects to be negligible:[3]

u+=

1
\kappa

lny++C+,

with

y+=

yu\tau
\nu

,

u
\tau=\sqrt{\tauw
\rho
} and

u+=

u
u\tau

where

y+

is the wall coordinate: the distance y to the wall, made dimensionless with the friction velocity uτ and kinematic viscosity ν,

u+

is the dimensionless velocity: the velocity u parallel to the wall as a function of y (distance from the wall), divided by the friction velocity uτ,

\tauw

is the wall shear stress,

\rho

is the fluid density,

u\tau

is called the friction velocity or shear velocity,

\kappa

is the Von Kármán constant,

C+

is a constant, and

ln

is the natural logarithm.
From experiments, the von Kármán constant is found to be

\kappa ≈ 0.41

and

C+ ≈ 5.0

for a smooth wall.[3]

With dimensions, the logarithmic law of the wall can be written as:[4]

{u}=

u\tau
\kappa

ln

y
y0

where y0 is the distance from the boundary at which the idealized velocity given by the law of the wall goes to zero. This is necessarily nonzero because the turbulent velocity profile defined by the law of the wall does not apply to the laminar sublayer. The distance from the wall at which it reaches zero is determined by comparing the thickness of the laminar sublayer with the roughness of the surface over which it is flowing. For a near-wall laminar sublayer of thickness

\delta\nu

and a characteristic roughness length-scale

ks

,[2]

ks<\delta\nu

hydraulically smooth flow,

ks\delta\nu

transitional flow,

ks>\delta\nu

hydraulically rough flow.

Intuitively, this means that if the roughness elements are hidden within the laminar sublayer, they have a much different effect on the turbulent law of the wall velocity profile than if they are sticking out into the main part of the flow.

This is also often more formally formulated in terms of a boundary Reynolds number,

Rew

, where
Re
w=u\tauks
\nu

.

The flow is hydraulically smooth for

Rew<3

, hydraulically rough for

Rew>100

, and transitional for intermediate values.[2]

Values for

y0

are given by:[2] [5]
y
0=\nu
9u\tau
 
for hydraulically smooth flow
y
0=ks
30
for hydraulically rough flow.

Intermediate values are generally given by the empirically derived Nikuradse diagram,[2] though analytical methods for solving for this range have also been proposed.

For channels with a granular boundary, such as natural river systems,

ks3.5D84,

where

D84

is the average diameter of the 84th largest percentile of the grains of the bed material.[6]

Power law solutions

Works by Barenblatt and others have shown that besides the logarithmic law of the wall — the limit for infinite Reynolds numbers — there exist power-law solutions, which are dependent on the Reynolds number.[7] [8] In 1996, Cipra submitted experimental evidence in support of these power-law descriptions. This evidence itself has not been fully accepted by other experts. In 2001, Oberlack claimed to have derived both the logarithmic law of the wall, as well as power laws, directly from the Reynolds-averaged Navier–Stokes equations, exploiting the symmetries in a Lie group approach.[3] However, in 2014, Frewer et al. refuted these results.

For scalars

For scalars (most notably temperature), the self-similar logarithmic law of the wall has been theorized (first formulated by B. A. Kader[9]) and observed in experimental and computational studies.[10] [11] [12] [13] In many cases, extensions to the original law of the wall formulation (usually through integral transformations) are generally needed to account for compressibility, variable-property and real fluid effects.

Near the wall

Below the region where the law of the wall is applicable, there are other estimations for friction velocity.[14]

Viscous sublayer

In the region known as the viscous sublayer, below 5 wall units, the variation of

u+

to

y+

is approximately 1:1, such that:

For

y+<5

u+=y+

where,

y+

is the wall coordinate: the distance y to the wall, made dimensionless with the friction velocity

u\tau

and kinematic viscosity

\nu

,

u+

is the dimensionless velocity: the velocity u parallel to the wall as a function of y (distance from the wall), divided by the friction velocity

u\tau

,
This approximation can be used farther than 5 wall units, but by

y+=12

the error is more than 25%.

Buffer layer

In the buffer layer, between 5 wall units and 30 wall units, neither law holds, such that:

For

5<y+<30

u+y+

u+

1
\kappa

lny++C+

with the largest variation from either law occurring approximately where the two equations intersect, at

y+=11

. That is, before 11 wall units the linear approximation is more accurate and after 11 wall units the logarithmic approximation should be used, though neither are relatively accurate at 11 wall units.

The mean streamwise velocity profile

u+

is improved for

y+<20

with an eddy viscosity formulation based on a near-wall turbulent kinetic energy

\kappa+

function and the van Driest mixing length equation. Comparisons with DNS data of fully developed turbulent channel flows for

109<Re\tau<2003

showed good agreement.

External links

Notes and References

  1. (also as: “Mechanical Similitude and Turbulence”, Tech. Mem. NACA, no. 611, 1931).
  2. Web site: Conservation of Mass and Momentum. Mohrig. David. 2004. 12.110: Sedimentary Geology, Fall 2004. MIT OCW. 2009-03-27.
  3. Schlichting & Gersten (2000) pp. 522–524.
  4. Schlichting & Gersten (2000) p. 530.
  5. Web site: Hydraulic Roughness. Whipple. Kelin. 2004. 12.163: Surface processes and landscape evolution. MIT OCW. 2009-03-27.
  6. Web site: Equivalent sand roughness of Nikuradse (ks) . Haws . Benjamin . 2009-03-27 .
  7. Web site: A flaw in the law . Lynn Yarris . Lawrence Berkeley National Laboratory, U.S. Department of Energy . Berkeley Lab: Highlights 97–98 .



  8. Kader . B. A. . 1981-09-01 . Temperature and concentration profiles in fully turbulent boundary layers . International Journal of Heat and Mass Transfer . en . 24 . 9 . 1541–1544 . 10.1016/0017-9310(81)90220-9 . 0017-9310.
  9. Simonich . J. C. . Bradshaw . P. . 1978-11-01 . Effect of Free-Stream Turbulence on Heat Transfer through a Turbulent Boundary Layer . Journal of Heat Transfer . 100 . 4 . 671–677 . 10.1115/1.3450875 . 0022-1481.
  10. Patel . Ashish . Boersma . Bendiks J. . Pecnik . Rene . 2017-08-21 . Scalar statistics in variable property turbulent channel flows . Physical Review Fluids . 2 . 8 . 084604 . 10.1103/PhysRevFluids.2.084604. 2017PhRvF...2h4604P .
  11. Toki . Takahiko . Teramoto . Susumu . Okamoto . Koji . 2020-01-01 . Velocity and Temperature Profiles in Turbulent Channel Flow at Supercritical Pressure . Journal of Propulsion and Power . 36 . 1 . 3–13 . 10.2514/1.B37381. 209963353 .
  12. Guo . J. . Yang . X. I. A. . Ihme . M. . March 2022 . Structure of the thermal boundary layer in turbulent channel flows at transcritical conditions . Journal of Fluid Mechanics . en . 934 . 10.1017/jfm.2021.1157 . 2022JFM...934A..45G . 0022-1120. free .
  13. Turbulent Flows (2000) pp. 273–274.