Ein Gas-Kinetic Scheme Ansatz zur Modellierung und Simulation von Feuer auf massiv paralleler Hardware

This work presents a simulation approach based on a Gas Kinetic Scheme (GKS) for the simulation of fire that is implemented on massively parallel hardware in terms of Graphics Processing Units (GPU) in the framework of General Purpose computing on Graphics Processing Units (GPGPU). Gas kinetic schemes belong to the class of kinetic methods because their governing equation is the mesoscopic Boltzmann equation, rather than the macroscopic Navier-Stokes equations. Formally, kinetic methods have the advantage of a linear advection term which simplifies discretization. GKS inherently contains the full energy equation which is required for compressible flows. GKS provides a flux formulation derived from kinetic theory and is usually implemented as a finite volume method on cell-centered grids. In this work, we consider an implementation on nested Cartesian grids. To that end, a coupling algorithm for uniform grids with varying resolution was developed and is presented in this work. The limitation to local uniform Cartesian grids allows an efficient implementation on GPUs, which belong to the class of many core processors, i.e. massively parallel hardware. Multi-GPU support is also implemented and efficiency is enhanced by communication hiding. The fluid solver is validated for several two- and three-dimensional test cases including natural convection, turbulent natural convection and turbulent decay. It is subsequently applied to a study of boundary layer stability of natural convection in a cavity with differentially heated walls and large temperature differences. The fluid solver is further augmented by a simple combustion model for non-premixed flames. It is validated by comparison to experimental data for two different fire plumes. The results are further compared to the industry standard for fire simulation, i.e. the Fire Dynamics Simulator (FDS). While the accuracy of GKS appears slightly reduced as compared to FDS, a substantial speedup in terms of time to solution is found. Finally, GKS is applied to the simulation of a compartment fire. This work shows that the GKS has a large potential for efficient high performance fire simulations.Diese Arbeit präsentiert einen Simulationsansatz basierend auf einer gaskinetischen Methode (eng. Gas Kinetic Scheme, GKS) zur Simulation von Bränden, welcher für massiv parallel Hardware im Sinne von Grafikprozessoren (eng. Graphics Processing Units, GPUs) implementiert wurde. GKS gehört zur Klasse der kinetischen Methoden, die nicht die makroskopischen Navier-Stokes Gleichungen, sondern die mesoskopische Boltzmann Gleichung lösen. Formal haben kinetische Methoden den Vorteil, dass der Advektionsterms linear ist. Dies vereinfacht die Diskretisierung. In GKS ist die vollständige Energiegleichung, die zur Lösung kompressibler Strömungen benötigt wird, enthalten. GKS formuliert den Fluss von Erhaltungsgrößen basierend auf der gaskinetischen Theorie und wird meistens im Rahmen der Finiten Volumen Methode umgesetzt. In dieser Arbeit betrachten wir eine Implementierung auf gleichmäßigen Kartesischen Gittern. Dazu wurde ein Kopplungsalgorithmus für die Kombination von Gittern unterschiedlicher Auflösung entwickelt. Die Einschränkung auf lokal gleichmäßige Gitter erlaubt eine effiziente Implementierung auf GPUs, welche zur Klasse der massiv parallelen Hardware gehören. Des Weiteren umfasst die Implementierung eine Unterstützung für Multi-GPU mit versteckter Kommunikation. Der Strömungslöser ist für zwei und dreidimensionale Testfälle validiert. Dabei reichen die Tests von natürlicher Konvektion über turbulente Konvektion bis hin zu turbulentem Zerfall. Anschließend wird der Löser genutzt um die Grenzschichtstabilität in natürlicher Konvektion bei großen Temperaturunterschieden zu untersuchen. Darüber hinaus umfasst der Löser ein einfaches Verbrennungsmodell für Diffusionsflammen. Dieses wird durch Vergleich mit experimentellen Feuern validiert. Außerdem werden die Ergebnisse mit dem gängigen Brandsimulationsprogramm FDS (eng. Fire Dynamics Simulator) verglichen. Die Qualität der Ergebnisse ist dabei vergleichbar, allerdings ist der in dieser Arbeit entwickelte Löser deutlich schneller. Anschließend wird das GKS noch für die Simulation eines Raumbrandes angewendet. Diese Arbeit zeigt, dass GKS ein großes Potential für die Hochleistungssimulation von Feuer hat

Compressible Shear Flow Transition and Turbulence: Enhancement of GKM Numerical Scheme and Simulation/Analysis of Pressure Effects on Flow Stabilization

Despite significant advancements in the understanding of fluid flows, combustion and material technologies, hypersonic flight still presents numerous technological challenges. In hypersonic vehicles turbulence is critical in controlling heat generation in the boundary layer, mixing inside the combustor, generation of acoustic noise, and mass flow in the intake. The study of turbulence in highly compressible flows is challenging compared to incompressible due to a drastic change in the behavior of pressure and a relaxation of the incompressibility constraint. In addition fluid flow inside a flight vehicle is complicated by wall-effects, heat generation and complex boundary conditions. Homogeneous shear flow contains most of the relevant physics of boundary and mixing layers without the aforementioned complicating effects. In this work we aim to understand and characterize the role of pressure, velocity-pressure interaction, velocity-thermodynamics interaction in the late-stage transition-to-turbulence regime in a high speed shear dominated flow by studying the evolution of perturbations in in a high Mach number homogeneous shear flow. We use a modal-analysis based approach towards understanding the statistical behavior of turbulence. Individual Fourier waves constituting the initial flow field are studied in isolation and in combination to understand collective statistical behavior. We demonstrate proof of concept of novel acoustic based strategies for controlling the onset of turbulence. Towards this goal we perform direct numerical simulations (DNS) in three studies: (a) development and evaluation of gas kinetic based numerical tool for DNS of compressible turbulence, and perform detailed evaluation of the efficacy of different interpolation schemes in capturing solenoidal and dilatational quantities, (b) modal investigation in the behavior of pressure and isolation of linear, non-linear, inertial and pressure actions, and (c) modal investigation in the possible acoustic based control strategies in homogeneously sheared compressible flows. The findings help to understand the manifestation of the effects of compressibility on transition and turbulence via the velocity-pressure interactions and the action of individual waves. The present study helps towards the design of control mechanisms for compressible turbulence and the development of physically consistent pressure strain correlation models

Gas-kinetic scheme for rarefied flow simulation

For increasingly rarefied flowfields, the predictions from continuum formulation, such as the Navier–Stokes equations lose accuracy. These inaccuracies are attributed primarily to the linear approximations of the stress and heat flux terms in the Navier– Stokes equations. The inclusion of higher order terms, such as Burnett, or high-order moment equations, could improve the predictive capabilities of such continuum formulations, but there has been limited success in the shock structure calculations, especially in the high Mach number case. Here, after reformulating the viscosity and heat conduction coefficients appropriate for the rarefied flow regime, we will show that the extended Navier–Stokes-type continuum formulation may still be properly used. The equations with generalized dissipative coefficients based on the closed solution of the Bhatnagar–Gross–Krook (BGK) model of the Boltzmann equation, are solved using the gas-kinetic numerical scheme. © 2006 IMACS. Published by Elsevier B.V. All rights reserved