Solutions Manual for Fluid Mechanics 5th Edition by Kundu

Full Download: [Link]

Fluid Mechanics, 5th Ed. Kundu, Cohen, and
Dowling

Exercise 1.1. Many centuries ago, a mariner poured 100 cm3 of water into the ocean. As time
passed, the action of currents, tides, and weather mixed the liquid uniformly throughout the
earth's oceans, lakes, and rivers. Ignoring salinity, estimate the probability that the next cup of
water you drink will contain at least one water molecule that was dumped by the mariner. Assess
your chances of ever drinking truly pristine water. [Some possibly useful facts: Mw for water is
18.0 kg per kg-mole, the radius of the earth is 6370 km, the mean depth of the oceans is
approximately 3.8 km and they cover 71% of the surface of the earth. One cup is ~240 ml.]

Solution 1.1. To get started, first list or determine the volumes involved:
d = volume of water dumped = 100 cm3, c = volume of a cup ≈ 240 cm3, and
V = volume of water in the oceans = 4R 2 D ,
where, R is the radius of the earth, D is the mean depth of the oceans, and  is the oceans'
coverage fraction. Here we've ignored the ocean volume occupied by salt and have assumed that
the oceans' depth is small compared to the earth's diameter. Putting in the numbers produces:

V  4 (6.37 106 m)2 (3.8 103 m)(0.71) 1.376 1018 m 3 .
For well-mixed oceans, the probability Po that any water molecule in the ocean came from the
dumped water is:
(100 cm3 of water)  d 1.0 104 m 3
 Po     7.27 1023 ,
