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

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

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

Chapter 5 Performance analysis I Steady level flight (Lectures 17 to 20) Keywords: Steady level flight equations of motion, minimum power required,

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

Chapter 3 Lecture 8. Drag polar 3. Topics. Chapter-3

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

Chapter 1 Lecture 2. Introduction 2. Topics. Chapter-1

Air Loads. Airfoil Geometry. Upper surface. Lower surface

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

Performance analysis II Steady climb, descent and glide 2

AIRFRAME NOISE MODELING APPROPRIATE FOR MULTIDISCIPLINARY DESIGN AND OPTIMIZATION

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

Chapter 8 Performance analysis IV Accelerated level flight and climb (Lecture 27) Keywords: Topics 8.1 Introduction 8.2 Accelerated level flight

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

Drag Computation (1)

Introduction to Aeronautics

Aircraft Structures Design Example

Chapter 4 Estimation of wing loading and thrust loading - 7 Lecture 15 Topics

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

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

AERO-STRUCTURAL MDO OF A SPAR-TANK-TYPE WING FOR AIRCRAFT CONCEPTUAL DESIGN

Flight Vehicle Terminology

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

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

Airfoils and Wings. Eugene M. Cliff

Masters in Mechanical Engineering Aerodynamics 1 st Semester 2015/16

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

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

Introduction to Aerospace Engineering

DEPARTMENT OF AEROSPACE ENGINEERING, IIT MADRAS M.Tech. Curriculum

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

The Importance of drag

ME 425: Aerodynamics

Performance analysis II Steady climb, descent and glide 3

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

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

Lecture-4. Flow Past Immersed Bodies

Research Article Aeroelastic Analysis of Wings in the Transonic Regime: Planform s Influence on the Dynamic Instability

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

Experimental Evaluation of Aerodynamics Characteristics of a Baseline Airfoil

MAE SUMMER 2015 HOMEWORK 1 SOLUTION

An Experimental Validation of Numerical Post-Stall Aerodynamic Characteristics of a Wing

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

Applied Fluid Mechanics

344 JAXA Special Publication JAXA-SP E 2. Prediction by the CFD Approach 2.1 Numerical Procedure The plane shape of the thin delta wing of the r

CHAPTER 3 ANALYSIS OF NACA 4 SERIES AIRFOILS

Chapter 1 Introduction (Lectures 1,2 and 3)

EVALUATION OF TRANSONIC BUFFETING ONSET BOUNDARY ESTIMATED BY TRAILING EDGE PRESSURE DIVERGENCE AND RMS DATA OF WING VIBRATION

Figure 1 UAV Geometry Base Case Design

Introduction to Aerospace Engineering

WALL ROUGHNESS EFFECTS ON SHOCK BOUNDARY LAYER INTERACTION FLOWS

Consider a wing of finite span with an elliptic circulation distribution:

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

Study of Preliminary Configuration Design of F-35 using simple CFD

THE EFFECT OF WING GEOMETRY ON LIFT AT SUPERSONIC SPEEDS

FREQUENCY DOMAIN FLUTTER ANALYSIS OF AIRCRAFT WING IN SUBSONIC FLOW

Configuration Aerodynamics

Supersonic Aerodynamics. Methods and Applications

Propeller theories. Blade element theory

1. Tasks of designing

Mestrado Integrado em Engenharia Mecânica Aerodynamics 1 st Semester 2012/13

AE 2020: Low Speed Aerodynamics. I. Introductory Remarks Read chapter 1 of Fundamentals of Aerodynamics by John D. Anderson

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

The E80 Wind Tunnel Experiment the experience will blow you away. by Professor Duron Spring 2012

Applied Aerodynamics - I

Chapter 8 Dynamic stability analysis II Longitudinal motion (Lectures 28 to 32)

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

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

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

Chapter 2 Earth s atmosphere (Lectures 4 and 5)

Flight and Orbital Mechanics

