Aerodynamics. Finite Wings Lifting line theory Glauert s method

Similar documents
Finite Wings Steady, incompressible flow

AERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY

NUMERICAL DIFFERENTIATION

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Appendix B. The Finite Difference Scheme

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

Numerical Heat and Mass Transfer

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

The Finite Element Method

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Perfect Competition and the Nash Bargaining Solution

Chapter 2 - The Simple Linear Regression Model S =0. e i is a random error. S β2 β. This is a minimization problem. Solution is a calculus exercise.

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation

A random variable is a function which associates a real number to each element of the sample space

Chapter 11: Simple Linear Regression and Correlation

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

Principles of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer

Here is the rationale: If X and y have a strong positive relationship to one another, then ( x x) will tend to be positive when ( y y)

Topic 5: Non-Linear Regression

CHAPTER 14 GENERAL PERTURBATION THEORY

ACTM State Calculus Competition Saturday April 30, 2011

e i is a random error

Georgia Tech PHYS 6124 Mathematical Methods of Physics I

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

AERODYNAMICS INVESTIGATION OF FINITE WINGS BY LIFTING-LINE MODELS USING THE TRACK METHOD

Lecture 16 Statistical Analysis in Biomaterials Research (Part II)

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Chapter 13: Multiple Regression

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models

Convexity preserving interpolation by splines of arbitrary degree

Indeterminate pin-jointed frames (trusses)

Introduction to Random Variables

Analysis of Discrete Time Queues (Section 4.6)

PART 8. Partial Differential Equations PDEs

Single Variable Optimization

Chapter 4. Velocity analysis

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models

ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

Econ Statistical Properties of the OLS estimator. Sanjaya DeSilva

Inductance Calculation for Conductors of Arbitrary Shape

One-sided finite-difference approximations suitable for use with Richardson extrapolation

DUE: WEDS FEB 21ST 2018

Interconnect Modeling

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

A Hybrid Variational Iteration Method for Blasius Equation

(c) (cos θ + i sin θ) 5 = cos 5 θ + 5 cos 4 θ (i sin θ) + 10 cos 3 θ(i sin θ) cos 2 θ(i sin θ) 3 + 5cos θ (i sin θ) 4 + (i sin θ) 5 (A1)

Linear Regression Analysis: Terminology and Notation

9.2 Seismic Loads Using ASCE Standard 7-93

Newton s Method for One - Dimensional Optimization - Theory

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.

In this section is given an overview of the common elasticity models.

Correlation and Regression. Correlation 9.1. Correlation. Chapter 9

Implicit Integration Henyey Method

Solving Fractional Nonlinear Fredholm Integro-differential Equations via Hybrid of Rationalized Haar Functions

Statistics for Managers Using Microsoft Excel/SPSS Chapter 13 The Simple Linear Regression Model and Correlation

Numerical Solution of Ordinary Differential Equations

Pressure Measurements Laboratory

Feb 14: Spatial analysis of data fields

Chapter 3. Two-Variable Regression Model: The Problem of Estimation

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

Lecture 21: Numerical methods for pricing American type derivatives

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

A Comparative Investigation into Aerodynamic Performances of Two Set Finned Bodies with Circular and Non Circular Cross Sections

Statistics for Business and Economics

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method

ORIGIN 1. PTC_CE_BSD_3.2_us_mp.mcdx. Mathcad Enabled Content 2011 Knovel Corp.

Polynomial Regression Models

#64. ΔS for Isothermal Mixing of Ideal Gases

APPENDIX A Some Linear Algebra

Mathematical Preparations

Laboratory 1c: Method of Least Squares

Error Probability for M Signals

6. Stochastic processes (2)

6. Stochastic processes (2)

Solutions Homework 4 March 5, 2018

Lecture 3: Probability Distributions

Department of Quantitative Methods & Information Systems. Time Series and Their Components QMIS 320. Chapter 6

Multigrid Methods and Applications in CFD

Lecture Notes on Linear Regression

Introduction to Statistical Methods

Some notes on Futaki invariant

Laboratory 3: Method of Least Squares

2 Finite difference basics

Primer on High-Order Moment Estimators

2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12

Lecture 2: Numerical Methods for Differentiations and Integrations

Grid Generation around a Cylinder by Complex Potential Functions

Aging model for a 40 V Nch MOS, based on an innovative approach F. Alagi, R. Stella, E. Viganò

Affine and Riemannian Connections

An identification algorithm of model kinetic parameters of the interfacial layer growth in fiber composites

ECONOMICS 351*-A Mid-Term Exam -- Fall Term 2000 Page 1 of 13 pages. QUEEN'S UNIVERSITY AT KINGSTON Department of Economics

P exp(tx) = 1 + t 2k M 2k. k N

