Kinetic Modelling of CH4/CO2 Dry Reforming in a Dielectric Barrier Discharge for Hydrogen Production: An Open, Sensitivity- Guided Validation Study

Technische Zusammenfassung: Kinetische Modellierung der CH4/CO2-Trockenreformierung in einer dielektrischen Barriereentladung

Problemstellung
Die Konversion von Methan (CH4) und Kohlendioxid (CO2) mittels plasmaunterstützter Trockenreformierung (DRM) in dielektrischen Barriereentladungen (DBD) bietet einen nicht-thermischen Pfad zur Produktion von Wasserstoff und Synthesegas. Kinematische Modellierungswerkzeuge für diesen Prozess stützen sich jedoch häufig auf spezialisierte, nicht-modulare Plattformen und verwenden oft illustrative Größenordnungs-Ratenkonstanten anstelle von streng validierten Daten. Zudem mangelt es an offenen Frameworks, die eine 0D-homogene Chemie explizit mit einem 1D-axialen Dispersions-Transport koppeln, was es schwierig macht, Betriebsparameter schnell zu prototypisieren oder Modelle mit quantitativer Präzision gegen experimentelle Benchmarks zu validieren.

Methodik
Die Autoren entwickelten ein offenes, modulares Python-Framework (0D/1D), um die CH4/CO2-Konversion in einer DBD zu simulieren.

• Kinetischer Mechanismus: Es wurde ein reduzierter Mechanismus verwendet, der aus 23 neutralen Spezies und 38 Reaktionen besteht. Dieser umfasst Elektronenstoß-Dissoziation, bimolekulare Neutral-Neutral-Reaktionen, Drei-Körper-Assoziationen sowie unimolekulare Zerfälle.
• Elektronenkinetik: Im Gegensatz zu früheren Iterationen, die auf analytischen Fits basierten, wurden die Elektronenstoß-Ratenkoeffizienten (E1–E5) aus einer echten Lösung der Zwei-Term-Elektronen-Boltzmann-Gleichung (unter Verwendung des BOLOS-Solvers) abgeleitet, die mit realen Elektronen-Neutral-Streuquerschnitten aus der LXCat/Morgan-Datenbank gespeist wurde.
• Plasma-Forcing: Die DBD wurde als periodische Pulsfolge modelliert, die filamentäre Mikroentladungen repräsentiert. Eine Zeitmittelungs-Approximation wurde angewandt, um instantane pulsaufgelöste Raten durch zyklusgemittelte Werte zu ersetzen, was validiert wurde, indem ein Fehler von <1 % bei der Reproduktion der pulsaufgelösten Ergebnisse erzielt wurde, während die Rechenkosten signifikant gesenkt wurden.
• Sensitivitätsanalyse: Eine globale Morris-Sensitivitätsanalyse wurde an den 33 Neutral-Neutral-Ratenkonstanten durchgeführt. Diese identifizierte sechs spezifische Reaktionen als die primären Treiber für die CO2-Konversion, das H2/CO-Verhältnis und die Energieeffizienz.
• Validierungsstrategie: Die sechs identifizierten kritischen Reaktionen (R2, R11–R14 und R33) wurden durch Werte aus der evaluierten Verbrennungs- und Atmosphärenkinetik-Literatur ersetzt. Die verbleibenden 27 Ratenkonstanten wurden als illustrative Größenordnungs-Parameter beibehalten.
• Modellstruktur: Das Framework koppelt ein 0D-Lagrange-Modell (das ein Gaspaket verfolgt) mit einem 1D-Euler-Axial-Dispersions-Pfropferohrreaktor-Modell, um räumliche Profile zu simulieren.

Hauptergebnisse

