Wind Tunnels
Wind Tunnels
Wind Tunnels
1. Introduction
The subject of fluid mechanics is filled with abstract concepts, mathematical methods, and
results. Historically, it has been a challenging subject for students, undergraduate and
graduate. In most institutions, the introductory course in fluid mechanics is accompanied
by a laboratory course. While institutional philosophy and orientation vary around the
world, the goal of that laboratory is to strengthen students’ understanding of fluid
mechanics using a variety of laboratory exercises (Feisel & Rosa, 2005).
The literature has identified six basic functions of experimental work. Indeed, the report
of the Laboratory Development Committee of the Commission on Engineering
Education identified six key functions and objectives of the instructional laboratory
(Ernest, 1983):
a. Familiarization
b. Model identification
c. Validation of assumptions
d. Prediction of the performance of complex systems
e. Testing for compliance with specifications
f. And exploration for new fundamental information.
The report states that “The role of the undergraduate instructional laboratory is to teach
student engineers to perform these six functions. Hence the primary goal of undergraduate
laboratories is to inculcate into the student the theory and practice of experimentation. This
includes instrumentation and measurement theory.” (Ernst, 1983).
The wind tunnel is one such instrument. This chapter focuses on the measurement theory on
which the wind tunnel is based and presents examples of its use in the undergraduate fluid
mechanics laboratory at Indiana University-Purdue University Fort Wayne, Fort Wayne,
Indiana, USA.
The remainder of the chapter is organized in the following manner:
1. Basic concepts discuss definitions, classifications, and various uses of wind tunnels.
2. Fundamental Equations present the equations that are used as foundations for the theory
and application of wind tunnels.
3. Applications of wind tunnels in teaching fluid mechanics present nine different examples
that are used in our laboratory to teach various aspects of fluid mechanics and its uses
in design, testing, model verification, and research.
4. References list all cited works in alphabetical order.
www.intechopen.com
236 Wind Tunnels and Experimental Fluid Dynamics Research
2. Basic concepts
2.1 Definition of a wind tunnel
A wind tunnel is a specially designed and protected space into which air is drawn, or
blown, by mechanical means in order to achieve a specified speed and predetermined flow
pattern at a given instant. The flow so achieved can be observed from outside the wind
tunnel through transparent windows that enclose the test section and flow characteristics
are measurable using specialized instruments. An object, such as a model, or some full-scale
engineering structure, typically a vehicle, or part of it, can be immersed into the established
flow, thereby disturbing it. The objectives of the immersion include being able to simulate,
visualize, observe, and/or measure how the flow around the immersed object affects the
immersed object.
2.2.1 Type 1 classification – The criterion for classification is the path followed by the
drawn air: Open- vs. closed-circuit wind tunnels
Open-circuit (open-return) wind tunnel. If the air is drawn directly from the surroundings
into the wind tunnel and rejected back into the surroundings, the wind tunnel is said to
have an open-air circuit. A diagram of such a wind tunnel is shown in Figure 1.
Fig. 1. Diagram of an open-circuit, also known as open-return, wind tunnel (from NASA)
An open-circuit wind tunnel is also called an open-return wind tunnel.
Closed-circuit, or closed-return, wind tunnel. If the same air is being circulated in such a
way that the wind tunnel does neither draw new air from the surrounding, nor return it into
www.intechopen.com
Wind Tunnels in Engineering Education 237
the surroundings, the wind tunnel is said to have a closed-air circuit. It is conventional to
call that a closed-circuit (closed-return ) wind tunnel. Figure 2 illustrates this configuration.
Fig. 2. Top view of a closed-circuit, also known as closed-return, wind tunnel ( NASA)
www.intechopen.com
238 Wind Tunnels and Experimental Fluid Dynamics Research
www.intechopen.com
Wind Tunnels in Engineering Education 239
forces and moments on airplane wings, airfoils, and tall buildings. A close-up view of a
model of an F-5 fighter plane mounted in the test section of a wind tunnel is shown in
Figure 4.
Fig. 4. Close-up of a tufted model of an F-5 fighter plane in the test section of a wind tunnel
(NASA)
www.intechopen.com
240 Wind Tunnels and Experimental Fluid Dynamics Research
2 ( p0 - p)
V= (1)
r
Where V is the speed of the fluid, P0 is the total, also called the stagnation, pressure at that
point of measurement, and p is the static pressure at the same point. This equation comes
from the application of Bernoulli’s equation for the steady flow of an incompressible and
inviscid fluid along a streamline. Bernoulli’s equation is typically obtained by integrating
Euler’s equations along a streamline. It will be recalled that Euler’s equations are a special
case of the Navier -Stokes equations, when the viscosity of the fluid has been neglected. The
Navier-Stokes’ equations, in turn, are obtained from Newton’s second law when it is
applied to a fluid for which the shear deformation follows Newton’s law of viscosity.
Accordingly, in order to establish the theoretical validity of this equation for use in
educational wind tunnels, it is important to review some basic results from the theory of
viscous and inviscid flows. For the interested reader, these are available in all introductory
textbooks of fluid mechanics (e,g. Pritchard, 2011). For this reason, the rest of this chapter
will emphasize applications of the results of fluid mechanics theory as they pertain to the
use of wind tunnels for instructional purposes.
www.intechopen.com
Wind Tunnels in Engineering Education 241
2 ( p0 - p)
V= (1)
r
Fig. 5. Sketch of the wind tunnel used (Courtesy of Joseph Thomas, 2006)
Students often wonder whether or not the use of a wall tap is correct; that is, if it can be
justified using analysis. And the answer is that it is and it can. The use of a wall tap is
allowable because the flow is presumed, and is in fact, essentially parallel. An illustration of
parallel flow is shown in Figure 7.
Under parallel-flow conditions, Eq.(2) , which is Euler’s equation, written along a coordinate
axis that is normal to the local streamline, indicates that the curvature of the local
streamlines is extremely large, which causes the pressure gradient in the direction
perpendicular to the streamlines to be zero, making the pressure constant in the direction
www.intechopen.com
242 Wind Tunnels and Experimental Fluid Dynamics Research
normal to the streamline. Therefore, the value of the static pressure measured by the wall
tap is the same as that which would have been measured at the tip of the stagnation probe.
Fig. 6. Two different ways to measure the total and static pressures inside the test section
(Pritchard, 2011).
Fig. 7. Illustration showing that the radius of curvature becomes very large inside the test
section (Pritchard, 2011).
æV 2 ÷ö ¶ p 1 ¶p V2
r ççç ÷÷ = = ;.
÷
çè R ø ¶ n r ¶n R
(2)
¶p
Parallel - flow : R ¥ , 0
¶n
www.intechopen.com
Wind Tunnels in Engineering Education 243
Fig. 8. Variation of wind speed in the test section with the size of the air gap.
4.2 Experimental verification of the existence of the boundary layer over a flat plate
4.2.1 Purpose
The purpose of this experiment is to learn how to use the wind tunnel to measure the
difference between the stagnation (total) pressure and the static pressure at a series of points
located on a vertical line selected in the flow field and use those differences to compute the
wind speeds at each such point using Bernoulli’s equation. The plotting of the resulting
velocity profile and its examination will be used to determine whether or not the existence
of the boundary layer can be detected.
Many viscous flows past solid bodies can be analyzed by dividing the flow region into two
subregions: one that is adjacent to the body and the other that covers the rest of the flow
field. The influence of viscosity is concentrated, and only important, in the first subregion,
that which is adjacent to the body. The effects of viscosity can be neglected in the second
subregion, that is, outside of the region adjacent to the body. The first region has been
called the boundary layer historically. This phrase is a translation of the German phrase
used by Prandtl , who introduced this concept. A big problem in fluid mechanics is locating
the line the demarcates the boundary between the two subregions. Locating this line is also
called determination of the boundary layer, or simply the boundary-layer problem.
The symbol used for the local thickness of this boundary layer is δ. It denotes the distance
between a point on the solid body and the point beyond which the effect of viscosity can be
considered to be negligible.
www.intechopen.com
244 Wind Tunnels and Experimental Fluid Dynamics Research
Blasius (Pritchard, 2011). In conformity with his work, the continuity equation and the
Navier-Stokes equation with the corresponding boundary conditions are ordinarily written
as shown below:
u v
0
x y
u u 2u
v
x y y2
u
(3)
u( y 0) 0; v( y 0) 0
u
u( y ) U ; ( y ) 0
y
Where u is the component of velocity along the plate and v is the component of velocity
perpendicular the plate. The origin of the coordinate system is at the leading edge of the
plate, with the x direction along the plate and the y direction perpendicular to it. The
magnitude of free-stream velocity, far from the plate is U.
Using similarity transformations, one introduces a change of variables as shown below. Let
y u
hµ = g(h );
d U
nx U
dµ h=y (4)
U nx
¶y ¶y y
u= ;v = - ; f (h ) =
¶y ¶x n xU
Applying this change of variables allows the second-order partial differential equation given
above to become a nonlinear, third-order, ordinary differential equation, with the associated
boundary conditions shown below:
d3 f d2 f
2 +f
dh 3 dh 2
df
f ( h = 0) = 0; (h = 0) = 0 (5)
dh
df
(h ¥) = 1.
dh
The solution to this equation is obtained numerically. From that numerical solution, it is
seen that, at = 5.0, u/U = 0.992. If the boundary layer thickness is defined as the value of y
for which u/U = 0.99, one gets
5.0 x Ux
d @= ; with Re x = (6)
Re x n
Using boundary-layer theory, a sketch of the velocity profile along a vertical line in the test
section of the wind tunnel is expected to look as shown below. In this application of the
www.intechopen.com
Wind Tunnels in Engineering Education 245
wind tunnel, one wishes to compare this profile to that obtained experimentally in the test
section of the wind tunnel (See Figure 9.).
δ
U u
Fig. 9. Graphical Representation of Boundary Layer Theory in Wind Tunnel Test Section
www.intechopen.com
246 Wind Tunnels and Experimental Fluid Dynamics Research
42.478
44.043
44.927
Height of point (m)
45.445
45.993
46.355
46.697
46.751
46.768
Air speed (m/s) 46.786
Fig. 10. Experimental velocity profile of the flow in the test section for a speed of 46 m/s
4.3 Determination and characterization of the boundary layer along a flat plate
4.3.1 Purpose
The purpose of this experiment is to learn how to use experimental data collected in a wind
tunnel to determine the thickness of the boundary layer and to characterize the type of
boundary layer that is represented by such data.
Ux
Re x = (7)
n
Where U is the freestream velocity, n is the kinematic viscosity, and x is the distance from
the leading edge of the plate to some point of interest. By convention, a boundary player
becomes turbulent when the Reynolds number of the flow exceeds 5 x 105.
For flow inside the test section, where, instead of inserting a plate, it is the bottom surface of
the test section that takes the role of the flat plate, the location of the leading edge of the
plate must be estimated. In the case of the data reported here, it was estimated in the
following way: the leading edge was defined as the line where the curved section that
constitutes the intake of the tunnel becomes horizontal, and hence, tangential to the inlet to
the test section. In the wind tunnel used for these tests, that line was at a distance of 0.356 m
< x < 0.457 m from the plane that passes through the geometric center (centroid) of the test
section. The free stream speed was U = 46 m/s; and using a kinematic viscosity of 15.68 x
www.intechopen.com
Wind Tunnels in Engineering Education 247
10-6 m2/s, corresponding to a room temperature of ( 200 C ) on the day of the experiment, the
Reynolds number was found to be
Ux
1.05 ´ 106 < Re x = < 1.40 ´ 106 (8)
n
Ux
From these results, 5 ´ 10 5 < 1.05 ´ 106 < Re x = < 1.40 ´ 106 , and the boundary layer is
n
turbulent. Therefore, the boundary layer thickness given by the Blasius solution is not
applicable. Instead, one can use the momentum- integral equation to estimate the thickness
of the boundary layer. That equation is given by
tw
r
=
d
dx
( )
U 2 q + d *U
dU
dx
(9)
Where t w is the shear stress at the wall, q is the momentum thickness, and d * is the
displacement thickness. These are defined as follows:
The displacement thickness: δ*
¥æ ö dæ uö
d* = ò ç1 - u ÷÷ dy = çç1 - ÷÷÷ dy
0 ççè Uø ÷ ò0 çè Uø
(10)
¥ u æç uö d uæ uö
q=ò ç1 - ÷÷÷ dy = ò0 ççç1 - ÷÷÷ dy
ç
(11)
0 Uè Uø Uè Uø
And the boundary-layer thickness is δ.
Since the expression for the flow inside the boundary layer is not known, one uses the
power-law formula for pipe flow but adjusted to boundary layer flows. It is
1/n
u æç y ÷ö
=ç ÷ = h 1/n (12)
U çè d ÷ø
Where n is unknown. Using different values of n allows one to construct a Table such as the
one shown below (Njock Libii, 2010).
u q d* d* d
b = C f (Re x )
1/5
(Re x )
1/5
= h 1/n H= aº
U d d q x
h 1/6 0.107143 0.142857 1. 333333 0.337345906 0.057830727
h 1/7
0.097222 0.125 1.285714 0.381143751 0.059289028
h 1/8
0.088889 0.111111 1.250000 0.423532215 0.060235693
h 1/9
0.0818181 0.10000 1.222222 0.464755563 0.060840728
h 1/10 0.0757575 0.9090909 1.200000 0.504990077 0.061210918
Table 2. Turbulent boundary layer over a flat plate at zero incidence: results
www.intechopen.com
248 Wind Tunnels and Experimental Fluid Dynamics Research
0.382 x Ux
d@ ; with Re x = (14)
(Re x )
1/5
n
In the case tested here, where U = 46 m/s , 0.356 m < x < 0.457 m, and for which
Ux
1.05 ´ 106 < Re x = < 1.40 ´ 106 , the thickness of the boundary layer is found to be
n
Pmax - Pmin
Fluctuation º ( )´100 (15a)
Pmax
www.intechopen.com
Wind Tunnels in Engineering Education 249
Fig. 11. Pressure fluctuations inside the test section at an air speed of 46 m/s.
éæ a2 ÷ö ù éæ a2 ÷ö ù
V = v êê ççç1 - 2 ÷÷ cos q úú r - v êêççç1 + 2 ÷÷ sin q úú q (16)
êë çè r ÷ø úû êëçè r ÷ø úû
www.intechopen.com
250 Wind Tunnels and Experimental Fluid Dynamics Research
Using Bernoulli’s equation, the pressure distribution on the surface of the cylinder is found
to be
Pc - P¥ v 2
r
=
2
(
1 - 4 sin 2 q ) (17)
PC - P¥
Cp = = 1 - 4 sin 2 q (18)
1 r v2
2
Unfortunately, an expression of Cp vs. similar to the one in equation (18) cannot be
obtained analytically for the steady flow of an incompressible viscous fluid. However, one
function of the position, . Such data can be plotted and compared with the graph of
can obtain experimental data by measuring the pressure on the surface of the cylinder as a
Fig. 12. Comparison of Cp for: inviscid and viscous flows (Pritchard, 2011).
www.intechopen.com
Wind Tunnels in Engineering Education 251
CL
FL
V 2 A
(19)
1
2
www.intechopen.com
252 Wind Tunnels and Experimental Fluid Dynamics Research
Where FL is the lift force, V is the average speed, and A is the reference area.
The drag coefficient, CD , is given by
CD
FD
V 2 A
(20)
1
2
Where FD is the drag force, V is the average speed, and A is the reference area.
www.intechopen.com
Wind Tunnels in Engineering Education 253
Fig. 14. Streamlines around a NACA 0012 airfoil at a moderate angle of attack.
Fig. 15. Schematic representation of lift, drag, thrust, and weight on an airfoil.
www.intechopen.com
254 Wind Tunnels and Experimental Fluid Dynamics Research
Fig. 16. Sample variation of lift and drag on an airfoil with the angle of attack.
(p r v2 d 2 )
FD = CD (21)
8
From Eq.(21), if a rough sphere and a smooth sphere of equal diameters are in the same
fluid, the differences in their drag forces would be borne by the drag coefficient, CD .
If a sphere is suspended in the test section of the wind tunnel in such a way as to create a
simple pendulum that makes and angle q with the vertical, and for which the sphere is the
www.intechopen.com
Wind Tunnels in Engineering Education 255
bob, then, neglecting the lift force, the drag force is related to the deflection of the pendulum
by
FD = mg tan q (22)
Using Eq. (22), the drag force per unit weight becomes
FD
= tan q (23)
mg
If one uses a rough sphere and a smooth sphere of equal weights, then, the effect of the
roughness on the drag force would be borne by the deflection angle.
www.intechopen.com
256 Wind Tunnels and Experimental Fluid Dynamics Research
assumes that the flow is entirely in the axial direction , and its velocity is only a function of
two variables, the radial position, r ,and the time, t. If one further assumes that the pipe is
– Kcos(t ) , then, the solution of the thus-simplified Navier-Stokes equations for the
horizontal with internal radius r0 and that the pressure gradient in the axial direction is
i
J 0 r
u e it 1
i
K
i
(24)
J 0 r0
Fig. 17. Deflection angle of the pendulum in the test section vs. mass of the sphere.
Where J0 denotes the Bessel function of the first kind and zero order, K is a constant, is
the mass density of the fluid and is the circular frequency of the oscillation.
One introduces non-dimensional quantities such that
r 2
r * ; * 0 ; u*
r u
(25)
r0 umax
that has a constant pressure gradient equal to – K . It is known that oscillating flows of this
Where umax is the maximum centerline velocity for a steady Poiseuille flow in the same pipe
type may become turbulent when * exceeds 2000. For large values of * a region of large
velocity is established near the wall of the pipe and the mean square velocity is given by
www.intechopen.com
Wind Tunnels in Engineering Education 257
e cos B *
e 2 B
1
u2 2 B
K / 2
2 2 (26)
r * r
where
B 1 r * *
2
(27)
u2
K / 2 2
When this equation is plotted as B vs. 2 , one observes that the velocity displays an
overshoot near the wall, indicating that the maximum velocity occurs near the wall and not
at the center of the pipe , as is the case for Poiseuille flow in a circular pipe. Such a plot is
shown in Figure 18.
Fig. 18. Overshoot of velocity near the wall, known as Richardson’s annular effect
A similar phenomenon was observed in the velocity profile of flow in the wind tunnel that
was shown in Figure 10.
www.intechopen.com
258 Wind Tunnels and Experimental Fluid Dynamics Research
performed calculations for low and high frequency regimes ( 1 h 20 ) in a rectangular
sectional dimensions a and h was solved numerically (Yakhot et al., 1998,1999). Yakhot et al.
www.intechopen.com
Wind Tunnels in Engineering Education 259
duct with different aspect ratios (a/h = 1 and 10). Their results clearly demonstrate the
existence of Richardson’s annular effect, Figure 20. The results in that Figure are very
similar to those that were obtained from the wind tunnel and are shown in Figure 19.
gradient variation with time. (b) Duct flow, h 8 : solid line, x/a =0.5; dashed, x/a = 0.25;
Fig. 20. Velocity profiles in pulsating flow at different instants of one period. (a) Pressure
dot-dashed, x/a =0.1. (c) Flow between two parallel plates (Yakhot et al.,1999).
5. References
Ernst, E. (1983). The Undergraduate Engineering Laboratory, Proceedings of an Engineering
Foundation Conference, Edward W. Ernst, Editor, ISBN: 0939204215, New England
College, Henniker, New Hampshire, July 1983
Feisel, L., Rosa, A. (2005) The role of the engineering laboratory, Journal of Engineering
Education, Vol. (January 2005), pp. (121-130), ISBN 0-309-08291-9.
Njock Libii, J. (2006) Design of an experiment to test the effect of dimples on the magnitude
of the drag force on a golf ball. World Transactions on Engineering and Technology
Education (WTE&TE), Vol.5, No. 3, (December 2006), pp. (477-480), ISSN 1446-2257
Njock Libii, J. (2007) Dimples and drag: Experimental demonstration of the aerodynamics of
golf balls. American Journal of Physics, Vol. 75, No. 8, (August 2007), pp. (764 -767),
ISSN 0002-9505.
www.intechopen.com
260 Wind Tunnels and Experimental Fluid Dynamics Research
Njock Libii, J. (2010) Laboratory exercises to study viscous boundary layers in the test
section of an open-circuit wind tunnel, World Transactions on Engineering and
Technology Education(WTE&TE), Vol. 8, No. 1, (March 2010), pp. (91-97), ISSN 1446-
2257.
Njock Libii, J. (2010) Using wind tunnel tests to study pressure distributions around a bluff
body: the case of a circular cylinder, World Transactions on Engineering and
Technology Education(WTE&TE), Vol. 8, No. 3, (December 2010), pp.( 361-367), ISSN
1446-2257
Njock Libii, J. & Russell R. Tobias, R. (2006). Lifting Bodies, In: USA in Space (Third Edition),
Russell R. Tobias and David G. Fisher, pp. (744-748), Salem Press, ISBN: 1-58765-
259-5, Pasadena, California
Pritchard, P. (2011). Fox and McDonald’s Introduction to Fluid Mechanics (Eighth edition ),
John Wiley & Sons, ISBN-13 9780470547557, ISBN-10 0470547553, New York, NY.
Rutgers, http://coewww.rutgers.edu/classes/mae/mae433/lab3.pdf, n.d.
Sonntag, R., Borgnakke, C., & Van Wylen, G. (1998). Fundamentals of Thermodynamics (Fifth
Edition), John Wiley & Sons, Inc., ISBN 0-471-18361-X, New York, 1998.
Yakhot, A., M. Arad, M., & Ben-Dor, G. (1998). Richardson's Annular Effect in Oscillating
Laminar Duct Flows. Journal of Fluids Engineering, Vol. 120, 1, (March 1998) pp.
(209-301), ISSN 0098-2202
Yakhot, A., M. Arad, M., & Ben-Dor, G. (1999). Numerical investigation of a laminar
pulsating flow in a rectangular duct. International Journal For Numerical Methods In
Fluids, Vol. 29, (1999) , pp. (935-950), ISSN 0271-2091
www.intechopen.com
Wind Tunnels and Experimental Fluid Dynamics Research
Edited by Prof. Jorge Colman Lerner
ISBN 978-953-307-623-2
Hard cover, 709 pages
Publisher InTech
Published online 27, July, 2011
Published in print edition July, 2011
The book “Wind Tunnels and Experimental Fluid Dynamics Research†is comprised of 33 chapters
divided in five sections. The first 12 chapters discuss wind tunnel facilities and experiments in incompressible
flow, while the next seven chapters deal with building dynamics, flow control and fluid mechanics. Third section
of the book is dedicated to chapters discussing aerodynamic field measurements and real full scale analysis
(chapters 20-22). Chapters in the last two sections deal with turbulent structure analysis (chapters 23-25) and
wind tunnels in compressible flow (chapters 26-33). Contributions from a large number of international experts
make this publication a highly valuable resource in wind tunnels and fluid dynamics field of research.
How to reference
In order to correctly reference this scholarly work, feel free to copy and paste the following:
Josué Njock Libii (2011). Wind Tunnels in Engineering Education, Wind Tunnels and Experimental Fluid
Dynamics Research, Prof. Jorge Colman Lerner (Ed.), ISBN: 978-953-307-623-2, InTech, Available from:
http://www.intechopen.com/books/wind-tunnels-and-experimental-fluid-dynamics-research/wind-tunnels-in-
engineering-education