(Ebook) Turbulence in the Atmosphere by John C. Wyngaard ISBN 9780521887694, 0521887690 All Chapters Instant Download
(Ebook) Turbulence in the Atmosphere by John C. Wyngaard ISBN 9780521887694, 0521887690 All Chapters Instant Download
(Ebook) Turbulence in the Atmosphere by John C. Wyngaard ISBN 9780521887694, 0521887690 All Chapters Instant Download
https://ebooknice.com/product/biota-grow-2c-gather-2c-cook-6661374
ebooknice.com
https://ebooknice.com/product/turbulence-structure-and-vortex-
dynamics-991848
ebooknice.com
https://ebooknice.com/product/essentials-of-meteorology-an-invitation-
to-the-atmosphere-4060912
ebooknice.com
https://ebooknice.com/product/essentials-of-meteorology-an-invitation-
to-the-atmosphere-7344156
ebooknice.com
(Ebook) Turbulence Modeling for CFD (Third Edition) by Wilcox, David
C. ISBN 9781928729082, 1928729088
https://ebooknice.com/product/turbulence-modeling-for-cfd-third-
edition-55969182
ebooknice.com
https://ebooknice.com/product/turbulence-modeling-for-cfd-third-
edition-2358078
ebooknice.com
https://ebooknice.com/product/essentials-of-meteorology-an-invitation-
to-the-atmosphere-11326010
ebooknice.com
https://ebooknice.com/product/large-scale-atmosphere-ocean-
dynamics-i-1774724
ebooknice.com
https://ebooknice.com/product/volatile-organic-compounds-in-the-
atmosphere-1558098
ebooknice.com
Turbulence in the Atmosphere 1st Edition John C.
Wyngaard Digital Instant Download
Author(s): John C. Wyngaard
ISBN(s): 9780521887694, 0521887690
Edition: 1
File Details: PDF, 3.99 MB
Year: 2010
Language: english
This page intentionally left blank
TURBULENCE IN THE ATMOSPHERE
John Wyngaard’s experience in turbulence research and teaching spans the Air
Force Cambridge Research Laboratories, the Wave Propagation Laboratory of the
National Oceanographic and Atmospheric Administration (NOAA) in Boulder, the
Atmospheric Analysis and Prediction Division of the National Center for Atmo-
spheric Research (NCAR), and the Department of Meteorology at Pennsylvania
State University, where he developed a sequence of courses on turbulence. This
book is based on those courses. He has published over 100 refereed journal papers
covering theoretical, observational, and numerical modeling aspects of engineering
and geophysical turbulence.
TURBULENCE IN THE
ATMOSPHERE
JOHN C. WYNGAARD
CAMBRIDGE UNIVERSITY PRESS
Cambridge, New York, Melbourne, Madrid, Cape Town, Singapore,
São Paulo, Delhi, Dubai, Tokyo
Published in the United States of America by Cambridge University Press, New York
www.cambridge.org
Information on this title: www.cambridge.org/9780521887694
© John C. Wyngaard 2010
Preface page ix
1 Introduction 3
1.1 Turbulence, its community, and our approach 3
1.2 The origins and nature of turbulence 4
1.3 Turbulence and surface fluxes 5
1.4 How do we study turbulence? 10
1.5 The equations of turbulence 11
1.6 Key properties of turbulence 14
1.7 Numerical modeling of turbulent flows 18
1.8 Physical modeling of turbulent flows 20
1.9 The impact of Kolmogorov 20
v
vi Contents
4 Turbulent fluxes 75
4.1 Introduction 75
4.2 Temperature flux in a boundary layer 76
4.3 Mass flux in scalar diffusion 78
4.4 Momentum flux in channel flow 83
4.5 The “mixture length” 84
4.6 Summary 85
Index 387
Preface
ix
x Preface
several years. Owen Coté and I, with the plotting and programming assistance
of Jack Izumi and Jean O’Donnell, mapped out perhaps the first observational
analyses of the conservation equations for stress and scalar flux in a turbulent
flow. Chandran Kaimal did his inimitable spectral analyses. Visitor Joost Businger
oversaw the analysis of the flux-mean profile relations over the wide range of
stability conditions. Henk Tennekes and I found that velocity derivative statistics
in the huge Reynolds number Kansas turbulence were off the older charts, but in
accord with the newer thinking of the Russian school. Those were heady times.
We returned to the field in 1973, in the very flat farming country of northwestern
Minnesota. In collaboration with a British Met Office group we reached deep into
the boundary layer with sensors on the tethering cable of a World War II surplus
barrage balloon. It stayed up for several weeks, despite the best late-night efforts
of rifle-toting cowboys, until it was brought down by a gust front. Part II discusses
some of the insights gained from these field programs.
Today’s main types of turbulent-flow models – large-eddy simulation and second-
order closure (Chapters 5 and 6) – were in their infancy in 1970. I remember
discovering with Owen Coté the myriad ways a second-order-closure model could
misbehave – negative variances, violations of Schwartz’s inequality, … What we
thought were obviously better closures gave poorer results. We saw the early hopes
of universal second-order closures dashed in buoyancy-dominated flows. We devel-
oped a wariness about turbulence modeling. As Ronald Reagan later said, “Trust,
but verify.”
Anyone who has developed models of the second-moment equations (Chapter 5),
discovered how poorly they can behave, and then in fatigue and discouragement
wondered how Nature keeps variances positive, can appreciate this story:
Some years ago, during the hall talk at a break in an NCAR meeting, a prominent senior sci-
entist became impatient with a mathematician’s fussing over obscure details of an equation.
“Hell,” he blurted, “in the atmospheric sciences we don’t even know what the equations are.”
The applied turbulence field seems different today. Numerical modeling of tur-
bulent flows is a dominant technology used by a second- or even third-generation
community. Programmers have ensured that the codes don’t misbehave like they
used to. Geophysical observations have not kept pace with the model predictions,
nor could they have; modeling and observational work have cruelly different time
scales. Now less likely to be rooted in personal experience, wariness of modeling
seems to be diminishing.
Recently I previewed a video of the EPA Fluid-Modeling Facility before showing
it in my class on atmospheric dispersion. The FMF, as it is called, is located in
Raleigh, NC, and contains low-speed wind tunnels, a stratified towing tank, and
a replica of the Deardorff–Willis convection tank (Chapter 11). The FMF is a
world-class facility put together largely by Bill Snyder beginning about 1970.
Preface xi
subrange. The final chapter covers the dissipative range, both as hypothesized in
1941 by Kolmogorov and more recently through dissipation-intermittency models,
and two-dimensional turbulence as described through the Kolmogorov-like notions
of Kraichnan and Batchelor.
Part II covers turbulence in the atmospheric boundary layer (ABL). The first
chapter generalizes the equations of Part I to a variable density environment in a
standard way, using a background-plus-deviation representation for density, tem-
perature, and pressure. The background state is hydrostatic, buoyancy is handled
through the Boussinesq approximation, and a conserved temperature is used that
in its ultimate form allows phase change. The four subsequent chapters survey the
structure and dynamics of the ABL, emphasizing for a non-meteorological audi-
ence those features that make its turbulence different from that in engineering flows.
They also cover turbulence in the surface layer, and discuss in depth the physics and
the efficacy of the Monin–Obukhov similarity hypothesis for its turbulence struc-
ture. There is a chapter on the convective boundary layer, whose turbulence physics
and structure have been extensively studied in the field and through large-eddy
simulation. The final chapter covers the stable ABL, which has some regimes in
which turbulence structure and dynamics have a reasonably simple interpretation.
Part III, “Statistical representation of turbulence,” includes a number of important
statistical tools and concepts – probability densities and distributions, covariances,
autocorrelations, spectra, and local isotropy – that are used in turbulence and other
stochastic problems. It has a number of illustrative examples of stochastic prob-
lems that can be solved analytically, including the wavenumber-space dynamics
of turbulence spectra; relating spectra in the plane to traditional spectra; and the
effects of spatial averaging, sensor separation, crosstalk, and probe-induced flow
distortion on turbulence measurements.
In the course of writing this book I received valuable input on technical mat-
ters from Bob Antonia, Bob Beare, Craig Bohren, Frank Bradley, Peter Bradshaw,
Jim Brasseur, Joost Businger, Steve Clifford, Steve Derbyshire, Diego Donzis,
Carl Friehe, Steve Hatlee, Reg Hill, Bert Holtslag, Tom Horst, Mark Kelly, Don
Lenschow, Charles Meneveau, Chin-Hoh Moeng, Parviz Moin, Ricardo Munoz,
Laurent Mydlarski, Bill Neff, Ray Shaw, K. R. Sreenivasan, Peter Sullivan, Dennis
Thomson, Chenning Tong, Zellman Warhaft, Jeff Weil, Keith Wilson, P. K. Yeung,
and Sergej Zilitinkevich, for which I am most grateful. I’d like to thank Lori Mattina
for expertly and patiently crafting the figures, Ned Patton for kindly setting up the
LATEX style files for me, and Peter Sullivan again for his generous and sustained
assistance with LATEX. I am grateful to the AFCRL group – Duane Haugen, Chan-
dran Kaimal, Owen Coté, Jack Izumi, Jim Newman, Jean O’Donnell, and Don
Stevens – for my once-in-a-lifetime experience in the 1968 Kansas experiment.
Finally, I thank John Lumley for inspiring my career in turbulence.
Part I
A grammar of turbulence
1
Introduction
3
4 Introduction
Figure 1.1 Instability of an axisymmetric jet. A laminar stream of air flows from a
circular tube at the left at Reynolds number 10 000 and is made visible by a smoke
wire. The edge of the jet develops axisymmetric oscillations, rolls up into vortex
rings, and then abruptly becomes turbulent. Photograph courtesy Robert Drubka
and Hassan Nagib. From Van Dyke (1982).
† For example, in Figure 1.1 U is the velocity averaged over the tube cross section and L is the tube diameter.
‡ To paraphrase Batchelor (1950), “eddy” does not refer to any specific local distribution of velocity; it is simply a
concise term for local turbulent motion with a certain length scale – an arbitrary local flow pattern characterized
by size alone. A turbulent flow has a spectrum of eddies of different size, determined by an analysis of the
velocity field into sinusoidal components of different wavelengths (Chapter 15).
1.3 Turbulence and surface fluxes 5
hr
f
uave D
Re=
n
Figure 1.2 The Moody chart, which shows the behavior of the Darcy friction factor
f , Eq. (1.5), in a circular pipe. In laminar flow f ∝ Re−1 , Eq. (1.6); f jumps to
larger values with the transition to turbulence at Re 2000, and in the region of
equilibrium turbulence past the critical zone f depends also on the wall-roughness
height hr relative to D. Adapted from Moody (1944).
where r is the radial coordinate, R = D/2 is the pipe radius, and umax is the
maximum (centerline) velocity. The velocity averaged over the cross section is
R
1 umax
uave = u(r)2π r dr = . (1.2)
π R2 0 2
with μ the dynamic viscosity of the fluid. Since ∂p/∂x does not depend on x
(Problem 1.1), we can write the axial force balance on a slug of fluid of length L
and diameter D as
∂P π D 2 ∂P
τwall π DL = − L , so that − D = 4τwall . (1.4)
∂x 4 ∂x
− ∂P
∂x D 4τwall
f ≡ 2
= . (1.5)
ρ(uave ) /2 ρ(uave )2 /2
64μuave 64
flam = 2
= . (1.6)
Dρ(uave ) Re
† The Fanning friction factor is the wall stress nondimensionalized with ρ(uave )2 /2. The Darcy friction factor,
Eq. (1.5), is larger by a factor of four.
1.3 Turbulence and surface fluxes 7
Re
Figure 1.3 The ratios of mean fluxes at the wall in turbulent and laminar flow
through smooth pipes. The momentum-flux ratio is Eq. (1.7) evaluated with f
data from Figure 1.2; the heat-flux ratio is Eq. (1.16) evaluated with Nu data from
Dittus and Boelter (1930), as summarized by Turns (2006).
This ratio is plotted for smooth pipes in Figure 1.3. It has very large values at large
Re, indicating the strong influence of turbulence on the wall stress.
Turns (2006) shows that a good fit to the classical mean-velocity measurements
of Nikuradse (1933) in turbulent pipe flow is
u(r) f 1/2 Re f 1/2 r
= √ 2.5ln √ 1− + 5.5 . (1.8)
uave 2 2 2 R
Figure 1.4 shows that this profile is much “flatter” in the core region than the laminar
profile (1.1). At large Re the mean-velocity gradient is significant only adjacent
to the wall, where it is much larger than in laminar flow of the same bulk fluid
velocity. The wall stress in turbulent flow is still defined by the velocity gradient at
the wall, Eq. (1.3), but that gradient, and therefore the wall shear stress, fluctuates
chaotically with time and with position. The mean value (which we designate by
an overbar) of the wall stress is
∂u ∂u
τ wall = −μ = −μ . (1.9)
∂r r=R ∂r r=R
We have used the property that the differentiation and averaging can be done in
either order (Problem 1.3). The sharp increase in the mean-velocity gradient at the
wall at transition (Figure 1.4) causes a sharp increase in wall stress (Figure 1.2).
8 Introduction
u
uave
Figure 1.4 Profiles of u/uave in turbulent pipe flow, Eq. (1.8), “flatten” as the
Reynolds number increases, making the mean shear and mean stress at the wall
much larger than in laminar flow with the same average velocity.
Pipes typically have some wall roughness, and Figure 1.2 indicates that the mean
wall stress increases with that roughness. The explanation (Kundu, 1990) is that
immediately adjacent to the wall in turbulent pipe flow is a laminar sublayer of
thickness δ ∼ 5ν/u∗ , with u∗ = (τ wall /ρ)1/2 the friction velocity. If the typ-
ical height hr of the individual “bumps” or roughness elements on the wall is
much less than δ, wall roughness has minimal effect and the mean wall stress is
the viscous one given by Eq. (1.9). But as hr approaches δ the roughness ele-
ments cause form drag through the pressure distribution on their surface, which
adds to the viscous drag and increases the friction factor f . When hr is large
enough this form drag dominates and f ceases to change with Re, as indicated in
Figure 1.2.
An analogous situation exists for the wall heat flux Hwall (watts m−2 ). It is carried
entirely by the molecular diffusion process called conduction heat transfer:
∂T
Hwall = −k , (1.10)
∂r r=R
with k the thermal conductivity (watts m−1 K−1 ). The heat flux is continuous at
the fluid–wall interface, but the temperature gradient there is discontinuous if k of
the wall material and the fluid differ. We shall consider the fluid side.
1.3 Turbulence and surface fluxes 9
The temperature profile in fully developed laminar pipe flow depends on the
temperature boundary conditions. We’ll consider the analytically simple case where
the fluid and wall temperatures vary linearly with x but their difference, and the
wall heat flux, are independent of x. Its temperature profile is (Problem 1.2)
R 2 uave ∂T r2 r2
T (0, x) − T (r, x) = − 1 − , (1.11)
α ∂x 2R 2 4R 2
with α = k/(ρcp ) the thermal diffusivity of the fluid. The relation between the wall
heat flux and ∂T /∂x is (Problem 1.2)
Duave ρcp ∂T
Hwall = − , (1.12)
4 ∂x
so the temperature profile (1.11) can be rewritten as
Hwall D r 2 r2
T (0, x) − T (r, x) = 1− . (1.13)
k 2R 2 4R 2
The wall heat flux made dimensionless with pipe diameter D, fluid thermal
conductivity k, and a temperature difference T is called a Nusselt number Nu:
Hwall D
Nu = . (1.14)
k T
T is defined through the wall temperature Tw (x) and the “bulk fluid temperature”
at that position, Tb (x):
R
u(r) T (r, x) 2π r dr Hwall D
Tb (x) = 0
, Nu = . (1.15)
2
π R uave k(Tb − Tw )
Turns (2006) shows that in the laminar case in this problem Nu 4.4.
In the turbulent case the heat flux at a point on the wall, like the stress there,
fluctuates chaotically in time. The turbulent mixing makes the mean temperature
gradient relatively small over most of the cross section; it is large only near the wall,
as for velocity (Figure 1.4). The mean wall heat flux, the product of the fluid thermal
conductivity k and the mean temperature gradient at the wall, is much larger than
in the laminar case.
From Eq. (1.14) in this problem we can write the ratio of wall heat fluxes as, for
given values of D, k, and T ,
H wall Nu Nu
= = . (1.16)
Hwall (laminar flow) Nu(laminar flow) 4.4
10 Introduction
Figure 1.3 shows the wall-flux ratios in the turbulent regime. The ratios for heat
and momentum differ by a constant factor of about 1.5–2 (a manifestation of the
Reynolds analogy between heat and momentum transfer) as they increase sharply
with Reynolds number.
If the flow over the earth’s surface were laminar, not turbulent, the environmental
effects would be profound. In clear summer weather, for example, the earth’s surface
temperature could routinely approach 100 ◦ C during the day and 0 ◦ C at night
(Chapter 9).
cited as the first LES. DNS and LES, which now use as many as 40963 6 × 1010
grid points, are leading research tools in turbulence today.
with p the pressure, gi the gravity vector and σij the viscous stress tensor. The left
side of Eq. (1.19) is density times the total acceleration following the motion, the
sum of local and advective accelerations; the right side is the sum of the pressure-
gradient, gravity, and viscous forces. Equation (1.19) requires a nonaccelerating
coordinate system, but our earth-based coordinate system is accelerating because
the earth rotates; as a result (1.19) also needs a Coriolis term. This can be important
in atmospheric turbulence (Part II) but we shall ignore it in Part I.
Batchelor (1967) points out that when ρ is uniform one can define a pressure
s
(p , say) whose gradient exactly balances the gravity force:
∂p s
0=− + ρgi . (1.20)
∂xi
p = ps + pm , (1.21)
with pm a modified pressure that is due to the fluid motion, then we can write
Eq. (1.19) as
∂ui ∂ui ∂p m ∂σij
ρ + uj =− + , (1.22)
∂t ∂xj ∂xi ∂xj
so the gravity term does not appear explicitly. In Part I we shall use the form (1.22)
and drop the superscript m on pressure.
Using Eq. (1.18) we can write the momentum equation (1.22) as
∂ui ∂p ∂
ρ =− + −ρui uj + σij . (1.23)
∂t ∂xi ∂xj
The final term in Eq. (1.23) is in flux form. It can be interpreted as the divergence of
the total flux of momentum, the sum of advective and viscous parts. Momentum is
mass times velocity; it is a vector. The momentum flux is the amount of momentum
passing through a unit area per unit time. It is a second-order tensor quantity; it
involves two directions, that of the unit normal to the area and that of the momentum.
Its units, density times velocity squared, are equivalent to (newtons/m2 ), or stress.
Thus, we can also interpret the final term in Eq. (1.23) as the divergence of a
generalized stress.
In a incompressible Newtonian fluid the viscous stress tensor σij is a linear
function of the strain-rate tensor sij . We write this as
∂ui ∂uj
σij = μ + = 2μsij , (1.24)
∂xj ∂xi
1.5 The equations of turbulence 13
It is usual to divide by density and use Eqs. (1.24) and (1.25) to write Eq. (1.19)
for a Newtonian fluid with constant viscosity as
This says that the total time derivative of vorticity is the sum of a term representing
the interaction of vorticity and the velocity gradient, which we will interpret shortly,
and a molecular-diffusion term.
The statement of mass conservation for a scalar c that has no sources or sinks
(such as the mass density of a nonreacting trace constituent in the fluid) is
∂c ∂cui ∂ 2c
+ =γ , (1.29)
∂t ∂xi ∂xi ∂xi
where γ is the molecular diffusivity of c in the fluid. We can also write Eq. (1.29)
in flux form,
∂c ∂ ∂c
=− cui − γ , (1.30)
∂t ∂xi ∂xi
which says that local time changes in c are due to the divergence of the total flux
of c, the sum of advective and molecular components.
In Part I we are considering constant-density fluids, in which the velocity
divergence vanishes, so Eq. (1.29) can also be written
Dc ∂c ∂c ∂ 2c
= + ui =γ . (1.31)
Dt ∂t ∂xi ∂xi ∂xi
14 Introduction
This says that following the fluid motion and neglecting molecular diffusion, c does
not change. We call such a scalar a conserved scalar.†
In flows where heating due to radiation, phase change, chemical reactions, and
viscous effects is negligible the thermal energy equation reduces to the same form
as Eq. (1.31),
DT ∂T ∂T ∂ 2T
= + ui =α , (1.32)
Dt ∂t ∂xi ∂xi ∂xi
where T is temperature and α is the thermal diffusivity of the fluid. Equation (1.32)
says that under these conditions temperature is a conserved variable, changing only
through conduction heat transfer.
If ∂u1 /∂x1 is positive (1.33) says the vortex is stretched in the x1 direction,
increasing the magnitude of ω1 . ∂u2 /∂x1 and ∂u3 /∂x1 can generate ω2 and ω3
from ω1 ; this is sometimes called vortex tilting.
In two-dimensional turbulence the velocity field is ui = [u1 (x, y), u2 (x, y), 0],
say. Then ωi = (0, 0, ω3 ) and the vortex-stretching term in Eq. (1.28) for ω3
is ω3 ∂u3 /∂x3 = 0. This demonstrates that three dimensionality is necessary for
vortex stretching.
A cascade of kinetic energy through eddies of diminishing size (Chapters 6
and 7) that terminates in viscous dissipation – the conversion of kinetic energy
into internal energy by viscous forces in the smallest eddies – is a defining feature
† Unfortunately, if the velocity divergence is nonzero the density of a mass-conserving species is not a conserved
scalar. The ratio of its density and the fluid density is a conserved scalar, however (Part II).
‡ Parts of this discussion are adapted from Lumley and Panofsky (1964).
1.6 Key properties of turbulence 15
† This distinction between random and stochastic is not always made in turbulence, but it is useful.
16 Introduction
far more rapidly in the mean chimney plume than in a laminar plume, but within any
realization local effluent concentrations can be far higher than in the mean plume.
† We use ∼ in the sense of Tennekes and Lumley (1972), the implied proportionality coefficient being between
1/5 and 5, say.
‡ This and other classic papers on turbulence are included in a collection by Friedlander and Topper (1961).
§ As we discuss in Chapter 7, Kolmogorov and others subsequently modified this hypothesis to allow for the
effects of dissipation intermittency.
1.6 Key properties of turbulence 17
for these are the only combinations of and ν that produce a length scale and
a velocity scale. The Reynolds number of the dissipative eddies is therefore
∼ υη/ν = 1, which confirms that they are strongly influenced by viscosity.
Using ∼ u3 / we can write, using Eq. (1.34)
3/4
1/4 u3/4 u 3/4
= 3/4 ∼ 1/4 3/4 = = Rt ,
η ν ν ν
(1.35)
u u u 1/4
= 1/4
∼ 3 = Rt .
υ ( ν) (u ν/)1/4
Equation (1.35) implies that in large-Rt turbulent flows the dissipative eddies are
quite weak and quite small compared to the energy-containing eddies. If, for exam-
ple, u ∼ 1 m s−1 and ∼ 103 m, as is typical in the atmospheric boundary layer,
then Rt ∼ 108 so that υ ∼ 10−2 m s−1 and η ∼ 10−3 m.
The ratio of vorticities typical of the dissipative and energy-containing eddies is
at large Rt the dissipative eddies contain essentially all the turbulent vorticity.
The smallness of η/ in large-Rt turbulence puts severe limits on the Rt that can
be reached in direct numerical simulations of turbulence (Problem 1.9). Although
few turbulent flows of practical importance have Rt values small enough to be
calculated in this way, the concept of “Reynolds number similarity” (Chapter 2)
does make them useful.
The Kolmogorov microscale η, which from Eq. (1.35) can be written as
ν 3/4 1/4
η∼ , (1.37)
u3/4
is seldom smaller than 10−4 m in engineering flows and about 10−3 m in the atmo-
sphere. This is almost always large enough to ensure the applicability of continuum
fluid mechanics.
cascade of kinetic energy from large scales to small (Chapters 6, 7). This cascade
dynamically “couples” all the eddies in a turbulent flow. Equation (1.35) shows
that the ratio of energy-containing and dissipative eddy sizes increases as the 3/4
power of the large-eddy Reynolds number Rt , so a turbulent flow of large Rt has
a huge number of interacting eddies. This has thwarted all attempts to solve the
turbulence equations analytically. Even at relatively small Rt the “bookkeeping”
for these interactions overwhelms even our largest computers.
∂ui ∂ui uj 1 ∂p
+ =− , (1.38)
∂t ∂xj ρ ∂xi
with the overbar denoting the average. We have assumed the averaged viscous term
is negligible (Problem 1.8). In writing (1.38) we have used the zero-divergence
property of uj to bring it into the derivative. But the averaging has created a ui uj
term in (1.38); in turbulent flow it differs from ui uj (Problem 1.17) and, hence, is
an unknown. If we write ui uj as
τij
ui uj = ui uj + ui uj − ui uj = ui uj − , (1.39)
ρ
Averaging the Navier–Stokes equation has produced Eq. (1.40) for the average
velocity field, but the equation contains a new term involving a turbulent stress τij .
We’ll see that this is a subtly different quantity for ensemble and space averages,
but in each case it is typically called the Reynolds stress. Experience has shown that
in turbulent flow it is always important.
1.7 Numerical modeling of turbulent flows 19
In the same way, averaging Eq. (1.29) for a conserved scalar and neglecting its
averaged molecular term produces
∂c ∂c ui ∂fi
+ =− , (1.41)
∂t ∂xi ∂xi
weighting of all contributions; in filtering the weighting function can vary with
position. The concept is to apply to the governing equations a spatial filter that
removes the smallest eddies and leaves the energy-containing ones unaffected,
producing equations for the larger-scale part of each variable. If L is the spa-
tial scale of the computational domain and f is the cutoff scale of the filter –
that is, the filter removes spatial variations of scale less than f – then of order
N = L/ f computational grid points are required in each direction to resolve the
filtered fields, of order N 3 grid points in all. For this reason the filtered variables are
often called resolvable-scale variables. Today N is typically in the range 30–300.
This modeling approach now goes by the apt name large-eddy simulation, or LES.
It was proposed by Lilly (1967) and first used successfully by Deardorff for turbu-
lent channel flow (Deardorff, 1970a) and then for the atmospheric boundary layer
(Deardorff, 1970b).
Earlier work (as summarized by Taylor (1935), for example) had shown that
a turbulent flow has large eddies that interact with the mean flow, contain most
of the kinetic energy of the turbulence, and lie at one end of a wide range of
eddies that interact nonlinearly through mechanisms such as vortex stretching. It
was understood that in equilibrium the rate of working of the fluid against the
viscous stresses in the smallest of these eddies dissipates kinetic energy at the same
mean rate it is extracted from the mean flow by the large eddies. But Taylor (1935)
misidentified the spatial scale of this dissipative microturbulence, as he called it.
He chose the scale (Chapter 7)
1/2
νu2
λ= , (1.42)
but as we discussed in Section 1.6.4, Kolmogorov argued (and we now accept) that
the dissipative-eddy scale is η = (ν 3 / )1/4 . They are related by
λ 1/4
∼ Rt (1.43)
η
(Problem 1.19), so that λ can be considerably larger. We now call λ the Taylor
microscale.
According to Batchelor (1996), Taylor’s “The spectrum of turbulence” (1938)
was his “last paper on turbulence before the needs of World War II took him away
from academic research.” In this paper Taylor connected the power spectral density,
or spectrum, to the two-point correlation function in physical space through the
Fourier transform (Part III) and, in a casual aside, introduced what is now known
as “Taylor’s hypothesis” for interpreting a time series at a spatial point as a spatial
record in the upstream direction.
But until Kolmogorov (1941) there was no unifying view of turbulence dynamics
across the entire scale range. He postulated that for scales beyond the energy-
containing range there are but two governing parameters: the mean rate of transfer of
kinetic energy per unit mass from larger eddies to smaller (the energy cascade rate,
Part III) and the fluid kinematic viscosity. Therefore (as we showed in Section 1.6.4)
this extensive scale range yields readily to dimensional analysis. The dimensional
analysis is even simpler in what we now call the inertial subrange, the larger-
scale end of this range. Here the local Reynolds number is large so the viscosity
is not important, and so the turbulence spectrum – the mean-squared amplitude
of velocity fluctuation as a function of spatial wavenumber (inverse scale) – is
then determined solely by the energy cascade rate. Its analytical form, the famous
Kolmogorov spectrum (Chapter 7), emerges directly from dimensional analysis.
In large-Rt turbulence this inertial subrange can be extensive; in the atmospheric
boundary layer it can span four decades.
22 Introduction
1.6 What old saying about snowflakes applies also to turbulent flows?
Explain.
1.7 Explain why most practically important turbulent flows cannot be calculated
numerically through their governing equations.
1.8 Explain how the Reynolds stress emerges when the equation of motion for
turbulent flow is averaged. Interpret the Reynolds stress physically.
1.9 What is a turbulence model? Why are they necessary? What are the two
broad types, and how do they differ?
1.10 Why do you think only the lowest part of the atmosphere is continuously
turbulent?
1.11 Explain the statement “turbulence is an unsolved problem.”
1.12 We say turbulence is random and stochastic. What do those terms mean?
1.13 Explain how the rate of viscous dissipation can be independent of the fluid
viscosity when that viscosity appears in its definition.
1.14 Explain why the density of a mass-conserving constituent in a turbulent flow
need not be a conserved variable.
1.15 Discuss what properties of the ordinary derivative are shared by the
substantial derivative. What property is not shared?
Problems
1.1 Consider steady, fully developed laminar flow in a circular pipe of diameter
D. The axial equation of motion is
1 ∂p(x) ν ∂ ∂u(r)
0=− + r .
ρ ∂x r ∂r ∂r
(a) Why does p not depend on r? u not depend on x? ∂p/∂x not depend
on x?
(b) Solve this differential equation for u(r).
(c) Find expressions for uave , Eq. (1.2), and the Darcy friction factor f ,
Eq. (1.5).
(d) Express the required pumping power per unit mass of fluid in terms of f .
1.2 Generalize Problem 1.1 to include conduction heat transfer in the radial
direction. The temperature equation is
∂T (x, r) α ∂ ∂T (x, r)
u(r) = r .
∂x r ∂r ∂r
Here α = k/(ρcp ) is the thermal diffusivity. Consider the case when ∂T /∂x
does not depend on x.
24 Introduction
1.14 Two geometrically identical turbulent flows, one using water and the other
air, have the same u and . By what factor do their total rates of dissipation
of kinetic energy differ?
1.15 If Sij and Aij are symmetric and antisymmetric tensors, respectively, show
that their contraction Sij Aij vanishes.
1.16 Resolve the paradox that is independent of ν but the expression for it
(Problem 1.4) involves ν.
1.17 Explain why, in connection with Eq. (1.40), ui uj is different from ui uj in a
turbulent flow.
1.18 Derive the vorticity equation (1.28) from the Navier–Stokes equation (1.26).
Hint: use the vector identities ω × u = (u · ∇)u − ∇(u · u)/2 and (u · ∇)u =
ω × u + ∇(u · u)/2.
1.19 Derive Eq. (1.43) for the ratio of the Taylor and Kolmogorov microscales.
1.20 Do Problem 1.4 for steady flow in a pipe, integrating between two cross
sections. Interpret your result.
1.21 Show that the total time derivative need not commute with other derivatives.
References
Batchelor, G. K., 1950: The application of the similarity theory of turbulence to
atmospheric diffusion. Quart. J. Roy. Meteor. Soc., 76, 133–146.
Batchelor, G. K., 1967: An Introduction to Fluid Dynamics. Cambridge University Press.
Batchelor, G. K., 1996: The Life and Legacy of G. I. Taylor. Cambridge University Press.
Bister, M., and K. A. Emanuel, 1998: Dissipative heating and hurricane intensity. Meteor.
Atmos. Phys., 65, 233–240.
Businger, S., and J. A. Businger, 2001: Viscous dissipation of turbulence kinetic energy in
storms. J. Atmos. Sci., 58, 3793–3796.
Carusi, E., and A. Favaro, Eds., 1924: Leonardo da Vinci’s Del Moto e Misura dell’
Acqua. Bologna.
Corrsin, S., 1961: Turbulent flow. Am. Scientist, 49, 300–325.
Deardorff, J. W., 1970a: A numerical study of three-dimensional turbulent channel flow at
large Reynolds numbers. J. Fluid Mech., 41, 453–480.
Deardorff, J. W., 1970b: A three-dimensional numerical investigation of the idealized
planetary boundary layer. Geophys. Fluid Dyn., 1, 377–410.
Deardorff, J. W., and G. E. Willis, 1975: A parameterization of diffusion into the mixed
layer. J. Appl. Meteorol., 14, 1451–1458.
Dittus, F. W., and L. M. K. Boelter, 1930: Heat transfer in automobile radiators of the
tubular type. Univ. Calif., Berkeley, Publ. Eng., 2, 443–461.
Friedlander, S. K., and L. Topper, Eds., 1961: Turbulence: Classical Papers on Statistical
Theory. New York: Interscience.
Frisch, U., 1995: Turbulence: The Legacy of A. N. Kolmogorov. Cambridge University
Press.
Gleick, J., 1987: Chaos. New York: Viking Penguin.
Holmes, P., J. L. Lumley, and G. Berkooz, 1996: Turbulence, Coherent Structures,
Dynamical Systems and Symmetry. Cambridge University Press.
26 Introduction
• Homogeneous turbulence has spatially uniform statistical properties (with the excep-
tion of mean pressure). A turbulent flow can be homogeneous in zero, one, two, or
three directions. A sphere wake is an example of the first. The turbulent boundary
layer near the leading edge of a flat plate, or downstream of a change in surface con-
ditions, can be homogeneous in one direction (the lateral) but is inhomogeneous in
the wall-normal and streamwise directions. The turbulent boundary layer over a uni-
form surface can be homogeneous in two directions, those in the plane parallel to
the surface, but is necessarily inhomogeneous in the normal direction. The grid tur-
bulence produced by a grating of bars spanning the cross section of a wind tunnel is
homogeneous in the cross-stream plane but inhomogeneous in the streamwise direction
because it decays as it goes downstream (Chapter 5). A carefully tailored homogeneous
shear flow in a wind tunnel (Tavoularis and Corrsin, 1981) is, to a good approximation,
homogeneous in all three directions.
• Steady (also called stationary) turbulence – turbulence in a laboratory flow
driven by a constant-speed blower, for example – has statistics that are independent of
time. Turbulence in the atmospheric boundary layer can be steady for up to a few hours,
but near sunrise and sunset in clear weather it is made unsteady (nonstationary) by the
changing surface energy budget (Chapter 9).
27
28 Getting to know turbulence
Figure 2.1 The turbulent wake of a bullet, which is several hundred wake diame-
ters to the left. This shadowgraph shows the strikingly sharp, irregular boundary
between the turbulent fluid and the almost motionless fluid outside. Photograph
courtesy Army Research Laboratory. From Van Dyke (1982).
• Isotropic turbulence has statistical properties that are independent of translation, rotation,
and reflection of the coordinate axes. It must decay in time, for the production mechanisms
that maintain stationary turbulence are anisotropic (Chapter 6). Isotropy of the smallest
spatial scales is called local isotropy (Chapter 14).
• The logarithmic profile (or law of the wall) and the constant-stress layer refer to the height
variation of mean velocity and the Reynolds shear stress, respectively, in the turbulent
boundary layer over a flat surface. Like the Reynolds stresses themselves (which G. I.
Taylor (1935) called “virtual mean stresses”) they exist only as averages.
• The Gaussian plume refers to a mean effluent plume in homogeneous turbulent flow, not
an instantaneous plume.
• A well-mixed state of a convective atmospheric boundary layer (Chapter 11) has mean,
not instantaneous, profiles of potential temperature and water-vapor mixing ratio that are
essentially uniform with height.
• A turbulent flux of a property is the mean, not instantaneous, amount of that property
flowing through unit area per unit time due to the turbulent velocity.
2.2 Averaging
All the dependent variables in a turbulent flow – velocity, vorticity, temperature (if
there is heat transfer), density of an advected constituent (if there is mass transfer),
pressure – are turbulent. At any instant they are distributed irregularly in space, at
any point in space they fluctuate chaotically in time, and at given point and a given
time they vary randomly from realization to realization. Since Osborne Reynolds’
2.2 Averaging 29
Figure 2.2 Upper: Profiles of temperature and ozone mixing ratio at the top of a
cloud-capped convective boundary layer. z is measured from the mean top. From
Lenschow et al. (1988). Lower: Airborne-lidar measurements of aerosol concen-
tration in a vertical plane of a clear convective boundary layer. The instantaneous
top is quite thin and variable in space. Courtesy C. Kiemle, DLR Germany, and
J. Grabon, Penn State.
time it has been traditional to separate a turbulent flow variable a(x, t) into mean
and fluctuating parts distinguished by the overbar and prime:
The prime can make this notation cumbersome, so we shall use instead the notation
of Tennekes and Lumley (1972), denoting the “full” turbulent flow variable with
30 Getting to know turbulence
a tilde, ã(x, t), and representing its mean and fluctuating parts with upper- and
lower-case symbols:
ã(x, t) = A(x, t) + a(x, t). (2.2)
Several types of averages have been used to define mean values in turbulence.
Reynolds (1895) used a volume average. Somewhat later (in the 1930s, accord-
ing to Monin and Yaglom (1971)), Kolmogorov and his school, and Kampé de
Fériet, brought the ensemble average of statistical physics to turbulence; it is
conceptually the most elegant. Tennekes and Lumley (1972) used a time aver-
age in steady conditions. A time average is almost always used with quasi-steady
observations, and space averages in homogeneous directions are convenient with
numerical-simulation results.
N
1
ã(x, t) ≡ A(x, t) ≡ lim ã(x, t; α). (2.3)
N →∞ N
α=1
As indicated in Eq. (2.3), the ensemble average can depend on both position and
time. Being linear, it commutes with other linear operations such as differentiation
and integration,
b b
∂ ã(x, t; α) ∂
= ã(x, t), ã(x, t; α) dt = ã(x, t) dt, (2.4)
∂t ∂t a a
• The average of a sum is the sum of the averages (the distributive property):
ã + b̃ = ã + b̃. (2.5)
† Reynolds (1895) used a volume average but assumed it follows the rules for the ensemble average, so in the
turbulence community they now bear his name. We’ll discuss the volume averaging rules in Chapter 3.
Discovering Diverse Content Through
Random Scribd Documents
— Moi je me couche avec les poules, mais pas les mêmes…
Vous les aimez, au moins ? on ne sait plus avec qui on a affaire.
Non, vous n’êtes pas de la bande ?
— Vous verrez bien.
— Tiens, oui.
Le visage de Françoise, bien que petit, n’était pas fin mais
semblait l’être, usé par une tendresse dévergondée. Et sous le
cotillon simple, d’un sans-façon affecté, un corps fluet, vif, aidé de
muscles minces, serrés.
— … mon quatrième fils…
Où diable a-t-elle pu les mettre ses quatre fils ? Drôle de petite
dégourdie. Ils se promenèrent dans le parc, elle sauta une barrière.
Elle était nue sous sa robe de flanelle, avec de longues chaussettes.
Mais son visage faisait rêver à l’entour une toilette bien plus
féminine, fraîche, vaporeuse, aux couleurs du matin.
— … Quand mon mari sera revenu de Paris…
— Vous aimez votre mari ?
— Bien sûr… Je l’ai aimé comme une bête pendant dix ans.
— Maintenant, c’est la onzième année.
— Dame oui… Quelquefois j’ai envie de m’en aller tout à fait de
l’autre côté de la terre avec un type tombé du ciel.
— Vous me montrerez vos chevaux ?
— Tout de suite, si vous voulez. Tenez, c’est cela. Laissons les
femmes et allons-y.
Ils revinrent en courant vers la maison et sautèrent dans la
voiture.
— Je vous le prends… Cinq minutes.
Son démarrage menaça un rosier et le coup de vent cassa leur
vague geste d’excuse.
— Comment vous appelez-vous ? lui demanda-t-elle un instant
après, avec un sourire d’une ironie mouillée et complice.
— Gille.
— Tiens, c’est drôle, c’est un peu niais, ça vous ressemble.
J’aime cela.
Elle se tenait droite à son volant, ses petites mains dans des
gants bien sales. Son vieux chapeau écrasant sur la légère
couperose de sa joue une mèche très blonde. Rien que ces artifices
campagnards : ni fard, ni poudre.
Gille sortit subrepticement du fond de son enfance le rêve d’une
châtelaine courageuse et pure.
Elle se jeta dans des chemins de traverse, parmi d’opulents
herbages. On toucha à un petit bois assez fourré.
— On va s’arrêter là. Il fait bon, on ira à pied à travers le taillis
jusqu’à ma ferme. Marchez derrière moi.
Aussitôt il y eut un jeu pour se protéger des branches l’un l’autre.
Ils se frôlèrent, se touchèrent ; leurs corps se heurtèrent, leurs mains
se pressèrent sur le même rameau.
Gille l’embrassa dans le cou, ce qui la renversa dans ses bras. Ils
mêlèrent aussitôt leurs bouches et leurs membres, tapis sous un
buisson. Leurs gestes étaient sûrs.
Après cela, que dire ? Et la châtelaine sur sa tour ? Gille ricanait
mais appréciait une bonne tenue. Elle le félicita avec des mots justes
de l’avoir contentée. Ses paroles, ses regards étaient infléchis par
cette douceur qu’il avait remarquée sur son visage.
Le début d’une aventure ouvrait une perspective naïve à Gille et
bien que tous les défauts de l’amoureuse le piquassent dès la
première minute, cette facilité d’illusion lui permettait d’enchanter
l’autre, un instant, comme lui-même. Encore avec Molly, il avait
passé toute la matinée à imaginer une prolongation assez
improbable de leur voisinage. Il la compara avec son nouveau
plaisir. Mais il l’avait beaucoup oubliée depuis quelques heures et il
ne retrouvait rien d’elle. Et pourtant son corps, plein de pulpe, lui
convenait mieux que celui de cette mère de famille émaciée, qui,
encore enfantin à trente-huit ans, faisait songer à une fillette meurtrie
par un stupre prématuré.
Bras dessus, bras dessous, ils s’en allèrent comme des
camarades qui ont joué, jusqu’à la ferme où ils burent du cidre. Il
arriva à Françoise, au moment de remonter dans la torpedo blanche
qui avait attendu dans l’herbe, narquoise, distillant légèrement son
essence parmi les fleurs des champs, de dire : « Je suis une grande
garce, tiens, pourtant j’avais un joli petit cœur », avec un rire un peu
fatigué, à moitié rentré et des yeux qui pétillaient d’astuce sous une
légère buée.
Comme ils rejoignaient la grande route, ils se trouvèrent nez à
nez avec la voiture de Finette qui, avec Luc et Molly, revenait de
raccompagner Lady Hyacinthia.
— Eh bien ! je vous le rends ; s’écria Françoise avec un regain de
bonne humeur. Il est gentil, vous savez, soignez-le. Et elle disparut.
Gille vint s’asseoir modestement entre Finette, confite dans la
plus hypocrite quiétude et Molly pincée, mais qui, pour ne pas mentir
à ses opinions libérales, parla d’autre chose, de façon assez fluente.
Le dîner remit tout le liant souhaitable entre eux, grâce à des
coquetailles et à la verve de Luc qui tenait l’excellent sujet de
l’Anglaise. Le rire eut vite détendu Molly, ce qui encouragea Gille à
lui jeter quelques regards sournois d’enfant prodigue. Elle y
répondait comme quelqu’un qui est décidé à tuer le veau gras le soir
même, ce qui ne laissa pas d’effrayer le jeune coureur. Mais il était
fort animé et suivant son penchant il ne put résister à l’envie de
ressaisir ce qu’il avait lâché. Pourtant après ce repas il revint à
Finette.