Numerical Modelling of High-Temperature Methane Pyrolysis Under Plasma-Inspired Thermal Gradients

Technische Zusammenfassung: Numerische Modellierung der Hochtemperatur-Methan-Pyrolyse unter plasma-inspirierten thermischen Gradienten

Problemstellung
Die vorliegende Arbeit befasst sich mit der Herausforderung, abgelegene Erdgasreserven (primär Methan) in transportierbare flüssige Kohlenwasserstoffe umzuwandeln. Während die Mikrowellenplasma-Technologie eine direkte Alternative zu den traditionellen, energieintensiven indirekten Routen (wie der Dampfreformierung von Methan mit anschließender Fischer-Tropsch-Synthese) bietet, steht sie vor erheblichen operativen Hürden. Insbesondere die gleichzeitige Bildung fester Rußpartikel in Hochtemperatur-Plasmazonen führt zu Verstopfungen im Reaktor und einer reduzierten Selektivität für die angestrebten flüssigen Edukte. Darüber hinaus vereinfachen bestehende numerische Simulationen die komplexen thermischen Gradienten in Plasma-Reaktoren oft durch die Annahme einheitlicher thermischer Zonen oder basaler eindimensionaler Approximationen. Diese Vereinfachungen versäumen es, die wahren chemischen Trajektorien der Moleküle zu erfassen, was zu ungenauen Vorhersagen der Produktausbeuten und der Rußbildungsraten führt. Der Kern des Problems ist das Fehlen eines rigorosen, prädiktiven Rahmens, der den hochtemperierten stationären Plasmakern effektiv von der schnellen thermischen Quench-Zone am Reaktorausgang entkoppelt und dabei recheneffizient genug für die mehrdimensionale Kopplung mit der Fluiddynamik bleibt.

Methodik
Die Studie verwendet ein thermisches Pyrolyse-Kinetik-Framework, um die Methankonversion in einer Reaktorumgebung unter atmosphärischem Druck zu modellieren, wobei explizite Elektronenkinetik zugunsten von Neutralradikal-Netzwerken unter aggressiven thermischen Gradienten ausgeschlossen wird.

• Erweiterung des kinetischen Modells: Der grundlegende Mechanismus basiert auf dem Dean-Mechanismus der Methanpyrolyse, erweitert auf das Spektrum von 300–4000 K, wie es für Mikrowellenplasma-Lichtbögen typisch ist. Eine entscheidende Ergänzung ist die Einbeziehung des kollisionsinduzierten Intersystem-Crossing (ISC) von elektronisch angeregtem Singulett-Methylen ($^1CH_2$) zu Triplett-Methylen ($^3CH_2$), das durch Kollisionen mit dem Argon-Badgas getrieben wird. Dieser quantenmechanische Übergang beeinflusst die Entwicklung des Radikalpools signifikant, insbesondere während des schnellen thermischen Quenchs.
• Thermische Profil-Entkopplung: Das Modell unterteilt den Reaktor in zwei distinkte Zonen:
  1. Aktiver Plasmakern: Eine kontinuierliche Heizzone mit einer Verweilzeit von etwa 15,54 ms (in Tabelle 3 als 155,4 s angegeben, vermutlich ein Formatierungsfehler im Originaltext), die zwischen 2500 K und 4000 K gehalten wird, wobei die elektromagnetische Energieabsorption das System in einen Quasi-Fließgleichgewichtszustand hochangeregter Radikalspezies treibt.
  2. Kühlungstrahl-Region (Cooling Jet Region): Eine diffusionslimitierte Quench-Zone, in der das Gas innerhalb von etwa 76 ms (spezifisch 76,9 ms laut Tabelle 3) von 2500 K auf 300 K abkühlt. Dieses Profil wurde durch Lösung einer radialen eindimensionalen Wärmeleitungsgleichung abgeleitet und gegen experimentelle Rayleigh-Streuungs- und OES-Daten validiert.
• Vernachlässigbarkeit der Elektronenkinetik: Eine nulldimensionale Boltzmann-Analyse bestätigt, dass bei 1 bar und 2500 K die thermischen Dissoziationsraten die Elektronenstoß-Dissoziationsraten um sieben Größenordnungen übersteigen. Folglich behandelt das Modell den Prozess als rein thermischen Pyrolyseprozess, der durch neutrale Radikalchemie gesteuert wird.
• Mechanismus-Reduktion: Um die Integration in mehrdimensionale numerische Strömungsmechanik (CFD) zu ermöglichen, wurde ein zweistufiges Reduktionsverfahren mittels KINALC-Software durchgeführt. Dies beinhaltete die skelettale Eliminierung redundanter Reaktionen sowie eine anschließende Hauptkomponentenanalyse (PCA), um die strukturelle Multikollinearität der Reaktionsraten zu adressieren.