Aerodynamic Classification of Swept-Wing Ice Accretion

Static Aeroelasticity Considerations in Aerodynamic Adaptation of Wings for Low Drag

AE Stability and Control of Aerospace Vehicles

Continuity Equation for Compressible Flow

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

Aircraft Performance Sensitivity to Icing Cloud Conditions

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

The S411, S412, and S413 Airfoils

for what specific application did Henri Pitot develop the Pitot tube? what was the name of NACA s (now NASA) first research laboratory?

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

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

Chapter 2 Lecture 7 Longitudinal stick fixed static stability and control 4 Topics

Lecture 7 Boundary Layer

Design and Computational Studies on Plain Flaps

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

Improved Method for Prediction of Attainable Wing Leading-Edge Thrust

Aerodynamic Investigation of a 2D Wing and Flows in Ground Effect

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

High Speed Aerodynamics. Copyright 2009 Narayanan Komerath

Nose Cone & Fin Optimization

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

INSTITUTE OF AERONAUTICAL ENGINEERING (Autonomous) Dundigal, Hyderabad

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

Aerospace Engineering

Optimization Framework for Design of Morphing Wings

Chapter 4 Estimation of wing loading and thrust loading - 8 Lecture 16 Topics

External Flow and Boundary Layer Concepts

Numerical and Experimental Investigation on Aerodynamic Characteristics of SMA Actuated Smart Wing Model

FUNDAMENTALS OF AERODYNAMICS

Transcription:

Chapter 5 Wing design - selection of wing parameters Lecture 0 Topics 5..4 Effects of geometric parameters, Reynolds number and roughness on aerodynamic characteristics of airfoils 5..5 Choice of airfoil camber 5..6 Choice of airfoil thickness ratio (t/c) 5.3 Selection of wing parameters 5.3.1 Choice of aspect ratio (A) 5..4 Effect of geometric parameters, Reynolds number and roughness on aerodynamic characteristics of airfoils The important aerodynamic characteristics of airfoil from the point of view of design are angle of zero lift ( α ol ), maximum lift coefficient ( C lmax ), stall pattern, minimum drag coefficient ( C ), lift coefficient corresponding to C dmin which is dmin also called optimum lift coefficient (C lopt ), extent of drag bucket for low drag airfoils, moment coefficient about aerodynamic centre (C mac ) and critical Mach number. At subsonic speeds these characteristics are affected by geometrical parameters viz. camber, thickness ratio (t/c), airfoil shape, Reynolds number and roughness. Various chapters in Refs.5.1 and 5. contain information about characteristics of NACA airfoils. These effects can be summarized as follows. (i) The camber decides α 0l, C lopt and C mac. For a given family of airfoils, with increase of camber, α 0l and C mac become more negative whereas C lopt increases. (ii) The thickness ratio influences C dmin and C lmax. For a given family of airfoils, the minimum drag coefficient (C dmin ) increases with (t/c). The maximum lift coefficient (C lmax ) is highest for (t/c) between 1 to 16%. The stall pattern is also gradual for these thickness ratios. (iii) The Reynolds number (Re) mainly influences C lmax and C dmin. The former (C lmax ) increases with Re and the latter generally decreases with Re. As noted in Dept. of Aerospace Engg., Indian Institute of Technology, Madras 1

the previous subsection, the extent of drag bucket indicated by the nomenclature of the airfoil is at Re = 9 x 10 6. (iv) The surface roughness influences C lmax and C dmin. With increase of roughness C lmax decreases and C dmin increases. (v)the critical Mach number, in connection with the airfoil, is defined as the Free stream Mach number at which the maximum Mach number on the airfoil is unity. This quantity can be obtained theoretically by calculating the pressure distribution on the airfoil, but cannot be determined experimentally. However, when the critical Mach number is exceeded, the drag coefficient starts to increase. Making use of this behavior, the term Drag divergence Mach number (M D ) is defined as the Mach number at which the drag coefficient shows an increase of 0.00 over the subsonic drag value (Fig.5.4). Some authors (Ref.4.3) define M D as the Mach number at which the slope of the C d vs. M curve has a value of 0.1 i.e. (dc d /dm) = 0.1 Fig.5.4 Drag divergence Mach number Dept. of Aerospace Engg., Indian Institute of Technology, Madras

