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

Résumé technique : Modélisation cinétique du reformage à sec de CH4/CO2 dans une décharge à barrière diélectrique

Énoncé du problème
La conversion du méthane (CH4) et du dioxyde de carbone (CO2) via le reformage à sec (DRM) assisté par plasma dans des décharges à barrière diélectrique (DBD) offre une voie non thermique pour produire de l'hydrogène et du gaz de synthèse. Cependant, les outils de modélisation cinétique pour ce processus reposent souvent sur des plateformes spécialisées et non modulaires, et utilisent fréquemment des constantes de vitesse illustratives de l'ordre de grandeur plutôt que des données rigoureusement validées. De plus, il existe un manque de cadres ouverts couplant explicitement la chimie homogène 0D avec le transport par dispersion axiale 1D, ce qui rend difficile le prototypage rapide des paramètres de fonctionnement ou la validation des modèles par rapport à des références expérimentales avec une précision quantitative.

Méthodologie
Les auteurs ont développé un cadre Python ouvert et modulaire (0D/1D) pour simuler la conversion de CH4/CO2 dans une DBD.

• Mécanisme cinétique : Un mécanisme réduit comprenant 23 espèces neutres et 38 réactions a été employé. Celui-ci inclut la dissociation par impact électronique, les réactions bimoléculaires neutre-neutre, les associations à trois corps et les décompositions unimoléculaires.
• Cinétique électronique : Contrare aux itérations précédentes utilisant des ajustements analytiques, les coefficients de vitesse d'impact électronique (E1–E5) ont été dérivés d'une véritable résolution de l'équation de Boltzmann électronique à deux termes (utilisant le solveur BOLOS) alimentée par de réelles sections efficaces de diffusion électron-neutre provenant de la base de données LXCat/Morgan.
• Forçage du plasma : La DBD a été modélisée comme un train d'impulsions périodiques représentant des microdécharges filamentaires. Une approximation de moyenne temporelle a été appliquée pour remplacer les taux instantanés résolus par impulsion par des valeurs moyennées par cycle, validée pour reproduire les résultats résolus par impulsion avec une erreur <1 % tout en réduisant considérablement le coût computationnel.
• Analyse de sensibilité : Une analyse de sensibilité globale de Morris a été réalisée sur les 33 constantes de vitesse neutre-neutre. Elle a identifié six réactions spécifiques comme étant les principaux moteurs de la conversion du CO2, du rapport H2/CO et de l'efficacité énergétique.
• Stratégie de validation : Les six réactions critiques identifiées (R2, R11–R14 et R33) ont été remplacées par des valeurs provenant de la littérature évaluée de la cinétique de combustion et atmosphérique. Les 27 autres constantes de vitesse ont été conservées comme paramètres illustratifs de l'ordre de grandeur.
• Structure du modèle : Le cadre couple un modèle Lagrangien 0D (suivant un paquet de gaz) avec un modèle de réacteur à flux piston à dispersion axiale 1D d'Euler pour simuler les profils spatiaux.

Résultats clés

• Validation du mécanisme : Le remplacement des suppositions illustratives par des valeurs de la littérature a corrigé des erreurs allant de 6 fois à plus de 3000 fois. Plus précisément, la réaction CO + O (+M) → CO2 (+M) (R33), qui est interdite par le spin et difficile à mesurer, a été ancrée à une valeur de la littérature à 200 K (la meilleure disponible, bien que le modèle opère à 400 K).
• Tendances de conversion : Avec le mécanisme validé, la conversion du CO2 respecte le plafond d'environ 25 % observé dans 19 sources de plasma à pression atmosphérique dans la littérature pour une énergie spécifique injectée (SEI) allant jusqu'à ~30 000 J/L. La conversion du CH4 augmente plus brusquement, approchant une conversion quasi complète à des SEI élevées.
• Sélectivité : Le rapport H2/CO se stabilise entre 1,1 et 2,1, reflétant un régime limité par le CH4 où le CH4 se dissocie plus rapidement que le CO2, ce qui est cohérent avec les observations expérimentales.
• Comparaison de référence :
  • À une SEI correspondante de 5,7 eV/molécule, le modèle a initialement surestimé la conversion du CH4 d'environ 5,8 fois et celle du CO2 de ~2,7 fois par rapport à la référence expérimentale d'Ozkan et al. [55].
  • L'application d'une SEI « corrigée par la puissance » (tenant compte de la puissance totale transférée plutôt que seulement du coût énergétique par impact électronique) a réduit ces écarts à environ 3,3 fois (CH4) et 1,6 fois (CO2).
  • Crucialement, à des valeurs de SEI plus basses et typiquement expérimentales (938–2 500 J/L), le modèle a montré un accord étroit avec l'expérience, sous-estimant légèrement la conversion du CH4 et correspondant étroitement au CO2. La divergence augmente uniquement à des SEI élevées, suggérant que le modèle est robuste dans la plage de SEI faible à modérée typique des DBD de laboratoire.
• Sources d'incertitude : L'analyse de sensibilité a confirmé que même après validation, R33 reste la source dominante d'incertitude pour la conversion du CO2 et le rapport H2/CO. La surestimation restante de la conversion du CH4 est attribuée à la chimie du réservoir de radicaux de CH4 non validée (R1, R3–R10) et à l'hypothèse de rapports de branchement résolus par canal pour la dissociation du CH4.

Signification et revendications
L'article présente une étape significative vers un outil ouvert et quantitativement validé pour le DRM par plasma. Ses principales contributions sont :

1. Cadre ouvert : Fournir un code Python reproductible qui couple la chimie 0D avec le transport 1D, accessible pour le prototypage rapide et l'enseignement.
2. Amélioration méthodologique : Passer d'ajustements analytiques arbitraires à de véritables coefficients de vitesse résolus par Boltzmann basés sur de réelles sections efficaces.
3. Validation ciblée : Démontrer que le remplacement par des valeurs de la littérature de seulement six constantes de vitesse clés améliore considérablement la fidélité du modèle, éliminant une source majeure d'erreur et rapprochant les prédictions de la réalité expérimentale.
4. Limites réalistes : Les auteurs affirment modestement que, bien que le modèle ne soit pas encore totalement quantitatif sur toutes les plages de SEI, il fournit une base « qualitative à semi-quantitative » pour le régime de SEI faible à modéré. Ils identifient explicitement la feuille de route restante pour une validation complète : obtenir des données à plus basse température pour R33, valider la chimie des radicaux de CH4 et calibrer les paramètres de forçage du plasma par rapport à des mesures électriques spécifiques du réacteur.

Le travail ne prétend pas avoir résolu quantitativement le problème du DRM, mais propose plutôt une base de calcul transparente et ouverte avec une feuille de route spécifique et resserrée pour parvenir à un modèle entièrement validé.