Ideal Folow

CIVL4160 2013/1 Advanced fluid mechanics
IDEAL-FLUID FLOW TUTORIALS
TUTORIAL 1
Attendance to tutorials is very strongly advised. Repeated absences by some individuals will be noted and
these would demonstrate some disappointing responsible behaviour.
Past course results demonstrated a very strong correlation between the performances at the end-of-semester
examination, the attendance of tutorials during the semester and the overall course result.
1. Pre-Requisite Knowledge - Tutorials
The first tutorial consists of basic pre-requisite knowledge.
1.1 Give the following fluid and physical properties(at 20 Celsius and standard pressure) with a 4-digit
accuracy.
Value Units
Air density:
Water density:
Air dynamic viscosity:
Water dynamic viscosity:
Gravity constant in Brisbane:
Surface tension (air and water) :
1.2 What is the definition of an ideal fluid ?

What is the dynamic viscosity of an ideal fluid ?
1.3 From what fundamental equation does the Navier-Stokes equation derive : (a) continuity, (b) momentum
equation, (c) energy equation, (d) other ?
From what fundamental principle derives the Bernoulli equation ?
1.4 Sketch the streamlines of the following two-dimensional flow situations :

A- A laminar flow past a circular cylinder,
B- A turbulent flow past a circular cylinder,
In each case, show the possible extent of the wake (if any). Indicate clearly in which regions the ideal fluid
flow assumptions are valid, and in which areas they are not.
Remember the CIVL3130 Fluid Mechanics experiment "Flow past a cylinder".
CIVL4160-1
2. Ideal Fluid Flow - Irrotational Flows
2.1 Quizz
- What is the definition of the velocity potential ?
- Is the velocity potential a scalar or a vector ?
- Units of the velocity potential ?
- What is definition of the stream function ? Is it a scalar or a vector ? Units of the stream function ?
For an ideal fluid with irrotational flow motion :
- Write the condition of irrotationality as a function of the velocity potential.
- Does the velocity potential exist for 1- an irrotational flow and 2- for a real fluid ?
- Write the continuity equation as a function of the velocity potential.
Further, answer the following questions :
- What is a stagnation point ?
- For a two-dimensional flow, write the stream function conditions.
- How are the streamlines at the stagnation point ?
Reference
CHANSON, H. (2009). "Applied Hydrodynamics: An Introduction to Ideal and Real Fluid Flows." CRC
Press, Taylor & Francis Group, Leiden, The Netherlands, 478 pages.
VALLENTINE, H.R. (1969). "Applied Hydrodynamics." Butterworths, London, UK, SI edition.
2.2 Basic applications

(1) Considering the following velocity field :
Vx = y  z  t
Vy = z  x  t
Vz = x  y  t
- Is the flow a possible flow of an incompressible fluid ?

- Is the motion irrotational ? If yes : what is the velocity potential ?
(2) Considering the following velocity field :

Vx = 2  x
Vy = -2  y
Is the motion irrotational ? In the affirmative, what is the velocity potential ?
(3) Draw the streamline pattern of the following stream functions :

(3.1)  = 2 × x
(3.2)  = 3 × y
(3.3)  = 3 × x - 4 × y
CIVL4160-2
(3.4)  = -1.5 × x2
(3.5)  = 4 × Ln(x) - 2/y
Remember: a streamline is curve along which  is constant
2.3 Two-dimensional flow

Considering a two-dimensional flow, find the velocity potential and the stream function for a two-
dimensional flow having the following velocity components :
2xy
Vx = - 2
(x2 + y2)
x2 - y 2
Vy = 2
(x2 + y2)
2.4 Applications
(a) Using the software 2DFlowPlus, investigate the flow field of a vortex (at origin, strength 2) superposed to
a sink (at origin, strength 1). Visualise the streamlines, the contour of equal velocity ad the contour of
constant pressure.
Repeat the same process for a vortex (at origin, strength 2) superposed to a sink (at x=-5, y=0, strength 1).
How would you describe the flow region surrounding the vortex.
(b) Investigate the superposition of a source (at origin, strength 1) and an uniform velocity field (horizontal
direction, V = 1). How many stagnation point do you observe ? What is the pressure at the stagnation point ?
What is the "half-Rankine" body thickness at x = +1 ? (You may do the calculations directly or use
2DFlowPlus to solve the flow field.)
(c) Using 2DFlowPlus, investigate the flow past a circular building (for an ideal fluid with irrotational flow
motion). How many stagnation points is there ? Compare the resulting flow pattern with real-fluid flow
pattern behind a circular bluff body (search Reference text in the library).
(d) Investigate the seepage flow to a sink (well) located close to a lake. What flow pattern would you use ?
Note : the software 2DFlowPlus is described in the textbook, Appendix E. A demonstration copy can be
downloaded from the course website {www.uq.edu.au/~2ehchans/civ4160.html}. It is installed in the
undergraduate laboratory network.
Reference
Press, Taylor & Francis Group, Leiden, The Netherlands, 478 pages.
CIVL4160-3
2.5 Basic equations (2)
Considering an two-dimensional irrotational flow of ideal fluid, which basic principle(s) is(are) used to
determine the pressure field ?
2.6 Laplace equation