• Mechanismus-Validierung: Das Ersetzen von illustrativen Schätzungen durch Literaturwerte korrigierte Fehler, die von dem 6-fachen bis zum über 3000-fachen reichten. Insbesondere wurde die Reaktion CO + O (+M) → CO2 (+M) (R33), die spin-verboten und schwer messbar ist, auf einen Literaturwert bei 200 K (der beste verfügbare, obwohl das Modell bei 400 K arbeitet) verankert.
• Konversionstrends: Mit dem validierten Mechanismus respektiert die CO2-Konversion die in der Literatur beobachtete Decke von ~25 % für spezifische Energiezufuhr (SEI) bis zu ~30.000 J/L über 19 atmosphärische Plasmaquellen hinweg. Die CH4-Konversion steigt steiler an und nähert sich bei hoher SEI einer nahezu vollständigen Konversion an.
• Selektivität: Das H2/CO-Verhältnis stabilisiert sich zwischen 1,1 und 2,1, was ein CH4-limitiertes Regime widerspiegelt, in dem CH4 schneller dissoziiert als CO2, was im Einklang mit experimentellen Beobachtungen steht.
• Benchmark-Vergleich:
  • Bei einer angeglichenen SEI von 5,7 eV/Molekül überschätzte das Modell die CH4-Konversion zunächst um das ~5,8-fache und die CO2-Konversion um das ~2,7-fache im Vergleich zum experimentellen Benchmark von Ozkan et al. [55].
  • Die Anwendung einer „leistungskorrigierten“ SEI (die die gesamte Leistungsübertragung statt nur der Elektronenstoß-Energiekosten berücksichtigt) reduzierte diese Abweichungen auf ~3,3-fach (CH4) und ~1,6-fach (CO2).
  • Entscheidend ist, dass das Modell bei niedrigeren, experimentell typischen SEI-Werten (938–2.500 J/L) eine enge Übereinstimmung mit dem Experiment zeigte, wobei es die CH4-Konversion leicht unterschätzte und die CO2-Konversion gut traf. Die Divergenz nimmt erst bei hoher SEI zu, was darauf hindeutet, dass das Modell im Labor-typischen Bereich niedriger bis moderater SEI robust ist.
• Unsicherheitsquellen: Die Sensitivitätsanalyse bestätigte, dass selbst nach der Validierung R33 die dominante Quelle der Unsicherheit für die CO2-Konversion und das H2/CO-Verhältnis bleibt. Die verbleibende Überschätzung der CH4-Konversion wird der unvalidierten CH4-Radikal-Pool-Chemie (R1, R3–R10) und der Annahme von kanalaufgelösten Verzweigungsverhältnissen für die CH4-Dissoziation zugeschrieben.

Bedeutung und Ansprüche
Das Paper stellt einen bedeutenden Schritt in Richtung eines quantitativ validierten Open-Source-Tools für das Plasma-DRM dar. Seine primären Beiträge sind:

1. Offenes Framework: Bereitstellung eines reproduzierbaren Python-Codebase, das 0D-Chemie mit 1D-Transport koppelt, zugänglich für schnelles Prototyping und Lehre.
2. Methodische Aufwertung: Der Übergang von willkürlichen analytischen Fits zu echten, durch den Boltzmann-Solver ermittelten Ratenkoeffizienten basierend auf realen Querschnitten.
3. Gezielte Validierung: Demonstration, dass der Sensitivitäts-gesteuerte Ersatz von nur sechs Schlüssel-Ratenkonstanten durch Literaturwerte die Modelltreue signifikant verbessert, die Fehlerquelle eliminiert und die Vorhersagen näher an die experimentelle Realität bringt.
4. Realistische Einschränkungen: Die Autoren behaupten bescheiden, dass das Modell zwar noch nicht vollständig quantitativ über alle SEI-Bereiche hinweg ist, aber eine „qualitative bis semi-quantitative“ Basis für den niedrigen bis moderaten SEI-Bereich bietet. Sie identifizieren explizit den verbleibenden Fahrplan für eine vollständige Validierung: die Gewinnung von Niedrigtemperaturdaten für R33, die Validierung der CH4-Radikalchemie und die Kalibrierung der Plasma-Forcing-Parameter gegen spezifische elektrische Messungen des Reaktors.

Die Arbeit beansprucht nicht, das DRM-Problem quantitativ gelöst zu haben, sondern bietet vielmehr eine transparente, offene rechnerische Basis mit einem spezifisch eingegrenzten Fahrplan zur Erreichung eines vollständig validierten Modells.