This preview shows page 1. Sign up to view the full content.
Unformatted text preview: Inviscid potential flow past an array of cylinders. The mathematics of potential theory, presented in this chapter, is both beautiful and manageable, but results may be unrealistic when
there are solid boundaries. See Figure 8.13b for the real (viscous) flow pattern. (Courtesy of
Tecquipment Ltd., Nottingham, England)  v v 214  eText Main Menu  Textbook Table of Contents  Study Guide Chapter 4
Differential Relations
for a Fluid Particle
Motivation. In analyzing fluid motion, we might take one of two paths: (1) seeking
an estimate of gross effects (mass flow, induced force, energy change) over a finite region or control volume or (2) seeking the pointbypoint details of a flow pattern by
analyzing an infinitesimal region of the flow. The former or grossaverage viewpoint
was the subject of Chap. 3.
This chapter treats the second in our trio of techniques for analyzing fluid motion,
smallscale, or differential, analysis. That is, we apply our four basic conservation laws
to an infinitesimally small control volume or, alternately, to an infinitesimal fluid system. In either case the results yield the basic differential equations of fluid motion. Appropriate boundary conditions are also developed.
In their most basic form, these differential equations of motion are quite difficult to
solve, and very little is known about their general mathematical properties. However,
certain things can be done which have great educational value. First, e.g., as shown in
Chap. 5, the equations (even if unsolved) reveal the basic dimensionless parameters
which govern fluid motion. Second, as shown in Chap. 6, a great number of useful solutions can be found if one makes two simplifying assumptions: (1) steady flow and
(2) incompressible flow. A third and rather drastic simplification, frictionless flow,
makes our old friend the Bernoulli equation valid and yields a wide variety of idealized, or perfectfluid, possible solutions. These idealized flows are treated in Chap. 8,
and we must be careful to ascertain whether such solutions are in fact realistic when
compared with actual fluid motion. Finally, even the difficult general differential equations now yield to the approximating technique known as numerical analysis, whereby
the derivatives are simulated by algebraic relations between a finite number of grid
points in the flow field which are then solved on a digital computer. Reference 1 is an
example of a textbook devoted entirely to numerical analysis of fluid motion. 4.1 The Acceleration Field
of a Fluid In Sec. 1.5 we established the cartesian vector form of a velocity field which varies in
space and time:
V(r, t) iu(x, y, z, t) j (x, y, z, t) kw(x, y, z, t) (1.4)  v v 215  eText Main Menu  Textbook Table of Contents  Study Guide Chapter 4 Differential Relations for a Fluid Particle This is the most important variable in fluid mechanics: Knowledge of the velocity vector field is nearly equivalent to solving a fluidflow problem. Our coordinates are fixed
in space, and we observe the fluid as it passes by—as if we had scribed a set of coordinate lines on a glass window in a wind tunnel. This is the eulerian frame of reference, as opposed to the lagrangian frame, which follows the moving position of individual particles.
To write Newton’s second law for an infinitesimal fluid system, we need to calculate the acceleration vector field a of the flow. Thus we compute the total time derivative of the velocity vector:
dV
dt a i du
dt j d
dt k dw
dt Since each scalar component (u, , w) is a function of the four variables (x, y, z, t), we
use the chain rule to obtain each scalar time derivative. For example,
du(x, y, z, t)
dt u
t u dx
x dt u dy
y dt u dz
z dt But, by definition, dx/dt is the local velocity component u, and dy/dt
w. The total derivative of u may thus be written in the compact form
du
dt u
t u
x u u
y u
t u
z w (V , and dz/dt )u (4.1) Exactly similar expressions, with u replaced by or w, hold for d /dt or dw/dt. Summing these into a vector, we obtain the total acceleration:
a dV
dt V
t V
x u Local V
y w V
z V
t (V )V (4.2) Convective The term V/ t is called the local acceleration, which vanishes if the flow is steady,
i.e., independent of time. The three terms in parentheses are called the convective acceleration, which arises when the particle moves through regions of spatially varying
velocity, as in a nozzle or diffuser. Flows which are nominally “steady” may have large
accelerations due to the convective terms.
Note our use of the compact dot product involving V and the gradient operator :
u x y w V z where i x j y k z The total time derivative—sometimes called the substantial or material derivative—
concept may be applied to any variable, e.g., the pressure:
dp
dt p
t u p
x p
y w p
z p
t (V )p (4.3) Wherever convective effects occur in the basic laws involving mass, momentum, or energy, the basic differential equations become nonlinear and are usually more complicated than flows which do not involve convective changes.
We emphasize that this total time derivative follows a particle of fixed identity, making it convenient for expressing laws of particle mechanics in the eulerian fluidfield  v v 216  eText Main Menu  Textbook Table of Contents  Study Guide 4.2 The Differential Equation of Mass Conservation 217 description. The operator d/dt is sometimes assigned a special symbol such as D/Dt as
a further reminder that it contains four terms and follows a fixed particle.
EXAMPLE 4.1
Given the eulerian velocityvector field
xzj 3ti ty2k xz V w find the acceleration of a particle. Solution
First note the specific given components
u 3t ty2 Then evaluate the vector derivatives required for Eq. (4.2)
V
t i V
x u
t zj j
V
y w
t k t 3i
V
z 2tyk y2k xj This could have been worse: There are only five terms in all, whereas there could have been as
many as twelve. Substitute directly into Eq. (4.2):
dV
dt (3i y2k) (3t)(zj) (xz)(2tyk) (ty2)(xj) Collect terms for the final result
dV
dt 3i (3tz txy2)j (2xyzt y2)k Ans. Assuming that V is valid everywhere as given, this acceleration applies to all positions and times
within the flow field. 4.2 The Differential Equation
of Mass Conservation All the basic differential equations can be derived by considering either an elemental
control volume or an elemental system. Here we choose an infinitesimal fixed control
volume (dx, dy, dz), as in Fig. 4.1, and use our basic controlvolume relations from
Chap. 3. The flow through each side of the element is approximately onedimensional,
and so the appropriate massconservation relation to use here is CV t d ( iAiVi)out ( iAiVi)in i 0 (3.22) i The element is so small that the volume integral simply reduces to a differential term  v v CV  eText Main Menu  t d Textbook Table of Contents t dx dy dz  Study Guide 218 Chapter 4 Differential Relations for a Fluid Particle
y
Control volume ρ u + ∂ (ρ u) d x dy dz
∂x ρ u dy dz
dy x
dz Fig. 4.1 Elemental cartesian fixed
control volume showing the inlet
and outlet mass flows on the x
faces. dx
z The massflow terms occur on all six faces, three inlets and three outlets. We make use
of the field or continuum concept from Chap. 1, where all fluid properties are considered to be uniformly varying functions of time and position, such as
(x, y, z, t).
Thus, if T is the temperature on the left face of the element in Fig. 4.1, the right face will
have a slightly different temperature T ( T/ x) dx. For mass conservation, if u is
known on the left face, the value of this product on the right face is u ( u/ x) dx.
Figure 4.1 shows only the mass flows on the x or left and right faces. The flows on
the y (bottom and top) and the z (back and front) faces have been omitted to avoid cluttering up the drawing. We can list all these six flows as follows: Face Inlet mass flow x u dy dz y dx dz z Outlet mass flow w dx dy u x
y w z ( u) dx dy dz
( ) dy dx dz
( w) dz dx dy Introduce these terms into Eq. (3.22) above and we have
t dx dy dz x ( u) dx dy dz y ( ) dx dy dz z ( w) dx dy dz 0 The element volume cancels out of all terms, leaving a partial differential equation involving the derivatives of density and velocity
t x ( u) y () z ( w) 0 (4.4)  v v This is the desired result: conservation of mass for an infinitesimal control volume. It
is often called the equation of continuity because it requires no assumptions except that
the density and velocity are continuum functions. That is, the flow may be either steady  eText Main Menu  Textbook Table of Contents  Study Guide 4.2 The Differential Equation of Mass Conservation 219 or unsteady, viscous or frictionless, compressible or incompressible.1 However, the
equation does not allow for any source or sink singularities within the element.
The vectorgradient operator
i j x k y z enables us to rewrite the equation of continuity in a compact form, not that it helps
much in finding a solution. The last three terms of Eq. (4.4) are equivalent to the divergence of the vector V
x ( u) () y z ( w) ( V) (4.5) so that the compact form of the continuity relation is
( V) t 0 (4.6) In this vector form the equation is still quite general and can readily be converted to
other than cartesian coordinate systems. Cylindrical Polar Coordinates The most common alternative to the cartesian system is the cylindrical polar coordinate system, sketched in Fig. 4.2. An arbitrary point P is defined by a distance z along
the axis, a radial distance r from the axis, and a rotation angle about the axis. The
three independent velocity components are an axial velocity z, a radial velocity r, and
a circumferential velocity , which is positive counterclockwise, i.e., in the direction
1 One case where Eq. (4.4) might need special care is twophase flow, where the density is discontinuous between the phases. For further details on this case, see, e.g., Ref. 2. υr υθ Typical point (r, θ , z) r θ υz
Base
line Typical
infinitesimal
element dr
dz r dθ Cyl ind rica l ax is  v v Fig. 4.2 Definition sketch for the
cylindrical coordinate system.  eText Main Menu z  Textbook Table of Contents  Study Guide 220 Chapter 4 Differential Relations for a Fluid Particle of increasing . In general, all components, as well as pressure and density and other
fluid properties, are continuous functions of r, , z, and t.
The divergence of any vector function A(r, , z, t) is found by making the transformation of coordinates
(x2 r y2)1/2 tan 1 y
x z z (4.7) (Az) (4.8) and the result is given here without proof2
1
r A 1
r (rAr) r (A ) z The general continuity equation (4.6) in cylindrical polar coordinates is thus
1
r t r (r 1
r r) ( ) ( z z) 0 (4.9) There are other orthogonal curvilinear coordinate systems, notably spherical polar coordinates, which occasionally merit use in a fluidmechanics problem. We shall not
treat these systems here except in Prob. 4.12.
There are also other ways to derive the basic continuity equation (4.6) which are interesting and instructive. Ask your instructor about these alternate approaches. Steady Compressible Flow If the flow is steady, / t
tion (4.6) reduces to
Cartesian:
Cylindrical: 0 and all properties are functions of position only. Equa x ( u) 1
(r
rr y () 1
r r) ( z
) ( w) z 0 ( z) 0 (4.10) Since density and velocity are both variables, these are still nonlinear and rather formidable, but a number of specialcase solutions have been found. Incompressible Flow A special case which affords great simplification is incompressible flow, where the
density changes are negligible. Then / t 0 regardless of whether the flow is steady
or unsteady, and the density can be slipped out of the divergence in Eq. (4.6) and divided out. The result is
V 0 (4.11) valid for steady or unsteady incompressible flow. The two coordinate forms are
u
x Cartesian: 1
(r r)
rr Cylindrical:  v v 2  y
1
r w
z
() 0 z (4.12a)
( z) 0 See, e.g., Ref. 3, p. 783. eText Main Menu  Textbook Table of Contents  Study Guide (4.12b) 4.2 The Differential Equation of Mass Conservation 221 These are linear differential equations, and a wide variety of solutions are known, as
discussed in Chaps. 6 to 8. Since no author or instructor can resist a wide variety of
solutions, it follows that a great deal of time is spent studying incompressible flows.
Fortunately, this is exactly what should be done, because most practical engineering
flows are approximately incompressible, the chief exception being the highspeed gas
flows treated in Chap. 9.
When is a given flow approximately incompressible? We can derive a nice criterion
by playing a little fast and loose with density approximations. In essence, we wish to
slip the density out of the divergence in Eq. (4.6) and approximate a typical term as,
e.g.,
u
x ( u) x (4.13) This is equivalent to the strong inequality u x x V V u or (4.14) As we shall see in Chap. 9, the pressure change is approximately proportional to the
density change and the square of the speed of sound a of the fluid
p a2 (4.15) Meanwhile, if elevation changes are negligible, the pressure is related to the velocity
change by Bernoulli’s equation (3.75)
p VV (4.16) Combining Eqs. (4.14) to (4.16), we obtain an explicit criterion for incompressible
flow:
V2
a2 Ma2 1 (4.17) where Ma V/a is the dimensionless Mach number of the flow. How small is small?
The commonly accepted limit is
Ma 0.3 (4.18) For air at standard conditions, a flow can thus be considered incompressible if the velocity is less than about 100 m/s (330 ft/s). This encompasses a wide variety of airflows: automobile and train motions, light aircraft, landing and takeoff of highspeed
aircraft, most pipe flows, and turbomachinery at moderate rotational speeds. Further,
it is clear that almost all liquid flows are incompressible, since flow velocities are small
and the speed of sound is very large.3  v v 3
An exception occurs in geophysical flows, where a density change is imposed thermally or mechanically rather than by the flow conditions themselves. An example is fresh water layered upon saltwater or
warm air layered upon cold air in the atmosphere. We say that the fluid is stratified, and we must account
for vertical density changes in Eq. (4.6) even if the velocities are small.  eText Main Menu  Textbook Table of Contents  Study Guide Chapter 4 Differential Relations for a Fluid Particle Before attempting to analyze the continuity equation, we shall proceed with the derivation of the momentum and energy equations, so that we can analyze them as a group.
A very clever device called the stream function can often make short work of the continuity equation, but we shall save it until Sec. 4.7.
One further remark is appropriate: The continuity equation is always important and
must always be satisfied for a rational analysis of a flow pattern. Any newly discovered momentum or energy “solution” will ultimately crash in flames when subjected
to critical analysis if it does not also satisfy the continuity equation.
EXAMPLE 4.2
Under what conditions does the velocity field
V
where a1, b1, etc. (a1x b 1y c1z)i (a2x b2y c2z)j (a3x b3y c3z)k const, represent an incompressible flow which conserves mass? Solution
Recalling that V ui
j wk, we see that u
(4.12a) for incompressible continuity, we obtain
x (a1x b 1y c1z) y or (a2x (a1x b2y
b2 a1 b1y c2z)
c3 z c1z), etc. Substituting into Eq. (a3x b3y c3z) 0 0
Ans. At least two of constants a1, b2, and c3 must have opposite signs. Continuity imposes no restrictions whatever on constants b1, c1, a2, c2, a3, and b3, which do not contribute to a mass increase or decrease of a differential element. EXAMPLE 4.3
An incompressible velocity field is given by
u a(x2 y2) unknown w b where a and b are constants. What must the form of the velocity component be? Solution
Again Eq. (4.12a) applies
x (ax2 ay2) or y y b
z 0 2ax (1) This is easily integrated partially with respect to y
(x, y, z, t)  v v 222  eText Main Menu  2axy f(x, z, t) Textbook Table of Contents Ans.  Study Guide 4.3 The Differential Equation of Linear Momentum 223
This is the only possible form for which satisfies the incompressible continuity equation. The
function of integration f is entirely arbitrary since it vanishes when is differentiated with respect to y.† EXAMPLE 4.4
A centrifugal impeller of 40cm diameter is used to pump hydrogen at 15°C and 1atm pressure.
What is the maximum allowable impeller rotational speed to avoid compressibility effects at the
blade tips? Solution
The speed of sound of hydrogen for these conditions is a 1300 m/s. Assume that the gas velocity leaving the impeller is approximately equal to the impellertip speed
V r 1
2 D Our rule of thumb, Eq. (4.18), neglects compressibility if
1
2 V
1
2 or (0.4 m) D 0.3a 390 m/s 390 m/s 1950 rad/s Thus we estimate the allowable speed to be quite large
310 r/s (18,600 r/min) Ans. An impeller moving at this speed in air would create shock waves at the tips but not in a light
gas like hydrogen. 4.3 The Differential Equation
of Linear Momentum Having done it once in Sec. 4.2 for mass conservation, we can move along a little faster
this time. We use the same elemental control volume as in Fig. 4.1, for which the appropriate form of the linearmomentum relation is
F Vd t (m iVi)out
˙ CV (m iVi)in
˙ (3.40) Again the element is so small that the volume integral simply reduces to a derivative
term
t Vd
CV t ( V) dx dy dz (4.19) The momentum fluxes occur on all six faces, three inlets and three outlets. Referring again to Fig. 4.1, we can form a table of momentum fluxes by exact analogy with
the discussion which led up to the equation for net mass flux:  v v †
This is a very realistic flow which simulates the turning of an inviscid fluid through a 60° angle; see
Examples 4.7 and 4.9.  eText Main Menu  Textbook Table of Contents  Study Guide Chapter 4 Differential Relations for a Fluid Particle Faces Inlet momentum flux x uV dy dz uV y V dx dz V z Outlet momentum flux wV dx dy wV x
y
z ( uV) dx dy dz
( V) dy dx dz
( wV) dz dx dy Introduce these terms and Eq. (4.19) into Eq. (3.40), and get the intermediate result
F dx dy dz t ( V) x ( uV) y ( V) z ( wV) (4.20) Note that this is a vector relation. A simplification occurs if we split up the term in
brackets as follows:
t ( V) x ( uV)
V y ( V) z ( wV) ( V) t V
t V
x u V
y w V
z (4.21) The term in brackets on the righthand side is seen to be the equation of continuity,
Eq. (4.6), which vanishes identically. The long term in parentheses on the righthand
side is seen from Eq. (4.2) to be the total acceleration of a particle which instantaneously occupies the control volume
V
t u V
x V
y w V
z dV
dt (4.2) Thus we have now reduced Eq. (4.20) to
F dV
dx dy dz
dt (4.22) It might be good for you to stop and rest now and think about what we have just done.
What is the relation between Eqs. (4.22) and (3.40) for an infinitesimal control volume? Could we have begun the analysis at Eq. (4.22)?
Equation (4.22) points out that the net force on the control volume must be of differential size and proportional to the element volume. These forces are of two types,
body forces and surface forces. Body forces are due to external fields (gravity, magnetism, electric potential) which act upon the entire mass within the element. The only
body force we shall consider in this book is gravity. The gravity force on the differential mass dx dy dz within the control volume is
dFgrav g dx dy dz (4.23) where g may in general have an arbitrary orientation with respect to the coordinate
system. In many applications, such as Bernoulli’s equation, we take z “up,” and
g
gk.  v v 224  eText Main Menu  Textbook Table of Contents  Study Guide 4.3 The Differential Equation of Linear Momentum 225 The surface forces are due to the stresses on the sides of the control surface. These
stresses, as discussed in Chap. 2, are the sum of hydrostatic pressure plus viscous
stresses ij which arise from motion with velocity gradients ij p xx yx p xy
xz zx
yy zy p yz zz (4.24) The subscript notation for stresses is given in Fig. 4.3.
It is not these stresses but their gradients, or differences, which cause a net force
on the differential control surface. This is seen by referring to Fig. 4.4, which shows
y σy y σy x σy z σx y σz y
σx x
σx z σz x x σz z σi j = Stress in j
direction on a face
normal to i axis z Fig. 4.3 Notation for stresses. (σyx + y ∂σyx
∂y dy) dx dz σz x dx dy
σx x d y dz (σx x +
dy ∂σx x
dx) dy dz
∂x x σy x dx dz
dz
dx  v v Fig. 4.4 Elemental cartesian fixed
control volume showing the surface
forces in the x direction only.  z
(σzx + eText Main Menu  ∂σzx
∂z dz) dx dy Textbook Table of Contents  Study Guide Chapter 4 Differential Relations for a Fluid Particle only the xdirected stresses to avoid cluttering up the drawing. For example, the leftward force xx dy dz on the left face is balanced by the rightward force xx dy dz on
the right face, leaving only the net rightward force ( xx/ x) dx dy dz on the right face.
The same thing happens on the other four faces, so that the net surface force in the x
direction is given by
dFx,surf ( x xx) y ( yx) z ( zx) dx dy dz (4.25) We see that this force is proportional to the element volume. Notice that the stress terms
are taken from the top row of the array in Eq. (4.24). Splitting this row into pressure
plus viscous stresses, we can rewrite Eq. (4.25) as
dFx
d p
x ( x xx) y ( yx) z ( zx) (4.26) In exactly similar manner, we can derive the y and z forces per unit volume on the
control surface
dFy
d p
y x dFz
d p
z x ( xy) ( xz) y
y ( yy) ( yz) z
z ( zy) ( zz) (4.27) Now we multiply Eqs. (4.26) and (4.27) by i, j, and k, respectively, and add to obtain
an expression for the net vector surface force
dF
d dF
d p
surf (4.28)
viscous where the viscous force has a total of nine terms:
dF
d i
viscous j
k xx yx zx x y z xy yy zy x y z xz yz zz x y z (4.29) Since each term in parentheses in (4.29) represents the divergence of a stresscomponent vector acting on the x, y, and z faces, respectively, Eq. (4.29) is sometimes expressed in divergence form
dF
d v   eText Main Menu  yx zx xy yy zy xz ij (4.30) ij
viscous
xx where v 226 yz zz Textbook Table of Contents (4.31)  Study Guide 4.3 The Differential Equation of Linear Momentum 227 is the viscousstress tensor acting on the element. The surface force is thus the sum of
the pressuregradient vector and the divergence of the viscousstress tensor. Substituting into Eq. (4.22) and utilizing Eq. (4.23), we have the basic differential momentum
equation for an infinitesimal element
g
dV
dt where V
t p
u dV
dt ij V
x V
y (4.32)
V
z w (4.33) We can also express Eq. (4.32) in words:
Gravity force per unit volume pressure force per unit volume viscous force per unit volume density acceleration (4.34) Equation (4.32) is so brief and compact that its inherent complexity is almost invisible. It is a vector equation, each of whose component equations contains nine terms.
Let us therefore write out the component equations in full to illustrate the mathematical difficulties inherent in the momentum equation:
gx p
x gy p
y gz p
z xz yz zz x y z xx yx y z xy yy zy x y u
t zx x z u
u t
w
t u u
x u
y x y w
x w
y w
w
w u
z
(4.35) z
w
z This is the differential momentum equation in its full glory, and it is valid for any fluid
in any general motion, particular fluids being characterized by particular viscousstress
terms. Note that the last three “convective” terms on the righthand side of each component equation in (4.35) are nonlinear, which complicates the general mathematical
analysis. Inviscid Flow: Euler’s Equation Equation (4.35) is not ready to use until we write the viscous stresses in terms of velocity components. The simplest assumption is frictionless flow ij 0, for which Eq.
(4.35) reduces to
g p dV
dt (4.36) This is Euler’s equation for inviscid flow. We show in Sec. 4.9 that Euler’s equation
can be integrated along a streamline to yield the frictionless Bernoulli equation, (3.75)
or (3.77). The complete analysis of inviscid flow fields, using continuity and the
Bernoulli relation, is given in Chap. 8.  v v Newtonian Fluid:
NavierStokes Equations  For a newtonian fluid, as discussed in Sec. 1.7, the viscous stresses are proportional to
the element strain rates and the coefficient of viscosity. For incompressible flow, the eText Main Menu  Textbook Table of Contents  Study Guide Chapter 4 Differential Relations for a Fluid Particle generalization of Eq. (1.23) to threedimensional viscous flow is4
u
x 2 xx xy yy u
y yx x yz 2 y
xz zy w
z 2 zz u
z w
x zx (4.37) w
y z where is the viscosity coefficient. Substitution into Eq. (4.35) gives the differential
momentum equation for a newtonian fluid with constant density and viscosity
p
x 2 u
x2 2 u
y2 2 gx u
z2 du
dt p
y 2 2 2 gy d
dt p
z 2 gz 2 x w
x2 2 z2 2 2 y w
y2 w
z2 (4.38) dw
dt These are the NavierStokes equations, named after C. L. M. H. Navier (1785 – 1836)
and Sir George G. Stokes (1819 – 1903), who are credited with their derivation. They
are secondorder nonlinear partial differential equations and are quite formidable, but
surprisingly many solutions have been found to a variety of interesting viscousflow
problems, some of which are discussed in Sec. 4.11 and in Chap. 6 (see also Refs. 4
and 5). For compressible flow, see eq. (2.29) of Ref. 5.
Equation (4.38) has four unknowns: p, u, , and w. It should be combined with the
incompressible continuity relation (4.12) to form four equations in these four unknowns.
We shall discuss this again in Sec. 4.6, which presents the appropriate boundary conditions for these equations.
EXAMPLE 4.5
Take the velocity field of Example 4.3, with b
2 u a(x 0 for algebraic convenience 2 y) 2axy w 0 and determine under what conditions it is a solution to the NavierStokes momentum equation
(4.38). Assuming that these conditions are met, determine the resulting pressure distribution when
z is “up” (gx 0, gy 0, gz
g). Solution
Make a direct substitution of u, , w into Eq. (4.38):
(0) p
x (2a 2a) 2a2 (x3 xy2) (1) 4
When compressibility is significant, additional small terms arise containing the element volume expansion rate and a second coefficient of viscosity; see Refs. 4 and 5 for details.  v v 228  eText Main Menu  Textbook Table of Contents  Study Guide 4.3 The Differential Equation of Linear Momentum 229 p
y (0) ( g) (0) p
z 2a2 (x2y
(0) (2) 0 (3) is not zero). Equation (3) can be integrated The viscous terms vanish identically (although
partially to obtain
p y3) gz f1(x, y) (4) i.e., the pressure is hydrostatic in the z direction, which follows anyway from the fact that the
flow is twodimensional (w 0). Now the question is: Do Eqs. (1) and (2) show that the given
velocity field is a solution? One way to find out is to form the mixed derivative 2p/( x y) from
(1) and (2) separately and then compare them.
Differentiate Eq. (1) with respect to y
2 p
xy 4a2 xy (5) Now differentiate Eq. (2) with respect to x
2 p
xy x [2a2 (x2y y3)] 4a2 xy (6) Since these are identical, the given velocity field is an exact solution to the NavierStokes
equation.
Ans.
To find the pressure distribution, substitute Eq. (4) into Eqs. (1) and (2), which will enable
us to find f1(x, y)
f1
x 2a2 (x3 xy2) (7) f1
y 2a2 (x2y y3) (8) Integrate Eq. (7) partially with respect to x
12
2 a (x4 f1 2x2y2) f2(y) (9) Differentiate this with respect to y and compare with Eq. (8)
f1
y 2a2 x2y f 2(y) (10) Comparing (8) and (10), we see they are equivalent if
2a2 y3 f 2(y)
or f2(y) 12
2 a y4 C (11) where C is a constant. Combine Eqs. (4), (9), and (11) to give the complete expression for pressure distribution
p(x, y, z) gz 12
2 a (x4 y4 2x2y2) C Ans. (12)  v v This is the desired solution. Do you recognize it? Not unless you go back to the beginning and
square the velocity components:  eText Main Menu  Textbook Table of Contents  Study Guide 230 Chapter 4 Differential Relations for a Fluid Particle
u2 2 w2 V2 a2(x4 y4 2x2y2) (13) Comparing with Eq. (12), we can rewrite the pressure distribution as
1
2 p V2 gz C (14) This is Bernoulli’s equation (3.77). That is no accident, because the velocity distribution given in
this problem is one of a family of flows which are solutions to the NavierStokes equation and
which satisfy Bernoulli’s incompressible equation everywhere in the flow field. They are called
irrotational flows, for which curl V
V
0. This subject is discussed again in Sec. 4.9. 4.4 The Differential Equation
of Angular Momentum Having now been through the same approach for both mass and linear momentum, we
can go rapidly through a derivation of the differential angularmomentum relation. The
appropriate form of the integral angularmomentum equation for a fixed control volume is
MO (r t V) d (r CV V) (V n) dA (3.55) CS We shall confine ourselves to an axis O which is parallel to the z axis and passes through
the centroid of the elemental control volume. This is shown in Fig. 4.5. Let be the
angle of rotation about O of the fluid within the control volume. The only stresses
which have moments about O are the shear stresses xy and yx. We can evaluate the
moments about O and the angularmomentum terms about O. A lot of algebra is involved, and we give here only the result
xy 1
2 yx x ( xy) dx 1
2 y ( yx) dy dx dy dz 1
(dx dy dz)(dx2
12 dy2) d2
dt2 (4.39) Assuming that the angular acceleration d2 /dt2 is not infinite, we can neglect all higherτ yx + ∂ (τ yx) d y
∂y θ = Rotation
angle
τ xy dy
Axis O  v v Fig. 4.5 Elemental cartesian fixed
control volume showing shear
stresses which may cause a net angular acceleration about axis O. τ x y + ∂ (τ x y) d x
∂x dx τ yx  eText Main Menu  Textbook Table of Contents  Study Guide 4.5 The Differential Equation of Energy 231 order differential terms, which leaves a finite and interesting result
xy (4.40) yx Had we summed moments about axes parallel to y or x, we would have obtained exactly analogous results
xz zx yz (4.41) zy There is no differential angularmomentum equation. Application of the integral theorem to a differential element gives the result, well known to students of stress analysis, that the shear stresses are symmetric: ij
ji. This is the only result of this section.5 There is no differential equation to remember, which leaves room in your brain
for the next topic, the differential energy equation. 4.5 The Differential Equation
of Energy6 We are now so used to this type of derivation that we can race through the energy equation at a bewildering pace. The appropriate integral relation for the fixed control volume of Fig. 4.1 is
˙
Q ˙
Ws ˙
W ed t e CV p (V n) dA (3.63) CS ˙
where Ws 0 because there can be no infinitesimal shaft protruding into the control
volume. By analogy with Eq. (4.20), the righthand side becomes, for this tiny element,
˙
Q ˙
W where
e
this becomes ( e) t x (u) y ( ) z ( w ) dx dy dz (4.42) p/ . When we use the continuity equation by analogy with Eq. (4.21),
˙
Q de
dt ˙
W V p dx dy dz (4.43) ˙
To evaluate Q, we neglect radiation and consider only heat conduction through the sides
of the element. The heat flow by conduction follows Fourier’s law from Chap. 1
q k T (1.29a) where k is the coefficient of thermal conductivity of the fluid. Figure 4.6 shows the
heat flow passing through the x faces, the y and z heat flows being omitted for clarity.
We can list these six heatflux terms:
Faces Inlet heat flux Outlet heat flux x qx dy dz qx y qy dx dz qy z qz dx dy qz x
y
z (qx) dx dy dz
(qy) dy dx dz
(qz) dz dx dy 5  v v We are neglecting the possibility of a finite couple being applied to the element by some powerful external force field. See, e.g., Ref. 6, p. 217.
6
This section may be omitted without loss of continuity.  eText Main Menu  Textbook Table of Contents  Study Guide 232 Chapter 4 Differential Relations for a Fluid Particle dx Heat flow per
unit area:
∂T
qx = – k
∂x qx + ∂ (qx ) d x
∂x dy
wx + ∂ (wx ) d x
∂x wx Fig. 4.6 Elemental cartesian control
volume showing heatflow and
viscousworkrate terms in the x
direction. Viscous
work rate
per unit
wx = –(uτ x x + υτ x y + w τx z)
area: dz By adding the inlet terms and subtracting the outlet terms, we obtain the net heat
added to the element
˙
Q x (qx) (qy) y (qz) dx dy dz z q dx dy dz (4.44) As expected, the heat flux is proportional to the element volume. Introducing Fourier’s
law from Eq. (1.29), we have
˙
Q (k T) dx dy dz (4.45) The rate of work done by viscous stresses equals the product of the stress component, its corresponding velocity component, and the area of the element face. Figure
4.6 shows the work rate on the left x face is
˙
Wυ,LF wx dy dz where wx (u υ xx w xy xz) (4.46) (where the subscript LF stands for left face) and a slightly different work on the right
face due to the gradient in wx. These work fluxes could be tabulated in exactly the same
manner as the heat fluxes in the previous table, with wx replacing qx, etc. After outlet
terms are subtracted from inlet terms, the net viscouswork rate becomes
˙
W x
z (u xx xy w xz) (u zx zy w zz) (V ij) y (u yx yy w yz) dx dy dz dx dy dz (4.47) We now substitute Eqs. (4.45) and (4.47) into Eq. (4.43) to obtain one form of the differential energy equation
de
dt V p (k T) (V ij) where e û 1
2 V2 gz (4.48) A more useful form is obtained if we split up the viscouswork term  v v (V  eText Main Menu  ij) V( ij) Textbook Table of Contents (4.49)  Study Guide 4.5 The Differential Equation of Energy 233 where is short for the viscousdissipation function.7 For a newtonian incompressible
viscous fluid, this function has the form
u
x 2 2 2 2 y 2 w
y 2 x u
y 2 2 w
x u
z z w
z 2 (4.50) Since all terms are quadratic, viscous dissipation is always positive, so that a viscous
flow always tends to lose its available energy due to dissipation, in accordance with
the second law of thermodynamics.
Now substitute Eq. (4.49) into Eq. (4.48), using the linearmomentum equation (4.32)
to eliminate
ij. This will cause the kinetic and potential energies to cancel, leaving a more customary form of the general differential energy equation
dû
dt p( V) (k T) (4.51) This equation is valid for a newtonian fluid under very general conditions of unsteady,
compressible, viscous, heatconducting flow, except that it neglects radiation heat transfer and internal sources of heat that might occur during a chemical or nuclear reaction.
Equation (4.51) is too difficult to analyze except on a digital computer [1]. It is customary to make the following approximations:
dû c dT c , , k, const (4.52) Equation (4.51) then takes the simpler form
c dT
dt k 2 T (4.53) which involves temperature T as the sole primary variable plus velocity as a secondary
variable through the total timederivative operator
dT
dt T
t T
x u T
y w T
z (4.54) A great many interesting solutions to Eq. (4.53) are known for various flow conditions,
and extended treatments are given in advanced books on viscous flow [4, 5] and books
on heat transfer [7, 8].
One wellknown special case of Eq. (4.53) occurs when the fluid is at rest or has
negligible velocity, where the dissipation
and convective terms become negligible
c T
t 2 k T (4.55) This is called the heatconduction equation in applied mathematics and is valid for
solids and fluids at rest. The solution to Eq. (4.55) for various conditions is a large part
of courses and books on heat transfer.
This completes the derivation of the basic differential equations of fluid motion.
7  v v For further details, see, e.g., Ref. 5, p. 72.  eText Main Menu  Textbook Table of Contents  Study Guide 234 Chapter 4 Differential Relations for a Fluid Particle There are three basic differential equations of fluid motion, just derived. Let us summarize them here: 4.6 Boundary Conditions for
the Basic Equations Continuity: ( V) t
dV
dt Momentum:
dû
dt Energy: g p( 0 (4.56) p V) (4.57) ij (k T) (4.58) where is given by Eq. (4.50). In general, the density is variable, so that these three
equations contain five unknowns, , V, p, û, and T. Therefore we need two additional
relations to complete the system of equations. These are provided by data or algebraic
expressions for the state relations of the thermodynamic properties
( p, T) û û( p, T) (4.59) For example, for a perfect gas with constant specific heats, we complete the system with
p
RT û c dT cT const (4.60) It is shown in advanced books [4, 5] that this system of equations (4.56) to (4.59) is
well posed and can be solved analytically or numerically, subject to the proper boundary conditions.
What are the proper boundary conditions? First, if the flow is unsteady, there must
be an initial condition or initial spatial distribution known for each variable:
At t , V, p, û, T 0: known f(x, y, z) (4.61) Thereafter, for all times t to be analyzed, we must know something about the variables
at each boundary enclosing the flow.
Figure 4.7 illustrates the three most common types of boundaries encountered in
fluidflow analysis: a solid wall, an inlet or outlet, a liquidgas interface.
First, for a solid, impermeable wall, there is no slip and no temperature jump in a
viscous heatconducting fluid
Vfluid Vwall Tfluid Twall solid wall (4.62) The only exception to Eq. (4.62) occurs in an extremely rarefied gas flow, where slippage can be present [5].
Second, at any inlet or outlet section of the flow, the complete distribution of velocity, pressure, and temperature must be known for all times:
Inlet or outlet: Known V, p, T (4.63) These inlet and outlet sections can be and often are at
, simulating a body immersed in an infinite expanse of fluid.
Finally, the most complex conditions occur at a liquidgas interface, or free surface,
as sketched in Fig. 4.7. Let us denote the interface by  v v Interface:  eText Main Menu z  (x, y, t) Textbook Table of Contents (4.64)  Study Guide 4.6 Boundary Conditions for the Basic Equations 235 Liquidgas interface z = η(x, y, t):
pliq = pgas – (R–1 + R–1)
x
y
dη
wliq = wgas =
dt
Equality of q and τ across interface Z Gas
Liquid
Inlet:
known V, p, T Outlet:
known V, p, T Solid contact:
( V, T )fluid = ( V, T )wall Fig. 4.7 Typical boundary conditions in a viscous heatconducting
fluidflow analysis. Solid impermeable wall Then there must be equality of vertical velocity across the interface, so that no holes
appear between liquid and gas:
wliq d
dt wgas u t x y (4.65) This is called the kinematic boundary condition.
There must be mechanical equilibrium across the interface. The viscousshear
stresses must balance
( zy)liq ( zy)gas ( zx)liq ( zx)gas (4.66) Neglecting the viscous normal stresses, the pressures must balance at the interface except for surfacetension effects
pliq pgas (R x 1 R y 1) (4.67)  v v which is equivalent to Eq. (1.34). The radii of curvature can be written in terms of the
freesurface position
x
Rx 1 Ry 1
x
1(
x)2 (
y)2
y
(4.68)
y
1(
x)2 (
y)2  eText Main Menu  Textbook Table of Contents  Study Guide 236 Chapter 4 Differential Relations for a Fluid Particle Finally, the heat transfer must be the same on both sides of the interface, since no
heat can be stored in the infinitesimally thin interface
(qz)liq (qz)gas (4.69) Neglecting radiation, this is equivalent to
k T
z T
z k
liq (4.70)
gas This is as much detail as we wish to give at this level of exposition. Further and even
more complicated details on fluidflow boundary conditions are given in Refs. 5 and 9. Simplified FreeSurface
Conditions In the introductory analyses given in this book, such as openchannel flows in Chap.
10, we shall back away from the exact conditions (4.65) to (4.69) and assume that the
upper fluid is an “atmosphere” which merely exerts pressure upon the lower fluid, with
shear and heat conduction negligible. We also neglect nonlinear terms involving the
slopes of the free surface. We then have a much simpler and linear set of conditions at
the surface
2 pliq pgas 2
2 V
z wliq y2 x T
z 0
liq t 0 (4.71) liq In many cases, such as openchannel flow, we can also neglect surface tension, so that
patm pliq (4.72) These are the types of approximations which will be used in Chap. 10. The nondimensional forms of these conditions will also be useful in Chap. 5. Incompressible Flow with
Constant Properties Flow with constant , , and k is a basic simplification which will be used, e.g., throughout Chap. 6. The basic equations of motion (4.56) to (4.58) reduce to:
Continuity: V
dV
dt Momentum:
Energy: c g
dT
dt 0
2 p
k (4.73)
V (4.74) 2 T (4.75) Since is constant, there are only three unknowns: p, V, and T. The system is closed.8
Not only that, the system splits apart: Continuity and momentum are independent of
T. Thus we can solve Eqs. (4.73) and (4.74) entirely separately for the pressure and
velocity, using such boundary conditions as
Solid surface: V Vwall (4.76) 8  v v For this system, what are the thermodynamic equivalents to Eq. (4.59)?  eText Main Menu  Textbook Table of Contents  Study Guide 4.6 Boundary Conditions for the Basic Equations Inlet or outlet: Known V, p Free surface: p pa w 237 (4.77)
(4.78) t Later, entirely at our leisure,9 we can solve for the temperature distribution from Eq.
(4.75), which depends upon velocity V through the dissipation and the total timederivative operator d/dt. InviscidFlow Approximations Chapter 8 assumes inviscid flow throughout, for which the viscosity
mentum equation (4.74) reduces to
dV
dt g 0. The mo p (4.79) This is Euler’s equation; it can be integrated along a streamline to obtain Bernoulli’s
equation (see Sec. 4.9). By neglecting viscosity we have lost the secondorder derivative of V in Eq. (4.74); therefore we must relax one boundary condition on velocity.
The only mathematically sensible condition to drop is the noslip condition at the wall.
We let the flow slip parallel to the wall but do not allow it to flow into the wall. The
proper inviscid condition is that the normal velocities must match at any solid surface:
Inviscid flow: (Vn)fluid (Vn)wall (4.80) In most cases the wall is fixed; therefore the proper inviscidflow condition is
Vn 0 (4.81) There is no condition whatever on the tangential velocity component at the wall in inviscid flow. The tangential velocity will be part of the solution, and the correct value
will appear after the analysis is completed (see Chap. 8).
EXAMPLE 4.6
For steady incompressible laminar flow through a long tube, the velocity distribution is given
by
z U1 r2
R2 r 0 where U is the maximum, or centerline, velocity and R is the tube radius. If the wall temperature is constant at Tw and the temperature T T(r) only, find T(r) for this flow. Solution
With T T(r), Eq. (4.75) reduces for steady flow to
c r dT
dr kd
dT
r
r dr
dr dz
dr 2 (1)  v v 9
Since temperature is entirely uncoupled by this assumption, we may never get around to solving for it
here and may ask you to wait until a course on heat transfer.  eText Main Menu  Textbook Table of Contents  Study Guide 238 Chapter 4 Differential Relations for a Fluid Particle But since
to obtain r 0 for this flow, the convective term on the left vanishes. Introduce
2 kd
dT
r
r dr
dr dz
dr z into Eq. (1) 4U2 r2
R4 2 (2) Multiply through by r/k and integrate once:
r U2r4
kR4 dT
dr (3) C1 Divide through by r and integrate once again:
U2r4
4kR4 T C1 ln r C2 (4) Now we are in position to apply our boundary conditions to evaluate C1 and C2.
First, since the logarithm of zero is
, the temperature at r 0 will be infinite unless
C1 0 (5) Thus we eliminate the possibility of a logarithmic singularity. The same thing will happen if we
apply the symmetry condition dT/dr 0 at r 0 to Eq. (3). The constant C2 is then found by
the walltemperature condition at r R
T
or U2
4k Tw
C2 Tw C2 U2
4k (6) The correct solution is thus
T(r) U2
1
4k Tw r4
R4 Ans. which is a fourthorder parabolic distribution with a maximum value T0
centerline. 4.7 The Stream Function Cylindrical: v v U2/(4k) at the We have seen in Sec. 4.6 that even if the temperature is uncoupled from our system of
equations of motion, we must solve the continuity and momentum equations simultaneously for pressure and velocity. The stream function is a clever device which allows us to wipe out the continuity equation and solve the momentum equation directly
for the single variable .
The streamfunction idea works only if the continuity equation (4.56) can be reduced to two terms. In general, we have four terms:
Cartesian:  Tw (7)  t ( u) 1
(r
rr t eText Main Menu x  r) y
1
r ()
( z ( w) ) Textbook Table of Contents ( z  0
z) (4.82a)
0 Study Guide (4.82b) 4.7 The Stream Function 239 First, let us eliminate unsteady flow, which is a peculiar and unrealistic application of
the streamfunction idea. Reduce either of Eqs. (4.82) to any two terms. The most common application is incompressible flow in the xy plane
u
x 0 y (4.83) This equation is satisfied identically if a function (x, y) is defined such that Eq. (4.83)
becomes
x y y 0 x (4.84) Comparison of (4.83) and (4.84) shows that this new function
that
u y or V must be defined such (4.85) x i j y x Is this legitimate? Yes, it is just a mathematical trick of replacing two variables (u and
) by a single higherorder function . The vorticity, or curl V, is an interesting function
2 curl V 2k k z 2 2 x2 y2 2 where (4.86) Thus, if we take the curl of the momentum equation (4.74) and utilize Eq. (4.86), we
obtain a single equation for
y x ( 2 ) x y ( 2 2 ) ( 2 ) (4.87) where
/ is the kinematic viscosity. This is partly a victory and partly a defeat:
Eq. (4.87) is scalar and has only one variable, , but it now contains fourthorder
derivatives and probably will require computer analysis. There will be four boundary
conditions required on . For example, for the flow of a uniform stream in the x direction past a solid body, the four conditions would be
At infinity: y At the body: U x y 0
(4.88) 0 x Many examples of numerical solution of Eqs. (4.87) and (4.88) are given in Ref. 1.
One important application is inviscid irrotational flow in the xy plane, where z 0.
Equations (4.86) and (4.87) reduce to
2 2 2
2 x y2 0 (4.89)  v v This is the secondorder Laplace equation (Chap. 8), for which many solutions and analytical techniques are known. Also, boundary conditions like Eq. (4.88) reduce to  eText Main Menu  Textbook Table of Contents  Study Guide 240 Chapter 4 Differential Relations for a Fluid Particle At infinity: Uy At the body: const (4.90) const It is well within our capability to find some useful solutions to Eqs. (4.89) and (4.90),
which we shall do in Chap. 8. Geometric Interpretation of The fancy mathematics above would serve by itself to make the stream function immortal and always useful to engineers. Even better, though, has a beautiful geometric interpretation: Lines of constant are streamlines of the flow. This can be shown
as follows. From Eq. (1.41) the definition of a streamline in twodimensional flow is
dy dx
u
or u dy dx 0 streamline (4.91) Introducing the stream function from Eq. (4.85), we have
x
Thus the change in dx y dy 0 (4.92) d is zero along a streamline, or
const along a streamline (4.93) Having found a given solution (x, y), we can plot lines of constant to give the
streamlines of the flow.
There is also a physical interpretation which relates to volume flow. From Fig.
4.8, we can compute the volume flow dQ through an element ds of control surface of
unit depth
dy
dx
dQ (V n) dA
i
j
i
j
ds(1)
y
x
ds
ds
x dx y dy d (4.94) Control surface
(unit depth
into paper) dQ = ( V • n) d A = dψ V = iu + jv dy
ds  v v Fig. 4.8 Geometric interpretation of
stream function: volume flow
through a differential portion of a
control surface. dx
n=  eText Main Menu  dy
dx
i–
j
ds
ds Textbook Table of Contents  Study Guide 4.7 The Stream Function
ψ 2 < ψ1 ψ 2 > ψ1 Flow Fig. 4.9 Sign convention for flow
in terms of change in stream function: (a) flow to the right if U is
greater; (b) flow to the left if L is
greater. 241 Flow ψ1 ψ1
(a) (b) Thus the change in across the element is numerically equal to the volume flow through
the element. The volume flow between any two points in the flow field is equal to the
change in stream function between those points:
2 Q1→2 2 (V n) dA d 1 2 1 1 (4.95) Further, the direction of the flow can be ascertained by noting whether increases or
decreases. As sketched in Fig. 4.9, the flow is to the right if U is greater than L,
where the subscripts stand for upper and lower, as before; otherwise the flow is to the
left.
Both the stream function and the velocity potential were invented by the French
mathematician Joseph Louis Lagrange and published in his treatise on fluid mechanics in 1781.
EXAMPLE 4.7
If a stream function exists for the velocity field of Example 4.5
u a(x2 y2) 2axy w 0 find it, plot it, and interpret it. Solution
Since this flow field was shown expressly in Example 4.3 to satisfy the equation of continuity,
we are pretty sure that a stream function does exist. We can check again to see if
u
x
Substitute: 2ax 0 y ( 2ax) 0 checks Therefore we are certain that a stream function exists. To find , we simply set
u ax2 y  v v x  eText Main Menu  Textbook Table of Contents ay2 (1) 2axy (2)  Study Guide Chapter 4 Differential Relations for a Fluid Particle
and work from either one toward the other. Integrate (1) partially
ax2y ay3
3 f(x) (3) Differentiate (3) with respect to x and compare with (2)
x
Therefore f (x) 0, or f 2axy f (x) To plot this, set C 2axy (4) constant. The complete stream function is thus found
y3
3 a x2y C Ans. (5) 0 for convenience and plot the function
3x2y y3 3
a (6) for constant values of . The result is shown in Fig. E4.7a to be six 60° wedges of circulating
motion, each with identical flow patterns except for the arrows. Once the streamlines are labeled,
the flow directions follow from the sign convention of Fig. 4.9. How can the flow be interpreted?
Since there is slip along all streamlines, no streamline can truly represent a solid surface in a
viscous flow. However, the flow could represent the impingement of three incoming streams at
60, 180, and 300°. This would be a rather unrealistic yet exact solution to the NavierStokes
equation, as we showed in Example 4.5. ψ = 2a a 0 –2a
–a ψ = 2a y
60° 60° 60° 60°
60° ψ = – 2a –a 0 a a Flow around a 60° corner
x –a
–2a The origin is a
stagnation point 2a Incoming stream impinging
against a 120° corner Flow around a
rounded 60° corner E4.7b E4.7a By allowing the flow to slip as a frictionless approximation, we could let any given streamline be a body shape. Some examples are shown in Fig. E4.7b. A stream function also exists in a variety of other physical situations where only
two coordinates are needed to define the flow. Three examples are illustrated here.  v v 242  eText Main Menu  Textbook Table of Contents  Study Guide 4.7 The Stream Function Steady Plane Compressible Flow Suppose now that the density is variable but that w
plane. Then the equation of continuity becomes
x ( u) y 243 0, so that the flow is in the xy () 0 (4.96) We see that this is in exactly the same form as Eq. (4.84). Therefore a compressibleflow stream function can be defined such that
u y (4.97) x Again lines of constant are streamlines of the flow, but the change in
to the mass flow, not the volume flow
dm
˙ (V n) dA is now equal d 2 or m1→2
˙ (V n) dA 2 1 (4.98) 1 The sign convention on flow direction is the same as in Fig. 4.9. This particular stream
function combines density with velocity and must be substituted into not only momentum but also the energy and state relations (4.58) and (4.59) with pressure and temperature as companion variables. Thus the compressible stream function is not a great
victory, and further assumptions must be made to effect an analytical solution to a typical problem (see, e.g., Ref. 5, chap. 7). Incompressible Plane Flow in
Polar Coordinates Suppose that the important coordinates are r and , with
is constant. Then Eq. (4.82b) reduces to
1
r r (r r) 1
r () z 0, and that the density 0 (4.99) After multiplying through by r, we see that this is the same as the analogous form of
Eq. (4.84)
r 0 r (4.100) By comparison of (4.99) and (4.100) we deduce the form of the incompressible polarcoordinate stream function
r 1
r r (4.101) Once again lines of constant are streamlines, and the change in is the volume flow
Q1→2
2
1. The sign convention is the same as in Fig. 4.9. This type of stream
function is very useful in analyzing flows with cylinders, vortices, sources, and sinks
(Chap. 8).  v v Incompressible Axisymmetric Flow  As a final example, suppose that the flow is threedimensional (υr, υz) but with no circumferential variations,
/
0 (see Fig. 4.2 for definition of coordinates). Such eText Main Menu  Textbook Table of Contents  Study Guide Chapter 4 Differential Relations for a Fluid Particle a flow is termed axisymmetric, and the flow pattern is the same when viewed on any
meridional plane through the axis of revolution z. For incompressible flow, Eq. (4.82b)
becomes
1
r r (r r) z ( z) 0 (4.102) This doesn’t seem to work: Can’t we get rid of the one r outside? But when we realize that r and z are independent coordinates, Eq. (4.102) can be rewritten as
r (r r) (r z) z 0 (4.103) By analogy with Eq. (4.84), this has the form
r z z 0 r (4.104) By comparing (4.103) and (4.104), we deduce the form of an incompressible axisymmetric stream function (r, z)
r 1
r 1
r z z (4.105) r Here again lines of constant are streamlines, but there is a factor (2 ) in the volume
flow: Q1→2 2 ( 2
1). The sign convention on flow is the same as in Fig. 4.9.
EXAMPLE 4.8
Investigate the stream function in polar coordinates
U sin r R2
r (1) where U and R are constants, a velocity and a length, respectively. Plot the streamlines. What
does the flow represent? Is it a realistic solution to the basic equations? Solution
The streamlines are lines of constant , which has units of square meters per second. Note that
/(UR) is dimensionless. Rewrite Eq. (1) in dimensionless form
UR sin η 1 η r
R (2) 0. From Eq. (1) or (2) this occurs when (a)
0
Of particular interest is the special line
or 180° and (b) r R. Case (a) is the xaxis, and case (b) is a circle of radius R, both of which
are plotted in Fig. E4.8.
For any other nonzero value of it is easiest to pick a value of r and solve for :
sin
In general, there will be two solutions for
ample take /(UR)
1.0:  v v 244  eText Main Menu  /(UR)
r/R R/r (3) because of the symmetry about the yaxis. For ex Textbook Table of Contents  Study Guide 4.8 Vorticity and Irrotationality 245
Streamlines converge,
highvelocity region ψ
= +1
UR
r=R –1 +1
2 0 0 0 0
0 +1 –1
2
–1 Singularity
at origin E4.8 Guess r/R 3.0 2.0 1.8 1.7 1.618 22°
158° Compute 2.5
28°
152° 42°
138° 54°
156° 64°
116° 90° This line is plotted in Fig. E4.8 and passes over the circle r R. You have to watch it, though,
because there is a second curve for /(UR)
1.0 for small r R below the xaxis:
Guess r/R 0.618 Compute 90° 0.6 0.5 0.4 0.3 0.2 0.1 70°
110° 42°
138° 28°
152° 19°
161° 12°
168° 6°
174° This second curve plots as a closed curve inside the circle r R. There is a singularity of infinite velocity and indeterminate flow direction at the origin. Figure E4.8 shows the full pattern.
The given stream function, Eq. (1), is an exact and classic solution to the momentum equation (4.38) for frictionless flow. Outside the circle r R it represents twodimensional inviscid
flow of a uniform stream past a circular cylinder (Sec. 8.3). Inside the circle it represents a rather
unrealistic trapped circulating motion of what is called a line doublet.  v v 4.8 Vorticity and Irrotationality  The assumption of zero fluid angular velocity, or irrotationality, is a very useful simplification. Here we show that angular velocity is associated with the curl of the localvelocity vector.
The differential relations for deformation of a fluid element can be derived by examining Fig. 4.10. Two fluid lines AB and BC, initially perpendicular at time t, move
and deform so that at t dt they have slightly different lengths A B and B C and are
slightly off the perpendicular by angles d and d . Such deformation occurs kinematically because A, B, and C have slightly different velocities when the velocity field V eText Main Menu  Textbook Table of Contents  Study Guide 246 Chapter 4 Differential Relations for a Fluid Particle
∂u d y d t
∂y
A′ dy + dβ ∂υ d y d t
∂y Time: t + dt
C′ B′
dx +
Time t dy B 0 ∂u dx dt
∂x V y Fig. 4.10 Angular velocity and
strain rate of two fluid lines deforming in the xy plane. ∂υ dx dt
∂x dα Line 2
A Line 1 dx C x has spatial gradients. All these differential changes in the motion of A, B, and C are
noted in Fig. 4.10.
We define the angular velocity z about the z axis as the average rate of counterclockwise turning of the two lines
1
2 ωz dβ
dt d
dt (4.106) But from Fig. 4.10, d and d are each directly related to velocity derivatives in the
limit of small dt
lim tan 1 dx dt→0 d ( 1 lim tan d dt→0 / x) dx dt
( u/ x) dx dt x ( u/ y) dy dt
dy ( / y) dy dt dt u
dt
y (4.107) Combining Eqs. (4.106) and (4.107) gives the desired result:
1
2 z u
y x (4.108) In exactly similar manner we determine the other two rates:  v v x  eText Main Menu  1
2 w
y z y 1
2 Textbook Table of Contents w
x u
z  Study Guide (4.109) 4.9 Frictionless Irrotational Flows 247 The vector i j x y k z is thus onehalf the curl of the velocity vector 1
2 1
(curl V)
2 i j k x y z u w (4.110) 1
2 Since the factor of is annoying, many workers prefer to use a vector twice as large,
called the vorticity:
2 curl V (4.111) Many flows have negligible or zero vorticity and are called irrotational
curl V 0 (4.112) The next section expands on this idea. Such flows can be incompressible or compressible, steady or unsteady.
We may also note that Fig. 4.10 demonstrates the shearstrain rate of the element,
which is defined as the rate of closure of the initially perpendicular lines
d
dt d
dt ˙ xy u
y x (4.113) When multiplied by viscosity , this equals the shear stress xy in a newtonian fluid,
as discussed earlier in Eqs. (4.37). Appendix E lists strainrate and vorticity components in cylindrical coordinates. 4.9 Frictionless Irrotational
Flows When a flow is both frictionless and irrotational, pleasant things happen. First, the momentum equation (4.38) reduces to Euler’s equation
dV
dt g p (4.114) Second, there is a great simplification in the acceleration term. Recall from Sec. 4.1
that acceleration has two terms
V
t dV
dt (V )V (4.2) A beautiful vector identity exists for the second term [11]:
(V )V ( 1 V2)
2 V (4.115) where
curl V from Eq. (4.111) is the fluid vorticity.
Now combine (4.114) and (4.115), divide by , and rearrange on the lefthand side.
Dot the entire equation into an arbitrary vector displacement dr:
V
t 12
V
2 V 1 p g dr 0 (4.116) Nothing works right unless we can get rid of the third term. We want  v v (  eText Main Menu  V) (dr) Textbook Table of Contents 0  Study Guide (4.117) 248 Chapter 4 Differential Relations for a Fluid Particle This will be true under various conditions:
1. V is zero; trivial, no flow (hydrostatics).
2.
is zero; irrotational flow.
3. dr is perpendicular to
V; this is rather specialized and rare.
4. dr is parallel to V; we integrate along a streamline (see Sec. 3.7).
Condition 4 is the common assumption. If we integrate along a streamline in frictionless compressible flow and take, for convenience, g
gk, Eq. (4.116) reduces to
V
dr
t d dp 12
V
2 g dz 0 (4.118) Except for the first term, these are exact differentials. Integrate between any two points
1 and 2 along the streamline:
2
1 2 V
ds
t dp 1 1
(V 2
2
2 2
V 1) g(z2 z1) 0 (4.119) where ds is the arc length along the streamline. Equation (4.119) is Bernoulli’s equation for frictionless unsteady flow along a streamline and is identical to Eq. (3.76). For
incompressible steady flow, it reduces to
12
V
2 p gz constant along streamline (4.120) The constant may vary from streamline to streamline unless the flow is also irrotational
(assumption 2). For irrotational flow
0, the offending term Eq. (4.117) vanishes
regardless of the direction of dr, and Eq. (4.120) then holds all over the flow field with
the same constant. Velocity Potential Irrotationality gives rise to a scalar function similar and complementary to the stream
function . From a theorem in vector analysis [11], a vector with zero curl must be the
gradient of a scalar function
If V 0 then V (4.121) where
(x, y, z, t) is called the velocity potential function. Knowledge of
immediately gives the velocity components
u x y w thus (4.122) z Lines of constant are called the potential lines of the flow.
Note that , unlike the stream function, is fully threedimensional and not limited
to two coordinates. It reduces a velocity problem with three unknowns u, , and w to
a single unknown potential ; many examples are given in Chap. 8 and Sec. 4.10. The
velocity potential also simplifies the unsteady Bernoulli equation (4.118) because if
exists, we obtain  v v V
dr
t  eText Main Menu  t ( ) dr d Textbook Table of Contents (4.123) t  Study Guide 4.9 Frictionless Irrotational Flows 249 Equation (4.118) then becomes a relation between
dp 1 2 t 2 and p
gz const (4.124) This is the unsteady irrotational Bernoulli equation. It is very important in the analysis of accelerating flow fields (see, e.g., Refs. 10 and 15), but the only application in
this text will be in Sec. 9.3 for steady flow. Orthogonality of Streamlines and
Potential Lines If a flow is both irrotational and described by only two coordinates, and both exist and the streamlines and potential lines are everywhere mutually perpendicular except at a stagnation point. For example, for incompressible flow in the xy plane, we
would have
u y (4.125) x (4.126) y x Can you tell by inspection not only that these relations imply orthogonality but also
that and satisfy Laplace’s equation?10 A line of constant would be such that the
change in is zero
d x dx y dy 0 u dx dy (4.127) Solving, we have
dy
dx u
const 1
(dy/dx) (4.128)
const Equation (4.128) is the mathematical condition that lines of constant and be mutually orthogonal. It may not be true at a stagnation point, where both u and are zero,
so that their ratio in Eq. (4.128) is indeterminate. Generation of Rotationality This is the second time we have discussed Bernoulli’s equation under different circumstances (the first was in Sec. 3.7). Such reinforcement is useful, since this is probably
the most widely used equation in fluid mechanics. It requires frictionless flow with no
shaft work or heat transfer between sections 1 and 2. The flow may or may not be irrotational, the latter being an easier condition, allowing a universal Bernoulli constant.
The only remaining question is: When is a flow irrotational? In other words, when
does a flow have negligible angular velocity? The exact analysis of fluid rotationality
under arbitrary conditions is a topic for advanced study, e.g., Ref. 10, sec. 8.5; Ref. 9,
sec. 5.2; and Ref. 5, sec. 2.10. We shall simply state those results here without proof.
A fluid flow which is initially irrotational may become rotational if
1. There are significant viscous forces induced by jets, wakes, or solid boundaries.
In this case Bernoulli’s equation will not be valid in such viscous regions.  v v 10
Equations (4.125) and (4.126) are called the CauchyRiemann equations and are studied in complexvariable theory.  eText Main Menu  Textbook Table of Contents  Study Guide 250 Chapter 4 Differential Relations for a Fluid Particle
Viscous regions where Bernoulli's equation fails:
Laminar
boundary
layer U Turbulent
boundary
layer Separated
flow Wake
flow Uniform
approach
flow
(irrotational)
(a)
Curved shock wave introduces rotationality
Viscous regions where Bernoulli is invalid: Fig. 4.11 Typical flow patterns illustrating viscous regions patched
onto nearly frictionless regions:
(a) low subsonic flow past a body
(U a); frictionless, irrotational
potential flow outside the boundary
layer (Bernoulli and Laplace equations valid); (b) supersonic flow
past a body (U a); frictionless,
rotational flow outside the boundary layer (Bernoulli equation valid,
potential flow invalid). Laminar
boundary
layer U Turbulent
boundary
layer Slight
separated
flow Wake
flow Uniform
supersonic
approach
(irrotational)
(b) 2. There are entropy gradients caused by curved shock waves (see Fig. 4.11b).
3. There are density gradients caused by stratification (uneven heating) rather than
by pressure gradients.
4. There are significant noninertial effects such as the earth’s rotation (the Coriolis
acceleration).  v v In cases 2 to 4, Bernoulli’s equation still holds along a streamline if friction is negligible. We shall not study cases 3 and 4 in this book. Case 2 will be treated briefly in
Chap. 9 on gas dynamics. Primarily we are concerned with case 1, where rotation is
induced by viscous stresses. This occurs near solid surfaces, where the noslip condition creates a boundary layer through which the stream velocity drops to zero, and in
jets and wakes, where streams of different velocities meet in a region of high shear.
Internal flows, such as pipes and ducts, are mostly viscous, and the wall layers grow
to meet in the core of the duct. Bernoulli’s equation does not hold in such flows unless it is modified for viscous losses.
External flows, such as a body immersed in a stream, are partly viscous and partly
inviscid, the two regions being patched together at the edge of the shear layer or boundary layer. Two examples are shown in Fig. 4.11. Figure 4.11a shows a lowspeed  eText Main Menu  Textbook Table of Contents  Study Guide 4.9 Frictionless Irrotational Flows 251 subsonic flow past a body. The approach stream is irrotational; i.e., the curl of a constant is zero, but viscous stresses create a rotational shear layer beside and downstream
of the body. Generally speaking (see Chap. 6), the shear layer is laminar, or smooth,
near the front of the body and turbulent, or disorderly, toward the rear. A separated, or
deadwater, region usually occurs near the trailing edge, followed by an unsteady turbulent wake extending far downstream. Some sort of laminar or turbulent viscous theory must be applied to these viscous regions; they are then patched onto the outer flow,
which is frictionless and irrotational. If the stream Mach number is less than about 0.3,
we can combine Eq. (4.122) with the incompressible continuity equation (4.73).
V ( ) 2
2 or 0 0 2
2 2
2 x y (4.129) z2 This is Laplace’s equation in three dimensions, there being no restraint on the number
of coordinates in potential flow. A great deal of Chap. 8 will be concerned with solving Eq. (4.129) for practical engineering problems; it holds in the entire region of Fig.
4.11a outside the shear layer.
Figure 4.11b shows a supersonic flow past a body. A curved shock wave generally
forms in front, and the flow downstream is rotational due to entropy gradients (case
2). We can use Euler’s equation (4.114) in this frictionless region but not potential theory. The shear layers have the same general character as in Fig. 4.11a except that the
separation zone is slight or often absent and the wake is usually thinner. Theory of separated flow is presently qualitative, but we can make quantitative estimates of laminar
and turbulent boundary layers and wakes. EXAMPLE 4.9
If a velocity potential exists for the velocity field of Example 4.5
a(x2 u y2) 2axy w 0 find it, plot it, and compare with Example 4.7. Solution
Since w 0, the curl of V has only one z component, and we must show that it is zero:
( V)z 2 z 2ay x
2ay u
y x 0 ( 2axy) y (ax2 checks ay2)
Ans. The flow is indeed irrotational. A potential exists.
To find (x, y), set
u ax2 x  v v y  eText Main Menu  Textbook Table of Contents ay2
2axy  Study Guide (1)
(2) 252 Chapter 4 Differential Relations for a Fluid Particle φ = –2 a
–a 2a
a
0 0 y –a a φ = –2 a 2a
x φ = 2a a E4.9 –2 a
0 –a Integrate (1)
ax3
3 axy2 2axy f (y) f(y) (3) Differentiate (3) and compare with (2)
y
Therefore f 0, or f 2axy (4) constant. The velocity potential is
ax3
3 axy2 C Ans. Letting C 0, we can plot the lines in the same fashion as in Example 4.7. The result is shown
in Fig. E4.9 (no arrows on ). For this particular problem, the lines form the same pattern as
the lines of Example 4.7 (which are shown here as dashed lines) but are displaced 30°. The
and lines are everywhere perpendicular except at the origin, a stagnation point, where they
are 30° apart. We expected trouble at the stagnation point, and there is no general rule for determining the behavior of the lines at that point. 4.10 Some Illustrative Plane
Potential Flows Chapter 8 is devoted entirely to a detailed study of inviscid incompressible flows, especially those which possess both a stream function and a velocity potential. As sketched
in Fig. 4.11a, inviscid flow is valid away from solid surfaces, and this inviscid pattern
is “patched” onto the nearwall viscous layers — an idea developed in Chap. 7. Various
body shapes can be simulated by the inviscidflow pattern. Here we discuss plane flows,
three of which are illustrated in Fig. 4.12. Uniform Stream in the x Direction A uniform stream V iU, as in Fig. 4.12a, possesses both a stream function and a velocity potential, which may be found as follows:  v v u  eText Main Menu U  x y 0 y Textbook Table of Contents x  Study Guide 4.10 Some Illustrative Plane Potential Flows 253 U m/r Fig. 4.12 Three elementary plane
potential flows. Solid lines are
streamlines; dashed lines are potential lines. K/r
(a) (b) (c) We may integrate each expression and discard the constants of integration, which do
not affect the velocities in the flow. The results are
Uniform stream iU: Uy Ux (4.130) The streamlines are horizontal straight lines (y const), and the potential lines are vertical (x const), i.e., orthogonal to the streamlines, as expected. Line Source or Sink at the Origin Suppose that the zaxis were a sort of thinpipe manifold through which fluid issued
at total rate Q uniformly along its length b. Looking at the xy plane, we would see a
cylindrical radial outflow or line source, as sketched in Fig. 4.12b. Plane polar coordinates are appropriate (see Fig. 4.2), and there is no circumferential velocity. At any
radius r, the velocity is
r Q
2 rb m
r 1
r 0 r r 1
r where we have used the polarcoordinate forms of the stream function and the velocity potential. Integrating and again discarding the constants of integration, we obtain
the proper functions for this simple radial flow:
Line source or sink: m m ln r (4.131) where m Q/(2 b) is a constant, positive for a source, negative for a sink. As shown
in Fig. 4.12b, the streamlines are radial spokes (constant ), and the potential lines are
circles (constant r). Line Irrotational Vortex A (twodimensional) line vortex is a purely circulating steady motion,
f(r) only,
0. This satisfies the continuity equation identically, as may be checked from Eq.
r
(4.12b). We may also note that a variety of velocity distributions (r) satisfy the
momentum equation of a viscous fluid, Eq. (E.6). We may show, as a problem exercise, that only one function (r) is irrotational, i.e., curl V 0, and that is
K/r,
where K is a constant. This is sometimes called a free vortex, for which the stream
function and velocity may be found:  v v r  eText Main Menu  0 1
r K
r r Textbook Table of Contents  r 1
r Study Guide 254 Chapter 4 Differential Relations for a Fluid Particle We may again integrate to determine the appropriate functions:
K ln r K (4.132) where K is a constant called the strength of the vortex. As shown in Fig. 4.12c, the
streamlines are circles (constant r), and the potential lines are radial spokes (constant
). Note the similarity between Eqs. (4.131) and (4.132). A free vortex is a sort of reversed image of a source. The “bathtub vortex,” formed when water drains through a
bottom hole in a tank, is a good approximation to the freevortex pattern. Superposition: Source Plus an
Equal Sink Each of the three elementary flow patterns in Fig. 4.12 is an incompressible irrotational
flow and therefore satisfies both plane “potential flow” equations 2
0 and
2
0. Since these are linear partial differential equations, any sum of such basic
solutions is also a solution. Some of these composite solutions are quite interesting and
useful.
For example, consider a source m at (x, y) ( a, 0), combined with a sink of
equal strength m, placed at ( a, 0), as in Fig. 4.13. The resulting stream function is
simply the sum of the two. In cartesian coordinates,
source sink m tan y 1 x a m tan y 1 x a Similarly, the composite velocity potential is
source sink 1
m ln [(x
2 a)2 y2] 1
m ln [(x
2 a)2  v v Fig. 4.13 Potential flow due to a
line source plus an equal line sink,
from Eq. (4.133). Solid lines are
streamlines; dashed lines are potential lines.  eText Main Menu  Textbook Table of Contents  Study Guide y2] 4.10 Some Illustrative Plane Potential Flows 255 By using trigonometric and logarithmic identities, these may be simplified to
Source plus sink: m tan 2ay
y2 a2 1
2 x (x
1
m ln
2
(x a)2
2 a) (4.133) y2
y2 These lines are plotted in Fig. 4.13 and are seen to be two families of orthogonal
circles, with the streamlines passing through the source and sink and the potential
lines encircling them. They are harmonic (laplacian) functions which are exactly
analogous in electromagnetic theory to the electriccurrent and electricpotential patterns of a magnet with poles at ( a, 0). Sink Plus a Vortex at the Origin An interesting flow pattern, approximated in nature, occurs by superposition of a sink
and a vortex, both centered at the origin. The composite stream function and velocity
potential are
Sink plus vortex: m K ln r m ln r K (4.134) When plotted, these form two orthogonal families of logarithmic spirals, as shown in
Fig. 4.14. This is a fairly realistic simulation of a tornado (where the sink flow moves
up the zaxis into the atmosphere) or a rapidly draining bathtub vortex. At the center
of a real (viscous) vortex, where Eq. (4.134) predicts infinite velocity, the actual circulating flow is highly rotational and approximates solidbody rotation ≈ Cr. y x  v v Fig. 4.14 Superposition of a sink
plus a vortex, Eq. (4.134), simulates a tornado.  eText Main Menu  Textbook Table of Contents  Study Guide 256 Chapter 4 Differential Relations for a Fluid Particle If we superimpose a uniform xdirected stream against an isolated source, a halfbody
shape appears. If the source is at the origin, the combined stream function is, in polar
coordinates, Uniform Stream Plus a Sink at
the Origin: The Rankine HalfBody Uniform stream plus source: Ur sin m (4.135) We can set this equal to various constants and plot the streamlines, as shown in Fig.
4.15. A curved, roughly elliptical, halfbody shape appears, which separates the source
flow from the stream flow. The body shape, which is named after the Scottish engineer W. J. M. Rankine (1820 – 1872), is formed by the particular streamlines
m. The halfwidth of the body far downstream is m/U. The upper surface may be
plotted from the relation
)
m(
U sin r (4.136) Us(max) = 1.26U∞ ψ = +π m
πa y
U∞ ψ= 0 x
a Fig. 4.15 Superposition of a source
plus a uniform stream forms a
Rankine halfbody. ψ = – πm It is not a true ellipse. The nose of the body, which is a “stagnation” point where V
0, stands at (x, y) ( a, 0), where a m/U. The streamline
0 also crosses this
point — recall that streamlines can cross only at a stagnation point.
The cartesian velocity components are found by differentiation:
u y U m
cos
r m
sin
r x Setting u
0, we find a single stagnation point at
180° and r
(x, y) ( m/U, 0), as stated. The resultant velocity at any point is
V2 u2 2 U2 1 a2
r2 2a
cos
r (4.137)
m/U, or (4.138)  v v where we have substituted m Ua. If we evaluate the velocities along the upper surface
m, we find a maximum value Us,max 1.26U at
63°. This point is labeled in Fig. 4.15 and, by Bernoulli’s equation, is the point of minimum pressure on  eText Main Menu  Textbook Table of Contents  Study Guide 4.10 Some Illustrative Plane Potential Flows 257 the body surface. After this point, the surface flow decelerates, the pressure rises, and
the viscous layer grows thicker and more susceptible to “flow separation,” as we shall
see in Chap. 7.
EXAMPLE 4.10
The bottom of a river has a 4mhigh bump which approximates a Rankine halfbody, as in
Fig. E4.10. The pressure at point B on the bottom is 130 kPa, and the river velocity is
2.5 m/s. Use inviscid theory to estimate the water pressure at point A on the bump, which is
2 m above point B. Water at 20°C
A
2.5 m/s 4m
2m E4.10 B 0 Solution
As in all inviscid theories, we ignore the lowvelocity boundary layers which form on solid surfaces due to the noslip condition. From Eq. (4.136) and Fig. 4.15, the downstream bump halfheight equals a. Therefore, for our case, a (4 m)/
1.27 m. We have to find the spot where
the bump height is half that much, h 2 m
a/2. From Eq. (4.136) we may compute
r a( hA )
2 sin a or 90° 2 Thus point A in Fig. E4.10 is directly above the (initially unknown) origin of coordinates (labeled O in Fig. E4.10) and is 1.27 m to the right of the nose of the bump. With r
a/2 and
/2 known, we compute the velocity at point A from Eq. (4.138):
2
VA or U2 1
VA a2
( a/2)2 1.185U 2a
cos
2
a/2 1.185(2.5 m/s) 2 1.405U2
2.96 m/s 3 For water at 20°C, take
998 kg/m and
9790 N/m . Now, since the velocity and elevation are known at point A, we are in a position to use Bernoulli’s inviscid, incompressibleflow
equation (4.120) to estimate pA from the known properties at point B (on the same streamline):
pA
pA
9790 N/m3 or V2
A
2g zA (2.96 m/s)2
2(9.81 m/s2) 2
VB
2g pB zB 130,000
9790 2m (2.5)2
2(9.81) 0 Solving, we find  v v pA  eText Main Menu  (13.60 2.45)(9790) Textbook Table of Contents  109,200 Pa Study Guide Ans. 258 Chapter 4 Differential Relations for a Fluid Particle
If the approach velocity is uniform, this should be a pretty good approximation, since water is
relatively inviscid and its boundary layers are thin. 4.11 Some Illustrative
Incompressible Viscous Flows The inviscid flows of Sec. 4.10 do not satisfy the noslip condition. They “slip” at the
wall but do not flow through the wall. To look at fully viscous noslip conditions, we
must attack the complete NavierStokes equation (4.74), and the result is usually not
at all irrotational, nor does a velocity potential exist. We look here at three cases: (1)
flow between parallel plates due to a moving upper wall, (2) flow between parallel
plates due to pressure gradient, and (3) flow between concentric cylinders when the inner one rotates. Other cases will be given as problem assignments or considered in
Chap. 6. Extensive solutions for viscous flows are discussed in Refs. 4 and 5. Couette Flow between a Fixed and
a Moving Plate Consider twodimensional incompressible plane ( / z 0) viscous flow between parallel plates a distance 2h apart, as shown in Fig. 4.16. We assume that the plates are
w 0.
very wide and very long, so that the flow is essentially axial, u 0 but
The present case is Fig. 4.16a, where the upper plate moves at velocity V but there is
no pressure gradient. Neglect gravity effects. We learn from the continuity equation
(4.73) that
u
x w
z y 0 u
x 0 0 or u u(y) only Thus there is a single nonzero axialvelocity component which varies only across the
channel. The flow is said to be fully developed (far downstream of the entrance). Substitute u u(y) into the xcomponent of the NavierStokes momentum equation (4.74)
for twodimensional (x, y) flow:
u
x u or p
x u
y (0 0) 0 0 2 2 u
x2 gx u
y2 d2u
dy2 0 (4.139) Fixed V
y = +h
y
u( y)
x umax  v v Fig. 4.16 Incompressible viscous
flow between parallel plates: (a) no
pressure gradient, upper plate moving; (b) pressure gradient p/ x
with both plates fixed. u( y)
y = –h
Fixed
(a)  Fixed
(b) eText Main Menu  Textbook Table of Contents  Study Guide 4.11 Some Illustrative Incompressible Viscous Flows 259 Most of the terms drop out, and the momentum equation simply reduces to
d2u
dy2 0 or u C1y C2 The two constants are found by applying the noslip condition at the upper and lower
plates:
At y h: u V C1h At y h: u 0 C1( h) or V
2h C1 C2 and C2
V
2 C2 Therefore the solution for this case (a), flow between plates with a moving upper wall, is
u V
y
2h V
2 h y h (4.140) This is Couette flow due to a moving wall: a linear velocity profile with noslip at each
wall, as anticipated and sketched in Fig. 4.16a. Note that the origin has been placed in
the center of the channel, for convenience in case (b) below.
What we have just presented is a rigorous derivation of the more informally discussed flow of Fig. 1.6 (where y and h were defined differently). Flow due to Pressure Gradient
between Two Fixed Plates Case (b) is sketched in Fig. 4.16b. Both plates are fixed (V 0), but the pressure varies
in the x direction. If
w 0, the continuity equation leads to the same conclusion
as case (a), namely, that u u(y) only. The xmomentum equation (4.138) changes
only because the pressure is variable:
d2u
dy2 p
x (4.141) w 0 and gravity is neglected, the y and zmomentum equations p
y Also, since
lead to 0 p
z and 0 or p p(x) only Thus the pressure gradient in Eq. (4.141) is the total and only gradient:
d2u
dy2 dp
dx const 0 (4.142)  v v Why did we add the fact that dp/dx is constant? Recall a useful conclusion from the
theory of separation of variables: If two quantities are equal and one varies only with
y and the other varies only with x, then they must both equal the same constant. Otherwise they would not be independent of each other.
Why did we state that the constant is negative? Physically, the pressure must decrease in the flow direction in order to drive the flow against resisting wall shear stress.
Thus the velocity profile u(y) must have negative curvature everywhere, as anticipated
and sketched in Fig. 4.16b.  eText Main Menu  Textbook Table of Contents  Study Guide Chapter 4 Differential Relations for a Fluid Particle The solution to Eq. (4.142) is accomplished by double integration:
1 dp y2
dx 2 u C1y C2 The constants are found from the noslip condition at each wall:
At y h: u 0 or C1 0 and dp h2
dx 2 C2 Thus the solution to case (b), flow in a channel due to pressure gradient, is
dp h2
dx 2 u 1 y2
h2 (4.143) The flow forms a Poiseuille parabola of constant negative curvature. The maximum velocity occurs at the centerline y 0:
dp h2
dx 2 umax (4.144) Other (laminar) flow parameters are computed in the following example. EXAMPLE 4.11
For case (b) above, flow between parallel plates due to the pressure gradient, compute (a) the
wall shear stress, (b) the stream function, (c) the vorticity, (d) the velocity potential, and (e) the
average velocity. Solution
All parameters can be computed from the basic solution, Eq. (4.143), by mathematical
manipulation.
(a) The wall shear follows from the definition of a newtonian fluid, Eq. (4.37):
w u
y xy wall x y y h dp
dx h2
2 1 y2
h2 y 2 umax
h dp
h
dx h Ans. (a) The wall shear has the same magnitude at each wall, but by our sign convention of Fig. 4.3,
the upper wall has negative shear stress.
(b) Since the flow is plane, steady, and incompressible, a stream function exists:
u y umax 1 y2
h2 umax y
At the walls, y v   eText Main Menu h and  0 x 0 at the centerline for convenience, we obtain Integrating and setting v 260 y3
3h2 Ans. (b) 2umaxh/3, respectively. Textbook Table of Contents  Study Guide 4.11 Some Illustrative Incompressible Viscous Flows 261
(c) In plane flow, there is only a single nonzero vorticity component:
(curl V)z z u
y x 2umax
y
h2 Ans. (c) The vorticity is highest at the wall and is positive (counterclockwise) in the upper half and
negative (clockwise) in the lower half of the fluid. Viscous flows are typically full of vorticity and are not at all irrotational.
(d) From part (c), the vorticity is finite. Therefore the flow is not irrotational, and the velocity
potential does not exist.
Ans. (d)
(e) The average velocity is defined as Vav Q/A, where Q
u dA over the cross section. For
our particular distribution u(y) from Eq. (4.143), we obtain
1
A Vav u dA 1
b(2h) y2
b dy
h2 h
h umax 1 2
umax
3 Ans. (e) In plane Poiseuille flow between parallel plates, the average velocity is twothirds of the
maximum (or centerline) value. This result could also have been obtained from the stream
function derived in part (b). From Eq. (4.95),
Qchannel upper lower 2umaxh
3 2umaxh
3 4
umaxh per unit width
3 whence Vav Q/Ab 1 (4umaxh/3)/(2h) 2umax/3, the same result.
This example illustrates a statement made earlier: Knowledge of the velocity vector V
[as in Eq. (4.143)] is essentially the solution to a fluidmechanics problem, since all other
flow properties can then be calculated. Flow between Long Concentric
Cylinders Consider a fluid of constant ( , ) between two concentric cylinders, as in Fig. 4.17.
/ z 0. Let the inner cylinder rotate at
There is no axial motion or end effect z
angular velocity i. Let the outer cylinder be fixed. There is circular symmetry, so the
velocity does not vary with and varies only with r.
Fixed ro
Ωi vθ
r
ri Fluid: ρ, µ  v v Fig. 4.17 Coordinate system for
incompressible viscous flow between a fixed outer cylinder and a
steadily rotating inner cylinder.  eText Main Menu  Textbook Table of Contents  Study Guide 262 Chapter 4 Differential Relations for a Fluid Particle The continuity equation for this problem is Eq. (D.2):
1
r r 1
r (r r) 1d
(r r)
r dr 0 or r r const Note that does not vary with . Since r 0 at both the inner and outer cylinders,
it follows that r 0 everywhere and the motion can only be purely circumferential,
(r). The momentum equation (D.6) becomes
(V p 1
r r ) r 2 g r2 For the conditions of the present problem, all terms are zero except the last. Therefore
the basic differential equation for flow between rotating cylinders is
2 d
1d
r
dr
r dr (4.145) r2 This is a linear secondorder ordinary differential equation with the solution
C2
C1r
r
The constants are found by the noslip condition at the inner and outer cylinders:
Outer, at r ro: Inner, at r 0 ri: C2
ro C1ro iri C2
ri C1ri The final solution for the velocity distribution is
Rotating inner cylinder: iri ro/r
ro/ri r/ro
ri/ro (4.146) The velocity profile closely resembles the sketch in Fig. 4.17. Variations of this case,
such as a rotating outer cylinder, are given in the problem assignments.
The classic Couetteflow solution11 of Eq. (4.146) describes a physically satisfying concave, twodimensional, laminarflow velocity profile as in Fig. 4.17. The solution is
mathematically exact for an incompressible fluid. However, it becomes unstable at a
relatively low rate of rotation of the inner cylinder, as shown in 1923 in a classic paper by G. I. Taylor [17]. At a critical value of what is now called the Taylor number,
denoted Ta, Instability of Rotating Inner
Cylinder Flow Tacrit ri(ro ri)3
2 2
i 1700 (4.147) the plane flow of Fig. 4.17 vanishes and is replaced by a laminar threedimensional
flow pattern consisting of rows of nearly square alternating toroidal vortices. An ex  v v 11
Named after M. Couette, whose pioneering paper in 1890 established rotating cylinders as a method,
still used today, for measuring the viscosity of fluids.  eText Main Menu  Textbook Table of Contents  Study Guide 4.11 Some Illustrative Incompressible Viscous Flows 263 (a) Fig. 4.18 Experimental verification
of the instability of flow between a
fixed outer and a rotating inner
cylinder. (a) Toroidal Taylor vortices exist at 1.16 times the critical
speed; (b) at 8.5 times the critical
speed, the vortices are doubly periodic. (After Koschmieder, Ref. 18.)
This instability does not occur if
only the outer cylinder rotates. (b) perimental demonstration of toroidal “Taylor vortices” is shown in Fig. 4.18a, measured at Ta 1.16 Tacrit by Koschmieder [18]. At higher Taylor numbers, the vortices
also develop a circumferential periodicity but are still laminar, as illustrated in Fig.
4.18b. At still higher Ta, turbulence ensues. This interesting instability reminds us that
the NavierStokes equations, being nonlinear, do admit to multiple (nonunique) laminar solutions in addition to the usual instabilities associated with turbulence and chaotic
dynamic systems. Summary  v v This chapter complements Chap. 3 by using an infinitesimal control volume to derive
the basic partial differential equations of mass, momentum, and energy for a fluid.
These equations, together with thermodynamic state relations for the fluid and appro  eText Main Menu  Textbook Table of Contents  Study Guide 264 Chapter 4 Differential Relations for a Fluid Particle priate boundary conditions, in principle can be solved for the complete flow field in
any given fluidmechanics problem. Except for Chap. 9, in most of the problems to be
studied here an incompressible fluid with constant viscosity is assumed.
In addition to deriving the basic equations of mass, momentum, and energy, this
chapter introduced some supplementary ideas — the stream function, vorticity, irrotationality, and the velocity potential — which will be useful in coming chapters, especially Chap. 8. Temperature and density variations will be neglected except in Chap.
9, where compressibility is studied.
This chapter ended by discussing a few classical solutions for inviscid flows (uniform stream, source, sink, vortex, halfbody) and for viscous flows (Couette flow due
to moving walls and Poiseuille flow due to pressure gradient). Whole books [11, 13]
are written on the basic equations of fluid mechanics. Whole books [4, 5, 15] are written on classical solutions to fluidflow problems. Reference 12 contains 360 solved
problems which relate fluid mechanics to the whole of continuum mechanics. This does
not mean that all problems can be readily solved mathematically, even with the modern digitalcomputer codes now available. Often the geometry and boundary conditions
are so complex that experimentation (Chap. 5) is a necessity. Problems
Most of the problems herein are fairly straightforward. More difficult or openended assignments are labeled with an asterisk. Problems labeled with an EES icon will benefit from the use of the Engineering Equation Solver (EES), while problems labeled with a
computer disk may require the use of a computer. The standard
endofchapter problems 4.1 to 4.91 (categorized in the problem
list below) are followed by word problems W4.1 to W4.10, fundamentals of engineering exam problems FE4.1 to FE4.3, and comprehensive problem C4.1.
Problem distribution
Section
4.1
4.2
4.3
4.4
4.5
4.6
4.7
4.8
4.9
4.10
4.11 Topic Problems The acceleration of a fluid
The continuity equation
Linear momentum: NavierStokes
Angular momentum: couple stresses
The differential energy equation
Boundary conditions
Stream function
Vorticity, irrotationality
Velocity potential
Plane potential flows
Incompressible viscous flows 4.1 – 4.8
4.9 – 4.25
4.26 – 4.37
4.38
4.39 – 4.42
4.43 – 4.46
4.47 – 4.55
4.56 – 4.60
4.61 – 4.67
4.68 – 4.78
4.79 – 4.91 the acceleration vector and (b) any unit vector normal to the
acceleration.
P4.2 Flow through the converging nozzle in Fig. P4.2 can be approximated by the onedimensional velocity distribution
u 2x
L V0 1 0 w 0 (a) Find a general expression for the fluid acceleration in
the nozzle. (b) For the specific case V0 10 ft/s and L
6 in, compute the acceleration, in g’s, at the entrance and at
the exit. V0 u = 3V0 x=L
x P4.2 x=0 P4.3 A twodimensional velocity field is given by
V P4.1 An idealized velocity field is given by the formula
V 4txi 2t 2yj 4xzk  v v Is this flow field steady or unsteady? Is it two or threedimensional? At the point (x, y, z) ( 1, 1, 0), compute (a)  eText Main Menu  (x2 y2 x)i (2xy y)j in arbitrary units. At (x, y) (1, 2), compute (a) the accelerations ax and ay, (b) the velocity component in the direction
40°, (c) the direction of maximum velocity, and (d)
the direction of maximum acceleration. Textbook Table of Contents  Study Guide Problems 265
P4.4 Suppose that the temperature field T 4x2 3y3, in arbitrary units, is associated with the velocity field of Prob. 4.3.
Compute the rate of change dT/dt at (x, y) (2, 1).
P4.5 The velocity field near a stagnation point (see Example 1.10)
may be written in the form
U0x
L u U0y
L U0 and L are constants (a) Show that the acceleration vector is purely radial. (b)
For the particular case L 1.5 m, if the acceleration at (x,
y) (1 m, 1 m) is 25 m/s2, what is the value of U0?
P4.6 Assume that flow in the converging nozzle of Fig. P4.2 has
the form V V0[1 (2x)/L]i. Compute (a) the fluid acceleration at x L and (b) the time required for a fluid particle to travel from x 0 to x L.
P4.7 Consider a sphere of radius R immersed in a uniform stream
U0, as shown in Fig. P4.7. According to the theory of Chap.
8, the fluid velocity along streamline AB is given by
V U0 1 ui P4.9 P4.10 P4.11 P4.12 the time for which the fluid acceleration at x L is zero.
Why does the fluid acceleration become negative after condition (b)?
A velocity field is given by V (3y2 3x2)i Cxyj 0k.
Determine the value of the constant C if the flow is to be
(a) incompressible and (b) irrotational.
Write the special cases of the equation of continuity for (a)
steady compressible flow in the yz plane, (b) unsteady incompressible flow in the xz plane, (c) unsteady compressible flow in the y direction only, (d) steady compressible
flow in plane polar coordinates.
Derive Eq. (4.12b) for cylindrical coordinates by considering the flux of an incompressible fluid in and out of the elemental control volume in Fig. 4.2.
Spherical polar coordinates (r, , ) are defined in Fig.
P4.12. The cartesian transformations are
x cos y r sin sin z R3
i
x3 r sin r cos Find (a) the position of maximum fluid acceleration along
AB and (b) the time required for a fluid particle to travel
from A to B. z υr
υφ y U0 P
B Sphere A
x = – 4R θ x r R υθ P4.7 y φ P4.8 When a valve is opened, fluid flows in the expansion duct
of Fig. 4.8 according to the approximation
iU 1 V x P4.12 x
Ut
tanh
2L
L Find (a) the fluid acceleration at (x, t) The cartesian incompressible continuity relation (4.12a) can
be transformed to the spherical polar form (L, L/U) and (b)
1
r2 r (r2 r) 1
r sin ( x=0 v   1
r sin Ky u
P4.8 sin ) () 0 What is the most general form of r when the flow is purely
radial, that is, and
are zero?
P4.13 A twodimensional velocity field is given by u (x, t) v r = constant x=L eText Main Menu  x2 Kx
y2 x2 y2 where K is constant. Does this field satisfy incompressible Textbook Table of Contents  Study Guide 266 Chapter 4 Differential Relations for a Fluid Particle P4.14 P4.15 P4.16 P4.17 u U continuity? Transform these velocities to polar components
. What might the flow represent?
r and
For incompressible polarcoordinate flow, what is the most
general form of a purely circulatory motion,
(r, , t)
and r 0, which satisfies continuity?
What is the most general form of a purely radial polarcoordinate incompressibleflow pattern, r
r(r, , t) and
0, which satisfies continuity?
An incompressible steadyflow pattern is given by u x3
2z2 and w y3 2yz. What is the most general form of the
third component, (x, y, z), which satisfies continuity?
A reasonable approximation for the twodimensional incompressible laminar boundary layer on the flat surface in
Fig. P4.17 is
y2 2y for y 2 Cx1/2, C where const (a) Assuming a noslip condition at the wall, find an expression for the velocity component (x, y) for y
. (b)
Then find the maximum value of at the station x 1 m,
for the particular case of airflow, when U 3 m/s and
1.1 cm.
Layer thickness δ (x)
U
U = constant y
U u (x, y) u (x, y) x 0 P4.17
P4.18 A piston compresses gas in a cylinder by moving at constant
speed , as in Fig. P4.18. Let the gas density and length at
t 0 be 0 and L0, respectively. Let the gas velocity vary linearly from u V at the piston face to u 0 at x L. If the
gas density varies only with time, find an expression for (t). V = constant u (x, t) ρ (t) P4.20 A twodimensional incompressible velocity field has u
K(1 e ay), for x L and 0 y
. What is the most
general form of (x, y) for which continuity is satisfied and
0? What are the proper dimensions for con0 at y
stants K and a?
P4.21 Air flows under steady, approximately onedimensional
conditions through the conical nozzle in Fig. P4.21. If the
speed of sound is approximately 340 m/s, what is the minimum nozzlediameter ratio De /D0 for which we can safely
neglect compressibility effects if V0 (a) 10 m/s and (b)
30 m/s? V0 Ve De P4.21 D0 P4.22 Air at a certain temperature and pressure flows through a
contracting nozzle of length L whose area decreases linearly,
A A0[1 x/(2L)]. The air average velocity increases
nearly linearly from 76 m/s at x 0 to 167 m/s at x L. If
the density at x 0 is 2.0 kg/m3, estimate the density at
x L.
P4.23 A tank volume contains gas at conditions ( 0, p0, T0). At
time t 0 it is punctured by a small hole of area A. According to the theory of Chap. 9, the mass flow out of such
a hole is approximately proportional to A and to the tank
pressure. If the tank temperature is assumed constant and
the gas is ideal, find an expression for the variation of density within the tank.
*P4.24 Reconsider Fig. P4.17 in the following general way. It is
known that the boundary layer thickness (x) increases monotonically and that there is no slip at the wall (y 0). Further, u(x, y) merges smoothly with the outer stream flow,
where u U constant outside the layer. Use these facts
to prove that (a) the component (x, y) is positive everywhere within the layer, (b) increases parabolically with y
very near the wall, and (c) is a maximum at y
.
P4.25 An incompressible flow in polar coordinates is given by x P4.18 x=0 r x = L (t)  v v P4.19 An incompressible flow field has the cylindrical components
Cr, z K(R2 r2), r 0, where C and K are constants and r R, z L. Does this flow satisfy continuity?
What might it represent physically?  eText Main Menu  K cos
K sin 1
1 b
r2
b
r2 Does this field satisfy continuity? For consistency, what Textbook Table of Contents  Study Guide Problems 267
should the dimensions of constants K and b be? Sketch the
surface where r 0 and interpret.
*P4.26 Curvilinear, or streamline, coordinates are defined in Fig.
P4.26, where n is normal to the streamline in the plane of
the radius of curvature R. Show that Euler’s frictionless momentum equation (4.36) in streamline coordinates becomes Stagnation
point
(u = 0) y a x 0 V
t V
V
s 1 p
s gs (1) V V2
R 1 p
n gn (2) t P4.31 Further show that the integral of Eq. (1) with respect to s is
none other than our old friend Bernoulli’s equation (3.76).
n s, V P4.32
z θ
y Streamline
R x P4.33 P4.26
P4.27 A frictionless, incompressible steadyflow field is given by
V P4.29 P4.30 P4.31
EES P4.34 y2j in arbitrary units. Let the density be 0 constant and neglect gravity. Find an expression for the pressure gradient
in the x direction.
If z is “up,” what are the conditions on constants a and b
for which the velocity field u ay,
bx, w 0 is an exact solution to the continuity and NavierStokes equations
for incompressible flow?
Consider a steady, twodimensional, incompressible flow of
a newtonian fluid in which the velocity field is known, i.e.,
u
2xy,
y2 x2, w 0. (a) Does this flow satisfy
conservation of mass? (b) Find the pressure field, p(x, y) if
the pressure at the point (x 0, y 0) is equal to pa.
Show that the twodimensional flow field of Example 1.10
is an exact solution to the incompressible NavierStokes
equations (4.38). Neglecting gravity, compute the pressure
field p(x, y) and relate it to the absolute velocity V 2 u2
2
. Interpret the result.
According to potential theory (Chap. 8) for the flow approaching a rounded twodimensional body, as in Fig. P4.31,
the velocity approaching the stagnation point is given by
u U(1 a2/x2), where a is the nose radius and U is the
velocity far upstream. Compute the value and position of
the maximum viscous normal stress along this streamline.  v v P4.28 2xyi  eText Main Menu  P4.35 Is this also the position of maximum fluid deceleration?
Evaluate the maximum viscous normal stress if the fluid is
SAE 30 oil at 20°C, with U 2 m/s and a 6 cm.
The answer to Prob. 4.14 is
f(r) only. Do not reveal
this to your friends if they are still working on Prob. 4.14.
Show that this flow field is an exact solution to the NavierStokes equations (4.38) for only two special cases of the
function f(r). Neglect gravity. Interpret these two cases physically.
From Prob. 4.15 the purely radial polarcoordinate flow
which satisfies continuity is r f( )/r, where f is an arbitrary function. Determine what particular forms of f( ) satisfy the full NavierStokes equations in polarcoordinate
form from Eqs. (D.5) and (D.6).
The fully developed laminarpipeflow solution of Prob.
3.53, z umax(1 r2/R2),
0, r 0, is an exact solution to the cylindrical NavierStokes equations (App. D).
Neglecting gravity, compute the pressure distribution in the
pipe p(r, z) and the shearstress distribution (r, z), using R,
umax, and as parameters. Why does the maximum shear
occur at the wall? Why does the density not appear as a parameter?
From the NavierStokes equations for incompressible flow
in polar coordinates (App. D for cylindrical coordinates),
find the most general case of purely circulating motion (r),
0, for flow with no slip between two fixed conr
z
centric cylinders, as in Fig. P4.35. υθ (r)
r No slip r=b P4.35 Textbook Table of Contents r=a  Study Guide 268 Chapter 4 Differential Relations for a Fluid Particle P4.36 A constantthickness film of viscous liquid flows in laminar motion down a plate inclined at angle , as in Fig. P4.36.
The velocity profile is
u Cy(2h y) w Tw
y=h 0 u(y) y Find the constant C in terms of the specific weight and viscosity and the angle . Find the volume flux Q per unit width
in terms of these parameters. x y=0
Tw P4.40 y
g 4umaxy(h
h2 u u(y) θ u x P4.36 *P4.37 A viscous liquid of constant and falls due to gravity between two plates a distance 2h apart, as in Fig. P4.37. The
flow is fully developed, with a single velocity component
w w(x). There are no applied pressure gradients, only
gravity. Solve the NavierStokes equation for the velocity
profile between the plates.
P4.42
h
x z, w P4.43
P4.37 v v P4.38 Reconsider the angularmomentum balance of Fig. 4.5 by
adding a concentrated body couple Cz about the z axis [6].
Determine a relation between the body couple and shear
stress for equilibrium. What are the proper dimensions for
Cz? (Body couples are important in continuous media with
microstructure, such as granular materials.)
P4.39 Problems involving viscous dissipation of energy are dependent on viscosity , thermal conductivity k, stream velocity U0,
and stream temperature T0. Group these parameters into the dimensionless Brinkman number, which is proportional to .
P4.40 As mentioned in Sec. 4.11, the velocity profile for laminar
flow between two plates, as in Fig. P4.40, is   y) w 0 If the wall temperature is Tw at both walls, use the incompressibleflow energy equation (4.75) to solve for the temperature distribution T(y) between the walls for steady flow.
*P4.41 The approximate velocity profile in Prob. 3.18 and Fig. P3.18
for steady laminar flow through a duct, was suggested as h h T(y) eText Main Menu  P4.44 P4.45 P4.46 y2
b2 umax 1 1 z2
h2 With
w 0, it satisfied the noslip condition and gave
a reasonable volumeflow estimate (which was the point of
Prob. 3.18). Show, however, that it does not satsify the xmomentum NavierStokes equation for duct flow with constant pressure gradient p/ x 0. For extra credit, explain
briefly how the actual exact solution to this problem is obtained (see, for example, Ref. 5, pp. 120 – 121).
In ductflow problems with heat transfer, one often defines
an average fluid temperature. Consider the duct flow of Fig.
P4.40 of width b into the paper. Using a controlvolume integral analysis with constant density and specific heat, derive an expression for the temperature arising if the entire
duct flow poured into a bucket and was stirred uniformly.
Assume arbitrary u(y) and T(y). This average is called the
cupmixing temperature of the flow.
For the draining liquid film of Fig. P4.36, what are the appropriate boundary conditions (a) at the bottom y 0 and
(b) at the surface y h?
Suppose that we wish to analyze the sudden pipeexpansion
flow of Fig. P3.59, using the full continuity and NavierStokes equations. What are the proper boundary conditions
to handle this problem?
Suppose that we wish to analyze the Utube oscillation flow
of Fig. P3.96, using the full continuity and NavierStokes
equations. What are the proper boundary conditions to handle this problem?
Fluid from a large reservoir at temperature T0 flows into a
circular pipe of radius R. The pipe walls are wound with an
electricresistance coil which delivers heat to the fluid at a
rate qw (energy per unit wall area). If we wish to analyze Textbook Table of Contents  Study Guide Problems 269 P4.47 P4.48 P4.49 P4.50 P4.51 P4.52 this problem by using the full continuity, NavierStokes, and
(2L, 0, 0), (2L, 0, b), (0, L, b), and (0, L, 0). Show the dienergy equations, what are the proper boundary conditions
rection of Q.
for the analysis?
*P4.55 In spherical polar coordinates, as in Fig. P4.12, the flow is
A twodimensional incompressible flow is given by the vecalled axisymmetric if
0 and /
0, so that r
locity field V 3yi 2xj, in arbitrary units. Does this flow
(r, ). Show that a stream function (r,
r(r, ) and
satisfy continuity? If so, find the stream function (x, y)
) exists for this case and is given by
and plot a few streamlines, with arrows.
1
1
Determine the incompressible twodimensional stream
r
r
r2 sin
r sin
function (x, y) which represents the flow field given in Example 1.10.
This is called the Stokes stream function [5, p. 204].
Investigate the stream function
K(x2 y2), K con P4.56 Investigate the velocity potential
Kxy, K constant.
stant. Plot the streamlines in the full xy plane, find any stagSketch the potential lines in the full xy plane, find any stagnation points, and interpret what the flow could represent.
nation points, and sketch in by eye the orthogonal streamInvestigate the polarcoordinate stream function
lines. What could the flow represent?
Kr1/2 sin 1 , K constant. Plot the streamlines in the full P4.57 Determine the incompressible twodimensional velocity po2
xy plane, find any stagnation points, and interpret.
tential (x, y) which represents the flow field given in ExInvestigate the polarcoordinate stream function
ample 1.10. Sketch a few potential lines and streamlines.
Kr2/3 sin (2 /3), K constant. Plot the streamlines in all ex P4.58 Show that the incompressible velocity potential in plane pocept the bottom right quadrant, and interpret.
lar coordinates (r, ) is such that
A twodimensional, incompressible, frictionless fluid is
1
guided by wedgeshaped walls into a small slot at the orir
r
r
gin, as in Fig. P4.52. The width into the paper is b, and the
volume flow rate is Q. At any given distance r from the slot,
Further show that the angular velocity about the zaxis in
the flow is radial inward, with constant velocity. Find an exsuch a flow would be given by
pression for the polarcoordinate stream function of this flow.
1
1
2z
(r )
( r)
rr
r
θ = π /4 vr
r
Slot θ=0 P4.52 as defined above satisfies Laplace’s
Finally show that
equation in polar coordinates for incompressible flow.
P4.59 Consider the simple flow defined by V xi yj, in arbitrary units. At t 0, consider the rectangular fluid element
defined by the lines x 2, x 3 and y 2, y 3. Determine, and draw to scale, the location of this fluid element
at t 0.5 unit. Relate this new element shape to whether
the flow is irrotational or incompressible.
P4.60 Liquid drains from a small hole in a tank, as shown in Fig.
P4.60, such that the velocity field set up is given by r 0,
0,
R2/r, where z H is the depth of the water
z P4.53 For the fully developed laminarpipeflow solution of Prob.
4.34, find the axisymmetric stream function (r, z). Use this
result to determine the average velocity V Q/A in the pipe
as a ratio of umax.
P4.54 An incompressible stream function is defined by
(x, y) U
(3x2y
L2 v v  patm
r
z=H y3) zC? where U and L are (positive) constants. Where in this chapter are the streamlines of this flow plotted? Use this stream
function to find the volume flow Q passing through the rectangular surface whose corners are defined by (x, y, z)  z eText Main Menu  z=0
r=R P4.60 Textbook Table of Contents  Study Guide 270 P4.61 P4.62 P4.63 P4.64 Chapter 4 Differential Relations for a Fluid Particle
far from the hole. Is this flow pattern rotational or irrotational? Find the depth zC of the water at the radius r R.
Investigate the polarcoordinate velocity potential
1
Kr1/2 cos 2 , K constant. Plot the potential lines in the
full xy plane, sketch in by eye the orthogonal streamlines,
and interpret.
Show that the linear Couette flow between plates in Fig.
1.6 has a stream function but no velocity potential. Why
is this so?
Find the twodimensional velocity potential (r, ) for the
polarcoordinate flow pattern r Q/r,
K/r, where Q
and K are constants.
Show that the velocity potential (r, z) in axisymmetric
cylindrical coordinates (see Fig. 4.2) is defined such that
r z r z Further show that for incompressible flow this potential
satisfies Laplace’s equation in (r, z) coordinates.
P4.65 A twodimensional incompressible flow is defined by
Ky u x2 Kx
y2 x2 y2 A doublet is formed in the limit as a goes to zero (the source
and sink come together) while at the same time their
strengths m and m go to infinity and minus infinity, respectively, with the product a m remaining constant. (a)
Find the limiting value of velocity potential for the doublet.
Hint: Expand the natural logarithm as an infinite series of
the form
ln K K ln r C and K const Find the velocity potential for this flow. Sketch some streamlines and potential lines, and interpret the flow pattern.
P4.68 Find the stream function and plot some streamlines for the
combination of a line source m at (x, y) (0, a) and an
equal line source placed at (0, a).
P4.69 Find the stream function and plot some streamlines for the
combination of a counterclockwise line vortex K at (x, y)
( a, 0) and an equal line vortex placed at ( a, 0).
*P4.70 Superposition of a source of strength m at ( a, 0) and a
sink (source of strength m) at (a, 0) was discussed briefly
in this chapter, where it was shown that the velocity potential function is  v v (x
1
m ln
(x
2  a)2
a)2 3 Manifold
Water
8m 25 cm/s const Find the stream function for this flow, sketch some streamlines and potential lines, and interpret the flow pattern.
P4.67 A stream function for a plane, irrotational, polarcoordinate flow is
C 2 as goes to zero. (b) Rewrite your result for doublet in cylindrical coordinates.
P4.71 Find the stream function and plot some streamlines for the
combination of a counterclockwise line vortex K at (x, y)
( a, 0) and an opposite (clockwise) line vortex placed at
( a, 0).
P4.72 A coastal power plant takes in cooling water through a vertical perforated manifold, as in Fig. P4.72. The total volume
flow intake is 110 m3/s. Currents of 25 cm/s flow past the
manifold, as shown. Estimate (a) how far downstream and
(b) how far normal to the paper the effects of the intake are
felt in the ambient 8mdeep waters. where K constant. Is this flow irrotational? If so, find
its velocity potential, sketch a few potential lines, and interpret the flow pattern.
P4.66 A plane polarcoordinate velocity potential is defined by
K cos
r 3 1
1 y2
y2 eText Main Menu  110 m3/s P4.72 P4.73 A twodimensional Rankine halfbody, 8 cm thick, is placed
in a water tunnel at 20°C. The water pressure far upstream
along the body centerline is 120 kPa. What is the nose radius of the halfbody? At what tunnel flow velocity will cavitation bubbles begin to form on the surface of the body?
P4.74 Find the stream function and plot some streamlines for the
combination of a uniform stream iU and a clockwise line
vortex K at the origin. Are there any stagnation points in
the flow field?
*P4.75 Find the stream function and plot some streamlines for the
combination of a line source 2m at (x, y) ( a, 0) and a
line source m at ( a, 0). Are there any stagnation points in
the flow field?
P4.76 Air flows at 1.2 m/s along a flat surface when it encounters
a jet of air issuing from the horizontal wall at point A, as in
Fig. 4.76. The jet volume flow is 0.4 m3/s per unit depth
into the paper. If the jet is approximated as an inviscid line
source, (a) locate the stagnation point S on the wall. (b) How
far vertically will the jet flow extend into the stream? Textbook Table of Contents  Study Guide Problems 271
0.4 m3/(s • m) 1.2 m/s
G
S
A P4.76
P4.77 A tornado is simulated by a line sink m
1000 m2/s plus
a line vortex K
1600 m2/s. Find the angle between any
streamline and a radial line, and show that it is independent
of both r and . If this tornado forms in sealevel standard
air, at what radius will the local pressure be equivalent to
29 inHg?
P4.78 The solution to Prob. 4.68 (do not reveal!) can represent a
line source m at (0, a) near a horizontal wall (y 0). [The
other source at (0, a) represents an “image” to create the
wall.] Find (a) the magnitude of the maxinun flow velocity
along the wall and (b) the point of minimum pressure along
the wall. Hint: Use Bernoulli’s equation.
*P4.79 Study the combined effect of the two viscous flows in Fig.
4.16. That is, find u(y) when the upper plate moves at speed
V and there is also a constant pressure gradient (dp/dx). Is
superposition possible? If so, explain why. Plot representative velocity profiles for (a) zero, (b) positive, and (c) negative pressure gradients for the same upperwall speed V.
*P4.80 Oil, of density and viscosity , drains steadily down the
side of a vertical plate, as in Fig. P4.80. After a development region near the top of the plate, the oil film will become independent of z and of constant thickness . Assume
that w w(x) only and that the atmosphere offers no shear
resistance to the surface of the film. (a) Solve the NavierStokes equation for w(x), and sketch its approximate shape.
(b) Suppose that film thickness and the slope of the velocity profile at the wall [ w/ x]wall are measured with a
laserDoppler anemometer (Chap. 6). Find an expression for
oil viscosity as a function of ( , , g, [ w/ x]wall). P4.81 Modify the analysis of Fig. 4.17 to find the velocity u when
the inner cylinder is fixed and the outer cylinder rotates at
angular velocity 0. May this solution be added to Eq.
(4.146) to represent the flow caused when both inner and
outer cylinders rotate? Explain your conclusion.
*P4.82 A solid circular cylinder of radius R rotates at angular velocity in a viscous incompressible fluid which is at rest
far from the cylinder, as in Fig. P4.82. Make simplifying
assumptions and derive the governing differential equation
and boundary conditions for the velocity field in the fluid.
Do not solve unless you are obsessed with this problem.
What is the steadystate flow field for this problem?
υθ (r, θ , t)
r θ
Ω
r=R P4.82 P4.83 The flow pattern in bearing lubrication can be illustrated by
Fig. P4.83, where a viscous oil ( , ) is forced into the gap
h(x) between a fixed slipper block and a wall moving at velocity U. If the gap is thin, h L, it can be shown that the
pressure and velocity distributions are of the form p p(x),
u u(y),
w 0. Neglecting gravity, reduce the NavierStokes equations (4.38) to a single differential equation for
u(y). What are the proper boundary conditions? Integrate and
show that
u 1 dp 2
(y
2 dx yh) U1 y
h where h h(x) may be an arbitrary slowly varying gap
width. (For further information on lubrication theory, see
Ref. 16.)
y Plate
Oil film Oil
inlet Fixed slipper
block Air h0 g h (x)
x z v v   eText Main Menu P4.83  Textbook Table of Contents u ( y) h1
U Moving wall
x P4.80 Oil
outlet  Study Guide 272 Chapter 4 Differential Relations for a Fluid Particle *P4.84 Consider a viscous film of liquid draining uniformly down
the side of a vertical rod of radius a, as in Fig. P4.84. At
some distance down the rod the film will approach a terminal or fully developed draining flow of constant outer radius b, with z
0. Assume that the atz(r),
r
mosphere offers no shear resistance to the film motion.
Derive a differential equation for z, state the proper boundary conditions, and solve for the film velocity distribution.
How does the film radius b relate to the total film volume
flow rate Q? z
pa
µa ≈ 0 Fully
developed
region a
b Film µ 6 cm Mercury 1m P4.86 υz ρ 8 mm Q P4.87 Suppose in Fig. 4.17 that neither cylinder is rotating. The
fluid has constant ( , , k, cp). What, then, is the steadyflow solution for (r)? For this condition, suppose that the
inner and outer cylinder surface temperatures are Ti and To,
respectively. Simplify the differential energy equation appropriately for this problem, state the boundary conditions,
and find the temperature distribution in the fluid. Neglect
gravity.
P4.88 The viscous oil in Fig. P4.88 is set into steady motion by a
concentric inner cylinder moving axially at velocity U inside a fixed outer cylinder. Assuming constant pressure and
density and a purely axial fluid motion, solve Eqs. (4.38)
for the fluid velocity distribution z(r). What are the proper
boundary conditions? r Q SAE 10
oil Fixed outer cylinder P4.84
vz(r)
a
U Oil:ρ, µ vz P4.88 y
Incompressible
viscous
fluid b r P4.85 A flat plate of essentially infinite width and breadth oscillates sinusoidally in its own plane beneath a viscous fluid,
as in Fig. P4.85. The fluid is at rest far above the plate. Making as many simplifying assumptions as you can, set up the
governing differential equation and boundary conditions for
finding the velocity field u in the fluid. Do not solve (if you
can solve it immediately, you might be able to get exempted
from the balance of this course with credit).  v v *P4.89 Modify Prob. 4.88 so that the outer cylinder also moves to
the left at constant speed V. Find the velocity distribution
z(r). For what ratio V/U will the wall shear stress be the
same at both cylinder surfaces?
P4.90 A 5cmdiameter rod is pulled steadily at 2 m/s through a
x
fixed cylinder whose clearance is filled with SAE 10 oil at
Plate velocity:
20°C, as in Fig. P4.90. Estimate the (steady) force required
to pull the inner rod.
P4.85
U0 sin ω t
*P4.91 Consider twodimensional, incompressible, steady Couette
flow (flow between two infinite parallel plates with the upP4.86 SAE 10 oil at 20°C flows between parallel plates 8 cm apart,
per plate moving at constant speed and the lower plate staas in Fig. P4.86. A mercury manometer, with wall pressure
EES
tionary, as in Fig. 4.16a). Let the fluid be nonnewtonian,
taps 1 m apart, registers a 6cm height, as shown. Estimate
with its viscous stresses given by
the flow rate of oil for this condition.
u (x, y, z, t)?  eText Main Menu  Textbook Table of Contents  Study Guide Fundamentals of Engineering Exam Problems xy yx u
x a xx 1
2 a c c a yy y c u
y xz x
yz zy 1
2 a zx z 1
2 w
y a c w
z a zz u
z w
x 273 1.5 m Fixed cylinder c SAE 10 oil 2 cm c 5 cm where a and c are constants of the fluid. Make all the same
assumptions as in the derivation of Eq. (4.140). (a) Find the
velocity profile u(y). (b) How does the velocity profile for
this case compare to that of a newtonian fluid? Moving rod 2 m/s 2 cm P4.90 Word Problems
W4.1 The total acceleration of a fluid particle is given by Eq.
(4.2) in the eulerian system, where V is a known function
of space and time. Explain how we might evaluate particle acceleration in the lagrangian frame, where particle position r is a known function of time and initial position,
r fcn(r0, t). Can you give an illustrative example?
W4.2 Is it true that the continuity relation, Eq. (4.6), is valid for
both viscous and inviscid, newtonian and nonnewtonian,
compressible and incompressible flow? If so, are there any
limitations on this equation?
W4.3 Consider a CD compact disk rotating at angular velocity
. Does it have vorticity in the sense of this chapter? If
so, how much vorticity?
W4.4 How much acceleration can fluids endure? Are fluids like
astronauts, who feel that 5g is severe? Perhaps use the flow
pattern of Example 4.8, at r R, to make some estimates
of fluidacceleration magnitudes.
W4.5 State the conditions (there are more than one) under which
the analysis of temperature distribution in a flow field can
be completely uncoupled, so that a separate analysis for
velocity and pressure is possible. Can we do this for both
laminar and turbulent flow? W4.6 Consider liquid flow over a dam or weir. How might the
boundary conditions and the flow pattern change when we
compare water flow over a large prototype to SAE 30 oil
flow over a tiny scale model?
W4.7 What is the difference between the stream function and
our method of finding the streamlines from Sec. 1.9? Or
are they essentially the same?
W4.8 Under what conditions do both the stream function and
the velocity potential exist for a flow field? When does
one exist but not the other?
W4.9 How might the remarkable threedimensional Taylor instability of Fig. 4.18 be predicted? Discuss a general procedure for examining the stability of a given flow pattern.
W4.10 Consider an irrotational, incompressible, axisymmetric
(/
0) flow in (r, z) coordinates. Does a stream function exist? If so, does it satisfy Laplace’s equation? Are
lines of constant equal to the flow streamlines? Does a
velocity potential exist? If so, does it satisfy Laplace’s
equation? Are lines of constant everywhere perpendicular to the lines? Fundamentals of Engineering Exam Problems  v v This chapter is not a favorite of the people who prepare the FE
Exam. Probably not a single problem from this chapter will appear
on the exam, but if some did, they might be like these.
FE4.1 Given the steady, incompressible velocity distribution V
3xi Cyj 0k, where C is a constant, if conservation of
mass is satisfied, the value of C should be
(a) 3,
(b) 3/2,
(c) 0,
(d) 3/2,
(e) 3
FE4.2 Given the steady velocity distribution V 3xi 0j Cyk,  eText Main Menu  where C is a constant, if the flow is irrotational, the value
of C should be
(a) 3,
(b) 3/2,
(c) 0,
(d) 3/2,
(e) 3
FE4.3 Given the steady, incompressible velocity distribution V
3xi Cyj 0k, where C is a constant, the shear stress xy
at the point (x, y, z) is given by
(a) 3 ,
( b ) (3 x
Cy) ,
( c ) 0,
(d) C ,
(e) (3 C) Textbook Table of Contents  Study Guide 274 Chapter 4 Differential Relations for a Fluid Particle Comprehensive Problem
ˆ
Estimate the wall shear stress, w as a function of Q, f, , b,
s, and time (t). (e) Finally, for the numbers given in the problem statement, estimate the amplitude of the wall shear stress,
ˆw, in N/m2. C4.1 In a certain medical application, water at room temperature
and pressure flows through a rectangular channel of length
L 10 cm, width s 1.0 cm, and gap thickness b 0.30
mm as in Fig. C4.1. The volume flow rate is sinusoidal with
ˆ
amplitude Q
0.50 mL/s and frequency f 20 Hz, i.e.,
ˆ
Q Q sin (2 f t).
(a) Calculate the maximum Reynolds number (Re Vb/ )
based on maximum average velocity and gap thickness.
Channel flow like this remains laminar for Re less than about
2000. If Re is greater than about 2000, the flow will be turbulent. Is this flow laminar or turbulent? (b) In this problem,
the frequency is low enough that at any given time, the flow
can be solved as if it were steady at the given flow rate. (This
is called a quasisteady assumption.) At any arbitrary instant
of time, find an expression for streamwise velocity u as a
function of y, , dp/dx, and b, where dp/dx is the pressure
gradient required to push the flow through the channel at volume flow rate Q. In addition, estimate the maximum magnitude of velocity component u. (c) At any instant of time,
find a relationship between volume flow rate Q and pressure
gradient dp/dx. Your answer should be given as an expression for Q as a function of dp/dx, s, b, and viscosity . (d) L
y x
z s Q
b C4.1 References 2.
3.
4.
5.
6.
7.
8.
9. D. A. Anderson, J. C. Tannehill, and R. H. Pletcher, Computational Fluid Mechanics and Heat Transfer, 2d ed.,
McGrawHill, New York, 1997.
G. B. Wallis, OneDimensional TwoPhase Flow, McGrawHill, New York, 1969.
S. M. Selby, CRC Handbook of Tables for Mathematics, 4th
ed., CRC Press Inc., Cleveland, OH, 1976.
H. Schlichting, Boundary Layer Theory, 7th ed., McGrawHill, New York, 1979.
F. M. White, Viscous Fluid Flow, 2d ed., McGrawHill, New
York, 1991.
L. E. Malvern, Introduction to Mechanics of a Continuous
Medium, PrenticeHall, Englewood Cliffs, NJ, 1969.
J. P. Holman, Heat Transfer, 8th ed., McGrawHill, New York,
1997.
W. M. Kays and M. E. Crawford, Convective Heat and Mass
Transfer, 3d ed., McGrawHill, New York, 1993.
G. K. Batchelor, An Introduction to Fluid Dynamics, Cambridge University Press, Cambridge, England, 1967.  v v 1.  eText Main Menu  10.
11.
12.
13. 14.
15.
16.
17. 18. L. Prandtl and O. G. Tietjens, Fundamentals of Hydro and
Aeromechanics, Dover, New York, 1957.
R. Aris, Vectors, Tensors, and the Basic Equations of Fluid
Mechanics, PrenticeHall, Englewood Cliffs, NJ, 1962.
G. E. Mase, Continuum Mechanics, Schaum’s Outline Series,
McGrawHill, New York, 1970.
W. F. Hughes and E. W. Gaylord, Basic Equations of Engineering Science, Schaum’s Outline Series, McGrawHill,
New York, 1964.
G. Astarita and G. Marrucci, Principles of NonNewtonian
Fluid Mechanics, McGrawHill, New York, 1974.
H. Lamb, Hydrodynamics, 6th ed., Dover, New York, 1945.
A. Szeri, Tribology: Friction, Lubrication, and Wear,
McGrawHill, New York, 1980.
G. I. Taylor, “Stability of a Viscous Liquid Contained between Two Rotating Cylinders,” Philos. Trans. Roy. Soc. London Ser. A, vol. 223, 1923, pp. 289 – 343.
E. L. Koschmieder, “Turbulent Taylor Vortex Flow,” J. Fluid
Mech., vol. 93, pt. 3, 1979, pp. 515 – 527. Textbook Table of Contents  Study Guide ...
View Full
Document
 Spring '08
 Sakar
 The Land

Click to edit the document details