- What is the Laplacian of a function ? Write the Laplacian of the scalar function  in Cartesian and polar
coordinates.
- Rewrite the definition of the Laplacian of a scalar function as a function of vector operators (e.g. grad, div,
curl).
Solution
 2  2  2 
 (x,y,z) =    (x,y,z) = div grad (x,y,z) = + +
x2 y2 z2
Laplacian of scalar
    
 F (x,y,z) =    F (x,y,z) = i  Fx + j  Fy + k  Fz Laplacian of vector
1   1 2  2 
 (r,,z) = r  r   + 2  2 + Polar coordinates
r r  r  z2
It yields:

 f = div grad f
     
(
 F = grad div F - curl curl F )
where f is a scalar.
Note the following operations:
 (f + g) =  f +  g
 (F + G ) =  F +  G
 
 (f  g) = g   f + f   g + 2  grad f  grad g
where f and g are scalars.

For a two-dimensional ideal fluid flow, write:
(a) the continuity equation,
(b) the streamline equation,
(c) the velocity potential and stream function,
(d) the condition of irrotationality and
CIVL4160-4
(e) the Laplace equation
in polar coordinates.
Solution
Continuity equation
 Vx  Vy 1  (rVr) 1  V
+ = 0  + r   = 0
x y r r
Momentum equation
 Vx  Vx  Vx  P
+ Vx  + Vy  = -  + g  z
t x y x  
 Vy  Vy  Vy  P
+ Vx  + Vy  = - + g  z
t x y y  
 Vr V  Vr V2
 Vr  P
+ Vr  + r  - r = -  + g  z
t r  r  
 V  V V  V Vr  V 1  P
+ Vr  + r  
t r 
+ r = - r    + g  z
Streamline equation
Vx  dy - Vy  dx = Vr  r  d - V  dr = 0
Velocity potential and stream function

   1 
Vx = - = - V = - = -r*
x y r r 
  1  
Vy = - = + V = - r  = +
y x  r
Q = Q = 
Condition of irrotationality
 Vy  Vx  V 1  Vr
- = 0 - r = 0
x y r 
Laplace equation
2  2  2  1 2 
x2
+
y2
= 0
r2
+ r2  2 = 0
2  2  2  1 2 
+ = 0 + * = 0
x2 y2 r2 r2 2

Considering an two-dimensional irrotational flow of ideal fluid:
- write the Navier-Stokes equation (assuming gravity forces),
- substitute the irrotational flow condition and the velocity potential,
CIVL4160-5
- integrate each equation with respect to x and y,
- what is the final integrated form of the three equations of motion ?
This equation is called the Bernoulli equation for unsteady flow.
- For a steady flow write the Bernoulli equation. When the velocity is known, how do you determine the
pressure ?
Solution
 Vx  Vx  Vx  P
1- + Vx  + Vy  = -  + g  z
t x y x  
 Vy  Vy  Vy  P
+ Vx  + Vy  = - + g  z
t x y y  
2- Substituting the irrotational flow conditions:

 Vx  Vy
=
y x
and the velocity potential:

 
Vx = - Vy = -
x y
the equations may be expressed as (STREETER 1948, p. 24):

   Vx2 Vy2 P 
- + 2 + 2 + + g  z = 0
x  t  
   Vx2 Vy2 P 
- + 2 + 2 + + g  z = 0
y  t  
3- Integrating with respect to x and y:

 V2 P
- + 2 + + g  z = Fx(y, t)
t 
 V2 P
- + 2 + + g  z = Fy(x, t)
t 
where V is defined as the magnitude of the velocity: V2 = Vx2 + Vx2. The integration of the three motion
equations are identical and the left-hand sides of the equations are the same:
Fx(y, t) = Fy(x, t)
4- The final integrated form of the three equations of motion is the Bernoulli equation for unsteady flow:
 V2 P
- + 2 + + g  z = F(t)
t 
 
In the general case of a volume force potential U (i.e. Fv = - grad U ):
 V2 P
- + 2 + + U = F(t)
t 
CIVL4160-6
More exercises in textbook pp. 24-25, 45-53.
Press, Taylor & Francis Group, Leiden, The Netherlands, 478 pages (ISBN: 978-0-415-49271-3
(Hardback); 978-0-203-87626-8 (eBook)) (PSE library: TC171 .C54 2009).
CIVL4160-7
TUTORIAL 2
3. Two-Dimensional Flows (1) Basic equations and flow analogies - Tutorials
3.1 Flow net (2006 examination paper)

