An alternative approach to integral equation method based on Treftz solution for inviscid incompressible flow
|
|
- Lucas Lynch
- 5 years ago
- Views:
Transcription
1 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, Covilhd 6200, Portugal Jan A. Kolodziej Poznan University of Technology, Institute ofapplied Mechanics, Poznan, Poland Abstract The paper deals with the numerical modelling of inviscid, incompressible flow about two-dimensional lifting surfaces at low speed. The problem is formulated in terms of the Laplace equation. For its solution two different approaches are discussed. The first one tackles the problem in terms of the velocity potential with normal derivative prescribed at the boundary. An integral equation method based on a distribution of vortices is utilised to solve the flow. As for the second approach the governing equation is written in terms of the stream function. Treftz method is then used to express the solution as a power series of the fundamental solutions of Laplace equation. 1 Introduction Modelling of inviscid incompressible flow about lifting surfaces of an arbitrary geometry may be conducted by means of integral or differential
2 58 Advances in Fluid Mechanics II representations. The problem is governed by a field equation which has originally been solved by numerical methods applied to a grid extended along the fluid domain. By making use of Green's theorem the solution of the problem may be reduced to solve an integral equation written for the physical boundaries of the flow. In this context the classical approach by surface singularity distributions has originated the well-known panel methods in two and three-dimensions - Hess [1]. Another classical approach, which is not so commonly applied to solve this class of problems, is the Treftz Method - Kolodziej [2], Herein we introduce this methodology as an alternative to the more traditional integral equation approach. In this last approach the problem is formulated in terms of the velocity potential with normal derivative prescribed at the boundary. This harmonic function is afterwards searched for as the velocity field generated by a distribution of vortices along the body contour. The resulting Fredholm integral equation, of the first kind, is subsequently solved for the tangential velocity at a discrete number of panels that are distributed along the body profile. The numerical results that were obtained by modelling a NACA 0012 airfoil have been subsequently compared to those experimentally obtained by the authors in a low-speed wind tunnel. 2 Potential flow equations Let us consider a two-dimensional lifting surface of an arbitrary shaped cross-section, like the one represented in Fig. 1. Its contour ( C ) is a closed line having continuous derivative everywhere except at the trailing edge, whose outward normal n points to the unbounded fluid domain (D). Here OXYZ is a stationary co-ordinate system and oxyz is a system of axes attached to the body. Fig. 1 Wing section moving with speed U and angle of attack a in fluid domain (D).
3 Advances in Fluid Mechanics II 59 The solid body advances with velocity U in air which is initially at rest; a is the angle between the chord c and the direction of motion of the wing. At low speed the flow induced in the vicinity of the body may be assumed as incompressible. Moreover, if we assume inviscid, irrotational flow, the fluid velocity V may be derived from an absolute potential. This function is the solution of the following system of equations: V^0(z) = 0, (zed) (1) - " ' m where z stands for the complex co-ordinate z=x+iy. The velocity potential is, therefore, the solution of Laplace equation (1) and the impermeability condition (2), a condition to be applied on the solid boundaries of the flow. At infinity a regularity condition applies which may be written in terms of the relative potential: = -U, (z-»±oo) (3) 3 Integral-equation approach The harmonic function formulated above is herein searched for as the real part of a complex velocity potential F(z). This velocity field will be generated in every point z of the fluid domain (D) by a distribution of Rankine vortices with density y(z') at a generic point along the contour (z' At this stage the velocity field associated with this function satisfies automatically to the governing equation (1) and must equally verify the boundary conditions stated above. As to condition (2) imposed on the surface of the body, this condition may be formulated in the integral form as follows:
4 60 Advances in Fluid Mechanics II = 0.5 (5) taking into account eq (3). In the previous equation the kernel of the integral may be easily determined after differentiation of the function defined in eq (4). The angle 0 represents here the argument of the unit-vector tangent to the contour of the body. Equation (5) is a Fredholm integral equation, of the first kind, which is to be solved in order to the vortex densities y, for a prescribed body motion with velocity U. Fig. 2 Discretized contour of the wing section. The solution of this equation for a body having a contour of arbitrary geometry requires the use of a numerical approach. A numerical treatment that is commonly adopted consists of sub-dividing the contour ( C ) in N linear elements Sj (j=l..,n), as it is shown in Fig. 2. Each of these elements Sj is the support of a uniform distribution of vortices with density yj. Sj is defined by its extremity points Zj and Zj+i, and its argument 8,; Z=z/c are non-dimensional complex co-ordinates. By enforcing the impermeability condition on N control points TJ taken as the middle point of each element of the contour, eq (5) assumes the discretized form:, = M,,U = i..,n) (6) The influence coefficients of the velocity By may be calculated after integration of the left-hand side of eq (5), taking yj constant along each linear boundary element.
5 Advances in Fluid Mechanics II 61 B,, =0 (i- j); As to the right-hand side of the linear system, n, is the scalar product of the unit-vector in the direction of the body motion to the normal unitvector at the middle point ij=0.5(zi+zi+i). The system of N equations (6) is to be solved for the N unknown normalised densities (yj/u). In order to get a finite value for the velocity at the wing trailing edge and the system to be determined it is, however, necessary to impose an additional condition at this point. Although Kutta condition may be imposed in few different ways, a convenient form for non-cusped airfoils is to ensure equal velocities on either side of the trailing edge, as suggested by Anderson [3]: yi+y#=0. This is achieved by simply replacing the last equation of system (6) by the former relationship. Finally the pressure coefficient is derived from Bernoulli's law: Following Katz & Plotkin [4], the tangential velocity within the fluid adjacent to the profile may be obtained directly from the vortices density that were calculated above, on each panel of the body. The lift coefficient of the airfoil may then be derived from Joukowski theorem as a function of the circulation F: C - F ^L- (9) 4 Treftz solution The governing equation is here written in terms of the stream function Y of the flow and the value of this function is prescribed at the boundary. In two-dimensional potential flow, Cauchy-Riemann relationships applied toeq(l) yell the following Laplace equation:
6 62 Advances in Fluid Mechanics II (10) The solution for this flow field is now assumed to be the superposition of a uniform stream of constant velocity U and incidence a, with a perturbation of the flow which is due to the presence of the body. k=\ i=\ (11) where A^ are unknown weighting coefficients of the truncated series that may be determined by means of a simple collocation procedure - Collatz [5] - or by using a Galerkin approximation; y\(x,y,xk,yd are trial functions that exactly fulfil eq (10). These fundamental solutions are sources forz = l, vortices for i=2, doublets for i=3 and multipoles for (12) where n -, X X i The co-ordinates (x^.y^) represent here the positions of these singularities placed inside the contour of the airfoil. Solution (11) fulfils exactly the governing equation (10) and must also satisfy the boundary condition required on the surface of the airfoil, where the stream function must be zero: (* \ ^-0). Keeping this in mind, normal and tangential velocities may be obtained after differentiation of equation (11),
7 Advances in Fluid Mechanics II 63 03) k = \ /=! (14) Here cosg and sing are direction cosines of the unit-vector tangent to the contour of the body, and: g\ 82=- X XL, 3 = X XL, i -4 -t, 5^, (15) /2 ~ ~ /3 = /, = -(z - 2)!r^ '^ cos[(/ - i = 4, 5,... In the numerical procedure used to determine the NxM coefficients A^ in eq (11), the boundary conditions which are to be imposed on the stream function and normal velocity will be fulfilled only approximately, at a finite number of collocation points along the body profile. In this way we are led to an over-determined system of equations
8 64 Advances in Fluid Mechanics II which must be solved by the least-squares method. Collocation of the stream function and Neumann conditions for (P-l) control points Tp(Xp,yp) defined along the boundary of the airfoil, yell the following system of equations: TV k = \ 7 = 1 AT M (16) Moreover, Kutta condition states that the backward stagnation point coincides with the trailing edge: TV k=\ 7=1 k=\ i=\ (18) k = \ 7 = 1 Two different collocation procedures, both in least squares sense, have been used. In the first one only the stream function condition is imposed along with Kutta condition; (P+2) linear equations are to be solved in this case. As for the second one, the collocation involves stream function and its tangential derivative (normal velocity), and Kutta condition is enforced as well. Equations (16), (17) and (18) constitute a linear system of (2P+1) equations to be solved for NxM A^j unknowns. Eq (14) is then used to calculate Cp from eq (8) and eq (9) delivers the lift coefficient.
9 Advances in Fluid Mechanics II 65 5 Numerical results and discussion The numerical results presented here have been obtained for a standard lifting surface (NACA 0012). Fig. 3 gives the pressure coefficient along the wing section, calculated by the integral equation method. The values of Cp are compared with the ones obtained by Moran [6] for the same profile. The lift coefficient is shown in Fig. 4 for different angles of attack a. The computational results are here compared to the results of Abbott and Doenhoff [7] for a Reynolds number equal to 6x10, and to our own experiments at Re =2.5x10. We have used for that purpose a low-speed wind tunnel Flint TE44 (0.45x0.45x0.572m), in which the velocity of air can be continuously varied till 30 m/sec. The model is mounted on an electronic three-component balance TE81/E, equipped with piezo-resistive load cells. A preliminary analysis of the uncertainty of the measured forces points out to a precision of 1% during the calibration procedure. The data acquisition is controlled by an A/D Lab PC card from National Instruments working in Lab VIEW environment. The signals are afterwards processed in order to evaluate the lift and drag coefficients, as well as the aerodynamic moment. The speed of the approaching stream is controlled at the main section of the tunnel by a Prandtl probe Cp x/c Fig. 3 Pressure coefficient along the upper and lower surfaces of the NACA 0012 airfoil. For an angle of attack till about 5 we have obtained a good correlation between the numerical predictions and the experimental results. For oc=10 the lift coefficient as predicted by the integral equation model, differs from the experimental value by about 25%. However, this fact is due to the separation of the flow and the
10 66 Advances in Fluid Mechanics II development of a wake behind the wing. Results comparable to the ones provided above are attained by using the concept of boundary collocation with adaptation - Golik & Kolodziej [8]. Fig. 4 Lift coefficient for the NACA 0012 wing section as a function of the angle of attack. The previous results were calculated for an average of N=200 panels distributed over the body surface. Fig. 5 shows the evolution of the values obtained for CL as we refine the discretization of the boundary. For the NACA 0012, a=8.3 is a particular angle of attack for which CL=!. By applying Joukowski theorem we have obtained a lift coefficient 2% close to this value (0^=0.98), using an average of N=100 panels. From direct integration of the pressure around the wing a more rough estimate of this coefficient is obtained (C[j=0.93), and this is certainly due to the numerical oscillations on the tangential velocity near the trailing edge. The effect of these unwanted oscillations on pressure at the vicinity of this point is visible in Fig. 3. Under these theoretical conditions the wing experiences a small drag Co=0.15, that can be seen as a residual error in fulfilling Kutta condition at the trailing edge of the wing section. This numerical difficulty may be overcome by reasoning on the field induced inside the body - Katz & Plotkin [4] - or by considering higher-order singularity distributions. In general we may conclude that the Integral Equation method reproduces correctly the flow properties. If vortices are used, a formulation based on the normal velocity leads inevitably to an integral equation of the first kind. Consequently, the matrix of the influence coefficients does not exhibit a dominant diagonal; the number of panels that we need to input in order to get an accurate solution is in fact high, if a direct elimination solver with full pivoting is used - Press et al [9]. The
11 Advances in Fluid Mechanics II 67 improvement gained by using a singular value decomposition had no significance in this case. Condition number based on the Euclidean norm is situated in our case between 10* and 10^, for N=10 to 200 boundary elements. Therefore, double precision is required in order to obtain increased accuracy of the numerical results. As it has been pointed out by Nathman [10] the evaluation of the relevant parameters using higher precision may have some impact on the final results for the case of very thin wings. However, the logarithmic nature of the velocity coefficients in eq (7) will certainly be an important limitation, with respect to the simulation of wings with small thickness N Fig. 5 Convergence of the lift coefficient with number of panels. As to the CPU time surveyed on calculations performed by integral equation method on a PC double Pentium processor at 200MHz, with 64 Mbyte RAM and symmetrical multi-processing, it is of the order of msec for N up to 900 and between 1 and 2 sec for N up to 3000 panels. After this value CPU time starts from 22 sec, as a result of limitations in memory space. Acknowledgements The research work that is presented in this paper was carried out within the framework of PRAXIS grant 3/3.1/CTAE/1940/95. The work reflects a close collaboration between Universidade da Beira Interior (Portugal) and The Technical University of Poznan, in Poland. The numerical and experimental results have been obtained at the Laboratory of Fluid Mechanics of Universidade da Beira Interior.
12 68 Advances in Fluid Mechanics II References [1] Hess, J. L., Panel methods in computational fluid dynamics, Ann. Rev. Fluid Mechanics, 22, pp , [2] Kolodziej, J. A., Review of application of boundary collocation method in mechanics of continuous media, Solid Mechanics Archives, 12, pp , [3] Anderson, J. D. Jr., Fundamentals of Aerodynamics, McGraw-Hill International, New York, pp , [4] Katz, J. & Plotkin, A., Low-Speed Aerodynamics: From Wing Theory to Panel Methods, McGraw Hill, New York, [5] Collatz, L., The Numerical Treatment of Differential Equations, Springer-Verlag, Berlin, [6] Moran, J., An Introduction to Theoretical and Computational Aerodynamics, John Wiley & Sons, New York, pp 241, [7] Abbott, I. & Doenhoff, A., Theory of Wing Sections, Dover Publications, pp , New York, [8] Golik W. L. & Kolodziej, J. A., An Adaptive Boundary Collocation Method for Linear PDEs, Numerical Methods for Partial Differential Equations, 11, pp , [9] Press, W. H., Flannery, B. P., Teukolsky, S. A. & Vetterling, W. T., Numerical Recipes, Cambridge University Press, Cambridge, pp 26-29, [10] Nathman, J., Precision Requirements for Potential Based Panel Methods, AIAA Journal, Vol. 32, 5, pp , 1993.
A comparison of velocity and potential based boundary element methods for the analysis of steady 2D flow around foils
A comparison of velocity and potential based boundary element methods for the analysis of steady 2D flow around foils G.B. Vaz, L. E a, J.A.C. Falcao de Campos Department of Mechanical Engineering, Institute
More informationNUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD
NUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD V. G. Guedes a, G. C. R. Bodstein b, and M. H. Hirata c a Centro de Pesquisas de Energia Elétrica Departamento de Tecnologias
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 informationInvestigation potential flow about swept back wing using panel method
INTERNATIONAL JOURNAL OF ENERGY AND ENVIRONMENT Volume 7, Issue 4, 2016 pp.317-326 Journal homepage: www.ijee.ieefoundation.org Investigation potential flow about swept back wing using panel method Wakkas
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 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 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 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 informationACD2503 Aircraft Aerodynamics
ACD2503 Aircraft Aerodynamics Session delivered by: Prof. M. D. Deshpande 1 Aims and Summary PEMP It is intended dto prepare students for participation i i in the design process of an aircraft and its
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 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 informationNumerical Simulation of Unsteady Aerodynamic Coefficients for Wing Moving Near Ground
ISSN -6 International Journal of Advances in Computer Science and Technology (IJACST), Vol., No., Pages : -7 Special Issue of ICCEeT - Held during -5 November, Dubai Numerical Simulation of Unsteady Aerodynamic
More informationLifting 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 informationAerodynamic Rotor Model for Unsteady Flow and Wake Impact
Aerodynamic Rotor Model for Unsteady Flow and Wake Impact N. Bampalas, J. M. R. Graham Department of Aeronautics, Imperial College London, Prince Consort Road, London, SW7 2AZ June 28 1 (Steady Kutta condition)
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 informationCorrection of Wind Tunnel Results for the Airfoils of ITA s Unmanned Aerial Vehicle
Correction of Wind Tunnel Results for the Airfoils of ITA s Unmanned Aerial Vehicle Carlos Diego Aguiar de Nogueira Gomes Empresa Brasileira de Aeronáutica S. A. (Embraer). Av. Brigadeiro Faria Lima, 2170,
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 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 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 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 informationIncompressible Flow Over Airfoils
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
More informationAerodynamic force analysis in high Reynolds number flows by Lamb vector integration
Aerodynamic force analysis in high Reynolds number flows by Lamb vector integration Claudio Marongiu, Renato Tognaccini 2 CIRA, Italian Center for Aerospace Research, Capua (CE), Italy E-mail: c.marongiu@cira.it
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 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 informationComputational Fluid Dynamics Study Of Fluid Flow And Aerodynamic Forces On An Airfoil S.Kandwal 1, Dr. S. Singh 2
Computational Fluid Dynamics Study Of Fluid Flow And Aerodynamic Forces On An Airfoil S.Kandwal 1, Dr. S. Singh 2 1 M. Tech Scholar, 2 Associate Professor Department of Mechanical Engineering, Bipin Tripathi
More informationExperimental Evaluation of Aerodynamics Characteristics of a Baseline Airfoil
Research Paper American Journal of Engineering Research (AJER) e-issn: 2320-0847 p-issn : 2320-0936 Volume-4, Issue-1, pp-91-96 www.ajer.org Open Access Experimental Evaluation of Aerodynamics Characteristics
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 informationAn Experimental Validation of Numerical Post-Stall Aerodynamic Characteristics of a Wing
Proceedings of ICTACEM 2017 International Conference on Theoretical, Applied, Computational and Experimental Mechanics December 28-30, 2017, IIT Kharagpur, India ICTACEM-2017/XXXX(paper No.) An Experimental
More informationUNSTEADY EFFECTS IN 2D MODELLING OF PARTIAL CAVITATION USING BEM
MÉTODOS COMPUTACIONAIS EM ENGENHARIA Lisboa, 31 de Maio - 2 de Junho, 24 APMTAC, Portugal 24 UNSTEADY EFFECTS IN 2D MODELLING OF PARTIAL CAVITATION USING BEM Guilherme Vaz, Johan Bosschers and J. A. C.
More informationUNSTEADY AERODYNAMIC ANALYSIS OF HELICOPTER ROTOR BY USING THE TIME-DOMAIN PANEL METHOD
6 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES UNSTEAD AERODNAMIC ANALSIS OF HELICOPTER ROTOR B USING THE TIME-DOMAIN PANEL METHOD Seawook Lee*, Leesang Cho*, Jinsoo Cho* *Hanyang University
More informationVORTEX METHOD APPLICATION FOR AERODYNAMIC LOADS ON ROTOR BLADES
EWEA 2013: Europe s Premier Wind Energy Event, Vienna, 4-7 February 2013 Figures 9, 10, 11, 12 and Table 1 corrected VORTEX METHOD APPLICATION FOR AERODYNAMIC LOADS ON ROTOR BLADES Hamidreza Abedi *, Lars
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 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 informationANALYSIS OF HORIZONTAL AXIS WIND TURBINES WITH LIFTING LINE THEORY
ANALYSIS OF HORIZONTAL AXIS WIND TURBINES WITH LIFTING LINE THEORY Daniela Brito Melo daniela.brito.melo@tecnico.ulisboa.pt Instituto Superior Técnico, Universidade de Lisboa, Portugal December, 2016 ABSTRACT
More informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
More informationAERODYNAMICS STUDY NOTES UNIT I REVIEW OF BASIC FLUID MECHANICS. Continuity, Momentum and Energy Equations. Applications of Bernouli s theorem
AERODYNAMICS STUDY NOTES UNIT I REVIEW OF BASIC FLUID MECHANICS. Continuity, Momentum and Energy Equations. Applications of Bernouli s theorem UNIT II TWO DIMENSIONAL FLOWS Complex Potential, Point Source
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 informationSimulation of Aeroelastic System with Aerodynamic Nonlinearity
Simulation of Aeroelastic System with Aerodynamic Nonlinearity Muhamad Khairil Hafizi Mohd Zorkipli School of Aerospace Engineering, Universiti Sains Malaysia, Penang, MALAYSIA Norizham Abdul Razak School
More informationAerodynamics. Professor: Luís Eça
Aerodynamics Professor: Luís Eça Program 1. Introduction Aerodynamical forces. Flow description. Dependent variables and physical principles that govern the flow Program 2. Incompressible, Viscous Flow
More informationEnclosure enhancement of flight performance
THEORETICAL & APPLIED MECHANICS LETTERS, 23 (21) Enclosure enhancement of flight performance Mehdi Ghommem, 1, a) Daniel Garcia, 2 Victor M. Calo 3 1) Center for Numerical Porous Media (NumPor), King Abdullah
More informationNumerical study of battle damaged two-dimensional wings
Advances in Fluid Mechanics IX 141 Numerical study of battle damaged two-dimensional wings S. Djellal, T. Azzam, M. Djellab & K. Lakkaichi Fluid Mechanics Laboratory Polytechnical School Bordj El Bahri,
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 informationMasters in Mechanical Engineering. Problems of incompressible viscous flow. 2µ dx y(y h)+ U h y 0 < y < h,
Masters in Mechanical Engineering Problems of incompressible viscous flow 1. Consider the laminar Couette flow between two infinite flat plates (lower plate (y = 0) with no velocity and top plate (y =
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 informationFundamentals of Aerodynamics
Fundamentals of Aerodynamics Fourth Edition John D. Anderson, Jr. Curator of Aerodynamics National Air and Space Museum Smithsonian Institution and Professor Emeritus University of Maryland Me Graw Hill
More informationActuator Surface Model for Wind Turbine Flow Computations
Actuator Surface Model for Wind Turbine Flow Computations Wen Zhong Shen* 1, Jens Nørkær Sørensen 1 and Jian Hui Zhang 1 Department of Mechanical Engineering, Technical University of Denmark, Building
More informationSelf-adaptive Trefftz procedure for harmonic problems
Computer Assisted Mechanics and Engineering Sciences, 4: 491 500, 1997. Copyright c 1997 by Polska Akademia Nauk Self-adaptive Trefftz procedure for harmonic problems Jan A. Ko lodziej and Roman Starosta
More informationAerodynamics. Professor: Luís Eça
Professor: Luís Eça 1. Introduction Aerodynamical forces. Flow description. Dependent variables and physical principles that govern the flow 2. Incompressible, Viscous Flow Analytical solutions of the
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 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 informationFundamentals of Aerodynamits
Fundamentals of Aerodynamits Fifth Edition in SI Units John D. Anderson, Jr. Curator of Aerodynamics National Air and Space Museum Smithsonian Institution and Professor Emeritus University of Maryland
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 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 informationAerodynamics. Lecture 1: Introduction - Equations of Motion G. Dimitriadis
Aerodynamics Lecture 1: Introduction - Equations of Motion G. Dimitriadis Definition Aerodynamics is the science that analyses the flow of air around solid bodies The basis of aerodynamics is fluid dynamics
More informationNUMERICAL INVESTIGATION OF VERTICAL AXIS WIND TURBINE WITH TWIST ANGLE IN BLADES
Eleventh International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 7-9 December 05 NUMERICAL INVESTIGATION OF VERTICAL AXIS WIND TURBINE WITH TWIST ANGLE IN BLADES
More informationINFLUENCE OF ACOUSTIC EXCITATION ON AIRFOIL PERFORMANCE AT LOW REYNOLDS NUMBERS
ICAS 2002 CONGRESS INFLUENCE OF ACOUSTIC EXCITATION ON AIRFOIL PERFORMANCE AT LOW REYNOLDS NUMBERS S. Yarusevych*, J.G. Kawall** and P. Sullivan* *Department of Mechanical and Industrial Engineering, University
More informationDepartment of Energy Sciences, LTH
Department of Energy Sciences, LTH MMV11 Fluid Mechanics LABORATION 1 Flow Around Bodies OBJECTIVES (1) To understand how body shape and surface finish influence the flow-related forces () To understand
More informationAerodynamics of Spinning Sphere in Ideal Flow R. C. Mehta Department of Aeronautical Engineering, Noorul Islam University, Kumaracoil , India
Scholars Journal of Engineering and Technology (SJET) Sch. J. Eng. Tech., 06; 4(5):5-9 Scholars Academic and Scientific Publisher (An International Publisher for Academic and Scientific Resources) www.saspublisher.com
More informationPart 3. Stability and Transition
Part 3 Stability and Transition 281 Overview T. Cebeci 1 Recent interest in the reduction of drag of underwater vehicles and aircraft components has rekindled research in the area of stability and transition.
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 informationActive Control of Separated Cascade Flow
Chapter 5 Active Control of Separated Cascade Flow In this chapter, the possibility of active control using a synthetic jet applied to an unconventional axial stator-rotor arrangement is investigated.
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
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 informationPEMP ACD2505. M.S. Ramaiah School of Advanced Studies, Bengaluru
Governing Equations of Fluid Flow Session delivered by: M. Sivapragasam 1 Session Objectives -- At the end of this session the delegate would have understood The principle of conservation laws Different
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 informationAirfoil Lift Measurement by Surface Pressure Distribution Lab 2 MAE 424
Airfoil Lift Measurement by Surface Pressure Distribution Lab 2 MAE 424 Evan Coleman April 29, 2013 Spring 2013 Dr. MacLean 1 Abstract The purpose of this experiment was to determine the lift coefficient,
More informationPropulsion Systems and Aerodynamics MODULE CODE LEVEL 6 CREDITS 20 Engineering and Mathematics Industrial Collaborative Engineering
TITLE Propulsion Systems and Aerodynamics MODULE CODE 55-6894 LEVEL 6 CREDITS 20 DEPARTMENT Engineering and Mathematics SUBJECT GROUP Industrial Collaborative Engineering MODULE LEADER Dr. Xinjun Cui DATE
More informationIntroduction to Aerospace Engineering
4. Basic Fluid (Aero) Dynamics Introduction to Aerospace Engineering Here, we will try and look at a few basic ideas from the complicated field of fluid dynamics. The general area includes studies of incompressible,
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 informationDesign and simulation of flow through subsonic axial compressor
Design and simulation of flow through subsonic axial compressor Mutahir Ahmed 1, Saeed Badshah 2, Rafi Ullah Khan 2, Muhammad sajjad 2, Sakhi Jan 2 1 NESCOM, Islamabad Pakistan 2 Department of Mechanical
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationUNIT IV BOUNDARY LAYER AND FLOW THROUGH PIPES Definition of boundary layer Thickness and classification Displacement and momentum thickness Development of laminar and turbulent flows in circular pipes
More informationUnsteady Subsonic Aerodynamic Characteristics of Wing in Fold Motion
Technical Paper DOI:10.5139/IJASS.2011.12.1.63 Unsteady Subsonic Aerodynamic Characteristics of Wing in Fold Motion Yoo-Yeon Jung* School of Mechanical and Aerospace Engineering, Seoul National University,
More informationDEVELOPMENT OF A THREE-DIMENSIONAL TIGHTLY COUPLED EULER/POTENTIAL FLOW SOLVER FOR TRANSONIC FLOW
DEVELOPMENT OF A THREE-DIMENSIONAL TIGHTLY COUPLED EULER/POTENTIAL FLOW SOLVER FOR TRANSONIC FLOW Yeongmin Jo*, Se Hwan Park*, Duck-Joo Lee*, and Seongim Choi *Korea Advanced Institute of Science and Technology,
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 informationMarine Hydrodynamics Prof.TrilochanSahoo Department of Ocean Engineering and Naval Architecture Indian Institute of Technology, Kharagpur
Marine Hydrodynamics Prof.TrilochanSahoo Department of Ocean Engineering and Naval Architecture Indian Institute of Technology, Kharagpur Lecture - 10 Source, Sink and Doublet Today is the tenth lecture
More informationDepartment of Mechanical Engineering
Department of Mechanical Engineering AMEE401 / AUTO400 Aerodynamics Instructor: Marios M. Fyrillas Email: eng.fm@fit.ac.cy HOMEWORK ASSIGNMENT #2 QUESTION 1 Clearly there are two mechanisms responsible
More informationWhy airplanes fly, and ships sail
Why airplanes fly, and ships sail A. Eremenko April 23, 2013 And windmills rotate and propellers pull, etc.... Denote z = x + iy and let v(z) = v 1 (z) + iv 2 (z) be the velocity field of a stationary
More informationComplex functions in the theory of 2D flow
Complex functions in the theory of D flow Martin Scholtz Institute of Theoretical Physics Charles University in Prague scholtz@utf.mff.cuni.cz Faculty of Transportation Sciences Czech Technical University
More informationLEE-SIDE FLOW SIMULATIONS OF CRUCIFORM WING- BODY CONFIGURATIONS AT INCOMPRESSIBLE MACH NUMBERS
LEE-SIDE FLOW SIMULATIONS OF CRUCIFORM WING- BODY CONFIGURATIONS AT INCOMPRESSIBLE MACH NUMBERS Janine Versteegh* ** *University of the Witwatersrand **Council for Scientific and Industrial Research (CSIR)
More informationEdwin van der Weide and Magnus Svärd. I. Background information for the SBP-SAT scheme
Edwin van der Weide and Magnus Svärd I. Background information for the SBP-SAT scheme As is well-known, stability of a numerical scheme is a key property for a robust and accurate numerical solution. Proving
More informationBLUFF-BODY AERODYNAMICS
International Advanced School on WIND-EXCITED AND AEROELASTIC VIBRATIONS OF STRUCTURES Genoa, Italy, June 12-16, 2000 BLUFF-BODY AERODYNAMICS Lecture Notes by Guido Buresti Department of Aerospace Engineering
More informationTheory of turbo machinery. Chapter 3
Theory of turbo machinery Chapter 3 D cascades Let us first understand the facts and then we may seek the causes. (Aristotle) D cascades High hub-tip ratio (of radii) negligible radial velocities D cascades
More informationAeroelastic Analysis Of Membrane Wings
Aeroelastic Analysis Of Membrane Wings Soumitra P. Banerjee and Mayuresh J. Patil Virginia Polytechnic Institute and State University, Blacksburg, Virginia 46-3 The physics of flapping is very important
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 informationImplementing a Partitioned Algorithm for Fluid-Structure Interaction of Flexible Flapping Wings within Overture
10 th Symposimum on Overset Composite Grids and Solution Technology, NASA Ames Research Center Moffett Field, California, USA 1 Implementing a Partitioned Algorithm for Fluid-Structure Interaction of Flexible
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 informationAerodynamic Investigation of a 2D Wing and Flows in Ground Effect
26 2 2009 3 CHINESE JOURNAL OF COMPUTATIONAL PHYSICS Vol. 26,No. 2 Mar., 2009 Article ID : 10012246 X(2009) 0220231210 Aerodynamic Investigation of a 2D Wing and Flows in Ground Effect YANG Wei, YANG Zhigang
More informationLarge-eddy simulations for wind turbine blade: rotational augmentation and dynamic stall
Large-eddy simulations for wind turbine blade: rotational augmentation and dynamic stall Y. Kim, I.P. Castro, and Z.T. Xie Introduction Wind turbines operate in the atmospheric boundary layer and their
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 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 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 informationChapter 6: Incompressible Inviscid Flow
Chapter 6: Incompressible Inviscid Flow 6-1 Introduction 6-2 Nondimensionalization of the NSE 6-3 Creeping Flow 6-4 Inviscid Regions of Flow 6-5 Irrotational Flow Approximation 6-6 Elementary Planar Irrotational
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 informationA simplified model for a small propeller with different airfoils along the blade
A simplified model for a small propeller with different airfoils along the blade Kamal A. R. Ismail 1) and *Célia V. A. G. Rosolen 2) 1), 2) State University of Campinas, Faculty of Mechanical Engineering,
More informationUnsteady Analysis of a Horizontal Axis Marine Current Turbine in Yawed Inflow Conditions With a Panel Method
First International Symposium on Marine Propulsors smp 09, Trondheim, Norway, June 2009 Unsteady Analysis of a Horizontal Axis Marine Current Turbine in Yawed Inflow Conditions With a J. Baltazar, J.A.C.
More informationNumerical Investigation of Laminar Flow over a Rotating Circular Cylinder
International Journal of Mechanical & Mechatronics Engineering IJMME-IJENS Vol:13 No:3 32 Numerical Investigation of Laminar Flow over a Rotating Circular Cylinder Ressan Faris Al-Maliky Department of
More informationVISCID/INVISCID INTERACTION ANALYSIS OF EJECTOR WINGS. M. Bevilaqua, C. J. Woan, and E. F. Schum
VISCID/INVISCID INTERACTION ANALYSIS OF EJECTOR WINGS OP. by M. Bevilaqua, C. J. Woan, and E. F. Schum 0 Rockwell International, North American Aircraft Division Columbus, Ohio - Paper Presented at Ejector
More informationCHAPTER 4 OPTIMIZATION OF COEFFICIENT OF LIFT, DRAG AND POWER - AN ITERATIVE APPROACH
82 CHAPTER 4 OPTIMIZATION OF COEFFICIENT OF LIFT, DRAG AND POWER - AN ITERATIVE APPROACH The coefficient of lift, drag and power for wind turbine rotor is optimized using an iterative approach. The coefficient
More informationSection III gives the authors intuitive version of their account. Taking their numbered points in order,
The classical theory of flight is one the most beautiful and subtle achievements of applied mathematics. However, it is no longer of interest to mathematicians, because they know that it has been a solved
More information