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

Technische Samenvatting: Numerieke Modellering van Methaanpyrolyse bij Hoge Temperaturen onder Plasma-geïnspireerde Thermische Gradiënten

Probleemstelling
Het artikel behandelt de uitdaging om afgelegen aardgasreserves (voornamelijk methaan) om te zetten in transporteerbare vloeibare koolwaterstoffen. Hoewel microgolfplasma-technologie een direct alternatief biedt voor traditionele, energie-intensieve indirecte routes (zoals stoomreforming van methaan gevolgd door Fischer-Tropsch-synthese), wordt het geconfronteerd met aanzienlijke operationele hindernissen. Specifiek leidt de gelijktijdige vorming van vaste roetdeeltjes in hoogtemperatuur-plasmazones tot verstopping van de reactor en een verminderde selectiviteit voor de beoogde vloeibare grondstoffen. Bovendien vereenvoudigen bestaande numerieke simulaties de complexe thermische gradiënten in plasma-reactoren vaak door uniforme thermische zones of basismatige eendimensionale benaderingen aan te nemen. Deze vereenvoudigingen falen in het vastleggen van de werkelijke chemische trajecten van moleculen, wat leidt tot onnauwkeurige voorspellingen van productopbrengsten en roetvormingssnelheden. Het kernprobleem is het gebrek aan een rigoureus, voorspellend kader dat de hoogtemperatuur-steady-state plasmakern effectief kan ontkoppelen van de snelle thermische quench-zone bij de reactoruitgang, terwijl het computationeel efficiënt genoeg blijft voor koppeling met meerdimensionale fluïdumdynamica.

Methodologie
De studie maakt gebruik van een thermisch pyrolyse-kinetisch kader om de methaanconversie in een reactoromgeving bij atmosferische druk te modelleren, waarbij expliciete elektronenkinetica bewust wordt uitgesloten ten gunste van neutrale radicaalnetwerken onder agressieve thermische gradiënten.

• Expansie van het Kinetisch Model: Het fundamentele mechanisme is gebaseerd op het Dean-mechanisme voor methaanpyrolyse, uitgebreid om het temperatuurspectrum van 300–4000 K te dekken dat typerend is voor microgolfplasma-bogen. Een cruciale toevoeging is de inclusie van collision-induced intersystem crossing (ISC) van elektronisch geëxciteerde singlet methyleen ($^1CH_2$) naar triplet methyleen ($^3CH_2$), gedreven door botsingen met het argon badgas. Deze kwantummechanische transitie beïnvloedt de ontwikkeling van de radicalenpool aanzienlijk, met name tijdens de snelle thermische quench.
• Decoupling van het Thermisch Profiel: Het model scheidt de reactor in twee afzonderlijke zones:
  1. Actieve Plasma Kern: Een continue verwarmingszone met een verblijftijd van ongeveer 15,54 ms (genoteerd als 155,4 s in Tabel 3, waarschijnlijk een formatteringsfout in de brontekst) die wordt onderhouden tussen 2500 K en 4000 K, waar elektromagnetische energieabsorptie het systeem naar een quasi-steady-state van sterk geëxciteerde radicale soorten drijft.
  2. Koelstraal Regio: Een diffusie-gelimiteerde quench-zone waarin het gas snel afkoelt van 2500 K naar 300 K over ongeveer 76 ms (specifiek 76,9 ms zoals vermeld in Tabel 3). Dit profiel is afgeleid door het oplossen van een 1D radiale warmtegeleidingsvergelijking, gevalideerd tegen experimentele Rayleigh-verstrooiing en OES-data.