Wesentliche Ergebnisse

• Konversion und Ausbeute: Das Modell sagt eine Nettomethankonversion von 98,7 % voraus. Unter der idealisierten Einphasen-Annahme (nur Gasphase) sind die Primärprodukte Wasserstoff und Acetylen (88,0 % Ausbeute), mit kleineren Anteilen an Ethylen (7,0 %), Propen (1,1 %) und Butadien (1,0 %).
• Sequentielle kinetische Erstarrung (Kinetic Freezing): Die Studie identifiziert distinkte thermische Schwellenwerte, die die Produktbildung während des Quenchs steuern:
  • Acetylen stabilisiert sich und erstarrt unterhalb von 3900 K.
  • Ethylen akkumuliert signifikant unterhalb von 1520 K.
  • Propen und Butadien entstehen bei 1060 K bzw. 870 K.
  • Kinetischer Lock: Alle molekularen Transformationen stoppen unterhalb von 500 K, wodurch die endgültige Produktverteilung fixiert wird.
• Auswirkung der Rußdeposition: Eine vergleichende Mehrstufen-Simulation, die die feste Kohlenstoffpräzipitation berücksichtigt, zeigt, dass Gasphasenmodelle die Acetylen-Ausbeute signifikant überbewerten. Wenn fester Kohlenstoff bei 2500 K präzipitieren darf, fällt 99 % der Kohlenstoffmasse als Ruß aus. In dieser kohlenstoffarmen, wasserstoffreichen Umgebung verschiebt sich die nachfolgende Quench-Chemie von ungesättigten Spezies (Acetylen/Ethylen) hin zu stabilem Ethan und unreagiertem Methan – ein Ergebnis, das spezifisch für das Mehrphasenmodell ist und im Gegensatz zur Gasphasen-Baseline steht.
• Recheneffizienz: Die Mechanismus-Reduktion komprimierte das detaillierte Schema (Hunderte von Spezies/Schritten) erfolgreich in einen schlanken Mechanismus aus 26 Spezies und 102 irreversiblen Reaktionen. Während das Abstract eine 61-fache Beschleunigung anführt, berichtet Abschnitt 3.5 explizit von einem 14,2-fachen Beschleunigungsfaktor für das Ausführungsprofil. Dieses reduzierte Modell reproduziert die Spezieskonzentrationen innerhalb einer Abweichung von 0,1 % gegenüber dem vollen Mechanismus.

Bedeutung und Ansprüche
Das Paper beansprucht, die prädiktive Kapazität von numerischen Plasma-Reaktor-Simulationen voranzutreiben, indem es eine präzisere Schätzung der komplexen Temperatur-Zeit-Profile etabliert, denen die reagierenden Gasgemische ausgesetzt sind. Durch die genaue Entkopplung des Hochtemperaturkerns von der schnellen Quench-Zone schafft der Autor einen physikalisch realistischen Rahmen für die detaillierte chemische Kinetikmodellierung.

Die primäre Bedeutung liegt in dem Nachweis, dass das kinetische Einfrieren (Kinetic Freezing) ein sequentieller Prozess ist, der durch spezifische thermische Schwellenwerte gesteuert wird, welche die endgültige Konfiguration der Produkte diktieren. Des Weiteren hebt die Studie hervor, dass heterogene Phasenwechsel (Rußbildung) die Gasphasen-Reaktionspfade fundamental verändern und die Produktverteilungen weg von den durch Einphasenmodelle vorhergesagten thermodynamischen Minima verschieben.

Schließlich bietet die Arbeit einen recheneffizienten Rahmen für zukünftige Ingenieursanwendungen. Der reduzierte Mechanismus, der gegen das vollständige kinetische Schema validiert wurde, ermöglicht hochauflösende, dreidimensionale Simulationen von industriellen Plasma-Reaktoren im vollen Maßstab und überwindet damit die prohibitiven Rechenkosten, die bisher mit solch detaillierter chemischer Kinetik in CFD-Codes verbunden waren.