Learning: Wind Tunnel Familiarization
Learning: Wind Tunnel Familiarization
LEARNING
MODULE 01:
Wind Tunnel
Familiarization
AE 324
AERONAUTICAL LABORATORY
TABLE OF CONTENTS
• Diffusers 4 18
• Driers 4 20
• Shock Tube 4 23
Assessment Activity 24
TABLE OF FIGURES
FIGURE REFERENCE
Tecquipment Academia. Subsonic Wind Tunnel 450 mm. Retrieved from:
1.1
https://www.tecquipment.com/subsonic-wind-tunnel-450mm
National Aeronautics and Space Administration. Wind Tunnel Design.
1.2 Retrieved from: https://www.grc.nasa.gov/www/k-
12/airplane/tunnozd.html
Lockheed Martin. Lockheed Martin Missiles and Fire Control’s High
1.3 Speed Wind Tunnel (HSWT). Retrieved from:
https://www.lockheedmartin.com/en-us/products/highspeedwindtunnel.html
Pope A., Goin K. High Speed Wind Tunnel Testing. New York: John
1.7-1.10
Wiley & Sons.
1.4-1.6 Clancy, L. J. (1975). Aerodynamics [E-book]. Pitman Publishing Limited.
1.11-1.17
TABLE OF REFERENCES
REFERENCES No.
Pope A., Goin K. High Speed Wind Tunnel Testing. New York: John Wiley &
1
Sons.
Lockheed Martin Missiles and Fire Control. High Speed Wind Tunnel and
Test Systems Design Handbook [document on the Internet]. 2002 May 06.
2
[cited 2009 November 15]. Retrieved from:
http://www.lockheedmartin.com/data/assets/13617.pdf
National Aeronautics and Space Administration. Wind Tunnel Index.
Retrieved October 1, 2020 from https://www.grc.nasa.gov/www/k- 3
12/airplane/shortt.html
Clancy, L. J. (1975). Aerodynamics [E-book]. Pitman Publishing Limited. 4
HONESTY CLAUSE
LEARNING OUTCOMES
I. INTRODUCTION
Wind tunnels are designed for a specific purpose and speed range and there
is a wide variety of wind tunnel types and model instrumentation. The model to be
tested in the wind tunnel is placed in the test section of the tunnel. The speed in
the test section is determined by the design of the tunnel. The choice of speed
range affects the design of the wind tunnel due to compressibility effects.
For subsonic flows, the air density remains nearly constant and decreasing
the cross-sectional area causes the flow to increase velocity and decrease
pressure. Similarly, increasing the area causes the velocity to decrease and the
pressure to increase. We want the highest possible velocity in the test section. For
a subsonic wind tunnel, the test section is placed at the end of the contraction
section and upstream of the diffuser. From a knowledge of the conservation of mass
for subsonic flows, we can design the test section to produce a desired velocity or
Mach number since the velocity is a function of the cross-sectional area. On the
figure, we note the changes in Mach number, velocity, and pressure through a
subsonic wind tunnel design. The plenum is the settling chamber on a closed return
tunnel, or the open room of an open return design.
For supersonic flows, the air density changes in the tunnel because of
compressibility. In fact, the density changes faster than the velocity by a factor of
the square of the Mach number. In a supersonic flow, decreasing the cross-
sectional area causes the flow to decrease in velocity and increase pressure.
Similarly, increasing the area causes the velocity to increase and the pressure to
decrease. This change in properties is exactly the opposite of the change that
occurs subsonically. In addition, compressible flows experience mass flow choking.
As a subsonic flow is contracted, the velocity and Mach number increase. When
the velocity reaches the speed of sound (M = 1), the flow chokes and the Mach
number cannot be increased beyond M = 1. We want the highest possible velocity
in the test section of the wind tunnel. For a supersonic wind tunnel, we contract the
flow until it chokes in the throat of a nozzle. We then diffuse the flow which increases
the speed supersonically. The test section of the supersonic tunnel is placed at the
end of the diffuser. From a consideration of conservation of mass for a compressible
flow, we can design the test section to produce a desired velocity or Mach based
on the area in the test section. On the figure we note the changes in Mach number,
velocity and pressure through a supersonic wind tunnel design.
High-speed tunnels are those whose usual operating speeds require the
inclusion of compressible flow effects. This, it turns out, means that in the high-
speed field we usually talk about “Mach number”—the ratio of a given velocity to
the speed of sound in the air about the body—as a more typical parameter than
velocity. A lower limit of “high speed” might be considered to be where the Mach
number is approximately 0.5—about 380 mph for standard sea level conditions.
The power to drive a low-speed wind tunnel varies as the cube of the velocity
in the wind tunnel. Although this rule does not hold into the high-speed regime, the
implication of rapidly increasing power requirements with increasing test speed is
correct. Because of the power requirements, high-speed wind tunnels are often of
the “intermittent” type, in which energy is stored in the form of pressure or vacuum
or both and is allowed to drive the tunnel only a few seconds out of each pumping
hour.
Figure 1.3 Lockheed Martin Missiles and Fire Control’s High-Speed Wind Tunnel
The problems of nozzle and diffuser design are common to all supersonic
wind tunnels. So is the problem of providing the necessary pressure ratio, and thus
the necessary power, to achieve supersonic flow. This latter problem is solved in
very different ways in different tunnels. The different designs include tunnels for
intermittent operation, and continuously running tunnels; tunnels in which the flow
goes straight through and return circuit tunnels: tunnels in which the flow is direct,
and induced flow tunnels. We shall discuss each type in turn.
A. Intermittent Operation
However, the running time is generally short, because of the limited amount
of air in the reservoir. Once the pressure in the reservoir drops below a certain
value, the pressure ratio across the tunnel will no longer be sufficient to give
supersonic flow in the working section. It is then necessary to stop the tunnel and
re-charge the reservoir, and the charging time between runs may be considerable.
A further disadvantage is that the stagnation temperature and pressure will vary
during a run, due to the expansion of the air in the reservoir, where the flow
originates.
C. Continuous Operation
RETURN CIRCUITS
CHOKING
The phenomenon of choking occurs when the local Mach number in a high-
speed tunnel reaches unity at some point in the flow other than the nozzle throat.
This may happen at points where the effective area of cross-section of the tunnel
is reduced as a result of boundary layer thickening, shock-induced separation, or
simply the presence of a model which is too thick, or at too high an incidence.
When choking occurs, the supersonic flow in the working section breaks down, and
this condition sets a limit to the range of operation of the tunnel-a limit which may
vary according to the model under test, and its attitude.
1. We are more in control of conditions and may usually return to a given flow
condition with more accuracy.
2. Since the “panic” of rapid testing is removed, check points are more easily
obtained.
3. Testing conditions can be held constant over a long period of time.
Some advantages of the indraft tunnel over the blowdown tunnel are:
The detail shape of the nozzle is a matter for very careful design, but the
theory on which it is based is beyond the scope of this book. The problem arises
from the fact that as the flow accelerates super sonically expansion waves
emanate from points on the tunnel walls. When these waves strike the opposite
walls, they may be reflected as shock waves, which would spoil the flow in the
working section. The shape of the liner walls must be calculated so as to avoid
this wave reflection,
Some liners are single-sided, i.e., one wall of the liner is plane, as
indicated in Fig. 1.12. This makes the replacement of one liner by another much
easier and cheaper, but it may be more difficult to achieve uniform flow. The
working section Mach number is only achieved downstream of an oblique shock
wave located across the tunnel and generated by the concavity of the liner wall.
(b) The visualization of the flow by optical methods involves the insertion of
plane windows in the side walls of the working section.
(c) When models are mounted in the working section, and tested, shock
waves are generated by the model. These shocks are inevitably reflected
by the tunnel walls, and it is imperative that the reflected shocks should
not strike the model. For this reason, at any given Mach number, the
working section must have a minimum height in relation to the model
length or chord. This would be impossible over the whole span of a two-
dimensional model if the section were circular.
• DIFFUSERS
The diffusion of the flow, i.e., the slowing down of the air after it has left
the working section, in the case of a supersonic tunnel has to be affected in two
stages. A supersonic diffuser is required in which the flow is slowed down from
supersonic to subsonic speeds, and this must be followed by a subsonic
diffuser in which the flow is reduced to the low speeds at which the air is
discharged to the atmosphere or sent on its return circuit. There is no reason
why the subsonic diffuser should be any different from one used in a low speed
tunnel, except that it will probably need to be longer, since it must affect a bigger
speed reduction. It is in the diffuser that most of the power losses occur, and
therefore the design of the diffuser, and particularly the supersonic stage, is a
vital feature in keeping down the power requirement. There are various types
of supersonic diffuser.
Figure 1.15 illustrates this concept and shows how the wedge design must take
account of shock reflection
• DRIERS
The design of such liners is based on a working section whose walls are
not solid but perforated. Slots or holes are made in the walls of the working
section, the area of the perforations increasing with distance along the length
of the section. The principle is that the excess air associated with boundary
layer growth or the presence of the model can escape through these
perforations, thus eliminating the effect of a narrowing channel, and preventing
choking. The detail positioning, size and shape of the holes or slots is often
crucial to the successful operation of the tunnel, and it is largely a matter of trial
and error. The perforated walls are surrounded by other, solid walls, and the air
which has escaped through the holes into the expansion chamber between,
which is known as the plenum chamber, is fed back into the tunnel downstream
of the working section Figure 1.16 illustrates this concept.
Figure 1.16
The shock tube is a device for producing a flow at very high Mach
number for a very short time. The principle is illustrated in Fig. 1.17. The device
consists simply of a long tube, in which a diaphragm separates a region in which
the air is compressed to a very high pressure, from an evacuated region. When
the pressure ratio across it reaches a certain value, the diaphragm bursts, (or
it may be punctured by a device incorporated in the tube) and air rushes at very
high speeds into the evacuated region. At the front of this body of air is a mixing
region, terminated by a normal shock, behind which is a region of uniform flow
at high Mach number. This passes over a model mounted at an appropriate
point in the initially evacuated region. When the wave front reaches the end of
the tube, it is reflected back again, and once it has reached the model position
the flow is spoilt. Thus, the time for which uniform flow is achieved is
exceedingly short, probably only a few micro-seconds. Measurement of the flow
properties is very difficult, and many of the problems created by the use of
shock tubes have consisted in the development of special instrumentation for
making such measurements on a very small-time scale. The tube itself is cheap
and easy to construct, and the power requirement is low. More complicated and
refined developments of the idea, such as the gun-tunnel, in which a charge of
air is fired along the tube, have been developed, and the end of the tube may
fork into two branches so as to give increased running time', but the above
account describes the basic principle of these devices.
The shock tube and gun-tunnel are experimental tools for the
investigation of flow at hypersonic speeds. Since hypersonic aero dynamics lies
outside the scope of this book, is not intended to deal further here with
ASSESSMENT ACTIVITY
Answer the following questions in your notebook.
AE 324 – AERONAUTICAL
LABORATORY II
1|P a ge
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
TABLE OF CONTENTS
Title Ref. no Page
Classification of Aerodynamic Flows 1 5
Equations Governing Incompressible Flow 1 5
Equations Governing Compressible Flow 1,6 9
Derivation of Speed of Sound 1 12
Derivation of Important Isentropic Relations and Area
1,6 14
Relations for High Speed Flows
Sample Problems 18
Assessment Activity
2|P a ge
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
TABLE OF REFERENCES
References No.
Anderson, Jr., J. D. (1989). Introduction to Flight (3rd Edition) [E-book].
1
McGraw-Hill Book Company
Jones, B. (1939). Elements of Practical Aerodynamics (2nd Edition) [E-
2
book]. John Wiley & Sons, Inc.
Clancy, L. J. (1975). Aerodynamics [E-book]. Pitman Publishing Limited. 3
Sadraey, M. H. (2013). Aircraft Design A Systems Engineering Approach.
4
John Wiley & Sons, Ltd.
FAA. (2016). Pilot’s Handbook of Aeronautical Knowledge. 5
Pope, A., & Goin, K. L. (1965). High Speed Wind Tunnel Testing [E-book].
6
John Wiley & Sons, Inc.
Figures and Tables No.
Figures 2.1-2.3 are taken from the reference above. 1
Note that tables presented are made for organizing purposes only.
3|P a ge
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
LEARNING OUTCOMES
Program Learning Outcomes (PLO)
Course Learning Outcomes (CLO)
PLO 1. Apply knowledge of Module Learning Outcomes (MLO)
mathematics, physical sciences,
CLO 2: Determine flow Topic Learning
engineering sciences to the
properties along surfaces using Outcomes (TLO)
practice of aeronautical MLO 1: Formulate equations
engineering. formulas taken from isentropic,
adiabatic and thermodynamics for high speed processes TLO 6: Review basic
PLO 5.Identify, formulate, and
solve aeronautical engineering for high speed flight. using isentropic, adiabatic thermodynamics using
problems. CLO 7: Commit to the and thermodynamic existing formulas
PLO 6. Understand the scientific process on solving properties. applicable for
professional and ethical problems by producing MLO 2: Solve flow aerodynamics.
responsibility solutions following the basic properties for high speed TLO 7: Explain
community and with society at steps of scientific method. isentropic properties
large
processes using derived
CLO 8: Display professional using formulas from
PLO 8.Understand the impact of working equations.
commitment to ethical practice thermodynamics.
aeronautical engineering solutions TLO 8: Derivation of
in global, economic, by showing sense of
responsibility in complying and working formulas for
environmental, and societal
context compiling academic high speed processes
PLO 9.Recognize the need for, requirement. using isentropic
and engage in life-long learning relations.
PLO 11.Use techniques, skills, TLO 9: Solve
and modern engineering tools problems involving
necessary for aeronautical high speed process
engineering practice through problem sets
PLO 13.Apply acquired and/or module
aeronautical engineering
activities.
knowledge and skills for national
development.
4|P a ge
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
CONTENT
In reality, all matter is compressible. It will change its volume but its mass will remain
the same. As an example, solids and liquids will change its volume when compressed
and its mass will stay the same but its density will change. For all practical purposes,
this change in density is insignificantly small thus the density is constant. Meanwhile,
for a gas, there will be a change in volume if compressed but the density can be
variable.
I. Classification of Aerodynamic Flows
1. Incompressible Flows
- Flows where the density of the fluid is always constant. In reality,
incompressible flows are a myth. They do not exist in nature but
they are classified as such because at the speeds where this fluid
flows the change in density is significantly small it is assumed
constant for practicality. Speeds classified under incompressible
flows are V < 100 m/s.
2. Compressible Flows
- Flows in which the change in density is not a constant. This
variability of density is particularly important in high speed flows
where the change in density is not small enough to be neglected.
By convention, if the speed exceeds 100 m/s it is considered
compressible.
NOTE: It is important that the student is able to differentiate the two flows
because it will be important for the computation process. If they
interchange the two type of flows, there will be a huge discrepancy in the
final answer.
II. Equations Governing Incompressible Flows
1. Equation of State
- At any point in the flow field, the flow properties such as pressure,
temperature and density are related by the equation of state given
as:
P= ρRT
Where:
P = Pressure
ρ = Density
T = Temperature
R = specific gas constant
2. Continuity Equation
- From the physical principle of conservation of mass which states
that “mass can neither be created nor destroyed.”
5|P a ge
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
Area A1 A2
Density ρ1 ρ2
Velocity V1 V2
Distance Covered V1 dt V2 dt
Volume Swept at a
A1 V1 dt A2 V2 dt
time interval dt
Differential mass dm = ρ1 (A1 V1 dt) dm = ρ2 (A2 V2 dt)
𝑑𝑚 𝑑𝑚
𝑚1 = = ρ1 A1 V1 𝑚2 = = ρ2 A2 V2
Mass Flow Formula 𝑑𝑡 𝑑𝑡
Continuity Equation
and from the
definition that mass 𝑚1 = 𝑚2 = 𝐶
can neither be ρ1 A1 V1 = ρ2 A2 V2
created nor
destroyed.
Continuity Equation A1 V1 = A2 V2 𝑠𝑖𝑛𝑐𝑒 𝑑𝑒𝑛𝑠𝑖𝑡𝑦 𝑖𝑠 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡
(Incompressible)
6|P a ge
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
3. Momentum Equation
- The momentum equation shows the relationship of pressure in the
aerodynamic flow. Indeed, the differences in pressure from one
point to another in the flow creates forces that act on the fluid
elements that cause them to move. Therefore, it is important to
determine the relation between pressure and velocity.
7|P a ge
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
𝑑𝑉 𝑑𝑉 𝑑𝑥 𝑑𝑉 𝑑𝑥
𝑎= = ( )= ( )
Taking acceleration and rewriting in 𝑑𝑡 𝑑𝑡 𝑑𝑥 𝑑𝑥 𝑑𝑡
𝑑𝑉
terms of Velocity (V) = 𝑉
𝑑𝑥
𝐹 = 𝑚𝑎
𝑑𝑃
𝐹=− 𝑑𝑥𝑑𝑦𝑑𝑧
𝑑𝑥
𝑚 = 𝜌𝑑𝑥𝑑𝑦𝑑𝑧
Going back to Newton’s Second 𝑑𝑉
Law of Motion 𝑎= 𝑉
𝑑𝑥
𝑑𝑃 𝑑𝑉
− 𝑑𝑥𝑑𝑦𝑑𝑧 = 𝜌𝑑𝑥𝑑𝑦𝑑𝑧 ∗ ( 𝑉)
𝑑𝑥 𝑑𝑥
𝑑𝑃 = −𝜌𝑉𝑑𝑉
Euler’s Equation/Momentum
𝑑𝑃 + 𝜌𝑉𝑑𝑉 = 0
Equation
2 2
Integrating from point 1 to point 2.
∫ 𝑑𝑃 + ∫ 𝜌𝑉𝑑𝑉 = 0
1 1
8|P a ge
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
𝜌
(𝑃2 − 𝑃1 ) + (𝑉22 − 𝑉12 ) = 0
The equation is constant along a 2
𝜌 2 𝜌
streamline. 𝑃2 + (𝑉2 ) = 𝑃1 + (𝑉12 ) = 𝐶
2 2
𝜌 2 𝜌
𝑃2 + (𝑉2 ) = 𝑃1 + (𝑉12 ) = 𝐶
Bernoulli’s Equation 2 2
2. Momentum Equation
- The momentum equation for compressible flow in a stream tube or
duct with no friction is given by:
𝑃2 + 𝜌2 (𝑉22 ) = 𝑃1 + 𝜌1 (𝑉12 ) = 𝐶
Where:
P = Pressure
ρ = Density
V = Velocity
3. Continuity Equation
- For a continuous flow in a duct or stream tube, the total mass flow
at any two points along the flow is given by:
ρ1 A1 V1 = ρ2 A2 V2
ρ1 ≠ ρ2
Where:
ρ = Density which is not constant along the flow
A = Area
V = Velocity
4. Isentropic Equation
- Defining different process from thermodynamics:
9|P a ge
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
STEPS EQUATIONS
𝛿𝑞 = 𝑑𝑢 + 𝑃𝑑𝜈 𝑏𝑢𝑡 𝛿𝑞 = 0
𝑑𝑢 = −𝑃𝑑𝜈 = 𝑐𝑣 𝑑𝑇
Since isentropic flow is an
adiabatic process 𝛿𝑞 = 𝑑ℎ − 𝜈𝑑𝑃 𝑏𝑢𝑡 𝛿𝑞 = 0
𝑑ℎ = 𝜈𝑑𝑃 = 𝑐𝑃 𝑑𝑇
𝑑𝑢 = −𝑃𝑑𝜈 = 𝑐𝑣 𝑑𝑇
Dividing du and dh. Note that v 𝑑ℎ = 𝜈𝑑𝑃 = 𝑐𝑃 𝑑𝑇
denotes Volume. 𝑑𝑃 − 𝑐𝑃 𝑑𝜈
=
𝑃 𝐶𝑣 𝜈
10 | P a g e
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
𝑃2 𝜈2
Integrating at from point 1 to 𝑙𝑛 = −𝑘 ln
𝑃1 𝜈1
point 2 along a streamline 𝑃2 𝜈2 −𝑘
𝑙𝑛 = ln ( )
𝑃1 𝜈1
𝑃2 𝜈2 −𝑘 1
= ( ) 𝐵𝑢𝑡 𝜈 =
𝑃1 𝜈1 𝜌
𝑘
𝑃2 𝜌2 𝑘 𝑇2 𝑘−1
Isentropic Flow Relationships = ( ) =( )
𝑃1 𝜌1 𝑇1
5. Energy Equation
- For high speed, compressible flows the physical principle
goes as follows: Energy can neither be created nor destroyed.
- It was mentioned that for incompressible flows, temperature along
the streamline is constant and therefore temperature is not
considered as an important parameter for equations used in such
flows. On the other hand, the density varies for compressible flows
and therefore, the temperature changes. The temperature changes
despite the process being adiabatic in nature.
STEPS EQUATIONS
Substitution 1
𝑑ℎ + 𝜈𝜌𝑉𝑑𝑉 = 0 𝑏𝑢𝑡 𝜈 =
NOTE THAT V denotes velocity 𝜌
and v denotes volume. 𝑑ℎ + 𝑉𝑑𝑉 = 0
11 | P a g e
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
𝑉22 𝑉12
(ℎ2 − ℎ1 ) + ( − )=0; ℎ
2 2
= 𝑐𝑃 𝑇
As the wave travels from the right side to the left side it imparts changes to
the properties originally ahead of it. Imagine a sound wave travelling the
air which have point properties such as Pressure, Temperature and
Density, as the air passes the wave a slight change occurs. An example is
the air speed a which moves away from the wave with a velocity of a + da.
STEPS EQUATIONS
𝜌1 𝐴1 𝑉1 = 𝜌2 𝐴2 𝑉2
𝐵𝑢𝑡 𝑉1 = 𝑎; 𝑉2 = 𝑎 + 𝑑𝑎; 𝜌2
= 𝜌 + 𝑑𝜌
𝜌𝑎 = (𝜌 + 𝑑𝜌)(𝑎 + 𝑑𝑎)
𝜌𝑎 = (𝜌 + 𝑑𝜌)(𝑎 + 𝑑𝑎)
Simplifying 𝜌𝑎 = (𝜌𝑎 + 𝑎𝑑𝜌) + 𝜌𝑑𝑎 + 𝑑𝜌𝑑𝑎
𝐵𝑈𝑇 𝑑𝜌𝑑𝑎 𝑖𝑠 𝑛𝑒𝑔𝑙𝑖𝑔𝑖𝑏𝑙𝑦 𝑠𝑚𝑎𝑙𝑙
12 | P a g e
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
𝜌𝑑𝑎
𝑎= −
𝑑𝜌
𝑑𝑃 = −𝜌𝑉𝑑𝑉 𝐵𝑢𝑡 𝑉 = 𝑎
𝑑𝑃 = −𝜌𝑎𝑑𝑎
From Euler’s Equation 𝑑𝑃
𝑑𝑎 = −
𝜌𝑎
𝜌𝑑𝑎
𝑎= −
𝑑𝜌
𝜌𝑑𝑃
Substitute da 𝑎=
𝜌𝑎𝑑𝜌
𝑑𝑃
𝑎2 =
𝑑𝜌
𝑑𝑃 𝑃
= 𝑘 (𝑘 ∗ 𝜌𝑘−1 )
𝑑𝜌 𝜌
Substituting the value of C 𝑑𝑃 𝑘𝑃
=
𝑑𝜌 𝜌
𝑑𝑃
𝑎=√
𝑑𝜌
Speed of sound
𝑘𝑃
𝑎=√ = √𝑘𝑅𝑇
𝜌
Note that speed of sound varies directly with the square root of
temperature, i.e. the speed of sound of a gas depends only on the
temperature of the gas. As defined before, temperature varies directly to
the mean average kinetic energy from the molecular collisions happening
in the gas. As a supporting evidence, the sound waves transfers some of
the molecular kinetic energy to the air and this results in a new average
13 | P a g e
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
𝑉 = 𝑀𝑎
𝑎 = √𝑘𝑅𝑇
𝑘𝑅
𝑐𝑝 =
𝑘−1
𝑘
𝑃2 𝜌2 𝑘 𝑇2 𝑘−1
= ( ) =( )
𝑃1 𝜌1 𝑇1
Applying stagnation conditions on the left hand side of the energy equation
where V = 0.
𝑉12 𝑉22
𝑐𝑝 𝑇1 + = 𝑐𝑝 𝑇2 +
2 2
𝑉12
𝑐𝑝 𝑇1 + = 𝑐𝑝 𝑇0 𝑏𝑢𝑡 𝑉 = 𝑀𝑎
2
𝑀12 𝑎12
𝑐𝑝 𝑇1 + = 𝑐𝑝 𝑇0 𝑏𝑢𝑡 𝑎 = √𝑘𝑅𝑇
2
𝑀12 (𝑘𝑅𝑇1 )
𝑐𝑝 𝑇1 + = 𝑐𝑝 𝑇0 𝑑𝑖𝑣𝑖𝑑𝑖𝑛𝑔 𝑏𝑜𝑡ℎ 𝑠𝑖𝑑𝑒𝑠 𝑏𝑦 𝑐𝑝 𝑇0
2
14 | P a g e
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
𝑐𝑝 𝑇1 𝑀12 (𝑘𝑅𝑇1 ) 𝑐𝑝 𝑇0
+ =
𝑐𝑝 𝑇0 2 (𝑐𝑝 𝑇0 ) 𝑐𝑝 𝑇0
𝑇1 𝑀12 (𝑘𝑅𝑇1 ) 𝑘𝑅
+ = 1 𝑏𝑢𝑡 𝑐𝑝 =
𝑇0 2 (𝑐𝑝 𝑇0 ) 𝑘−1
𝑇1 𝑀12 (𝑘𝑅𝑇1 )
+ =1
𝑇0 2 ( 𝑘𝑅 ∗ 𝑇 )
𝑘−1 0
𝑇1 (𝑘 − 1)𝑀12 (𝑇1 )
+ =1
𝑇0 2 (𝑇0 )
𝑇1 (𝑘 − 1)𝑀12
(1 + )=1
𝑇0 2
𝑇1 1
=
𝑇0 (𝑘 − 1)𝑀12
[1 + ]
2
−1
𝑇1 (𝑘 − 1)𝑀12
[1
= + ]
𝑇0 2
𝑘
𝑃1 (𝑘 − 1)𝑀12 𝑘−1
= [1 + ]
𝑃0 2
𝑘
−1∗
𝜌0 𝑘 (𝑘 − 1)𝑀12 𝑘−1
( ) = [1 + ]
𝜌1 2
𝑘−1
𝜌1 (𝑘 − 1)𝑀12
= [1 + ]
𝜌0 2
15 | P a g e
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
At another point along the stream line consider point two with similar
isentropic relations:
−1
𝑇2 (𝑘 − 1)𝑀22
= [1 + ]
𝑇0 2
𝑘
𝑃2 (𝑘 − 1)𝑀22 𝑘−1
= [1 + ]
𝑃0 2
𝑘−1
𝜌2 (𝑘 − 1)𝑀22
= [1 + ]
𝜌0 2
𝑇1
𝑇01 =
(𝑘 − 1)𝑀12
[1 + ]
2
𝑇2
𝑇02 =
(𝑘 − 1)𝑀22
[1 + ]
2
𝑇1 𝑇2
−1 = −1
(𝑘 − 1)𝑀12 (𝑘 − 1)𝑀22
[1 + ] [1 + ]
2 2
(𝑘 − 1)𝑀22
𝑇1 [1 + ]
2
=
𝑇2 (𝑘 − 1)𝑀12
[1 + ]
2
16 | P a g e
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
1
(𝑘 − 1)𝑀22 𝑘−1
𝜌1 [1 + ]
2
={ }
𝜌2 (𝑘 − 1)𝑀12
[1 + ]
2
A1 ρ2 V2
=
A2 ρ1 𝑉1
A1 ρ2 M2 𝑎2
=
A2 ρ1 𝑀1 𝑎1
A1 ρ2 M2 √𝑘𝑅𝑇2
=
A2 ρ1 𝑀1 √𝑘𝑅𝑇1
A1 M2 ρ2 𝑇2 1/2
= ( )
A2 𝑀1 ρ1 𝑇1
1
− −1/2
(𝑘 − 1)𝑀22 𝑘−1 (𝑘 − 1)𝑀22
A1 M2 [1 + ] [1 + ]
2 2
= { } ∗{ }
A2 𝑀1 (𝑘 − 1)𝑀12 (𝑘 − 1)𝑀12
[1 + ] [1 + ]
2 2
𝑘+1
−
(𝑘 − 1)𝑀22 2(𝑘−1)
A1 M2 [1 + ]
2
= { }
A2 𝑀1 (𝑘 − 1)𝑀12
[1 + ]
2
𝑘+1
−
(𝑘 − 1) 2(𝑘−1)
A1 1 [1 + ]
2
= { }
A∗ 𝑀1 (𝑘 − 1)𝑀12
[1 + ]
2
17 | P a g e
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
𝑘+1
(𝑘 − 1)𝑀12 2(𝑘−1)
A1 1 [1 + ]
2
= { }
A∗ 𝑀1 (𝑘 − 1)
[1 + 2 ]
The equation presented are used for supersonic nozzles whose throat
have a Mach = 1.
18 | P a g e
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
19 | P a g e
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
ANSWER KEY:
1. Consider a convergent duct with an inlet area A1 = 5 m2. Air enters this duct
with a velocity V1 = 10 m/s and leaves the duct exit with a velocity V2 = 30
m/s. What is the area of the duct exit?
Solution:
Since the flow velocities are less than 100 m/s, then we can assume that the
flow is incompressible.
From the continuity equation:
𝜌1 𝐴1 𝑉1 = 𝜌2 𝐴2 𝑉2 𝑏𝑢𝑡 𝜌 = 𝐶
𝐴1 𝑉1 = 𝐴2 𝑉2
𝐴1 𝑉1
𝐴2 =
𝑉2
𝑚
5 𝑚2 ∗ 10 𝑠
𝐴2 = 𝑚
30 𝑠
𝐴2 = 1.67 𝑚2
2. Consider an airfoil in a flow of air, where far ahead of the airfoil, the pressure,
velocity, and density are 2116 psf, 100 mph, and 0.002377 slug/ft. 3,
respectively. At a given point A on the airfoil, the pressure is 2070 psf. What is
the velocity at that point?
Solution:
88 𝑓𝑡/𝑠
𝑉1 = 100 𝑚𝑝ℎ ∗ = 146.67 𝑓𝑡/𝑠
60 𝑚𝑝ℎ
Hence, the flow is incompressible.
From the Bernoulli’s equation, we have:
𝜌1 𝑉12 𝜌𝐴 𝑉𝐴2
𝑃1 + = 𝑃𝐴 + 𝑏𝑢𝑡 𝜌 = 𝐶
2 2
𝜌
𝑃1 − 𝑃𝐴 = [𝑉𝐴2 − 𝑉12 ]
2
1/2
2 2
[ ]
𝑉𝐴 = { 𝑃1 − 𝑃𝐴 + 𝑉1 }
𝜌
1/2
2
𝑉𝐴 = { [2116 𝑝𝑠𝑓 − 2070 𝑝𝑠𝑓] + (146.67 𝑓𝑡/𝑠)2 }
0.002377
𝑉𝐴 = 245.39 𝑓𝑡/𝑠
20 | P a g e
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
3. Consider the same convergent duct and conditions as in problem #1. If the air
pressure and temperature at the inlet are 1.2x10 5 Pa and 330 K, respectively.
Calculate the pressure at the exit.
Solution:
For the inlet conditions, we have the equation of state for the density:
𝑃 1.2 𝑥 105 𝑃𝑎 𝑘𝑔
𝜌= = = 1.27 3
𝑅𝑇 287.08 𝐽 ∗ 330 𝐾 𝑚
𝑘𝑔 𝐾
𝑃2 = 61629.51 𝑃𝑎
5. A jet transport is flying at a standard altitude of 30,000 ft. with a velocity of 550
mph. What is its Mach number?
Solution:
At 30000 ft, we have:
𝑇 = 𝑇𝑜 + 𝑎ℎ
21 | P a g e
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
𝑅
𝑇 = 519 𝑅 − 0.003566 (30000 𝑓𝑡) = 412.02 𝑅
𝑓𝑡
𝑓𝑡
𝑎 = √𝑘𝑅𝑇 = √1.4 ∗ 1716 ∗ 412.02 = 994.91
𝑠
𝑓𝑡
88 𝑠 𝑓𝑡
𝑉 = 550 𝑚𝑝ℎ ∗ = 806.67
60 𝑚𝑝ℎ 𝑠
𝑓𝑡
𝑉 806.67 𝑠
𝑀= = = 0.811
𝑎 994.91 𝑓𝑡
𝑠
6. In a low speed subsonic wind tunnel, one side of a mercury manometer is
connected to the settling chamber (reservoir) and the other side is connected
to the test section. The contraction ratio of the nozzle A2/A1 equals 1/15. The
reservoir pressure and temperature are 1.1 atm and 300 K, respectively.
When the tunnel is running, the height difference between the two columns of
mercury is 10 cm. The density of liquid mercury is 1.36 x10 4 kg/m3. Calculate
the airflow velocity in the test section V2.
Solution:
𝑃1 – 𝑃2 = 𝜌𝑔∆ℎ
𝑘𝑔 𝑚 1𝑚
𝑃1 – 𝑃2 = 1.36𝑥104 3 ∗ 9.81 2 ∗ (10 𝑐𝑚 ∗ )
𝑚 𝑠 100 𝑐𝑚
𝑃1 – 𝑃2 = 13341.6 𝑃𝑎
From the equation of state for the reservoir:
101325 𝑃𝑎
𝑃1 1.1 atm ∗ 1 𝑎𝑡𝑚
𝜌1 = = = 1.29 𝑘𝑔/𝑚3
𝑅𝑇1 (287.08 𝐽 )(300 𝐾)
𝑘𝑔 𝐾
Since we are dealing with low speed subsonic flow, we can assume that it is
an incompressible flow. From the Bernoulli’s Equation, we have: Designating
point 2 as the area corresponding to the test section.
𝜌1 𝑉12 𝜌2 𝑉22
𝑃1 + = 𝑃2 +
2 2
𝜌2 𝑉22 𝜌1 𝑉12
𝑃1 − 𝑃2 = −
2 2
𝜌 2
𝑃1 − 𝑃2 = [𝑉2 − 𝑉12 ]
2
2
𝑉22 − 𝑉12 = (𝑃1 − 𝑃2 )
𝜌
From the continuity equation, we have:
𝐴1 𝑉1 = 𝐴2 𝑉2
22 | P a g e
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
𝐴2
𝑉1 = 𝑉
𝐴1 2
Substituting, we have:
2
𝑉22 − 𝑉12 = (𝑃 − 𝑃2 )
𝜌 1
2
𝐴2 2
𝑉22 − ( 𝑉2 ) = (𝑃1 − 𝑃2 )
𝐴1 𝜌
2
𝐴2 2
𝑉22 (1 − ( ) ) = (𝑃1 − 𝑃2 )
𝐴1 𝜌
(𝑃1 − 𝑃2 )
𝑉2 =
√ 𝐴 2
𝜌(1 − (𝐴2 ) )
1
(13341.6 𝑃𝑎)
𝑉2 = √
𝑘𝑔 1 2
1.29 3 (1 − (15) )
𝑚
𝑉2 = 144.14 𝑚/𝑠
Note that the answer for velocity is above 100 m/s where we designated that
flow starts to be compressible. Therefore, our assumption for incompressible
flow is inaccurate of about 8%.
7. An experimental rocket-powered aircraft is flying with a velocity of 3000 mph
at an altitude where the ambient pressure and temperature are 151 psf and
390 ⁰R, respectively. A pitot tube is mounted in the nose of the aircraft. What
is the pressure measured by the Pitot tube?
Solution:
Is the flow subsonic or supersonic?
𝑓𝑡
88 𝑠
𝑉 = 3000 𝑚𝑝ℎ ∗ = 4400 𝑓𝑡/𝑠
60 𝑚𝑝ℎ
𝑎 = √𝑘𝑅𝑇 = √1.4 ∗ 1716 ∗ 390 = 967.95 𝑓𝑡/𝑠
𝑓𝑡
𝑉 4400 𝑠
𝑀= = = 4.55
𝑎 967.95 𝑓𝑡
𝑠
For pressure measured by the pitot tube:
𝑘
𝑃𝑂2 (𝑘 + 1)2 𝑀12 𝑘−1 1 − 𝑘 + 2𝑘𝑀 2
1
=[ 2 ] [ ]
𝑃1 4𝑘𝑀1 − 2(𝑘 − 1) 𝑘+1
1.4
𝑃𝑂2 (1.4 + 1)2 𝑀12 1.4−1 1 − 1.4 + (2 ∗ 1.4 ∗ (4.55)2 )
=[ ] [ ]
𝑃1 4 ∗ 1.4 ∗ (4.55)2 − 2(1.4 − 1) 1.4 + 1
23 | P a g e
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 2: Isentropic, Adiabatic, and other Thermodynamics for High Speed
Processes
24 | P a g e
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 03: Airflow and Compressible Flow Visualization (Prelim)
LEARNING
MODULE 03:
Airflow and Compressible
Flow Visualization
AE 324
AERONAUTICAL LABORATORY
Presented by:
AERONAUTICAL ENGINEERING DEPARTMENT
TABLE OF CONTENTS
TABLE OF FIGURES
FIGURE REFERENCE
Akbiyik, Hurrem & Yavuz, Hakan & Akansu, Yahya. (2016). Investigation
of the effect of the plasma actuators vertically placed on the spanwise of
3.1 NACA0015 airfoil. Retrieved from
https://www.researchgate.net/figure/The-schematic-view-of-the-smoke-
wire-flow-visualization-technique_fig3_301899199
Akbiyik, Hurrem & Yavuz, Hakan & Akansu, Yahya. (2016). Investigation
of the effect of the plasma actuators vertically placed on the spanwise of
3.2 NACA0015 airfoil. Retrieved from
https://www.researchgate.net/figure/Schematic-view-of-the-experimental-
setup_fig1_301899199
Retrieved from https://broadtechengineering.com/cfd-consultancy/cfd-
3.3 study-of-the-effects-of-boundary-layer-suction-on-transonic-airfoil-
performance/
3.4, 3.6 Pope A., Goin K. High Speed Wind Tunnel Testing. New York: John
3.15-3.22 Wiley & Sons.
Atish Kumar. Flow Visualization Techniques and their use. Lovely
3.5, 3.7- Professional University, India.
3.8 Retrieved from https://www.scribd.com/document/42487468/Flow-
visualization-techniques-and-their-use-RH5001A53-Atish-Kumar-CIV-208
Ueyama, Atsushi (2019). Basics of Flow. Retrieved from
3.9-3.14
https://www.cradle-cfd.com/media/column/a69
Phillip Keane, Tutorial: Performing Flow Simulation of an Aerofoil.
3.15-3.37 Retrieved from https://www.engineersrule.com/tutorial-performing-flow-
simulation-aerofoil/
TABLE OF REFERENCES
REFERENCES No.
Pope A., Goin K. High Speed Wind Tunnel Testing. New York: John Wiley &
1
Sons.
National Aeronautics and Space Administration. Wind Tunnel Index.
Retrieved October 1, 2020 from https://www.grc.nasa.gov/www/k- 2
12/airplane/shortt.html
Atish Kumar. Flow Visualization Techniques and their use. Lovely
Professional University, India.
3
Retrieved from https://www.scribd.com/document/42487468/Flow-
visualization-techniques-and-their-use-RH5001A53-Atish-Kumar-CIV-208
Ueyama, Atsushi (2019). Basics of Flow. Retrieved from https://www.cradle-
4
cfd.com/media/column/a69
Phillip Keane, Tutorial: Performing Flow Simulation of an Aerofoil. Retrieved
from https://www.engineersrule.com/tutorial-performing-flow-simulation- 5
aerofoil/
AE 324 – AERONAUTICAL LABORATORY 3|Page
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 03: Airflow and Compressible Flow Visualization (Prelim)
HONESTY CLAUSE
LEARNING OUTCOMES
The principle requirement of a smoke tunnel is for uniform flow with low
turbulence. If this is not achieved, the smoke will quickly become diffused. The
speed must not be too low, or gravity will affect the paths of the smoke
particles, curving these paths downwards. On the other hand, if the speed is
too high, the smoke filaments may become attenuated, even when large
quantities of smoke are introduced. Good lighting is needed if the smoke
filaments are to be clearly seen, and an open circuit tunnel is required in order
to get rid of the smoke. The main disadvantage is that usually only isolated
spanwise positions can be examined, and not the flow over the whole wing.
This technique involves using a heated wire inserted into the airstream
to produce a filament of air of different density from the mainstream. The path
of this filament can then be followed by using one of several optical systems,
in the context of their use with high speed wind tunnels.
Similar to the hot wire in principle is the spark technique. Small volumes
of air are heated by the discharge of a series of electric sparks. Again, an
optical system is used to trace paths of such ‘particles.’ Such systems have
been used to photograph velocity distributions in boundary layers.
For many studies of fluid flow, the main information required about the
boundary layer consists of a knowledge of the nature of the layer, whether
laminar or turbulent, the position of the transition point and the position of the
separation point, if any. Visual methods are very useful in obtaining this
information, and several such methods are common in use.
One such method involves the use of smoke. A very thin filament of
smoke is introduced very close to the surface of the model. In the laminar
boundary layer, the filament remains clearly defined. If transition occurs, the
filament is broken up suddenly, and behind the transition point there is a region
diffused smoke due to the turbulent layer. If, on the other hand, the laminar
boundary layer separates, the filament is suddenly lifted well clear of the
surface, but remains well defined for some distance downstream, leaving just
a few wisps of smoke in the dead air region. The separation and transition
points are thus easily and accurately located and distinguished from one
another. The main difficulty is to ensure that there are no disturbances at the
point where the smoke is introduced into the boundary layer. Good lighting is
essential. Again, the investigation is confined to isolated spanwise stations.
Perhaps the most widely used method for locating transition is the china
clay method. The surface of the model is thinly sprayed with kaolin (china
clay). When dry, this is a white, crystalline solid, which can be rubbed down to
a high degree of smoothness. This deposit is then sprayed with a volatile
liquid, whose refractive index is approximately the same as that of the kaolin,
so that the resultant deposit, while wet, is transparent. The model is then
placed in an airstream. In the turbulent boundary layer, the rate of evaporation
is higher than in the laminar region, so that the surface dries first behind the
transition point. When it does so, the kaolin again appears as a white, opaque
deposit, and the transition point is located. The fluids commonly used for this
purpose are ethyl and methyl salicylate. There are two special advantages of
AE 324 – AERONAUTICAL LABORATORY 7|Page
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 03: Airflow and Compressible Flow Visualization (Prelim)
this technique. After one test run, the model can be re-sprayed with fluid and
used again; the kaolin coating lasts for a considerable time. Also, the
investigation covers the whole surface of the model, and not just a few isolated
spanwise positions.
A simpler alternative to the china clay method is the liquid film method.
The model is simply sprayed with a volatile oil, and when it is placed in an
airstream, the region in which the boundary layer is turbulent dries more
quickly than the laminar region. A comparison of the wet and dry regions then
enables the transition point to be located. However, the visual contrast
between the laminar and turbulent regions is less good than with the china
clay method, and it does not last as long.
a. OPTICAL METHODS
When light passes through a region of a gas in which the density varies,
the changes in the refractive index of the gas associated with these density
changes cause the direction of the rays to be altered. The amount of the
deflection is proportional to the density gradient, and the sense of the
deflection is determined by the sign of this gradient. In a high-speed flow, the
density of the fluid changes throughout the field of flow, and the above
principle may be utilized in order to show up these variations. There are three
commonly used methods of doing so, the shadowgraph, schlieren and
interferometer systems. They are all based on the same principle, but,
whereas the first two are mainly useful simply for flow visualization, the third
may also be used for quantitative measurement, though it is much more
difficult to interpret than the other two. The shadowgraph system is the
simplest, but the schlieren system is more effective, and is the most commonly
used of all three methods.
Other systems may employ optical means, e.g., prisms, of splitting light
from a white source into different colors. Color schlieren systems give
attractive and easily interpreted pictures, but the definition and contrast are
usually not so good, nor so easily controlled, as in black and white.
• The Interferometer
For a flow which does not change with time, the streamline, streakline,
and pathline are the same line. A flow which does not change with time is
called a steady-state flow . On the other hand, a flow which varies with time is
called a transient flow. For transient flows, the streamline, streakline, and
pathline are all different lines.
A streakline represents all the points that have passed through a certain
location. For this situation, where the streakline can be likened to represent
the smoke trail from a chimney, the streakline goes to the south during the first
10 seconds because all the smoke is going south. Then, when the wind shifts
to the east, all the smoke particles that were initially heading south (emitted
for Time < 10 seconds) start to be offset to the east. The newer smoke
particles (emitted for Time > 10 seconds) head directly east. After 20 seconds
the streakline is at a right angle as shown on the right in Figure 3.13.
A pathline goes to the south during the first 10 seconds just as the
streakline. The pathline can be thought of as the path traced by a balloon
floating in the air. When the wind shifts to the east, the balloon starts moving
east. The pathline goes to the east from the point where the wind direction has
been changed. As a result, after 20 seconds the pathline bends at a right angle
as shown on the right in Figure 3.14.
As can be seen in this chapter, the analysis results for transient flows
will be different depending on the method used to visualize the flow.
Understanding this difference is important when visualizing analysis and
experimental results.
tube is also crucial. The flow speeds up round the nose, and if the holes are
too far forward, we may get reduced pressure readings as a result. The holes
are therefore placed some distance back from the nose, where the effect of
the stem or mounting helps to compensate for this effect.
As the Mach number increases, the local accelerations in the flow may
give rise to regions of supersonic flow with shock waves and Consequent
errors in reading. To minimize these effects, we need to use a very thin and
fairly pointed tube. At supersonic speeds, it is necessary to have a sharp nose
to ensure shock attachment, but, in the absence of special design, there will
still be errors resulting from the presence of the shock. There are various
possible designs for a static tube for use at supersonic speeds. These include:
(a) A very thin, single-sided wedge, such as that shown in Fig. 3.7. In
the design attitude, there is no flow deviation on to the lower surface, and the
pressure recorded there is the free stream pressure. Such a device is very
sensitive to incidence changes, because they give rise to shock or expansion
waves whose intensity does vary with incidence.
(b) A long, thin cone, with four holes located well behind the shoulder,
as shown in Fig. 3.8. This device records values of the pressure which are
practically the same as the free stream values, because the expansion round
the shoulder is equal and opposite to the compression at the nose. It is less
sensitive to incidence changes than the wedge described above.
Like the airspeed, the Mach number can be determined, either in wind
tunnel experiments or in an aircraft in flight, from pitot and static tube readings.
A Machmeter uses the pressures taken from a pitot and a static tube to
record the ratio of total pressure minus static pressure to static pressure, and
this ratio determines the Mach number uniquely. In supersonic flow, however,
the total head is necessarily measured behind a normal shock, so that its value
is less than the free stream value. Nevertheless, for given conditions upstream
of the shock, the total head behind it is determined, and, if 𝑝𝑇 , is the total head
behind the shock, i.e., the total head recorded by the pitot tube, and H, is the
free stream total head, then we may write 𝑝𝑇2 /𝑝𝑇1 , as a function of Mach
number. Also, 𝑝𝑇1 /p1, is a function of Mach number, where p1, is the free
stream static pressure, so that we may write 𝑝𝑇2 /p1, as a function of Mach
number. Therefore,
𝑝𝑇2 − 𝑝1
= 𝑓2 (𝑀)
𝑝1
Thus, the same pressure ratio as in the case of subsonic flow again
determines the Mach number uniquely, though the functional relationship is
not the same.
• MEASUREMENT OF AIRSPEED
In the case of supersonic flow, the problems are a little more difficult.
The same airspeed indicator is used, which consists of pitot and static tubes
which now ideally record the pressure difference 𝑝𝑇 -P, where 𝑝𝑇 , is the total
head behind the shock, and p, the static pressure ahead of it, and the airspeed
indicator reading, as always, is based on this pressure difference. However,
we have seen above that (𝑝𝑇𝑎 – p)/p, is a function of Mach number, so that
(𝑝𝑇 - P) is a function of Mach number and p. Thus, a correction for Mach
number will have to be applied to the airspeed indicator reading, so that the
instrument will always need to be read in con junction with the Machmeter.
Further, the correction will vary with altitude, since p, depends on altitude.
b. YAWMETERS
Figure 3.19 depicts two examples of bent tube yawmeters. The arms
are pitot tubes inclined at 45̊ to the axis of the instrument, since pitot tubes are
very sensitive to yaw at an angle of 45 ̊. When the pressure difference between
the two arms is zero, the flow is parallel to the axis. As the inclination of the
flow to the axis of the yawmeter increases, so also does the pressure
difference between the two arms.
AE 324 – AERONAUTICAL LABORATORY 22 | P a g e
PHILIPPINE STATE COLLEGE OF AERONAUTICS
INSTITUTE OF ENGINEERING AND TECHNOLOGY
AERONAUTICAL ENGINEERING DEPARTMENT
Learning Module 03: Airflow and Compressible Flow Visualization (Prelim)
V. SPECIAL FACILITIES
• FLIGHT TESTS
Assessment Activity
Answer the following questions in your notebook.
1. Differentiate the different flow patterns in visualizing flow.
2. Enumerate the various systems used in flow visualization. Explain how each of
the operates.
3. Explain the importance of visualizing the flow in aerodynamics; in aviation.
4. Classify and analyze the measurement of fluid velocity, total head, static
pressure, Mach number and airspeed
5. In your own words, how do we visualize flow?
I. AIRFOIL MODELLING
4) Click OK, and you will see the airfoil curve appear in the design window, as
seen in the image below.
5) Of course, being a curve, it is still not useful for creating geometry, so select
the Front Plane from the design tree and click Sketch from the Sketch tab.
6) Now click Convert Entities from the Sketch tab, and in the Convert Entities
panel, select the airfoil curve from the design window.
7) Next, we want to make a centerline from the trailing edge to just beneath the
leading edge. This will represent the chord length of the airfoil, and once we
have constrained it we can alter the chord length at will.
8) After the chord line is sketched, we need to put another line connected to the
last line near the leading edge. This new line needs to be tangential to the
airfoil, as shown below.
Figure 3.20
9) Then, we can select both the chord line and the tangent line and constrain them
so that they are perpendicular to each other. Why? Because when we rotate
the sketch or extend the chord length, we want it to retain shape, and the
perpendicular constraint will ensure that the whole thing remains airfoil-shaped.
10) Now that the sketch is constrained, we can just double-click on the chord
line and enter a value for how long we want the chord to be. In this case, let’s
set it to 1.6 meters.
Figure 3.21
11) Congrats! You have now converted your airfoil curve into a sketch entity.
Now we can model our solid airfoil.
c. Airfoil 2D to 3D
1) This part is easy. Simply select the airfoil sketch and extrude it to 4 meters.
This will provide us with a basic constant-chord (i.e., non-tapered),
rectangular wing. This type of wing, incidentally, is referred to colloquially as
a “Hershey Bar”.
Figure 3.22
2) And there it is. Our Hershey Bar wing is now ready for some flow simulation!
1) Load up the Flow Simulation add-in by clicking Tools > Add-ins and checking
the SOLIDWORKS Flow Simulation box. Once it is loaded, select the Flow
Simulation tab and click the Wizard button to start the Flow Simulation
Wizard.
Figure 3.23
2) On the first page of the wizard (Project Name), name your project and click
Next.
3) On the second page (Unit System), select your preferred unit system. For
consistency, we will select SI units here (m-kg-s). Then click Next.
4) On page three (Analysis Type), we can select Internal or External study.
Internal studies are for simulating flows that are constrained by some kind of
vessel, such as a pipe, and external studies are for simulating flows around
the outside of a body such as a truck or an airfoil. So, we click External, and
then press Next to advance to the next page.
5) The next page (Default Fluid) allows us to select the fluid in our study. This
is an aerodynamic study, so we select Air from the top list and click Add.
Once the default fluid has been added, we can click Next.
Figure 3.24
6) We can skip over the next page (Wall Conditions) by clicking Next.
7) The final page that we need to deal with in the wizard is the Initial and
Ambient Conditions page. This is where we set the temperature and
pressure of the environment and the velocity of the flow in the x-direction.
We have set the temperature and pressure to SSL (standard sea level)
values and the velocity in x-direction to 55m/s (about 200km/h).
Figure 3.25
8) That’s all we need to worry about with the wizard. Click Finish and the wizard
will close.
Figure 3.26
9) You will notice that the wizard has created a box around the wing. This is our
Computational Domain, where all the magic happens. Think of it as the inside
of a wind tunnel. Everything inside it is part of the simulation, and everything
outside it is irrelevant.
10) Note that a larger Computational Domain requires more processing.
11) Click on Computational Domain on the left-hand panel (as seen below) and
you will notice six handles appear on the box. Drag these handles until the
domain box fits just around the wing model. Be sure to leave enough room at
the fore and aft of the wing so we can get some sweet visualization of the fluid
flow as it passes around the wing.
Figure 3.27
1) First, we want to select the faces of the wing that we want included in the
study. In the Surface Goals panel, click the blue Selection area to activate it and
click all of the faces of the wing model.
Figure 3.28
2) Next, go down the list and check the minimum, maximum and average for
the following parameters:
Note that we have selected Velocity (X), because this is the direction
that the flow will be travelling in.
Figure 3.29
5) Check the iterations box and ensure it is set to 100 iterations. It may be that
your simulation requires less, or even more. But for now, 100 iterations is
fine. This should be enough for the goals to reach convergence. More
iterations will generally give a better result, but after a point, the trade-off
between accuracy and time-taken simply isn’t worth it. You can run the
simulation all day long and the gains to accuracy will become very modest.
So, 100 is fine in this case. Click OK to exit.
6) Now that our simulation is set up, we can run it. You can find the Run button
in the top ribbon (as seen below). Click it and you will see the solver screen
appear, informing you of how many iterations are left.
Figure 3.30
1) Now the calculations have finished, we can go into the study panel on the
left and expand the Results section to show us a selection of graphs and
plots. Right clicking any of these plots will allow you to insert the plot into the
main window.
Figure 3.31
d. Cut Plots
1) The first plots we will look at are cut plots. This type of plot will display a 2D
slice (a plane) of the model, and you can drag the green arrow to move the
slice along any part of the 3D model.
2) Right click Cut Plot and select Insert.
3) In this instance, I select Front Plane, then I select Contours to show a
contour plot. In the Contours section, you can see that the default parameter
should be Velocity (X). We would like to see the pressure contours here, so
we can click the parameters box and select Pressure.
Figure 3.32
4) Click the green check mark and you will see your plot appear in the main
design window. You can move the slice along the length of the wing by using
the green drag handle and you can rotate the plot as you would do your 3D
model. The image below shows an isometric view and a side view. The color
code shows how the colors relate to differences in pressure.
Figure 3.33
5) Because the results are already loaded into your computer, you can easily
switch between data types by clicking the parameter name just beneath the
colour scale and selecting new results to display.
Figure 3.34
6) So, if I want to change from a pressure contour plot to a velocity contour plot,
I simply click Pressure beneath the colored scale (as seen above) and switch
it to velocity. The main plot will change accordingly.Note, if you want to see
the slice scan along the entire length of the wing, you can right click on Cut
Plot and select Play for a little animation.
Figure 3.35
e. Flow Trajectories
1) Cut plots are nice, but they don’t show the holistic view of what is going on;
they simply show a 2D slice of the 3D whole. The trajectory plot is more
useful for showing behavior over the full length of the wing at any given time.
This is more like a wind tunnel with smoke injected into the chamber, which
you may be more familiar with from university.
2) Right click on the Flow Trajectory option in the study pane and select Insert.
This will open the Flow Trajectories pane.
3) In this pane, select the faces that we want to be a part of the study as we did
for the Cut Plot.
4) In the Number of Points box, type 15 and set the Spacing to 0.03m.
5) In the Appearance section, we select Static option, and then we select the
appearance of the trajectory. In this instance we select Pipes, but feel free to
play around here and experiment with different appearances.
Figure 3.36
6) Again, in this plot we will be looking at Pressure, so select that from the
Appearance section, and then click the green check icon. The plot will appear
in the main window, as you can see below.
Figure 3.37
7) Here we can see the variations of pressure as the air flows over the aerofoil,
and also we can get some idea of the turbulence/vortices created by the wing
tip.