• Verwaarloosbaarheid van Elektronenkinetica: Een zero-dimensionale Boltzmann-analyse bevestigt dat bij 1 bar en 2500 K de thermische dissociatiesnelheden de elektron-impact dissociatiesnelheden met zeven ordes van grootte overtreffen. Bijgevolg behandelt het model het systeem als een puur thermisch pyrolyseproces dat wordt beheerst door neutrale radicalenchemie.
• Mechanisme-reductie: Om integratie in meerdimensionale Computationele Fluïdumdynamica (CFD) mogelijk te maken, werd een tweestaps-reductieprocedure uitgevoerd met KINALC-software. Dit omvatte de skeletale eliminatie van redundante reacties, gevolgd door Principal Component Analysis (PCA) om structurele multi-collineariteit in reactiesnelheden aan te pakken.

Belangrijkste Resultaten

• Conversie en Opbrengst: Het model voorspelt een netto methaanconversie van 98,7%. Onder de geïdealiseerde single-phase (gas-alleen) aanname zijn de primaire producten waterstof en acetylen (88,0% opbrengst), met kleinere fracties etheen (7,0%), propeen (1,1%) en butadieen (1,0%).
• Sequentiële Kinetische Bevriezing: De studie identificeert onderscheidende thermische drempelwaarden die de productvorming tijdens de quench beheersen:
  • Acetylen stabiliseert en bevriest onder 3900 K.
  • Etheen accumuleert significant onder 1520 K.
  • Propeen en butadieen ontstaan bij respectievelijk 1060 K en 870 K.
  • Kinetische Lock: Alle moleculaire transformaties stoppen onder 500 K, waardoor de uiteindelijke productdistributie wordt vastgelegd.
• Impact van Roetdepositie: Een vergelijkende meerfasen-simulatie die de precipitatie van vast koolstof meenemen, laat zien dat gasfase-modellen de acetylenopbrengst aanzienlijk overschatten. Wanneer vast koolstof bij 2500 K kan precipiteren, valt 99% van de koolstofmassa uit als roet. In deze koolstof-arme, waterstofrijke omgeving verschuift de daaropvolgende quench-chemie de uiteindelijke producten van onverzadigde soorten (acetylen/etheen) naar stabiel ethaan en niet-gereageerd methaan, een resultaat dat specifiek is voor het meerfasenmodel en contrasteert met de gasfase-baseline.
• Computationele Efficiëntie: De reductie van het mechanisme slaagde erin het gedetailleerde schema (honderden soorten/stappen) te comprimeren tot een slank mechanisme van 26 soorten en 102 onomkeerbare reacties. Terwijl de abstract een 61-voudige versnelling vermeldt, rapporteert Sectie 3.5 expliciet een 14,2 keer snelheidstoename voor het executieprofiel. Dit gereduceerde model reproduceert de concentraties van soorten binnen 0,1% van het volledige mechanisme.

Betekenis en Claims
Het artikel beweert de voorspellende capaciteit van numerieke plasmasimulaties te verbeteren door een striktere schatting te geven van de complexe temperatuur-tijdprofielen die gereageerde gasmengsels ervaren. Door de hoogtemperatuurkern accuraat te ontkoppelen van de snelle quench-zone, biedt de auteur een fysiek realistisch kader voor gedetailleerde chemische kinetische modellering.

De primaire betekenis ligt in de demonstratie dat kinetische bevriezing een sequentieel proces is dat wordt beheerst door specifieke thermische drempelwaarden, wat de uiteindelijke productconfiguratie bepaalt. Verder benadrukt de studie dat heterogene faseveranderingen (roetvorming) de gasfase-reactiepaden fundamenteel veranderen, waardoor de productdistributies verschuiven van de thermodynamische minima die door single-phase modellen worden voorspeld.

Ten slotte biedt het werk een computationeel efficiënt kader voor toekomstige engineeringtoepassingen. Het gereduceerde mechanisme, gevalideerd tegen het volledige kinetische schema, maakt hoogwaardige, driedimensionale simulaties van industriële plasma-reactoren op volledige schaal mogelijk, waarbij de prohibitieve computationele kosten die voorheen geassocieerd werden met dergelijke gedetailleerde chemische kinetiek in fluïdumdynamica-codes, worden overwonnen.