Chapter 13 Compressible Thin Airfoil Theory

Total Page:16

File Type:pdf, Size:1020Kb

Chapter 13 Compressible Thin Airfoil Theory CHAPTER 13 COMPRESSIBLE THIN AIRFOIL THEORY 13.1 COMPRESSIBLE POTENTIAL FLOW 13.1.1 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 from an inviscid, irrotational model of the flow. Recall the equations developed in Chapter 6 governing steady, irrotational, homentropic (¢s0= ) flow in the absence of body forces. ¢l()U = 0 UU ¢P ¢£¥--------------- +0-------- = ¤¦2 l . (13.1) P l a ------ = £¥------ ¤¦ P0 l0 The gradient of the isentropic relation is 2 ¢Pa= ¢l. (13.2) Recall from the development in Chapter 6 that P a – 1 ¢P ¢£¥--- = £¥------------ -------- (13.3) ¤¦l ¤¦a l Using (13.3) the momentum equation becomes a P UU ¢£¥£¥------------ --- + --------------- = 0 (13.4) ¤¦¤¦a – 1 l 2 Substitute (13.2) into the continuity equation and use (13.3). The continuity equa- tion becomes 2 2 U ¢a + ()a – 1 a ¢ U = 0 . (13.5) Equate the Bernoulli integral to free stream conditions. bjc 13.1 5/31/13 Compressible potential flow 2 2 2 2 a UU a' U' a' a – 1 2 ------------ + --------------- ==------------ + -------- ------------£¥1 + ------------M =C T (13.6) a – 1 2 a – 1 2 a – 1¤¦2 ' p t Note that the momentum equation is essentially equivalent to the statement that the stagnation temperature T t is constant throughout the flow. Using (13.6) we can write 2 a UU £¥------------ = h – --------------- (13.7) ¤¦a – 1 t 2 The continuity equation finally becomes UU UU ()a – 1 £¥h – --------------- ¢ UU– ¢£¥--------------- = 0 (13.8) ¤¦t 2 ¤¦2 The equations governing compressible, steady, inviscid, irrotational motion reduce to a single equation for the velocity vector U . The irrotationality condition ¢ × U = 0 permits the introduction of a velocity potential. U = ¢\ (13.9) and (13.8) becomes ¢\ ¢\ 2 ¢\ ¢\ ()a – 1 £¥h – ------------------------ ¢ \¢\¢– £¥------------------------ = 0 . (13.10) ¤¦t 2 ¤¦2 For complex body shapes numerical methods are normally used to solve for \ . However the equation is of relatively limited applicability. If the flow is over a thick airfoil or a bluff body for instance then the equation only applies to the sub- sonic Mach number regime at Mach numbers below the range where shocks begin to appear on the body. At high subsonic and supersonic Mach numbers where there are shocks then the homentropic assumption (13.2) breaks down. Equation (13.8) also applies to internal flows without shocks such as fully expanded nozzle flow. 13.1.2 THE NONLINEAR SMALL DISTURBANCE APPROXIMATION In the case of a thin airfoil that only slightly disturbs the flow, equation (13.8) can be simplified using small disturbance theory. 5/31/13 13.2 bjc Compressible potential flow Consider the flow past a thin 3-D airfoil shown below. v y=f(x,z) U' + u U' y x z The velocity field consists of a freestream flow plus a small disturbance UU= ' + u Vv= (13.11) Ww= where . (13.12) uU ' « 1 , v U' « 1 , w U' « 1 Similarly the state variables deviate only slightly from freestream values. PP= ' + P' TT= ' + T' (13.13) ll= ' + l' and . (13.14) aa= ' + a' This decomposition of variables is substituted into equation (13.8). Various terms are 2 2 2 2 UU U' u v w --------------- = --------uU++++----- ----- ------ 2 2 ' 2 2 2 (13.15) ¢ U = ux ++vy wz and bjc 13.3 5/31/13 Compressible potential flow UU ¢£¥--------------- = (u U +++uu vv ww , ¤¦2 x ' x x x (13.16) uyU' +++uuy vvy wwy, uzU' +++uuz vvz wwz ) as well as UU ()a – 1 £¥h – --------------- ¢ U = ¤¦t 2 2 2 2 2 £¥£¥U' u v w ()a – 1 ²´ht – ²´-------- ++++uU' ----- ----- ------ ux + ¤¦¤¦2 2 2 2 (13.17) 2 2 2 2 £¥£¥U' u v w ()a – 1 ²´ht – ²´-------- ++++uU' ----- ----- ------ vy + ¤¦¤¦2 2 2 2 2 2 2 2 £¥£¥U' u v w ()a – 1 ²´ht – ²´-------- ++++uU' ----- ----- ------ wz ¤¦¤¦2 2 2 2 and, finally UU 2 U ¢£¥--------------- = u U +++uu U vv U ww U + ¤¦2 x ' x ' x ' x ' 2 uu U +++u u uvv uww + x ' x x x . (13.18) 2 vuyU' +++vuuy v vy vwwy + 2 wuzU' +++wuuz wvvz w wz Neglect terms in (13.17) and (13.18) that are of third order in the disturbance velocities. Now 5/31/13 13.4 bjc Compressible potential flow UU UU ()a – 1 £¥h – --------------- ¢ UU– ¢£¥--------------- ¤¦t 2 ¤¦2 2 £¥£¥U' ()a – 1 ²´ht – ²´-------- + uU' ux + ¤¦¤¦2 2 £¥£¥U' (13.19) ()a – 1 ²´ht – ²´-------- + uU' vy + ¤¦¤¦2 2 £¥£¥U' 2 ()a – 1 ²´ht – ²´-------- + uU' wz – ()uxU' +++uuxU' vvxU' wwxU' – ¤¦¤¦2 uuxU' – vuyU' – wuzU' or UU UU ()a – 1 £¥h – --------------- ¢ UU– ¢£¥--------------- ¤¦t 2 ¤¦2 2 (13.20) ()a – 1 ()h' – uU' ()ux ++vy wz – ()uxU' +++uuxU' vvxU' wwxU' –uuxU' – vuyU' – wuzU' 2 Recall that ()a – 1 h' = a' . Equation (13.20) can be rearranged to read UU UU ()a – 1 £¥h – --------------- ¢ UU– ¢£¥--------------- ¤¦t 2 ¤¦2 2 2 (13.21) a'()ux ++vy wz – U'ux – ()a + 1 uuxU' – ()vvxU' + wwxU' ()a – 1 ()uvyU' + uwzU' – vuyU' – wuzU' 2 Divide through by a' . bjc 13.5 5/31/13 Compressible potential flow UU UU ()a – 1 £¥h – --------------- ¢ UU– ¢£¥--------------- ¤¦t 2 ¤¦2 2 ()a + 1 M' ()1M– ' ux ++vy wz – --------------------------- uux – (13.22) a' M' ---------()()a – 1 ()uvy + uwz ++vuy wuz ++vvx wwx a' This equation contains both linear and quadratic terms in the velocity disturbances and one might expect to be able to neglect the quadratic terms. But note that the first term becomes very small near M' = 1 . Thus in order to maintain the small disturbance approximation at transonic Mach numbers the uux term must be retained. The remaining quadratic terms are small at all Mach numbers and can be dropped. Finally the small disturbance equation is 2 ()a + 1 M' . (13.23) ()1M– ' ux ++vy wz – --------------------------- uux = 0 a' The velocity potential is written in terms of a freestream potential and a distur- bance potential . (13.24) \ = U'x + q()xyz,, The small disturbance equation in terms of the disturbance potential becomes 2 M' (13.25) ()1M– ' qxx ++qyy qzz = ()a + 1 ---------qxqxx a' Equation (13.25) is valid over the whole range of subsonic, transonic and super- sonic Mach numbers. 13.1.3 LINEARIZED POTENTIAL FLOW If we restrict our attention to subsonic and supersonic flow, staying away from Mach numbers close to one, the nonlinear term on the right side of (13.25) can be dropped and the small disturbance potential equation reduces to the linear wave equation. 5/31/13 13.6 bjc Compressible potential flow 2 .` qxx – ()qyy + qzz = 0 (13.26) 2 where ` = M' – 1 . In two dimensions 2 ` qxx – qyy = 0 . (13.27) The general solution of (13.27) can be expressed as a sum of two arbitrary functions q()xy, = Fx()– `y + Gx()+ `y . (13.28) 2 Note that if M' < 1 the 2-D linearized potential equation (13.27) is an elliptic equation that can be rescaled to form Laplace’s equation and (13.28) expresses the solution in terms of conjugate complex variables. In this case the subsonic flow can be analyzed using the methods of complex analysis. Presently we will restrict our attention to the supersonic case. The subsonic case is treated later in the chapter. 2 If M' > 1 then (13.27) is the 2-D wave equation and has solutions of hyperbolic type. Supersonic flow is analyzed using the fact that the properties of the flow are constant along the characteristic lines x ± `y = constan t. The figure below illustrates supersonic flow past a thin airfoil with several characteristics shown. Notice that in the linear approximation the characteristics are all parallel to one another and lie at the Mach angle µ' of the free stream. Information about the flow is carried in the value of the potential assigned to a given characteristic and in the spacing between characteristics for a given flow change. Right-leaning characteristics carry the information about the flow on the upper surface of the wing and left-leaning characteristics carry information about the flow on the lower surface. bjc 13.7 5/31/13 Compressible potential flow x – `y = constant y µ' x M' > 1 µ -1 ' = Sin (1/M') x +const`y = tan All properties of the flow, velocity, pressure, temperature, etc. are constant along the characteristics. Since disturbances only propagate along downstream running characteristics we can write the velocity potential for the upper and lower surfaces as q()xy, = Fx()– `y y>0 . (13.29) q()xy, = Gx()+ `y y<0 Let yfx= () define the coordinates of the upper surface of the wing and ygx= () define the lower surface. The full nonlinear boundary condition on the upper surface is v df ---- = ------ . (13.30) U dx yf= In the spirit of the thin airfoil approximation this boundary condition can be approximated by the linearized form v df -------- = ------ (13.31) U' dx y0= which we can write as 5/31/13 13.8 bjc Compressible potential flow ,q()xy, df -------------------- = U £¥------ (13.32) ,y '¤¦dx y0= or U' df F'()x = –-------- £¥------ . (13.33) ` ¤¦dx On the lower surface the boundary condition is U' dg G'()x = -------- £¥------ . (13.34) ` ¤¦dx In the thin airfoil approximation the airfoil itself is, in effect, collapsed to a line along the x-axis, the velocity potential is extended to the line y0= and the sur- face boundary condition is applied at y0= instead of at the physical airfoil surface. The entire effect of the airfoil on the flow is accounted for by the vertical velocity perturbation generated by the local slope of the wing. The linearized boundary condition is valid on 2-D thin wings and on 3-D wings of that are of “thin planar form”.
Recommended publications
  • Che 253M Experiment No. 2 COMPRESSIBLE GAS FLOW
    Rev. 8/15 AD/GW ChE 253M Experiment No. 2 COMPRESSIBLE GAS FLOW The objective of this experiment is to familiarize the student with methods for measurement of compressible gas flow and to study gas flow under subsonic and supersonic flow conditions. The experiment is divided into three distinct parts: (1) Calibration and determination of the critical pressure ratio for a critical flow nozzle under supersonic flow conditions (2) Calculation of the discharge coefficient and Reynolds number for an orifice under subsonic (non- choked) flow conditions and (3) Determination of the orifice constants and mass discharge from a pressurized tank in a dynamic bleed down experiment under (choked) flow conditions. The experimental set up consists of a 100 psig air source branched into two manifolds: the first used for parts (1) and (2) and the second for part (3). The first manifold contains a critical flow nozzle, a NIST-calibrated in-line digital mass flow meter, and an orifice meter, all connected in series with copper piping. The second manifold contains a strain-gauge pressure transducer and a stainless steel tank, which can be pressurized and subsequently bled via a number of attached orifices. A number of NIST-calibrated digital hand held manometers are also used for measuring pressure in all 3 parts of this experiment. Assorted pressure regulators, manual valves, and pressure gauges are present on both manifolds and you are expected to familiarize yourself with the process flow, and know how to operate them to carry out the experiment. A process flow diagram plus handouts outlining the theory of operation of these devices are attached.
    [Show full text]
  • Laws of Similarity in Fluid Mechanics 21
    Laws of similarity in fluid mechanics B. Weigand1 & V. Simon2 1Institut für Thermodynamik der Luft- und Raumfahrt (ITLR), Universität Stuttgart, Germany. 2Isringhausen GmbH & Co KG, Lemgo, Germany. Abstract All processes, in nature as well as in technical systems, can be described by fundamental equations—the conservation equations. These equations can be derived using conservation princi- ples and have to be solved for the situation under consideration. This can be done without explicitly investigating the dimensions of the quantities involved. However, an important consideration in all equations used in fluid mechanics and thermodynamics is dimensional homogeneity. One can use the idea of dimensional consistency in order to group variables together into dimensionless parameters which are less numerous than the original variables. This method is known as dimen- sional analysis. This paper starts with a discussion on dimensions and about the pi theorem of Buckingham. This theorem relates the number of quantities with dimensions to the number of dimensionless groups needed to describe a situation. After establishing this basic relationship between quantities with dimensions and dimensionless groups, the conservation equations for processes in fluid mechanics (Cauchy and Navier–Stokes equations, continuity equation, energy equation) are explained. By non-dimensionalizing these equations, certain dimensionless groups appear (e.g. Reynolds number, Froude number, Grashof number, Weber number, Prandtl number). The physical significance and importance of these groups are explained and the simplifications of the underlying equations for large or small dimensionless parameters are described. Finally, some examples for selected processes in nature and engineering are given to illustrate the method. 1 Introduction If we compare a small leaf with a large one, or a child with its parents, we have the feeling that a ‘similarity’ of some sort exists.
    [Show full text]
  • Aerodynamics Material - Taylor & Francis
    CopyrightAerodynamics material - Taylor & Francis ______________________________________________________________________ 257 Aerodynamics Symbol List Symbol Definition Units a speed of sound ⁄ a speed of sound at sea level ⁄ A area aspect ratio ‐‐‐‐‐‐‐‐ b wing span c chord length c Copyrightmean aerodynamic material chord- Taylor & Francis specific heat at constant pressure of air · root chord tip chord specific heat at constant volume of air · / quarter chord total drag coefficient ‐‐‐‐‐‐‐‐ , induced drag coefficient ‐‐‐‐‐‐‐‐ , parasite drag coefficient ‐‐‐‐‐‐‐‐ , wave drag coefficient ‐‐‐‐‐‐‐‐ local skin friction coefficient ‐‐‐‐‐‐‐‐ lift coefficient ‐‐‐‐‐‐‐‐ , compressible lift coefficient ‐‐‐‐‐‐‐‐ compressible moment ‐‐‐‐‐‐‐‐ , coefficient , pitching moment coefficient ‐‐‐‐‐‐‐‐ , rolling moment coefficient ‐‐‐‐‐‐‐‐ , yawing moment coefficient ‐‐‐‐‐‐‐‐ ______________________________________________________________________ 258 Aerodynamics Aerodynamics Symbol List (cont.) Symbol Definition Units pressure coefficient ‐‐‐‐‐‐‐‐ compressible pressure ‐‐‐‐‐‐‐‐ , coefficient , critical pressure coefficient ‐‐‐‐‐‐‐‐ , supersonic pressure coefficient ‐‐‐‐‐‐‐‐ D total drag induced drag Copyright material - Taylor & Francis parasite drag e span efficiency factor ‐‐‐‐‐‐‐‐ L lift pitching moment · rolling moment · yawing moment · M mach number ‐‐‐‐‐‐‐‐ critical mach number ‐‐‐‐‐‐‐‐ free stream mach number ‐‐‐‐‐‐‐‐ P static pressure ⁄ total pressure ⁄ free stream pressure ⁄ q dynamic pressure ⁄ R
    [Show full text]
  • A Semi-Hydrostatic Theory of Gravity-Dominated Compressible Flow
    Generated using version 3.2 of the official AMS LATEX template 1 A semi-hydrostatic theory of gravity-dominated compressible flow 2 Thomas Dubos ∗ IPSL-Laboratoire de Météorologie Dynamique, Ecole Polytechnique, Palaiseau, France Fabrice voitus CNRM-Groupe d’étude de l’Atmosphere Metéorologique, Météo-France, Toulouse, France ∗Corresponding author address: Thomas Dubos, LMD, École Polytechnique, 91128 Palaiseau, France. E-mail: [email protected] 1 3 Abstract 4 From Hamilton’s least action principle, compressible equations of motion with density diag- 5 nosed from potential temperature through hydrostatic balance are derived. Slaving density 6 to potential temperature suppresses the degrees of freedom supporting the propagation of 7 acoustic waves and results in a sound-proof system. The linear normal modes and dispersion 8 relationship for an isothermal state of rest on f- and β- planes are accurate from hydrostatic 9 to non-hydrostatic scales, except for deep internal gravity waves. Especially the Lamb wave 10 and long Rossby waves are not distorted, unlike with anelastic or pseudo-incompressible 11 systems. 12 Compared to similar equations derived by Arakawa and Konor (2009), the semi-hydrostatic 13 system derived here possesses an additional term in the horizontal momentum budget. This 14 term is an apparent force resulting from the vertical coordinate not being the actual height 15 of an air parcel, but its hydrostatic height, i.e. the hypothetical height it would have after 16 the atmospheric column it belongs to has reached hydrostatic balance through adiabatic 17 vertical displacements of air parcels. The Lagrange multiplier λ introduced in Hamilton’s 18 principle to slave density to potential temperature is identified as the non-hydrostatic ver- 19 tical displacement, i.e.
    [Show full text]
  • Brazilian 14-Xs Hypersonic Unpowered Scramjet Aerospace Vehicle
    ISSN 2176-5480 22nd International Congress of Mechanical Engineering (COBEM 2013) November 3-7, 2013, Ribeirão Preto, SP, Brazil Copyright © 2013 by ABCM BRAZILIAN 14-X S HYPERSONIC UNPOWERED SCRAMJET AEROSPACE VEHICLE STRUCTURAL ANALYSIS AT MACH NUMBER 7 Álvaro Francisco Santos Pivetta Universidade do Vale do Paraíba/UNIVAP, Campus Urbanova Av. Shishima Hifumi, no 2911 Urbanova CEP. 12244-000 São José dos Campos, SP - Brasil [email protected] David Romanelli Pinto ETEP Faculdades, Av. Barão do Rio Branco, no 882 Jardim Esplanada CEP. 12.242-800 São José dos Campos, SP, Brasil [email protected] Giannino Ponchio Camillo Instituto de Estudos Avançados/IEAv - Trevo Coronel Aviador José Alberto Albano do Amarante, nº1 Putim CEP. 12.228-001 São José dos Campos, SP - Brasil. [email protected] Felipe Jean da Costa Instituto Tecnológico de Aeronáutica/ITA - Praça Marechal Eduardo Gomes, nº 50 Vila das Acácias CEP. 12.228-900 São José dos Campos, SP - Brasil [email protected] Paulo Gilberto de Paula Toro Instituto de Estudos Avançados/IEAv - Trevo Coronel Aviador José Alberto Albano do Amarante, nº1 Putim CEP. 12.228-001 São José dos Campos, SP - Brasil. [email protected] Abstract. The Brazilian VHA 14-X S is a technological demonstrator of a hypersonic airbreathing propulsion system based on supersonic combustion (scramjet) to fly at Earth’s atmosphere at 30km altitude at Mach number 7, designed at the Prof. Henry T. Nagamatsu Laboratory of Aerothermodynamics and Hypersonics, at the Institute for Advanced Studies. Basically, scramjet is a fully integrated airbreathing aeronautical engine that uses the oblique/conical shock waves generated during the hypersonic flight, to promote compression and deceleration of freestream atmospheric air at the inlet of the scramjet.
    [Show full text]
  • No Acoustical Change” for Propeller-Driven Small Airplanes and Commuter Category Airplanes
    4/15/03 AC 36-4C Appendix 4 Appendix 4 EQUIVALENT PROCEDURES AND DEMONSTRATING "NO ACOUSTICAL CHANGE” FOR PROPELLER-DRIVEN SMALL AIRPLANES AND COMMUTER CATEGORY AIRPLANES 1. Equivalent Procedures Equivalent Procedures, as referred to in this AC, are aircraft measurement, flight test, analytical or evaluation methods that differ from the methods specified in the text of part 36 Appendices A and B, but yield essentially the same noise levels. Equivalent procedures must be approved by the FAA. Equivalent procedures provide some flexibility for the applicant in conducting noise certification, and may be approved for the convenience of an applicant in conducting measurements that are not strictly in accordance with the 14 CFR part 36 procedures, or when a departure from the specifics of part 36 is necessitated by field conditions. The FAA’s Office of Environment and Energy (AEE) must approve all new equivalent procedures. Subsequent use of previously approved equivalent procedures such as flight intercept typically do not need FAA approval. 2. Acoustical Changes An acoustical change in the type design of an airplane is defined in 14 CFR section 21.93(b) as any voluntary change in the type design of an airplane which may increase its noise level; note that a change in design that decreases its noise level is not an acoustical change in terms of the rule. This definition in section 21.93(b) differs from an earlier definition that applied to propeller-driven small airplanes certificated under 14 CFR part 36 Appendix F. In the earlier definition, acoustical changes were restricted to (i) any change or removal of a muffler or other component of an exhaust system designed for noise control, or (ii) any change to an engine or propeller installation which would increase maximum continuous power or propeller tip speed.
    [Show full text]
  • Potential Flow Theory
    2.016 Hydrodynamics Reading #4 2.016 Hydrodynamics Prof. A.H. Techet Potential Flow Theory “When a flow is both frictionless and irrotational, pleasant things happen.” –F.M. White, Fluid Mechanics 4th ed. We can treat external flows around bodies as invicid (i.e. frictionless) and irrotational (i.e. the fluid particles are not rotating). This is because the viscous effects are limited to a thin layer next to the body called the boundary layer. In graduate classes like 2.25, you’ll learn how to solve for the invicid flow and then correct this within the boundary layer by considering viscosity. For now, let’s just learn how to solve for the invicid flow. We can define a potential function,!(x, z,t) , as a continuous function that satisfies the basic laws of fluid mechanics: conservation of mass and momentum, assuming incompressible, inviscid and irrotational flow. There is a vector identity (prove it for yourself!) that states for any scalar, ", " # "$ = 0 By definition, for irrotational flow, r ! " #V = 0 Therefore ! r V = "# ! where ! = !(x, y, z,t) is the velocity potential function. Such that the components of velocity in Cartesian coordinates, as functions of space and time, are ! "! "! "! u = , v = and w = (4.1) dx dy dz version 1.0 updated 9/22/2005 -1- ©2005 A. Techet 2.016 Hydrodynamics Reading #4 Laplace Equation The velocity must still satisfy the conservation of mass equation. We can substitute in the relationship between potential and velocity and arrive at the Laplace Equation, which we will revisit in our discussion on linear waves.
    [Show full text]
  • Chapter 5 Dimensional Analysis and Similarity
    Chapter 5 Dimensional Analysis and Similarity Motivation. In this chapter we discuss the planning, presentation, and interpretation of experimental data. We shall try to convince you that such data are best presented in dimensionless form. Experiments which might result in tables of output, or even mul- tiple volumes of tables, might be reduced to a single set of curves—or even a single curve—when suitably nondimensionalized. The technique for doing this is dimensional analysis. Chapter 3 presented gross control-volume balances of mass, momentum, and en- ergy which led to estimates of global parameters: mass flow, force, torque, total heat transfer. Chapter 4 presented infinitesimal balances which led to the basic partial dif- ferential equations of fluid flow and some particular solutions. These two chapters cov- ered analytical techniques, which are limited to fairly simple geometries and well- defined boundary conditions. Probably one-third of fluid-flow problems can be attacked in this analytical or theoretical manner. The other two-thirds of all fluid problems are too complex, both geometrically and physically, to be solved analytically. They must be tested by experiment. Their behav- ior is reported as experimental data. Such data are much more useful if they are ex- pressed in compact, economic form. Graphs are especially useful, since tabulated data cannot be absorbed, nor can the trends and rates of change be observed, by most en- gineering eyes. These are the motivations for dimensional analysis. The technique is traditional in fluid mechanics and is useful in all engineering and physical sciences, with notable uses also seen in the biological and social sciences.
    [Show full text]
  • Cavitation Experimental Investigation of Cavitation Regimes in a Converging-Diverging Nozzle Willian Hogendoorn
    Cavitation Experimental investigation of cavitation regimes in a converging-diverging nozzle Willian Hogendoorn Technische Universiteit Delft Draft Cavitation Experimental investigation of cavitation regimes in a converging-diverging nozzle by Willian Hogendoorn to obtain the degree of Master of Science at the Delft University of Technology, to be defended publicly on Wednesday May 3, 2017 at 14:00. Student number: 4223616 P&E report number: 2817 Project duration: September 6, 2016 – May 3, 2017 Thesis committee: Prof. dr. ir. C. Poelma, TU Delft, supervisor Prof. dr. ir. T. van Terwisga, MARIN Dr. R. Delfos, TU Delft MSc. S. Jahangir, TU Delft, daily supervisor An electronic version of this thesis is available at http://repository.tudelft.nl/. Contents Preface v Abstract vii Nomenclature ix 1 Introduction and outline 1 1.1 Introduction of main research themes . 1 1.2 Outline of report . 1 2 Literature study and theoretical background 3 2.1 Introduction to cavitation . 3 2.2 Relevant fluid parameters . 5 2.3 Current state of the art . 8 3 Experimental setup 13 3.1 Experimental apparatus . 13 3.2 Venturi. 15 3.3 Highspeed imaging . 16 3.4 Centrifugal pump . 17 4 Experimental procedure 19 4.1 Venturi calibration . 19 4.2 Camera settings . 21 4.3 Systematic data recording . 21 5 Data and data processing 23 5.1 Videodata . 23 5.2 LabView data . 28 6 Results and Discussion 29 6.1 Analysis of starting cloud cavitation shedding. 29 6.2 Flow blockage through cavity formation . 30 6.3 Cavity shedding frequency . 32 6.4 Cavitation dynamics . 34 6.5 Re-entrant jet dynamics .
    [Show full text]
  • Introduction to Compressible Computational Fluid Dynamics James S
    Introduction to Compressible Computational Fluid Dynamics James S. Sochacki Department of Mathematics James Madison University [email protected] Abstract This document is intended as an introduction to computational fluid dynamics at the upper undergraduate level. It is assumed that the student has had courses through three dimensional calculus and some computer programming experience with numer- ical algorithms. A course in differential equations is recommended. This document is intended to be used by undergraduate instructors and students to gain an under- standing of computational fluid dynamics. The document can be used in a classroom or research environment at the undergraduate level. The idea of this work is to have the students use the modules to discover properties of the equations and then relate this to the physics of fluid dynamics. Many issues, such as rarefactions and shocks are left out of the discussion because the intent is to have the students discover these concepts and then study them with the instructor. The document is used in part of the undergraduate MATH 365 - Computation Fluid Dynamics course at James Madi- son University (JMU) and is part of the joint NSF Grant between JMU and North Carolina Central University (NCCU): A Collaborative Computational Sciences Pro- gram. This document introduces the full three-dimensional Navier Stokes equations. As- sumptions to these equations are made to derive equations that are accessible to un- dergraduates with the above prerequisites. These equations are approximated using finite difference methods. The development of the equations and finite difference methods are contained in this document. Software modules and their corresponding documentation in Fortran 90, Maple and Matlab can be downloaded from the web- site: http://www.math.jmu.edu/~jim/compressible.html.
    [Show full text]
  • A Novel Full-Euler Low Mach Number IMEX Splitting 1 Introduction
    Commun. Comput. Phys. Vol. 27, No. 1, pp. 292-320 doi: 10.4208/cicp.OA-2018-0270 January 2020 A Novel Full-Euler Low Mach Number IMEX Splitting 1, 2 2,3 1 Jonas Zeifang ∗, Jochen Sch ¨utz , Klaus Kaiser , Andrea Beck , Maria Luk´aˇcov´a-Medvid’ov´a4 and Sebastian Noelle3 1 IAG, Universit¨at Stuttgart, Pfaffenwaldring 21, DE-70569 Stuttgart, Germany. 2 Faculty of Sciences, Hasselt University, Agoralaan Gebouw D, BE-3590 Diepenbeek, Belgium. 3 IGPM, RWTH Aachen University, Templergraben 55, DE-52062 Aachen, Germany. 4 Institut f¨ur Mathematik, Johannes Gutenberg-Universit¨at Mainz, Staudingerweg 9, DE-55128 Mainz, Germany. Received 15 October 2018; Accepted (in revised version) 23 January 2019 Abstract. In this paper, we introduce an extension of a splitting method for singularly perturbed equations, the so-called RS-IMEX splitting [Kaiser et al., Journal of Scientific Computing, 70(3), 1390–1407], to deal with the fully compressible Euler equations. The straightforward application of the splitting yields sub-equations that are, due to the occurrence of complex eigenvalues, not hyperbolic. A modification, slightly changing the convective flux, is introduced that overcomes this issue. It is shown that the split- ting gives rise to a discretization that respects the low-Mach number limit of the Euler equations; numerical results using finite volume and discontinuous Galerkin schemes show the potential of the discretization. AMS subject classifications: 35L81, 65M08, 65M60, 76M45 Key words: Euler equations, low-Mach, IMEX Runge-Kutta, RS-IMEX. 1 Introduction The modeling of many processes in fluid dynamics requires the solution of the com- pressible Euler equations.
    [Show full text]
  • INCOMPRESSIBLE FLOW AERODYNAMICS 3 0 0 3 Course Category: Programme Core A
    COURSE CODE COURSE TITLE L T P C 1151AE107 INCOMPRESSIBLE FLOW AERODYNAMICS 3 0 0 3 Course Category: Programme core a. Preamble : The primary objective of this course is to teach students how to determine aerodynamic lift and drag over an airfoil and wing at incompressible flow regime by analytical methods. b. Prerequisite Courses: Fluid Mechanics c. Related Courses: Airplane Performance Compressible flow Aerodynamics Aero elasticity Flapping wing dynamics Industrial aerodynamics Transonic Aerodynamics d. Course Educational Objectives: To introduce the concepts of mass, momentum and energy conservation relating to aerodynamics. To make the student understand the concept of vorticity, irrotationality, theory of air foils and wing sections. To introduce the basics of viscous flow. e. Course Outcomes: Upon the successful completion of the course, students will be able to: Knowledge Level CO Course Outcomes (Based on revised Nos. Bloom’s Taxonomy) Apply the physical principles to formulate the governing CO1 K3 aerodynamics equations Find the solution for two dimensional incompressible inviscid CO2 K3 flows Apply conformal transformation to find the solution for flow over CO3 airfoils and also find the solutions using classical thin airfoil K3 theory Apply Prandtl’s lifting-line theory to find the aerodynamic CO4 K3 characteristics of finite wing Find the solution for incompressible flow over a flat plate using CO5 K3 viscous flow concepts f. Correlation of COs with POs: COs PO1 PO2 PO3 PO4 PO5 PO6 PO7 PO8 PO9 PO10 PO11 PO12 CO1 H L H L M H H CO2 H L H L M H H CO3 H L H L M H H CO4 H L H L M H H CO5 H L H L M H H H- High; M-Medium; L-Low g.
    [Show full text]