Aerodynamics. Airfoils Kármán-Treftz mapping. The exponent k defines the internal angle of the trailing edge, τ, from. ( k)

Similar documents
Introduction to Aeronautics

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

CHAPTER 3 ANALYSIS OF NACA 4 SERIES AIRFOILS

Masters in Mechanical Engineering Aerodynamics 1 st Semester 2015/16

An Internet Book on Fluid Dynamics. Joukowski Airfoils

Incompressible Flow Over Airfoils

1. Fluid Dynamics Around Airfoils

Copyright 2007 N. Komerath. Other rights may be specified with individual items. All rights reserved.

Aerodynamics. High-Lift Devices

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

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

AIRFOIL DESIGN PROCEDURE A MODIFIED THEODORSEN E-FUNCTION. by Raymond L. Barger. Langley Research Center NASA TECHNICAL NOTE NASA TN D-7741

Flight Vehicle Terminology

Theory of turbo machinery. Chapter 3

COMPUTATIONAL SIMULATION OF THE FLOW PAST AN AIRFOIL FOR AN UNMANNED AERIAL VEHICLE

ν δ - 1 -

Airfoils and Wings. Eugene M. Cliff

An Approach to the Constrained Design of Natural Laminar Flow Airfoils

Lifting Surfaces. Lifting Surfaces

INSTITUTE OF AERONAUTICAL ENGINEERING (Autonomous) Dundigal, Hyderabad

Wings and Bodies in Compressible Flows

Lifting Surfaces. Lifting Surfaces

NUMERICAL INVESTIGATION OF VERTICAL AXIS WIND TURBINE WITH TWIST ANGLE IN BLADES

Chapter three. Two-dimensional Cascades. Laith Batarseh

Lift Enhancement on Unconventional Airfoils

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

Syllabus for AE3610, Aerodynamics I

ENGR 4011 Resistance & Propulsion of Ships Assignment 4: 2017

Active Control of Separated Cascade Flow

Propeller theories. Blade element theory

A Study of Transonic Flow and Airfoils. Presented by: Huiliang Lui 30 th April 2007

Explicit algebraic Reynolds stress models for internal flows

FIXED BEAMS CONTINUOUS BEAMS

Exercise 9: Model of a Submarine

Lecture 5: Pure bending of beams

HIGHER SCHOOL CERTIFICATE EXAMINATION MATHEMATICS 4 UNIT (ADDITIONAL) Time allowed Three hours (Plus 5 minutes reading time)

FLOW AROUND A SYMMETRIC OBSTACLE

MATH 566: FINAL PROJECT

Lab Reports Due on Monday, 11/24/2014

Lecture 7 Boundary Layer

COURSE ON VEHICLE AERODYNAMICS Prof. Tamás Lajos University of Rome La Sapienza 1999

Optimization of Unstalled Pitching and Plunging Motion of an Airfoil

Aero III/IV Conformal Mapping

Coupled Fluid and Heat Flow Analysis Around NACA Aerofoil Profiles at Various Mach Numbers

Correction of Wind Tunnel Results for the Airfoils of ITA s Unmanned Aerial Vehicle

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

a) Derive general expressions for the stream function Ψ and the velocity potential function φ for the combined flow. [12 Marks]

Shape Optimization of Low Speed Airfoils using MATLAB and Automatic Differentiation. Christian Wauquiez

High Speed Aerodynamics. Copyright 2009 Narayanan Komerath

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

HEAT-3 APPLICATION OF DERIVATIVES BY ABHIJIT KUMAR JHA MAX-MARKS-(112(3)+20(5)=436)

Laminar Flow. Chapter ZERO PRESSURE GRADIENT

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

Section 4.3: Quadratic Formula

Limit Cycle Oscillations of a Typical Airfoil in Transonic Flow

Restriction/ Classification Cancelled

Aerodynamics. Basic Aerodynamics. Continuity equation (mass conserved) Some thermodynamics. Energy equation (energy conserved)

2-D and 3-D Assessment of Cambered and Symmetric Airfoils: A CFD Study

FLUID MECHANICS. Lecture 7 Exact solutions

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

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

Fluid structure interaction of a parked wind turbine airfoil

