Danny Messig

Flame structure analysis and flamelet progress variable modelling of strained coal flames

Strained two-phase pulverised coal flames in a counterflow configuration are investigated numerically. Three operating conditions with different coal-to-primary-air ratios and inlet velocities were evaluated in order to establish different flame regimes. At first, the two-phase flow of the fully resolved reference cases is calculated solving the transport equation for the species and directly evaluating the reaction rates. Different flame structures are identified using the heat release rate and the chemical explosive mode as markers, showing that complex structures with a combination of lean premixed and non-premixed flames can be observed in strained counterflow coal flames. In addition to the fully resolved simulation, the suitability of the Flamelet-Progress Variable (FPV) model is investigated. Both premixed and non-premixed tables are employed. At first, the suitability of the look-up tables is evaluated by means of an a priori analysis, using the fully resolved simulations as reference solutions, showing that the non-premixed flamelet table correctly predicts the structure of the strained coal flames, while the premixed table shows sensible deviations in terms of temperature and species, especially at rich conditions. Finally, the a posteriori analysis shows that the fully coupled FPV model with a non-premixed flamelet look-up table can accurately predict strained coal flames.

An abstraction layer for efficient memory management of tabulated chemistry and flamelet solutions

A large number of methods for simulating reactive flows exist, some of them, for example, directly use detailed chemical kinetics or use precomputed and tabulated flame solutions. Both approaches couple the research fields computational fluid dynamics and chemistry tightly together using either an online or offline approach to solve the chemistry domain. The offline approach usually involves a method of generating databases or so-called Lookup-Tables (LUTs). As these LUTs are extended to not only contain material properties but interactions between chemistry and turbulent flow, the number of parameters and thus dimensions increases. Given a reasonable discretisation, file sizes can increase drastically. The main goal of this work is to provide methods that handle large database files efficiently. A Memory Abstraction Layer (MAL) has been developed that handles requested LUT entries efficiently by splitting the database file into several smaller blocks. It keeps the total memory usage at a minimum using thin allocation methods and compression to minimise filesystem operations. The MAL has been evaluated using three different test cases. The first rather generic one is a sequential reading operation on an LUT to evaluate the runtime behaviour as well as the memory consumption of the MAL. The second test case is a simulation of a non-premixed turbulent flame, the so-called HM1 flame, which is a well-known test case in the turbulent combustion community. The third test case is a simulation of a non-premixed laminar flame as described by McEnally in 1996 and Bennett in 2000. Using the previously developed solver ‘flameletFoam’ in conjunction with the MAL, memory consumption and the performance penalty introduced were studied. The total memory used while running a parallel simulation was reduced significantly while the CPU time overhead associated with the MAL remained low.

Application of Algorithmic Differentiation for Exact Jacobians to the Universal Laminar Flame Solver

We introduce algorithmic differentiation (AD) to the C++ Universal Laminar Flame (ULF) solver code. ULF is used for solving generic laminar flame configurations in the field of combustion engineering. We describe in detail the required code changes based on the operator overloading-based AD tool CoDiPack. In particular, we introduce a global alias for the scalar type in ULF and generic data structure using templates. To interface with external solvers, template-based functions which handle data conversion and type casts through specialization for the AD type are introduced. The differentiated ULF code is numerically verified and performance is measured by solving two canonical models in the field of chemically reacting flows, a homogeneous reactor and a freely propagating flame. The models stiff set of equations is solved with Newtons method. The required Jacobians, calculated with AD, are compared with the existing finite differences (FD) implementation. We observe improvements of AD over FD. The resulting code is more modular, can easily be adapted to new chemistry and transport models, and enables future sensitivity studies for arbitrary model parameters.