A hydrodynamic study of the flow past a flat plate is conducted in water. The plate length is L = 0.2 m and
the plate width equals the test section width (0.5 m). For a particular angle of incidence, dye injection in a
water tunnel gives the flow pattern shown in Figure E3-1 on the next page. The mean flow is horizontal.
(a) Complete the flow net by drawing the suitable equipotentials. (Draw the complete flow net on the
examination paper.)
(b) The upstream velocity is 6.5 m/s. Estimate the drag force and the lift force acting on the plate. Indicate
clearly the sign convention. You may have to draw more streamline and equipotentials to obtain a good
accuracy.
Assume water at 20 Celsius ( = 998.2 kg/m3,  = 1.005 E-5 Pa.s,  = 0.0736 N/m).
CIVL4160-8
Fig. E3-1 - Flow visualisation and streamline patterns around an inclined flat plate (undistorted scale)
CIVL4160-9
3.2 Flow net
Let us consider the water flow past a circular arc sketched in Figure E3-2. The chord is L = 0.35 m and the
plate width equals 1 m. For a particular angle of incidence shown in Figure E3-2, conduct a graphical
analysis of the flow field. The mean flow is horizontal.
(a) Draw the flow net by drawing the suitable streamlines and equipotentials. (Draw the complete flow net
with sufficient details in the vicinity of the circular arc.)
(b) The upstream velocity is 15.2 m/s. Estimate the drag force and the lift force acting on the cambered plate.
Indicate clearly the sign convention..
Assume water at 20 Celsius ( = 998.2 kg/m3,  = 1.005 E-5 Pa.s,  = 0.0736 N/m).
CIVL4160-10
Fig. E3-2 - Flow past a circular arc (undistorted scale)
CIVL4160-11
3.3 Flow net beneath a cutoff wall
For a two-dimensional seepage under a impervious structure with a cutoff wall (Fig. E3-3), the boundary
conditions are : H = 6 m, K = 2.0 m/day.
(a) What is the hydraulic conductivity in m/s ? What type of soil is it ?
(b) Using the flow net, estimate the seepage flow par meter width of dam.
Fig. E3-3 - Flow net beneath an impervious dam
3.4 Flow net under a sheetpiling

For a two-dimensional seepage under sheetpiling with a permeable foundation (Fig. E3-4), the boundary
conditions are : a = 9.4 m, b = 4.7 m, c = 4 m, H = 2.5 m, K = 2.0 10-3 cm/s.
Using a dimensioned flow net (with 5 stream tubes or more), estimate the seepage flow par meter width of
dam.
Use graph paper.
CIVL4160-12
Fig. E3-4 - Flow net beneath a sheet pile
3.5 Flow net under an impervious dam

For a two-dimensional seepage under an impervious dam with apron (textbook, CHANSON 2009, Fig. 3-
1B), the boundary conditions are : H = 80 m, K = 5 10-5 m/s.
(A) Calculate the seepage flow rate in presence of an apron and a cutoff wall (Fig. 3-1B).
(B) In absence of the cutoff wall, sketch the flow net and determine the seepage flow. Determine the pressure
distribution along the base of the dam and beneath the apron. Calculate the uplift forces on the dam
foundation and on the apron.
3.6 Flow net under an impervious dam with cutoff wall

For a two-dimensional seepage under an impervious dam with a cutoff wall (Fig. E3-5), the boundary
conditions are : H1 = 60 m, H2 = 5 m, a = 60 m, b = 100 m, L = 70 m, K = 1 10-5 m/s.
(A) In absence of cutoff wall (i.e. a = 0 m), sketch the flow net; determine the seepage flow; determine the
pressure distribution along the base of the dam; calculate the uplift force.
(B) With the cutoff wall (i.e. a = 60 m): same questions : sketch the flow net; determine the seepage flow;
determine the pressure distribution along the base of the dam; calculate the uplift force
(C) Comparison and discuss the results.
Use graph paper.
CIVL4160-13
Fig. E3-5 - Flow net beneath a dam with cutoff wall
3.7 Pressure distribution above some bed forms

Figure E3-6 illustrates an open channel flow above some quasi-sinusoidal bed forms. The wave length is 2.2
m and the average water depth is 3.3 m. Neglecting sediment motion, calculate the pressure distributions at
the bed form crest and trough.
On a flat bed, the pressure distribution is hydrostatic. Determine the deviations from the hydrostatic
pressure distribution.
CIVL4160-14
Fig. E3-6 - Streamlines past a wavy bed
3.8 Cauchy-Riemann equations

- Determine the stream function for parallel flow with a velocity V inclined at an angle  to the x-axis.
- Find the velocity potential using the Cauchy-Riemann equations.
3.9 Flow analogy