Digital Signal Processing

Transcription:

α ( y) l Γ( y) r ( y) V c( y) β b 4 V Glauert s method b ( y) + r dy dγ y y dy Soluton procedure that transforms the lftng lne ntegro-dfferental equaton nto a system of algebrac equatons - Restrcted to symmetrc wngs wth no dhedral or sweep angles

Glauert s method Substtuton of ndependent varable y cos( θ ) b bsen dy ( θ ) θ d α ( θ ) l Γ( θ ) r ( θ ) V c( θ ) β ( θ ) + r 4 V 0 b ( cos( θ ) cos( θ )) dγ dθ dθ

The crculaton must be zero at the wng tps θ 0 and θ and symmetrc Γ ( θ ) Γ( θ ) The Fourer seres expanson of Γ( θ ) that satsfes these condtons s gven by Γ Glauert s method ( θ ) Γ sen( nθ ) n n,3,... The cosne terms are elmnated by the tp condtons and symmetry requres the elmnaton of all terms wth even n

Glauert s method The frst term of the seres corresponds to an ellptcal crculaton dstrbuton Γ ( θ ) Γ sen( θ ) Γ cos ( θ ) or Γ y b Γ Γ y + b

ω 4 Glauert s method Determnaton of the descendng velocty (downwash) 0 b( cos( θ ) cos( θ )) n,3,... nγ n cos ( nθ ) d θ ω b sen For a term seres ( θ ) n,3,... nγ n sen ( nθ ) ω Γ b

ftng lne equaton Glauert s method α α eff + α α ( θ ) l Γ n n,3,... sen r ( nθ ) ( θ ) V c( θ ) β ( θ ) + n,3,... nγ r bv n sen sen ( nθ ) ( θ ) Γ sen ( nθ ) r ( θ ) V c( θ ) n n,3,... l bv + ( nθ ) ( θ ) nsen r sen α ( θ ) + β ( θ )

hoosng a fnte number (N) of terms for the seres expanson leads to a system of algebrac equatons that satsfes + ( N ) ( nθ ) n ( nθ ) sen sen Γ n α θ ( ) ( ) ( ) + β θ r + r n l θ V c θ bv θ,3,... sen for N ponts of the nterval 0 < θ wth θ, wth,,..., N N Glauert s method ( ) ( )

Example for N3 Γ ( θ ) Γ sen( θ ) + Γ sen( 3θ ) + Γ sen( 5θ ) A A A θ 3 A A A 3 Glauert s method, 6 A A A 3 3 33 θ Γ Γ Γ 3, 3 3 5 θ α α α 3 5 ( θ) + β ( θ) ( θ ) + β ( θ ) ( ) ( ) θ3 + β θ3

A j Example for N3 sen Γ l (( + ( j ) ) θ ) ( + ( j ) ) sen( ( + ( j ) ) θ ) r + r ( θ ) V c( θ ) bv sen( θ ) sen( 5 6) 5sen( 5 6) A3 r + r ( 6) V c( 6) bv sen( 6) ( θ ) Γ sen( θ ) + Γ sen( 3θ ) + Γ sen( 5θ ) θ, 6 l Glauert s method θ 3, 3 θ 3 5

r ρ V 0 sen 0 Glauert s method Wth the knowledge of the seres coeffcents, Γ n, t s possble to calculate the lft force,, and the nduced drag force, D Γ n n,3,... ( mθ ) sen( nθ ) sen dθ b ( nθ ) sen( θ ) 0 m m dθ n n

Glauert s method Wth the knowledge of the seres coeffcents, Γ n, t s possble to calculate the lft force,, and the nduced drag force, D r ρbv Γ 4 The lft force depends only on the frst term of the seres. Obvously, ths does not mean that term s enough to calculate the exact value of

Glauert s method Wth the knowledge of the seres coeffcents, Γ n, t s possble to calculate the lft force,, and the nduced drag force, D D 4 ( nθ ) Γ sen( nθ ) ρ nγ n sen 0 n,3,... n,3,... n dθ D ρ nγn 8 n,3,...

Glauert s method Wth the knowledge of the seres coeffcents, Γ n, t s possble to calculate the lft force,, and the nduced drag force, D D δ ργ 8 ( + δ ) Γ n n n 3,5,... Γ

Glauert s method Wth the knowledge of the seres coeffcents, Γ n, t s possble to calculate the lft force,, and the nduced drag force, D Γ 4 r ρbv D r ρ V b ( + δ )

ft and nduced drag coeffcents ( ) δ ρ ρ + Γ D S V D bv S V r r r Glauert s method

Glauert s method D D r ρ V S ( + δ ) The nduced drag coeffcent s proportonal to the square of the lft coeffcent and t goes to zero the aspect rato,, tends to nfnty,.e. b

