AE 451 Aeronautical Engineering Design I Aerodynamics. Prof. Dr. Serkan Özgen Dept. Aerospace Engineering December 2015

Similar documents
AE 451 Aeronautical Engineering Design I Aerodynamics. Prof. Dr. Serkan Özgen Dept. Aerospace Engineering December 2017

The polar for uncambered and cambered airfoils. The angel of attack α

Drag (2) Induced Drag Friction Drag Form Drag Wave Drag

To highlight the change in drag with lift: Drag = Zero-Lift Drag + Lift-Dependent Drag + Compressibility Drag

Definitions. Temperature: Property of the atmosphere (τ). Function of altitude. Pressure: Property of the atmosphere (p). Function of altitude.

APPENDIX C DRAG POLAR, STABILITY DERIVATIVES AND CHARACTERISTIC ROOTS OF A JET AIRPLANE (Lectures 37 to 40)

Chapter 5 Wing design - selection of wing parameters 2 Lecture 20 Topics

Semi-Empirical Prediction of Aircraft Low-Speed Aerodynamic Characteristics. Erik D. Olson. NASA Langley Research Center, Hampton, VA 23681

MDTS 5705 : Aerodynamics & Propulsion Lecture 2 : Missile lift and drag. G. Leng, MDTS, NUS

and K becoming functions of Mach number i.e.: (3.49)

The Importance of drag

Experimental Aerodynamics. Experimental Aerodynamics

Air Loads. Airfoil Geometry. Upper surface. Lower surface

Airfoils and Wings. Eugene M. Cliff

ν δ - 1 -

Mechanics of Flight. Warren F. Phillips. John Wiley & Sons, Inc. Professor Mechanical and Aerospace Engineering Utah State University WILEY

Supersonic Aerodynamics. Methods and Applications

Introduction to Aerospace Engineering

Drag Computation (1)

Flight Vehicle Terminology

PERFORMANCE ANALYSIS OF A PISTON ENGINED AIRPLANE - PIPER CHEROKEE PA E.G.TULAPURKARA S.ANANTH TEJAS M. KULKARNI REPORT NO: AE TR

Induced Drag and High-Speed Aerodynamics Robert Stengel, Aircraft Flight Dynamics, MAE 331, 2018

Masters in Mechanical Engineering Aerodynamics 1 st Semester 2015/16

AIRFRAME NOISE MODELING APPROPRIATE FOR MULTIDISCIPLINARY DESIGN AND OPTIMIZATION

Brenda M. Kulfan, John E. Bussoletti, and Craig L. Hilmes Boeing Commercial Airplane Group, Seattle, Washington, 98124

Based on the Aircraft Design Lecture Notes and Master Thesis of Roberto Segovia García. Hamburg University of Applied Sciences.

Given the water behaves as shown above, which direction will the cylinder rotate?

SPC Aerodynamics Course Assignment Due Date Monday 28 May 2018 at 11:30

Nonlinear Aerodynamic Predictions Of Aircraft and Missiles Employing Trailing-Edge Flaps

Drag Analysis of a Supermarine. Spitfire Mk V at Cruise Conditions

A Novel Airfoil Circulation Augment Flow Control Method Using Co-Flow Jet

Flight Dynamics and Control. Lecture 3: Longitudinal stability Derivatives G. Dimitriadis University of Liege

Nose Cone & Fin Optimization

ADVERSE REYNOLDS NUMBER EFFECT ON MAXIMUM LIFT OF TWO DIMENSIONAL AIRFOILS

Application of larger wings to create more lift force in aeroplane Lokendra Bhati 1,Rajat Malviya 2 1 B.E.in mechanical engineering, I.I.ST.

Compressible Potential Flow: The Full Potential Equation. Copyright 2009 Narayanan Komerath

BOUNDARY LAYER FLOWS HINCHEY

In-Flight Lift-Drag Characteristics for a Forward-Swept Wing Aircraft (and Comparisions With Contemporary Aircraft)

F-35: A case study. Eugene Heim Leifur Thor Leifsson Evan Neblett. AOE 4124 Configuration Aerodynamics Virginia Tech April 2003

Degree Project in Aeronautics

Lecture 7 Boundary Layer

Reynolds Number Effects on the Performance of Lateral Control Devices

An Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond

Development of a Wave Drag Prediction Tool for the Conceptual Design Phase

INSTITUTE OF AERONAUTICAL ENGINEERING (Autonomous) Dundigal, Hyderabad

Stability and Control Some Characteristics of Lifting Surfaces, and Pitch-Moments

RECENT near-sonic and low-sonic boom transport aircraft

Fundamentals of Airplane Flight Mechanics

9.4 Miscellaneous topics flight limitations, operating envelop and V-n diagram Flight limitations Operating envelop 9.4.