A Hele-Shaw cell apparatus is built with a gap between plates  = 1.5 mm filled with Glycerol ( = 1260
kg/m3,  = 1.4 Pa.s). Estimate the hydraulic conductivity of the apparatus.
CIVL4160-15
More exercises in textbook pp. 74-80.
Exercise Solutions
Exercise 3.1
The problem is solved by drawing the equipotentials and completing the flow net. Close to the foil,
additional streamlines and equipotentials may be drawn to improve the estimates of the velocities next to the
extrados and intrados of the foil.
The pressure field is derived from the Bernoulli equation and the integration of the pressure distributions
yields the lift and drag forces.
Remarks
- The flow net method and technique are presented in the textbook (CHANSON 2009).
- The complete theory of lift and drag on airfoils, wings and hydrofoils is developed in the textbook
(CHANSON 2009, chapter I-6).
Exercise 3.2
The problem may be solved analytically using the Joukowski transformation and theorem of Kutta-
Joukowski. The complete theory of lift and drag on airfoils, wings and hydrofoils is developed in the
textbook (CHANSON 2009, chapter I-6).
The flow net solution can be compared with the theory of lift and drag.
Exercise 3.3
Solution
q = 4.6 m2/day
Remark
See Textbook (CHANSON 2009), pp. 75-76.
Exercise 3.4
Solution
q = 2.7 m2/day
Discussion
CIVL4160-16
The flow pattern may be analysed analytically using a finite line source (for the sheet pile) and the theory of
images. The resulting streamlines are the equipotentials of the sheet-pile flow.
Remember : A velocity potential can be found for each stream function. If the stream function satisfies the
Laplace equation the velocity potential also satisfies it. Hence the velocity potential may be considered as
stream function for another flow case. The velocity potential  and the stream function  are called
"conjugate functions" (Chapter I-2).
Exercise 3.5
Solution
(A) The problem is similar to the flow net sketched in Figure 3-1B (Textbook, pp. 58-60).
(B) In absence of cutoff wall, the streamlines are shorter and the seepage flow rate is greater.
The pressure distribution beneath the dam foundation and apron may be deduced from the equipotential
lines, since  = KH where K is the hydraulic conductivity and H is the piezometric head (Chapter I-3,
paragraph 3.3).
The uplift pressure on the dam foundation and apron are very significant. The apron structure would be
subjected to high risks of uplift and damage.
Remark
See Textbook (CHANSON 2009), pp. 60.
Exercise 3.6
See Textbook (CHANSON 2009), pp. 77-78.
Exercise 3.7
Solve graphically the problem : (a) Complete the equipotential lines; (b) Calculate the velocity magnitude at
all vertical elevation; (c) Apply the Bernoulli principle.
Remember that the pressure gradient is hydrostatic far away upstream.
Application
This flow pattern is typical of the flow above standing wave bed forms, although these tend to be
significantly larger (KENNEDY 1963, CHANSON 2000). The height of the wall undulations is closer to
large ripples and small dunes (CHANSON 2204, pp. 151-155 & 223-231).
KENNEDY, J.F. (1963). "The Mechanics of Dunes and Antidunes in Erodible-Bed Channels." Jl of Fluid
Mech., Vol. 16, No. 4, pp. 521-544 (& 2 plates).
CHANSON, H. (2000). "Boundary Shear Stress Measurements in Undular Flows : Application to Standing
Wave Bed Forms." Water Res. Res., Vol. 36, No. 10, pp. 3063-3076 (ISSN 0043-1397).
{http://espace.library.uq.edu.au/view.php?pid=UQ:11117}
CIVL4160-17
CHANSON, H. (2004). "The Hydraulics of Open Channel Flow : An Introduction." Butterworth-
Heinemann, Oxford, UK, 2nd edition, 630 pages (ISBN 978 0 7506 5978 9).
Exercise 3.8
Solution
By definition:

= - Vx = - V  cos = - V  cos  y + f1(x)
y

= Vy = V  sin = V  sin  x + f2(y)
x
Hence:
 = V  sin  x - V  cos  y + constant
A streamline equation is:


y = tan  x -
V  cos
which is a straight line equation.

The Cauchy-Riemann equations are:
 
=  = - V  cos  x + g1(y)
x y
 
= -  = - V  sin  y + g2(x)
y x
Hence:
 = - V  sin  y - V  cos  x + constant
Vx = Vo  cos
Vy = Vo  sin
 = - Vo  (y  cos - x  sin)
Exercise 3.9
Solution
In laminar flows between two parallel plates, the application of the momentum principle in its integral form
yields an expression of the longitudinal head loss:
H f V2
= -D 
x H 2g
where DH is the equivalent pipe diameter (or hydraulic diameter), V is the cross-sectional averaged velocity
and f is the Darcy-Weisbach friction factor:
64  
f =
  V  DH
Since DH = 2×, it yields:
CIVL4160-18
   2 H
V 
8   x
By analogy with the Darcy law, this gives a hydraulic conductivity K:

  2
K = 2.5×10-4 m/s
8
CIVL4160-19
TUTORIAL 3
4. Two-Dimensional Flows (2) Basic flow patterns - Tutorials
4.1 Doublet in uniform flow (1)

Select the strength of doublet needed to portray an uniform flow of ideal fluid with a 20 m/s velocity around
a cylinder of radius 2 m.
4.2 Source and sink

A source discharging 0.72 m2/s is located at (-1, 0) and a sink of twice the strength is located at (+2, 0). For a
remote pressure (far away) of 7.2 kPa,  = 1,240 kg/m3, find the velocity and pressure at (0, 1) and (1, 1).
Note : When some measurements are conducted with a Prandtl-Pitot tube, the pressure tapping at the leading
edge of the tube gives the dynamic pressure, while the pressure tappings on the side give the piezometric
pressure. Remember that, at the leading edge of the tube, stagnation occurs.
Remarks
The Pitot tube is named after the Frenchman Henri PITOT. The first presentation of the concept of the Pitot tube
was made in 1732 at the French Academy of Sciences by Henri PITOT. The original Pitot tube included basically a
total head reading. Ludwig PRANDTL improved the device by introducing a pressure (or piezometric head)
reading. The modified Pitot tube is sometimes called a Pitot-Prandtl tube.
For many years, aeroplanes used Prandtl-Pitot tubes to estimate their relative velocity.
4.3 Flow pattern (2)

In two-dimensional flow we now consider a source, a sink and an uniform stream. For the pattern resulting
from the combinations of a source (located at (-L, 0)) and sink (located at (+L, 0)) of equal strength Q in
uniform flow (velocity +Vo parallel to the x-axis) :
(a) Sketch streamlines and equipotential lines;
(b) Give the velocity potential and the stream function.
CIVL4160-20
This flow pattern is called the flow past a Rankine body. W.J.M. RANKINE (1820-1872) was a Scottish
engineer and physicist who developed the theory of sources and sinks. The shape of the body may be altered
by varying the distance between source and sink (i.e. 2L) or by varying the strength of the source and sink.
Other shapes may be obtained by the introduction of additional sources and sinks and RANKINE developed
ship contours in this way.
(c) What is the profile of the Rankine body (i.e. find the streamline that defines the shape of the body)?
(d) What is the length and height of the body ?
(e) Explain how the flow past a cylinder can be regarded as a Rankine body. Give the radius of the cylinder
as a function of the Rankine body parameter.
4.4 Flow pattern (3)

In two-dimensional flow we consider again a source, a sink and an uniform stream. But. the source is located
at (+L, 0) and the sink is located at (-L, 0) (i.e. opposite to a Rankine body flow pattern). They are of equal
strength q in an uniform flow (velocity +Vo parallel to the x-axis).
Derive the relationship between the discharge q, the length L and the flow velocity such that no flow injected
at the source becomes trapped into the sink.
4.5 Doublet in uniform flow (2)

We consider the air flow (Vo = 9 m/s, standard conditions) past a suspension bridge cable (Ø = 20 mm),
(a) Select the strength of doublet needed to portray the uniform flow of ideal fluid around the cylindrical
cable.
(b) In real fluid flow, calculate the hydrodynamic frequency of the vortex shedding.
4.6 Flow past buildings - 2003 exam paper

Let us consider a new architectural landmark to be built at the Mt Cootha Lookout. The structure consists of
three circular cylinders (Height : 25 m - Diameter : 2, 3 and 5 m). The landmark will be facing North-East,
while the dominant winds are Easterlies (Fig. 1).
You will assume that the wind flow around the structure is a two-dimensional irrotational flow of ideal fluid.
The atmospheric conditions are : P = Patm = 105 Pa ; T = 25 Celsius.
CIVL4160-21
Fig. 1 - Sketch of the Mt Cootha landmark - View in elevation
(1) On graph paper, sketch the flow net with the landmark for a 25 m/s Easterly wind. Indicate clearly on the
graph the discharge between two streamlines, the x-axis and y-axis, their direction, and use the centre of the
5-m diameter cylinder as the origin of your system of coordinates (with x in the South-East direction and y in
the North-East direction).
(2) The Brisbane City Council is concerned about wind velocities between the buildings that may blow down
tourists and damage cars.
(a) From your flow net, compute the wind velocity and the pressure at :
x = 1.5 m, y = 6 m
x = 5 m, y =4.5 m
x = 2.1 m, y = 2.1 m
These locations would be typical of tourists standing in front of the vertical cylinders.
CIVL4160-22
(b) Where is located the point of maximum velocity and minimum pressure ?
(c) What is the maximum velocity and minimum pressure between the cylinders ? Indicate that location on
your flow net.
(d) Discuss your results. Do you think that this result is realistic ? Why ?
(3) Explain what standard flow patterns you would use to describe the flow around these three buildings.
(4) Write the stream function and the velocity potential as a function of the wind speed Vo (25 m/s) and the
two cylinder diameters D1 (5 m), D2 (3 m) and D3 (2 m).
Do not use numbers. Express the results as functions of the above symbols.
(5) For a real fluid flow, what is (are) the drag force(s) on each cylinder (Height: 25 m) ?
4.7 Magnus effect (1)

Two 15-m high rotors 3 m in diameter are used to propel a ship. Estimate the total longitudinal force exerted
upon the rotors when the relative wind velocity is 25 knots, the angular velocity of the rotors is 220
revolutions per minute and the wind direction is at 60º from the bow of the ship.
Perform the calculations for (a) an ideal fluid with irrotational motion and (b) a real fluid.
(c) What orientation of the vector of relative wind velocity would yield the greatest propulsion force upon
the rotorship ? Calculate the magnitude of this force.
Assume a real fluid flow. The result is trivial for ideal fluid with irrotational flow motion.
(d) Determine how nearly into the wind the rotorship could sail. That is, at wind angle would the resultant
propulsion force be zero ignoring the wind effect on the ship itself.
4.8 Magnus effect (2)

A infinite rotating cylinder (R = 1.1 m) is placed in a free-stream flow (Vo = 15 knots) of water. The cylinder
is rotating at 45 rpm and the ambient pressure far away is 1.1 E+5 Pa.
(a) Calculate and plot the pressure distribution on the cylinder surface.
(b) Find the maximum and minimum pressures on the cylinder surface.
(c) Find the location  where the pressure on the cylinder surface is the fluid pressure far away from the
cylinder.
4.9 Whirlpools
Whirlpools may be approximated by a series of vortices of same signs advected into an uniform flow.
(a) Consider two vortices of equal strength K = +1 located at (-2, ) and (+2, 0). Estimate how far away the
effect of the vortices is perceived to be that of an unique vortex. What would be the strength of that vortex ?
CIVL4160-23
(b) Consider two vortices of equal strength K = +1 located at (-2, ) and (+2, 0) in a horizontal uniform flow
V = +0.03. What are the stream function and velocity potential of the resulting flow motion.
4.10 Magnus effect aircraft

Two rotating cylinders are used instead of conventional wings to provide the lift to an aircraft. Calculate the
length of each cylinder wing for the following design conditions :
Cruise speed : 320 km/h
Cruise altitude : 2,000 m
Aircraft mass : 8 E+6 kg
Cylinder radius : 1.8 m
Cylinder rotation speed : 500 rpm
4.11 Wind force on a Nissen hut

A 45 km/h wind flows over a Nissen hut which has a 3.5 m radius and is 54.9 m long (Fig. E4-1). The
upstream pressure and temperature are 1.013 E+5 Pa and 288.2 K respectively, and they are equal to the
pressure and temperature inside the Nissen hut.
(a)`Calculate the lift and drag forces on the building.
(b) Find the location  on the building roof where the pressure is 1.013 E+5 Pa.
Assume an irrotational flow motion of ideal fluid in parts (a) and (b).
(c) Calculate the drag force for a real fluid flow.
Notes: The Nissen hut is a building made from a semi-circle of corrugated steel. A variant was the Quonset
hut used extensively during World War 2 by the Commonwealth and US military for army camps and air
bases. The design was named after Major Peter Norman NISSEN, 29th Company Royal Engineers who
experimented with hut design and constructed three prototype semi-circular huts in April 1916. The building
was not only economical but also portable.
CIVL4160-24
Fig. E4-1 - Wind flow past a Nissen hut
4.12 Wind flow past columns (2008 Exam paper)

The facade of a 150 m tall building is hindered by two cylindrical columns (a service column and a structural
column) (Fig. E4-2). There are concerns about the wind flow around the columns during storm conditions
(Vo = 35 m/s).
The column dimensions are:
Column 1 Column 2 Units
D 0.8 2.0 m
x (centre) 0 2.8 m
y (centre) 1.0 2.0 m
You will assume that the wind flow around the building is a two-dimensional irrotational flow of ideal fluid.
The atmospheric conditions are: P = Patm = 105 Pa; T = 20 Celsius.
(a) On graph paper, sketch the flow net for a 35 m/s wind flow. Indicate clearly on the graph the discharge
between two streamlines, the x-axis and y-axis, and their direction.
(b) The developer is concerned about the maximum wind speeds next to the building facade that may
damage the cladding and service column (column 1).
(b1) From your flow net, compute the wind velocity and the pressure at:
x = 0 m, y = 0.4 m
x = 2.8 m, y = 1.0 m
x = 1.4 m, y = 0.8 m
(b2) Where is located the point of maximum wind velocity and minimum ambient pressure ? Show
that location on your flow net.
(b3) Discuss your results. Do you think that these results are realistic ? Why ?
(c) Explain what standard flow patterns you would use to describe the flow around these two columns.
(d) Write the stream function and the velocity potential as a function of the wind speed Vo (35 m/s), the two
cylinder diameters D1 and D2, and their locations (x1, y1) and (x2, y2). Do not use numbers. Express the
results as functions of the above symbols.
(e) For a real fluid flow, what is (are) the drag force(s) on each cylinder (Height: 150 m) neglecting
interactions between the columns?
CIVL4160-25
Fig. E4-2 - View in elevation of the building facade and cylindrical columns
More exercises in textbook pp. 74-80.
Exercise Solutions
Exercise 4.1
Solution
(a) A doublet and uniform flow is analogous to the flow past a cylinder of radius :
-
R = Vo
where  is the strength of the doublet. Hence :

 = - Vo  R2 = 80 m3/s
Remark
See Textbook (CHANSON 2009), Chapter I-4.
Exercise 4.3
Solution
The flow past a Rankine body is the pattern resulting from the combinations of a source and sink of equal
strength in uniform flow (velocity +Vo parallel to the x-axis) :
 = - Vo  r  cos -  +  Lnr1
q r
 2  2 
 = - Vo  r  sin -  +  (1 - 2)

q
 2 
CIVL4160-26
where the subscript 1 refers to the source, the subscript 2 to the sink and q is positive for the source located at
(-L, 0) and the sink located at (+L, 0).
The profile of the Rankine body is the streamline  = 0 :
q
 = - Vo  r  sin +  (1 - 2) = 0
2
q  (1 - 2)
r =
2    Vo  sin
The length of the body equals the distance between the stagnation points where :

q q q 1 1 
V = Vo + - = Vo + - = 0
2    r1 2    r2 2   rs - L rs + L
and hence :
q
Lbody = 2  rs = 2  L  1 +
  L  Vo
The half-width of the body h is deduced from the profile equation at the point (h, /2) :
q  (1 - 2)
h =
2    Vo
where : 1 =  and 2 =  -  and hence :

   h  Vo
 = 2 - q
But also :
h
tan = L
So the half-width of the body is the solution of the equation :

V
h = L  cot q o  h
 
Remark
Exercise 4.4
Exercise 4.5
Solution
(a) A doublet and uniform flow is analogous to the flow past a cylinder of radius :
-
R = Vo
where  is the strength of the doublet. Hence :

 = - Vo  R2 = - 9 E-4 m3/s
CIVL4160-27
(b) The Reynolds number of the flow is 1.1 E-4. For that range of Reynolds number, the vortex shedding
behind the cable is characterised by a well-defined von Karman street of vortex. The hydrodynamic
frequency satisfies :
2R
St = Vo ~ 0.2
It yields :  = 90 Hz. If the hydrodynamic frequency happens to coincide with the natural frequency of the
structure, the effects may be devastating : e.g., Tacoma Narrows bridge failure on 7 November 1940.
4.6 2003 exam paper

(2)
(a1) x = 1.5 m, y = 6 m V ~ 23.2 m/s P - Patm = +52 Pa
(a2) x = 5 m, y =4.5 m V ~ 28.5 m/s P - Patm = -112 Pa
(a3) x = 2.1 m, y = 2.1 m V ~ 10.5 m/s P - Patm = +309 Pa
(b) between the 5-m and 2-m diameter cylinders
(c) V ~ 60 m/s, P - Patm = -1785 Pa

Such a large maximum wind velocity may cause a potential hazard for pedestrians and cyclists.
(d) The flow is turbulent : VD/ ~ 3.2 E+6 (D = 2 m). Separation is likely to occur behind the cylinder.
However, since the location of maximum velocity is likely to be outside of a wake region, the above results
are very likely representative.
(3) Use 3 doublet patterns ( = 100, 225, 625 m3/s)
(5)
D (m)= 2 3 5
Re = 3.2 E+6 4.8 E+6 8 E+6
CD = 0.7 0.75 0.75
Drag (N) = 1.9 E+4 2.1 E+4 3.5 E+4
Exercise 4.7
Solution
The tangential velocity of the rotors is :
R   = 2  220 / 60  1.5 = 34.56 m/s
The relative velocity of the wind is :
Vo = 25  1852 / 3600 = 12.86 m/s
Note : if the wind comes from starboard, the rotation of the rotor masts must be in the trigonometric positive
direction to propel the ship forward.
CIVL4160-28
(a) For an ideal fluid with irrotational motion :
Total Lift = 2  15  1.2  12.86  13.57 = 150.8 kN
In the direction of flow motion, the total force is :
Total force = Total Lift  cos30º = 131 kN
(b) The total lift and drag forces are :
Total Lift ~ 62.5 kN
Total Drag ~ 17.9 kN
In the direction of flow motion, the total force is :
Total force = 45.2 kN
Exercise 4.8
Solution
(a) The flow pattern has two stagnation points.
(b) The minimum and maximum pressures at the cylinder surface are respectively -7.24 E+4 and +1.40 E+5
Pa for  = 270º and 20º(and 160) respectively.
Notes : (1) The minimum pressure is sub-atmospheric and may lead to some cavitation. (2) There are two
locations the pressure is maximum which correspond both the location of a stagnation point.
Exercise 4.9
Solution
Use 2D Flow Plus to asses the flow pattern.
A whirlpool is a vortex of vertical axis, with a downward velocity component near its centre. A good
example is the bathtub vortex. VAN DYKE (1982, p. 59) presented a superb illustration. See also the
Queensland Science Museum. A related example is the vortex dropshaft design. In coastal zones, whirlpools
are produced by the interaction of rising and falling tides. They are often observed at the edges of straits with
large tidal currents. (At Naruto, currents of up to 9 knots were observed.) The vortex (whirlpool) is a
coherent structure typical of shear flows where there is a velocity difference across the shear layer. It affect
the surrounding flow and water can be seen going back and forth across the shear layer between vortices.
Notable oceanic whirlpools include those of Garofalo along the coast of Calabria in southern Italy, and of
Messina in the strait between Sicily and peninsular Italy, the Maelstrøm (from Dutch for "whirling stream")
located near the Lofoten Islands off the coast of Norway. Whirlpools near the Hebrides and Orkney islands,
and in the Naruto strait between Awaji and Shikoku islands, are also well known.
See Textbook (CHANSON 2009), App. F.

Also : {http://www.uq.edu.au/~e2hchans/whirlpl.html}.
CIVL4160-29
Exercise 4.10
Solution
The air density at 2,000 m altitude is about 1.1 kg/m3.
The ideal fluid flow pattern is the superposition of an uniform flow (Vo = 88.9 m/s), a doublet (strength )
and a vortex strength K.
The cylinder radius and doublet strength are linked as :

R
Vo
while the vortex strength and rotation speed ( = 52.36 rad/s) satisfy :
K K  Vo
 = =
2    R2 2
The lift force per unit length of wing is

Lift = -   Vo  K
The right wing must rotate in the anti-clockwise direction as seen by the pilot and  must be negative to
provide a positive lift.
The calculations for cruise conditions imply that each wing is 376 m long (!).
Note : The Magnus effect lift force decreases with decreasing speed and it becomes small at take-off and
landing conditions.
Exercise 4.11
Solution
(a) Using the method of images, the flow field is that a doublet in an uniform flow. The pressure on the roof
is
1
P = Po + 2    Vo2  (1 - 4  sin2)
where Po = 1.013 E+5 Pa, Vo = 12.5 m/s and  = 1.225 kg/m3.

The flow direction is positive from left to right, and the doublet is centred at the origin. The lift and drag
forces on the building are respectively:

 (P - Po)  R  sin  d
Lift = - 
0
 (P - Po)  R  cos  d
Drag = - 
0
The integration yields:

Lift = + 2    Vo2  R = 1.3 kN/m
Drag = 0
The total lift force on the 54.9 m long Nissen hut is 72 kN.
CIVL4160-30
Notes
(1 - 4sin2)  d = sin(2) - 


 2 1
(1 - 4sin )  sin  d = 2  cos() - 3  cos(3  )
 2 1
(1 - 4sin )  cos  d = 3  sin(3  )
(b) The pressure on the building roof is 1.013 E+5 Pa for  = +30º (/6) and +150º (5/6).
(c) For a real fluid flow, separation takes place on the downstream of the building. The Reynolds number 
VoD/ equals 5.6 E+6, and the corresponding drag coefficient is about 0.8 (Fig. 4-13). The total drag force
equals hence:
1 1
Drag force = 2  CD  2    Vo2  D  L = 14 kN
where D = 2R.
CIVL4160-31

Ideal Folow

Uploaded by

Copyright:

Available Formats

Ideal Folow

Uploaded by

Document Information

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Ideal Folow

Uploaded by

Copyright:

Available Formats

CIVL4160 2013/1 Advanced fluid mechanics

IDEAL-FLUID FLOW TUTORIALS

1. Pre-Requisite Knowledge - Tutorials

The first tutorial consists of basic pre-requisite knowledge.

1.2 What is the definition of an ideal fluid ?

1.4 Sketch the streamlines of the following two-dimensional flow situations :

2.2 Basic applications

- Is the flow a possible flow of an incompressible fluid ?

(2) Considering the following velocity field :

Is the motion irrotational ? In the affirmative, what is the velocity potential ?

(3) Draw the streamline pattern of the following stream functions :

2.3 Two-dimensional flow

2.6 Laplace equation

where f and g are scalars.

2.7 Basic equations (2)

Velocity potential and stream function

2.8 Basic equations (3)

2- Substituting the irrotational flow conditions:

and the velocity potential:

the equations may be expressed as (STREETER 1948, p. 24):

3- Integrating with respect to x and y:

3. Two-Dimensional Flows (1) Basic equations and flow analogies - Tutorials

3.1 Flow net (2006 examination paper)

Fig. E3-3 - Flow net beneath an impervious dam

3.4 Flow net under a sheetpiling

3.5 Flow net under an impervious dam

3.6 Flow net under an impervious dam with cutoff wall

3.7 Pressure distribution above some bed forms

3.8 Cauchy-Riemann equations

3.9 Flow analogy

A streamline equation is:

which is a straight line equation.

Since DH = 2×, it yields:

By analogy with the Darcy law, this gives a hydraulic conductivity K:

4. Two-Dimensional Flows (2) Basic flow patterns - Tutorials

4.1 Doublet in uniform flow (1)

4.2 Source and sink

4.3 Flow pattern (2)

4.4 Flow pattern (3)

4.5 Doublet in uniform flow (2)

4.6 Flow past buildings - 2003 exam paper

4.7 Magnus effect (1)

4.8 Magnus effect (2)

4.10 Magnus effect aircraft

4.11 Wind force on a Nissen hut

4.12 Wind flow past columns (2008 Exam paper)

More exercises in textbook pp. 74-80.

where  is the strength of the doublet. Hence :

 = - Vo  r  sin -  +  (1 - 2)

where : 1 =  and 2 =  -  and hence :

So the half-width of the body is the solution of the equation :

where  is the strength of the doublet. Hence :

4.6 2003 exam paper

(b) between the 5-m and 2-m diameter cylinders

(c) V ~ 60 m/s, P - Patm = -1785 Pa

(3) Use 3 doublet patterns ( = 100, 225, 625 m3/s)

See Textbook (CHANSON 2009), App. F.

The lift force per unit length of wing is

where Po = 1.013 E+5 Pa, Vo = 12.5 m/s and  = 1.225 kg/m3.

The integration yields: