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

Technische Samenvatting: Kinetische Modellering van CH4/CO2 Droge Reforming in een Dielektrische Barrièreontlading

Probleemstelling
De conversie van methaan (CH4) en kooldioxide (CO2) via plasma-ondersteunde droge reforming (DRM) in dielektrische barrièreontladingen (DBD's) biedt een niet-thermisch pad voor de productie van waterstof en syngas. Kinetische modelleringsinstrumenten voor dit proces vertrouwen echter vaak op gespecialiseerde, niet-modulaire platforms en maken veelal gebruik van illustratieve grootteorde-orde van grootte reactieconstanten in plaats van rigoureus gevalideerde gegevens. Bovendien is er een gebrek aan open frameworks die 0D homogene chemie expliciet koppelen aan 1D axiale dispersietransport, wat het lastig maakt om operationele parameters snel te prototypen of modellen te valideren tegen experimentele benchmarks met kwantitatieve precisie.

Methodologie
De auteurs hebben een open, modulair Python-framework (0D/1D) ontwikkeld om de CH4/CO2-conversie in een DBD te simuleren.

• Kinetisch Mechanisme: Een gereduceerd mechanisme bestaande uit 23 neutrale soorten en 38 reacties werd gehanteerd. Dit omvat elektron-impact dissociatie, bimoleculaire neutraal-neutraal reacties, drie-lichaam associaties en unimoleculaire ontledingen.
• Elektronenkinetiek: In tegen tegenstelling tot eerdere iteraties die gebruikmaakten van analytische fits, werden de elektron-impact reactiecoëfficiënten (E1–E5) afgeleid van een echte twee-term elektron Boltzmann-vergelijking oplossing (met behulp van de BOLOS-solver), gevoed met echte elektron-neutraal verstrooiingsdoorsneden uit de LXCat/Morgan-database.
• Plasma-forcing: De DBD werd gemodelleerd als een periodieke pulsreeks die filamentaire microontladingen representeert. Een tijdgemiddelde benadering werd toegepast om instantane, puls-opgeloste snelheden te vervangen door cyclus-gemiddelde waarden, wat gevalideerd is door de puls-opgeloste resultaten met <1% fout weer te geven, terwijl de computationele kosten aanzienlijk werden verminderd.
• Gevoeligheidsanalyse: Er werd een globale Morris-gevoeligheidsanalyse uitgevoerd op de 33 neutraal-neutraal reactieconstanten. Dit identificeerde zes specifieke reacties als de primaire drijvers van de CO2-conversie, de H2/CO-ratio en de energie-efficiëntie.
• Validatiestrategie: De zes geïdentificeerde kritische reacties (R2, R11–R14 en R33) werden vervangen door waarden uit geëvalueerde verbrandings- en atmosferische kinetiekliteratuur. De resterende 27 reactieconstanten werden behouden als illustratieve grootteorde-parameters.
• Modelstructuur: Het framework koppelt een 0D Lagrangiaans model (dat een gaspakketje volgt) aan een 1D Euleriaans axiaal-dispersie plugflow-reactor model om ruimtelijke profielen te simuleren.

Belangrijkste Resultaten

• Mechanisme Validatie: Het vervangen van illustratieve gissingen door literatuurwaarden corrigeerde fouten variërend van 6-voudig tot meer dan 3000-voudig. Specifiek werd de reactie CO + O (+M) → CO2 (+M) (R33), die spin-verboden is en moeilijk te meten is, verankerd aan een literatuurwaarde bij 200 K (de beste beschikbare, hoewel het model bij 400 K opereert).
• Conversietrends: Met het gevalideerde mechanisme respecteert de CO2-conversie de plafondwaarde van ~25% die wordt waargenomen bij 19 atmosferische druk plasma-bronnen in de literatuur voor een Specifieke Energie-input (SEI) tot ~30.000 J/L. De CH4-conversie stijgt steiler en nadert bijna volledige conversie bij een hoge SEI.
• Selectiviteit: De H2/CO-ratio stabiliseert tussen 1,1 en 2,1, wat een CH4-beperkt regime reflecteert waarin CH4 sneller dissocieert dan CO2, consistent met experimentele waarnemingen.
• Benchmark Vergelijking:
  • Bij een overeenkomstige SEI van 5,7 eV/molecuul overschatte het model aanvankelijk de CH4-conversie met ~5,8-voudig en de CO2-conversie met ~2,7-voudig vergeleken met de Ozkan et al. [55] experimentele benchmark.
  • Het toepassen van een "kracht-gecorrigeerde" SEI (rekening houdend met de totale energieoverdracht in plaats van alleen de elektron-impact energiekosten) verkleinde deze kloven tot ~3,3-voudig (CH4) en ~1,6-voudig (CO2).
  • Cruciaal is dat bij lagere, experimenteel typische SEI-waarden (938–2.500 J/L) het model nauwe overeenstemming vertoonde met het experiment, waarbij de CH4-conversie licht werd onderschat en de CO2-conversie nauwkeurig werd benaderd. De divergentie neemt alleen toe bij een hoge SEI, wat suggereert dat het model robuust is in het laag-tot-matig SEI-regime dat typerend is voor laboratorium DBD's.
• Onzekerheidsbronnen: De gevoeligheidsanalyse bevestigde dat zelfs na validatie R33 de dominante bron van onzekerheid blijft voor de CO2-conversie en de H2/CO-ratio. De resterende overschatting van de CH4-conversie wordt toegeschreven aan ongevalideerde CH4-radicaal-pool chemie (R1, R3–R10) en de aanname van kanaal-opgeloste vertakkingsratio's voor de CH4-dissociatie.

Betekenis en Claims
Het artikel vormt een belangrijke stap naar een kwantitatief gevalideerde, open-source tool voor plasma DRM. De belangrijkste bijdragen zijn:

1. Open Framework: Het bieden van een reproduceerbare Python-codebase die 0D chemie koppelt aan 1D transport, toegankelijk voor snelle prototyping en onderwijs.
2. Methodologische Upgrade: De overstap van willekeurige analytische fits naar echte Boltzmann-opgeloste reactiecoëfficiënten gebaseerd op echte doorsneden.
3. Gerichte Validatie: Het aantonen dat het vervangen van slechts zes sleutelreactieconstanten met literatuurwaarden de modelgetrouwheid aanzienlijk verbetert, waardoor een belangrijke foutbron wordt geëlimineerd en voorspellingen dichter bij de experimentele realiteit worden gebracht.
4. Realistische Beperkingen: De auteurs claimen bescheiden dat hoewel het model nog niet volledig kwantitatief is over alle SEI-bereiken, het een "kwalitatieve-tot-semi-kwantitatieve" basis biedt voor het laag-tot-matig SEI-regime. Ze identificeren expliciet de resterende roadmap voor volledige validatie: het verkrijgen van data bij lagere temperaturen voor R33, het valideren van de CH4-radicaalchemie en het kalibreren van plasma-forcing parameters tegen specifieke reactor elektrische metingen.

Het werk claimt niet het DRM-probleem kwantitatief te hebben opgelost, maar biedt eerder een transparante, open computationele basis met een specifieke, verfijnde roadmap om een volledig gevalideerd model te bereiken.