(oceans' volume) V 1.376 10 m 18 3

Denote the probability that at least one molecule from the dumped water is part of your next cup
as P1 (this is the answer to the question). Without a lot of combinatorial analysis, P1 is not easy
to calculate directly. It is easier to proceed by determining the probability P2 that all the
moleculesin your cup are not from the dumped water. With these definitions, P can be
1
determined from: P1 = 1 – P2. Here, we can calculate P2 from:
P2 = (the probability that a molecule was not in the dumped water)[number of molecules in a cup].
The number of molecules, Nc, in one cup of water is
1.00g gmole molecules
N c  240cm 3  3
  6.023 10 23  8.03 10 24 molecules
cm 18.0g gmole
24
Thus, P2  (1 Po ) N c  (1 7.27 1023 ) 8.0310 . Unfortunately, electronic calculators and modern
computer math programs cannot evaluate this expression, so analytical techniques are required.
First, take the natural log of both sides, i.e.

ln(P2 )  Nc ln(1 Po )  8.031024 ln(1 7.27 1023 )

then expand the natural logarithm using ln(1–) ≈ – (the first term of a standard Taylor series
for   0 )
ln(P2 )  Nc  Po  8.031024  7.27 1023  584 ,

and exponentiate to find:
P2  e584 10254 ... (!)

Therefore, P1 = 1 – P2 is very-very close to unity, so there is a virtual certainty that the next cup

of water you drink will have at least one molecule in it from the 100 cm3 of water dumped many
years ago. So, if one considers the rate at which they themselves and everyone else on the planet

uses water it is essentially impossible to get a truly fresh cup to drink.

Full download all chapters instantly please go to Solutions Manual, Test Bank site: [Link]
 Fluid Mechanics, 5th Ed. Kundu, Cohen, and
Dowling

Exercise 1.2. An adult human expels approximately 500 ml of air with each breath during
ordinary breathing. Imagining that two people exchanged greetings (one breath each) many
centuries ago, and that their breath subsequently has been mixed uniformly throughout the
atmosphere, estimate the probability that the next breath you take will contain at least one air
molecule from that age-old verbal exchange. Assess your chances of ever getting a truly fresh
breath of air. For this problem, assume that air is composed of identical molecules having Mw =
29.0 kg per kg-mole and that the average atmospheric pressure on the surface of the earth is 100
kPa. Use 6370 km for the radius of the earth and 1.20 kg/m3 for the density of air at room
temperature and pressure.

Solution 1.2. To get started, first determine the masses involved.

m = mass of air in one breath = density x volume = 1.20kg /m 3 0.5 103 m 3  = 0.60 103 kg

M = mass of air in the atmosphere = 4R 2
 (z)dz
z 0


Here, R is the radius of the earth, z is the elevation above the surface of
the earth, and (z) is the
air density as function of elevation. From the law for static pressure in a gravitational field,
z 
dP dz  g , the surface pressure,
 Ps, on the earth is determined from Ps  P   (z)gdz so
z 0
P P
that: M  4R s   4 (6.37 10 6 m) 2 (10 5 Pa)  5.2 1018 kg .
2

g
 where the pressure (vacuum) in outer space = P∞ = 0, and g is assumed constant throughout the

atmosphere. For a well-mixed atmosphere, the probability Po that any molecule in the
atmosphere came from the age-old verbal exchange is

2  (mass of one breath) 2m 1.2 103 kg
Po     2.311022 ,
(mass of the whole atmosphere) M 5.2 1018 kg
where the factor of two comes from one breath for each person. Denote the probability that at
least one molecule from the age-old verbal exchange is part of your next breath as P1 (this is the
answer to the question). Without a lot of combinatorial analysis, P1 is not easy to calculate
 It is easier to proceed by determining the probability P that all the molecules in your
directly. 2
next breath are not from the age-old verbal exchange. With these definitions, P1 can be
determined from: P1 = 1 – P2. Here, we can calculate P2 from:
P2 = (the probability that a molecule was not in the verbal exchange)[number of molecules in a breath].
The number of molecules, Nb, involved in one breath is
0.6 103 kg 10 3 g molecules
Nb    6.023 10 23  1.25 10 22 molecules
29.0g /gmole kg gmole
22
Thus, P2  (1 Po ) N b  (1 2.311022 )1.2510 . Unfortunately, electronic calculators and modern
computer math programs cannot evaluate this expression, so analytical techniques are required.
First,
 take the natural log of both sides, i.e.
ln(P2 )  Nb ln(1 Po ) 1.25 1022 ln(1 2.311022 )

then expand the natural logarithm using ln(1–) ≈ – (the first term of a standard Taylor series
for   0 )
ln(P2 )  Nb  Po  1.25 1022  2.311022  2.89 ,

and exponentiate to find:


Fluid Mechanics, 5th Ed. Kundu, Cohen, and
Dowling

P2  e2.89  0.056 .
Therefore, P1 = 1 – P2 = 0.944 so there is a better than 94% chance that the next breath you take
will have at least one molecule in it from the age-old verbal exchange. So, if one considers how
often they themselves and everyone else breathes, it is essentially impossible to get a breath of
truly fresh air. 
 Fluid Mechanics, 5th Ed. Kundu, Cohen, and
Dowling

Exercise 1.3. In Cartesian coordinates, the Maxwell probability distribution, f(u) = f(u1,u2,u3), of
molecular velocities in a gas flow with average velocity U = (U1,U2,U3) is
 m 3 2  m 2 
f (u)    exp  u  U 
2kB T   2k B T 
where n is the number of gas molecules in volume V, m is the molecular mass, kB is Boltzmann’s
constant and T is the absolute temperature.
a) Verify that U is the average molecular velocity, and determine the standard deviations (1,
 12
 
2,3) of each component of U using  i   (ui  U i ) f (u)d u for i = 1, 2, and 3.
2 3

all u 
b) Using the molecular version of perfect gas law (1.21), determine n/V at room temperature T =
295 K and atmospheric pressure p = 101.3 kPa.
c) Determine n for volumes V = (10 m)3, 1 m3, and (0.1 m)3.

d) For the ith velocity component, the standard deviation of the average, a,i, over n molecules
isa,i =  i n when n >> 1. For an airflow at U = (1.0 ms–1, 0, 0), compute the relative
uncertainty, 2 a ,1 U1 , at the 95% confidence level for the average velocity for the three volumes
listed in part c).
e)For the conditions specified in parts b) and d), what is the smallest volume of gas that ensures

