Modellierung der partiellen Oxidation von Methan

Motivation

Die Verbrennung eines Brenngases bestehend aus reinen Kohlenwasserstoffen und dem Oxidationsmittel Sauerstoff erlaubt über die Betrachtung der Stöchiometrie eine Einteilung in zwei Bereiche: Die Verbrennung von Brenngasen unterstöchiometrischer oder überstöchiometrischer Zusammensetzung. Den Grenzfall bildet dabei die Verbrennung eines Brenngases in stöchiometrischer Zusammensetzung, wobei, im Fall von reinen Kohlenwasserstoffen, die Sauerstoffmenge im Brenngas exakt ausreicht, sodass bei vollständigem Umsatz nur die Produkte Wasser und Kohlenstoffdioxid entstehen. Bei der Verbrennung eines Brenngases überstöchiometrischer Zusammensetzung ist der Sauerstoffanteil größer als im Fall stöchiometrischer Zusammensetzung und bei der Verbrennung eines Brenngases unterstöchiometrischer Zusammensetzung entsprechend geringer. Ein Charakteristikum dieser Einteilung besteht in den Produkten der Verbrennung. Während bei der Verbrennung von Brenngasen mit überstöchiometrischer Zusammensetzung der Sauerstoffüberschuss einen vergleichsweise geringen Einfluss auf die Produktzusammensetzung zeigt, treten bei der Verbrennung von Brenngasen mit unterstöchiometrischer Zusammensetzung, neben Wasser und Kohlenstoffdioxid, nennenswerte Mengen weiterer Produkte wie Wasserstoff, Kohlenstoffmonoxid, höhere Kohlenwasserstoffe bis hin zu Festoffen wie Ruß auf.

Bei “klassischer“ Verwendung der Verbrennung zur Bereitstellung hoher Temperaturen, beispielsweise bei der Energieumwandlung, werden Brenngase in überstöchiometrischer Zusammensetzung verbrannt, da die Produkte der Verbrennung von Brenngasen unterstöchiometrischer Zusammensetzung, beispielsweise, als Schadstoffe für die Umwelt deklariert sind.
In der chemischen Industrie werden die Produkte der Verbrennung eines Brenngases mit unterstöchiometrischer Zusammensetzung als Ausgangsstoffe für die Produktion von Düngemitteln, Kunststoffen, Elektroden, Schmierstoffen, Spezialchemikalien und vielem mehr eingesetzt. Die genaue Spezieszusammensetzung der Abgase aus der Verbrennung eines Brenngases unterstöchiometrischer Zusammensetzung in Abhängigkeit der technischen Randbedingungen ist Gegenstand dieser Forschung am Institut für Verbrennungstechnik. Die numerische Simulation soll dabei als Werkzeug zur Abbildung von technischen Systemen zum Einsatz kommen, da solche Systeme nur unter hohem Aufwand experimentell untersucht werden können und in vielen Fällen keine Messtechnik zur Untersuchung unter den vorliegenden Bedingungen existiert.

Die partielle Oxidation

Die Verbrennung eines Brenngases mit unterstöchiometrischer Zusammensetzung wird aufgrund der anteiligen chemischen Umsetzung des Brennstoffes zu Wasser und Kohlenstoffdioxid auch als partielle Oxidation bezeichnet. Es wird zwischen zwei technischen Verfahren zur partiellen Oxidation unterschieden: Der katalytischen partiellen Oxidation und der thermischen partiellen Oxidation. Bei der katalytischen partiellen Oxidation wird das Brenngas in Kontakt mit einem Katalysator gebracht, wodurch die Reaktion zu den jeweiligen Produkten energetisch begünstigt wird und daher bei niedrigen Temperaturen mit hoher Selektivität abläuft. Bei der thermischen Oxidation unter autothermen Bedingungen bildet sich eine Flamme aus. Diese Flamme unterscheidet sich in ihrerer Struktur signifikant von einer Flamme bei der ein Brenngas mit überstöchiometrischer Zusammensetzung verbrannt wird.