Lifting Airfoils in Incompressible Irrotational Flow. AA210b Lecture 3 January 13, AA210b - Fundamentals of Compressible Flow II 1

Numerical study of battle damaged two-dimensional wings

chapter 1 vector geometry solutions V Consider the parallelogram shown alongside. Which of the following statements are true?

Review: Properties of Exponents (Allow students to come up with these on their own.) m n m n. a a a. n n n m. a a a. a b a

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

Drag Computation (1)

FLUID MECHANICS. Atmosphere, Ocean. Aerodynamics. Energy conversion. Transport of heat/other. Numerous industrial processes

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

Large-eddy simulations for wind turbine blade: rotational augmentation and dynamic stall

Department of Energy Sciences, LTH

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

Numerical Study on Performance of Curved Wind Turbine Blade for Loads Reduction

2016 SEC 4 ADDITIONAL MATHEMATICS CW & HW

Fluid Mechanics Prof. T. I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay

Q. 1 Q. 25 carry one mark each.

4-1 Graphing Quadratic Functions

Design and Computational Studies on Plain Flaps

ANALYSIS OF CURRENT CLOSE TO THE SURFACE OF NET STRUCTURES Mathias Paschen, University of Rostock, Germany

Thin airfoil theory. Chapter Compressible potential flow The full potential equation

Performance of an Axial Cascade

Precalculus Summer Packet

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

A2 Assignment lambda Cover Sheet. Ready. Done BP. Question. Aa C4 Integration 1 1. C4 Integration 3

Aerodynamic Rotor Model for Unsteady Flow and Wake Impact

Implementing a Partitioned Algorithm for Fluid-Structure Interaction of Flexible Flapping Wings within Overture

Correction of the wall interference effects in wind tunnel experiments

Information Knowledge


2.00AJ / 16.00AJ Exploring Sea, Space, & Earth: Fundamentals of Engineering Design Spring 2009

Nonlinear Frequency Domain Methods Applied to the Euler and Navier-Stokes Equations p.1/50

Chapter 6: Incompressible Inviscid Flow

Aircraft Design I Tail loads

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

NATIONAL SENIOR CERTIFICATE GRADE 11

Aeroelastic modelling of vertical axis wind turbines

SOUTH AFRICAN TERTIARY MATHEMATICS OLYMPIAD

Aerodynamic Performance 1. Figure 1: Flowfield of a Wind Turbine and Actuator disc. Table 1: Properties of the actuator disk.

Inviscid & Incompressible flow

Transcription:

z = b Aerodynamics Kármán-Treftz mapping ( ζ + b) + ( ζ b) ( ζ + b) ( ζ b) The eponent defines the internal angle of the trailing edge, τ, from ( ) τ = π = = corresponds to the Jouowsi mapping z z + b b τ π ζ b = ζ + b

z = b Aerodynamics Kármán-Treftz mapping ( ζ + b) + ( ζ b) ( ζ + b) ( ζ b) = corresponds to a dividing streamline that leaves the trailing edge with tangential continuity. The trailing edge is not a stagnation point z z + b b ζ b = ζ + b τ = 0

z = b Aerodynamics Kármán-Treftz mapping ( ζ + b) + ( ζ b) ( ζ + b) ( ζ b) =1,95 corresponds to a dividing streamline that forms a dihedral at the trailing edge with an angle smaller than π. The trailing edge is a stagnation point z z + b b o τ = 9 ζ b = ζ + b

Kármán-Treftz mapping 1. Cylinder centre at the origin η z = b ( ζ + b) + ( ζ b) ( ζ + b) ( ζ b) y b ξ -b b =1,95, c = 3, 9b

η z b Aerodynamics Kármán-Treftz mapping. Cylinder centre on the imaginary ais = ( ζ + b) + ( ζ b) ( ζ + b) ( ζ b) y b ξ -b b =1,95, c = 3, 9b

η Aerodynamics Kármán-Treftz mapping 3. Cylinder centre on the real negative ais z = b ( ζ + b) + ( ζ b) ( ζ + b) ( ζ b) y b ξ -b- b = b ε ( 1+ ε ) ε = 1,95, c = 3, 9b +

