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


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
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
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
Microstructure, local dynamics, and flow behavior of colloidal suspensions with weak attractive interactions.
Weis, C.; Oelschlaeger, C.; Dijkstra, D.; Ranft, M.; Willenbacher, N.
2016. Scientific reports, 6, Art.Nr.:33498. doi:10.1038/srep33498
Enhancing 4D PC-MRI in an aortic phantom considering numerical simulations.
Kratzke, J.; Schoch, N.; Weis, C.; Müller-Eschner, M.; Speidel, S.; Farag, M.; Beller, C.; Heuvelinea, V.
2015. SPIE Medical Imaging 2015: Physics of Medical Imaging, 22 - 25 February 2015, Orlando, Florida, United States. Ed.: C.Hoeschen, Art.Nr. 94121E, Society of Photo-optical Instrumentation Engineers (SPIE). doi:10.1117/12.2082483
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
Effect of weak attractive interactions on flow behavior of highly concentrated crystalline suspensions.
Weis, C.; Natalia, I.; Willenbacher, N.
2014. Journal of Rheology, 58 (5), 1583–1597. doi:10.1122/1.4882955