Von Kármán Vortex Streets

Total Page:16

File Type:pdf, Size:1020Kb

Von Kármán Vortex Streets Von Karm´ an´ Vortex Streets: Simulating Natural Phenomena at Scale By: Sean Wang April 25th, 2019 1 Introduction The turn of the 19th century ushered in a golden era of fluid and aerodynamics research, with many scientists attacking the great questions that not only arose from advent of the powered airplane, but also those posed by the founders of fluid dynamics decades or centuries before. Possibly the greatest fluid dynamicist of all time, Ludwig Prandtl sought to answer of these questions, particularly as they related to flight forces. In his exploration of bluff body drag in near potential flows, Prandtl observed a natural phenomenon, an oscillatory sheet of shed vortices. This paper will explore this phenomenon as well as its outreaching applications in modelling at scale. 2 The Karm´ an´ Vortex Street The original understanding behind vortex shedding comes from a combination of Prandtl’s boundary layer theory (1904) and the existing potential flow theory. In potential flow theory, for an inviscid, incompressible (∇·U~ ), and irrotational (∇×U~ ) fluid, d’Alembert’s paradox states that there is zero drag. Experimentally, this is never true, with the closest examples being flows of ex- tremely low Reynolds numbers paired with a relatively slow free stream velocity (gen- erally M < :3). In the case of bluff bodies I.E. a cylinder in a flow, there is a strong adverse pressure gradient past the half chord, which causes almost immediate flow separation and a subsequent low pressure wake, integrating across the imbalance of pressures gives the resulting the drag. An additional consequence of separation is the phenomenon discussed earlier, alternating vortices caused by imbalances in the rotat- ing, separating wake; this is what is known as a Karm´ an´ street, after Theodore Van Karm´ an.´ Figure 1: A low pressure wake is visible on a simulated vortex street, a result of separation past the half chord. Re = 800, i = 6000 Von Karm´ an,´ Prandtl’s graduate student proposed idealizations for bluff bodies to predict vortex shedding in 1911, when he assumed that vortices were 2D, equidistant along their respective vortex sheets, and of equal strength whose circulation was given by: I G = Ud~ ~r (1) 1 Von Karm´ an´ went on to further generalize the case of symmetrically staggered vortices of equal circulation (or vorticity). Then the ratio of the distance between vortex sheets b and the longitudinal distance between vortices l can be described: b 1 p = cosh−1 2 ' :2806 (2) l p Von Karm´ an’s´ result is helpful both in predicting and verifying a pattern between vortices for bluff body vortex shedding, but is not an ideal foundation for simulation of such a vortex streets. Instead Navier-Stokes or Lattice-Boltzmann are used to do continuum and/or particle flow simulations. 3 The Lattice-Boltzmann Method The Lattice-Boltzmann method is a subset of fluid dynamics techniques typically used to calculate fluid flows at the microscopic scale. As opposed to macroscopic techniques like Navier-Stokes, which considers the aggregate quantities of a particular flow (den- sity, temperature, velocity), microscopic techniques approach flows as a conglomera- tion of particles instead of a continuum. Figure 2: An example of an octagonal lattice used in LBM. Adapted f rom Bao & Meskas 2011 Lattice-Boltzmann, more specifically, deals with fluids in a discrete lattice, with interactions occurring either at the boundaries or at nodes where particles meet. In this way Lattice-Boltzmann operates very much like gas/fluid lattice automata, with the key factor in vortex street simulation being the streaming and collision processes that occur in the model, defined by: 2 Collision Streaming z }| eq { z }| { [ fi(~x;t) − fi (~x;t)] fi(~x + c~eiDt;t + Dt) − fi(~x;t) = − (3) w eq ~ Where i represents some iterative value, w is a relaxation variable, and fi (x;t) is defined as, eq fi (~x;t) = wir + rsi(~u(~x;t)) (4) where wi are respective weights (0 at the walls: i = 1;2;3), and, e ·~u 9 (e ·~u)2 3 ~u ·~u s (~u) = w (3 i + i − ) (5) i i c 2 c2 2 c2 For a simulation of the Karm´ an´ vortex street, the relaxation factor is taken to be, 1 w = (6) 3nLB + :5 where U rcyl n = ¥ (7) LB Re With the foundations of Lattice-Boltzmann, in equations (3) through (7), a proper simulation can be constructed for the Karm´ an´ vortex street. The following is a code snippet defining non-iterative parameters considered in the vortex street simulation, 1 ###### Flow definition ################################## 2 iterPerFrame = 1000 ’’’Number of iterations before a new 3 frame is saved’’’ 4 maxIter = 50000 # Total number of time iterations 5 Re = 666.0 # Reynolds number :: start @ 220 6 nx = 800; ny = 180; ly=ny −1.0; q = 9 ’’’Lattice dimensions 7 and populations ’’’ 8 cx = nx/4; cy=ny/2; r=ny/9; # Coordinates of the cylinder. 9 uLB = 0 . 0 8 # Velocity in lattice units. 10 nulb = uLB∗ r/Re; omega = 1.0 / (3. ∗ nulb + 0 . 5 ) ; ’ ’ ’ Relax 11 ation parameter’’’ Next, the lattice constants are defined, and the code moves to iterating both stream- ing and colliding over maxIter. Velocities are subsequently stored in an array, then plotted for time intervals based on iterPerFrame (see Page 4, Appendix). 3 Figure 3: The development of the Karm´ an´ vortex street from separation wake, to instability, to periodic vortex street. Re = 800, ULB = :08, iterations ' [6000;46000] 4 4 Non-Dimensionalization for Flow over Volcanic Is- lands In terms of natural phenomena, volcanic islands are particularly interesting because of their particular penchant for shedding Karm´ an´ vortex streets. This should probably be no surprise; volcanic islands are usually isolated enough to encounter relatively uniform free stream flow, and their calderas can operate, at scale, like roughly circular bluff bodies. The Lattice-Boltzmann method can then be used model the vortex street from a given volcanic island. In order to do so, a scaling of non-dimensionalized parameters in the program code is necessary in order to obtain what is expected to be an accurate simulation of these shed vortices. Because Lattice-Boltzmann is sensitive to both high Re and compress- ibility, namely for higher free stream velocities, conservative scaling estimates are pref- erentially chosen to achieve model stability. Important scalings for modelling scale eddies, particularly as they relate to Reynolds number, include the following, r r p = alt (8) rm rSL T U m = ( m )2 (9) Tp Up L r T m = p ( m ):2 (10) Lp rm Tp and Reynold’s number, r U L r U L p p p = m m m (11) mp mm However, for the bulk of volcanic islands, particularly those of interest (alt ' 1800m), altitude variations in r, T, and m are of negligible order of magnitude dif- ference for vortex simulation. The respective proportional scalings are then assumed to be 1. Then comparing absolute scalings, assuming U¥ = 10m=s and LGranCanaria ' 47km, the simulation results can be scaled and compared to that of real observations. Conve- nient Re and ULB are chosen based upon the respective scalings. Then, Re m ULp ReULscale 10 s · 47km 3 = = m ' 587 · 10 (12) ReULm ReULsim (:04) s · 20m The simulation can be run and compare with the expected results. Note from the future: scalings were highly incorrect 5 5 Model Accuracy For Re = 350; ULB = :04, the simulated vortex street is as follows, Figure 4: The development of the Karm´ an´ vortex street from separation wake, to instability, to periodic vortex street. Re = 350, ULB = :04, times scales [10000;18000;20000;49000] 6 The particular island of interest for this simulation is Gran Canaria, in the Canary Island Chain, and is the island for which the scalings were calculated. Composing the simulated vortex street (Re = 350, ULB = :04, i = 49000, see last frame above) and the NASA satellite imagery side by side and marking distances, it is clear that the simulation models the scale vortices of Gran Canaria, quite well. Figure 5: The development of the Karm´ an´ vortex street from separation wake, to instability, to periodic vortex street. Re = 350, ULB = :04, times scales [10000;18000;20000;49000] A further discussion of results of Lattice-Boltzmann is continued on the next page. 7 LBM: Good LBM: Bad Given a correct scaling factor, accu- The calculated scaling factor was com- rately predicts the longitudinal spac- pletely incorrect for length L (probably ing between vortices down the Karm´ an´ my error) street – to a less extent the spacing of Does not accurately predict the long vortices across the street as well region of smooth trailing flow off Gran Simple and relatively fast (though not Canaria before the eddies on Python) Does not give information on relative It’s very nice vortex strengths (useful for anticipat- ing flight hazards) Does not follow von Karm´ an’s´ calcu- b lated l ratio (probably also my error) This particular simulation falls apart at high Re and ULB, a limitation of the code Table 1: A discussion of simulated results and empirical data. 6 Conclusion The phenomenon of Karm´ an´ vortex sheets as a model for scale phenomena, namely eddies off volcanic islands, is not only highly interesting, but carries implications into other applications. It is a theory that could govern or contribute intuition into other such phenomena, an example of which might be trailing tip vortices on commercial airliners, or perhaps the meteorological effects these streets may have. It is an area of fluid/aerodynamics still worthy of exploration today. In simulating these vortex streets, we found the Lattice-Boltzmann Method for fluid flows to be very effective in accurately predicting the geometries of these vortices.
Recommended publications
  • Kármán Vortex Street Energy Harvester for Picoscale Applications
    Kármán Vortex Street Energy Harvester for Picoscale Applications 22 March 2018 Team Members: James Doty Christopher Mayforth Nicholas Pratt Advisor: Professor Brian Savilonis A Major Qualifying Project submitted to the Faculty of WORCESTER POLYTECHNIC INSTITUTE in partial fulfilment of the requirements for the degree of Bachelor of Science This report represents work of WPI undergraduate students submitted to the faculty as evidence of a degree requirement. WPI routinely publishes these reports on its web site without editorial or peer review. For more information about the projects program at WPI, see http://www.wpi.edu/Academics/Projects. Cover Picture Credit: [1] Abstract The Kármán Vortex Street, a phenomenon produced by fluid flow over a bluff body, has the potential to serve as a low-impact, economically viable alternative power source for remote water-based electrical applications. This project focused on creating a self-contained device utilizing thin-film piezoelectric transducers to generate hydropower on a pico-scale level. A system capable of generating specific-frequency vortex streets at certain water velocities was developed with SOLIDWORKS modelling and Flow Simulation software. The final prototype nozzle’s velocity profile was verified through testing to produce a velocity increase from the free stream velocity. Piezoelectric testing resulted in a wide range of measured dominant frequencies, with corresponding average power outputs of up to 100 nanowatts. The output frequencies were inconsistent with predicted values, likely due to an unreliable testing environment and the complexity of the underlying theory. A more stable testing environment, better verification of the nozzle velocity profile, and fine-tuning the piezoelectric circuit would allow for a higher, more consistent power output.
    [Show full text]
  • A Concept of the Vortex Lift of Sharp-Edge Delta Wings Based on a Leading-Edge-Suction Analogy Tech Library Kafb, Nm
    I A CONCEPT OF THE VORTEX LIFT OF SHARP-EDGE DELTA WINGS BASED ON A LEADING-EDGE-SUCTION ANALOGY TECH LIBRARY KAFB, NM OL3042b NASA TN D-3767 A CONCEPT OF THE VORTEX LIFT OF SHARP-EDGE DELTA WINGS BASED ON A LEADING-EDGE-SUCTION ANALOGY By Edward C. Polhamus Langley Research Center Langley Station, Hampton, Va. NATIONAL AERONAUTICS AND SPACE ADMINISTRATION For sale by the Clearinghouse for Federal Scientific and Technical Information Springfield, Virginia 22151 - Price $1.00 A CONCEPT OF THE VORTEX LIFT OF SHARP-EDGE DELTA WINGS BASED ON A LEADING-EDGE-SUCTION ANALOGY By Edward C. Polhamus Langley Research Center SUMMARY A concept for the calculation of the vortex lift of sharp-edge delta wings is pre­ sented and compared with experimental data. The concept is based on an analogy between the vortex lift and the leading-edge suction associated with the potential flow about the leading edge. This concept, when combined with potential-flow theory modified to include the nonlinearities associated with the exact boundary condition and the loss of the lift component of the leading-edge suction, provides excellent prediction of the total lift for a wide range of delta wings up to angles of attack of 20° or greater. INTRODUCTION The aerodynamic characteristics of thin sharp-edge delta wings are of interest for supersonic aircraft and have been the subject of theoretical and experimental studies for many years in both the subsonic and supersonic speed ranges. Of particular interest at subsonic speeds has been the formation and influence of the leading-edge separation vor­ tex that occurs on wings having sharp, highly swept leading edges.
    [Show full text]
  • Three Types of Horizontal Vortices Observed in Wildland Mass And
    1624 JOURNAL OF CLIMATE AND APPLIED METEOROLOGY VOLUME26 Three Types of Horizontal Vortices Observed in Wildland Mas~ and Crown Fires DoNALD A. HAINES U.S. Department ofAgriculture, Forest Service, North Central Forest Experiment Station, East Lansing, Ml 48823 MAHLON C. SMITH Department ofMechanical Engineering, Michigan State University, East Lansing, Ml 48824 (Manuscript received 25 October 1986, in final form 4 May 1987) ABSTRACT Observation shows that three types of horizontal vortices may form during intense wildland fires. Two of these vortices are longitudinal relative to the ambient wind and the third is transverse. One of the longitudinal types, a vortex pair, occurs with extreme heat and low to moderate wind speeds. It may be a somewhat common structure on the flanks of intense crown fires when burning is concentrated along the fire's perimeter. The second longitudinal type, a single vortex, occurs with high winds and can dominate the entire fire. The third type, the transverse vortex, occurs on the upstream side of the convection column during intense burning and relatively low winds. These vortices are important because they contribute to fire spread and are a threat to fire fighter safety. This paper documents field observations of the vortices and supplies supportive meteorological and fuel data. The discussion includes applicable laboratory and conceptual studies in fluid flow and heat transfer that may apply to vortex formation. 1. Introduction experiments showed that when air flowed parallel to a heated metal ribbon that simulated the flank of a crown The occurrence of vertical vortices in wildland fires fire, a thin, buoyant plume capped with a vortex pair has been well documented as well as mathematically developed above the ribbon along its length.
    [Show full text]
  • Hurricane Vortex Dynamics During Atlantic Extratropical Transition
    714 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 65 Hurricane Vortex Dynamics during Atlantic Extratropical Transition CHRISTOPHER A. DAVIS National Center for Atmospheric Research,* Boulder, Colorado ϩ SARAH C. JONES AND MICHAEL RIEMER Universität Karlsruhe, Forschungszentrum Karlsruhe, Karlsruhe, Germany (Manuscript received 2 April 2007, in final form 6 July 2007) ABSTRACT Simulations of six Atlantic hurricanes are diagnosed to understand the behavior of realistic vortices in varying environments during the process of extratropical transition (ET). The simulations were performed in real time using the Advanced Research Weather Research and Forecasting (WRF) model (ARW), using a moving, storm-centered nest of either 4- or 1.33-km grid spacing. The six simulations, ranging from 45 to 96 h in length, provide realistic evolution of asymmetric precipitation structures, implying control by the synoptic scale, primarily through the vertical wind shear. The authors find that, as expected, the magnitude of the vortex tilt increases with increasing shear, but it is not until the shear approaches 20 m sϪ1 that the total vortex circulation decreases. Furthermore, the total vertical mass flux is proportional to the shear for shears less than about 20–25 m sϪ1, and therefore maximizes, not in the tropical phase, but rather during ET. This has important implications for predicting hurricane-induced perturbations of the midlatitude jet and its consequences on downstream predictability. Hurricane vortices in the sample resist shear by either adjusting their vertical structure through preces- sion (Helene 2006), forming an entirely new center (Irene 2005), or rapidly developing into a baroclinic cyclone in the presence of a favorable upper-tropospheric disturbance (Maria 2005).
    [Show full text]
  • 04 Delta Wings
    ExperimentalExperimental AerodynamicsAerodynamics Lecture 4: Delta wing experiments G. Dimitriadis Experimental Aerodynamics Introduction •! In this course we will demonstrate the use of several different experimental aerodynamic methodologies •! The particular application will be the aerodynamics of Delta wings at low airspeeds. •! Delta wings are of particular interest because of their lift generation mechanism. Experimental Aerodynamics Delta wing history •! Until the 1930s the vast majority of aircraft featured rectangular, trapezoidal or elliptical wings. •! Delta wings started being studied in the 1930s by Alexander Lippisch in Germany. •! Lippisch wanted to create tail-less aircraft, and Delta wings were one of the solutions he proposed. Experimental Aerodynamics Delta Lippisch DM-1 Designed as an interceptor jet but never produced. The photos show a glider prototype version. Experimental Aerodynamics High speed flight •! After the war, the potential of Delta wings for supersonic flight was recognized both in the US and the USSR. MiG-21 Convair XF-92 Experimental Aerodynamics Low speed performance •! Although Delta wings are designed for high speeds, they still have to take off and land at small airspeeds. •! It is important to determine the aerodynamic forces acting on Delta wings at low speed. •! The lift generated by such wings are low speeds can be split into two contributions: –! Potential flow lift –! Vortex lift Experimental Aerodynamics Delta wing geometry cb Wing surface: S = 2 2b Aspect ratio: AR = "! c c! b AR Sweep angle: tan ! = = 2c 4 b/2! Experimental Aerodynamics Potential flow lift •! Slender wing theory •! The wind is discretized into transverse segments. •! The flow around each segment is modeled as a 2D flow past a flat plate perpendicular to the free stream Experimental Aerodynamics Slender wing theory •! The problem of calculating the flow around the wing becomes equivalent to calculating the flow around each 2D segment.
    [Show full text]
  • Fluids – Lecture 16 Notes 1
    Fluids – Lecture 16 Notes 1. Vortex 2. Lifting flow about circular cylinder Reading: Anderson 3.14 – 3.16 Vortex Flowfield Definition A vortex flow has the following radial and tangential velocity components C V = 0 , V = r θ r where C is a scaling constant. The circulation around any closed circuit is computed as θ2 C Γ ≡ − V~ · d~s = − Vθ r dθ = − r dθ = −C (θ2 − θ1) I I Zθ1 r y y V r dθ ds dθ x x The integration range θ2−θ1 =2π if the circuit encircles the origin, but is zero otherwise. −2πC , (circuit encircles origin) Γ = ( 0 , (circuit doesn’t encircle origin) y y θ1 θ2 θ1 θ2 x x In lieu of C, it is convenient to redefine the vortex velocity field directly in terms of the circulation of any circuit which encloses the vortex origin. Γ V = − θ 2πr 1 A positive Γ corresponds to clockwise flow, while a negative Γ corresponds to counterclock- wise flow. Cartesian representation The cartesian velocity components of the vortex are Γ y u(x, y) = 2π x2 + y2 Γ x v(x, y) = − 2π x2 + y2 and the corresponding potential and stream functions are as follows. Γ Γ φ(x, y) = − arctan(y/x) = − θ 2π 2π Γ Γ ψ(x, y) = ln x2 + y2 = ln r 2π 2π q Singularity As with the source and doublet, the origin location (0, 0) is called a singular point of the vortex flow. The magnitude of the tangential velocity tends to infinity as 1 V ∼ θ r Hence, the singular point must be located outside the flow region of interest.
    [Show full text]
  • The Multiple Vortex Structure of a Tornado
    THE MULTIPLE VORTEX STRUCTURE OF A TORNADO Joshua Wurman University of Oklahoma 1. Introduction Conceptual models (Davies-Jones 1976), and 1 km computer (Rotunno 1984, Lewellen 1993, Lewellen et al 1997, Fiedler 1998) and laboratory simulations (Ward 1972, Church and Snow 1993, Church et al. 1979) of tornado structure predict that, under certain conditions, a primary tornado vortex will break down into several sub-tornadic scale multiple vortices (hereinafter “multiple-vortices”). Multiple vortices have been observed visually and in patterns of damage for decades (Fujita 1970, Pauley and Snow 1988). Direct radar evidence of sub-tornado- scale wind maxima probably associated with multiple-vortices was first obtained in a large tornado Figure 1.Doppler Velocity and Power in Spencer SD that destroyed much of the small town of Spencer, SD in tornado illustrating multiple vortices (black ovals). 1998 (Wurman 1999, Fig. 1). The tornado was observed W W ’ ’ at ranges of 1.7-5 km (to center). The vortices were Fig 2. Path of torn near Mulhall4 3 2 2 ° ° Tornado Core Region OK. Core flow diameter shown.7 7 weaker than the tornadic flow, with perturbations of ~ 9 9 20-30 ms-1 on a parent flow of ~ 85-95 ms-1. The 25m and vortices caused a degradation of the typical clear eye 37.5m non- 03:17:38 structure observed frequently in received power data oversampled 03:16:38 obtained in tornadoes not exhibiting strong multiple range vortices (Wurman and Gill 2000, Burgess et al 2001. resolution. 03:15:27 3x° xx’ N Sometimes, several “eyes” were evident.
    [Show full text]
  • A Vortex Model for Forces and Moments on Low-Aspect-Ratio Wings in Side-Slip with Experimental Validation Adam Devoria, Kamran Mohseni
    A vortex model for forces and moments on low-aspect-ratio wings in side-slip with experimental validation Adam Devoria, Kamran Mohseni To cite this version: Adam Devoria, Kamran Mohseni. A vortex model for forces and moments on low-aspect-ratio wings in side-slip with experimental validation. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, Royal Society, The, 2017, 473 (2198), 10.1098/rspa.2016.0760. hal-01918593 HAL Id: hal-01918593 https://hal.archives-ouvertes.fr/hal-01918593 Submitted on 11 Nov 2018 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. A vortex model for forces and rspa.royalsocietypublishing.org moments on low-aspect-ratio wings in side-slip with experimental validation Research Adam C. DeVoria1 and Kamran Mohseni1,2 Cite this article: DeVoria AC, Mohseni K. 2017 A vortex model for forces and moments on 1Department of Mechanical and Aerospace Engineering, and low-aspect-ratio wings in side-slip with 2Department of Electrical and Computer Engineering, experimental validation. Proc.R.Soc.A473: University of Florida, Gainesville, FL 32611, USA 20160760. http://dx.doi.org/10.1098/rspa.2016.0760 KM, 0000-0002-1382-221X This paper studies low-aspect-ratio (A) rectangular Received: 7 October 2016 wings at high incidence and in side-slip.
    [Show full text]
  • A Revised Tornado Definition and Changes in Tornado Taxonomy
    1256 WEATHER AND FORECASTING VOLUME 29 A Revised Tornado Definition and Changes in Tornado Taxonomy ERNEST M. AGEE Department of Earth, Atmospheric, and Planetary Sciences, Purdue University, West Lafayette, Indiana (Manuscript received 4 June 2014, in final form 30 July 2014) ABSTRACT The tornado taxonomy presented by Agee and Jones is revised to account for the new definition of a tor- nado provided by the American Meteorological Society (AMS) in October 2013, resulting in the elimination of shear-driven vortices from the taxonomy, such as gustnadoes and vortices in the eyewall of hurricanes. Other relevant research findings since the initial issuance of the taxonomy are also considered and in- corporated, where appropriate, to help improve the classification system. Multiple misoscale shear-driven vortices in a single tornado event, when resulting from an inertial instability, are also viewed to not meet the definition of a tornado. 1. Introduction and considerations from a cumuliform cloud, and often visible as a funnel cloud and/or circulating debris/dust at the ground.’’ In The first proposed tornado taxonomy was presented view of the latest definition, a few changes are warranted by Agee and Jones (2009, hereafter AJ) consisting of in the AJ taxonomy. Considering the roles played by three types and 15 species, ranging from the type I buoyancy and shear on a variety of spatial and temporal (potentially strong and violent) tornadoes produced by scales (from miso to meso to synoptic), coupled with the the classic supercell, to the more benign type III con- requirement in the latest definition that a tornado must vective and shear-driven vortices such as landspouts and be pendant from a cumuliform cloud, it is necessary to gustnadoes.
    [Show full text]
  • Vortex Dynamics Resulting from the Interaction Between Two NACA 23 012 Airfoils Thierry M
    Vortex dynamics resulting from the interaction between two NACA 23 012 airfoils Thierry M. Faure, Laurent Hétru, Olivier Montagnier To cite this version: Thierry M. Faure, Laurent Hétru, Olivier Montagnier. Vortex dynamics resulting from the interaction between two NACA 23 012 airfoils. 50th 3AF International Conference on Applied Aerodynamics, Association Aéronautique et Astronautique de France - 3AF, Mar 2015, Toulouse, France. hal- 01157857 HAL Id: hal-01157857 https://hal.archives-ouvertes.fr/hal-01157857 Submitted on 28 May 2015 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. 50h 3AF International Conference on Applied Aerodynamics 29-30 March – 01 April 2015, Toulouse - France FP10-2015-faure Vortex dynamics resulting from the interaction between two NACA 23 012 airfoils Thierry M. Faure (1,2), Laurent Hétru (1), Olivier Montagnier (1,3) (1) Centre de Recherche de l’Armée de l’air, École de l’Air, B.A. 701, 13661 Salon-de-Provence, France, Email: [email protected], [email protected], [email protected] (2) Laboratoire d’Informatique pour la Mécanique et les Sciences de l’Ingénieur, bâtiment 508, rue J. Von Neumann 91403 Orsay, France (3) Laboratoire de Mécanique et d’Acoustique, 31 chemin Joseph Aiguier, 13402 Marseille, France ABSTRACT An experimental study of the interaction between two airfoils, corresponding to a T-tail aircraft configuration, is implemented in a wind tunnel for a range of medium Reynolds numbers.
    [Show full text]
  • NUMERICAL INVESTIGATIONS of a TORNADO VORTEX USING VORTICITY CONFINEMENT HOLLY C. HASSENZAHL a Thesis Submitted in Partial
    NUMERICAL INVESTIGATIONS OF A TORNADO VORTEX USING VORTICITY CONFINEMENT by HOLLY C. HASSENZAHL A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science (Atmospheric and Oceanic Sciences) at the UNIVERSITY OF WISCONSIN MADISON 2007 i ABSTRACT The numerical simulation of a realistically strong tornado vortex and its associated condensation funnel has proven to be very difficult to resolve in atmospheric modeling. Many have attributed this failure to insufficient resolution of the models being used. Others have conjectured that the problem lies in the fact that strong gradients are eroded by numerical diffusion, thus prohibiting the formation of strong vortices. This latter hypothesis led engineers Steinhoff and Underhill (1994) to conceive the Vorticity Confinement (VC) technique, in an effort to restore the vorticity gradients lost to diffusion. In this study, the University of Wisconsin Non-Hydrostatic Modeling System (UW-NMS) is used to investigate the aforementioned hypotheses on a three-dimensional extension of the Wicker and Wilhelmson (1995) tornado vortex. These idealized simulations are carried out with two-way interactive nested grids at horizontal resolutions of 24 m and 12 m. Simulations without the VC technique do produce tornado vortices at both resolutions, however they are too weak to form condensation funnels extending to the surface. Comparisons with simulations employing the VC technique show that a realistically strong tornado vortex is resolved at the 24 m resolution, producing a beautiful condensation funnel that descends to the ground. However, when extended to a resolution of 12 m, the VC technique fails to converge to the real solution.
    [Show full text]
  • Circulation and Vorticity
    Lecture 4: Circulation and Vorticity • Circulation •Bjerknes Circulation Theorem • Vorticity • Potential Vorticity • Conservation of Potential Vorticity ESS227 Prof. Jin-Yi Yu Measurement of Rotation • Circulation and vorticity are the two primary measures of rotation in a fluid. • Circul ati on, whi ch i s a scal ar i nt egral quantit y, i s a macroscopic measure of rotation for a finite area of the fluid. • Vorticity, however, is a vector field that gives a microscopic measure of the rotation at any point in the fluid. ESS227 Prof. Jin-Yi Yu Circulation • The circulation, C, about a closed contour in a fluid is defined as thliihe line integral eval uated dl along th e contour of fh the component of the velocity vector that is locally tangent to the contour. C > 0 Î Counterclockwise C < 0 Î Clockwise ESS227 Prof. Jin-Yi Yu Example • That circulation is a measure of rotation is demonstrated readily by considering a circular ring of fluid of radius R in solid-body rotation at angular velocity Ω about the z axis . • In this case, U = Ω × R, where R is the distance from the axis of rotation to the ring of fluid. Thus the circulation about the ring is given by: • In this case the circulation is just 2π times the angular momentum of the fluid ring about the axis of rotation. Alternatively, note that C/(πR2) = 2Ω so that the circulation divided by the area enclosed by the loop is just twice thlhe angular spee dfifhid of rotation of the ring. • Unlike angular momentum or angular velocity, circulation can be computed without reference to an axis of rotation; it can thus be used to characterize fluid rotation in situations where “angular velocity” is not defined easily.
    [Show full text]