Law of the wall

(Redirected from Logarithmic law of the wall)
law of the wall, horizontal velocity near the wall with mixing length model

In fluid dynamics, the 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 by Theodore von Kármán, in 1930.[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^+ = \frac{1}{\kappa} \ln\, y^+ + C^+,$   with   $y^+ = \frac{y\, u_\tau}{\nu},$   $u_\tau=\sqrt{\frac{\tau_w}{\rho}}$   and   $u^+ = \frac{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τ, τw is the wall shear stress, ρ is the fluid density, uτ is called the friction velocity or shear velocity κ 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 κ≈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} = \frac{u_\tau}{\kappa} \ln\, \frac{y}{y_0}\$

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 δν and a characteristic roughness length-scale ks,[2]

 $k_s < \delta_\nu\,$ : hydraulically smooth flow, $k_s \approx \delta_\nu\,$ : transitional flow, $k_s > \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=\frac{u_\tau k_s}{\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=\frac{\nu}{9 u_\tau}$ for hydraulically smooth flow $y_0=\frac{k_s}{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.[6]

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

$k_s \approx 3.5 D_{84},\$

where D84 is the average diameter of the 84th largest percentile of the grains of the bed material.[7]

Power law solutions

Work by Barenblatt and others has 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.[8][9] In 1996, Cipra submitted experimental evidence in support of these power-law descriptions.[10] This evidence itself has not been fully accepted by other experts.[11]

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][12] All scaling laws derived in [12] should represent the leading order solutions of a corresponding averaged Navier-Stokes boundary value problem for sufficiently high Reynolds numbers, and hence should all be valid in spatial regimes where the influence of the boundaries can be ignored. This derivation from first principles however lacks mathematical rigor. It is based on several critical assumptions, mainly to justify the use of the pure temporal scaling symmetry of the inviscid Euler equations, a symmetry not admitted by the viscous Navier-Stokes equations, not even to leading order in any strict perturbative sense for small viscosities. The well-known challenging issues related to the vanishing viscosity limit in transiting from the viscous Navier-Stokes to the inviscid Euler equations,[13] particularly for bounded flows, also resides in every Lie group analysis for variable transformations, in that the limit is singularly unstable for symmetry transformations and non-unique for equivalence transformations.[14][15] These crucial issues were not addressed in the derivation, making this derivation from first principles thus highly questionable. The difference here between a symmetry and an equivalence transformation is that the former defines viscosity as an arbitrary but fixed parameter, while the latter contrarily defines viscosity as an own additional Lie group variable next to all coordinates and flow variables.[12] Eventually an ultimate test for these scaling laws allegedly emerging from first principles would be to show if the proposed Lie group methodology allows for a consistent unified treatment in agreement with experiment or DNS data when higher-order statistical quantities are included, in particular when involving the quantity of dissipation with its anomalous behavior in the vanishing viscosity limit. Otherwise the assumptions made in order to derive these proposed wall scaling laws in [12] distorts the conclusion of dealing with a rigorous derivation from first principles.

Near the Wall

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

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τ 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τ,

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
$u^+ \neq y^+$
$u^+ \neq \frac{1}{\kappa} \ln\, y^+ + C^+$

with the largest variation from either law occurring approximately where the two equations intercept, 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 k^+ function and the van Driest mixing length equation. Comparisons with DNS data of fully-developed turbulent channel flows for 109<Reτ<2003 showed good agreement.[17]

Notes

1. ^ von Kármán, Th. (1930), "Mechanische Ähnlichkeit und Turbulenz", Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen, Fachgruppe 1 (Mathematik) 5: 58–76 (also as: “Mechanical Similitude and Turbulence”, Tech. Mem. NACA, no. 611, 1931).
2. Mohrig, David (2004). "Conservation of Mass and Momentum". 12.110: Sedimentary Geology, Fall 2004. MIT OCW. Retrieved 2009-03-27.
3. ^ a b c Schlichting & Gersten (2000) pp. 522–524.
4. ^ Schlichting & Gersten (2000) p. 530.
5. ^ Whipple, Kelin (2004). "Hydraulic Roughness". 12.163: Surface processes and landscape evolution. MIT OCW. Retrieved 2009-03-27.
6. ^ Le Roux, J.P. (2004), "An integrated law of the wall for hydrodynamically transitional flow over plane beds", Sedimentary Geology 163 (3–4): 311–321, Bibcode:2004SedG..163..311L, doi:10.1016/j.sedgeo.2003.07.005
7. ^ Haws, Benjamin. "Equivalent sand roughness of Nikuradse (ks)". Retrieved 2009-03-27.
8. ^ Lynn Yarris. "A flaw in the law". Berkeley Lab: Highlights 97–98. Lawrence Berkeley National Laboratory, U.S. Department of Energy.
9. ^ Barenblatt, G.I. (1993), "Scaling laws for fully developed turbulent shear flows. Part 1. Basic hypotheses and analysis", Journal of Fluid Mechanics 248: 513–520, Bibcode:1993JFM...248..513B, doi:10.1017/S0022112093000874
Barenblatt, G.I.; Prostokishin, V.M. (1993), "Scaling laws for fully developed turbulent shear flows. Part 2. Processing of experimental data", Journal of Fluid Mechanics 248: 521–529, Bibcode:1993JFM...248..521B, doi:10.1017/S0022112093000886
Barenblatt, G.I.; Goldenfeld, N. (1995), "Does fully developed turbulence exist? Reynolds number independence versus asymptotic covariance", Physics of Fluids 7 (12): 3078–3084, arXiv:cond-mat/9507132, Bibcode:1995PhFl....7.3078B, doi:10.1063/1.868685
Barenblatt, G.I.; Chorin, A.J. (May 1996), Scaling laws and zero viscosity limits for wall-bounded shear flows and for local structure in developed turbulence, University of California at Berkeley Internal Report, PAM-678, Center for Pure and Applied Mathematics
10. ^ Cipra, B. (May 1996), "A New Theory of Turbulence Causes a Stir Among Experts", Science 272 (5264): 951, Bibcode:1996Sci...272..951C, doi:10.1126/science.272.5264.951
11. ^ Zagarola, M.V.; Perry, A.E.; Smits, A.J. (1997), "Log Laws or Power Laws: The Scaling in the Overlap Region", Physics of Fluids 9 (7): 2094–2100, Bibcode:1997PhFl....9.2094Z, doi:10.1063/1.869328
12. ^ a b c d Oberlack, Martin (2001), "A unified approach for symmetries in plane parallel turbulent shear flows", Journal of Fluid Mechanics 427: 299–328, doi:10.1017/S0022112000002408
13. ^ Constantin, Peter (2007), "On the Euler equations of incompressible fluids", Bulletin of the A. Math. Soc., 44(4): 603–621
14. ^ Ibragimov, N.H.; Kovalev, V.F. (2009), Approximate and Renormgroup Symmetries, Nonlinear Physical Science, Springer-Verlag
15. ^ Andreev, V.K.; Kaptsov, O.V.; Pukhnachev, V.V.; Rodionov, A.A. (1998), Applications of Group-Theoretical Methods in Hydrodynamics, Mathematics and its Applications, Springer-Verlag
16. ^ Turbulent Flows (2000) pp. 273–274.Pope, Stephen (2000), Turbulent Flows (1st revised ed.), Cambridge University Press, ISBN 0-521-59125-2
17. ^ Absi, Rafik (2009), "A simple eddy viscosity formulation for turbulent boundary layers near smooth walls", Comptes Rendus Mécanique 337 (3): 158–165, arXiv:1106.0985, Bibcode:2009CRMec.337..158A, doi:10.1016/j.crme.2009.03.010