01 FlightPrinciples New
01 FlightPrinciples New
S. Altmeyer
Samuel Langley’s
“Aerodrome” (1903)
Leonardo da Vinci (≈ 1500) Advance in modern times
Outline
History
Let’s get Airborne
Environmental Model
The Earth
The Atmosphere
Mechanical Model
Reference Frames
Equations: Newton Laws for a Rigid Body
Aerodynamical Model
Fluid Dynamics Basics
Aerofoil Aerodynamic Forces
Wing Aerodynamic Forces
Aircraft Aerodynamic Forces
Propulsive Model
Simple Flight Equilibria
The Aircraft (Basic components)
The Aircraft (Basic components)
wing
vertical
stabiliser
horizontal
stabiliser
fuselage
engine
nacelle
The Aircraft (Basic components)
wing placement
wing planform
fuselage - basic structure
principle lifting force
tail assembly
horizontal & vertical
stabilizer
landing gear
→ stability
undercarriage
Hypotheses:
▶ Constant mass
▶ Rigid body
▶ Reflection-symmetric
The Wing - Basics
extrados Parts:
intrados ▶ Leading Edge
▶ Trailing Edge
▶ Intrados / Lower Surface
▶ Extrados / Upper Surface
Naming:
▶ S: projected surface
▶ b: span
ge ▶ λ = b 2 /S: aspect ratio
ed
d ing
lea ▶ ltip /lroot : taper
e ▶ Λ: sweep
edg
ling
trai ▶ Γ: dihedral
▶ θi : twist (incidence)
The Wing - control surfaces
1 winglet
2,3 ailerons (outboard, inboard)
4 flap deployment system
5,6 slats (inboard, outboard)
7,8 flaps (inboard, outboard)
9,10 spoilers (outboard, inboard)
The Aerofoil
Parts:
▶ leading edge
leading upp
chor camber er s
urfa
edge d lin line
e ce ▶ trailing edge
trailing
lower edge
surface ▶ intrados / lower surface
▶ extrados / upper surface
Flight parameters:
▶ α: angle of attack Geometry parameters:
▶ l: chord
▶ γ: slope
▶ t: thickness
▶ θ = α + γ: pitch angle
▶ c: camber
▶ V : aerodynamic speed
Outline
History
Let’s get Airborne
Environmental Model
The Earth
The Atmosphere
Mechanical Model
Reference Frames
Equations: Newton Laws for a Rigid Body
Aerodynamical Model
Fluid Dynamics Basics
Aerofoil Aerodynamic Forces
Wing Aerodynamic Forces
Aircraft Aerodynamic Forces
Propulsive Model
Simple Flight Equilibria
The Earth
Hypotheses: (M = 2, h = 20km)
▶ Steady earth: (0.8%)
acoriolis 2Ωe V
= ≪1
agravity g
acentripetal V2
= ≪1
agravity Re g
g0 − gh 1
=1− ≪1
g0 (1 + h/Re )2
Air
Properties (sea level values ISA):
▶ p: Pressure (p0 = 101325Pa)
▶ T : Temperature (T0 = 288.15K )
▶ ρ: Density (ρ0 = 1.225kg /m3 )
▶ µ: Dynamic viscosity (µ0 = 1.78 · 10−5 Pa · s )
▶ R = R/M: Air constant (287.14J/(Kg K ))
Composition of dry air ▶ Cp , Cv : specific heat capacities (J/(Kg K ))
(volume fraction): Derived Properties:
▶ 78.1% Nitrogen (N2 ) ▶ γ = Cp /Cv : Adiabatic index (1.4)
√
▶ 20.9% Oxigen (O2 ) ▶ a = γRT : Speed of sound (m/s)
▶ 0.9% Argon (Ar) ▶ ν = µ/ρ: Kinematic viscosity (m2 /s)
▶ 0.1% Other gases Flow parameters:
Water vapor (H2 O): 0.4% ▶ M = V /a: Mach number
(1-4% at sea level) ▶ Re = VL/ν: Reynolds number
The International Standard Atmosphere
ISA (ISA online calculator) Hypotheses:
▶ Fixed atmosphere
▶ Ideal gas: p/ρ = RT
▶ Laplace’s law: dp = −ρgdh
▶ Temperature lapse rate: dT = Th dh
▶ Sutherland’s law for µ
Earth Axes:
▶ ze : Vertical axis. Oriented by local gravitational acceleration vector.
▶ xe : Latitude axis. South-North direction.
▶ ye : Longitude axis. West-East direction.
Body Axes:
▶ x: Fuselage (roll) axis. Forward pointing, aligned with the fuselage.
▶ z: Belly (yaw) axis. In the airplane symmetry plane.
▶ y: Right wing (pitch) axis. Completing a direct orthonormal basis,
aligned with the right wing.
Aerodynamic Axes:
▶ xa : Oriented by the aerodynamic velocity vector.
▶ za : In the airplane symmetry plane.
▶ ya : Completes a direct orthonormal basis.
Coordinate Changes
Body Axes (Aircraft attitude)
lateral vertical
Aerodynamic Axes
plane plane
symmetry
plane
horizontal
plane
"wings"
plane
symmetry
plane
The full aircraft equations of motion reflect a rather complicated relationship between
the forces and moments on the aircraft, and resulting aircraft motion
Newton’s second law for translational F = m a and rotational L = I ω = r × p motions
→ 6 second-order equations:
▶ Linear Momentum Equation
d 2 r⃗G X
m 2
= F⃗ext = F⃗ + T
⃗ +W
⃗
dt
d X
⃗ G ,ext = M⃗F + M⃗T
(IG ω
⃗) = M
dt
along a streamline:
pstatic + pdynamic = ptotal
▶ restrictions:
inviscid, steady, incompressible
(small v ), no heat addition,
negligible change in height
Generation of Lift – Theories
But how?
▶ inviscid (& incompressible) potential flow:
▶ flow leaves unperturbed
▶ no net (drag) force: D’Alembert’s paradox
▶ non-physical effects at trailing edge
▶ inviscid flow + circulation:
▶ corrects for non-physical effects (Kutta)
▶ provides lift (Kutta-Zhoukowski)
▶ But how is this accomplished?
Generation of Lift (and drag)
▶ inviscid flow + circulation:
▶ viscous flow:
▶ fluid particles interact at the molecular level
▶ adjacent fluid layers with different velocity:
▶ generate vorticity (circulation) −→ LIFT
▶ produce shear (friction) −→ DRAG
▶ Where?
▶ in a thin layer around solid obstacles: BOUNDARY LAYER
Generation of Lift
Aerodynamic forces:
▶ Pressure forces act normal
(perpendicular) to surface
▶ Force on the body is the vector
sum of pressure * area around the
entirePsolid body
F⃗ = surface p⃗
H
nA = p⃗ ndA
Lift = Fnormal Drag = Fstream
For a lifting airfoil:
(ps + pdyn = ptotal = const.)
▶ the surface static pressure varies from
top to bottom and front to back
▶ the average pressure on the upper
surface is lower than the average
pressure on the lower surface
▶ The pressure differende produces Lift!
Generation of Lift
Question:
Why does the surface static pressure
vary top to bottom and front to back?
▶ The ideal surface velocity varies
from top to bottom and front to
back.
▶ The surface of the foil is a
streamline, so Bernoulli’s equation
relates the surface velocity to
surface pressure.
▶ When we include viscosity, the
surface velocity is zero, but the
local velocity varies at the
boundary layer edge.
Generation of Lift
The flow around the airfoil must satisfy the conservation laws:
▶ Conservation of mass (continuity)
▶ Conservation of momentum (2 or 3 components)
▶ Conservation of energy
=⇒ LIFT
without viscosity NO lift
Lift (velocity & pressure distribution)
Ideal fluid flow about an airfoil Real fluid flow vs ideal fluid
Boundary layer
Ideal fluid
▶ slip over the surface
▶ NO drag (d’Alembert)
Real fluid
▶ laminar
▶ no-slip condition
▶ → boundary layer (0 to V∞ )
▶ relative velocities → internal
friction
▶ → drag
▶ turbulent
▶ random motion in the BL
▶ BL has same size as laminar
▶ velocity builds up more
quickly from the wall
▶ fluid further away reenergize
Boundarylayer flows: the slower moving fluid near
ideal fluid →← real fluid the wall
Flow regimes
Shock-wave formation
Useful simplifications in aerodynamics:
▶ Inviscid / viscous regions decoupling:
Boundary layer hypothesis:
▶ Inviscid flow in most of the
domain −→ Euler Equations
Re ⇈
▶ viscosity effects confined to
narrow layer close to solid
obstacles −→ Boundary Layer
Equations
▶ Low-pass filtering + turbulence Flow regimes for viscous fluids:
modelling ▶ Re −→ 0: Stokes flow (no inertial
▶ ideal gas (simple state equations) effects)
▶ Incompressible flow (solenoidal flow) ▶ moderate Re: laminar flow
M⇊ ▶ high Re: turbulent flow
▶ Stationary flow (permanent regime) ▶ increased apparent viscosity
▶ Geometric simplifications (2D flow, (friction)
axisymmetric flow...) ▶ higher momentum transport
(later bl detachment)
Flow regimes
Flow regimes for compressible fluids: Shock-wave formation
▶ Low Subsonic / Incompressible:
M ≲ 0.3
▶ Subsonic Compressible:
0.3 ≲ M ≲ 0.75
▶ Transonic: 0.75 ≲ M ≲ 1.2
▶ Supersonic: M ≳ 1.2
▶ Hypersonic: M ≳ 5
→ spanwise flow
L/D with M Double delta wing
Stall fences and
vortex generators
Sonic Boom - shock-wave formation
aerofoil 3D body
CL & CD functions of:
surface A l S
shape, α, M = V /a, Re = Vl/ν...
units [F ]/[L] [F ]
α Angle of Attack (AoA): angle between plane of the wing (airfoil chord)
and direction of motion (free stream velocity V∞ )
Lift (Aerofoil)
1 Plain Aerofoil
2 Hinged Flap
3 Split Flap
4 Leading-Edge Slat
5 Single-Slotted flap
6 Double-Slotted flap
7 Double-Slotted flap + Leading-Edge Slat
8 Boundary-Layer Suction
Lift (Aerofoil + High-lift devices)
Aerodynamic Centre (Aerofoil)
Centre of Lift (pressure) (xL ):
Flap deployment:
▶ CD0 increases
▶ k decreases
Consequence:
▶ Deploy flaps when
high CL is required
▶ Stow flaps for low
CL demand
Go to Example
Drag (Wing + High-lift devices)
▶ increase camber
→ increase lift
▶ boundary layer control →
increase max lift and αstall
Flaps deployed:
▶ CL increases
▶ CD0 increases
Aerodynamic Forces (Aircraft)
Aerodynamic Reference Frame:
a
−CD Cl
1 1 1 1
F⃗ = ρSV 2 C⃗F = ρSV 2 CQ ⃗ = ρlSV 2 C⃗M = ρlSV 2 C a
M m
2 2 2 2
−CL Cna
Body Axes:
T T
C⃗F = (−CA , CY , −CN ) M⃗F = (Cl , Cm , Cn )
Lift-to-Drag Ratio (E or f ):
L CL 1
E (α, M, Re) = = =
D CD tan γ
Depends on:
▶ Aerofoil characteristics
▶ Wing planform and configuration
▶ Angle of attack α
▶ Mach number M
▶ Reynolds number Re
Maximum Efficiency / Optimal Lift:
r
CD0 1
CLopt = Emax = p
k 2 CD0 k
Go to Example
Outline
History
Let’s get Airborne
Environmental Model
The Earth
The Atmosphere
Mechanical Model
Reference Frames
Equations: Newton Laws for a Rigid Body
Aerodynamical Model
Fluid Dynamics Basics
Aerofoil Aerodynamic Forces
Wing Aerodynamic Forces
Aircraft Aerodynamic Forces
Propulsive Model
Simple Flight Equilibria
Propulsive Model
Simple model:
T = kT ρV λT δx
▶ kT : some constant
▶ λT : engine-dependent exponent
≃ −1 propeller / turbofan
≃0 turbojet
≃1 turbojet + afterburner
≃2 ram-jet
▶ δx : thrust lever position
Outline
History
Let’s get Airborne
Environmental Model
The Earth
The Atmosphere
Mechanical Model
Reference Frames
Equations: Newton Laws for a Rigid Body
Aerodynamical Model
Fluid Dynamics Basics
Aerofoil Aerodynamic Forces
Wing Aerodynamic Forces
Aircraft Aerodynamic Forces
Propulsive Model
Simple Flight Equilibria
Longitudinal Equilibrium
Banked Turn:
▶ x-equilibrium: T = D
▶ y -equilibrium:
L sin(ϕ) = mV 2 /Rt
▶ z-equilibrium:
L cos(ϕ) = W = mg
Load Factor: n = L/W = 1/ cos(ϕ)
Stall Speed: √ p
Vsn̸=1 = Vsn=1 n = Vsn=1 / cos(ϕ)
Go to Example
Load Factor & Stall Speed
Go to Example
Example
◦
2. Its lift coefficient is CL = 0.4 + 0.2α, with αs = 11
2.1 Zero-lift angle of attack (α0 )
2.2 α = 0 lift coefficient (CL0 )
2.3 Maximum lift coefficient (CLmax ) back
4. With deployed flaps: CLf = 0.6 + 0.2α and CDf = 0.02 + 0.03CL2f
4.1 Optimal lift coefficient for configuration change (CLopt )
4.2 Optimal angle of attack (αopt ) back
Example
Airplane characteristics: m = 9000kg , S = 32m2 , b = 20m
5. Straight & level flight (SL and h = 30000ft)
5.1 Speed for α = 0 (V (0))
5.2 Speed for maximum efficiency (V |Emax )
5.3 Speed for optimal flap removal (Vfr ) back
7. Stall speed
7.1 SL straight-flight (VEASs )
7.2 h = 30000ft straight-flight (Vs )
7.3 h = 30000ft banked turn with previous ϕSL (Vs ) back