η Aerodynamics Kármán-Treftz mapping 4. Cylinder centre on the nd quadrant z = b ( ζ + b) + ( ζ b) ( ζ + b) ( ζ b) y b ξ -b- b = b ε ( 1+ ε ) ε = 1,95, c = 3, 9b +

Kármán-Treftz mapping 1. Cylinder centre at the origin 1 -C p 0.75 0.5 0.5 0-0.5 Etradorso, α=3 o Intradorso, α=3 o Etradorso, α=10 o Intradorso, α=10 o C p = p p r 1 ρ V -0.5-0.75-1 0 0.5 0.5 0.75 1 /c

Kármán-Treftz mapping. Cylinder centre on the imaginary ais 3 -C p.5 1.5 1 0.5 Etradorso, α=0 o Intradorso, α=0 o Etradorso, α=10 o Intradorso, α=10 o Etradorso, α=-10 o Intradorso, α=-10 o C p = p p r 1 ρ V 0-0.5-1 0 0.5 0.5 0.75 1 /c

Kármán-Treftz mapping 3. Cylinder centre on the real negative ais 4 3.5 3 Etradorso, α=0 o Intradorso, α=0 o Etradorso, α=10 o Intradorso, α=10 o -C p.5 1.5 1 0.5 C p = p p r 1 ρ V 0-0.5-1 0 0.5 0.5 0.75 1 /c

Kármán-Treftz mapping 4. Cylinder centre on the nd quadrant 4.5 -C p 4 3.5 3.5 1.5 1 Etradorso, α=0 o Intradorso, α=0 o Etradorso, α=10 o Intradorso, α=10 o Etradorso, α=-10 o Intradorso, α=-10 o C p = p p r 1 ρ V 0.5 0-0.5-1 0 0.5 0.5 0.75 1 /c

Generalization of the conformal mapping The Jouowsi transform is one of the mappings defined by z = ζ + n= 1 a n n ζ For the Jouowsi case: a1 = b an = 0 para n

Generalization of the conformal mapping z = ζ + n= 1 a n n ζ In general, the coefficients a n are comple Theoretically, the generalized mapping is able to obtain any airfoil from a cicular cylinder

Generalization of the conformal mapping z = ζ + n= 1 At small angles of attac, the lift coefficient, C l, of an airfoil in steady, irrotacional and incompressible flow is given by C For a Jouowsi airfoil a n n ζ d = π + at + c l 1 ( α β ) a t 0,77

M For an airfoil Q=0 Aerodynamics Pitching moment relative to the airfoil centre QΓ = ρ π + R Μ e 0 M 0 [ ] iα i πρu [ ] iα i πρ U = R e Μ Μ is the term proportional to z- of the comple velocity at large distances from the airfoil Assuming that the a 1 coefficient is given by a 1 iλ = b e Μ = a r r iα iα 1 V e a V e

Pitching moment relative to the airfoil centre For small values of α M c r 4πρ b V + ( α λ) Assuming a chord c 4b C M c = c M r ρ V 1 c π ( α + λ)

NACA airfoils The development of the NACA airfoils started in 1933 at the National Advisory Committee for Aeronautics (NACA), designated National Aeronautics and Space Administration (NASA) in nowadays NACA airfoils are obtained from the addition of a thicness distribution to a camber line. Several series of airfoils have been developed along the years

NACA airfoils 4 digits series, NACA ABCD A Maimum relative camber in percentage, f/c B Coordinate of maimum camber, m, given by 10 m /c CD Maimum relative thicness in percentage, d/c

NACA airfoils 4 digits series, NACA ABCD y - Thicness distribution d = 0, ( ) 3 4 0,969 0,16 0,3516 + 0,843 0.1015 - Camber lines parabola arcs matching at the point of maimum camber, located at m

NACA airfoils 4 digits series, NACA ABCD y - Thicness distribution d y = 0, m f = m 1 m ( ) 3 4 0,969 0,16 0,3516 + 0,843 0.1015 - Camber lines f ( ) m ( ) ( 1 )( 1+ ) m m m > m A f = 100 m B = 10

NACA airfoils 4 digits series - NACA 001 Symmetric airfoil with 1% of maimum relative thicness