The drag divergence Mach number (M D ) depends on airfoil shape, thickness ratio, and lift coefficient. For a given airfoil M D is highest near C lopt. It decreases with thickness ratio. Supercritical airfoil For airplanes flying at high subsonic speeds the lift coefficient under cruising condition (C Lcr ) is around 0.5. At this value of lift coefficient, the older NACA airfoils have drag divergence Mach number (M D ) of around 0.68 for a thickness ratio (t/c) of around 15%. With the advancements in computational fluid dynamics (CFD) it was possible, in 1970 s to compute transonic flow past airfoils. This enabled design of improved airfoils, called supercritical airfoils, which have M D around 0.75 for t/c of 15% (Ref.1.1 part II, chapter 6). For comparison, the shapes of older airfoil (NACA 66-15) and a supercritical airfoil are shown in Fig.5.3 d and g. Note the flat upper surface of the supercritical airfoil. Refer chapter 3 of Ref.3.4 and Ref.5.4, for additional information. Remarks : (i) To illustrate the effects of Reynolds number, roughness, camber and thickness and ratio on C l, C d and C mac, the experimentally obtained variations (Ref.5.5) are presented in Figs.5.5 a to e. They related to NASA MS(1)-0317, MS(1)-0313 and LS(1)-0417 airfoils. It may be mentioned that the airfoil LS(1)-0417 is 17% thick airfoil with C ldesign of 0.4. It is designed specifically for low speed airplanes. Later NASA MS(1)-0317 with, thickness ratio of 17% was designed for applications to medium speed airplanes (M 0.7 ). The value of C ldesign is 0.3. The airfoil NASA MS(1)-0313 is similar to NASA MS(01)-0317, but has t/c of 13% Figure 5.5a shows the effect of varying Reynolds number from x10 6 to 1 x 10 6 on lift characteristics of NASA MS (01)-0317 airfoil. It is observed that C lmax increases from about 1.6 to.0. Note that α ol is -3 o. Figure 5.5b shows the effect of varying Reynolds number on C d vs C l curve of the same airfoil. It is seen that C dmin occurs around C l = 0.3 but is almost constant between C l = 0.1 to 0.5, effect of Reynolds number on C dmin is not very clear, Dept. of Aerospace Engg., Indian Institute of Technology, Madras 3

near C l = 0.3 but at higher values of C l (C l > 0.75) the values of C d decrease as Re increases. Figure 5.5c shows the C d vs C l curves with Re as parameter for the same airfoil but with rough surface details of roughness see Ref.5.5 comparing Figs.5.5 b and c. It is seen that C dmin is significantly higher for the rough airfoil as compared to the smooth one. The value of C dmin decreases with Re. Figure 5.5d compares the C l vs α, C d vs C l and C c vs C l curves for NASA m 4 LS(1)-0417 and MS(1)-0.317 airfoils. The cambers of the two airfoils are different, being higher for LS(1)-0.417. It is seen that α ol and C c are more m 4 negative for the LS(1)-0417. Figure 5.5 e compares C l vs α, C d vs C l and C c vs C l curves for NASA MS(1)- m 4 0317 and MS(1)-0313 airfoils. It is observed that the thinner airfoil has slightly lower value of C dmin. (ii)appendix F of Ref.1.0 gives the designations of airfoils used on many airplanes. Dept. of Aerospace Engg., Indian Institute of Technology, Madras 4

Fig.5.5a Effect of Reynolds number on C l vs α curve Airfoil : NASA MS(1)-0317; M = 0.15 ; smooth surface (Adapted from Ref.5.5) Dept. of Aerospace Engg., Indian Institute of Technology, Madras 5

