Fluid Waves 1st Edition Manasseh Download PDF
Fluid Waves 1st Edition Manasseh Download PDF
Fluid Waves 1st Edition Manasseh Download PDF
com
https://ebookmeta.com/product/fluid-waves-1st-
edition-manasseh/
OR CLICK BUTTON
DOWLOAD EBOOK
https://ebookmeta.com/product/breaking-waves-1st-edition-barry-
litherland/
https://ebookmeta.com/product/waves-and-optics-1st-edition-
harish-parthasarathy/
https://ebookmeta.com/product/waves-of-global-terrorism-1st-
edition-david-rapoport/
https://ebookmeta.com/product/waves-1st-edition-julie-anne-
addicott-ashley-lane/
Understanding Gravitational Waves (Astronomers'
Universe) 1st Edition Kitchin
https://ebookmeta.com/product/understanding-gravitational-waves-
astronomers-universe-1st-edition-kitchin/
https://ebookmeta.com/product/theory-of-elastic-waves-1st-
edition-peijun-wei/
https://ebookmeta.com/product/sound-waves-and-acoustic-
emission-1st-edition-claudia-barile/
https://ebookmeta.com/product/nonlinear-dispersive-waves-and-
fluids-1st-edition-shijun-zheng/
https://ebookmeta.com/product/in-the-waves-lennon-stella-maisy-
stella/
Fluid Waves
Fluid Waves
Richard Manasseh
First edition published 2022
by CRC Press
6000 Broken Sound Parkway NW, Suite 300, Boca Raton, FL 33487-2742
Reasonable efforts have been made to publish reliable data and information, but the author and
publisher cannot assume responsibility for the validity of all materials or the consequences of their use.
The authors and publishers have attempted to trace the copyright holders of all material reproduced
in this publication and apologize to copyright holders if permission to publish in this form has not
been obtained. If any copyright material has not been acknowledged please write and let us know so
we may rectify in any future reprint.
Except as permitted under U.S. Copyright Law, no part of this book may be reprinted, reproduced,
transmitted, or utilized in any form by any electronic, mechanical, or other means, now known or
hereafter invented, including photocopying, microfilming, and recording, or in any information
storage or retrieval system, without written permission from the publishers.
For permission to photocopy or use material electronically from this work, access www.copyright.com
or contact the Copyright Clearance Center, Inc. (CCC), 222 Rosewood Drive, Danvers, MA 01923,
978-750-8400. For works that are not available on CCC please contact mpkbookspermissions@tandf.
co.uk
Trademark notice: Product or corporate names may be trademarks or registered trademarks and are
used only for identification and explanation without intent to infringe.
DOI: 10.1201/9780429295263
Typeset in CMR10
by KnowledgeWorks Global Ltd.
To Irena and Dylan
Contents
Nomenclature . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . xi
Preface xvii
vi
Contents vii
2 Water-surface waves 47
2.1 Summary of key points . . . . . . . . . . . . . . . . . . . . . . 47
2.2 An example . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
2.3 Linear water-wave theory . . . . . . . . . . . . . . . . . . . . 51
2.3.1 The waves we see . . . . . . . . . . . . . . . . . . . . . 51
2.3.2 Potential flow . . . . . . . . . . . . . . . . . . . . . . . 52
2.3.2.1 Physical assumptions that lead to potential
flow . . . . . . . . . . . . . . . . . . . . . . . 52
2.3.3 Laplace’s equation . . . . . . . . . . . . . . . . . . . . 56
2.3.4 Boundary conditions for water waves . . . . . . . . . . 56
2.3.5 Airy’s solution for surface gravity waves . . . . . . . . 59
2.3.5.1 Separation of variables solution . . . . . . . . 59
2.3.5.2 Applying the boundary conditions . . . . . . 60
2.3.5.3 Dispersion relation . . . . . . . . . . . . . . . 61
2.3.5.4 Ripples . . . . . . . . . . . . . . . . . . . . . 62
2.3.5.5 Phase speed . . . . . . . . . . . . . . . . . . 63
2.3.5.6 Velocity field . . . . . . . . . . . . . . . . . . 65
2.3.6 Surface elevation . . . . . . . . . . . . . . . . . . . . . 66
2.3.7 Particle trajectories . . . . . . . . . . . . . . . . . . . 69
2.3.8 Group velocity . . . . . . . . . . . . . . . . . . . . . . 70
2.3.9 Deep-water approximation . . . . . . . . . . . . . . . . 75
2.3.10 Consequences of deep water . . . . . . . . . . . . . . . 77
2.3.10.1 Maximum wavelength of ocean swell . . . . . 77
2.3.10.2 V-shaped wakes in deep water . . . . . . . . 77
2.3.10.3 Deep-water wave focusing . . . . . . . . . . . 80
2.3.11 Shallow-water approximation . . . . . . . . . . . . . . 80
2.3.12 Consequences of shallow water . . . . . . . . . . . . . 82
2.4 Problems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
3 Sound waves 89
3.1 Summary of key points . . . . . . . . . . . . . . . . . . . . . . 89
3.2 An example . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
3.3 Linear sound-wave theory . . . . . . . . . . . . . . . . . . . . 92
3.3.1 Use and control of sound . . . . . . . . . . . . . . . . 92
3.3.2 The wave equation for sound waves . . . . . . . . . . . 94
3.3.3 Solution of the wave equation . . . . . . . . . . . . . . 98
3.3.4 Relation to shallow-water waves . . . . . . . . . . . . 100
3.3.5 Refraction . . . . . . . . . . . . . . . . . . . . . . . . . 101
3.3.6 Acoustic impedance . . . . . . . . . . . . . . . . . . . 102
3.3.7 Reflection, scattering and transmission . . . . . . . . . 103
3.3.8 Representation and measurement of sound . . . . . . . 105
3.3.8.1 Spectral representation of sound . . . . . . . 105
3.3.8.2 Sound-measurement instruments . . . . . . . 107
3.3.9 Geometrical spreading . . . . . . . . . . . . . . . . . . 110
3.3.10 Doppler effect . . . . . . . . . . . . . . . . . . . . . . . 111
viii Contents
Bibliography 275
Index 287
Nomenclature xi
Nomenclature
Number type R: real; C: complex. MKS: metre-kilogram-second units; SI: standard
Système International units. To save space with multiple units, slashes are used
instead of the power notation used elsewhere in this book, so kg/m/s2 is kg m−1 s−2 .
Dashes (-) indicate dimensionless quantities. A few symbols, e.g. m and T , have two
meanings; if so, the meanings are quite different and should be apparent from the
context. Symbols only used briefly during a derivation or explanation are not below,
since their definition is on the same page or very close.
Roman letters
Type Units
MKS SI
n Integer number R -
N Buoyancy frequency R rad/s
N Integer number, usually maximum in series R -
P Pressure, total (or ‘absolute’) R kg/m/s2 Pa
P0 Pressure, total, in some initial steady state R kg/m/s2 Pa
P∅ Pressure, ambient R kg/m/s2 Pa
Patm Pressure, atmospheric R kg/m/s2 Pa
Ph Pressure, hydrostatic R kg/m/s2 Pa
Ps Pressure, engineering static R kg/m/s2 Pa
P Function of spatial variables (e.g. x, y, z) only R kg/m/s2 Pa
P̄η̂ Power, cycle-averaged, per unit crest length of
of single-frequency wave R kg m/s3 W/m
P̄H Power, cycle-averaged, per unit crest length of
of ocean-wave spectrum R kg m/s3 W/m
P̄P Power, cycle-averaged, engineering system R kg m2 /s3 W
p Pressure, dynamic C kg/m/s2 Pa
p̂ Pressure amplitude R kg/m/s2 Pa
pR Real part of p R kg/m/s2 Pa
p∅ Pressure, dynamic, reference (acoustics) R kg/m/s2 Pa
p∞ Pressure from external driver (acoustics) C kg/m/s2 Pa
pσ Pressure, surface tension R kg/m/s2 Pa
Q Pressure constant, usually a pressure scale R kg/m/s2 Pa
R Radius, time-varying, of an object R m
R0 Radius constant R m
R Function of r only R -
<(•) Real part of •
Re Reynolds number R -
r Distance in the radial direction R m
St Strouhal number R -
s Distance in an arbitrary direction R m
T Time constant, typically one period (cycle) R s
T Temperature, absolute R K K
T Function of t only C -
t Time R s
t Dimensionless time R -
xiv Nomenclature
Greek letters
α Angle R rad
αT Thermal diffusivity R m2 /s
β Angle R rad
γ Adiabatic index (ratio of specific heats) R -
∆ Prefix indicating a difference in a variable
δ Perturbation in bubble radius C m
δω Stokes-layer thickness R m
Dimensionless small parameter in derivations R -
ζ Damping ratio R -
ζP Damping ratio due to useful-power extraction R -
ζr Damping ratio due to wave radiation R -
ζT Damping ratio due to thermal losses (bubbles) R -
ζµ Damping ratio, linearised fluid-dynamical
(wave-energy conversion) or viscosity (bubbles) R -
η Surface elevation of surface gravity waves C m
ηR Real part of η R m
η̂ Amplitude of surface gravity waves; if waves
are nonlinear, defined as Ĥ/2 R m
θ Angle, or latitude of the Earth R rad
κ Wavenumber in x-direction (Cartesian);
in r-direction (spherical polars) R m−1
κp Polytropic index R -
Λ Wavenumber in r-direction (cylindrical polars) R m−1
λ Wavelength R m
µ Dynamic viscosity R kg/m/s Pa s
ν Kinematic viscosity R m2 /s
ξ Displacement of a fluid particle (real) or of
mechanical oscillator in a fluid (complex) C m
ρ Density R kg/m3
ρ0 Time-averaged, constant density (may vary R kg/m3
slowly in z)
% Small variation in density from ρ0 C kg/m3
σ Surface tension coefficient R kg/s2 N/m
τ Stress component C kg/m/s2 Pa
τ Stress tensor C kg/m/s2 Pa
Φ Phase angle R rad
φ Velocity potential C m2 /s
ϕ Azimuthal angle R rad
Ω Angular rotation rate R rad/s
Ω⊕ Angular rotation rate of the Earth R rad/s
ω Frequency R rad/s
ω0 Frequency, natural (undamped) R rad/s
ω0ζ Damped natural frequency R rad/s
xvi Nomenclature
Operators
xvii
xviii Preface
Structure
Part I addresses Aims 1 and 2. It comprises selected introductory fundamental
topics and results on fluid waves. These are topics thought to be of perpetual
relevance, such as surface gravity waves and sound waves. It is not compre-
hensive, omitting some significant topics such as waves due to instabilities in
fluid flows, many classes of geophysical waves affected by stratification and
rotation, and shock waves in compressible flows. Within each chapter, the
introductory nature of this book has also meant that many sub-topics are
omitted. The references provided cover the missing material. Meanwhile, a
few topics are subjected to full, if rather lengthy, mathematical derivations so
that there are few gaps left as ‘exercises for the student’. Thus, the reader is
made aware of the many steps required to address a typical problem in this
field and is fully equipped to tackle similar problems. Inevitably, however,
many detailed derivations are not undertaken.
Part II addresses Aim 3. It is intended to showcase the breadth of rele-
vance of wave phenomena in fluids. Thus, Part II is rather like a collection
of specialised review papers, which in most cases refer to sources right up to
the date of publication of this book. Some chapters involve significant deriva-
tions, but most are explanatory. In some cases, the applications, which are as
diverse as renewable energy, cancer treatment, climate change, and the search
for habitable exoplanets, are pursued to a depth such that the relevance to
fluid waves may appear tangential; nevertheless, without the ability to cal-
culate fluid-wave problems, rigorous pursuit of the applications would not be
possible.
P∅ = Patm + Ph + Ps ,
P = P∅ + p .
• The Ideal Gas Law, relating the absolute pressure P and the volume
per unit mass, V, to their initial values, P0 and V0 , is given by (1.5) on
page 10,
P V κp = P0 V0 κp ,
where κp is the polytropic index (κp = γ for an adiabatic or isentropic
process, where the adiabatic index γ is very close to 1.4 for air).
• The Viscosity Equation (Newton’s Law of Viscosity), (1.8) on page 13,
relates the shear stress τ to the gradient in velocity by
du
τ =µ ,
dy
DOI: 10.1201/9780429295263-1 3
4 Fundamentals
Du 1
= − ∇P + ν∇2 u + g, ,
Dt ρ0
ω
f= ,
2π
where f is the frequency in Hertz (or cycles per second) and ω is the fre-
quency in radians per second (often just called the ‘frequency’); secondly,
Summary of key points 5
UL
Reynolds Re Nonlinear-inertia to
ν viscous
U
Froude Fr √ Nonlinear-inertia to
gL gravity
L
Strouhal St Oscillatory- to
UT nonlinear-inertia
UT
Keulegan– KC Nonlinear-inertia to
Carpenter L oscillatory
p
Womersley Wo L √ω/ν Oscillatory-inertia to
or 2 L/δω viscous
U
Mach Ma Nonlinear-inertia to
c bulk-stiffness (flow
speed to sound speed)
6 Fundamentals
distinct surface bounding a finite volume, whereas a given mass of gas ex-
pands to fill whatever container it is housed in. Despite this difference, the
same fluid-dynamical laws - and almost all of the same fluid-wave phenomena
- apply equally to liquids and to gases.
Fluid mechanics is the study of forces and motions in fluids, and of the
forces and motions caused by fluids on solid objects and vice versa. Together
with solid mechanics, fluid mechanics forms the branch of classical physics
called continuum mechanics. The term fluid dynamics is often used inter-
changeably with ‘fluid mechanics’, but is usually reserved for those fluid me-
chanical situations where there is some motion. For example, the calculation
of the force on a dam wall due to the pressure of motionless water on the wall
is a fluid mechanics problem, but not a fluid dynamics problem. Since fluid
waves inherently involve motion, they are fluid-dynamical phenomena.
Continuum mechanics, and thus fluid mechanics and fluid dynamics, is
only valid in a continuum. This is a substance that is definitely composed
of indivisible ‘particles’, such as molecules or atoms, but whose behaviour
only matters to us for volumes very much larger than the individual particles.
Thus, it is possible to average over the very many particles inside the smallest
volume of interest to us and ignore the reality that the particles exist. Thus,
in fluid dynamics, individual molecules or atoms, which move in a statistical
fashion, confer their averaged quantities such as pressure, velocity, density or
temperature to continuum equations, allowing us to use the rules of calculus
to determine the behaviour of the continuum. It is worth noting that if the
smallest volume of interest is large enough, a continuum can also be considered
to be composed of ‘particles’ that are much larger than atoms or molecules.
For example, sand or grains of wheat can be observed to flow rather like a fluid,
and subject to suitable approximations, the flow of such granular materials
could be treated as a fluid dynamics problem. Meanwhile, the particles of
dust that compose the accretion disk of a forming solar system (discussed in
§11.2.1), even though there are great distances between the particles, can be
treated as a continuum when the smallest volume of interest is the size of a
planet, so there are still very many particles in the smallest volume.
Stress is force per unit area and has units of kg m−1 s−2 , with the SI name
Pascal, abbreviated Pa. It is divided into two types, normal stress and shear
stress. The normal stresses create the pressure in the fluid, though the normal
stresses are not uniquely related to the pressure unless the fluid is motionless
or incompressible. Stress, in general, is usually denoted with the symbol τ .
Since there are three dimensions of space (in Cartesian coordinates, the x,
y and z directions), forces could be acting in each of the three directions,
but the force in each direction could be acting on areas normal to all of the
x, y and z directions, giving nine possible combinations. Hence, τ is a 3 × 3
tensor. The three normal stresses are τxx , τyy and τzz , whereas the six shear
stresses are τxy , τxz , τyx , τyz , τzx and τzy . We will return to the shear stresses
in §1.2.2.4 when considering viscosity, but first, the normal stresses creating
pressure must be understood.
A bewildering variety of quantities are called ‘pressure’ in fluid dynam-
ics, engineering, physics and medicine. It is common for students to confuse
these and get the wrong result, so it is helpful to begin with a precise set of
definitions. The hydrostatic pressure is the pressure due to the weight of the
fluid above. If the fluid above is a liquid and thus almost incompressible, the
hydrostatic pressure Ph is given by
Ph = ρ0 gh , (1.1)
P∅ = Patm + Ph + Ps . (1.2)
In a motionless fluid, p = 0, so the total and ambient pressures are equal and
are equal to the average of the three normal stresses (and the normal stresses
are negative since the total pressure is positive), so P = P∅ = − 31 (τxx + τyy +
τzz ). If the fluid is moving but is incompressible, only the total pressure P is
given by this average, so P = − 13 (τxx + τyy + τzz ). If the fluid is moving and
is compressible, some part of the normal stresses is due to viscous resistance
to the rate of change of volume with time, which will be outlined in §1.2.2.4
below.
1.2.2.2 Strain
Linear strain is the ratio of one of the lengths of the substance when nor-
mal stress is applied along that length, to its length when not under stress.
Engineers typically study linear strain when learning how metal stretches
under load, for example, but for fluid flows, variations in the volume are
paramount.
Volumetric strain is the ratio of the volume, V , of the substance when
under some additional pressure to its volume when not under such pressure,
V∅ . If the contributions to normal stresses from fluid motion are all negative (so
that the dynamic pressure is positive), the substance is said to be compressed
by the motion. If these contributions are all positive (so that the dynamic
pressure is negative), the motion is said to put the substance under rarefaction
or sometimes ‘expansion’. Note that the ‘pressure’ we have just referred to is
the dynamic pressure p that is due to fluid motion, and thus the compression or
rarefaction is relative to the inherent compression due to the ambient pressure.
If a fluid is compressed, its density, ρ, will be higher than its ambient density
ρ0 , and if a fluid is rarefied, ρ will be lower than ρ0 . For example, in the
propagation of sound waves, ρ varies cyclically about ρ0 , being higher than
ρ0 during the compression half of the oscillation and less than ρ0 during the
rarefaction half of the oscillation.
Shear strain is the angle by which a substance is distorted owing to shear
stress.
10 Fundamentals
in which the volumetric strain can be identified once we realise that ρ/∆ρ =
∆V /V . It can be thought of as the ‘stiffness’ of the fluid since the higher
the bulk modulus, the more pressure is needed to cause a given volumetric
strain. In general, the bulk modulus is a function of temperature and fluid
composition. For liquids, the relation is often extracted from empirical data.
The bulk modulus will be used in §3.3.2 in the derivation of the speed of
sound, and in practice, the bulk modulus is actually determined from speed-
of-sound measurements (Fine and Millero, 1973). For an incompressible liquid,
which of course exists only as a theoretical approximation (albeit a very useful
approximation that we will employ many times), the bulk modulus is infinite.
For water Ev ' 2.2 × 109 Pa.
For gases, however, a much more convenient relation is available, and,
moreover, it can be derived from the fundamental theory of the mechanics of
gas molecules, called the kinetic theory of gases. The Ideal Gas Law relates
any value of total pressure, P , and the specific volume or volume per unit
mass, V, to any other ‘initial’ values of the total pressure and volume per unit
mass of the same gas, P0 and V0 . Here, care has been taken to introduce a
symbol for the initial state of the pressure, P0 , that is different to that of the
ambient pressure, P∅ , since some unchanging or ‘steady’ fluid flow might be
occurring (steady flows are defined in §1.2.4 below), causing some background
dynamic pressure; but if there is no such background flow, P0 = P∅ . Note that
V = 1/ρ and V0 = 1/ρ0 . The Ideal Gas Law is given by
P V κp = P0 V0 κp , (1.5)
where κp is the polytropic index whose value depends on the nature of the
heat transfer occurring during compression or expansion of the gas. If the
compression or expansion is isothermal, meaning that the temperature of the
gas does not alter during the compression or expansion, κp = 1. An isothermal
compression or expansion requires that heat flows out of or into the gas across
the boundaries of whatever is containing the gas without any restriction, so
that the temperature can stay constant. The bulk modulus for a gas is given
by Ev = κp P0 and, as just-noted, it will be used in §3.3.2 in the derivation
of the speed of sound. For isothermal air at atmospheric pressure and 20◦ C,
Ev ' 105 Pa.
If the compression or expansion is adiabatic or isentropic, meaning that
there is zero heat flow across the boundaries, κp = γ, where γ is the adiabatic
index. Adiabatic conditions need not imply the gas is in some container, like
Basic fluid dynamics 11
FIGURE 1.1
A set of oceanographic probes and sample bottles is winched into the ocean
from a research ship to a depth of over 1000 m. At such depths, water cannot
be considered incompressible, so the bulk modulus (1.4) is needed. For seawa-
ter, the bulk-modulus relation between stress and volumetric strain requires
an empirical relation, determined in practice by speed-of-sound measurements;
see Fine and Millero (1973) or Gill (1982). Photograph by Richard Manasseh.
the cylinder of a bicycle pump. If the expansions and contractions are very
rapid and are very small in magnitude, which is typical of most sound waves,
there is insufficient time for significant heat to flow one way on the crest of the
wave before it has to flow the opposite way on the trough, so in practice, the
heat flow is extremely small, if not exactly zero. The two cases of isothermal
and adiabatic behaviour represent the two idealised limiting cases for heat
transfer, so that 1 < κp < γ.
Now, the adiabatic index is related to the molecular structure of the gas
and is given by the formula
2
γ =1+ , (1.6)
nf
12 Fundamentals
∆u∆t
Infinitesimally
small element
∆θ
∆y under shear stress,
at time t + ∆t
FIGURE 1.2
The definition of the rate of shear strain
where nf is the number of degrees of freedom of the gas molecule. This is the
number of ways the atoms in the molecule can move and also move relative to
each other. In moving, the molecule possesses kinetic energy and thus heat.
In principle, the number of degrees of freedom of the molecule is purely a
geometric property of the molecule, but in practice, some degrees of freedom
only occur at high temperatures. For a gas such as argon (Ar) consisting of just
one atom, nf = 3 since the atom can move in the three dimensions of space.
For a diatomic molecule, for example, nitrogen (N2 ) or oxygen (O2 ), two added
degrees of freedom are active at room temperature, so nf = 3 + 2 = 5. Hence,
for both N2 and O2 , γ = 1 + 2/5 = 1.4, and since the Earth’s atmosphere
is about 78% N2 and 21% O2 , the air has a value of γ very close to 1.4.
Meanwhile, carbon dioxide (CO2 ) is a molecule with three atoms arranged in
a line. It has a value of γ of about 1.3 at room temperature. The adiabatic
index is also is equal to the ratio of specific heats of the gas, CP /CV , sometimes
called the ‘heat capacity ratio’.
1.2.2.4 Relation between shear stress and the rate of shear strain
If a fluid is flowing in the x-direction with speed u, and this speed u varies
in the y-direction (at right angles to the x-direction), the fluid is in a state of
shear strain. Imagine that u is increasing in the y-direction. The shear strain is
then constantly increasing with time as the faster fluid at higher y constantly
overtakes the slower fluid. The ability of a fluid to continually suffer an ever-
increasing shear strain is the attribute that defines a fluid; as noted in §1.2.1,
it deforms indefinitely in response to shear stress.
The definition of the rate of change of shear strain can be understood with
the aid of figure 1.2. The rate of change of shear strain is the rate at which
Basic fluid dynamics 13
the ‘tilt angle’, ∆θ, changes with time, i.e. ∆θ/∆t, where for small ∆θ,
∆u∆t
∆θ ' tan(∆θ) =
∆y
∆θ ∆u
⇒ ' (1.7)
∆t ∆y
Therefore, when the usual limits are taken as ∆y and ∆u tend to zero, the
rate of change of shear strain becomes du/dy.
However, the fluid does not experience stress without resistance. There
are chemical and physio-chemical forces between the particles of the fluid; for
example, there are attractive forces between the molecules of H2 0 in flowing
water, which transfer momentum from the faster-flowing portions of fluid to
the slower portions, slowing down the faster flows. Similarly, as faster-flowing
portions of gas overtake the slower portions, the molecules may collide, like-
wise transferring momentum. Owing to the random orientation of molecules,
these microscopic interactions inevitably transfer energy from the kinetic en-
ergy of the bulk flow to tiny and random molecular motions, and since the
amount of vibration of atoms and molecules is related to the temperature,
the fluid is thus heated. If the fluid is to continue deforming indefinitely, there
must be a shear stress continually applied to it that balances the shear stress
due to these inter-particle interactions.
The bulk fluid property that controls this transfer of momentum across
fluid layers is called dynamic viscosity (or sometimes ‘molecular viscosity’),
and is denoted with the symbol µ. The relation between the rate of shear
strain, du/dy, and the shear stress due to viscosity, τ , was first understood
by Isaac Newton; it is
du
τ =µ , (1.8)
dy
which is called the Viscosity Equation and traditionally called Newton’s Law
of Viscosity. For water, µ = 1.0016 × 10−3 Pa s at 20◦ C, and for air, µ =
1.81 × 10−5 Pa s at 20◦ C (Streeter and Wylie, 1979).
Now, a fluid flowing only in the x-direction with speed u could have gra-
dients in u in the z-direction as well as the y-direction, and indeed this is
occurring as fluid flows in any conduit (pipe, duct or channel); whatever the
value of the speed u in the middle of the conduit, u must be zero for those
fluid molecules physio-chemically bound to the walls. Thus the τ in (1.8) is
more properly written τyx , the shear stress in the x-direction acting on a plane
normal to the y-direction. Including the other directions gives
∂u
τyx = µ ,
∂y
∂u
τzx =µ . (1.9)
∂z
14 Fundamentals
FIGURE 1.3
A thin film of water flows down a glass feature window, forming waves affected
by surface tension. Photograph by Richard Manasseh.
Furthermore, the speed u could also vary in the x-direction, and as the flow
tries to compress (or expand) the fluid in the x-direction, energy could also be
transferred to heat, giving a form of normal stress in addition to that due to
pressure, the volume-viscous stress, with the corresponding parameter being
called the dilational viscosity or sometimes the ‘second coefficient of viscosity’.
The fluid could also be flowing in any and all of the x-, y- and z-directions
with velocity given by the vector u, where u = (u, v, w). Once flows in the
other directions are possible, for small strains, considering the geometry of
the deformed fluid element, (1.9) becomes
1 ∂u ∂v
τyx = µ + ,
2 ∂y ∂x
1 ∂u ∂w
τzx = µ + , (1.10)
2 ∂z ∂x
and similarly for τxy , τzy , τxz and τyz . This gives a total of six shear stresses.
Plus, in addition to the variation of u in the x-direction, possible variations
of v in the y-direction and w in the z-direction add to the volume viscous
stress. Thus, as noted in §1.2.2.1, there are nine stresses in the stress tensor,
τ , comprising the parts of the three normal stresses due to pressure plus the
three normal stresses due to volume viscosity, and the six shear stresses due
to Newtonian viscosity.
liquid together, whereas, for gas, the molecules or atoms move comparatively
freely relative to each other and constantly bounce off each other. At the liq-
uid surface, the attractive force is defined to have a value of σ Newton per
metre, so that an imaginary line of length L ‘drawn’ on the surface would have
a force of σL Newtons acting at right angles to it. For a water-air surface,
σ = 0.07275 ± 0.00036 N m−1 at 20◦ C (Vargaftik et al., 1983), or, slightly
more precisely, σ = 0.07268 ± 0.00018 N m−1 at 22◦ C (Berry et al., 2015),
and decreases with temperature (see Vargaftik et al., 1983, for values at var-
ious temperatures). If the surface is curved, it can be shown in a few lines
(Batchelor, 1973) that these forces develop a component normal to the sur-
face, causing an increase in pressure, pσ , on the concave side of the curved
surface given by
1 1
pσ = σ + , (1.11)
r1 r2
where r1 and r2 are the radii of curvature of the surface in the two planes at
right angles to the surface. Clearly, if length scales are small, so r1 , or r2 , or
both are small, surface tension becomes important; conversely, in flows with
large length scales, surface tension is negligible. Examples of such cases are in
§2.3.5.4 and §9.3.1.
Infinitesimally
u − ∆u/2 small element, u + ∆u/2
Inflow - cross-sectional - Outflow
ρ − ∆ρ/2 area ρ + ∆ρ/2
A×
∆x -
FIGURE 1.4
Transport of mass in one dimension for an infinitesimally small element of
fluid.
that mass. However, the principle of mass and energy equivalence, elucidated
by Albert Einstein in 1905, and encapsulated by the relation E = mc2 where
E is energy, m is mass and c is the speed of light, effectively unites the Law
of Conservation of Mass with the Law of Conservation of Energy. Excluding
nuclear reactions, however, it is convenient to express the conservation of mass
and energy with separate laws.
The conservation of mass in a fluid-dynamics context is easily expressed
by first imagining a one-dimensional flow in which fluid with density ρ is
transported by a flow with speed u along a channel in the x-direction, as
shown in figure 1.4. The channel has a cross-sectional area A× , say. After
travelling an infinitesimally small distance along the channel, ∆x say, imagine
that the speed has changed from u − ∆u/2 at the entry to the infinitesimally
small element, to u + ∆u/2 at the exit, and the density has likewise changed
from ρ − ∆ρ/2 to ρ + ∆ρ/2. The volume of this short segment is A× ∆x.
Therefore, mass is exiting this short segment of the channel of length ∆x at
a rate (ρ + ∆ρ/2)(u + ∆u/2)A× while it is entering at a rate (ρ − ∆ρ/2)(u −
∆u/2)A× , and the net change of mass in transiting the short segment would
be [(ρ − ∆ρ/2)(u − ∆u/2) − (ρ + ∆ρ/2)(u + ∆u/2)]A× .
Since ∆ρ and ∆u are by definition small, their product is negligible, so
the net change of mass becomes −(∆ρu + ρ∆u)A× . In general, ∆ρ and ∆u
could be positive or negative. Meanwhile, it is possible that the density inside
the short segment is changing with time, so the rate of net increase in mass
inside the short segment is (∆ρ/∆t)A× ∆x. Since mass can neither be created
nor destroyed, it can be neither gained nor lost in its journey along the short
segment of the channel, so that
Dividing (1.12) by the volume A× ∆x of the small segment and taking the
limits as the small quantities tend to zero gives
∂ρ ∂ρ ∂u
=− u−ρ ,
∂t ∂x ∂x
∂ρ ∂
⇒ = − (ρu). (1.13)
∂t ∂x
This reasoning is an application of the Reynolds Transport Theorem to mass.
Now, imagine that rather than being a segment of a channel, the small
volume is in the middle of a flowing mass of fluid, and, adopting Cartesian
coordinates, imagine that mass is able to flow into and out of the volume in
the y and z directions as well, with velocity components u, v and w in the
x, y and z directions respectively. The analysis is the same: dividing by the
volume ∆x∆y∆z at the equivalent of the step that led to (1.13) gives
∂ρ ∂ ∂ ∂
= − (ρu) − (ρv) − (ρw), (1.14)
∂t ∂x ∂y ∂z
or, using standard vector-calculus notation, use the operator ∇ for the deriva-
tives in space, which has the advantage that the coordinate system does not
need to be specified. Then (1.14) has the compact form
∂ρ
= −∇ · (ρu). (1.15)
∂t
F = ma, (1.16)
where F is the force vector, m is mass, and a is the acceleration vector. The
right-hand side of (1.16) is the rate of change of momentum with time. To
prepare for fluid-dynamical calculations, first, make a trivial re-arrangement
of (1.16),
F
a= . (1.17)
m
The purpose of this re-arrangement is to highlight that in fluid dynamics, the
greatest challenge is often posed by the determination of the motion of the
fluid for a given applied force. The calculation of the velocity resulting from
the acceleration can be extremely difficult, as we will see shortly.
18 Fundamentals
Infinitesimally
u − 12 ∆u small element, u + 12 ∆u
Inflow - cross-sectional - Outflow
1 1
ρu − 2
∆(ρu) area ρu + 2
∆(ρu)
A×
∆x -
FIGURE 1.5
Transport of momentum in one dimension for an infinitesimally small element
of fluid.
1.D. The copyright laws of the place where you are located also
govern what you can do with this work. Copyright laws in most
countries are in a constant state of change. If you are outside the
United States, check the laws of your country in addition to the terms
of this agreement before downloading, copying, displaying,
performing, distributing or creating derivative works based on this
work or any other Project Gutenberg™ work. The Foundation makes
no representations concerning the copyright status of any work in
any country other than the United States.
• You pay a royalty fee of 20% of the gross profits you derive from
the use of Project Gutenberg™ works calculated using the
method you already use to calculate your applicable taxes. The
fee is owed to the owner of the Project Gutenberg™ trademark,
but he has agreed to donate royalties under this paragraph to
the Project Gutenberg Literary Archive Foundation. Royalty
payments must be paid within 60 days following each date on
which you prepare (or are legally required to prepare) your
periodic tax returns. Royalty payments should be clearly marked
as such and sent to the Project Gutenberg Literary Archive
Foundation at the address specified in Section 4, “Information
about donations to the Project Gutenberg Literary Archive
Foundation.”
• You comply with all other terms of this agreement for free
distribution of Project Gutenberg™ works.
1.F.
1.F.4. Except for the limited right of replacement or refund set forth in
paragraph 1.F.3, this work is provided to you ‘AS-IS’, WITH NO
OTHER WARRANTIES OF ANY KIND, EXPRESS OR IMPLIED,
INCLUDING BUT NOT LIMITED TO WARRANTIES OF
MERCHANTABILITY OR FITNESS FOR ANY PURPOSE.
Please check the Project Gutenberg web pages for current donation
methods and addresses. Donations are accepted in a number of
other ways including checks, online payments and credit card
donations. To donate, please visit: www.gutenberg.org/donate.
Most people start at our website which has the main PG search
facility: www.gutenberg.org.