a relative uncertainty in U of less than one percent?

Solution 1.3. a) Use the given distribution, and the definition of an average:
 m 3 2 + + +  m 2 
uave   u f (u)d u  
3
    u exp  u  U d 3 u .
all u 2k B T  – – –  2k B T 
Consider the first component of u, and separate out the integrations in the "2" and "3" directions.
 m 3 2 + + +  m 
(u1 ) ave       u1 exp 
2k B T  – – –  2k B T
(u1  U1 ) 2  (u2  U 2 ) 2  (u3  U 3 ) 2 du1du2 du3

32
 m  +  m(u1  U1 ) 2  +  m(u2  U 2 ) 2  +  m(u3  U 3 ) 2 
    1u exp 
 
du 1  exp 
 du 2  exp  du3
2k B T  –  2k B T  –  2kB T  –  2k B T 
 The integrations in the "2" and "3" directions are equal to: 2kB T m , so
12

 m 1 2 +  m(u1  U1 ) 2 
 (u )     1u exp  du1
2k B T  –
1 ave
 2k B T 
The change of integration variable to   (u  U  )m 2k T  changes this integral to:
12
1 1 B

1  2k B T  
+ 12
1
(u
1 ) ave   
 
 –  m 
  U 1 
exp
   2
d  0 

U1   U1 ,


where the first term of the integrand is an odd function integrated on an even interval so its
contribution is zero. This procedure is readily repeated for the other directions to find (u2)ave =
U2, and (u3)ave = U3. Using the same simplifications and change of integration variables

produces:
 m 3 2 + + +  m 
1  
2
    (u1  U1 ) exp 
2k B T  – – –
2

 2k B T
(u1  U1 ) 2  (u2  U 2 ) 2  (u3  U 3 ) 2 du1du2 du3



 Fluid Mechanics, 5th Ed. Kundu, Cohen, and
Dowling

 m 1 2 +  m(u1  U1 ) 2  1 2kB T + 2

    
     exp d .
2 2
 (u U ) exp du 
2 k T
 B  –
1 1
 2k B T 
1
  m –

The final integral over  is:  2 , so the standard deviations of molecular speed are
1  kB T m   2   3 ,
12

 where the second two equalities follow from repeating this calculation for the second and third
directions. 
b) From (1.21), n V  p kB T  (101.2kPa) [1.3811023 J /K  295K]  2.487 1025 m3

c) From n/V from part b): n  2.487 1010 for V = 103 m3 = 10–15 m3
n  2.487 107 for V = 1.0 m3 = 10–18 m3
 n  2.487 104 for V = 0.001 m3 = 10–21 m3
d) From (1.22), the gas constant is R = (kB/m), and R = 287 m2/s2K for air. Compute:
1 2
2 a,1 U1  2kB T m n  1m / s  2RT n 1m /s  2287  295 n   582 n . Thus,
12 12

for V = 10–15 m3
: 2 a,1 U1 = 0.00369,
–18 3
V = 10 m : 2 a,1 U1 = 0.117, and
 V = 10–21 m3 : 2 a,1 U1 = 3.69.
e) To achieve
a relative uncertainty of 1% we need n ≈ (582/0.01)2 = 3.39  109, and this
corresponds to a volume of 1.36  10 m which is a cube with side dimension ≈ 5 m.
-16 3




 Fluid Mechanics, 5th Ed. Kundu, Cohen, and
Dowling

Exercise 1.4. Using the Maxwell molecular velocity distribution given in Exercise 1.3 with U =
 1 2
v   u f (u)d u
2
0, determine the average molecular speed = 3
and compare it with c = speed of
all u 
sound in a perfect gas under the same conditions.

Solution 1.4. Use the specified formfor v and the Maxwell distribution
 m 3 2 + + + 2  m 
v   u f (u)d u       u1  u2  u3 exp  u12  u22  u32 du1du2 du3 .
2 2 3 2 2

