Airfoil Characteristics
Airfoil Characteristics
Drag – component of force parallel to the free-stream direction ahead of the wing.
Axial force – is that component of force parallel to the chord or axis of the wing or body and is equal to the form
drag plus skin friction at zero degree angle of attack.
Moment – is the resultant of the moments of all forces on a body about some chordwise reference point.
Resultant force – is the force representing the resultant of the addition of all local aerodynamic and viscous
forces on the wing from which lift and drag force are resolved.
Skin Friction – is the force created by the tendency of the layer of air next to the surface of the body to
cling to the layer next to it until the free-stream velocity is reached. This shearing action creates a drag on
the surface which is a function of the viscosity, velocity, density and the type of flow (laminar or turbulent)
of the air, as well as the airwise surface dimension.
Drag due to lift or “drag due to normal force” – it is the component of the normal force which is parallel
to the direction of the free-stream. It exists only when the airfoil is in an attitude to supply lift and hence is usually
given the name shown at the heading. This force is called induced drag in the subsonic case but is not so called in the
supersonic case because the type of flow over the lifting surface is of a different character. The term “induced drag” is
sometimes used to describe the drag created in the region of the tips of rectangular wings producing lift in supersonic
flow, however.
ENGR. JIM JOSHUA PANARES 2021
BASIC AERODYNAMIC CHARACTERISTICS OF THE WING
𝐴𝑥
𝐿 = 𝑁𝑦 − 𝐴𝑦 𝑁𝑥 𝛼 𝐴𝑦 𝑁𝑦
𝑁𝑦 𝐴
𝑐𝑜𝑠𝛼 =
𝑁 𝑁 𝐿 𝑁 𝑅
𝑁𝑦 = 𝑁𝑐𝑜𝑠𝛼 𝑁𝑦
𝐴𝑦 𝛼
𝑠𝑖𝑛𝛼 = 𝛼
𝐴
𝐴𝑦 = 𝐴𝑠𝑖𝑛𝛼 𝑀
𝐿 = 𝑁𝑐𝑜𝑠𝛼 − 𝐴𝑠𝑖𝑛𝛼
𝛼
Freestream Velocity, 𝑉 𝑁𝑥 𝐴𝑥
𝐴𝑦 𝛼 𝐷
𝐴
ENGR. JIM JOSHUA PANARES 2021
BASIC AERODYNAMIC CHARACTERISTICS OF THE WING
𝐴𝑥
𝐷 = 𝑁𝑥 + 𝐴𝑥 𝑁𝑥 𝛼 𝐴𝑦 𝑁𝑦
𝑁𝑥 𝐴
𝑠𝑖𝑛𝛼 =
𝑁 𝑁 𝐿 𝑁 𝑅
𝑁𝑥 = 𝑁𝑠𝑖𝑛𝛼 𝑁𝑦
𝐴𝑥 𝛼
𝑐𝑜𝑠𝛼 = 𝛼
𝐴
𝐴𝑥 = 𝐴𝑐𝑜𝑠𝛼 𝑀
𝐷 = 𝑁𝑠𝑖𝑛𝛼 + 𝐴𝑐𝑜𝑠𝛼
𝛼
𝑁𝑥 𝐴𝑥
Freestream Velocity, 𝑉 𝐴𝑦 𝛼 𝐷
𝐴
ENGR. JIM JOSHUA PANARES 2021
AERODYNAMIC CHARACTERISTICS OF A TWO-DIMENSIONAL
SYMMETRICAL DOUBLE-WEDGE AIRFOIL
Recalling Approximation Theory (Second Order):
∆𝑃
= 𝐶1 𝜃 + 𝐶2 𝜃 2
𝑞ത
2
𝐶1 =
𝑀1 2 − 1
① ②
2
𝛾𝑀1 4 + 𝑀1 2 − 2 𝛽
𝐶2 = 2
2 𝑀1 2 − 1
𝛼 ③ ④
Where:
∆𝑃
=pressure change across an OSW or EW
𝑞ത Note: 𝜃 =(+) towards
𝜃 =local angle of attack in radians between the surface and freestream Mach number the surface, (-) away
𝑀1 =freestream Mach number from the surface
ENGR. JIM JOSHUA PANARES 2021
AERODYNAMIC CHARACTERISTICS OF A TWO-DIMENSIONAL
SYMMETRICAL DOUBLE-WEDGE AIRFOIL
The values of 𝜃 to be used in the pressure equation
for each surface are as follows:
𝜃1 = −𝛼 + 𝛽
𝜃2 = −𝛼 − 𝛽
𝜃3 = 𝛼 + 𝛽 ① ②
𝛽
𝜃4 = 𝛼 − 𝛽
Where: 𝛼 ③ ④
𝛼 = angle between the freestream
direction and chord line of the wing
𝛽 = semi-vertex angle of leading and
trailing edge
ENGR. JIM JOSHUA PANARES 2021
PRESSURE DIFFERENCE OVER THE FRONT AND REAR HALVES OF
THE AIRFOIL AT ANGLE OF ATTACK
∆𝑃 ∆𝑃 ∆𝑃 Front
= −
𝑞ത 𝐹
𝑞ത 3
𝑞ത 1
①
Rear
∆𝑃
= 𝐶1 𝜃3 + 𝐶2 𝜃3 2 ②
𝑞ത 3 ③
∆𝑃 2
= 𝐶1 𝛼 + 𝛽 + 𝐶2 𝛼 + 𝛽
𝑞ത 3 ④
∆𝑃 Upper surface
= 𝛼𝐶1 + 𝛽𝐶1 + 𝛼 2 𝐶2 + 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2
𝑞ത 3
∆𝑃 ∆𝑃 −
= 𝐶1 𝜃1 + 𝐶2 𝜃1 2 Lifting surface
𝑞ത 1 𝑞ത +
∆𝑃 2
= 𝐶1 −𝛼 + 𝛽 + 𝐶2 −𝛼 + 𝛽
𝑞ത 1
Lower surface
∆𝑃
= −𝛼𝐶1 + 𝛽𝐶1 + 𝛼 2 𝐶2 − 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2
𝑞ത 1
PRESSURE DIFFERENCE OVER THE FRONT AND REAR HALVES OF
THE AIRFOIL AT ANGLE OF ATTACK
∆𝑃 ∆𝑃 ∆𝑃 Front
= − ①
𝑞ത 𝐹
𝑞ത 3
𝑞ത 1 Rear
∆𝑃 ②
= 𝛼𝐶1 + 𝛽𝐶1 + 𝛼 2 𝐶2 + 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2
𝑞ത 3 ③
∆𝑃
= −𝛼𝐶1 + 𝛽𝐶1 + 𝛼 2 𝐶2 − 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2
𝑞ത 1
④
Upper surface
∆𝑃
= 𝛼𝐶1 + 𝛽𝐶1 + 𝛼 2 𝐶2 + 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2
𝑞ത 𝐹
− −𝛼𝐶1 + 𝛽𝐶1 + 𝛼 2 𝐶2 − 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2 ∆𝑃 −
Lifting surface
∆𝑃 𝑞ത +
= 𝛼𝐶1 + 𝛽𝐶1 + 𝛼 2 𝐶2 + 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2
𝑞ത 𝐹
+𝛼𝐶1 − 𝛽𝐶1 − 𝛼 2 𝐶2 + 2𝛼𝛽𝐶2 − 𝛽 2 𝐶2
Lower surface
∆𝑃
= 2𝛼𝐶1 + 4𝛼𝛽𝐶2
𝑞ത 𝐹 ENGR. JIM JOSHUA PANARES 2021
PRESSURE DIFFERENCE OVER THE FRONT AND REAR HALVES OF
THE AIRFOIL AT ANGLE OF ATTACK
∆𝑃 ∆𝑃 ∆𝑃 Front
= −
𝑞ത 𝑅
𝑞ത 4
𝑞ത 2 ①
Rear
∆𝑃
= 𝐶1 𝜃4 + 𝐶2 𝜃4 2 ②
𝑞ത 4 ③
∆𝑃 2
= 𝐶1 𝛼 − 𝛽 + 𝐶2 𝛼 − 𝛽
𝑞ത 4 ④
∆𝑃 Upper surface
= 𝛼𝐶1 − 𝛽𝐶1 + 𝛼 2 𝐶2 − 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2
𝑞ത 4
∆𝑃 ∆𝑃 −
= 𝐶1 𝜃2 + 𝐶2 𝜃2 2 Lifting surface
𝑞ത 2 𝑞ത +
∆𝑃 2
= 𝐶1 −𝛼 − 𝛽 + 𝐶2 −𝛼 − 𝛽
𝑞ത 2
Lower surface
∆𝑃
= −𝛼𝐶1 − 𝛽𝐶1 + 𝛼 2 𝐶2 + 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2
𝑞ത 2
PRESSURE DIFFERENCE OVER THE FRONT AND REAR HALVES OF
THE AIRFOIL AT ANGLE OF ATTACK
∆𝑃 ∆𝑃 ∆𝑃 Front
= − ①
𝑞ത 𝑅
𝑞ത 4
𝑞ത 2 Rear
∆𝑃 ②
= 𝛼𝐶1 − 𝛽𝐶1 + 𝛼 2 𝐶2 − 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2
𝑞ത 4 ③
∆𝑃
= −𝛼𝐶1 − 𝛽𝐶1 + 𝛼 2 𝐶2 + 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2
𝑞ത 2
④
Upper surface
∆𝑃
= 𝛼𝐶1 − 𝛽𝐶1 + 𝛼 2 𝐶2 − 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2
𝑞ത 𝑅
− −𝛼𝐶1 − 𝛽𝐶1 + 𝛼 2 𝐶2 + 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2 ∆𝑃 −
Lifting surface
∆𝑃 𝑞ത +
= 𝛼𝐶1 − 𝛽𝐶1 + 𝛼 2 𝐶2 − 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2
𝑞ത 𝑅
+𝛼𝐶1 + 𝛽𝐶1 − 𝛼 2 𝐶2 − 2𝛼𝛽𝐶2 − 𝛽 2 𝐶2
Lower surface
∆𝑃
= 2𝛼𝐶1 − 4𝛼𝛽𝐶2
𝑞ത 𝑅 ENGR. JIM JOSHUA PANARES 2021
NORMAL FORCE COEFFICIENT
𝑁
The normal force coefficient 𝐶𝑁 = 𝑞𝑆
ത
will be the total difference between the upper and lower pressure difference
(divided by 2 to retain the coefficient form based on the total plan-form area).
𝑁 1 ∆𝑃 1 ∆𝑃
𝐶𝑁 = = +
𝑞ത 𝑆 2 𝑞ത 𝐹
2 𝑞ത 𝑅
1 1
𝐶𝑁 = 2𝛼𝐶1 + 4𝛼𝛽𝐶2 + 2𝛼𝐶1 − 4𝛼𝛽𝐶2
2 2
𝐶𝑁 = 𝛼𝐶1 + 2𝛼𝛽𝐶2 + 𝛼𝐶1 − 2𝛼𝛽𝐶2
4𝛼
𝐶𝑁 = 2𝛼𝐶1 =
𝑀1 2 − 1
Where:
2
𝐶1 =
𝑀1 2 − 1 ENGR. JIM JOSHUA PANARES 2021
PRESSURE DIFFERENCE OVER THE UPPER AND LOWER HALVES OF
THE AIRFOIL AT ANGLE OF ATTACK
∆𝑃 ∆𝑃 ∆𝑃 Front
= −
𝑞ത 𝑈
𝑞ത 1
𝑞ത 2
①
Rear
∆𝑃
= −𝛼𝐶1 + 𝛽𝐶1 + 𝛼 2 𝐶2 − 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2 ②
𝑞ത 1 ③
∆𝑃
= −𝛼𝐶1 − 𝛽𝐶1 + 𝛼 2 𝐶2 + 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2
𝑞ത 2 ④
∆𝑃 Upper surface
= −𝛼𝐶1 + 𝛽𝐶1 + 𝛼 2 𝐶2 − 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2
𝑞ത 𝑈
− −𝛼𝐶1 − 𝛽𝐶1 + 𝛼 2 𝐶2 + 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2 −
∆𝑃
∆𝑃 Lifting surface
= −𝛼𝐶1 + 𝛽𝐶1 + 𝛼 2 𝐶2 − 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2 𝑞ത +
𝑞ത 𝑈
+𝛼𝐶1 + 𝛽𝐶1 − 𝛼 2 𝐶2 − 2𝛼𝛽𝐶2 − 𝛽 2 𝐶2
Lower surface
∆𝑃
= 2𝛽𝐶1 − 4𝛼𝛽𝐶2
𝑞ത 𝑈
PRESSURE DIFFERENCE OVER THE UPPER AND LOWER HALVES OF
THE AIRFOIL AT ANGLE OF ATTACK
∆𝑃 ∆𝑃 ∆𝑃 Front
= −
𝑞ത 𝐿
𝑞ത 3
𝑞ത 4
①
Rear
∆𝑃
= 𝛼𝐶1 + 𝛽𝐶1 + 𝛼 2 𝐶2 + 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2 ②
𝑞ത 3 ③
∆𝑃
= 𝛼𝐶1 − 𝛽𝐶1 + 𝛼 2 𝐶2 − 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2
𝑞ത 4 ④
∆𝑃 Upper surface
= 𝛼𝐶1 + 𝛽𝐶1 + 𝛼 2 𝐶2 + 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2
𝑞ത 𝐿
− 𝛼𝐶1 − 𝛽𝐶1 + 𝛼 2 𝐶2 − 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2 −
∆𝑃
∆𝑃 Lifting surface
= 𝛼𝐶1 + 𝛽𝐶1 + 𝛼 2 𝐶2 + 2𝛼𝛽𝐶2 + 𝛽 2 𝐶2 𝑞ത +
𝑞ത 𝐿
−𝛼𝐶1 + 𝛽𝐶1 − 𝛼 2 𝐶2 + 2𝛼𝛽𝐶2 − 𝛽 2 𝐶2
∆𝑃 Lower surface
= 2𝛽𝐶1 + 4𝛼𝛽𝐶2
𝑞ത 𝐿
FORM DRAG COEFFICIENT (SYMMETRICAL DOUBLE-WEDGE
AIRFOIL)
1 ∆𝑃 1 ∆𝑃 For small 𝛽 in radian, 𝑡𝑎𝑛𝛽 ≈ 𝛽:
𝐶𝐷𝐹 ′ = +
2 𝑞ത 𝑈
2 𝑞ത 𝐿 4β2
𝐶𝐷𝐹 = 𝐶𝐷𝐹 ′β =
1 1
𝐶𝐷𝐹 ′ = 2𝛽𝐶1 + 4𝛼𝛽𝐶2 + 2𝛽𝐶1 − 4𝛼𝛽𝐶2 𝑀1 2 − 1
2 2
𝐶𝐷𝐹 ′ = 𝛽𝐶1 + 2𝛼𝛽𝐶2 + 𝛽𝐶1 − 2𝛼𝛽𝐶2 Where:
4𝛽 2
𝐶𝐷𝐹 ′ = 2𝛽𝐶1 = 𝐶1 =
2
𝑀1 − 1 𝑀1 2 − 1
4𝛼 2
𝐶𝐷𝐿 =
𝑀1 2 − 1
4β2 4𝛼 2
𝐶𝐷 = + + 𝐶𝐷𝑓
𝑀1 2 − 1 𝑀1 2 − 1
4 β2 + 𝛼 2
𝐶𝐷 = + 𝐶𝐷𝑓
𝑀1 2 − 1
𝑁 𝐷
The optimum angle of attack for maximum is found by differentiating the expression to determine the slope as a
𝐷 𝑁
𝐷
function of angle of attack and setting this expression equal to zero because the optimum angle occurs where is a
𝑁
minimum.
𝑁
The expression for optimum angle of attack is related to the maximum 𝐷 ratio by the relation:
1
𝐶𝑁
2 2
𝐾1 𝜏 + 𝐶𝐷𝑓 𝑀1 − 1 𝐶𝐷 𝑀𝐴𝑋
𝛼𝑜𝑝𝑡. (𝑟𝑎𝑑. ) = 𝛼𝑜𝑝𝑡. (𝑟𝑎𝑑. ) =
4 2
1
𝛼𝑜𝑝𝑡. (𝑟𝑎𝑑. ) =
𝐶𝑁
𝐾1 𝜏2 2
+ 𝐶𝐷𝑓 𝑀1 − 1 2
𝐶𝐷 𝑀𝐴𝑋
𝛼𝑜𝑝𝑡. (𝑟𝑎𝑑. ) =
2
ENGR. JIM JOSHUA PANARES 2021
MOMENT AND CENTER OF PRESSURE (OTHER CROSS-SECTIONAL
SHAPES)
𝐶𝑚0.5 = 2𝐶2 𝐴′ 𝛼 𝐴𝑐𝑠
Type 𝐴′ = 2
Where: 𝑐
𝐴𝑐𝑠 𝑡𝑐 𝛽 𝜏
′
𝐴 = 2 = 2= Double Wedge
𝑐 2𝑐 2 2
2𝜏
𝛽 = 2𝐴′ Biconvex
3
𝐶2 𝐴′ 𝑐 1 2𝜏
𝐶. 𝑃.0.5 = Modified Double Wedge 𝑎∗ = 3
𝐶1 3
(4) 0.06 2
Required: 𝐶𝐴 = + 0.0053
𝐶𝑁 , 𝐶𝐷 , 𝐶𝐿 , 𝐶𝑚0.5 , 𝐶. 𝑃.0.5 (1.8)2 −1
𝐶𝐴 = 0.01492140471
𝐶𝐿 = 𝐶𝑁 cos α − 𝐶𝐴 sin 𝛼
𝐶𝐿 = 0.3731670907 cos 8° − 0.01492140471 sin 8°
ENGR. JIM JOSHUA PANARES 2021
𝐶𝐿 = 0.3674587961
EXAMPLE PROBLEM NO. 1
Find 𝐶𝑁 , 𝐶𝐷 , 𝐶𝐿 , 𝐶𝑚0.5 , and 𝐶. 𝑃.0.5 of a symmetrical double-wedge airfoil of 6% thickness-to-chord ratio at 8-degree
angle of attack in a supersonic stream of Mach 1.8. Total skin friction drag coefficient is 0.0053.
Given: Solution:
Symmetrical Double-Wedge Airfoil 2 𝐶𝑚0.5 = 𝛼𝛽𝐶2
𝜏 = 6% = 0.06 𝐶1 =
𝜋 𝑟𝑎𝑑.
𝛼 = 8° 𝑀1 2 − 1 𝐶𝑚0.5 = 8° 0.06 1.617729592
𝑀1 = 1.8 180°
2
𝐶𝐷𝑓 = 0.0053 𝐶1 = 𝐶𝑚0.5 = 0.01355265974
(1.8)2 −1
𝛽𝐶2 ∙ 𝑐
𝐶1 = 1.33630621 𝐶. 𝑃.0.5 =
Required: 2 2𝐶1
𝐶𝑁 , 𝐶𝐷 , 𝐶𝐿 , 𝐶𝑚0.5 , 𝐶. 𝑃.0.5 𝛾𝑀1 4 + 𝑀1 2 − 2 0.06 1.617729592 𝑐
𝐶2 = 2 𝐶. 𝑃.0.5 =
2 𝑀1 2 − 1 2 1.33630621
1.4 (1.8)4 + 1.8 2 − 2 2 𝐶. 𝑃.0.5 0.03631793926𝑐
𝐶2 = 𝐶. 𝑃.0.5 = 3.631793926% c
2 1.8 2 − 1 2
𝐶2 = 1.617729592
ENGR. JIM JOSHUA PANARES 2021
EXAMPLE PROBLEM NO. 2
For the given symmetrical supersonic airfoil (𝜏=6%) below, determine the form drag coefficient. The freestream flow
has a Mach number of 2.
Given: Solution:
Modified Double Wedge 𝑥 1.08 𝑐𝑚 𝐾1 𝜏 2
𝜏 = 6% = 0.06
𝑎∗ =
𝑐 0.06 − 8 𝑐𝑚 𝐶𝐷𝐹 =
𝑡
𝑎∗ = 2 𝑀1 2 − 1
𝑀1 = 2 From 𝜏 = 𝑐 and 𝑐 = 8 + 2𝑥: 1.08 𝑐𝑚
𝑐 − 8 𝑐𝑚 0.06 7.2 0.06 2
Required: 2 ∗
𝑎 = 0.2777777778 𝐶𝐷𝐹 =
𝑎∗ = 2 2−1
𝐶𝐷𝐹 𝑐
𝑡 2 𝐶𝐷𝐹 = 0.01496491898
𝜏 − 8 𝑐𝑚 𝐾1 = ∗
2 𝑎
𝑎∗ = 𝑡 2
𝐾1 =
𝜏 0.2777777778
𝐾1 = 7.2