Conservative LES-CMC modelling for turbulent jet flames

Abstract

The objective of the present study is to analyze turbulent non-premixed flames by utilizing a combined large eddy simulation - conditional moment closure (LES-CMC) method. LES is based on spatial filtering, and it resolves large scales of turbulent motion while modelling the small turbulent structures using a subgrid model, here the Smagorinsky model. CMC is a conserved scalar method where fluctuations of the reactive scalar variables can be associated with fluctuations of the mixture fraction. Therefore, CMC is applied to turbulent combustion modelling in this work using mixture fraction as the conditioning variable. In the last decade, computations using a non-conservative LES-CMC formulation have provided good predictions of major and minor species for different flames. However, inaccurate predictions occur in CMC cells which have large temporal variations of the mixture fraction field. A lack of weighting the convective term by a filtered density function (FDF) ratio in non-conservative CMC is believed to be a major reason for these inaccurate predictions. In contrast to non-conservative LES-CMC, the present conservative formulation is inherently mass conserving. It considers weighting the convective term by an FDF ratio so that improved predictions of local conditional scalars can be obtained.

Investigations of turbulent jet flames (Sandia Flames D, E and F) are performed using the conservative LES-CMC approach. Flame D is used as the first test case to validate the numerical results in comparison with well-established experimental data. A study of the flow and mixing parameters is carried out first to establish the parameters for Flames E and F. Results from these studies show that the optimal values of Schmidt number, Sc, and turbulent Schmidt number, Sct are 0.7 and 0.4, respectively. The appropriate value of the subgrid-scale variance modelling constant is 0.2. A sensitivity analysis of the results demonstrates that inflow velocity variance levels corresponding to 2/3 , 1/3 and 2/9 of the measured variances at z/D = 0.14 are suitable inflow conditions for Flames D, E and F, respectively.

Subsequently, parametric studies of the combustion model are performed for all test cases. The statistical predictions of scalars compared with measurements show that the LES with the conservative CMC formulation is better than the one based on non-conservative equations. However, similar predictions are obtained from two different flux approximation methods (computing the CMC convective fluxes based on the LES cells located at the CMC faces or based on the CMC cell centres). A parametric study of the CMC grid resolution shows that a resolution of 8*8*80 cells in x-, y- and z-directions yields appropriate predictions within a reasonable computational time. Simulation results of Flames E and F show that the CMC simulations presented here cannot capture local extinction and re-ignition accurately. This is partially due to the averaging effect of the large CMC cells on the modelled conditional dissipation. Much finer CMC cells of the order of the LES cell size will capture more of the fluctuation of scalar dissipation rates and may lead to a more accurate prediction of the local extinction events. The parametric study of CMC grid resolution for Flames E and F shows that a finer CMC grid (16*16*80) predicts slightly better results than the reference grid (8*8*80), but predictions could still be improved. It is possible that the problem is associated with the accuracy of the chemical source term. Hence, some possible solutions, such as second-order closure and doubly conditional reaction source terms, are suggested for future works.

Das Ziel der vorliegenden Arbeit ist die Analyse turbulenter, nicht-vorgemischter Flammen mittels einer kombinierten Large Eddy Simulation - Conditional Moment Closure (LES-CMC) Methode. LES beruht auf räumlicher Filterung und löst die großen turbulenten Skalen auf, wohingegen die kleinen turbulenten Strukturen durch ein Feinstrukturmodell beschrieben werden, in diesem Fall durch das Smagorinsky-Modell. CMC basiert auf einer skalaren Erhaltungsgröße, dem Mischungsbruch, und es wird angenommen, dass die Fluktuationen der reaktiven Skalare mit den Fluktuationen des Mischungsbruchs korrelieren. Deshalb wird CMC mit Konditionierung auf den Mischungsbruch in dieser Arbeit eingesetzt, um turbulente Flammen vorherzusagen. Im letzten Jahrzehnt hat eine Reihe an Simulationen basierend auf der nicht-konservativen LES-CMC Formulierung zu guten Vorhersagen der chemischen Spezies in unterschiedlichen Flammen geführt. Ungenaue Vorhersagen ergeben sich allerdings in CMC-Zellen, die starke zeitliche Schwankungen des Mischungsbruchsfelds aufweisen. Ein wichtiger Grund für diese schwachen Vorhersagen liegt im Fehlen einer Gewichtung des konvektiven Terms mit dem Verhältnis der filtered density functions (FDFs) in den nicht-konservativen CMC-Gleichungen. Im Gegensatz zur nicht-konservativen LES-CMC erhält die hier vorgestellte konservative Formulierung streng die Masse. Die Gewichtung des konvektiven Terms mit dem FDF-Verhältnis wird berücksichtigt, so dass verbesserte Vorhersagen der lokalen, konditionierten Skalare erzielt werden können.