Methodik

Eine Möglichkeit zur Simulation turbulenter vorgemischter technischer Verbrennungsysteme zur partiellen Oxidation besteht in der Anwendung einer CFD-Simulation. Dabei werden partielle Differentialgleichungen welche den Transport diverser Modellgrößen beschreiben numerisch gelöst. Die Komplexität solcher Simulationen hängt einerseits von der geometrischen Dimension des betrachteten Rechengebietes und der angestrebten oder erforderlichen Auflösung ab, und andererseits von der Anzahl der zu lösenden Differentialgleichungen sowie deren Abhängigkeiten untereinander.
Bei der fluiddynamischen Simulation eines mehrdimensionalen Rechengebietes werden für die Beschreibung des Impulstransportes die Navier-Stokes-Gleichungen gelöst. Beim Auftreten turbulenter Strömungen sind die auftretenden Turbulenzelemente aufzulösen oder teilweise bzw. vollständig über Modelle zu beschreiben. Wenn das durchströmende Medium, bei den im Rechengebiet auftretenden Druckänderungen, kompressible Effekte aufweist, oder sich die Dichte durch beispielsweise chemische Reaktionen ändert wird eine weitere partielle Differentialgleichung für die Enthalpie oder innere Energie gelöst. Im Falle der chemischen Reaktion werden für die jeweiligen Transportgleichungen der Spezies entsprechend ihrer Anzahl gelöst. Die Quellterme dieser Spezieserhaltungsgleichungen sind über die chemischen Reaktionen der Spezies untereinander gekoppelt. Das hat zur Folge, dass neben den Transporttermen ein gekoppeltes nichtlineares Gleichungssystem zu lösen ist. Bei der überstöchiometrischen Verbrennung von Methan werden, beispielsweise im detaillierten GRI3.0 Mechanismus, 53 Spezies betrachtet. Diese Spezies sind über 325 Reaktionen untereinander gekoppelt. Dieser Umstand führt zur Notwendigkeit der Reduktion des Reaktionsystemes, bei gleichzeitiger Erhaltung der Informationen über die Modellgrößen welche im Fokus der Simulation stehen.
Die in dieser Forschungsarbeit angewendete Methode basiert auf der Reduktion der Spezieserhaltungsgleichungen auf eine für die Reaktion charakteristische Fortschrittvariable und einen die Zusammensetzung beschreibenden Mischungsbruch. Für diese Größen wird dann jeweils eine partielle Differentialgleichung, welche den Transport der jeweiligen Größe beschreibt, gelöst. Der Quellterm der Transportgleichung für die Fortschrittsvariable wird, bei dem zur Anwendung kommenden Modell, einer mehrdimensionalen Tabelle entnommen. In dieser Tabelle sind für jeden Zustand, im Raum von Fortschrittsvariable und Mischungsbruch, die Zusammensetzung, thermophysikalische Daten sowie der Quellterm der Fortschrittsvariablen gespeichert. Zur Erstellung einer solchen Tabelle wird das Reaktormodell einer adiabaten laminaren eindimensionalen Vormischflamme, unter Verwendung eines detaillierten chemischen Reaktionsmechanismus für verschiedene Mischungsbrüche gelöst. Die berechneten Reaktionsverläufe werden über einer definierten Fortschrittsvariable aufgetragen und die Quellterme berechnet. Um den Einfluss der Turbulenz auf die Reaktion abzubilden wird diese Tabelle unter der Verwendung einer Wahrscheinlichkeitsdichtefunktion integriert.

Ziele