all u 2k B T  – – –  2k B T 
This can be re-arranged and expanded  into a total of nine one-variable integrations:
32
 m  +  mu12  +  mu22  +  mu32 
v 2     1 
  1    2   du3
2
u exp du exp du exp
 2k B T  –  2k B T  –  2k B T  –  2k B T 
 m 3 2 +  mu12  + 2  mu22  +  mu32 
   exp  du1  u2 exp  du2  exp  du3
2kB T  –  2k B T  –  2k B T  –  2k B T 
  m 3 2 +  mu12  +  mu22  + 2  mu32 
   exp  du1  exp  du2  u3 exp  du3 .
2kB T  –  2k B T  –  2k B T  –  2k B T 
arrangement, the six off-diagonal integrals are equal to 2k T m1 2 and the three on-
In this B

diagonal integrals are equal to 2kB T m

32

 2 . Thus, 

 m  2kB T 2k B T 3 2 2k B T 3 2 2k B T 3 2 
32
3k T
v  
2
           , or v 2  B .
2k B T   m  m   m   m   2 m
 so v  3RT and this speed has the same temperature dependence but is
From (1.22), R = (kB/m)
a factor of 3  larger than the speed of sound in a perfect gas: c  RT .
 

 
 Fluid Mechanics, 5th Ed. Kundu, Cohen, and
Dowling

Exercise 1.5. By considering the volume swept out by a moving molecule, estimate how the
mean-free path, l, depends on the average molecular cross section dimension d and the
molecular number density n˜ for nominally spherical molecules. Find a formula for ln˜ 1 3 (= the
ratio of the mean-free path to the mean intermolecular spacing) in terms of the molecular volume
( d 3 ) and the available volume per molecule ( 1 n˜ ). Is this ratio typically bigger
 or smaller than
one?  

 Solution 1.5. The combined collision cross section for two spherical molecules having diameter

d is d 2 . The mean free path l is the average distance traveled by a molecule between collisions.
Thus, the average molecule should experience one collision when
sweeping a volume equal to d 2l . If the molecular number density is
n˜ , then the volume per molecule is n˜ 1, and the mean intermolecular
  d
spacing is n˜ 1 3 . Assuming that the swept volume necessary to
produce one collision is proportional to the volume per molecule

produces:
 
d 2l  C n˜ or l  C n˜ d 2 , 

where C is a dimensionless constant presumed to be of order unity. The dimensionless version of
this equation is:
mean free path l
   1 3  l n˜ 1 3
mean intermolecular spacing n˜
C C n˜ 1 2 3 volume per molecule 2 3
   C 3   C  ,
n˜ 2 3d 2 n˜d 
3 23
d   molecular volume 

where all numerical constants like  have been combined into C. Under ordinary conditions in
not tightly packed so l  n˜ 1 3 . In liquids, the molecules are tightly
gases, the molecules are
packed so l ~ n˜ 1 3 .



 Fluid Mechanics, 5th Ed. Kundu, Cohen, and
Dowling

Exercise 1.6. In a gas, the molecular momentum flux (MFij) in the j-coordinate direction that
crosses a flat surface of unit area with coordinate normal direction i is:
n
MFij   mui u j f (u)d 3 u where f(u) is the Maxwell distribution given in Exercise 1.3. For a
V all u
perfect gas that is not moving on average (i.e. U = 0), show that MFij = p, the pressure, when i =
j, and that MFij = 0, when i ≠ j.
 Solution 1.6. Start from the given equation using the Maxwell distribution:
32
n nm  m  + + +  m 
MFij   mui u j f (u)d u 
V all u
3
     ui u j exp 
V 2kB T  – – –  2kB T
u12  u22  u32 du1du2 du3

and first consider i = j = 1, and recognize  = nm/V as the gas density (see (1.22)).
 m 3 2 + + + 2  m 

MF11       u1 exp 
2k B T  – – –  2k B T
u12  u22  u32 du1du2 du3

32
 m  + 2  mu12  +  mu22  +  mu32 
    u1 exp  du1  exp  du2  exp  du3