Wngs wth mnmum nduced drag The nduced drag coeffcent s gven by D D ρ V r ( + δ ) wth S δ Γ n n n 3,5,... Γ By defnton δ 0, therefore, the mnmum nduced drag corresponds to δ 0,.e. Γ 0 for n > n

Wngs wth mnmum nduced drag One term seres corresponds to an ellptc crculaton dstrbuton y Γ( θ ) Γ sen( θ ) Γ cos ( θ ) Γ b The nduced velocty, ω, s constant along the span and equal to ω Γ b

ft and nduced drag coeffcents gven by Γ Γ ρ ρ,, 4 D b V D bv bv r r r Wngs wth mnmum nduced drag

Wngs wth mnmum nduced drag onstant nduced angle of attack along the span α Γ bv r If the wng has no twst the geometrc angle of attack does not change along the span, ths leads to a constant α α eff α

Wngs wth mnmum nduced drag The lft coeffcent,, may be calculated from b b b l dy For a wng wth ellptcal crculaton dstrbuton ( e α ct ) and no twst ( e α ct ) and the same arfol secton along the span b l dy l b b

Wngs wth mnmum nduced drag For l, becomes or α α α eff + α l l β + + ( α + β )

Wngs wth mnmum nduced drag The rato between the slopes of the lnes that gve the lft coeffcent as a functon of the angle of attack for the wng and ts secton (arfol) s gven by l l + If we assume l then l +

l 0.9 0.8 0.7 Wngs wth mnmum nduced drag l + 0.6 0.5 0.4 0.3 0. 0. 0 0 4 6 8 0

α Wngs wth mnmum nduced drag An ellptcal crculaton dstrbuton may be obtaned for a wng wth no twst and the same secton (arfol) along the span sen( θ ) r β + V c( θ ) ( θ ) c sen( θ ) Γ Γ r l bv c r where the root chord, c r, s related to the requred lft and to the geometrc propertes of the wng secton (arfol thckness and camber)

Wngs wth mnmum nduced drag

Wngs wth mnmum nduced drag

Wngs wth mnmum nduced drag

Wngs wth mnmum nduced drag A wng wth ellptcal planform s harder to buld than a rectangular wng A tapered wng may produce a nearly ellptcal crculaton dstrbuton wth a smpler constructon than an ellptcal wng Taper rato s the rato between the tp chord, c t and the root chord, c r

Wngs wth mnmum nduced drag AR e δ Γ n n n 3,5,... Γ

Wngs wth mnmum nduced drag

Introducton of vscous effects The fnte wng lft coeffcent depends on the aspect rato and on the lfy coeffcent of the wng sectons (arfols) l ( y) ( y) α ( y) + β ( y) l b b b ( ) eff ( y) If l and β are determned for a vscous flow the effects of vscosty are ncluded n the determnaton of the wng lft coeffcent l dy

Introducton of vscous effects For a vscous flud, the drag coeffcent of the wng sectons (arfols) s not zero. Therefore, D d + D d + + + δ arfol frcton d pressure arfol wth ( ) b d d ( y) dy frcton arfol b b τw b ( ) b d d ( y) pressure ( ) ( ) arfol b δ * dy

In a real flow, there s only frcton and pressure drag (resstance) D where pressure + D D frcton D pressure ( ) b ( y) b ( ) b + ( + δ ) ( y) dy + ( + δ ) d Introducton of vscous effects pressure D frcton arfol d frcton arfol b b b d δ * d τw dy

Prandtl transformaton formulae Prandtl s transformaton formulae allow the converson of the functons D ( ) and (α) obtaned (expermentally or numercally) for a gven aspect rato, nto functons vald for wngs wth the same secton but dfferent aspect ratos - It s assumed that the tp effects depend only on the aspect rato,, and that the nduced drag and nduced angle of attack correspond to the values obtaned for an ellptcal crculaton dstrbuton

Assumptons: - Induced resstance coeffcent gven by - Induced angle of attack equal to - Drag coeffcent s equal to the sum of the wng secton (arfol) drag and nduced drag coeffcents -, β and are ndependent of l darfol Prandtl transformaton formulae - The lft coeffcent of the wng sectons s constant along the span and equal to

+ + d D d D perfl perfl + + + + D D D d D d D arfol arfol Prandtl transformaton formulae D as a functon of : For two wngs wth dfferent aspect ratos, and, and the same

Prandtl transformaton formulae D as a functon of :

α as a functon of : For two wngs wth dfferent aspect ratos, and, and the same + + β α α α α l eff + + + + α α α β α β α l l Prandtl transformaton formulae

Prandtl transformaton formulae α as a functon of :