Das Ziel dieser Forschungsarbeit ist die Simulation turbulenter vorgemischter technischer Verbrennungssysteme zur partiellen Oxidation. Dafür werden die folgenden Schritte unternommen:

 

  • Identifikation eines geeigneten Reaktionsmechanismus zur Beschreibung der laminaren Vormischflamme im Bereich der partiellen Oxidation
  • Recherche unterschiedlicher Reaktionsmechanismen und Berechnung der laminaren Flammengeschwindigkeit
  • Messung der laminaren Flammengeschwindigkeit als Validierungsparameter
  • Bestimmung einer geeigneten eindeutigen Fortschrittsvariablen als Reaktionsfortschritt
  • Tabellierung laminarer Vormischflammen und Integration mit einer Wahrscheinlichkeitsdichtefunktion
  • Berechnung der isothermen Strömungsfelder eines technischen Systems
  • Berechnung der reaktiven adiabaten Strömungsfelder

 

©

Publikationen


2022
3D direct pore level simulations of radiant porous burners
Wieland, C.; Weis, C.; Habisreuther, P.; Trimis, D.
2022. Combustion and Flame, 245, Art.-Nr.: 112370. doi:10.1016/j.combustflame.2022.112370
Experimental investigation of synthesis gas production in fuel-rich oxy-fuel methane flames
Sentko, M. M.; Schulz, S.; Weis, C.; Stelzner, B.; Anderlohr, C.; Vicari, M.; Trimis, D.
2022. Fuel, 317, Art.-Nr.: 123452. doi:10.1016/j.fuel.2022.123452
2021
Combustion Characterization of Solvents used in Coil Coating Processes: Experiments and Kinetic Modelling
Malliotakis, Z.; Weinbrecht, P.; Cameron, F.; Katsourinis, D.; Vourliotakis, G.; Weis, C. G.; Skevis, G.; Trimis, D.; Founti, M.
2021. 10th European Combustion Meeting (ECM 2021), April 14-15, 2021, Virtual Edition. Proceedings Volume, 1361–1366, MCM
Large Eddy Simulation of the Two-Phase Reactive Flow Field in a Single Sector Laboratory Scale Rich-Quench-Lean Combustion Chamber
Kaddar, D.; Langenthal, T. von; Weis, C. G.; Galeazzo, F. C. C.; Habisreuther, P.; Zarzalis, N.
2021. 10th European Combustion Meeting (ECM 2021), April 14-15, 2021, Virtual Edition. Proceedings Volume, 1121–1126, MCM
Detailed Numerical Simulation for Optimization of Radiation Efficiency of Porous Burners
Wieland, C.; Weinbrecht, P.; Weis, C.; Habisreuther, P.; Trimis, D.
2021. 30. Deutscher Flammentag, Deutsche Sektion des Combustion Institutes und DVV/VDI-Gesellschaft Energie und Umwelt Dinkelacker, Friedrich; Pitsch Heinz; Scherer, Viktor (eds.), 1066–1074
2020
Energy Efficient Coil Coating Process - ECCO
Schneider, M.; Weinbrecht, P.; Wieland, C.; Weis, C.; Trimis, D.
2020. 12th European Conference on Industrial Furnaces and Boilers (INFUB-12) : 10th and 11th November 2020
2019
Determination of a correlation for predicting lean blow off limits of gaseous fueled, premixed turbulent jet flame arrays enclosed in a hexagonal dump combustor
Weis, C. G.; Schwagerus, A.; Faller, S.; Bhagwan, R.; Habisreuther, P.; Zarzalis, N.
2019. 9th European Combustion Meeting (ECM), Lissabon, Portugal, 14 - 17 April 2019
Combustion characteristics of natural gas fueled, premixed turbulent jet flame arrays confined in a hexagonal combustor
Bhagwan, R.; Huth, M.; Schwagerus, A.; Koestlin, B.; Weis, C.; Dederichs, S.; Habisreuther, P.; Zarzalis, N.
2019. ASME Conference Proceedings: Turbo Expo: Power for Land, Sea, and Air. Volume 4A: Combustion, Fuels, and Emissions, The American Society of Mechanical Engineers (ASME). doi:10.1115/GT2019-90286
Determination of dispersion coefficients of heat and mass for porous media by detailed numerical simulation [in press]
Wieland, C.; Zhou, J.; Weis, C.; Habisreuther, P.; Trimis, D.
2019. Seventeenth International Conference on Numerical Combustion, 06-08 May, Aachen, Germany, MS3–199
Entwicklung einer Pilotstrecke zur Demonstration eines innovativen Ofenkonzeptes für die Stahlbandbeschichtung im Rahmen des europäischen Verbundprojektes: Energy Efficient Coil Coating - ECCO [in press]
Weis, C.; Cameron, F.; Weinbrecht, P.; Trimis, D.
2019. Jahrestreffen der ProcessNet Fachgruppe Hochtemperaturtechnik, KIT Karlsruhe, 2.-3. April 2019, Karlsruher Institut für Technologie (KIT)
Entwicklung eines Porenbrenners mit breitem Leistungsspektrum für niederkalorische, wasserstoffhaltige Gase
Weinbrecht, P.; Wieland, C.; Weis, C.; Stelzner, B.; Trimis, D.
2019. 29. Deutscher Flammentag (2019), Bochum, Deutschland, 17.–18. September 2019
Development of a porous burner for low calorific gaseous fuels offering a wide operating range [in press]
Wieland, C.; Weinbrecht, P.; Weis, C.; Habisreuther, P.; Trimis, D.
2019. Proceedings of the European Combustion Meeting – 2019, April 14-17, Lisboa, Portugal, S3_AIII_53
2016
Super-adiabatic flame temperatures in premixed methane flames: A comparison between oxy-fuel and conventional air combustion
Stelzner, B.; Weis, C.; Habisreuther, P.; Zarzalis, N.; Trimis, D.
2016. 1st International Workshop on Oxy-Fuel Combustion, Montabaur, 10.-11.02.2016
2015
Super-adiabatic flame temperatures in premixed methane-oxygen flames
Stelzner, B.; Weis, C.; Habisreuther, P.; Zarzalis, N.; Trimis, D.
2015. Proceedings of the European Combustion Meeting - 2015, March 30-April 2, 2015, Budapest, Hungary : ECM 2015
Bestimmung der laminaren Flammengeschwindigkeit von Methan/Sauerstoff-Gemischen im POX Bereich
Sentko, M. M.; Weis, C.; Habisreuther, P.; Zarzalis, N.; Trimis, D.
2015. Verbrennung und Feuerung : 27.Deutscher Flammentag, Clausthal-Zellerfeld, 16.-17.September 2015, 491–500, VDI Verlag
Temperature dependency of the laminar burning velocity of fuel-rich methane oxygen measurements
Sentko, M. M.; Weis, C.; Harth, S.; Habisreuther, P.; Zarzalis, N.; Trimis, D.
2015. 7th European Combustion Meeting (ECM2015), Budapest, H, March 30 - April 2, 2015
Numerische Untersuchung von super-adiabaten Flammentemperaturen in vorgemischten Methan-Sauerstoff-Flammen
Stelzner, B.; Weis, C.; Habisreuther, P.; Zarzalis, N.; Trimis, D.
2015. Verbrennung und Feuerung : 27.Deutscher Flammentag, Clausthal-Zellerfeld, 16.-17.September 2015, 491–500, VDI Verlag
Super-adiabatic flame temperatures in premixed methane-oxygen flames
Stelzner, B.; Weis, C.; Habisreuther, P.; Zarzalis, N.; Trimis, D.
2015. 7th European Combustion Meeting (ECM2015), Budapest, H, March 30 - April 2, 2015
2014
Measurement of laminar burning velocities of fuel-rich methane-oxygen mixtures
Weis, C.; Sentko, M. M.; Habisreuther, P.; Zarzalis, N.; Trimis, D.
2014. 3rd Heat Flux Burner Workshop, Berlin, September 25th, 2014. Innerhalb von: Flame Structure 2014, 8th International Seminar on Flame Structure, 21.-24. Sept. 2014, Berlin