F24AeroII cn01
F24AeroII cn01
Aerospace vehicles & flight missions (Flight Mechanics & Control, Aerodynamics, Solid Mechanics, Propulsion)
Design and Analysis
Objective of aerodynamics and the basic problem (including ground coordinate system and body C.S.)
WHAT/WHY/HOW of aerodynamics:
WHAT: (vehicle as a solid boundary) Physical objectThe Concept of continuum
Fluid element and flow field.
WHY: the science of aerodynamics (Physical variables, governing equations, see Appendix A below).
HOW: apply the analytical knowledge, design aerospace vehicles.
Aerodynamic Models:
Bernoulli Equation (Chapter 3, Anderson)
Low-Speed Thin-Airfoil Theory (Chapters 3 and 4)
Lifting-Line Theory (Chapter 5)
Comments:
(1) Presentation of core technical content: (i) verbal reasoning; (ii) formula—equations, derivations,
computations.
(2) Degree of difficulty: high (seems abstract; easily get confused).
F24AeroII—CN#01: Page 1 of 23
(3) Warning: materials covered in weeks 3 & 4 are critical.
(4) Suggestion: AERODYNAMIC MODELS are the key.
(5) Philosophy: WHAT/WHY/HOW, and we shall cover:
(6) The object of interest and (i) its motion in fluids, or (ii) fluids going through,
(7) Applying fundamental governing equations for fluid motion in specific cases.
(8) Design the objects that meet certain requirements.
We are mainly interested in the changes of flow in the main flow direction, x. The cross-sectional area is
↘ ↗ ↘ Nozzle
↗ ↘ ↗ Diffuser
A typical problem—Given: ρ, A1, V1, p1; Given: A2; Find: V2, p2;
Diffuser with given ρ, A1, V1, p1: as A2 goes very large, V2→0; p2→p0; Kinetic energy → Internal energy
Nozzle with given ρ, A1, V1, p1: as A2 goes very small, V2→∞; p2→-∞; Kinetic energy ← Internal energy
One may argue: p2 can only go as low as 0. As p2→0, can density stay as a constant?
F24AeroII—CN#01: Page 2 of 23
Volta Conference (1935): see Section 7.1, Anderson
Discussion: airplane at 800 km/h (Air speed record in 1934: 700km/h, see Appendix B below)
Myth: sound barrier; Swept wing design (A. Busemann)
Importance of variable-density aerodynamics; Transonic and supersonic wind tunnels
High speed High kinetic energy Heating Temperature increase
The field of compressible flow theory (subsonic, transonic, supersonic, and hypersonic)
We may follow the thinking of the leading scientists and engineers and ask:
Question 0101: Why do we need to consider changes in density when flight speed is high?
Question 0102: What does “high” speed mean quantitatively?
Question 0103: What exactly causes density to change?
Example: Messerschmitt Me 262 – German Luftwaffe (the world's first operational jet-powered fighter
aircraft.) (http://en.wikipedia.org/wiki/Messerschmitt_Me_262)
F24AeroII—CN#01: Page 3 of 23
Example: Lockheed Martin F35 Lightning II
(http://en.wikipedia.org/wiki/Lockheed_Martin_F-35_Lightning_II)
Bell X-1
F24AeroII—CN#01: Page 4 of 23
2. Flight Speed and Mach Number: An Example
Before we reveal the answers to the questioned listed under the section on the Volta Conference, we first
take a look at an example we have studied before.
Question: An aircraft flies at a constant speed V in the air. Take the air as the standard sea-level atmosphere.
Convert the problem to a coordinate system fixed on the aircraft (why do we want to convert to such a
coordinate system?). Find the Mach number of the incoming free stream at various values of V.
Solution: Standard sea-level air has the following thermodynamic property:
Temperature: T=288K; Pressure: p=1.0125x105Pa (density ρ=1.225kg/m3); Adiabatic coefficient: γ=1.4;
Gas constant for air: R=287.14
We may find the sonic speed (to be explained and derived later) to be,
𝑐= 𝛾𝑅𝑇 = 340𝑚/𝑠
Case 1 2 3 4 5 6 7 8
V (km/h) 10 100 500 1000 1,224 2,500 6,000 12,000
V (m/s) 2.78 27.8 139 278 340 694 1,667 3,333
M=V/c 0.008 0.08 0.41 0.82 1 2.04 4.9 9.8
Comments:
For cases 1 and 2, we can use the low-speed theory covered in Aerodynamics I to solve for the lift
coefficient of the airplane (or to be exact, airfoil sections and the wing). Case 3 is already a stretch for us
regarding the low-speed aerodynamic models. Cases 4 to 8 are definitely beyond our present knowledge.
The Mach numbers are presented in the table above. But we have not dealt with it in a rigorous way. It
is just simply stated that the Mach number of 0.3 is a dividing point for low- and high-speed flows. Naturally,
we need to find out what the sonic speed really is, and why shall we compare flow speed with the sonic
speed, before we can understand compressible flows.
We shall come back to solving for the high-speed flows around an airfoil (or a wing) in Chapters 11 and
12. You may find out that the materials we cover from now on until these chapters are mostly for the so-called
quasi-one-dimensional (Q1D) flows. These simple flows, along with some 2-D cases, would actually provide
us with almost all the essence of inviscid compressible flows, which we shall appreciate when the time comes
for us to review the entire course.
where v is the specific volume, and p the pressure. Thus, from definition, τ is the relative change in specific
volume per unit change in pressure. Or, for 1 Pa increase in pressure, how much change (decrease) in the
volume occupied by 1kg mass of fluid has occurred. For positive dp (compression), dv is negative, so we have
F24AeroII—CN#01: Page 5 of 23
the minus sign in the definition to make τ a positive parameter. We may replace the expression with the density,
as it is far more common to use density in aerodynamics rather than specific volume,
𝜏≡ (1a)
Of course, τ is a physical parameter of the fluid. We may measure it under various conditions/processes.
Generally, the value of τ depends on the precise process of the compression, not just the starting and ending
thermodynamic states. We can further define the isothermal compressibility and the isentropic compressibility
respectively as,
𝜏 ≡ 𝜏 ≡ (1b)
which are more precise in thermodynamic sense. (Here, we have dealt with the compressibility of a fluid in a
static way. What happens to it when the fluid is in motion, which is the case we really care about in
aerodynamics?)
4. The Sonic Speed (Disturbance Generated by A Plane Source. Section 8.3, Anderson)
Fig. 1. Left: a fire. Center: a beep of disturbance. Right: a beep sending out disturbance waves.
An explosion: chemical energy is released as heat into the air, which gains internal energy.
A ‘beep’: a very small explosion which releases a tiny amount of energy into the air.
Molecular point of view: outwards from the center of the beep, energy is transferred to the air (say, originally
standard sea-level, and stationary). One may reasonably argue: the process happens layer-by-layer (spherical),
reaching larger and larger radius. If the process happens at a finite speed, we may view it as a spherical wave.
Let’s look at an example: for stationary air everywhere with constant thermodynamic state of p, ρ, T,
suppose there is a disturbance at point A where a tiny solid ball suddenly explodes. From a physics point of
view, the solid at A will have interaction (pushing the air outward) with the air molecules around it.
F24AeroII—CN#01: Page 6 of 23
Fig. 2. A plane disturbance of infinitesimal strength.
From a fluid mechanics point of view, there are pressure waves propagating from point A. The air is been
pushed, its pressure would increase, albeit by a very small amount. Similar changes happen to ρ and T.
From fluid mechanics theory, we can find out the speed at which the wave travels!
Considering the symmetry of the case, the pressure wave is of spherical shape, which would make the
solution process a little bit more complex than the case shown below and it is left as a homework assignment.
We shall proceed instead to an easier case next.
Consider fluid in a duct with constant area A: V=0, p, ρ, T. A very small wave passes through. The wave
speed is at a constant c (For now, let us just assume this speed is a constant).
After the wave passes, parameters changed to p+dp, ρ+dρ, T+dT, c-dV. Let the viewer travels with the
wave, i.e., in a different frame of reference (F.O.R). Now, we have three unknowns: dV, dp, dρ.
𝑐= (2)
So, somehow, the change in density caused by change in pressure is quantified by a velocity measure
(m/s). We can see from this, at least vaguely, that flow speed is related to the compressibility of a fluid. The
compressibility becomes significant when the flow reaches high enough speed. The isentropic compressibility
can be expressed as,
𝜏 ≡ = 𝑐= (2a)
Definition: Speed of sound (sonic speed) is the propagation speed of small perturbations in a medium.
For a perfect gas such as air, the isentropic relation gives,
𝑝 𝑑𝑝 𝑑𝜌 𝑑𝑝 𝛾𝑝
= 𝑐𝑜𝑛𝑠𝑡. => =𝛾 => = => 𝑐 = 𝛾𝑅𝑇
𝜌 𝑝 𝜌 𝑑𝜌 𝜌
𝜏 = = = (2b)
F24AeroII—CN#01: Page 7 of 23
on the temperature for a specific perfect gas.
(2) The sound speed, even with the word speed in its name, is actually a thermodynamic parameter; it does
not depend on the flow velocity! It is not natural at all to compare the sonic speed and the flow velocity.
(3) Flow speed can be of an arbitrary value in an arbitrary coordinate system. Talking about an arbitrary
relative value of flow speed is nonsense. Comparing the flow speed to the sound speed should be in the sense
that the flow speed is relative to the source of the disturbance. A high flow speed should mean that the flow
speed is large relative to the source of disturbance (such as a solid wall, which forces the flow to change
direction).
(4) For Air at 15oC: c is about 340m/s=1,224 km/h (Note: see the Standard Atmosphere Model below); For
liquids: a typical value of c is about 1500m/s=5,400 km/h. For incompressible flows, c is infinity, meaning any
perturbation would be felt throughout the flow field immediately.
(5) For a general flow field, we may view the thermodynamic parameters as fields of space-time. Thus, Eq.
2 can also be viewed as a field, i.e., c is a thermodynamic field.
Case A: the disturbance source is moving while generating a wave at time intervals of 1 sec.
∆t=1s, c=300m/s, V=200m/s. For t=0, 1, 2, 3, and 4s.
Observations: (1) A wave would contain all the later-generated waves. (2) The disturbance source is contained
inside all the waves. (3) As enough time passes by, any faraway place would feel the disturbance.
F24AeroII—CN#01: Page 8 of 23
Case B: conditions the same as in Case A, except that V=600m/s
Observations: (1) Waves can cross each other (not all). (2) The disturbance source is ‘leading’ all the waves.
(3) There are always places ahead of the source that do not feel the disturbance. (4) All the waves are
contained in a0 wedge/cone (for 2D/3D, respectively).
The results show that the pictures differ qualitatively for cases with V being less than, equal to, or larger
than c, respectively. The patterns of the wave family would certainly be different from each other. When V is
less than c, a larger wave (larger radius, generated less recently) contains a smaller wave (smaller radius,
generated more recently). When V is the same as c, a larger wave contains a smaller wave, but they share a
common point, which is also the point the source currently resides. When V is the larger than c, a larger wave
may not contain a smaller wave and the source is always to the left of all waves. Also, there exists a wedge
(cone for the case with a point source) with its tip at the source and contains all the waves. More precisely, the
wedge is tangent to all the waves, which certainly seems more interesting. The wedge/cone seems to be
special thing. And it has a name: the Mach cone (or Mach wedge). Of course, the shape of the wedge/cone
can be represented entirely by its half-angle which satisfies:
𝒄 𝟏 𝟏
𝒔𝒊𝒏 𝝁 = 𝑽 = 𝑽/𝒄 ≡ 𝑴 (3)
We call this angle μ the Mach angle. For it to be meaningful, M=V/c must be no less than 1.
F24AeroII—CN#01: Page 9 of 23
Mach Number: M=V/c
M < 1: flow is subsonic; M > 1: flow is supersonic; M < 1 and > 1 co-exist: transonic.
What’s the meaning and significance of M? M represents the ratio of macro kinetic energy to internal energy;
so when M is small, V has small influence on the thermodynamic state of the fluid. M also represents the
compressibility of the fluid in motion. M is a very important parameter for compressible flows
F24AeroII—CN#01: Page 10 of 23
(4) For aeronautical engineers, the most important parameter here is the altitude. Common sense tells us
the higher the altitude, the cooler the air is and the lower the air pressure is.
(5) People have devised a term called the standard atmosphere, which gives the air properties for altitude
from zero (sea-level) up to the outer edge. It is based on empirical data and is season-averaged.
(6) Now, we usually use the U.S. Standard Atmosphere (1976) as our source for getting the air properties.
Please refer to “US_Standard_Air1976.m” and “US_Standard_Air1976_test.m” for getting the air properties.
Below is a plot of the pressure, density, and temperature of the standard air at altitudes between 0 and up to
80km. You will need to use these Matlab files in your future study, so learn how they work.
4 U.S. Standard Atmosphere, 1976
x 10
9
p/p0
8 rho/rho
0
T/T
7 0
6
Altitude, m
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
Short Review
(1) Pursuing higher speed
(2) Density change (compressible) is significant with high-speed flight
(3) “High-speed” is relative to “sonic speed,” which is the propagation speed of small disturbances (pressure
wave)Mach number & Mach Cone“Regions of Dependence” and “Regions of Influence”
(4) Density change caused by pressure change is the main concern (essence of compressible flows)
Looking Ahead……
(1) Focus on steady, inviscid high-speed flows
(2) 1D2D3D; InternalExternal;
(3) Aerodynamic Models (Q1D Isentropic Model, Normal Shock Model, Shock-Expansion Model)
F24AeroII—CN#01: Page 11 of 23
Appendix A: Review of the Fundamental Governing Equations
Fig. A1. A portrait of the great Leonhard Euler,1707-1783 (pic. taken from Wikipedia).
Euler is accredited with applying the Newton mechanics to fluid flows and building the so-called Euler
equations. His approach (the Eulerian Approach) of constructing an arbitrary control volume that is fixed in
space in the chosen coordinate system, contrary to a more direct/natural approach where a fluid element or a
group of fluid elements are followed and tracked all the time, was a work of genius. He should be regarded as
the greatest founder of fluid mechanics/aerodynamics as well as the most important scientist in the field of
analytical mechanics. Note: to many learners of fluid mechanics, failing to switch to the Eulerian Approach
frame of mind is the main cause of a lot of confusion and misunderstandings. It is highly recommended that a
student should have a quite thorough grasp of the Eulerian Approach before learning the Lagrangian Approach.
F24AeroII—CN#01: Page 12 of 23
Can we take the multi-body approach? It turns out that we should take fluid as a continuum and concern
with the macroscopic properties of fluid
Example: define the density field. We will not use the approach of analyzing a “fluid particle”, but rather
a small sum of fluid particles – “fluid element”
The important forces are applying on surfaces, not on specific particles. This is the main difference
between the study of classic mechanics and continuum mechanics
Summary:
1. The smallest mass of interest is called a fluid element, which is infinitesimally smalla point in space
2. The physical properties of a fluid element are macroscopic properties of many fluid molecules
3. Any portion of fluid consists of continuously distributed fluid elements
Thinking Process:
Interaction between fluid elementsForcesMotionNewton’s Second Law
F24AeroII—CN#01: Page 13 of 23
The power of differentiation!
We define aerodynamic variables at physical points. BUT, we can not just look at points
1. ACTUALLY, we must look at differentiation (infinitesimal)
2. Of course, we will also do the opposite operation—integrate
3. Finally, we will have things in differential form that are easy to use
𝑚(𝑡) = 𝜌(𝑥, 𝑦, 𝑧, 𝑡) 𝑑∀
∀
Continuity Equation
a. Control Volume and Control Surface
An arbitrary control volume is shown in Fig. A2. The volume is enclosed by a closed surface, S, which is
called the control surface. In the Eulerian approach, we let it have a fixed location in space. As time goes by,
fluid comes in/out of the CS.
Control
Volume, V
Control Surface, S
F24AeroII—CN#01: Page 14 of 23
water
n
Vn
θ
V (x1,y1,z1,t)
point A
(x1,y1,z1)
dS Vt
Figure A4: A small piece of the control surface with a velocity vector
As shown in Fig. A4, on dS, we can define a local velocity vector, V, of the fluid there. Now, we are
interested in the amount of mass flow out/in of the CV through dS. Geometry dictates that,
𝑉 = 𝑉⃗ • 𝑛⃗ = 𝑉 𝑐𝑜𝑠 𝜃 𝑉 = 𝑉 𝑠𝑖𝑛 𝜃
The parallel component, Vt, does not bring mass in/out of the CV through dS. It just flows by dS. The
mass that flows in/out of the CV through dS from time t to t+Δt is,
𝑀𝑎𝑠𝑠 = 𝑑𝑒𝑛𝑠𝑖𝑡𝑦 × 𝑣𝑜𝑙𝑢𝑚𝑒 = 𝑑𝑒𝑛𝑠𝑖𝑡𝑦 × 𝑎𝑟𝑒𝑎 × 𝑙𝑒𝑛𝑔𝑡ℎ
𝛥𝑀 = −𝜌(𝑥 , 𝑦 , 𝑧 , 𝑡)(𝑑𝑆)(𝑉 𝛥𝑡)
where the minus sign is due to the convention that the mass flows into the CV is marked as positive and that
flows out is regarded as negative. The mass flow through dS per unit time is,
𝛥𝑀
𝑑𝑚̇ = = −𝜌(𝑥 , 𝑦 , 𝑧 , 𝑡)𝑉 (𝑥 , 𝑦 , 𝑧 , 𝑡)𝑑𝑆 = −𝜌𝑉⃗ • 𝑛⃗𝑑𝑆 = −𝜌𝑉⃗ • 𝑑𝑆⃗
𝛥𝑡
where we use lower case m(dot) for mass flow rate (kg/s)
c. Mass Flow through the Entire Control Surface
The total mass flow rate through the entire CS is just the summation of all the little pieces of dS,
𝑚̇ = − ∑ 𝜌𝑉⃗ • 𝑑𝑆⃗ = − ∑ 𝜌𝑉 𝑑𝑆 (A1)
which is just the following integration in the limit of each dS having infinitesimal area,
𝑚̇ = 𝑙𝑖𝑚(𝑑𝑆 → 0) ∑ 𝜌𝑉⃗ • 𝑑𝑆⃗ = ∯ 𝜌𝑉⃗ • 𝑑𝑆⃗ (A2)
F24AeroII—CN#01: Page 15 of 23
d. Conservation of Mass
Now go back to building the continuity equation. Conservation of mass simply states: mass cannot be
destroyed nor created. Applied to the situation shown in Fig. A2, we have,
(1) The time rate of change of the total mass inside the CV
plus (2) the mass flow rate (going out of the CV) through the CS
must equal zero.
We have calculated the latter value (Eq. A2). Now, let us calculate the first part.
e. Mass inside an Infinitesimal Fluid Element and a Control Volume
Consider first an infinitesimal fluid element around a point B at (x2, y2, z2) inside the CV, the mass of
which is,
𝑀𝑎𝑠𝑠 = 𝑑𝑒𝑛𝑠𝑖𝑡𝑦 × 𝑣𝑜𝑙𝑢𝑚𝑒 𝛥𝑀(𝐹. 𝐸. ) = −𝜌(𝑥 , 𝑦 , 𝑧 , 𝑡)𝑑∀
The total mass inside the CV is the limit of the summation of all the small elements when d(Vol)0 in
size.
The time rate of changes in the total mass inside the CV is,
( ) ( , , , )
= ∰∀ 𝜌(𝑥 , 𝑦 , 𝑧 , 𝑡)𝑑∀ = ∰∀ 𝑑∀ (A3)
Note that the CV does not change with time, so we can take the partial derivative of t inside the integral.
Eq. A4 is the integral form of the continuity equation. We have not assumed anything about the shape, the size,
nor the location of the control volume. The only thing that is required is that the CS is a closed surface. Thus,
Eq. A4 is valid for all shapes and sizes of CV in the flow field as one wish to define. Note: from a strict calculus
point of view, there are requirements placed on the CS and the differential properties of the physical variables.
Among them are continuity of the variables, differentiability, etc. Please check relevant books if you are
interested in these. For now, we will just assume that Eq. A4 can be written and any additional requirements
are implicitly satisfied.
Eq. A4 can be put into a single integration via the Gauss Theorem. The Gauss Theorem says,
𝐴⃗ • 𝑑𝑆⃗ = 𝛻 • 𝐴⃗ 𝑑∀
∀
which means that, for any vector field A, its (directional) integral over a closed surface, S, equals the volume
integral of the divergence of A on the volume enclosed by S. Replace the vector with ρV (which is indeed a
vector field), we get the second term in (A4) as,
F24AeroII—CN#01: Page 16 of 23
𝜕𝜌
𝑑∀ + 𝛻 • 𝜌𝑉⃗ 𝑑∀ = 0
∀ 𝜕𝑡 ∀
∰∀ + 𝛻 • 𝜌𝑉⃗ 𝑑∀ = 0 (A5)
which is the differential form of the continuity equation. In Cartesian coordinate system, we have,
+ + + =0 (A7)
where, u, v, and w are the x, y, and z component of the velocity vector, respectively.
Eqs. A6 and A7 can be rewritten as,
𝜕𝜌 𝜕𝜌 𝜕𝜌 𝜕𝜌 𝜕𝜌 𝜕𝑢 𝜕𝑢 𝜕𝑢
+ 𝑉⃗ • 𝛻𝜌 + 𝜌𝛻 • 𝑉⃗ = 0 +𝑢 +𝑣 +𝑤 +𝜌 + + =0
𝜕𝑡 𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧
We shall see the power of the deferential form later on, although, it should be critically noted that the
integral form should always be regarded as the more fundamental form.
𝛻 • 𝜌𝑉⃗ 𝑑∀ = 0 → 𝛻 • 𝜌𝑉⃗ = 0 → 𝑉⃗ • 𝛻𝜌 + 𝜌𝛻 • 𝑉⃗ = 0
∀
𝜕𝜌𝑢 𝜕𝜌𝑣 𝜕𝜌𝑤 𝜕𝜌 𝜕𝜌 𝜕𝜌 𝜕𝑢 𝜕𝑢 𝜕𝑢
+ + =0 𝑢 +𝑣 +𝑤 +𝜌 + + =0
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧
Definition: A flow is regarded as incompressible if the density of all fluid elements does not change with
time. Incompressible does not necessarily mean ρ(x, y, z, t)=constant. We can have two elements with different
density values, which do not change with time. We shall look at this later on. For uniform incompressible fluid
flow, ρ(x, y, z, t)=constant, we get,
𝜕𝑢 𝜕𝑣 𝜕𝑤
𝛻 • 𝑉⃗𝑑∀ = 0 𝛻 • 𝑉⃗ = 0 + + =0
∀ 𝜕𝑥 𝜕𝑦 𝜕𝑧
which is true for steady and unsteady flows.
F24AeroII—CN#01: Page 17 of 23
Momentum Equation
Referring to Fig. A2 again, Newton’s second law dictates that,
(1) The time rate of change in momentum for a CV
plus (2) the momentum flow rate (going out of the CV) through the CS
equal (3) the total force acting on the CV.
Control
Volume, V
Control Surface, S
The time rate of changes in the total momentum inside the CV is,
𝜕𝑀 (𝐶𝑉) 𝜕 𝜕𝜌𝑉⃗
= 𝜌𝑉⃗𝑑∀ = 𝑑∀
𝜕𝑡 𝜕𝑡 ∀ ∀ 𝜕𝑡
Note again, the CV does not change with time, so we can take the partial derivative inside the integral.
Through the CS,
𝑀𝑜𝑚𝑒𝑛𝑡𝑢𝑚 𝑓𝑙𝑜𝑤 𝑟𝑎𝑡𝑒 = 𝑚𝑎𝑠𝑠 𝑓𝑙𝑜𝑤 𝑟𝑎𝑡𝑒 × 𝑣𝑒𝑙𝑜𝑐𝑖𝑡𝑦
𝑀𝑎𝑠𝑠 = 𝑑𝑒𝑛𝑠𝑖𝑡𝑦 × 𝑣𝑜𝑙𝑢𝑚𝑒 = 𝑑𝑒𝑛𝑠𝑖𝑡𝑦 × 𝑎𝑟𝑒𝑎 × 𝑙𝑒𝑛𝑔𝑡ℎ 𝛥𝑀 = −𝜌(𝑑𝑆)(𝑉 𝛥𝑡)
𝛥𝑀 (𝛥𝑀)𝑉⃗
= = −𝜌𝑉 𝑑𝑆𝑉⃗ 𝑀 (𝐶𝑆) = − 𝜌𝑉⃗ 𝑉 𝑑𝑆 = − 𝜌𝑉⃗ 𝑉⃗ • 𝑑𝑆⃗
𝛥𝑡 𝛥𝑡
The forces acting on the CV can be divided into two groups. First, there are forces that act on each of
the fluid element inside CV without making direct contact. We call these “body forces”. For example, gravity is
a body force. We use F for body force per unit mass.
𝐵𝑜𝑑𝑦 𝐹𝑜𝑟𝑐𝑒 = 𝑚𝑎𝑠𝑠 × 𝑓𝑜𝑟𝑐𝑒 𝑝𝑒𝑟 𝑢𝑛𝑖𝑡 𝑚𝑎𝑠𝑠 𝛥𝐹⃗ = (𝜌𝑑∀)𝐹⃗
Secondly, there are forces acting on the CS. We have talked about the pressure and shear forces. These
forces are molecular in nature. They act on the fluid elements on the CS. The source of the forces is from
neighboring fluid elements or solid body. Here, let’s consider only the pressure forces. By neglecting the shear
forces, i.e., neglecting the influence of viscosity, we confine out study to inviscid flows. As shown below, the
pressure force acting on a small dS is always in –n direction.
F24AeroII—CN#01: Page 18 of 23
n
p (x1,y1,z1,t)
point A
(x1,y1,z1)
dS
Figure A5: A small piece of the control surface with a pressure force
𝐹𝑜𝑟𝑐𝑒 = 𝑝𝑟𝑒𝑠𝑠𝑢𝑟𝑒 × 𝑎𝑟𝑒𝑎 𝛥𝐹⃗ = −𝑝(𝑥 , 𝑦 , 𝑧 , 𝑡)𝑛⃗𝑑𝑆 = −𝑝(𝑥 , 𝑦 , 𝑧 , 𝑡)𝑑𝑆⃗
𝐹⃗ = − 𝑝𝑑𝑆⃗ = − 𝑝𝑑 𝑆⃗
which is the momentum equation for inviscid flows. We call it the Euler equation.
Employing the so-called Gradient Theorem,
𝑝𝑑 𝑆⃗ = (𝛻𝑝)𝑑∀
∀
𝜕𝜌𝑉⃗
𝑑∀ + 𝛻 • 𝜌𝑉⃗𝑉⃗ 𝑑∀= 𝜌𝐹⃗ 𝑑∀ − (𝛻𝑝)𝑑∀ ==>
∀ 𝜕𝑡 ∀ ∀ ∀
⃗
∰∀ + 𝛻 • 𝜌𝑉⃗𝑉⃗ 𝑑∀= ∰∀ 𝜌𝐹⃗ − 𝛻𝑝 𝑑∀ (A8a)
+ + + = 𝜌𝐹 − (A8c)
+ + =− (A8d)
F24AeroII—CN#01: Page 19 of 23
Energy Equation
The first law of thermodynamics dictates that,
(1) The time rate of change in energy for a CV
+ (2) the rate of energy (going out of the CV) through the CS
equal (3) the work rate done by body forces acting on the CV
+ (4) the work rate done by surface forces acting on the CS
We will follow the same procedure that we employed in getting the continuity and momentum equations.
The total energy inside a CV is,
𝐸𝑛𝑒𝑟𝑔𝑦 = 𝑚𝑎𝑠𝑠 × (𝑒𝑛𝑒𝑟𝑔𝑦 𝑝𝑒𝑟 𝑢𝑛𝑖𝑡 𝑚𝑎𝑠𝑠)
𝑉 𝑉
𝛥𝐸 = 𝜌(𝑥 , 𝑦 , 𝑧 , 𝑡)𝑑∀ 𝑒 + 𝐸(𝐶𝑉) = 𝜌 𝑒+ 𝑑∀
2 ∀ 2
where e is the internal energy (per unit mass).
The time rate of changes in the total energy inside the CV is,
𝜕𝐸(𝐶𝑉) 𝜕 𝑉 𝜕 𝑉
= 𝜌 𝑒+ 𝑑∀ = 𝜌 𝑒+ 𝑑∀
𝜕𝑡 𝜕𝑡 ∀ 2 ∀ 𝜕𝑡 2
Note again, the CV does not change with time, so we can take the partial derivative inside the integral. Through
the CS,
𝐸𝑛𝑒𝑟𝑔𝑦 𝐼𝑛 = 𝑀𝑎𝑠𝑠 𝐼𝑛 × (𝐸𝑛𝑒𝑟𝑔𝑦 𝑝𝑒𝑟 𝑢𝑛𝑖𝑡 𝑚𝑎𝑠𝑠)
𝑀𝑎𝑠𝑠 = 𝑑𝑒𝑛𝑠𝑖𝑡𝑦 × 𝑣𝑜𝑙𝑢𝑚𝑒 = 𝑑𝑒𝑛𝑠𝑖𝑡𝑦 × 𝑎𝑟𝑒𝑎 × 𝑙𝑒𝑛𝑔𝑡ℎ
𝑉 𝑉
𝛥𝐸 = −𝜌(𝑑𝑆)(𝑉 𝛥𝑡) 𝑒 + 𝐸(𝐶𝑆) = − 𝜌 𝑒+ 𝑉⃗ • 𝑑𝑆⃗
2 2
The work rate done by the body forces,
𝑓𝑜𝑟𝑐𝑒 × 𝐷𝑖𝑠 𝑡𝑎𝑛 𝑐 𝑒
𝐵𝑜𝑑𝑦𝐹𝑜𝑟𝑐𝑒𝑊𝑜𝑟𝑘𝑅𝑎𝑡𝑒 = = 𝑓𝑜𝑟𝑐𝑒 × 𝑣𝑒𝑙𝑜𝑐𝑖𝑡𝑦
𝑇𝑖𝑚𝑒𝑃𝑒𝑟𝑖𝑜𝑑
Putting the four terms above together, we have the inviscid energy equation,
F24AeroII—CN#01: Page 20 of 23
Appendix B: Air-Speed Record – Major Timeline (http://en.wikipedia.org/wiki/Flight_airspeed_record)
Airspeed
Date Pilot Aircraft Location
mph km/h
1903 Wilbur Wright 6.82 10.98 Wright Flyer Kitty Hawk, NC, USA
1947 Chuck Yeager 670.0 1078.00 Bell X-1- Rocket powered* Muroc, California, USA
From this point records are set at altitude and not sea level.
1955/10 Horace A. Hanes 822.1 1,323.0 N. American F-100C Super Sabre Palmdale, USA
1958/05 Cap. WW Irwin, USAF 1,404 2,259.5 Lockheed F-104C Starfighter Edwards AF Base, USA
1961/11/22 Robert G. Robinson 1606.3 2585.1 Modified M-D F-4 Phantom II Edwards AF Base, USA
Airspeed
Date Pilot Aircraft Location
mph km/h
1903 Wilbur Wright 6.82 10.98 Wright Flyer Kitty Hawk, NC, USA
F24AeroII—CN#01: Page 21 of 23
1911/06/12 Alfred Leblanc 77.640 125.00 Blériot
1918/08 Roland Rohlfs 163 262.30 Curtiss Wasp Not officially recognized.
F24AeroII—CN#01: Page 22 of 23
Rocket powered Me 163B
1944 Heini Dittmar 702 1130.00 Lagerlechfeld
V18*
1947/08/20 Cmdr. Turner Caldwell 640.663 1031.05 Douglas Skystreak Muroc, California, USA
1947/08/25 Major Marion E. Carl 650.796 1047.36 Douglas Skystreak Muroc, California, USA
1947 Chuck Yeager 670.0 1078.00 Bell X-1- Rocket powered* Muroc, California, USA
1952 J. Slade Nash 698.5 1,124.1 North American F-86D Sabre Salton Sea, USA
1953 William Barnes 715.7 1,151.8 North American F-86D Sabre Salton Sea, USA
1953/09/26 Mike Lithgow 735.7 1,184.0 Supermarine Swift F4 Castel Idris, Tripoli, Libya
1953/10/03 James B. Verdin, USN 752.9 1,211.5 Douglas F4D Skyray Salton Sea, USA
From this point records are set at altitude and not sea level.
1957/12 USAF 1,207.6 1,943.5 McDonnell F-101A Voodoo Edwards AF Base, USA
1958/05 Cap. WW Irwin, USAF 1,404 2,259.5 Lockheed F-104C Starfighter Edwards AF Base, USA
1959/10/31 Col. Georgii Mosolov 1,484 2,388.0 Ye-66 (proto MiG-21) USSR
1959/12 Maj. Joseph Rogers 1,525.9 2,455.7 Convair F-106 Delta Dart Edwards AF Base, USA
1961/11/22 Robert G. Robinson 1606.3 2585.1 Modified M-D F-4 Phantom II Edwards AF Base, USA
1962/07/07 Col. Georgii Mosolov 1665.9 2681.0 Mikoyan Gurevich Ye-166 USSR
F24AeroII—CN#01: Page 23 of 23