Fig.5.5b Effect of Reynolds number on C d vs C l curve Airfoil : NASA MS(1)-0317; M = 0.15 ; smooth surface (Adapted from Ref.5.5) Dept. of Aerospace Engg., Indian Institute of Technology, Madras 6

Fig.5.5c Effect of Reynolds number on C d vs C l curve Airfoil : NASA MS(1)-0317; M = 0.15 ; Rough surface (Adapted from Ref.5.5) Dept. of Aerospace Engg., Indian Institute of Technology, Madras 7

Fig.5.5d C l vs α, C l vs C d, C l vs C mc/4 curves for NASA LS(1)-0417 and NASA MS(1)-0317 airfoils; Re = 6 x 10 6 ; M = 0.15 ; Rough surface (Adapted from Ref.5.5) Fig.5.5e C l vs α, C l vs C d, C l vs C mc/4 curves for NASA MS(1)-0317 and NASA MS(1)-0313 airfoils; Re = 6 x 10 6 ; M = 0.15 ; Rough surface (Adapted from Ref.5.5) Dept. of Aerospace Engg., Indian Institute of Technology, Madras 8

5..5 Choice of airfoil camber The choice of the airfoil for the airplane wing involves the selection of camber, thickness ratio and shape of the airfoil. The camber decides the C lopt of the airfoil and the thickness ratio decides the characteristics like C lmax, C dmin, drag divergence Mach number (M D ), weight of the wing and the stall pattern. For a good design, the camber should be chosen such that C lopt of the airfoil is close to the lift coefficient of the aircraft (C L ) in the flight corresponding to the mission of the airplane. This lift coefficient is called design lift coefficient (C Ldesign ). In most of the cases, this would correspond to the cruise flight condition. 1 Assuming L = W = ρvscl C Ldesign = Remark: W 1 ρvs ; ρ and V correspond to mission of the airplane e.g cruise The camber of the airfoil is chosen such that C lopt approximately equals C Ldesign. 5..6 Choice of airfoil thickness ratio (t/c) The thickness ratio (t/c) affects C dmin, C lmax, stall pattern, wing structural weight and M D. The influence of (t/c) on C dmin, C lmax and stall pattern has been dealt with in subsection 5..4. The following may be noted to understand the effect of thickness ratio (t/c) on the structural weight of the wing. The wing structure consists of spars (front and rear), stingers and skin (see Airbus 380 cut-away section in Appendix 1.1 and cut away drawing of airplanes in Ref.1.1).The spars are like I section beams. The flanges of the I section take the bending moment and the web takes the shear. If the wing section is thicker, then the spar flanges will be away from the centroidal axis of the section. Now, the bending moment resisted by an I section beam is proportional to the product of the area of the flange and the distance of flange from centroidal axis.thus, for a given bending moment, a thicker I beam would require lower Dept. of Aerospace Engg., Indian Institute of Technology, Madras 9