NACA airfoils 4 digits series - NACA 441 Airfoil with maimum relative camber of 4% located at m =0.4c (40% from the leading edge) and maimum relative thicness of 1%

5 digits series, NACA ABCDE A Approimate value of BC Coordinate of maimum camber, m, given by 00 m /c DE Maimum relative thicness in percentage, d/c ( C l ) proj NACA airfoils 10 3 is the project lift coefficient ( C l ) proj

NACA airfoils 5 digits series, NACA ABCDE y - Thicness distribution identical to 4 digits series d = 0, ( ) 3 4 0,969 0,16 0,3516 + 0,843 0.1015 - Camber line polynomials with decreasing curvature from starting at the leading edge. Straight line to the right of the maimum camber point, located at m

NACA airfoils 5 digits series, NACA ABCDE - Thicness distribution identical to 4 digits series d y = 0, - Camber line 1 6 = 1 6 y ( ) 3 4 0,969 0,16 0,3516 + 0,843 0.1015 m 1 1 ( 3 m) ( 3 3m + m ) m 3 ( 1 ) m > m m

Aerodynamics 5 digits series, NACA ABCDE Designação m NACA airfoils - Camber line Values of 1, m e m for C l = ( ) 0, 3 proj m 1 10 0,05 0,0580 361,4 0 0,10 0,160 51,64 30 0,15 0,05 15,957 40 0,0 0.900 6,643 50 0,5 0,3910 3,30

NACA airfoils 5 digits series - NACA 3015 ( ) 3 Airfoil with C l = 0,, maimum camber located proj at m =0.15c (15% from the leading edge) and 15% of maimum relative thicness

NACA airfoils 6 digits series, NACA 6A,B-CDE, a=a o A Value of p 10c. p is the horizontal coordinate of the suction pea of the corresponding symmetric foil at zero degrees angle of attac. p is related to the thicness distribution of the airfoil

NACA airfoils 6 digits series, NACA 6A,B-CDE, a=a o B Range of C l values (multiplied by 10) above and below the value of C l for project conditions, which correspond to favourable or close to zero pressure gradients on the upper and lower sides of the airfoil. 10 ( C l ) proj ( C l ) proj C. is the lift coefficient for project conditions

NACA airfoils 6 digits series, NACA 6A,B-CDE, a=a o DE Maimum relative thicness in percentage a o Maimum horizontal coordinate of the region that presents an approimately constant loading (pressure difference between upper and lower sides of the airfoil). For >a o the loading decreases linearly. a o is related to the camber line NACA 66-09 (a=1)

y M V r BA α Pitching moment relative to the leading edge L α M c - Moment propagation M C M BA = M = C c M + L cos + C l ( α ) cos c M 1 c c + L ( α ) CM + Cl BA c c 1

y M V r BA α Pitching moment relative to the leading edge L α M c - For an airfoil at small values of α π C l π C M + c so π π CM ( β λ) + ( α + β ) BA π Cl CM γ + com γ = β λ BA 4 ( α + β ) ( α λ)

Aerodynamic centre L α y M - V r α M c Aerodynamic centre is the location that ehibits a pitching moment independent of the angle of attac, α Pitching moment relative to a point CM CM C c l c

Aerodynamic centre L α y M - V r α M c The location of the aerodynamic centre is given by dcm dc M = 0 = 0 so dα dcl dcm dc dc ca ca M M 1 c c c = 0 = = dc c c dc dα dc l l l dα

Aerodynamic centre L α y M - V r α M c For an airfoil at small values of α dcm c π dcl π dα dα so ca c π 1 = π 1 4

Aerodynamic centre L α y M - V r α M c The aerodynamic centre is approimately located 5% of the chord downstream of the leading edge The pitching moment relative to the aerodynamic centre is given by dcm dc π c CM = ca dα dα M c ( λ β ) = γ γ

Pressure centre L α y M - V r α M c Pressure centre is the point along the line that contains the chord that ehibits a zero pitching moment CM CM Cl = 0 c c

V r α Pressure centre L α y M c M cp C CM CM Cl = 0 = c c c C For an airfoil at small values of α π C l π C M + c so cp 1 1 γ 1 πγ = + = + = c 4 4 α + β 4 C ( α + β ) ( α λ) l M l c 1 4 + C - M C l ca