Basic Low Speed Aerodynamics

Relaminerization of a Highly Accelerated Flow on a Convex Curvature

ROAD MAP... D-1: Aerodynamics of 3-D Wings D-2: Boundary Layer and Viscous Effects D-3: XFLR (Aerodynamics Analysis Tool)

Aircraft Design Studies Based on the ATR 72

Improved Method for Prediction of Attainable Wing Leading-Edge Thrust

PEMP ACD2501. M.S. Ramaiah School of Advanced Studies, Bengaluru

Part 3. Stability and Transition

AOE 3104 Problem Sheet 10 (ans)

AEROSPACE ENGINEERING DEPARTMENT. Second Year - Second Term ( ) Fluid Mechanics & Gas Dynamics

Chapter 5 Lecture 19. Performance analysis I Steady level flight 3. Topics. Chapter-5

MACH NUMBERS BETWEEN 0.8 AND 2.2 AS DETERMINED FROM THE

The Wright Stuff. Engineering Analysis. Aaron Lostutter Adam Nelessen Brandon Perez Zev Vallance Jacob Vincent. November 2012

= 390 ( R). Then, for a given edge Mach number, M e., and ratio of wall temperature to adiabatic wall temperature T W = T W. T AW = T e 2 M e 2

Configuration Aerodynamics

LONGITUDINAL STABILITY AND TRIM OF AN ARIANE 5 FLY-BACK BOOSTER

( ) (where v = pr ) v V

ACTIVE SEPARATION CONTROL ON A SLATLESS 2D HIGH-LIFT WING SECTION

Applied Fluid Mechanics

COLD WEATHER PROTECTION GENERAL

Project 2 Conceptual Design

Reynolds number influence on high agility aircraft vortical flows

Rainer Buffo Eike Stumpf Institute of Aerospace Systems (ILR) RWTH Aachen University, Germany

AAE 251 Formulas. Standard Atmosphere. Compiled Fall 2016 by Nicholas D. Turo-Shields, student at Purdue University. Gradient Layer.

An Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond

Design and Computational Studies on Plain Flaps

Aerodynamic Classification of Swept-Wing Ice Accretion

Design of a Droopnose Configuration for a Coanda Active Flap Application. Marco Burnazzi and Rolf Radespiel

Given a stream function for a cylinder in a uniform flow with circulation: a) Sketch the flow pattern in terms of streamlines.

Wings and Bodies in Compressible Flows

RESEARCH MEMORANDUM FUGHT CHARACTERSSTICS OF A WINGIZSS ROCKET-POWERED MODEL WITH FOUR EXTERNALLY MOUNTED AIR-TO-APR MISSILES,$ Langley.Field, Va.

Continuity Equation for Compressible Flow

AIAA Investigation of Reynolds Number Effects on a Generic Fighter Configuration in the National Transonic Facility (Invited)

TUTORIAL SHEET 1. magnitude of P and the values of ø and θ. Ans: ø =74 0 and θ= 53 0

General Remarks and Instructions

Aerodynamics. High-Lift Devices

Introduction to Atmospheric Flight. Dr. Guven Aerospace Engineer (P.hD)

PEMP ACD2505. Finite Wing Theory. M.S. Ramaiah School of Advanced Studies, Bengaluru

Flight Testing of Laminar Flow Control in High-Speed Boundary Layers

A CONCEPTUAL STUDY OF AIRFOIL PERFORMANCE ENHANCEMENTS USING CFD. A Thesis by. Armin Ghoddoussi. Bachelor of Science, Sojo University, 1998

Test-to-Test Repeatability of Results From a Subsonic Wing-Body Configuration in the National Transonic Facility

University of California at Berkeley Department of Mechanical Engineering ME 163 ENGINEERING AERODYNAMICS FINAL EXAM, 13TH DECEMBER 2005

A Numerical Study of Circulation Control on a Flapless UAV

Transonic Aerodynamics Wind Tunnel Testing Considerations. W.H. Mason Configuration Aerodynamics Class

AN ENGINEERING LEVEL PREDICTION METHOD FOR NORMAL-FORCE INCREASE DUE TO WEDGE SECTIONS

What is the crack propagation rate for 7075-T6 aluminium alloy.

AFWAL-TM & AIRFRAME BRANCH JULY UNLIMITED

External Flow and Boundary Layer Concepts

Fundamentals of Aerodynamits

Learning. Goals. Page 1 of 6. Document : V1.1

Masters in Mechanical Engineering. Problems of incompressible viscous flow. 2µ dx y(y h)+ U h y 0 < y < h,

Mach number, relative thickness, sweep and lift coefficient of the wing - An empirical investigation of parameters and equations

Transcription:

