Scribd
Upload
4 views3 pages

Lift and Drag Calculations for Airfoil

Given the Mach number and airfoil geometry, the document calculates the lift and drag per unit span. It uses oblique shock equations to find the pressure below the airfoil (pL) and isentropic flow equations to find the pressure above (pU). Plugging the calculated pressures into the lift and drag equations yields a lift of 110.0 kN/m and a drag of 35.7 kN/m.

Uploaded by

Raissa Benjamin
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
4 views3 pages

Lift and Drag Calculations for Airfoil

Given the Mach number and airfoil geometry, the document calculates the lift and drag per unit span. It uses oblique shock equations to find the pressure below the airfoil (pL) and isentropic flow equations to find the pressure above (pU). Plugging the calculated pressures into the lift and drag equations yields a lift of 110.0 kN/m and a drag of 35.7 kN/m.

Uploaded by

Raissa Benjamin
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

Problem 13.

204 [Difficulty: 4]

Given: Mach number and airfoil geometry

Find: Lift and drag per unit span

Solution:

The given or available data is: R = 286.9 J/kg.K


k = 1.4
p1 = 50 kPa
M1 = 1.75
= 18 o

c = 1 m

Equations and Computations:

The net force per unit span is F = (p L - p U)c

Hence, the lift force per unit span is


L = (p L - p U)c cos() (1)

The drag force per unit span is


D = (p L - p U)c sin() (2)

For the lower surface (oblique shock):

We need to find M 1n

The deflection angle is = 

= 18 o

From M 1 and , and Eq. 13.49


(using built-in function Theta (M , ,k ))

(13.49)

For = 18.0 o

= 62.9 o

(Use Goal Seek to vary  so that  is correct)


From M 1 and  M 1n = 1.56

From M 1n and p 1, and Eq. 13.48d


(using built-in function NormpfromM (M ,k ))

(13.48d)

p2 = 133.2 kPa

pL = p2

pL = 133.2 kPa

For the upper surface (isentropic expansion wave):

For isentropic flow p 0 = constant

p 02 = p 01

For p 01 we use Eq. 13.7a


(using built-in function Isenp (M , k ))

(13.7a)

p 01 = 266 kPa

p 02 = 266 kPa

For the deflection =  (Compression )

= 18.0 o

We use Eq. 13.55

(13.55)

and
Deflection = 2 - 1 = (M 2) - (M 1) (3)

From M 1 and Eq. 13.55 (using built-in function Omega (M , k ))

1 = 19.3 o

Applying Eq. 3 2 = 1 + 

2 = 37.3 o
From 2, and Eq. 13.55 (using built-in function Omega (M , k ))

For 2 = 37.3 o

M2 = 2.42

(Use Goal Seek to vary M 2 so that 2 is correct)

Hence for p 2 we use Eq. 13.7a


(using built-in function Isenp (M , k ))

p 2 = p 02/(p 02/p 2)

p2 = 17.6 kPa

pU = p2

pU = 17.6 kPa

From Eq. 1 L = 110.0 kN/m

From Eq. 2 D = 35.7 kN/m

You might also like