area of flange. Consequently, it would be lighter. Thus, a thicker wing will result in lighter wing. Reference 1.9, chapter 8; Ref.1.1, Pt.V, ch.9; Ref.1.15, chapter 6; Ref.1.18, chapter 15; Ref.1.19, chapter 8 and Ref.1.0, chapter 0 give formulae for the weight of wing (W W ) in terms of the geometrical parameters of the wing. Based on these, W w can be expressed as : c d e =CS A t/c 1+λ cos (5.5) a b W W w where, S W = wing area, A = aspect ratio, t/c = thickness ratio, λ = taper ratio and = sweep and C is a constant. The exponents a to e depend on the type of airplane. Their values lie in the following ranges. a = 0.6 to 0.76; b = 0.5 to 0.79; C = -0.3 to -0.4 ; d = 0.05 to 0.1 and e = -1. Remarks: (i) The final selection of the airfoil involve trade-off studies. It is seen that an increase in (t/c) results in increase of C lmax, decrease in wing weight and increase in C dmin. The trade-off studies would involve selecting different (t/c) values and examining which value gives minimum weight or fuel required etc. At the preliminary design stage the guidelines are obtained from the airfoils used on similar airplanes. Low speed airplanes have thickness ratio between 15 to 18%. NASA LS(1) 0417 is being used on low speed airplanes. NASA MS (01)-031 is being used on medium speed airplanes with turboprop engines. The high subsonic airplanes use supercritical airfoils of camber which would give C lopt = C ldesign and (t/c) around 14%. At supersonic speeds, C dmin is proportional to (t/c). These airplanes have (t/c) between 3 to 5%. Concorde airplane had biconvex airfoil of t/c = 0.035. (ii) Sometimes the (t/c) of the airfoil at the wing root is larger than the (t/c) of airfoil near wing tip. This is a compromise between the conflicting effects of increase of (t/c) on C dmin and the wing weight. Values of (t/c) = 0.18 at root and 0.15 at tip have been used. Dept. of Aerospace Engg., Indian Institute of Technology, Madras 10

5.3 Selection of wing parameters In this section, the selection of aspect ratio (A), sweep and taper ratioλare considered. 5.3.1 Choice of aspect ratio(a) Aspect ratio affects the slope of the lift curve structural weight of the wing and the wing span. a)effect of aspect ratio on slope of the lift curve C Lα, the induced drag C Di The slope of lift curve of a wing in subsonic flow for A > 4, is given by : (Ref.5.6, section 3.) π A C L = (5.6) tan Λ 1 A β + 4 + 1 + η β where β = 1 - M, η = C / ( π), lα Λ 1 = sweep of the half chord line, C lα is the slope of lift curve of the airfoil used on wing. Equation (5.6) shows that C Lα b)effect of aspect ratio on induced drag, decreases as aspect ratio decreases. The induced drag coefficient (C Di ) of a subsonic airplane is given by : CL C Di = 1+δ (5.7) πa where, δ depends on wing geometry i.e. aspect ratio, taper ratio and sweep. c)effect of aspect ratio on structural weight, the Equation (5.5) shows that the wing weight increases as A b where b = 0.5 to 0.79. The reason for this is as follows. As the aspect ratio increases the wingspan(b) increases (b= AS ). An increase in the span would increase the bending moment at the wing root. This would require higher moment of inertia of the spar and hence higher weight. Dept. of Aerospace Engg., Indian Institute of Technology, Madras 11

d) Effect of aspect ratio on span For a chosen wing area, the aspect ratio decides the span of the wing 1/ b= A S. In turn the span decides the hanger space needed for the airplane. Hence, for personal airplanes, a moderate aspect ratio of 6 to 7 is generally chosen. Further, the ride in turbulent weather is poor for a high aspect ratio wing. Hence, agricultural and other airplanes, which fly in proximity of ground, are subjected to air turbulence and have moderate aspect ratio of 6 to 7. Remark : The final choice of the aspect ratio would be arrived at after the trade-off studies which would involve selecting various values of aspect ratio and examining their effect on the criterion for the design of the particular airplane. At the preliminary design stage guidelines are obtained from the aspect ratios used on similar airplanes. Low speed airplanes of earlier designs had aspect ratio between 6 to 7.5, but the current trend is to choose between 7.5 to 8.5. The medium speed airplane, using turboprop engines, of earlier design had aspect ratio between 9 to 11. The current trend is the aspect ratio between 11 to 13. The high subsonic jet transport of earlier designs had aspect ratio between 7 to 8. The current trend is between 8.5 to 10.0. The trend towards higher aspect ratio appears to be due to availability of carbon epoxi material for fabrication of wing. This material is lighter than aluminium and has more stiffness. Reference 1.18, chapter 4 be referred for guidelines to select aspect ratios of other types of airplanes. Dept. of Aerospace Engg., Indian Institute of Technology, Madras 1