AE 451 Aeronautical Engineering Design I Aerodynamics Prof. Dr. Serkan Özgen Dept. Aerospace Engineering December 2015

Lift curve 2

Lift curve slope 3

Subsonic lift curve slope C Lα = 2 + 4 + AR2 β 2 η 2 2πAR 1 + tan2 Λ max,t β 2 Valid until M dd, fairly accurate until M=1. β 2 = 1 M 2 S exposed S F η: airfoil efficiency, = 0.95 for most airfoils. F = 1.07 1 + d/b 2, fuselage lift factor. 4

Subsonic lift curve slope AR eff = AR(1 + 1.9 h b); effective AR with endplates, h: height of the endplate. AR eff 1.2AR; effective AR with winglets. 5

Supersonic lift curve slope Theory: C Lα = 4 β Practice: use the charts valid for trapezoidal wings. Correct the values read with S exposed S F 6

Supersonic lift curve slope 7

Maximum lift (clean) For moderate to high aspect ratio wings with moderate sweep and high leading edge radius: C L,max = 0.9c l,max cosλ c/4 If a wing has low AR or high sweep and a sharp leading edge, maximum lift will increase due to leading edge vortices. This is a function of the shape of the upper surface of the leading edge: y = y 0.06c y 0.015c 8

Maximum lift (clean) 9

Maximum lift (clean) 10

Maximum lift (clean) For high aspect ratio wings: C L,max = c l,max C L,max c l,max + C L,max, Correct for F, S exposed. s α CL,max = C L,max C L,α + α 0L + α CL,max 11

Maximum lift (clean) 12

Maximum lift (clean) 13

Maximum lift (clean) 14

Maximum lift (clean) A wing has low AR if: AR 3 C 1 +1 cos Λ LE. C L,max = C L,max,base + ΔC L,max. α CL,max = α CL,max,base + α CL,max 15

Maximum lift (clean) 16

Maximum lift (clean) 17

Maximum lift (clean) 18

Maximum lift (clean) 19

Maximum lift (clean) 20

Maximum lift (clean) At transonic speeds, maximum lift is limited by structural buffeting and controllability considerations rather than aerodynamics. 21

Maximum lift (with high lift devices) 22

Maximum lift (with high lift devices) 23

Maximum lift (with high lift devices) 24

Maximum lift (with high lift devices) Trailing edge devices decrease the stall angle of attack by increasing the pressure drop over the top of the airfoil promoting flow separation. In order to increase α stall a leading edge device must be used. ΔC L,max = 0.9Δc l,max S flapped S cosλ HL S flapped Δα 0L = Δα 0L,airfoil cosλ S HL HL: hinge line of the high lift device 25

Maximum lift (with high lift devices) For takeoff, the increments of about 60-80% of the increment calculated above should be used. Maximum lift occurs at a flap setting of about 40 o -45 o. Δ α0l,airfoil 15 o (landing setting), Δ α0l,airfoil 10 o (takeoff setting), 26

Maximum lift (with high lift devices) 27

Maximum lift (with high lift devices) 28

Maximum lift (with high lift devices) Leading edge devices increase lift by: Increasing camber, Increasing wing area, Delaying separation. Leading edge devices are particularly useful at high α. During takeoff and landing, they are useful when in combination with trailing edge devices as they prevent stall. 29

Estimation of C Do, equivalent skin friction method C Do = S wet S C fe C fe : equivalent skin friction coefficient is a function of the Reynolds number, Re. 30

Equivalent skin friction coefficients 31

Wetted area ratio 32

Estimation of C Do, component build-up method Total parasite drag coefficient: ) C Do subsonic = C fcff c Q c S wet,c S C fc : flat plate skin friction coefficient, C fc = C fc Re, M, k ; k: skin roughness. + C D,misc + C D,L&P FF c : form factor, estimates pressure drag due to separation, Q: interference factor. 33

Estimation of C Do, component build-up method C D,misc : drag of flaps, landing gears, upswept aft fuselage, base area. C D,L&P : drag of leakages and protuberances. 34

Flat plate skin friction coefficient Laminar flow: C f = 1.328 Re, Re = ρ V l, μ l: characteristic length. Turbulent flow: C f = 0.455 log Re 2.58 1 + 0.144M 2 0.65 35

Flat plate skin friction coefficient 36

Flat plate skin friction coefficient If the surface is rough, the skin friction coefficient will be higher. The smaller of the cut-off Reynolds number and the actual reynolds number shall be used. Subsonic flow: Re cutoff = 38.21 l k 1.053, Transonic or Supersonic flow: Re cutoff = 44.62 l k 1.053 M 1.16 37

Flat plate skin friction coefficient 38