2kB T  –  2k B T  –  2k B T  –  2k B T 
 
The first integral is equal to 2k T m3 2  2 while the second two integrals are each equal to
B

2kB T m
12
. Thus:

 m 3 2 2k B T 3 2  2k B T 1 2 2k B T 1 2 k T
MF11            B  RT  p

2k B T   m  2  m   m  m
 where kB/m = R from (1.22). This analysis may be repeated with i = j = 2, and i = j = 3 to find:
MF22 = MF33 = p, as well.
Now consider the case i ≠ j. First note that MFij = MFji because the velocity product
under the triple integral may be written in either order u u = u u , so there are only three cases of
i j j i
interest. Start with i = 1, and j = 2 to find:
 m 3 2 + + +  m 
MF12       u1u2 exp 
2k B T  – – –  2k B T
u12  u22  u32 du1du2 du3

32
 m  +  mu1 
2 +  mu2  +
2  mu32 
    1 u exp  
du1  2
u exp  du 2  exp  du3
2kB T  –  2k B T  –  2k B T  –  2k B T 
 Here we need only consider the first integral. The integrand of this integral is an odd function
because it is product of an odd function, u1, and an even function, expmu12 2kB T. The
 integral of an odd function on an even interval [–∞,+∞] is zero, so MF12 = 0. And, this analysis
may be repeated for i = 1 and j = 3, and i = 2 and j = 3 to find MF13 = MF23 = 0.

Fluid Mechanics, 5th Ed. Kundu, Cohen, and
Dowling

Exercise 1.7. Consider the viscous flow in a channel of width 2b. The channel is aligned in the
x-direction, and the velocity u in the x-direction at a distance y from the channel centerline is
 
given by the parabolic distribution u(y) = U 0 1 y b . Calculate the shear stress  as a
2

function y, , b, and Uo. What is the shear stress at y = 0?

du d  y 2  y

Solution 1.7. Start from (1.3):      U o 1    –2U o 2 . At y = 0 (the location of
dy dy  b   b
  m2Uo b .
maximum velocity)  = 0. At At y = ±b (the locations of zero velocity), 



Fluid Mechanics, 5th Ed. Kundu, Cohen, and
Dowling

Exercise 1.8. Estimate the height to which water at 20 C will rise in a capillary glass tube 3 mm
in diameter that is exposed to the atmosphere. For water in contact with glass the wetting angle is
nearly 90. At 20 C, the surface tension of an water-air interface is  = 0.073 N/m. (Answer: h =
0.99 cm.)

Solution 1.8. Start from the result of Example 1.1.

2 sin  2(0.073N /m)sin(90)
h   9.92mm
gR (10 kg /m 3 )(9.81m /s2 )(1.5 103 m)
3


Fluid Mechanics, 5th Ed. Kundu, Cohen, and
Dowling

Exercise 1.9. A manometer is a U-shaped tube containing mercury of density m. Manometers
are used as pressure measuring devices. If the fluid in the tank A has a pressure p and density ,
then show that the gauge pressure in the tank is: p  patm = mgh  ga. Note that the last term on
the right side is negligible if  << m. (Hint: Equate the pressures at X and Y.)

Solution 1.9. Start by equating the pressures at X and Y.

pX = p + ga = patm + mgh = pY.
Rearrange to find:
p – patm = mgh – ga.
 Fluid Mechanics, 5th Ed. Kundu, Cohen, and
Dowling

Exercise 1.10. Prove that if e(T, ) = e(T) only and if h(T, p) = h(T) only, then the (thermal)
equation of state is (1.22) or p = kT.

Solution 1.10. Start with the first member of (1.18): de = Tds – pd, and rearrange it:
1 p s  s 
ds  de  d    de    d ,
T T e   e
where the second equality holds assuming the entropy depends on e and . Here we see that:
1 s  p s 
   , and    .
 T e  T  e
  s     s  
Equality of the crossed second derivatives of s,          , implies:
 e  e e  e 
  1T  p T 
     .
  e  e 
  1 T  1 T 