Simulationen turbulenter Strahlflammen (Sandia Flammen D, E und F) mittels der konservativen LES-CMC Methode werden durchgeführt. Flamme D wird als ein erster Testfall zur Validierung der numerischen Ergebnisse anhand etablierter experimenteller Daten verwendet. Eine Studie des Strömungs- und Mischungsfeldes der Flamme D dient zur Festlegung der relevanten Parameter für die Flammen E und F. Die Studie ergibt optimale Werte für die laminare und turbulente Schmidt-Zahl von 0,7 und 0,4. Ein geeigneter Wert für die Konstante im Modell für die Feinstruktur-Varianz ergibt sich zu 0,2. Eine Sensitivitätsstudie des Varianzlevels am Einstromrand zeigt, dass Varianzwerte von 2/3 , 1/3 und 2/9 der gemessenen Geschwindigkeitsvarianz bei z/D = 0,14 geeignete Einstrom-Randbedingungen für die Flammen D, E und F darstellen.

Nachfolgend wird eine Studie der Parameter des Verbrennungsmodells für alle Testfälle durchgeführt. Der Vergleich der Skalarvorhersagen mit Messungen zeigt, dass LES kombiniert mit der konservativen CMC-Formulierung bessere Ergebnisse erzielt, als auf Basis der nicht-konservativen Gleichungen. Allerdings ergeben sich ähnliche Vorhersagen beim Vergleich zweier Approximationsmodelle für den konvektiven Fluss in CMC (Berechnung der konvektiven CMC-Flüsse basierend auf den LES-Zellen am CMC-Zellrand oder basierend auf dem Wert am CMC-Zellmittelpunkt). Eine Parameterstudie der CMC-Gitterauflösung zeigt, dass 8*8*80 Zellen in x-, y- und z-Richtung akzeptable Vorhersagen innerhalb einer vertretbaren Rechenzeit ergeben. Die Simulationsergebnisse der Flammen E und F zeigen, dass die vorliegenden CMC-Simulationen lokales Verlöschen und Wiederentzünden nicht ausreichend abbilden. Dies liegt zum Teil am Mittelungseffekt der großen CMC-Zellen auf die konditionierte Dissipation. Wesentlich kleinere CMC-Zellen in der Größenordnung der LES-Filterweite könnten Extremwerte hoher Dissipation deutlich besser abbilden und somit zu genaueren Vorhersagen des lokalen Verlöschens führen. Eine Parameterstudie der CMC-Gitterauflösung für die Flammen E und F ergibt etwas bessere Ergebnisse auf einem feineren CMC-Gitter (16*16*80) im Vergleich zum Referenzgitter (8*8*80), wobei diese Vorhersagen immer noch verbesserungswürdig sind. Das Problem liegt möglicherweise in der Genauigkeit der Beschreibung des chemischen Quellterms. Aus diesem Grund werden mögliche Vorschläge für zukünftige Arbeiten gemacht, wie eine Schließung zweiter Ordnung und zweifach-konditionierte chemische Quellterme.