Incompressible Flow Over Airfoils
|
|
- Giles Cummings
- 5 years ago
- Views:
Transcription
1 Chapter 7 Incompressible Flow Over Airfoils Aerodynamics of wings: -D sectional characteristics of the airfoil; Finite wing characteristics (How to relate -D characteristics to 3-D characteristics) How to obtain -D characteristics? (a)experimental methods (NACA airfoils) (b)analytical methods (c)numerical methods Airfoil characteristics: () C l α () C d α (3) C m,c/4 α Characteristics of C l α curve: (a)α l= (b)a = dc l dα (c)c lmax 6
2 Copyright c 5 6
3 7.. VORTEX FILAMENT At low to moderate angles of attack C l α curve is linear. The flow moves slowly over the airfoil and is attached over most of the surface. At high angles of attack, the flow tends to separate from the top surface. C l,max occurs prior to stall C l,max is dependent on R e = ρvc µ C m,c/4 is independent of R e except for large α C d is dependent on R e The linear portion of the C l α curve is independent of R e and can be predicted using analytical methods. Theoretical/analytical methods to evaluate -D characteristics: Recall from potential flow that a spinning cylinder produces lift L = ρv Γ. clockwise): ψ = Γ lnr For a -D Vortex (spinning φ = Γ θ v θ = Γ r v r = 7. Vortex Filament Consider -D/point vortices of same strength duplicated in every plane parallel to the z-x plane along the y-axis from to. The flow is -D and is irrotational everywhere except the y-axis. y-axis is the straight vortex filament and may be defined as a line. Definition: A vortex filament is a straight or curved line in a fluid which coincides with the axis of rotation of successive fluid elements. Helmholt s vortex theorems:. The strength of a vortex filament is constant along its length. Proof: A vortex filament induces a velocity field that is irrotational at every point excluding the filament. Enclose a vortex filament with a sheath from which a slit has been removed. The vorticity at every point on the surface=. Evaluate the circulation for the sheath. Circulation = C V d s = S ( V ) d A. Sheath is irrotational. Thus V =, C V d s = or V d s = b V d s + d V d s + a V d s + V d s = a b c d However, a V d s + V d s = b d Copyright c 5 6
4 7.. VORTEX FILAMENT as it constitutes the integral across the slit. Thus, b d V d s + V d s = a c b a d V d s = c V d s = d V d s = Γ. A vortex filament cannot end in a fluid; it must extend to the boundaries of the fluid or form a closed path. 3. In the absence of rotational external flow, a fluid that is irrotational remains irrotational. 4. In the absence of rotational extenal force, if the circulation around a path enclosing a definite group of particles is initially zero, it will remain zero. 5. In the absence of rotational extenal force, the circulation around a path that encloses a tagged group of elements is invariant. Copyright c 5 63
5 7.. VORTEX FILAMENT Vortex Sheet or vortex surface An infinite number of straight vortex filaments placed side by side form a vortex sheet. Each vortex filament has an infinitesimal strength γ(s). ds v θ v r r θ vortex sheet γ(s) is the strength of vortex sheet per unit length along s. v θ = Γ r for -D (point vortex). A small portion of the vortex sheet of strength γds induces an infinitesimally small velocity dv at a field point P (r, θ). So, v θ vortexfilament = γds r... dv P = γds r. Circulation Γ around a point vortex is equal to the strength of the vortex. Similarly, the circulation around the vortex sheet is the sum of the strengths of the elemental vortices. Therefore, the circulation Γ for a finite length from point a to point b on the vortex sheet is given by: Γ = b a γ(s)ds Across a vortex sheet, there is a discontinuous change in the tangential component of velocity and the normal component of velocity is preserved. z u w n w s u x Copyright c 5 64
6 7.. VORTEX FILAMENT As n, we get Γ = Γ = v d l Box v d l = [w n u s w n + u s] γ s = (u u ) s + (w w ) n γ s = (u u ) s or γ = (u u ) γ = (u u ) states that the local jump in tangential velocity across the vortex sheet is equal to the local sheet strength. Kutta Condition: For a given airfoil at a given angle of attack, the value of Γ around the airfoil is such that the flow leaves the trailing edge smoothly. The Kutta condition tells us how to find Γ; it is based on experimental observation. A body with finite angle TE in relative motion through a fluid will create about itself a circulation of sufficient strength to hold the rear stagnation point at the TE. If the TE has a zero angle, the Kutta Condition requires that the velocity of fluid leaving upper and lower surfaces at the TE be equal and non-zero. A body with a finite TE angle will have crossing streamlines at the TE unless the TE is a stagnation point. The Kutta Condition eliminates the crossing streamlines. Consider the TE as a vortex sheet: γ(t E) = V u V l If the TE has a finite angle V u = V l = because TE is a stagnation point. Or γ(t E) =. If the TE is a cusp (zero angle), V u = V l and hence γ(t E) = V u V l =. Thus, Kutta Condition expressed mathematically in terms of vortex strength is γ(t E) =. 7.. Bound Vortex and Starting Vortex The question might arise: Does a real airfoil flying in a real fluid give rise to a circulation about itself? The answer is yes. When a wing section with a sharp T.E is put into motion, the fluid has a tendency to go around the sharp T.E from the lower to the upper surface. As the airfoil moves along vortices are shed from the T.E which form a vortex sheet. Copyright c 5 65
7 7.. VORTEX FILAMENT Helmholtz s theorem: If Γ = originally in a flow it remains zero. Kelvin s theorem: Circulation around a closed curve formed by a set of continuous fluid elements remains constant as DΓ the fluid elements move through the flow, Dt =. Substantial derivative gives the time rate of change following a given fluid element. The circulation about the airfoil is replaced by a vortex of equal strength. This is the bound vortex since it remains bound to the airfoil. A true vortex remains attaced to the same fluid particles and moves with the general flow. Thus, as far as resultant forces are concerned, a bound vortex of proper strength in a uniform flow is equivalent to a body with circulation in a uniform flow. Both the theorems are satisfied by the starting vortex and bound vortex system. In the beginning, Γ = when the flow is started within the contour C. When the flow over the airfoil is developed, Γ within C is still zero which includes the starting vortex Γ 3 and the bound vortex Γ 4 which are equal and opposite. Copyright c 5 66
8 7.. VORTEX FILAMENT Copyright c 5 67
9 7.. FUNDAMENTAL EQUATION OF THIN AIRFOIL THEORY 7. Fundamental Equation of Thin Airfoil Theory Principle: Mean camber line is a streamline of the flow. Velocity induced by a -D vortex is V = v θ ê θ = Γ r êθ where Γ is the strength of the -D vortex. Similarly the velocity induced by the vortex sheet of infinitesimal length ds is given by d v P = γ(s)ds r êθ To force the mean camber line to be a streamline, the sum of all velocity components normal to the mcl must be equal to zero. Consider the flow induced by an elemental vortex sheet ds at a point P on the vortex sheet. It is perpendicular to the line connecting the center of ds to the point P given by d v P = γ(s)ds r êθ Thus dw P the velocity normal to the mcl is: dw P γ(s) cos βds = dv P cos β = r where β is the angle made by dv p to the normal at P, and r is the distance from the center of ds to the point P. The induced velocity due to the vortex sheet representing the entire mcl is given by; w P (s) = T E γ(s) cos β ds r Now determine the component of the freestream velocity normal to the mcl. Copyright c 5 68
10 7.3. FLAT PLATE AT AN ANG OF ATTACK V,n = V sin(α + ɛ) where α is the angle of attack and ɛ is the angle made by the tangent at point P to the x-axis. The slope of the tangent line at point P is given by: or dz = tan( ɛ) = tan ɛ dx ɛ = tan ( dz dx ) V,n = V (sin α + tan ( dz dx )) In order that the mcl is a streamline. w P (s) + V,n =, or T E γ(s) cos β ds + V (sin α + tan ( dz r dx )) = within thin airfoil theory approximation s x, ds dx, cos β = and r (x x), where x varies from to c, and x refers to the point P. After changing these variables and making the small angle approximation for sin and tan, and upon rearrangement we get: c γ(x) x x dx = V (α dz dx ) 7.3 Flat Plate at an Angle of Attack The following analysis is an exact solution to the flat plate or an approximate solution to the symmetric airfoil. The mean camber line becomes the chord and hence: dz dx = c γ(x) (x x) dx = V α In order to facilitate analytic solution, we do a variable transformation such that: x = c ( cos θ) x = c ( cos θ ) θ = at and θ = at TE and θ increases in CW, dx = c sin θdθ γ(θ) c sin θ c [( cos θ ) ( cos θ)] dθ = V α Copyright c 5 69
11 7.4. -D LIFT COEFFICIENT FOR A THIN/SYMMETRICAL AIRFOIL Here we simply state a rigorous solution for γ(θ) as: γ(θ) sin θ (cos θ cos θ ) dθ = V α γ(θ) = αv + cosθ sin θ We can verify this solution by substitution as follows: γ(θ) sin θ (cos θ cos θ ) dθ = V α We now use the following result to evaluate the above integral. V α Thus, it satisfies the equation: ( + cos θ) (cos θ cos θ ) dθ = V α (cos nθ) (cos θ cos θ ) dθ = sin nθ sin θ ( + cos θ) (cos θ cos θ ) dθ dθ + cos θ cos θ cos θ dθ (7.) cos θ cos θ = V α ( + ) = V α (7.) γ(θ) sin θ (cos θ cos θ ) dθ = V α In addition, the solution for γ also satisfies the Kutta condition. When θ =, By using L Hospital rule, we get Thus it satifies the Kutta condition. γ() = V α γ() = V α sin cos = 7.4 -D lift coefficient for a thin/symmetrical airfoil Where s is along the mcl. By using thin airfoil approximation: L L ρv = ρ V Γ = ρv T E T E γ(x)dx = ρv γ(s)ds γ(x)dx Copyright c 5 7
12 7.4. -D LIFT COEFFICIENT FOR A THIN/SYMMETRICAL AIRFOIL Using the transformation x = c ( cosθ) Substituting the solution: L ρ V γ(x)dx = ρ V c γ(θ) sin θdθ L γ(θ) = αv ( + cos θ) sin θ ρ V c L αρ V c L V α( + cos θ)dθ cαρ V ( + cos θ)dθ dc l dα L α( ρ V )(c ) L αq S or C l = L q = α, and dc l S dα = = shows that lift curve is linearly proportional to the angle of attack Calculation of Moment Coefficient M = = = = ρ V = ρ V x(dl ) (7.3) x(ρ V dγ) (7.4) x(ρ V γ(x)dx) (7.5) = αρ V c 4 = αρ V c ( = q c α ) γ(x)xdx (7.6) αv ( + cos θ) sin θ ( ) c c ( cos θ) sin θdθ (7.7) ( cos θ)dθ (7.8) (7.9) (7.) Copyright c 5 7
13 7.5. THIN AIRFOIL THEORY FOR CAMBERED AIRFOIL C m, = M q Sc = M q c = α C l = α C m, = C l 4 M = M c/4 L c/4 C m, = C m,c/4 C l /4 C m, = C l /4 C m,c/4 is equal to zero for all values of α. C m,c/4 = c/4 is the aerodynamic center. Aerodynamic center is that point on an airfoil where moments are independant of angle of attack. 7.5 Thin Airfoil Theory for Cambered Airfoil γ(x)dx x x = V (α dz dx ) where dz dx is the slope of mcl at x. For symmetric airfoil, mcl is a straight line and hence dz dx = everywhere. On the other hand, for a cambered airfoil dz dx varies from point to point. As before, we do a variable transformation given by: Equation (A) becomes: x = c ( cos θ) dx = c sin θdθ (A) c γ(θ) sin θdθ (cos θ cos θ ) = V (α dz dx ) (B) has a solution of the form: Sub. solution in equation (B) γ(θ) = V [A ( + cos θ sin θ ) + A n sin nθ] n= A ( + cos θ) (cos θ cos θ ) dθ + n= A n sin nθ sin θ dz dθ = (α ) (cos θ cos θ ) dx x Using the integral cos nθ (cos θ cos θ ) dθ = sin nθ sin θ Copyright c 5 7
14 7.5. THIN AIRFOIL THEORY FOR CAMBERED AIRFOIL The first term becomes: Using the integral The second term becomes: = = A A ( + cos θ) (cos θ cos θ ) dθ n= A (cos θ cos θ ) dθ + A cos θ (cos θ cos θ ) dθ sin nθ sin θ (cos θ cos θ ) dθ = cos nθ A n sin nθ sin θ (cos θ cos θ ) dθ = A n cos nθ n= Therefore Equation (B) becomes: A n= A n cos nθ = α dz dx x Upon rearrangement the slope at a point P on the mcl is given by: dz dx = (α A ) + A n cos nθ n= From Fourier serious Where, f(θ) = B + B = B n = n =,,..., B n cos nθ n= f(θ)dθ (α A ) = B = A n = B n f(θ) cos nθdθ ( dz dx )dθ Evaluation of Γ From thin airfoil theroy Γ = γ(x)dx = c γ(θ) sin θdθ γ(θ) = V [ A ( + cos θ) + A n sin nθ] sin θ n= Γ = c γ(θ) sin θdθ = c V A ( + cos θ)dθ + c = cv [A (θ + sin θ) + A n sin nθ sin θdθ] n= V n= A n sin nθ sin θdθ Copyright c 5 73
15 7.5. THIN AIRFOIL THEORY FOR CAMBERED AIRFOIL Using A n sin nθ sin θdθ n= = for n = = for n Γ = cv [A + A ] L = ρ V Γ = ρv c[a + A ] C l = L q s = L q c = [A + A ] C l is normalized by the α as seen by the chord connecting the and TE of the mcl. c is the chord connecting the and TE of the mcl. C l = A + A ( = α [ = α + as is the case for symmetric airfoil. Also, Determination of moment coefficient ( ) dz dx ( dz dx ) ( dθ + ) (cos θ )dθ dc l dα = C l = dc l dα (α α L=) = (α α L= )... α L= = ] (cos θ ) dz dx dθ M = ρ V xγ(x)dx C m, = M q sc = V c As before we do a variable transformation from x to θ. Thus, x = c ( cos θ) xγ(x)dx ( ) ) dz cos θdθ dx dx = c sin θdθ γ(θ) = V [ A ( + cos θ) + A n sin nθ] sin θ n= C m, = A ( cos θ)dθ A n ( cos θ) sin θ sin nθdθ n= Copyright c 5 74
16 7.5. THIN AIRFOIL THEORY FOR CAMBERED AIRFOIL Using the following definite integrals: cos θdθ = sin θdθ = sin θ sin nθdθ = for n = = for n =,..., sin θ cos θ sin nθdθ = for n = = 4 for n = = for n = 3,..., C m, = A dθ + A cos θdθ (A n sin θ sin nθdθ A n n= = A + A A + A 4 = A A + A 4 C l = A + A C m, = [ ] A + A A ) sin θ cos θ sin nθdθ [ Cl = 4 + ] 4 (A A ) M = M c 4 L c 4 C m, = c m, c 4 C l 4 = C m, c 4 = 4 (A A ) C m, = x cp C l c [ Cl 4 + ] 4 (A A ) x cp = [ c 4 + c (A A )] = c [c + (A A )] 4C l 4 C l Relationship between pressure on mcl and γ () () (3) dl = (p l p u ) cos η(ds ) L = L = T E T E (p l p u ) cos ηds (A) ρv γ(s)ds Equating (A) and (B) and setting cos η =, we get T E (p l p u )ds = T E (B) ρv γ(s)ds Copyright c 5 75
17 7.5. THIN AIRFOIL THEORY FOR CAMBERED AIRFOIL or Using Bernoullis equation: (p l p u ) = ρv γ(s) () From vortex sheet theory: From (), () and (3) p l + ρ(u l + w ) = p u + ρ(u u + w ) p l p u = ρ (u u + u l )(u u u l ) () u u u l = γ(s) (3) V = u u + u l i.e., within thin airfoil approximation, the average of top and bottom surface velocities at any point on the mcl is equal to the freestream velocity. c p,l c p,u = (p l p ) q (p u p ) q = p l p u q = ρ(u u + u l )(u u u l ) q = γ(s)(v ) v = γ(s) V Copyright c 5 76
Lifting Airfoils in Incompressible Irrotational Flow. AA210b Lecture 3 January 13, AA210b - Fundamentals of Compressible Flow II 1
Lifting Airfoils in Incompressible Irrotational Flow AA21b Lecture 3 January 13, 28 AA21b - Fundamentals of Compressible Flow II 1 Governing Equations For an incompressible fluid, the continuity equation
More information1. Fluid Dynamics Around Airfoils
1. Fluid Dynamics Around Airfoils Two-dimensional flow around a streamlined shape Foces on an airfoil Distribution of pressue coefficient over an airfoil The variation of the lift coefficient with the
More informationCopyright 2007 N. Komerath. Other rights may be specified with individual items. All rights reserved.
Low Speed Aerodynamics Notes 5: Potential ti Flow Method Objective: Get a method to describe flow velocity fields and relate them to surface shapes consistently. Strategy: Describe the flow field as the
More informationGiven a stream function for a cylinder in a uniform flow with circulation: a) Sketch the flow pattern in terms of streamlines.
Question Given a stream function for a cylinder in a uniform flow with circulation: R Γ r ψ = U r sinθ + ln r π R a) Sketch the flow pattern in terms of streamlines. b) Derive an expression for the angular
More informationASTR 320: Solutions to Problem Set 3
ASTR 30: Solutions to Problem Set 3 Problem : The Venturi Meter The venturi meter is used to measure the flow speed in a pipe. An example is shown in Fig., where the venturi meter (indicated by the dashed
More informationFluid Mechanics Prof. T. I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay
Fluid Mechanics Prof. T. I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay Lecture No. # 35 Boundary Layer Theory and Applications Welcome back to the video course on fluid
More informationPEMP ACD2505. M.S. Ramaiah School of Advanced Studies, Bengaluru
Two-Dimensional Potential Flow Session delivered by: Prof. M. D. Deshpande 1 Session Objectives -- At the end of this session the delegate would have understood PEMP The potential theory and its application
More informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
More informationContinuum Mechanics Lecture 7 Theory of 2D potential flows
Continuum Mechanics ecture 7 Theory of 2D potential flows Prof. http://www.itp.uzh.ch/~teyssier Outline - velocity potential and stream function - complex potential - elementary solutions - flow past a
More informationGiven the water behaves as shown above, which direction will the cylinder rotate?
water stream fixed but free to rotate Given the water behaves as shown above, which direction will the cylinder rotate? ) Clockwise 2) Counter-clockwise 3) Not enough information F y U 0 U F x V=0 V=0
More informationTHE VORTEX PANEL METHOD
THE VORTEX PANEL METHOD y j m α V 4 3 2 panel 1 a) Approimate the contour of the airfoil by an inscribed polygon with m sides, called panels. Number the panels clockwise with panel #1 starting on the lower
More informationAerodynamics. Airfoils Kármán-Treftz mapping. The exponent k defines the internal angle of the trailing edge, τ, from. ( k)
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
More informationAerodynamics. High-Lift Devices
High-Lift Devices Devices to increase the lift coefficient by geometry changes (camber and/or chord) and/or boundary-layer control (avoid flow separation - Flaps, trailing edge devices - Slats, leading
More informationOffshore Hydromechanics Module 1
Offshore Hydromechanics Module 1 Dr. ir. Pepijn de Jong 4. Potential Flows part 2 Introduction Topics of Module 1 Problems of interest Chapter 1 Hydrostatics Chapter 2 Floating stability Chapter 2 Constant
More informationComputing potential flows around Joukowski airfoils using FFTs
AB CD EF GH Computing potential flows around Joukowski airfoils using FFTs Frank Brontsema Institute for Mathematics and Computing Science AB CD EF GH Bachelor thesis Computing potential flows around
More informationINSTITUTE OF AERONAUTICAL ENGINEERING (Autonomous) Dundigal, Hyderabad
INSTITUTE OF AERONAUTICAL ENGINEERING (Autonomous) Dundigal, Hyderabad - 500 043 AERONAUTICAL ENGINEERING TUTORIAL QUESTION BANK Course Name : LOW SPEED AERODYNAMICS Course Code : AAE004 Regulation : IARE
More informationExercise 9: Model of a Submarine
Fluid Mechanics, SG4, HT3 October 4, 3 Eample : Submarine Eercise 9: Model of a Submarine The flow around a submarine moving at a velocity V can be described by the flow caused by a source and a sink with
More informationAn Internet Book on Fluid Dynamics. Joukowski Airfoils
An Internet Book on Fluid Dynamics Joukowski Airfoils One of the more important potential flow results obtained using conformal mapping are the solutions of the potential flows past a family of airfoil
More informationLab Reports Due on Monday, 11/24/2014
AE 3610 Aerodynamics I Wind Tunnel Laboratory: Lab 4 - Pressure distribution on the surface of a rotating circular cylinder Lab Reports Due on Monday, 11/24/2014 Objective In this lab, students will be
More informationIntroduction to Aeronautics
Introduction to Aeronautics ARO 101 Sections 03 & 04 Sep 30, 2015 thru Dec 9, 2015 Instructor: Raymond A. Hudson Week #8 Lecture Material 1 Topics For Week #8 Airfoil Geometry & Nomenclature Identify the
More informationBernoulli's equation: 1 p h t p t. near the far from plate the plate. p u
UNSTEADY FLOW let s re-visit the Kutta condition when the flow is unsteady: p u p l Bernoulli's equation: 2 φ v 1 + + = () = + p h t p t 2 2 near the far from plate the plate as a statement of the Kutta
More informationFinal Review of AerE 243 Class
Final Review of AeE 4 Class Content of Aeodynamics I I Chapte : Review of Multivaiable Calculus Chapte : Review of Vectos Chapte : Review of Fluid Mechanics Chapte 4: Consevation Equations Chapte 5: Simplifications
More informationExercise 9, Ex. 6.3 ( submarine )
Exercise 9, Ex. 6.3 ( submarine The flow around a submarine moving at at velocity V can be described by the flow caused by a source and a sink with strength Q at a distance a from each other. V x Submarine
More informationAll that begins... peace be upon you
All that begins... peace be upon you Faculty of Mechanical Engineering Department of Thermo Fluids SKMM 2323 Mechanics of Fluids 2 «An excerpt (mostly) from White (2011)» ibn Abdullah May 2017 Outline
More informationAE301 Aerodynamics I UNIT B: Theory of Aerodynamics
AE301 Aerodynamics I UNIT B: Theory of Aerodynamics ROAD MAP... B-1: Mathematics for Aerodynamics B-: Flow Field Representations B-3: Potential Flow Analysis B-4: Applications of Potential Flow Analysis
More informationHomework Two. Abstract: Liu. Solutions for Homework Problems Two: (50 points total). Collected by Junyu
Homework Two Abstract: Liu. Solutions for Homework Problems Two: (50 points total). Collected by Junyu Contents 1 BT Problem 13.15 (8 points) (by Nick Hunter-Jones) 1 2 BT Problem 14.2 (12 points: 3+3+3+3)
More informationHigh Speed Aerodynamics. Copyright 2009 Narayanan Komerath
Welcome to High Speed Aerodynamics 1 Lift, drag and pitching moment? Linearized Potential Flow Transformations Compressible Boundary Layer WHAT IS HIGH SPEED AERODYNAMICS? Airfoil section? Thin airfoil
More informationSPC Aerodynamics Course Assignment Due Date Monday 28 May 2018 at 11:30
SPC 307 - Aerodynamics Course Assignment Due Date Monday 28 May 2018 at 11:30 1. The maximum velocity at which an aircraft can cruise occurs when the thrust available with the engines operating with the
More informationThe Divergence Theorem Stokes Theorem Applications of Vector Calculus. Calculus. Vector Calculus (III)
Calculus Vector Calculus (III) Outline 1 The Divergence Theorem 2 Stokes Theorem 3 Applications of Vector Calculus The Divergence Theorem (I) Recall that at the end of section 12.5, we had rewritten Green
More informationWeek 2 Notes, Math 865, Tanveer
Week 2 Notes, Math 865, Tanveer 1. Incompressible constant density equations in different forms Recall we derived the Navier-Stokes equation for incompressible constant density, i.e. homogeneous flows:
More informationMAE SUMMER 2015 HOMEWORK 1 SOLUTION
MAE 04 - SUMMER 205 HOMEWORK SOLUTION Problem :. Plot the drag coefficient, c D, the aerodynamic efficiency, AE, and the center of pressure, x cp, of the airplane as a function of the angle of attack.
More informationFlight Vehicle Terminology
Flight Vehicle Terminology 1.0 Axes Systems There are 3 axes systems which can be used in Aeronautics, Aerodynamics & Flight Mechanics: Ground Axes G(x 0, y 0, z 0 ) Body Axes G(x, y, z) Aerodynamic Axes
More informationSome Basic Plane Potential Flows
Some Basic Plane Potential Flows Uniform Stream in the x Direction A uniform stream V = iu, as in the Fig. (Solid lines are streamlines and dashed lines are potential lines), possesses both a stream function
More informationAirfoils and Wings. Eugene M. Cliff
Airfoils and Wings Eugene M. Cliff 1 Introduction The primary purpose of these notes is to supplement the text material related to aerodynamic forces. We are mainly interested in the forces on wings and
More informationMAE 101A. Homework 7 - Solutions 3/12/2018
MAE 101A Homework 7 - Solutions 3/12/2018 Munson 6.31: The stream function for a two-dimensional, nonviscous, incompressible flow field is given by the expression ψ = 2(x y) where the stream function has
More information3.1 Definition Physical meaning Streamfunction and vorticity The Rankine vortex Circulation...
Chapter 3 Vorticity Contents 3.1 Definition.................................. 19 3.2 Physical meaning............................. 19 3.3 Streamfunction and vorticity...................... 21 3.4 The Rankine
More informationPART II: 2D Potential Flow
AERO301:Spring2011 II(a):EulerEqn.& ω = 0 Page1 PART II: 2D Potential Flow II(a): Euler s Equation& Irrotational Flow We have now completed our tour through the fundamental conservation laws that apply
More informationSimplifications to Conservation Equations
Chater 5 Simlifications to Conservation Equations 5.1 Steady Flow If fluid roerties at a oint in a field do not change with time, then they are a function of sace only. They are reresented by: ϕ = ϕq 1,
More informationPart A Fluid Dynamics & Waves Draft date: 17 February Conformal mapping
Part A Fluid Dynamics & Waves Draft date: 17 February 4 3 1 3 Conformal mapping 3.1 Wedges and channels 3.1.1 The basic idea Suppose we wish to find the flow due to some given singularities (sources, vortices,
More informationNumerical Modeling of Lifting Flows in the Presence of a Free Surface
University of New Orleans ScholarWorks@UNO University of New Orleans Theses and Dissertations Dissertations and Theses Spring 5-18-2012 Numerical Modeling of Lifting Flows in the Presence of a Free Surface
More informationFluid mechanics, topology, and complex analysis
Fluid mechanics, topology, and complex analysis Takehito Yokoyama Department of Physics, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo 152-8551, Japan (Dated: April 30, 2013 OMPLEX POTENTIAL
More informationInviscid & Incompressible flow
< 3.1. Introduction and Road Map > Basic aspects of inviscid, incompressible flow Bernoulli s Equation Laplaces s Equation Some Elementary flows Some simple applications 1.Venturi 2. Low-speed wind tunnel
More informationTHE AERODYNAMICS OF A RING AIRFOIL*
136 THE AERODYNAMICS OF A RING AIRFOIL* H. J. STEWART California Institute of Technology Abstract. The downwash required to produce a given vorticity distribution is computed for a ring airfoil and the
More informationChapter 9 Flow over Immersed Bodies
57:00 Mechanics of Fluids and Transport Processes Chapter 9 Professor Fred Stern Fall 009 1 Chapter 9 Flow over Immersed Bodies Fluid flows are broadly categorized: 1. Internal flows such as ducts/pipes,
More informationMestrado Integrado em Engenharia Mecânica Aerodynamics 1 st Semester 2012/13
Mestrado Integrado em Engenharia Mecânica Aerodynamics 1 st Semester 212/13 Exam 2ª época, 2 February 213 Name : Time : 8: Number: Duration : 3 hours 1 st Part : No textbooks/notes allowed 2 nd Part :
More informationContinuum Mechanics Lecture 5 Ideal fluids
Continuum Mechanics Lecture 5 Ideal fluids Prof. http://www.itp.uzh.ch/~teyssier Outline - Helmholtz decomposition - Divergence and curl theorem - Kelvin s circulation theorem - The vorticity equation
More information18.325: Vortex Dynamics
8.35: Vortex Dynamics Problem Sheet. Fluid is barotropic which means p = p(. The Euler equation, in presence of a conservative body force, is Du Dt = p χ. This can be written, on use of a vector identity,
More informationContents. MATH 32B-2 (18W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables. 1 Multiple Integrals 3. 2 Vector Fields 9
MATH 32B-2 (8W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables Contents Multiple Integrals 3 2 Vector Fields 9 3 Line and Surface Integrals 5 4 The Classical Integral Theorems 9 MATH 32B-2 (8W)
More informationThe vorticity field. A dust devil
The vorticity field The vector ω = u curl u is twice the local angular velocity in the flow, and is called the vorticity of the flow (from Latin for a whirlpool). Vortex lines are everywhere in the direction
More informationSteady waves in compressible flow
Chapter Steady waves in compressible flow. Oblique shock waves Figure. shows an oblique shock wave produced when a supersonic flow is deflected by an angle. Figure.: Flow geometry near a plane oblique
More informationSyllabus for AE3610, Aerodynamics I
Syllabus for AE3610, Aerodynamics I Current Catalog Data: AE 3610 Aerodynamics I Credit: 4 hours A study of incompressible aerodynamics of flight vehicles with emphasis on combined application of theory
More informationMATH 566: FINAL PROJECT
MATH 566: FINAL PROJECT December, 010 JAN E.M. FEYS Complex analysis is a standard part of any math curriculum. Less known is the intense connection between the pure complex analysis and fluid dynamics.
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationAerodynamic Performance 1. Figure 1: Flowfield of a Wind Turbine and Actuator disc. Table 1: Properties of the actuator disk.
Aerodynamic Performance 1 1 Momentum Theory Figure 1: Flowfield of a Wind Turbine and Actuator disc. Table 1: Properties of the actuator disk. 1. The flow is perfect fluid, steady, and incompressible.
More informationWater is sloshing back and forth between two infinite vertical walls separated by a distance L: h(x,t) Water L
ME9a. SOLUTIONS. Nov., 29. Due Nov. 7 PROBLEM 2 Water is sloshing back and forth between two infinite vertical walls separated by a distance L: y Surface Water L h(x,t x Tank The flow is assumed to be
More informationMath 575-Lecture Failure of ideal fluid; Vanishing viscosity. 1.1 Drawbacks of ideal fluids. 1.2 vanishing viscosity
Math 575-Lecture 12 In this lecture, we investigate why the ideal fluid is not suitable sometimes; try to explain why the negative circulation appears in the airfoil and introduce the vortical wake to
More information2.00AJ / 16.00AJ Exploring Sea, Space, & Earth: Fundamentals of Engineering Design Spring 2009
MIT OpenCourseWare http://ocw.mit.edu 2.00AJ / 16.00AJ Exploring Sea, Space, & Earth: Fundamentals of Engineering Design Spring 2009 For information about citing these materials or our Terms of Use, visit:
More informationGeneral Solution of the Incompressible, Potential Flow Equations
CHAPTER 3 General Solution of the Incompressible, Potential Flow Equations Developing the basic methodology for obtaining the elementary solutions to potential flow problem. Linear nature of the potential
More informationMasters in Mechanical Engineering Aerodynamics 1 st Semester 2015/16
Masters in Mechanical Engineering Aerodynamics st Semester 05/6 Exam st season, 8 January 06 Name : Time : 8:30 Number: Duration : 3 hours st Part : No textbooks/notes allowed nd Part : Textbooks allowed
More informationII. Ideal fluid flow
II. Ideal fluid flow Ideal fluids are Inviscid Incompressible The only ones decently understood mathematically Governing equations u=0 Continuity u 1 +( u ) u= ρ p+ f t Euler Boundary conditions Normal
More informationA Discrete-Vortex Method for Studying the Wing Rock of Delta Wings
NASA/TP-2002-211965 A Discrete-Vortex Method for Studying the Wing Rock of Delta Wings Thomas G. Gainer Langley Research Center, Hampton, Virginia December 2002 The NASA STI Program Office... in Profile
More informationIntroduction to Aerospace Engineering
Introduction to Aerosace Engineering Lecture slides hallenge the future Introduction to Aerosace Engineering Aerodynamics & Prof. H. Bijl ir. N. Timmer &. Airfoils and finite wings Anderson 5.9 end of
More informationUnsteady Aerodynamics and Vortex-sheet Formation of A Two-dimensional Airfoil arxiv: v1 [physics.flu-dyn] 17 Nov 2016
arxiv version Unsteady Aerodynamics and Vortex-sheet Formation of A Two-dimensional Airfoil arxiv:6.5729v [physics.flu-dyn] 7 Nov 26 X. X I A A N D K. M O H S E N I,2 Department of Mechanical and Aerospace
More informationChapter 5 Wing design - selection of wing parameters 2 Lecture 20 Topics
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
More informationMethod of Images
. - Marine Hdrodnamics, Spring 5 Lecture 11. - Marine Hdrodnamics Lecture 11 3.11 - Method of Images m Potential for single source: φ = ln + π m ( ) Potential for source near a wall: φ = m ln +( ) +ln
More informationWings and Bodies in Compressible Flows
Wings and Bodies in Compressible Flows Prandtl-Glauert-Goethert Transformation Potential equation: 1 If we choose and Laplace eqn. The transformation has stretched the x co-ordinate by 2 Values of at corresponding
More informationFUNDAMENTALS OF AERODYNAMICS
*A \ FUNDAMENTALS OF AERODYNAMICS Second Edition John D. Anderson, Jr. Professor of Aerospace Engineering University of Maryland H ' McGraw-Hill, Inc. New York St. Louis San Francisco Auckland Bogota Caracas
More information1. Introduction - Tutorials
1. Introduction - Tutorials 1.1 Physical properties of fluids Give the following fluid and physical properties(at 20 Celsius and standard pressure) with a 4-digit accuracy. Air density : Water density
More informationAerodynamics. Basic Aerodynamics. Continuity equation (mass conserved) Some thermodynamics. Energy equation (energy conserved)
Flow with no friction (inviscid) Aerodynamics Basic Aerodynamics Continuity equation (mass conserved) Flow with friction (viscous) Momentum equation (F = ma) 1. Euler s equation 2. Bernoulli s equation
More informationThin airfoil theory. Chapter Compressible potential flow The full potential equation
hapter 4 Thin airfoil theory 4. ompressible potential flow 4.. The full potential equation In compressible flow, both the lift and drag of a thin airfoil can be determined to a reasonable level of accuracy
More informationLecture Notes Fluid Mechanics of Turbomachines II
Lecture Notes Fluid Mechanics of Turbomachines II N.P. Kruyt 999-2009 N.P. Kruyt Turbomachinery Laboratory Engineering Fluid Dynamics Department of Mechanical Engineering University of Twente The Netherlands
More informationA Study of Transonic Flow and Airfoils. Presented by: Huiliang Lui 30 th April 2007
A Study of Transonic Flow and Airfoils Presented by: Huiliang Lui 3 th April 7 Contents Background Aims Theory Conservation Laws Irrotational Flow Self-Similarity Characteristics Numerical Modeling Conclusion
More information2.25 Advanced Fluid Mechanics Fall 2013
.5 Advanced Fluid Mechanics Fall 013 Solution to Problem 1-Final Exam- Fall 013 r j g u j ρ, µ,σ,u j u r 1 h(r) r = R Figure 1: Viscous Savart Sheet. Image courtesy: Villermaux et. al. [1]. This kind of
More informationFLOW SEPARATION. Aerodynamics Bridge-Pier Design Combustion Chambers Human Blood Flow Building Design Etc.
FLOW SEPARATION Aerodynamics Bridge-Pier Design Combustion Chambers Human Blood Flow Building Design Etc. (Form Drag, Pressure Distribution, Forces and Moments, Heat And Mass Transfer, Vortex Shedding)
More informationa) Derive general expressions for the stream function Ψ and the velocity potential function φ for the combined flow. [12 Marks]
Question 1 A horizontal irrotational flow system results from the combination of a free vortex, rotating anticlockwise, of strength K=πv θ r, located with its centre at the origin, with a uniform flow
More informationConsider a wing of finite span with an elliptic circulation distribution:
Question 1 (a) onsider a wing of finite span with an elliptic circulation distribution: Γ( y) Γo y + b = 1, - s y s where s=b/ denotes the wing semi-span. Use this equation, in conjunction with the Kutta-Joukowsky
More informationDetailed Outline, M E 521: Foundations of Fluid Mechanics I
Detailed Outline, M E 521: Foundations of Fluid Mechanics I I. Introduction and Review A. Notation 1. Vectors 2. Second-order tensors 3. Volume vs. velocity 4. Del operator B. Chapter 1: Review of Basic
More informationReview of fluid dynamics
Chapter 2 Review of fluid dynamics 2.1 Preliminaries ome basic concepts: A fluid is a substance that deforms continuously under stress. A Material olume is a tagged region that moves with the fluid. Hence
More informationOffshore Hydromechanics
Offshore Hydromechanics Module 1 : Hydrostatics Constant Flows Surface Waves OE4620 Offshore Hydromechanics Ir. W.E. de Vries Offshore Engineering Today First hour: Schedule for remainder of hydromechanics
More informationLift Enhancement on Unconventional Airfoils
Lift Enhancement on nconventional Airfoils W.W.H. Yeung School of Mechanical and Aerospace Engineering Nanang Technological niversit, North Spine (N3), Level 2 50 Nanang Avenue, Singapore 639798 mwheung@ntu.edu.sg
More informationThe purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and engineering.
Lecture 16 Applications of Conformal Mapping MATH-GA 451.001 Complex Variables The purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and
More informationMATH 228: Calculus III (FALL 2016) Sample Problems for FINAL EXAM SOLUTIONS
MATH 228: Calculus III (FALL 216) Sample Problems for FINAL EXAM SOLUTIONS MATH 228 Page 2 Problem 1. (2pts) Evaluate the line integral C xy dx + (x + y) dy along the parabola y x2 from ( 1, 1) to (2,
More informationAn alternative approach to integral equation method based on Treftz solution for inviscid incompressible flow
An alternative approach to integral equation method based on Treftz solution for inviscid incompressible flow Antonio C. Mendes, Jose C. Pascoa Universidade da Beira Interior, Laboratory of Fluid Mechanics,
More informationChapter 2 Dynamics of Perfect Fluids
hapter 2 Dynamics of Perfect Fluids As discussed in the previous chapter, the viscosity of fluids induces tangential stresses in relatively moving fluids. A familiar example is water being poured into
More informationMAE 222 Mechanics of Fluids Final Exam with Answers January 13, Give succinct answers to the following word questions.
MAE 222 Mechanics of Fluids Final Exam with Answers January 13, 1994 Closed Book Only, three hours: 1:30PM to 4:30PM 1. Give succinct answers to the following word questions. (a) Why is dimensional analysis
More informationVorticity Equation Marine Hydrodynamics Lecture 9. Return to viscous incompressible flow. N-S equation: v. Now: v = v + = 0 incompressible
13.01 Marine Hydrodynamics, Fall 004 Lecture 9 Copyright c 004 MIT - Department of Ocean Engineering, All rights reserved. Vorticity Equation 13.01 - Marine Hydrodynamics Lecture 9 Return to viscous incompressible
More informationComplex Analysis MATH 6300 Fall 2013 Homework 4
Complex Analysis MATH 6300 Fall 2013 Homework 4 Due Wednesday, December 11 at 5 PM Note that to get full credit on any problem in this class, you must solve the problems in an efficient and elegant manner,
More informationAERODYNAMIC ANALYSIS OF THE HELICOPTER ROTOR USING THE TIME-DOMAIN PANEL METHOD
7 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES AERODYNAMIC ANALYSIS OF THE HELICOPTER ROTOR USING THE TIME-DOMAIN PANEL METHOD Seawook Lee*, Hyunmin Choi*, Leesang Cho*, Jinsoo Cho** * Department
More informationApplied Fluid Mechanics
Applied Fluid Mechanics 1. The Nature of Fluid and the Study of Fluid Mechanics 2. Viscosity of Fluid 3. Pressure Measurement 4. Forces Due to Static Fluid 5. Buoyancy and Stability 6. Flow of Fluid and
More informationAA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow
AA210A Fundamentals of Compressible Flow Chapter 1 - Introduction to fluid flow 1 1.2 Conservation of mass Mass flux in the x-direction [ ρu ] = M L 3 L T = M L 2 T Momentum per unit volume Mass per unit
More informationAerothermodynamics of high speed flows
Aerothermodynamics of high speed flows AERO 0033 1 Lecture 6: D potential flow, method of characteristics Thierry Magin, Greg Dimitriadis, and Johan Boutet Thierry.Magin@vki.ac.be Aeronautics and Aerospace
More informationThe E80 Wind Tunnel Experiment the experience will blow you away. by Professor Duron Spring 2012
The E80 Wind Tunnel Experiment the experience will blow you away by Professor Duron Spring 2012 Objectives To familiarize the student with the basic operation and instrumentation of the HMC wind tunnel
More informationi.e. the conservation of mass, the conservation of linear momentum, the conservation of energy.
04/04/2017 LECTURE 33 Geometric Interpretation of Stream Function: In the last class, you came to know about the different types of boundary conditions that needs to be applied to solve the governing equations
More informationADVERSE REYNOLDS NUMBER EFFECT ON MAXIMUM LIFT OF TWO DIMENSIONAL AIRFOILS
ICAS 2 CONGRESS ADVERSE REYNOLDS NUMBER EFFECT ON MAXIMUM LIFT OF TWO DIMENSIONAL AIRFOILS Kenji YOSHIDA, Masayoshi NOGUCHI Advanced Technology Aircraft Project Center NATIONAL AEROSPACE LABORATORY 6-
More informationShape Optimization of Low Speed Airfoils using MATLAB and Automatic Differentiation. Christian Wauquiez
Shape Optimization of Low Speed Airfoils using MATLAB and Automatic Differentiation Christian Wauquiez Stockholm 2000 Licentiate s Thesis Royal Institute of Technology Department of Numerical Analysis
More informationModel to Evaluate the Aerodynamic Energy Requirements of Active Materials in Morphing Wings
Model to Evaluate the Aerodynamic Energy Requirements of Active Materials in Morphing Wings Gregory W. Pettit Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University
More information1. Introduction, tensors, kinematics
1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and
More information= π + sin π = π + 0 = π, so the object is moving at a speed of π feet per second after π seconds. (c) How far does it go in π seconds?
Mathematics 115 Professor Alan H. Stein April 18, 005 SOLUTIONS 1. Define what is meant by an antiderivative or indefinite integral of a function f(x). Solution: An antiderivative or indefinite integral
More informationHW6. 1. Book problems 8.5, 8.6, 8.9, 8.23, 8.31
HW6 1. Book problems 8.5, 8.6, 8.9, 8.3, 8.31. Add an equal strength sink and a source separated by a small distance, dx, and take the limit of dx approaching zero to obtain the following equations for
More information3 Generation and diffusion of vorticity
Version date: March 22, 21 1 3 Generation and diffusion of vorticity 3.1 The vorticity equation We start from Navier Stokes: u t + u u = 1 ρ p + ν 2 u 1) where we have not included a term describing a
More information