However, if e depends only on T, then (∂/∂)e = (∂/∂)T, thus       0 , so
  e   T
 p T  
   0 , and this can be integrated once to find: p/T = f1(), where f1 is an undetermined
 e 
function. 
Now repeat this procedure using the second member of (1.18), dh = Tds + dp.
1  s  s 
 ds  dh  dp    dh    dp .
T T h p p h
1 s   s 
Here equality of the coefficients of the differentials implies:    , and     .
T h p T p h
  1 T    T 
So, equality of the crossed second derivatives implies:      .
 p h  h p
  1 T  1 T 

Yet, if h depends only on T, then (∂/∂p)h = (∂/∂p)T, thus       0 , so
 p h  p T
  T  
   0 , and this can be integrated once to find: /T = f2(p), where f2 is an undetermined
 h p
function. 
Collecting the two results involving f1 and f2, and solving for T produces:
p 
 T  or pf 2 (p)  f1( )  k ,
f1 ( ) f 2 ( p)
where k must be is a constant since p and  are independent thermodynamic variables.
Eliminating f1 or f2 from either equation of the left, produces p = kT.
And finally, using both versions  of (1.18) we can write: dh – de = dp + pd = d(p).

When e and h only depend on T, then dh = CpdT and de = CvdT, so
dh – de = (Cp – Cv)dT = d(p) = kdT , thus k = Cp – Cv = R,
where R is the gas constant. Thus, the final result is the perfect gas law: p = kT/ = RT.
 Fluid Mechanics, 5th Ed. Kundu, Cohen, and
Dowling

Exercise 1.11. Starting from the property relationships (1.18) prove (1.25) and (1.26) for a
reversible adiabatic process when the specific heats Cp and Cv are constant.

Solution 1.11. For an isentropic process: de = Tds – pd = –pd, and dh = Tds + dp = +dp.
Equations (1.25) and (1.26) apply to a perfect gas so the definition of the specific heat capacities
(1.14), and (1.15) for a perfect gas, dh = CpdT, and de = CvdT , can be used to form the ratio
dh/de:
dh Cp dT Cp dp d d dp
    or    .
de Cv dT Cv pd   p
The final equality integrates to: ln(p) = ln() + const which can be exponentiated to find:
p = const.,
which is (1.25). The constant may be evaluated at a reference condition po and o to find:

p po   o  and this may be inverted to put the density ratio on the left
1
 o  p po  ,
which is the second member of (1.26). The remaining relationship involving the temperature is
 found by using the perfect gas law, p = RT, to eliminate  = p/RT:
1 1 
 p RT  pTo  p  T p  p   p ( 1) 
     or       ,
 o po RTo poT po  To po po  po 
which is the first member of (1.26).

 
Fluid Mechanics, 5th Ed. Kundu, Cohen, and
Dowling

Exercise 1.12. A cylinder contains 2 kg of air at 50 C and a pressure of 3 bars. The air is
compressed until its pressure rises to 8 bars. What is the initial volume? Find the final volume for
both isothermal compression and isentropic compression.

Solution 1.12. Use the perfect gas law but explicitly separate the mass M of the air and the
volume V it occupies via the substitution  = M/V:
p = RT = (M/V)RT.
Solve for V at the initial time:
Vi = initial volume = MRT/pi = (2 kg)(287 m2/s2K)(273 + 50°)/(300 kPa) = 0.618 m3.
For an isothermal process:
Vf = final volume = MRT/pf = (2 kg)(287 m2/s2K)(273 + 50°)/(800 kPa) = 0.232 m3.
For an isentropic process:
V f  Vi pi p f  0.618m 3 300kPa 800kPa  0.307m 3 .
1 1 1.4


 Fluid Mechanics, 5th Ed. Kundu, Cohen, and
Dowling

Exercise 1.13. Derive (1.29) starting from the arguments provided at the beginning of Section
1.10 and Figure 1.8.

