AE 429 - Aircraft Performance and Flight Mechanics
Aerodynamics of the Airplane
Airfoils
An airfoil is a section of a wing (or a fin, or a stabilizer, or a propeller, etc.)
Cambered
Symmetrical Laminar Flow
Reflexed
Supercritical
Airfoil Nomenclature
chord line
straight line connecting the LE and TE
mean camber line camber
locus of points halfway between upper and lower surfaces maximum distance between mean camber line and chord line
More definitions
relative wind
direction of v
angle-of-attack drag lift
angle between relative wind and chord line component of resultant aerodynamic force parallel to the relative wind component of aerodynamic force perpendicular to relative wind pressure distribution also produces rotational torque
moment
M c / 4 Moment about the wing quarterchord point
Review of Aircraft Metrics
Wing Area = S Wing Span = b Mean Chord = c = S b = ( ct + cr ) 2 Root Chord = c0 = cr Tip Chord = ct Taper Ratio = = ct cr 2 Aspect Ratio = AR = b c = b S
The source of aerodynamic forces Force due to pressure and force due to friction
R = pndS + kdS
S S
Velocity and pressure around an airfoil
Streamline upper surface
Pressure and velocity distributions = 0
p-p ( > 0)
p-p ( = 0) Pressure and velocity distributions > 0
Pressure distributions = 0 Symmetric airfoil
Ideal and real fluid flow about an airfoil
Center of Pressure
Center of Pressure c.p.
the resultant forces (lift and drag) acting at the center of pressure produce no moment (use centroid rule to find R location) since the pressure distribution over the airfoil changes with angleof-attack, the location of the center of pressure varies with angleof-attack
Aerodynamic Center a.c.
the a.c. is the point on the airfoil about which moments do not vary with angle of attack, assuming v is constant if L = 0, the moment is a pure couple equal to the moment about the aerodynamic center, ma.c. simple airfoil theory places the a.c.
at the quarter chord point for low speed airfoils. at the half chord point for supersonic airfoils
Aerodynamic Center and Center of Pressure
Dimensional analysis
Lift, drag, and pitching moment
Are functions of several variables:
D = f2 (V , , ,S , a ) L = f1 (V , , ,S , a )
Because it is not feasible to experimentally vary each of these independent variables, we seek to group parameters and thus reduce our effort Allows us to intelligently group the variables Is an application of the Buckingham Pi Theorem Assume: a b e f L = Z (V , , S d , a , )
Z , a, b, d , e, f
Dimensionless constants
M = f3 (V , , ,S , a )
Dimensional analysis
Dimensional analysis
Principle: dimensions on both sides of the equations must be identical
Fundamental units: m, l, t
Are related to physical quantities
ml t2 Equating the dimensions on left and right of the lift force equation L
ml l m 2 = l t2 t l3
a b
For example:
( )
l m t lt
Equating mass exponents Equating length exponents Equating time exponents
1= b+f 1 = a 3 b + 2d + e f
2 = a e f
Dimensional analysis
Continuing our analysis
solving the preceding three equations for a, b, and d (in terms of e and f) 1 f 2e f 1 e f L = Z (V ) f , S 2 , a ,
a 2 L = ZV S V V S
e f
a 1 and = V M our characteristic length
Noting that
has units of length, we choose c as
V c
f
Then we can replace
V S
with
Now, the lift equation becomes
e
1 1 L = ZV2 S M Rec
Force/moment coefficients
Now we define the airfoils section lift coefficient
Or we could have simply defined lift coefficient as
Dynamic pressure L 1 2 q S q = V 2 Notice that cl is dimensionless It is represented as a function of Mach number and Reynolds number, but since the dimensional analysis was carried out for a given angle of attack, an airfoils section lift coefficient depends on all three of these variables
1 cl Z M 2
1 Re c
L=
1 2 V Scl 2
cl
cl = f (, M , Re )
Force/moment coefficients
A similar analysis gives
Drag coefficient Moment coefficient
D = qScd
M = qSccm
Notice that this coefficient has c explicitly included This term accounts for the force x length units
Summarizing
cl
Where
L q S
cd
D q S
cm
M qSc
cl = f1 (, M , Re )
cd = f2 (, M , Re )
cm = f3 (, M , Re )
Generic lift coefficient variation with
Generic moment coefficient variation with
How to find the Aerodynamic Center xa.c. ?
M a.c. = Lxa.c. + M c / 4
M a .c . L xa.c. M c / 4 = + q Sc q S c q Sc
L
Mc/4
a .c .
c/4
xa.c.
cma.c. = cl dcma.c. d
xa.c. + cmc / 4 c dc x dc = l a.c. + mc / 4 d c d
xa.c. + m0 c
0 = a0
xa.c. m = 0 = const c a0
10
Generic lift coefficient variation with Mach number
Generic drag coefficient variation with Mach number
11
Cl - Experimental data are essential to aircraft design
NACA/NASA data camber changes L=0
Cl varies linearly with
This linear relationship breaks down when stall occurs
Cl - At high , the boundary layer tends to separate
Lift decreases Drag increases Moment becomes nose down
12
Pressure coefficient
Pressure coefficient
Let us define another coefficient that describes the pressure distribution over an airfoil surface
Cp p p p p = 1 q 2 V 2
The sketch shows how cp varies over both upper and lower surfaces cp can be measured experimentally in the wind tunnel
Pressure coefficient
Pressure coefficient versus mach
Prandtl-Glauert rule
Cp =
Cp ,0
2 1 M
13
Obtaining cl from cp
Lift per unit span
LE (leading edge) Ds cos = dx
ds dx
TE (trailing edge)
Lift is the upward force due to the pressure distribution on the lower surface minus the down-ward force due to the pressure distribution on the upper surface
L=
TE p LE l
cos ds
c
TE p cos ds LE u
Or, since ds cos = dx
L = pl dx pu dx
0 0 c
Obtaining cl from cp
Lift per unit span
Adding and subtracting p
L = (pl p )dx (pu p )dx 0 0
c c
Recalling the definition of lift coefficient
cl L L L = = qS qc(1 ) qc
Combining these two equations
cl =
Note:
1 c pl p 1 c pu p 0 q dx c 0 q dx c
C p ,l
It follows that:
pl p p p and C p ,u u q q
cl =
1 c C p ,l C p ,u dx c 0
14
Obtaining cl from cp Summary
If plots of pressure coefficient data over the upper and lower surfaces vs. chordwise distance (x) are available
The lift coefficient can be found as the net area between the upper and lower pressure coefficient curves divided by the chord length The section lift coefficient (cl ) equation is a good approximation only for small angles of attack
cl = 1 c C p ,l C p,u dx c 0
IF 0
Pressure coefficients defined
C p,l Pressure coefficient on the lower surface C p,u Pressure coefficient on the upper surface
Lift, Drag and Moment Coefficients
Lift, Wing and Drag coefficients are written as Cab Where a describes the type of coefficient (lift, drag or moment) and b describes what the coefficient is a function of (reference coefficient, angle of attack, sideslip) b also describes the units of the coefficient For example if b is 0 (reference coefficient) then there are no units. If b is (angle of attack) then the units are per degree or per radian
15
Lift Coefficients
Lift Coefficient vs Angle of Attack
Cl Stall
Cl0 0L
dCl = Cl d
0L = zero lift angle of attack Slope of the linear region gives the infinite wing lift coefficient The reference lift coefficient is given by the point where the angle of attack is zero
Angle of Attack (Degrees)
Moment Coefficients
Moment Coefficient vs Angle of Attack
Cm
Cm0
dCm = Cm d
Slope of the linear region gives the infinite wing moment coefficient The reference moment coefficient is given by the point where the angle of attack is zero
Angle of Attack (Degrees)
16
Drag Coefficients
Drag Coefficient vs Angle of Attack
Cd
Cd0
Angle of Attack (Degrees)
The reference drag coefficient is given by the point where the angle of attack is zero The other drag coefficients can be determined using Excel, Matlab etc to perform a quadratic regression
References:
Pinkerton, R. M. and Greenberg, H, "Aerodynamic Characteristics of a Large Number of Airfoils Tested in the VariableDensity Wind Tunnel", NACA Report No. 628 Abbott, I. H. and Von Doenhoff, A. E., "Theory of Wing Sections",Dover Publications, New York, 1959
17
NACA Airfoils
NACA Airfoils
18
NACA Airfoils
Boeing Airfoils
http://www.engr.utk.edu/~rbond/airfoil.html
19
NACA Airfoil data
AIRFOIL CLASSIFICATIONS
NACA 2412
NACA 23012
NACA 63-210
Some NACA Airfoils
NACA 0012: four digit airfoil. First two digits indicate no camber. Last two digits indicate max t/c=12 percent NACA 6412: This airfoil combines a 0012 thickness (four digit) with a twodigit 64 camber line. A 64 camber line has 6% max camber at 40% chord.
NACA 16-015: This airfoil is identical to a 0015-45 airfoil. The modified fourdigit thickness has a leading edge index of 4 and the maximum thickness is at 50% chord NACA 23012: This airfoil combines a 230 mean line (three-digit) with a 0012 thickness (four digit). A 230 mean line has CL=0.3 and maximum camber at 15% chord. NACA 63A010: First three characters: 63A series thickness. Fourth
character: no camber (CL design=0) Last two digits: 10 percent t/c
NACA 63A409: First three characters: 63A series thickness. Fourth
character: CL design=0.4 (63A airfoils always use 6-series modified mean line) Last two digits: 9 percent t/c
20
