Highly-resolved numerical simulation of turbulent premixed and stratified combustion under adiabatic and non-adiabatic conditions with tabulated chemistry
In this thesis, premixed and stratified turbulent combustion processes under adiabatic and non-adiabatic conditions are investigated with highly resolved large eddy simulation (LES) and direct numerical simulation (DNS). The flamelet generated manifolds (FGM) approach is used to model the combustion chemistry. The simulations are performed with the ‘PsiPhi’ CFD solver (originally developed by Kempf), which is improved, extended, accelerated and made suitable for efficient computations on ten thousands of processor cores.
The FGM model is extended to account for the effect of stratified mixture and heat losses. The models are validated by LES with the artificial thickened flame model (ATF) of a bluff-body burner investigated at Cambridge/Sandia and a model gas turbine combustor investigated at the DLR Stuttgart, respectively. Radial statistics of the velocity components, temperature and major species mass fractions obtained from simulation and measurement are compared and a good agreement is found. The combustion process is investigated closer with probability density functions (PDFs) and joint probability density functions (JPFs). For the Cambridge/Sandia burner, it is shown that the flame burns in a premixed regime near the burner exit and in a stratified regime further downstream and that the combustion process can be well described by an ensemble of premixed flamelets burning at the local equivalence ratio. For the DLR combustor, it is shown that the flame is exposed to massive heat loss throughout the whole combustor, a newly introduced model for the inclusion of heat loss in the one-dimensional FGM model flames yields comparable results to established models while allowing for a large range of heat losses in the FGM table generation. A further introduced model for a correction factor depending on the local strain rate improves the temperature and species predictions further downstream.
Afterwards, the Cambridge/Sandia burner is simulated with a grid resolution of 100 μm, which is fine enough to resolve all relevant scales directly. The computational domain consists of 1.6 billion cells, requiring 10 million core hours on 64,000 processor cores. The turbulent flame structure is analyzed with PDFs and JPFs, where a conditioning on the distance from the burner and the progress of the reaction is necessary, as the location of the local flame front in the combustion regime diagram depends on these quantities. The analysis is carried out with respect to previous findings from DNS of simpler flame configurations, it is demonstrated that most of the findings are comparable. The turbulent flame thickness is found to be mostly slightly smaller than the laminar one due to predominantly compressive flame normal strain rate. The displacement speed of the turbulent flame front is shown to be mainly determined by the curvature of the progress variable field on the unburned side of the flame and the progress variable gradient on the burned side, partially negative values also occur. The results are filtered at different filter widths and used for a-priori and, together with LES results with corresponding cell sizes, a-posteriori investigation of sub-filter closures for the ATF model. The established sub-filter model by Charlette is used with a static and a dynamic model parameter, together with a newly introduced dynamic power-law model. The static model over predicts the sub-filter wrinkling, whereas both dynamic models are able to predict the sub-filter wrinkling well except for a certain under prediction for the largest filter width. It is found that this under prediction is compensated in the LES by larger local gradients of the progress variable field. The LES results converge against the DNS results with grid refinement for the dynamic models.
In dieser Arbeit werden vorgemischte und partiell vorgemischte Verbrennungsprozesse unter adiabaten und nicht-adiabaten Bedingungen mittels hochaufgelöster Large Eddy Simulation(LES) und direkter numerischer Simulation (DNS) untersucht. Die Verbrennungs-Chemie wird mit dem ‘flamelet generated manifolds’ (FGM) Ansatzes modelliert. Die Simulationen werden mit dem ‘PsiPhi’ CFD Löser (ursprüglich von Kempf entwickelten) durchgeführt, der verbessert, erweitert, beschleunigt und für effiziente Berechnungen auf zehntausenden Prozessorkernen anwendbar gemacht wird.
Das FGM Modell wird für partiell vorgemischte und nicht-adiabate Bedingungen erweitert. Die Modelle werden durch LES eines in Cambridge/Sandia untersuchten Staukörper-Brenners sowie eines bei der DLR Stuttgart untersuchten Gasturbinen-Brennkammer-Modells validiert, dabei wird das ‘artificial thickened flame’ (ATF) Modell verwendet. Radiale Statistiken von Geschwindigkeits-Komponenten, Temperatur und Massenbrüche der Haupt-Spezies aus Simulation und Experiment zeigen eine gute Ubereinstimmung. Der Verbrennungs-Prozess wird mittels Wahrscheinlichkeitsdichtefunktionen (PDFs) und Verbund-PDFs (JPDFs) genauer untersucht. Es wird gezeigt, dass die Flamme beim Cambridge/Sandia Brenner nahe des Brenneraustritts vorgemischt brennt und weiter stromabwärts partiell vorgemischt, was sich durch ein Ensemble kleiner vorgemischter Flammen mit dem lokalen Verbrennungsluftverhältnis gut beschreiben lässt. Beim DLR Brenner ist die Flamme starken Wärmeverlusten ausgesetzt, wobei ein neu eingeführtes Modell f ̈ur die Einbringung von Wärmeverlusten in die eindimensionalen FGM-Modellflammen vergleichbare Ergebnisse mit etablierten Modellen ergibt und dabei einen größeren maximalen Wärmeverlust bei der Tabellenerzeugung erlaubt. Ein ebenfalls neu eingeführter Korrekturfaktor abhängig von der lokalen Verformungsgeschwindigkeit führt zu einer Verbesserung der Temperatur- und Spezies-Vorhersagen weiter stromabwärts.
Im Anschluß wird der Cambridge/Sandia Brenner mit einer Netzauflösung von 100 μm simuliert, welche fein genug ist um alle relevanten Skalen direkt aufzulösen. Das Rechengebiet setzt sich aus 1,6 Milliarden Zellen zusammen, es werden 10 Millionen Rechenstunden auf 64.000 Prozessorkernen benötigt. Die Struktur der turbulenten Flamme wird mittels PDFs und JPDFs analysiert, wobei eine Konditionierung auf den Abstand vom Brenneraustritt sowie den Reaktionsfortschritt erforderlich ist, da die Positionierung der lokalen Flammenfront im Verbrennungs-Regime-Diagramm von diesen Faktoren abhängt. Die Resultate werden mit früheren DNS von einfacheren Flammen-Konfiguration verglichen, es wird gezeigt dass viele der Erkenntnisse vergleichbar sind. Die Dicke der turbulenten Flamme ist etwas geringer als die der laminaren Flamme, der Grund dafür ist eine überwiegende Druckverformung in Flammen-Normalenrichtung. Die Ausbreitungsgeschwindigkeit der turbulenten Flammenfront wird auf der unverbrannten Seite der Flamme durch die Hauptkrümmung des Fortschrittsvariablenfelds bestimmt, auf der verbrannten Seite hingegen durch den Gradienten des Fortschrittsvariablenfelds, dabei treten lokal auch negative Werte auf. Die Ergebnisse werden mit unterschiedlich
großen Filtern gefiltert und für a-priori sowie, in Verbindung mit entsprechenden LES Ergebnissen, a-posteriori Untersuchungen von Schließungsansätzen für das ATF-Modell verwendet. Dabei wird das etablierte Modell von Charlette et al. mit einem statischen und einem dynamischen Modellparameter sowie ein neu eingeführtes dynamische Potenzgesetz verwendet. Das statische Model überschätzt dabei die nicht-aufgelöste Flammenkrümmung, wohingegen beide dynamischen Modelle diese passend vorhersagen, mit Ausnahme einer leicht zu niedrigen Vorhersage für die grösste Filterweite. Es zeigt sich, dass diese zu niedrige Vorhersage in der LES durch größere lokale Gradienten des Fortschrittsvariablenfelds kompensiert wird. Die LES Ergebnisse nähern sich dabei mit feiner werdendem Netz an die DNS Ergebnisse an.
Vorschau
Zitieren
Zitierform:
Zitierform konnte nicht geladen werden.