Solution 1.13. Take the z axis vertical, and consider a small fluid element m of fluid having
volume V that starts at height z0 in a stratified fluid medium having a vertical density profile =
(z), and a vertical pressure profile p(z). Without any vertical displacement, the small mass and
its volume are related by m = (z0)V. If the small mass is displaced vertically a small distance
 via an isentropic process, its density will change isentropically according to:
a (z0   )  (z0 )  da dz  ...
where da/dz is the isentropic density at z0. For a constant m, the volume of the fluid element
will be:
m m m  1 da 
V    1   ...
a (z0 )  da dz  ... (z0 )  (z0 ) dz 
The background density at z0 +  is:
(z0   )  (z0 )  d dz  ...
If g is the acceleration of gravity, the (upward) buoyant force on the element at the vertically

displaced location will be g(z0 + )V, while the (downward) weight of the fluid element at any
vertical location is gm. Thus, a vertical application Newton's second law implies:

d 2 m  1 da 
m 2  g(z0   )V  gm  g(z0 )  d dz  ... 1   ... gm ,
dt (z0 )  (z0 ) dz 
where the second equality follows from substituting for (z0 + ) and V from the above
equations. Multiplying out the terms in (,)-parentheses and dropping second order terms
produces:

d 2 gm d gm da gm d da 
m 2  gm     ... gm    
dt (z0 ) dz (z0 ) dz (z0 ) dz dz 
Dividing by m and moving all the terms to the right side of the equation produces:
d 2 g d da 
     0
 dt 2 (z0 ) dz dz 
Thus, for oscillatory motion at frequency N, we must have
g d da 
N2     ,
(z0 ) dz dz 

which is (1.29).


Fluid Mechanics, 5th Ed. Kundu, Cohen, and
Dowling

Exercise 1.14. Starting with the hydrostatic pressure law (1.8), prove (1.30) without using
perfect gas relationships.

Solution 1.14. The adiabatic temperature gradient dTa/dz, can be written terms of the pressure
gradient:
dTa T  dp T 
    g 
dz p s dz p s
where the hydrostatic law dp/dz = –g has been used to reach the second equality. Here, the final
partial derivative can be exchanged for one involving  = 1/ and s, by considering:
h  h 
 dh    ds    dp  Tds  dp .
s p p s
  h    h  
Equality of the crossed second derivatives of h,          , implies:
p s p s s p s p
 T      T    s 
               ,
p s s p T p s p T p T p
where the second two equalities are mathematical
 manipulations that allow the introduction of
1     h  s 
        , and Cp     T  .

 T p T p T p T p
Thus,
dTa T    s  C  gT
 g   g      g  p   .
 dz  p 
s
T
p 
 T 
p  T  C p


Fluid Mechanics, 5th Ed. Kundu, Cohen, and
Dowling

Exercise 1.15. Assume that the temperature of the atmosphere varies with height z as T = T0 
g/KR
 T0 
Kz where K is a constant. Show that the pressure varies with height as p = p0   , where
T0  Kz 
g is the acceleration of gravity and R is the gas constant for the atmospheric gas.

Solution 1.15. Start with the hydrostatic and perfect gas laws, dp/dz = –g, and p = RT,

eliminate the density, and substitute in the given temperature profile to find:
dp p p dp g dz
  g   g g or  .
dz RT R(T0  Kz) p R (T0  Kz)
The final form may be integrated to find:
g
ln p   lnT0  Kz  const.
RK
At z = 0, 
the pressure must be p0, therefore: 
g
ln p0   lnT0   const.
RK
Subtracting this from theequation above and invoking the properties of logarithms produces:
 p  g T0  Kz 
ln   ln 
 p 0  RK  T0 
Exponentiating produces:
g/KR g/KR
p T0  Kz   T0 
=   , which is the same as: p = p0   .
p0   T0  T0  Kz 

 
Fluid Mechanics, 5th Ed. Kundu, Cohen, and
Dowling

Exercise 1.16. Suppose the atmospheric temperature varies according to: T = 15  0.001z, where
T is in degrees Celsius and height z is in meters. Is this atmosphere stable?

