A Study of Solving Navier-Stokes Equations with a Finite Volume Method Based on Polygonal Unstructured Grids and the Computation

Total Page:16

File Type:pdf, Size:1020Kb

A Study of Solving Navier-Stokes Equations with a Finite Volume Method Based on Polygonal Unstructured Grids and the Computation A STUDY OF SOLVING NAVIER-STOKES EQUATIONS WITH A FINITE VOLUME METHOD BASED ON POLYGONAL UNSTRUCTURED GRIDS AND THE COMPUTATIONAL ANALYSIS OF GROUND VEHICLE AERO- DYNAMICS by JIANNAN TAN, B.S. A Dissertation In MECHNICAL ENGINEERING Submitted to the Graduate Faculty of Texas Tech University in Partial Fulfillment of the Requirements for the Degree of DOCTOR OF PHILOSOPHY Approved Siva Parameswaran, Ph.D. Chair of the Committee Jingzhou Yang, Ph.D. Zhaoming He, Ph.D. Sukalyan Bhattacharya, Ph.D. Xinzhong Chen, Ph.D. Fred Hartmeister, Dean of the Graduate School December, 2010 Copyright 2010, Jiannan Tan Texas Tech University, Jiannan Tan, December2010 ACKNOWLEDGMENTS Although the following dissertation is an individual work, I could never complete it without the help from a lot of people. Firstly I’d like to express my sincere gratitude to my advisor, Dr. Siva Parameswaran for his consistent guidance and encouragement in the past five years. His knowledge, expertise and enthusiasm on CFD have always been the driving force throughout my graduate research career at Texas Tech Univer- sity. I want to thank my committee professors, Dr. Jingzhou Yang, Dr. Zhaoming He, Dr. Sukalyan Bhattacharya and Dr. Xinzhong Chen for their kind support and help in the past few years. My thanks also extend to my lab mates: Dongdae Lee, Zixi Chen, Tengxiao Liu, Di- vya R. and Suranga Dharmarathne. Working with them has always been a pleasure, which makes my research progressive and my life cheerful and colorful. Finally, my very special thanks go to my family, especially my wife, Ou Zhang. This piece of work would never exist without her supervision and support on a daily basis. ii Texas Tech University, Jiannan Tan, December2010 TABLE OF CONTENTS ACKNOWLEDGMENTS………………………………………………………………………... ii ABSTRACT……………………………………………………………………………………….. v LIST OF TABLES……………………………………………………………………………….. vi LIST OF FIGURES……………………………………………………………………………... vii 1. INTRODUCTION .................................................................................................1 1.1 The Development of Computational Fluid Mechanics .............................1 1.2 Objectives ...............................................................................................2 1.3 Contents of Dissertation ..........................................................................3 2. THE DISCRETIZATION OF N-S EQUATIONS ......................................................4 2.1 Derivation and Description of the Navier-Stokes (N-S) Equations ..........4 2.2 N-S Equations under Cartesian Coordinating System in A 2-D Case .......5 2.3 Finite Volume Method ............................................................................5 2.4 The Discretization of 2-D N-S Equations ................................................6 2.4.1 Gauss’s Theorem in a 2-D Polygonal Cell Case ............................................... 6 2.4.2 Discretization of 2-D N-S Equations by casting Gauss’s Theorem ..................... 7 2.4.3 Important Variables of the Discretized N-S Equations ....................................... 8 3. GRID STRUCTURE AND GRID GENERATION ....................................................9 3.1 Structured and Unstructured Grids/Mesh ................................................9 3.2 Triangular, Quadrilateral and Polygonal Grids/Mesh ............................ 10 3.3 Cartesian Grids ..................................................................................... 11 3.4 Grid Generation .................................................................................... 12 3.5 The Data Structure of the Mesh ............................................................ 12 3.6 An Example: The Data Structure of a 2-D Mesh File from ANSYS® FLUENT® ............................................................................................................. 14 3.6.1 Point Information Block .................................................................................. 14 3.6.2 Face Information Block ................................................................................... 16 3.6.3 Cell Information Block .................................................................................... 18 3.6.4 User Defined Information Block ...................................................................... 18 3.6.5 The Internal Code and the Corresponding Boundary Condition ....................... 19 3.6.6 The Demonstrated Data structured of the Mesh in Figure 3.5 ......................... 19 iii Texas Tech University, Jiannan Tan, December2010 4. SOLVING THE DISCRETIZED N-S EQUATIONS ............................................... 21 4.1 The Goals of Solving the N-S Equations ............................................... 21 4.2 The Positions of the Flow Variables ...................................................... 21 4.3 The Iterative Method ............................................................................ 22 4.4 Handling Diffusive Term, Convective Term and Pressure Term ........... 23 4.4.1 The Diffusive Term ......................................................................................... 24 4.4.2 The Pressure Term and Least Square Method ............................................... 25 4.4.3 Convective Term and Momentum Interpolation for Face Velocity .................... 27 4.4.4 A New Momentum Interpolation Method ......................................................... 27 4.5 Conclusion............................................................................................ 39 5. CFD VALIDATION ............................................................................................ 40 5.1 Internal Flow in a Straight Channel ....................................................... 40 5.2 Internal Flow in a Z-Pipe ...................................................................... 42 5.3 Internal Flow In a Jet ............................................................................ 45 5.4 Laminar Flow Past a Square Cylinder ................................................... 49 5.5 Conclusion............................................................................................ 54 6. EFFECTS OF CROSS WIND ON SPORT UTILITY VEHICLES (SUV): A COMPUTATIONAL STUDY....................................................................................... 56 6.1 INTRODUCTION ................................................................................ 56 6.2 Aerodynamic Forces and coefficients.................................................... 57 6.2.1 Aerodynamic forces and coefficients .............................................................. 57 6.2.2 Method of determining yaw angle and the resultant wind velocity ................... 58 6.3 CAD models and cfd simulation ........................................................... 59 6.3.1 Vehicle CAD model ........................................................................................ 59 6.3.2 The wind tunnel geometry and boundary configuration ................................... 59 6.3.3 Blockage ratio ................................................................................................ 60 6.3.4 Method of dealing with different angles of attack ............................................. 61 6.3.5 Mesh generation and simulation ..................................................................... 62 6.4 Results and discussion .......................................................................... 63 6.5 Concluding remarks .............................................................................. 65 REFERENCES ........................................................................................................... 68 iv Texas Tech University, Jiannan Tan, December2010 ABSTRACT Navier-Stokes (N-S) equations describe the motion of fluid flow in the nature and they are called the governing equations of fluid flows. Solving Navier-Stokes equations is of great interest to the scientists and researchers. Due to the high nonlinearity, achiev- ing the analytical solutions for the N-S equations is extremely difficult, if not impossi- ble. Thus, people have to switch to numerical solutions with putting on certain restric- tions on the N-S equations. This leads to the development of Computational Fluid Dy- namics (CFD). This dissertation contains two major sections. The first section is about theoretical study of CFD. We go through the whole process that a CFD analysis normally re- quires: generating mesh, setting boundary conditions and achieving numerical solu- tions of N-S equations, and post-processing to achieve flow field plots. An in-house 2- D CFD code based on unstructured polygonal mesh is presented, in which a new mo- mentum interpolation method is developed and implemented to calculate the flow flux on the cell faces. The 2-D code is also validated by comparing the numerical results with widely-known analytical results, if available, or by benchmarking with the results produced by commercial CFD software packages. The second section of this disserta- tion is about one of the applications of CFD in modern auto industry – ground vehicle aerodynamics. The cross wind effect on a sport utility vehicle (SUV) is studied and analyzed using CFD methods and compared with available wind tunnel experimental results. The first section contributes to the philosophy of mechanical physics and the second section aims to fulfill the purpose of engineering. v Texas Tech University, Jiannan Tan, December2010 LIST OF TABLES Table 3.1 Table of Internal Code
Recommended publications
  • Sensitivity Analysis of Non-Linear Steep Waves Using VOF Method
    Tenth International Conference on ICCFD10-268 Computational Fluid Dynamics (ICCFD10), Barcelona, Spain, July 9-13, 2018 Sensitivity Analysis of Non-linear Steep Waves using VOF Method A. Khaware*, V. Gupta*, K. Srikanth *, and P. Sharkey ** Corresponding author: [email protected] * ANSYS Software Pvt Ltd, Pune, India. ** ANSYS UK Ltd, Milton Park, UK Abstract: The analysis and prediction of non-linear waves is a crucial part of ocean hydrodynamics. Sea waves are typically non-linear in nature, and whilst models exist to predict their behavior, limits exist in their applicability. In practice, as the waves become increasingly steeper, they approach a point beyond which the wave integrity cannot be maintained, and they 'break'. Understanding the limits of available models as waves approach these break conditions can significantly help to improve the accuracy of their potential impact in the field. Moreover, inaccurate modeling of wave kinematics can result in erroneous hydrodynamic forces being predicted. This paper investigates the sensitivity of non-linear wave modeling from both an analytical and a numerical perspective. Using a Volume of Fluid (VOF) method, coupled with the Open Channel Flow module in ANSYS Fluent, sensitivity studies are performed for a variety of non-linear wave scenarios with high steepness and high relative height. These scenarios are intended to mimic the near-break conditions of the wave. 5th order solitary wave models are applied to shallow wave scenarios with high relative heights, and 5th order Stokes wave models are applied to short gravity waves with high wave steepness. Stokes waves are further applied in the shallow regime at high wave steepness to examine the wave sensitivity under extreme conditions.
    [Show full text]
  • Waves and Structures
    WAVES AND STRUCTURES By Dr M C Deo Professor of Civil Engineering Indian Institute of Technology Bombay Powai, Mumbai 400 076 Contact: [email protected]; (+91) 22 2572 2377 (Please refer as follows, if you use any part of this book: Deo M C (2013): Waves and Structures, http://www.civil.iitb.ac.in/~mcdeo/waves.html) (Suggestions to improve/modify contents are welcome) 1 Content Chapter 1: Introduction 4 Chapter 2: Wave Theories 18 Chapter 3: Random Waves 47 Chapter 4: Wave Propagation 80 Chapter 5: Numerical Modeling of Waves 110 Chapter 6: Design Water Depth 115 Chapter 7: Wave Forces on Shore-Based Structures 132 Chapter 8: Wave Force On Small Diameter Members 150 Chapter 9: Maximum Wave Force on the Entire Structure 173 Chapter 10: Wave Forces on Large Diameter Members 187 Chapter 11: Spectral and Statistical Analysis of Wave Forces 209 Chapter 12: Wave Run Up 221 Chapter 13: Pipeline Hydrodynamics 234 Chapter 14: Statics of Floating Bodies 241 Chapter 15: Vibrations 268 Chapter 16: Motions of Freely Floating Bodies 283 Chapter 17: Motion Response of Compliant Structures 315 2 Notations 338 References 342 3 CHAPTER 1 INTRODUCTION 1.1 Introduction The knowledge of magnitude and behavior of ocean waves at site is an essential prerequisite for almost all activities in the ocean including planning, design, construction and operation related to harbor, coastal and structures. The waves of major concern to a harbor engineer are generated by the action of wind. The wind creates a disturbance in the sea which is restored to its calm equilibrium position by the action of gravity and hence resulting waves are called wind generated gravity waves.
    [Show full text]
  • THERMODYNAMICS, HEAT TRANSFER, and FLUID FLOW, Module 3 Fluid Flow Blank Fluid Flow TABLE of CONTENTS
    Department of Energy Fundamentals Handbook THERMODYNAMICS, HEAT TRANSFER, AND FLUID FLOW, Module 3 Fluid Flow blank Fluid Flow TABLE OF CONTENTS TABLE OF CONTENTS LIST OF FIGURES .................................................. iv LIST OF TABLES ................................................... v REFERENCES ..................................................... vi OBJECTIVES ..................................................... vii CONTINUITY EQUATION ............................................ 1 Introduction .................................................. 1 Properties of Fluids ............................................. 2 Buoyancy .................................................... 2 Compressibility ................................................ 3 Relationship Between Depth and Pressure ............................. 3 Pascal’s Law .................................................. 7 Control Volume ............................................... 8 Volumetric Flow Rate ........................................... 9 Mass Flow Rate ............................................... 9 Conservation of Mass ........................................... 10 Steady-State Flow ............................................. 10 Continuity Equation ............................................ 11 Summary ................................................... 16 LAMINAR AND TURBULENT FLOW ................................... 17 Flow Regimes ................................................ 17 Laminar Flow ...............................................
    [Show full text]
  • POTENTIAL FLOW: in Which the Vorticity Is Zero
    Chapter 1 Potential flow: 1.1 General formulation Inviscid and irrotational flows in the limit of high Reynolds number are referred to as ‘potential’ or ‘ideal’ flows. The term ‘inviscid’ refers to flows where viscous forces are small compared to inertial forces, so that the fluid viscosity can be neglected in comparison to fluid inertia. ‘Potential’ or ‘ideal’ flows are a class of inviscid flows in which the vorticity ω, which is the curl of the velocity vector, is zero, i. e. ∂uk ωi = ǫijk = 0 (1.1) ∂xj Since the curl of the velocity is zero, the velocity can be expressed as the gradient of a potential φ, ∂φ ui = (1.2) ∂xi and hence the name ‘potential flow’. Using equation 1.2 for the velocity field, the Navier-Stokes mass and mo- mentum equations can be written in terms of the velocity potential. The mass conservation equation, expressed in terms of the velocity potential, is, 2 ∂ui ∂ φ = 2 = 0 (1.3) ∂xi ∂xi Thus, the mass conservation equation simply states that the velocity potential satisfies the Laplace equation, 2φ = 0. Therefore, for potential flow we can use all the techiques developed∇ for solving the Laplace equation earlier. The momentum conservation equation an inviscid flow is, ∂ui ∂ui 1 ∂p fi + uj = + (1.4) ∂t ∂xj −ρ ∂xi ρ where fi is the body force per unit volume acting on the fluid. The second term on the right side of the above equation can be simplified for an irrotational flow 1 2 CHAPTER 1. POTENTIAL FLOW: in which the vorticity is zero.
    [Show full text]
  • A New Class of Exact Solutions of the Navier–Stokes Equations for Swirling Flows in Porous and Rotating Pipes
    Advances in Fluid Mechanics VIII 67 A new class of exact solutions of the Navier–Stokes equations for swirling flows in porous and rotating pipes A. Fatsis1, J. Statharas2, A. Panoutsopoulou3 & N. Vlachakis1 1Technological University of Chalkis, Department of Mechanical Engineering, Greece 2Technological University of Chalkis, Department of Aeronautical Engineering, Greece 3Hellenic Defence Systems, Greece Abstract Flow field analysis through porous boundaries is of great importance, both in engineering and bio-physical fields, such as transpiration cooling, soil mechanics, food preservation, blood flow and artificial dialysis. A new family of exact solution of the Navier–Stokes equations for unsteady laminar flow inside rotating systems of porous walls is presented in this study. The analytical solution of the Navier–Stokes equations is based on the use of the Bessel functions of the first kind. To resolve these equations analytically, it is assumed that the effect of the body force by mass transfer phenomena is the ‘porosity’ of the porous boundary in which the fluid moves. In the present study the effect of porous boundaries on unsteady viscous flow is examined for two different cases. The first one examines the flow between two rotated porous cylinders and the second one discusses the swirl flow in a rotated porous pipe. The results obtained reveal the predominant features of the unsteady flows examined. The developed solutions are of general application and can be applied to any swirling flow in porous axisymmetric rotating geometries.
    [Show full text]
  • An Essay on Lagrangian and Eulerian Kinematics of Fluid Flow
    Lagrangian and Eulerian Representations of Fluid Flow: Kinematics and the Equations of Motion James F. Price Woods Hole Oceanographic Institution, Woods Hole, MA, 02543 [email protected], http://www.whoi.edu/science/PO/people/jprice June 7, 2006 Summary: This essay introduces the two methods that are widely used to observe and analyze fluid flows, either by observing the trajectories of specific fluid parcels, which yields what is commonly termed a Lagrangian representation, or by observing the fluid velocity at fixed positions, which yields an Eulerian representation. Lagrangian methods are often the most efficient way to sample a fluid flow and the physical conservation laws are inherently Lagrangian since they apply to moving fluid volumes rather than to the fluid that happens to be present at some fixed point in space. Nevertheless, the Lagrangian equations of motion applied to a three-dimensional continuum are quite difficult in most applications, and thus almost all of the theory (forward calculation) in fluid mechanics is developed within the Eulerian system. Lagrangian and Eulerian concepts and methods are thus used side-by-side in many investigations, and the premise of this essay is that an understanding of both systems and the relationships between them can help form the framework for a study of fluid mechanics. 1 The transformation of the conservation laws from a Lagrangian to an Eulerian system can be envisaged in three steps. (1) The first is dubbed the Fundamental Principle of Kinematics; the fluid velocity at a given time and fixed position (the Eulerian velocity) is equal to the velocity of the fluid parcel (the Lagrangian velocity) that is present at that position at that instant.
    [Show full text]
  • Electromagnetic Field Orientation and Dynamics Governs Advection Characteristics Within Pendent Droplets
    Electromagnetic field orientation and dynamics governs advection characteristics within pendent droplets Purbarun Dhar a, *, Vivek Jaiswal a and A R Harikrishnan b a Department of Mechanical Engineering, Indian Institute of Technology Ropar, Rupnagar–14001, India b Department of Mechanical Engineering, Indian Institute of Technology Madras, Chennai–600036, India * Corresponding author: E–mail: [email protected] Phone: +91-1881-24-2173 Abstract The article reports the domineering governing role played by the direction of electric and magnetic fields on the internal advection pattern and strength within salt solution pendant droplets. Literature shows that solutal advection drives circulation cells within salt based droplets. Flow visualization and velocimetry reveals that the direction of the applied field governs the enhancement/reduction in circulation velocity and the directionality of circulation inside the droplet. Further, it is noted that while magnetic fields augment the circulation velocity, the electric field leads to deterioration of the same. The concepts of electro and magnetohydrodynamics are appealed to and a Stokesian stream function based mathematical model to deduce the field mediated velocities has been proposed. The model is found to reveal the roles of and degree of dependence on the governing Hartmann, Stuart, Reynolds and Masuda numbers. The theoretical predictions are observed to be in good agreement with experimental average spatio-temporal velocities. The present findings may have strong implications in microscale electro and/or magnetohydrodynamics. Keywords: Droplet; Particle Image Velocimetry; Electrohydrodynamics; Magnetohydrodynamics; Hartmann number; Masuda number Introduction – Fluid dynamics1, heat transfer2, 3 and species transport4, 5 dynamics in droplets, both pendant6, 7 as well as sessile8, 9, has been an area of immense academic importance within the research community.
    [Show full text]
  • On Vortex Shedding and Prediction of Vortex-Induced Vibrations of Circular Cylinders
    INEI-NO NTNU Trondheim Norges teknisk-naturvitenskapelige universitet — , Doktor ingenieravhandling 1997:44 1049 Institutt for marine konstraksjoner MTA-rapport 1997:117 cylinders On Vortex Shedding and Prediction of Vortex-Induced Vibrations of Circular Cylinders A thesis submitted in partial fulfillment of the requirements for the degree of doktor ingenipr by Karl H. Halse Trondheim, February 1997 Department of Marine Structures Faculty of Marine Technology Norwegian University of Science and Technology DISCLAIMER Portions of this document may be illegible in electronic image products. Images are produced from the best available original document. Abstract This work addresses the vortex shedding phenomenon and the problem of predicting vortex- induced vibrations of slender marine structures. This high-frequency response form may con ­ tribute significantly to the fatigue damage of a structure. As the development of hydrocarbons moves to deeper waters, the importance of accurately predicting the vortex-induced response has increased. Thus the need for proper response prediction methods is large. The problem of vortex-induced vibrations is a classical example of a fluid-structure interaction problem, which involves both fluid dynamics and structural mechanics. In this work an extensive review of existing research publications about vortex shedding from circular cylinders and the vortex-induced vibrations of cylinders is presented. Furthermore an overview of the different numerical approaches to modelling the fluid flow is given. A comparison study of existing response prediction methods revealed very large deviations be ­ tween the response predicted by the different methods. Differences were observed both in re­ sponse shapes and in vibration amplitudes. Most of the existing prediction methods are based upon simple force expressions, where the flow velocity is the only parameter representing the flow.
    [Show full text]
  • Simulations of Vortex Evolution and Phase Slip in Oscillatory Potential Flow of the Superfluid Component of 4He Through an Aperture
    Simulations of Vortex Evolution and Phase Slip in Oscillatory Potential Flow of the Superfluid Component of 4He Through an Aperture J. A. Flaten,1 C. T. Borden,2 C. A. Lindensmith,3 and W. Zimmermann, Jr.4 1Department of Physics, Luther College, Decorah, IA 52101, USA 2Michigan Technological University, Houghton, MI 49931, USA 3Jet Propulsion Laboratory, Pasadena, CA 91109, USA 4Tate Laboratory of Physics, University of Minnesota, Minneapolis, MN 55455, USA Abstract. The evolution of semicircular quantum vortex loops in oscillating potential flow emerging from an aperture is simulated in some highly symmetrical cases. As the frequency of potential flow oscillation increases, vortex loops that are evolving so as eventually to cross all of the streamlines of potential flow are drawn back toward the aperture when the flow reverses. As a result, the escape size of the vortex loops, and hence the net energy transferred from potential flow to vortex flow in such 2π phase-slip events, decreases as the oscillation frequency increases. Above some aperture-dependent and flow-dependent threshold frequency, vortex loops are drawn back into the aperture. Simulations are preformed using both radial potential flow and oblate-spheroidal potential flow. PACS numbers: 67.40.Vs, 47.15.Ki, 47.32.Cc I. INTRODUCTION When superfluid 4He potential flow through an aperture reaches the critical velocity associated with that aperture, energy begins to be dissipated from the flow. This energy dissipation is believed to be due to the motion of quantum vortices in the vicinity of the aperture. When a vortex core evolves in such a way as to cross streamlines of the potential flow through the aperture, energy is transferred from the potential flow field to the flow field of the vortex itself.
    [Show full text]
  • Fluid Mechanics Fluid Kinematics Velocity Field Flow Patterns
    www.getmyuni.com FLUID MECHANICS FLUID KINEMATICS Fluid Kinematics gives the geometry of fluid motion. It is a branch of fluid mechanics, which describes the fluid motion, and it’s consequences without consideration of the nature of forces causing the motion. Fluid kinematics is the study of velocity as a function of space and time in the flow field. From velocity, pressure variations and hence, forces acting on the fluid can be determined. VELOCITY FIELD Velocity at a given point is defined as the instantaneous velocity of the fluid particle, which at a given instant is passing through the point. It is represented by V=V(x,y,z,t). Vectorially, V=ui+vj+wk where u,v,w are three scalar components of velocity in x,y and z directions and (t) is the time. Velocity is a vector quantity and velocity field is a vector field. FLOW PATTERNS Fig. Flow Patterns Fluid Mechanics is a visual subject. Patterns of flow can be visualized in several ways. Basic types of line patterns used to visualize flow are streamline, path line, streak line and time line. (a) Stream line is a line, which is everywhere tangent to the velocity vector at a given instant. (b) Path line is the actual path traversed by a given particle. (c) Streak line is the locus of particles that have earlier passed through a prescribed point. (d) Time line is a set of fluid particles that form a line at a given instant. 1 Streamline is convenient to calculate mathematically. Other three lines are www.getmyuni.comeasier to obtain experimentally.
    [Show full text]
  • Navier–Stokes Equations
    Navier–Stokes equations In physics, the Navier–Stokes equations (/nævˈjeɪ stoʊks/), named after French engineer and physicist Claude-Louis Navier and British physicist and mathematician George Gabriel Stokes, describe the motion of viscous fluid substances. These balance equations arise from applying Isaac Newton's second law to fluid motion, together with the assumption that the stress in the fluid is the sum of a diffusing viscous term (proportional to the gradient of velocity) and a pressure term—hence describing viscous flow. The main difference between them and the simpler Euler equations for inviscid flow is that Navier–Stokes equations also factor in the Froude limit (no external field) and are not conservation equations, but rather a dissipative system, in the sense that they cannot be put into the quasilinear homogeneous form: Claude-Louis Navier George Gabriel Stokes The Navier–Stokes equations are useful because they describe the physics of many phenomena of scientific and engineering interest. They may be used to model the weather, ocean currents, water flow in a pipe and air flow around a wing. The Navier–Stokes equations, in their full and simplified forms, help with the design of aircraft and cars, the study of blood flow, the design of power stations, the analysis of pollution, and many other things. Coupled with Maxwell's equations, they can be used to model and study magnetohydrodynamics. The Navier–Stokes equations are also of great interest in a purely mathematical sense. Despite their wide range of practical uses, it has not yet been proven whether solutions always exist in three dimensions and, if they do exist, whether they are smooth – i.e.
    [Show full text]
  • A Review of Vortex Methods and Their Applications: from Creation to Recent Advances
    fluids Review A Review of Vortex Methods and Their Applications: From Creation to Recent Advances Chloé Mimeau * and Iraj Mortazavi Laboratory M2N, CNAM, 75003 Paris, France; [email protected] * Correspondence: [email protected]; Tel.: +33-0140-272-283 Abstract: This review paper presents an overview of Vortex Methods for flow simulation and their different sub-approaches, from their creation to the present. Particle methods distinguish themselves by their intuitive and natural description of the fluid flow as well as their low numerical dissipation and their stability. Vortex methods belong to Lagrangian approaches and allow us to solve the incompressible Navier-Stokes equations in their velocity-vorticity formulation. In the last three decades, the wide range of research works performed on these methods allowed us to highlight their robustness and accuracy while providing efficient computational algorithms and a solid mathematical framework. On the other hand, many efforts have been devoted to overcoming their main intrinsic difficulties, mostly relying on the treatment of the boundary conditions and the distortion of particle distribution. The present review aims to describe the Vortex methods by following their chronological evolution and provides for each step of their development the mathematical framework, the strengths and limits as well as references to applications and numerical simulations. The paper ends with a presentation of some challenging and very recent works based on Vortex methods and successfully applied to problems such as hydrodynamics, turbulent wake dynamics, sediment or porous flows. Keywords: hydrodynamics; incompressible flows; particle methods; remeshing; semi-Lagrangian vortex methods; vortex methods; Vortex-Particle-Mesh methods; vorticity; wake dynamics Citation: Mimeau, C.; Mortazavi, I.
    [Show full text]