Component form factors Wing, tail, strut and pylon: FF = 1 + 0.6 x/c m t c + 100 t c 4 1.34M 0.18 cos Λ m 0.28 x c m : chordwise location of the maximum thickness point, Λ m : sweep angle at the same location 39

Component form factors Fuselage and smooth canopy: FF = 1 + 60 f 3 + f 400 f = l d = l 4 π A max : fineness ratio. Nacelle and external stores: FF = 1 + 0.35 f 40

Component form factors For a tail surface with a hinged control surface: +10% A square sided fuselage: +40% For a two piece canopy: +40% For an external boundary-layer diverter for a fuselage mounted inlet: Double wedge: FF = 1 + d l, Single wedge: FF = 1 + 2d l. 41

Component form factors 42

Component interference factors Nacelle or external store mounted on wing or fuselage: Q=1.5. Nacelle or external store mounted on wing or fuselage: Q=1.3 (if mounted less than one diameter away). Nacelle or external store mounted on wing or fuselage: Q=1.1 (if mounted more than one diameter away). Wingtip mounted missiles: Q=1.25. High-wing, mid-wing or a well-filleted low-wing: Q=1.0. Unfilleted low-wing: Q=1.1-1.4. Fuselage: Q=1.0. Tail surfaces: Q=1.03 (V-tail), 1.08 (H-Tail), 1.04-1.05 (conventional tail). 43

Miscellaneous drag Upswept aft fuselage: D q = 3.83θ2.5 A max. 44

Miscellaneous drag Landing gear: summation of the drags of the wheels, struts, and other gear components. Q=1.2, *1.07 for retractable landing gears accounting for the hollow landing gear well. 45

Miscellaneous drag Flaps: ΔC Do,flap = F flap c flap c S flapped S δ flap 10 o. F flap = 0.0144: plain flaps, F flap = 0.0074: slotted flaps. Speed brakes: Fuselage mounted: D q = 1.0A frontal Wing mounted: D q = 1.6A frontal 46

Miscellaneous drag Canopies (transport and light aircraft): D q = 0.50A frontal,wind shield Cannon port: D q = 0.2 ft2. 47

Leakage and protuberance drag Antennas, lights, door edges, fuel vents, control surface external hinges, actuator fairings, rivets, rough or misaligned panels Jet transports and bombers: 2-5% parasite drag, Propeller aircraft: 5-10%, Fighters: 10-15% (old), 5-10% (new). 48

Supersonic Wave Drag For supersonic skin friction drag Q = FF = 1. C D0 ) s.sonic = C fcs wet + C S D,misc + C D,L&P + C D,wave Leakage and protuberance drag percentages apply only to skin-friction drag. For preliminary wave drag analysis (M 1.2): 0.77 D = E q wd 1 0.386(M 1.2) 0.57 1 πλ LE,deg wave 100 D q Sears Haack 49

Sears-Haack body 50

Supersonic wave drag E wd : wave drag efficiency factor. D =1.0 for a perfect Sears-Haack body, =1.2 for a smooth volume distribution, blended delta wing, =1.8-2.2 for a supersonic fighter, bomber. q Sears Haack = 9π 2 A max l 2 ; subtract inlet capture area. l: aircraft length length with constant cross sectional area. 51

Transonic wave drag Boeing formulation: M DD = M DD,L=0 LF DD 0.05C L,design wing. 52

Transonic wave drag 53

Transonic wave drag 54

Transonic wave drag L n : length of fuselage from nose to the location where fuselage cross section becomes constant. d: equivalent diameter of the fuselage there. 55

Transonic wave drag For initial analysis: 56

Transonic wave drag M 1.2: use supersonic wave drag expression. C D,wave M = 1.05 = C D,wave M = 1.2. C D,wave M = 1.0 = C D,wave M=1.05 2 M cr = M DD 0.08.. C D M DD = C D M cr + 0.002. 57

Complete drag build-up Subsonic drag: skin friction drag (including form factor and interference) + miscellaneous drag + Leakage & protuberance drag Supersonic drag: skin friction drag + miscellaneous drag + leakage and protuberance drag + wave drag. 58

Complete drag build-up 59

Complete drag build-up 60

Drag due to lift (induced drag) Induced drag coefficient: K = 1 πare Straight-winged airplane: e = 1.78 1 0.045AR 0.68 0.64 Swept winged airplane: e = 4.61 1 0.045AR 0.68 cos Λ 0.15 LE 3.1 (Λ LE < 30 o ) At supersonic speeds: K = AR(M2 1) 4AR M 2 1 2 cos Λ LE 61

Drag due to lift (induced drag) Flap effect on induced drag: ΔC Di = K f 2 (ΔC L,flap ) 2 cos Λ c 4, K f = 0.14: full span flaps, K f = 0.28: partial span flaps. 62