Solution 1.16. Compute the temperature gradient:

dT d C C
 (15  0.001z)  0.001  1.0 .
dz dz m km
For air in the earth's gravitational field, the adiabatic temperature gradient is:
dTa gT (9.81m /s2 )(1/T)T C
   9.8 .
 dz Cp 1004m /s C
2 2
km
Thus, the given temperature profile is stable because the magnitude of its gradient is less than
the magnitude of the adiabatic temperature gradient.


 Fluid Mechanics, 5th Ed. Kundu, Cohen, and
Dowling

Exercise 1.17. Consider the case of a pure gas planet where the hydrostatic law is:
dp dz  (z)Gm(z) z 2 . Here G is the gravitational constant, and m(z)  4  o ( ) 2 d is the
z

planetary mass up to distance z from the center of the planet. If the planetary gas is perfect with
gas constant R, determine (z) and p(z) if this atmosphere is isothermal at temperature T. Are
 these vertical profiles of  and p valid as z increases without
bound?
Solution 1.17. Start with the given relationship for m(z), differentiate it with respect to z, and use
the perfect gas law, p = RT to replace the  with p.
dm d  z  p(z)
 4   ( ) 2 d  4z 2 (z)  4 z 2 .
dz dz  0  RT
Now use this and the hydrostatic law to obtain a differential equation for m(z),
dp Gm(z) d  RT dm   1 dm Gm(z)
  (z)   2    .
dz z 2
dz 4z dz  4 z 2 dz  z 2

After recognizing T as a constant, the nonlinear second-order differential equation for m(z)
simplifies to:
RT d  1 dm  1 dm
  2   4 m .
G dz z dz  z dz
This equation can be solved by assuming a power law: m(z) = Azn. When substituted in, this trial
solution produces:
RT d 2 RT
 z Anz  n  3Anz n4  z4 A2nz 2n1.
n1

G dz G
Matching exponents of z across the last equality produces: n – 4 = 2n – 5, and this requires n = 1.
For this value of n, the remainder of the equation is:
RT RT
 2Az 3  z4 A 2z1 , which reduces to: A  2 .
G G
Thus, we have m(z) = 2RTz/G, and this leads to:
2RT 1 2R 2T 2 1
(z)  , and p(z)  .
 G 4 z 2
 G 4 z 2
Unfortunately, these profiles are not valid as z increases without bound, because this leads to an
unbounded planetary mass.

 
Solutions Manual for Fluid Mechanics 5th Edition by Kundu
Full Download: [Link]
Fluid Mechanics, 5th Ed. Kundu, Cohen, and
Dowling

Exercise 1.18. Consider a heat-insulated enclosure that is separated into two compartments of
volumes V1 and V2, containing perfect gases with pressures and temperatures of p1, p2, and T1, T2,
respectively. The compartments are separated by an impermeable membrane that conducts heat
(but not mass). Calculate the final steady-state temperature assuming each gas has constant
specific heats.

Solution 1.18. Since no work is done and no heat is transferred out of the enclosure, the final
energy Ef is the sum of the energies, E1 and E2, in the two compartments.
E1 + E2 = Ef implies 1V1Cv1T1 + 2V2Cv2T2 = (1V1Cv1 + 2V2Cv2)Tf,
where the Cv's are the specific heats at constant volume for the two gases. The perfect gas law
can be used to find the densities: 1 = p1/R1T1 and 2 = p2/R2T2, so
p1V1Cv1/R1 + p2V2Cv2/R2 = (p1V1Cv1/R1T1 + p2V2Cv2/R2T2)Tf.
A little more simplification is possible, Cv1/R1 = 1/(1 – 1) and Cv2/R1 = 1/(2 – 1). Thus, the final
temperature is:
p1V1 (1 1)  p2V2 ( 2 1)
Tf  .
p1V1 (1 1)T1   p2V2 ( 2 1)T2 



Full download all chapters instantly please go to Solutions Manual, Test Bank site: [Link]