Turbulent Transition Simulation and Particulate Capture Modeling with an Incompressible Lattice Boltzmann Method

Total Page:16

File Type:pdf, Size:1020Kb

Turbulent Transition Simulation and Particulate Capture Modeling with an Incompressible Lattice Boltzmann Method Michigan Technological University Digital Commons @ Michigan Tech Dissertations, Master's Theses and Master's Reports 2017 TURBULENT TRANSITION SIMULATION AND PARTICULATE CAPTURE MODELING WITH AN INCOMPRESSIBLE LATTICE BOLTZMANN METHOD John R. Murdock Michigan Technological University, [email protected] Copyright 2017 John R. Murdock Recommended Citation Murdock, John R., "TURBULENT TRANSITION SIMULATION AND PARTICULATE CAPTURE MODELING WITH AN INCOMPRESSIBLE LATTICE BOLTZMANN METHOD", Open Access Dissertation, Michigan Technological University, 2017. https://doi.org/10.37099/mtu.dc.etdr/525 Follow this and additional works at: https://digitalcommons.mtu.edu/etdr Part of the Applied Mechanics Commons, Computer-Aided Engineering and Design Commons, and the Heat Transfer, Combustion Commons TURBULENT TRANSITION SIMULATION AND PARTICULATE CAPTURE MODELING WITH AN INCOMPRESSIBLE LATTICE BOLTZMANN METHOD By John R. Murdock A DISSERTATION Submitted in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY In Mechanical Engineering-Engineering Mechanics MICHIGAN TECHNOLOGICAL UNIVERSITY 2017 © 2017 John R. Murdock This dissertation has been approved in partial fulfillment of the requirements for the Degree of DOCTOR OF PHILOSOPHY in Mechanical Engineering-Engineering Mechanics. Department of Mechanical Engineering-Engineering Mechanics Dissertation Advisor: Dr. Song-Lin Yang Committee Member: Dr. Franz X. Tanner Committee Member: Dr. Kazuya Tajiri Committee Member: Dr. Youngchul Ra Department Chair: Dr. William W. Predebon Dedication To my wife who provided the utmost support for this endeavor from the outset - without which I would neither be who I am nor would this work be what it is today. Contents List of Figures ................................. xv List of Tables .................................. xxiii Preface ...................................... xxv Acknowledgments ............................... xxvii Nomenclature .................................. xxix Abstract ..................................... xxxv 1 Introduction ................................. 1 1.1 Transition . 3 1.2 Filtration . 4 1.3 Research Objectives and Structure . 8 1.3.1 Objectives . 8 1.3.2 Document Structure . 9 2 Literature Review ............................. 11 vii 2.1 Incompressible Lattice Boltzmann Method . 11 2.2 Turbulence Simulation . 13 2.3 Particulate Matter Transport . 15 2.3.1 LBE-Lagrange . 16 2.3.2 LBE-LGA . 16 2.3.3 LBE-Euler for Particle Concentration . 18 2.4 Dominant Mechanisms . 19 2.5 Heat Transfer . 19 2.5.1 DDF-LBE . 20 2.5.2 Conjugate Heat Transfer . 21 2.6 Porous and Disordered Media . 22 2.6.1 REV Scale . 23 2.6.2 Pore Scale . 24 2.7 Aftertreatment Synthesis . 25 2.8 Summary . 27 3 Methodology ................................ 29 3.1 The Lattice Boltzmann Method . 29 3.1.1 Origins . 30 3.1.2 The Lattice Boltzmann Equation . 32 3.1.3 Lattices . 35 3.1.4 Boundary Conditions . 36 viii 3.1.5 Multiple Relaxation Time Collision Operator . 40 3.1.6 Heat Transfer . 41 3.1.7 Particle Concentration Transport . 44 3.2 Turbulence and Transition Simulation . 46 3.2.1 Modeling and Simulation of Turbulence . 46 3.2.2 Direct Numerical Simulation and Large Eddy Simulation . 47 3.2.3 Evaluation Methodology . 49 3.3 Substrate Approximation . 50 3.4 Computational Language and Parallelization . 52 4 The Incompressible LBE ......................... 55 4.1 Introduction . 56 4.2 Derivation of the D2Q9 Incompressible LBE . 59 4.2.1 Standard LBE Form . 59 4.2.2 An Incompressible LBE Derivation . 60 4.2.3 Discussion and Comparison . 67 4.3 Numerical Results . 68 4.3.1 Developing Channel Flow . 69 4.3.2 Lid Driven Cavity . 70 4.3.3 Womersley Flow . 73 4.3.4 Backward Facing Step . 77 4.3.5 3D Lid Driven Cavity . 78 ix 4.4 Conclusions . 80 4.5 MRT-iLBE for the D2Q9 Form . 81 5 2D DNS and 3D LES with the iLBE ................. 85 5.1 Transition Flow with an Incompressible Lattice Boltzmann Equation (2D DNS) . 87 5.1.1 Introduction . 88 5.1.2 Incompressible Multiple Relaxation Time Lattice Boltzmann Equation . 91 5.1.3 Simulation Parameters . 95 5.1.3.1 Post Processing . 97 5.1.4 Results . 100 5.1.4.1 Verifying MRT-iLBM with steady LDC Flow . 100 5.1.4.2 The First Hopf Bifurcation and Transition to Chaos 102 5.1.4.3 Post-Transition Statistics . 105 5.1.5 Conclusions . 112 5.2 3D LES of a Driven Cavity . 114 5.2.1 Simulation Parameters . 115 5.2.2 Vreman Sub-grid Scale Model . 119 5.2.3 Statistical Processing and Post Processing . 121 5.2.4 Numerical Results . 122 5.2.4.1 Verifying D3Q15 MRT-iLBE with Steady Flow . 122 x 5.2.4.2 The Hopf Bifurcation . 125 5.2.5 Conclusions . 129 6 Substrate Generation and Verification ................ 131 6.1 Introduction . 132 6.2 Substrate Generation and Measurement Algorithm and Pseudocode 134 6.3 The Lattice Boltzmann Method . 137 6.3.1 The Incompressible Lattice Boltzmann Equation . 137 6.3.2 The Lattice . 141 6.3.3 Boundary Conditions . 142 6.3.4 Multiphysics Capability . 143 6.4 Verification Methodology . 144 6.4.1 Darcy's Law . 144 6.4.2 Ergun Equation . 145 6.4.3 Koponen Model . 147 6.4.4 A New Semi-Empirical Model . 148 6.4.4.1 The Conduit-Model: Gebart Model . 149 6.4.4.2 Flow Past a Submerged Body Model and the Cell Model . 150 6.4.4.3 Tortuosity Factor . 151 6.4.4.4 Blending functions/Basis-functions . 152 6.4.4.5 The Model Equations . 153 xi 6.5 Results . 155 6.5.1 Verification by Darcy's Law . 157 6.5.2 Verification by Ergun Equation . 158 6.5.3 Verification by Porosity-Permeability Correlations . 159 6.5.3.1 Verification by Koponen Model . 160 6.5.3.2 Verification by Gebart-Brinkmann (GB) Model . 161 6.5.3.3 Verification by Gebart-Kubawara (GK) Model . 161 6.5.4 Experimental Verification . 163 6.6 Conclusions . 164 7 Multiphysics and PM Filtration .................... 167 7.1 Heat Transfer . 168 7.1.1 Recovery of the Advection-Diffusion Equation . 168 7.1.2 Numerical Results . 168 7.1.2.1 Thermal Lid Driven Cavity . 169 7.1.2.2 Conjugate Heat Transfer-Channel Flow . 170 7.2 Particle Transport and Deposition . 173 7.2.1 Diffusion, Advection, and Deposition Bases . 173 7.2.2 Implementation on a Lattice . 174 7.2.3 Numerical Results . 176 7.2.3.1 Backward Facing Step . 176 7.2.3.2 Cylinder in Flow . 179 xii 7.3 Full Coupling in a Filtration Model . 186 7.3.1 Low Porosity and High Sticking Probability . 188 7.3.2 High Porosity and Moderate Sticking Probability . 196 7.3.3 Low Porosity and Low Sticking Probability . 201 8 Conclusions ................................. 211 8.1 Contributions to the Incompressible Lattice Boltzmann Method and its Applications . 211 8.2 Summary of Findings . 212 8.2.1 Incompressible Lattice Boltzmann Equation . 212 8.2.2 Turbulent Transition . 214 8.2.3 Multiphysics Filtration . 215 8.3 Future Work . 219 Bibliography .................................. 225 A iLBE Derivations .............................. 237 A.1 Chapman-Enskog Expansion for the iLBE . 237 A.2 MRT-iLBE . 243 A.2.1 D2Q9 . 243 A.2.2 D3Q15 . 249 B Thermal D2Q4 LBE Derivations .................... 253 B.1 Equilibrium Distribution Form: Ansatz Method . 253 xiii B.2 Chapman-Enskog Expansion for the Thermal D2Q4 LBE . 255 C Units ..................................... 259 D Filter Pseudocode ............................. 265 E Copyright Agreements and Permissions ............... 271 xiv List of Figures 1.1 Schematic of DPF flow. 5 3.1 PDF of He-4 particle velocities at different temperatures. 32 3.2 Some valid lattices of the LBM. 36 3.3 Bounce-back schematics over a single time step. 39 3.4 Speedup of LDC OpenMP code on Superior. 54 4.1 D2Q9 lattice velocities, a. ....................... 61 4.2 Developing channel flow convergence of umax error with node count. 71 4.3 Developing channel flow residual history of standard and incompress- ible LBEs. 72 4.4 Comparison of x and y-velocities along the centerlines in the lid driven cavity flow for Re 1000 with Marchi et al.. 74 4.5 Lid driven cavity flow streamlines at Re 1000. 75 4.6 Womersley flow centerline velocity as a function of time and y. 76 4.7 Backward facing step streamlines for Re 800. 77 4.8 D3Q19 Lattice. ..
Recommended publications
  • Large Eddy Simulation of Turbulent Flow Over a Backward-Facing Step
    IT 18 009 Examensarbete 45 hp Mars 2018 Large Eddy Simulation of Turbulent Flow over a Backward-Facing Step Edmond Shehadi Institutionen för informationsteknologi Department of Information Technology Abstract Large Eddy Simulation of Turbulent Flow over a Backward-Facing Step Edmond Shehadi Teknisk- naturvetenskaplig fakultet UTH-enheten This work studies the effect of grid resolution and subgrid scale modeling on the predictive accuracy of Large Eddy Simulation. In particular, the problems considered Besöksadress: are: turbulent flow over a backward-facing step and fully-developed turbulent channel Ångströmlaboratoriet Lägerhyddsvägen 1 flow. A SGS-free model and four subgrid scale models are used: Smagorsinky, Hus 4, Plan 0 k-equation, dynamic k-equation and WALE. Different combinations of streamwise and spanwise grid resolutions are considered along with different strategies for the Postadress: cell-size distribution in the wall-normal direction. Box 536 751 21 Uppsala First, a detailed study pertaining to fully-developed turbulent channel flow is Telefon: presented for target friction Reynolds numbers 180 and 300. A symmetric 018 – 471 30 03 three-layered wall-normal meshing strategy in conjunction with a SGS-free model is Telefax: shown to be the best compromise between computational efficiency and agreement 018 – 471 30 00 with benchmark data. Overall, the accuracy of the results was observed to be most sensitive to the spanwise resolution of the grid. Cell sizes in wall units ranging in Hemsida: between 25-28 and 12-20 in the streamwise and spanwise directions, respectively, http://www.teknat.uu.se/student yielded excellent agreement with reference data. Second, a detailed study of turbulent flow over a backward-facing step at a step-height Reynolds number of 5100 is conducted.
    [Show full text]
  • Large Eddy Simulation of the Pulsating, Non-Reacting Flow in Combustion Chambers
    Large Eddy Simulation of the Pulsating, Non-Reacting Flow in Combustion Chambers Balázs Pritz Karlsruhe 2009 II Cover page: Iso-surfaces of the Q-criterion at 5·104 s-2 in the model combustion chamber III IV Large Eddy Simulation of the Pulsating, Non-Reacting Flow in Combustion Chambers Zur Erlangung des akademischen Grades Doktor der Ingenieurwissenschaften der Fakultät für Maschinenbau Universität Karlsruhe (TH) genehmigte Dissertation von Dipl.-Ing. Balázs Pritz aus Budapest (Ungarn) Tag der mündlichen Prüfung: 03.08.2009 Hauptreferent: Prof. Dr.-Ing. M. Gabi Korreferenten: Prof. Dr.-Ing. habil. H. Büchner Dr. L. Kullmann V VI Abstract It is well known that in order to fulfil the stringent demands for low emissions of NOx, the lean premixed combustion concept is commonly used. However, lean premixed combustors are susceptible to thermo-acoustic instabilities driven by the combustion process and possibly sustained by a resonant feedback mechanism coupling pressure and heat release. This resonant feedback mechanism creates pulsations typically in the frequency range of several hundred Hz, which reach high amplitudes so that the system has to be shut down or it is even damaged. Although the research activities of the recent years have contributed to a better understanding of this phenomenon, the underlying mechanisms are still not understood well enough. For the prediction of the stability of technical combustion systems regarding the development and maintaining of self-sustained combustion instabilities the knowledge of the periodic-non- stationary mixing and reacting behaviour of the applied flame type and a quantitative description of the resonance characteristics of the gas volumes in the combustion chamber is conclusively needed.
    [Show full text]
  • General Methodology Involved in Solving Problems of Fluid Mechanics
    International Journal of Advances in Engineering Research http://www.ijaer.com/ (IJAER) 2011, Vol. No. 1, Issue No. V, May ISSN: 2231-5152 METHODS INVOLVED IN SOLVING PROBLEMS OF FLUID MECHANICS USING COMPUTATIONAL FLUID DYNAMICS- A STUDY Anuj Kumar Sharma, Research Scholar, Manav Bharti University Solan, HP Parveen Sharma, Asst. Prof.,GNIT Girls Institute of Technology, Greater Noida U.K.Singh, Prof., Kamla Nehru Institute of Technology,Sultanpur ABSTRACT This research paper offers an overview of the techniques used to solve problems in fluid mechanics using computational fluid dynamics (CFD) on computers and describes in detail those most often used in practice. Included are advanced methods in computational fluid dynamics, like direct and large-eddy simulation of turbulence, Probability density function (PDF) methods and vortex methods. It also contains a section dealing with grid quality and an extended description of discretization methods like finite element method (FEM), finite volume method (FVM), finite difference method (FDM) and boundary element method. This paper shows common roots and basic principles for many different methods. Keywords: Fluid Dynamics, Navier–Stokes equations, Hydrodynamics, Boundary layer flow, Computational Fluid Dynamics. INTRODUCTION Computational fluid dynamics (CFD) is a branch of fluid mechanics that uses numerical methods and algorithms to solve and analyze problems that involve fluid flows. Computers are used to perform the calculations required to simulate the interaction of liquids and gases with surfaces defined by boundary conditions. With high-speed supercomputers, better solutions can be achieved. Ongoing research, however, yields software that improves the accuracy and speed of complex simulation scenarios such as transonic or turbulent flows.
    [Show full text]
  • On the Reliability of Large-Eddy Simulation for Dispersed Two-Phase Flows
    On the Reliability of Large-Eddy Simulation for Dispersed Two-Phase Flows Vom Fachbereich Maschinenbau an der Technischen Universität Darmstadt zur Erlangung des Grades eines Doktor-Ingenieurs (Dr.-Ing.) genehmigte Dissertation vorgelegt von Dipl.-Ing. Desislava Nikolova Dimitrova aus Haskovo (Bulgarien) Berichterstatter: Prof. Dr.–Ing. J. Janicka Mitberichterstatter: Prof. Dr.–Ing. P. Stephan Tag der Einreichung: 04. Mai 2010 Tag der mündlichen Prüfung: 07. Juni 2010 Darmstadt 2011 D 17 Abstract "On the Reliability of Large-Eddy Simulation for Dispersed Two-Phase Flows" Desislava Dimitrova This work deals with the numerical description of dilute polydispersed and turbulent two- phase flows. The method used to describe the two-phase flow is based on the EULERIAN- LAGRANGIAN approach. The carrier phase turbulent flow behavior is modeled utilizing Large Eddy Simulation (LES). Most of the recent DNS investigations attempting to provide insight into particle-fluid interac- tion mechanisms are performed on simple, mostly homogeneous turbulent flows with very low REYNOLDS numbers. In contrast to this, LES is an opportunity to extend this studies to flows with significantly higher REYNOLDS numbers and to more complex geometries. Moreover, it is less restrictive to flows containing particles, larger than the KOLMOGOROV length scale. The main objective of this work is to assess the state of the art capabilities of the two-phase LES and to identify weak points. Recommendations for additional model development for increasing the predictive ability are derived. The reliability of transfer of findings from one configuration to others is assessed. The conclusions are based on a systematic variation of relevant flow pa- rameters, such as the REYNOLDS number and the STOKES number, so that a wide range of applications is covered.
    [Show full text]
  • Large Eddy Simulation of Turbulent Flows Using Finite Volume Methods with Structured, Unstructured, and Zonal Embedded Grids Todd Allen Simons Iowa State University
    Iowa State University Capstones, Theses and Retrospective Theses and Dissertations Dissertations 1998 Large eddy simulation of turbulent flows using finite volume methods with structured, unstructured, and zonal embedded grids Todd Allen Simons Iowa State University Follow this and additional works at: https://lib.dr.iastate.edu/rtd Part of the Mechanical Engineering Commons Recommended Citation Simons, Todd Allen, "Large eddy simulation of turbulent flows using finite volume methods with structured, unstructured, and zonal embedded grids " (1998). Retrospective Theses and Dissertations. 11892. https://lib.dr.iastate.edu/rtd/11892 This Dissertation is brought to you for free and open access by the Iowa State University Capstones, Theses and Dissertations at Iowa State University Digital Repository. It has been accepted for inclusion in Retrospective Theses and Dissertations by an authorized administrator of Iowa State University Digital Repository. For more information, please contact [email protected]. INFORMATION TO USERS This manuscript has been reproduced from the microfilm master. UMI films the text directly from the original or copy submitted. Thus, some thesis and dissertation copies are in typewriter face, while others may be from any type of computer printer. The quality of this reproduction is dependent upon the quality of the copy submitted. Broken or indistinct print, colored or poor quality illustrations and photographs, print bleedthrough, substandard margins, and improper alignment can adversely affect reproduction. In the unlikely event that the author did not send UMI a complete manuscript and there are missing pages, these will be noted. Also, if unauthorized copyright material had to be removed, a note will indicate the deletion.
    [Show full text]
  • Large Eddy Simulations of Turbulent and Buoyant Flows in Urban and Complex Terrain Areas Using the Aeolus Model
    atmosphere Article Large Eddy Simulations of Turbulent and Buoyant Flows in Urban and Complex Terrain Areas Using the Aeolus Model Akshay A. Gowardhan *, Dana L. McGuffin , Donald D. Lucas , Stephanie J. Neuscamman, Otto Alvarez and Lee G. Glascoe Lawrence Livermore National Laboratory, Livermore, CA 94551, USA; [email protected] (D.L.M.); [email protected] (D.D.L.); [email protected] (S.J.N.); [email protected] (O.A.); [email protected] (L.G.G.) * Correspondence: [email protected] Abstract: Fast and accurate predictions of the flow and transport of materials in urban and complex terrain areas are challenging because of the heterogeneity of buildings and land features of different shapes and sizes connected by canyons and channels, which results in complex patterns of turbulence that can enhance material concentrations in certain regions. To address this challenge, we have developed an efficient three-dimensional computational fluid dynamics (CFD) code called Aeolus that is based on first principles for predicting transport and dispersion of materials in complex terrain and urban areas. The model can be run in a very efficient Reynolds average Navier–Stokes (RANS) mode or a detailed large eddy simulation (LES) mode. The RANS version of Aeolus was previously validated against field data for tracer gas and radiological dispersal releases. As a part of this work, we have validated the Aeolus model in LES mode against two different sets of data: (1) turbulence quantities measured in complex terrain at Askervein Hill; and (2) wind and tracer data from the Citation: Gowardhan, A.A.; Joint Urban 2003 field campaign for urban topography.
    [Show full text]
  • Arxiv:2009.02411V1 [Physics.Flu-Dyn] 4 Sep 2020 Ρ = Density, Kg/M3 P = Pressure, Pa
    Heat transfer in three-dimensional intersecting shock-wave/turbulent boundary-layer interactions with wall-modeled large-eddy simulations Lin Fu*, Sanjeeb Bose†, Parviz Moin‡ Stanford University, Stanford CA 94305-3024, USA The accurate prediction of aerothermal surface loading is of paramount importance for the design of high speed flight vehicles. In this work, we consider the numerical solution of hypersonic flow over a double-finned geometry, representative of the inlet of an air-breathing flight vehicle, charac- terized by three-dimensional intersecting shock-wave/turbulent boundary-layer interaction at Mach 6 8.3. High Reynolds numbers (ReL ≈ 11:6 × 10 based on free-stream conditions) and the pres- ence of cold walls (Tw=To ≈ 0:27) leading to large near-wall temperature gradients necessitate the use of wall-modeled large-eddy simulation (WMLES) in order to make calculations computation- ally tractable. The comparison of the WMLES results with experimental measurements shows good agreement in the time-averaged surface heat flux and wall pressure distributions, and the WMLES predictions show reduced errors with respect to the experimental measurements than prior RANS calculations. The favorable comparisons are obtained using an LES wall model based on equilib- rium boundary layer approximations despite the presence of numerous non-equilibrium conditions including three dimensionality, shock-boundary layer interactions, and flow separation. Lastly, it is also demonstrated that the use of semi-local eddy viscosity scaling (in lieu of the
    [Show full text]
  • Large-Eddy Simulation: Current Capabilities, Recommended Practices, and Future Research
    NASA/TM—2009-215616 AIAA–2009–948 Large-Eddy Simulation: Current Capabilities, Recommended Practices, and Future Research Nicholas J. Georgiadis Glenn Research Center, Cleveland, Ohio Donald P. Rizzetta Air Force Research Laboratory, Wright-Patterson Air Force Base, Dayton, Ohio Christer Fureby FOI, Swedish Defense Research Agency, Stockholm, Sweden July 2009 NASA STI Program . in Profile Since its founding, NASA has been dedicated to the • CONFERENCE PUBLICATION. Collected advancement of aeronautics and space science. The papers from scientific and technical NASA Scientific and Technical Information (STI) conferences, symposia, seminars, or other program plays a key part in helping NASA maintain meetings sponsored or cosponsored by NASA. this important role. • SPECIAL PUBLICATION. Scientific, The NASA STI Program operates under the auspices technical, or historical information from of the Agency Chief Information Officer. It collects, NASA programs, projects, and missions, often organizes, provides for archiving, and disseminates concerned with subjects having substantial NASA’s STI. The NASA STI program provides access public interest. to the NASA Aeronautics and Space Database and its public interface, the NASA Technical Reports • TECHNICAL TRANSLATION. English- Server, thus providing one of the largest collections language translations of foreign scientific and of aeronautical and space science STI in the world. technical material pertinent to NASA’s mission. Results are published in both non-NASA channels and by NASA in the NASA STI Report Series, which Specialized services also include creating custom includes the following report types: thesauri, building customized databases, organizing and publishing research results. TECHNICAL PUBLICATION. Reports of completed research or a major significant phase For more information about the NASA STI of research that present the results of NASA program, see the following: programs and include extensive data or theoretical analysis.
    [Show full text]
  • A Definition for Large Eddy Simulation Approximations of the Navier-Stokes Equations
    A DEFINITION FOR LARGE EDDY SIMULATION APPROXIMATIONS OF THE NAVIER-STOKES EQUATIONS J.-L. GUERMOND1,3 AND S. PRUDHOMME2 Abstract. This paper proposes a mathematical definition of Large Eddy Sim- ulation (LES) for three-dimensional turbulent incompressible viscous flows. We propose to call LES approximation of the Navier–Stokes equations any se- quence of finite-dimensional approximations of the velocity and pressure fields that converge (possibly up to subsequences) in the Leray class to a suitable weak solution. 1. Introduction 1.1. Introductory comments. At the present time, computer predictions of tur- bulence phenomena by so-called Direct Numerical Simulation (DNS) of the Navier- Stokes equations are still considered a formidable task for Reynolds numbers larger than a few thousands. Since the times of Boussinesq and Reynolds, numerous turbulence models based on time-averaged or space-averaged quantities (Reynolds Averaged Navier–Stokes models, k- models, etc.) have been developed and then used in engineering applications as a means of overcoming, though often with lim- ited success, the lack of sufficient computer resources required by DNS. During the past forty years a new class of turbulence models, collectively known as Large Eddy Simulation (LES) models [23], has emerged in the literature. These models are founded on the observation that representing the whole range of flow scales may not be important in many engineering applications as one is generally inter- ested only in the large scale features of the flow. With this objective in mind, LES modelers have devised artifacts for representing the interactions between the unreachable small scales and the large ones.
    [Show full text]
  • Large Eddy Simulation Applications
    VKI Lecture Series Numerical Investigations in Turbomachinery: A State of the Art Large Eddy Simulation Applications G. Dufour, N. Gourdain, F. Duchaine, O. Vermorel, L.Y.M Gicquel, J.-F. Boussuge and T. Poinsot Notes prepared for the von Karman Institute for Fluid Dynamics September, 21-25 2009 Contact: [email protected] http://www.cerfacs.fr Copyright 2009 by CERFACS Contents 1 Introduction 5 2 Literature overview 7 2.1 Idealized configurations . .8 2.2 Basic configurations . 11 2.3 Complex geometries . 14 2.4 Synthesis . 17 3 Case studies 19 3.1 High-pressure turbine stator . 19 3.2 Conjugate heat transfer for a cooled turbine blade . 22 4 Conclusion 29 Bibliography 31 3 4 CONTENTS Chapter 1 Introduction Since Computational Fluid Dynamics (CFD) simulations have become central in the design process of turbomachinery components, the need for as accurate predictions as possible has followed the constant increase in computing power (see Horlock and Denton [57] for an historical review of design practices). Early on and since the pioneering approaches relying on single element/single passage (see Refs. [51, 4] for instance), CFD simulations have evolved to include multiple rows, using the mixing plane approach under the steady flow assumption [30, 28]. With increasing stage loading and reduced weight, strong interactions between rows have fostered the emergence of fully unsteady simulations [100, 47, 53]. Now, state-of-the-art unsteady simulations can be performed on the entire geometry of a multi-stage compressor, including all the passages of all the blades [141, 44]. All the published literature for such applications presents approaches within the (Unsteady) Reynolds-Averaged Navier-Stokes (U)RANS equations framework.
    [Show full text]
  • Large-Eddy Simulation of Atmospheric Flow Over Complex Terrain
    RIS0 Large-Eddy Simulation of Atmospheric Flow over Complex Terrain Andreas Bechmann Ris0-PhD-28(EN) Ris0 National Laboratory Technical University of Denmark Roskilde, Denmark November 2006 Author: Andreas Bechmann Rise PhD 28(EN) Title: Large-Eddy Simulation of Atmospheric Flow over Complex November 2006 Terrain Department: Wind Energy Department Tins thesis is submitted in partial fulfilment of the requirements for the Ph D. degree at the Technical University of Demnark (DTU) Abstract (max. 2000 char.): ISBN 87-550-3556-6 The present report describes the development and validation of a turbulence model designed for atmospheric flows based on the concept of Large-Eddy Simulation (EES). The background for the work is the high Reynolds number k - s model, which has been implemented on a finite-volume code of the incompressible Reynolds-averaged Navier-Stokes equations (RANS). The k - s model is traditionally used for Contract no.: RANS computations, but is here developed to also enable EES. EES is able to provide detailed descriptions of a wide range of engineering flows at low Reynolds numbers. For Group's own reg. no. : atmospheric flows, however, the high Reynolds numbers and the rough surface of the earth provide difficulties normally not compatible with EES. Since these issues are most severe Sponsorship: near the surface they are addressed by handling the near surface region with RANS and only use EES above this region. Using this method, the developed turbulence model is Cover : able to handle both engineering and atmospheric flows and can be run in both RANS or EES mode. For EES simulations a time-dependent wind field that accurately represents the turbulent structures of a wind environment must be prescribed at the computational inlet.
    [Show full text]
  • Direct Numerical Simulation and Wall-Resolved Large Eddy Simulation in Nuclear Thermal Hydraulics Iztok Tiselj, Cédric Flageul, Jure Oder
    Direct Numerical Simulation and Wall-Resolved Large Eddy Simulation in Nuclear Thermal Hydraulics Iztok Tiselj, Cédric Flageul, Jure Oder To cite this version: Iztok Tiselj, Cédric Flageul, Jure Oder. Direct Numerical Simulation and Wall-Resolved Large Eddy Simulation in Nuclear Thermal Hydraulics. Nuclear Technology, American Nuclear Society, 2019, pp.1-15. 10.1080/00295450.2019.1614381. hal-02154872 HAL Id: hal-02154872 https://hal.archives-ouvertes.fr/hal-02154872 Submitted on 5 Jul 2019 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. NUCLEAR TECHNOLOGY © 2019 American Nuclear Society DOI: https://doi.org/10.1080/00295450.2019.1614381 Critical Review Direct Numerical Simulation and Wall-Resolved Large Eddy Simulation in Nuclear Thermal Hydraulics Iztok Tiselj, * Cedric Flageul, and Jure Oder Jožef Stefan Institute, Jamova 39, 1000 Ljubljana, Slovenia Received March 15, 2019 Accepted for Publication April 30, 2019 Abstract— The critical review discusses the most accurate methods for description of turbulent flows: the computationally very expensive direct numerical simulation (DNS) and slightly less accurate and slightly less expensive large eddy simulation (LES) methods. Both methods have found their way into nuclear thermal hydraulics as tools for studies of the fundamental mechanisms of turbulence and turbulent heat transfer.
    [Show full text]