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

तकनीकी सारांश: प्लाज्मा-प्रेरित तापीय प्रवणता के तहत उच्च-तापमान मीथेन पायरोलिसिस का संख्यात्मक मॉडलिंग

समस्या विवरण
यह शोध पत्र दूरस्थ प्राकृतिक गैस भंडारों (मुख्य रूप से मीथेन) को परिवहन योग्य तरल हाइड्रोकार्बन में परिवर्तित करने की चुनौती को संबोधित करता है। जबकि माइक्रोवेव प्लाज्मा तकनीक पारंपरिक, ऊर्जा-गहन अप्रत्यक्ष मार्गों (जैसे स्टीम मीथेन रिफॉर्मिंग के बाद फिशर-ट्रॉप्स संश्लेषण) के लिए एक सीधा विकल्प प्रदान करती है, यह महत्वपूर्ण परिचालन बाधाओं का सामना करती है। विशेष रूप से, उच्च-तापमान प्लाज्मा क्षेत्रों में ठोस कालिख (सूट) कणों का एक साथ बनना रिएक्टर में रुकावट पैदा करता है और लक्षित तरल फीडस्टॉक के लिए चयनात्मकता को कम करता है। इसके अलावा, मौजूदा संख्यात्मक सिमुलेशन अक्सर प्लाज्मा रिएक्टरों में पाए जाने वाले जटिल तापीय प्रवणता (thermal gradients) का सरलीकरण करते हैं, जो समान तापीय क्षेत्रों या बुनियादी एक-आयामी सन्निकटन को मानते हैं। ये सरलीकरण अणुओं के वास्तविक रासायनिक पथों को पकड़ने में विफल रहते हैं, जिससे उत्पाद प्राप्ति (yield) और कालिख निर्माण दरों की गलत भविष्यवाणियां होती हैं। मुख्य समस्या एक कठोर, भविष्य कहनेवाला ढांचे के अभाव की है जो रिएक्टर निकास पर उच्च-तापमान स्थिर-अवस्था प्लाज्मा कोर और तीव्र तापीय शमन (quenching) क्षेत्र को सटीक रूप से अलग कर सके, और साथ ही बहु-आयामी द्रव गतिकी युग्मन (fluid dynamics coupling) के लिए पर्याप्त रूप से गणनात्मक रूप से कुशल भी रहे।

कार्यप्रणाली
यह अध्ययन वायुमंडलीय-दबाव रिएक्टर वातावरण में मीथेन रूपांतरण को मॉडल करने के लिए एक तापीय पायरोलिसिस गतिज ढांचे (kinetic framework) का उपयोग करता है, जिसमें स्पष्ट इलेक्ट्रॉन गतिकी के बजाय न्यूट्रल रेडिकल नेटवर्क को शामिल किया गया है।

• गतिज मॉडल विस्तार: आधारभूत तंत्र मीथेन पायरोलिसिस के डीन तंत्र (Dean mechanism) पर आधारित है, जिसे माइक्रोवेव प्लाज्मा आर्क के विशिष्ट 300–4000 K तापमान स्पेक्ट्रम को कवर करने के लिए विस्तारित किया गया है। एक महत्वपूर्ण जुड़ाव आर्गन बाथ गैस के साथ टकराव से प्रेरित, इलेक्ट्रॉनिक रूप से उत्तेजित सिंगलेट मेथिलीन ($^1CH_2$) के ट्रिपलेट मेथिलीन ($^3CH_2$) में टकराव-प्रेरित इंटरसिस्टम क्रॉसिंग (ISC) का समावेश है। यह क्वांटम मैकेनिकल संक्रमण रेडिकल पूल के विकास को महत्वपूर्ण रूप से प्रभावित करता है, विशेष रूप से तीव्र तापीय शमन के दौरान।
• तापीय प्रोफाइल का पृथक्करण: मॉडल रिएक्टर को दो विशिष्ट क्षेत्रों में विभाजित करता है:
  1. सक्रिय प्लाज्मा कोर: एक निरंतर तापन क्षेत्र जिसका निवास समय लगभग 15.54 ms (तालिका 3 में 155.4 s के रूप में नोट किया गया है, जो संभवतः मूल पाठ में स्वरूपण त्रुटि है) है, जिसे 2500 K और 4000 K के बीच बनाए रखा जाता है, जहाँ विद्युत चुम्बकीय ऊर्जा अवशोषण प्रणाली को अत्यधिक उत्तेजित रेडिकल प्रजातियों की एक अर्ध-स्थिर अवस्था (quasi-steady state) की ओर ले जाता है।
  2. कूलिंग जेट क्षेत्र: एक विसरण-सीमित शमन क्षेत्र जहाँ गैस लगभग 76 ms (विशेष रूप से तालिका 3 में सूचीबद्ध 76.9 ms) में 2500 K से 300 K तक तेजी से ठंडी होती है। यह प्रोफाइल एक 1D रेडियल ऊष्मा चालन समीकरण को हल करके प्राप्त किया गया है, जिसे रेले स्कैटरिंग और OES डेटा के विरुद्ध सत्यापित किया गया है।
• इलेक्ट्रॉन गतिकी की उपेक्षा: एक शून्य-आयामी बोल्ट्ज़मैन विश्लेषण पुष्टि करता है कि 1 बार और 2500 K पर, तापीय वियोजन दरें इलेक्ट्रॉन-इम्पैक्ट वियोजन दरों से सात गुना अधिक हैं। फलस्वरूप, मॉडल प्रणाली को न्यूट्रल रेडिकल केमिस्ट्री द्वारा नियंत्रित एक शुद्ध तापीय पायरोलिसिस प्रक्रिया के रूप में मानता है।
• तंत्र न्यूनीकरण (Mechanism Reduction): बहु-आयामी कम्प्यूटेशनल फ्लूइड डायनेमिक्स (CFD) में एकीकरण को सक्षम करने के लिए, KINALC सॉफ्टवेयर का उपयोग करके एक दो-चरणीय न्यूनीकरण प्रक्रिया निष्पादित की गई। इसमें रेडंडेंट प्रतिक्रियाओं के कंकाल उन्मूलन (skeletal elimination) के बाद संरचनात्मक बहु-सहसंबंध (multi-collinearity) को संबोधित करने के लिए प्रिंसिपल कंपोनेंट एनालिसिस (PCA) शामिल था।

प्रमुख परिणाम

• रूपांतरण और प्राप्ति: मॉडल 98.7% शुद्ध मीथेन रूपांतरण की भविष्यवाणी करता है। आदर्श एकल-चरण (केवल गैस) धारणा के तहत, प्राथमिक उत्पाद हाइड्रोजन और एसिटिलीन (88.0% प्राप्ति) हैं, जिसमें एथिलीन (7.0%), प्रोपीन (1.1%) और ब्यूटाडिएन (1.0%) के छोटे अंश शामिल हैं।
• अनुक्रमिक गतिज फ्रीजिंग (Sequential Kinetic Freezing): अध्ययन शमन के दौरान उत्पाद निर्माण को नियंत्रित करने वाले विशिष्ट तापीय थ्रेसहोल्ड की पहचान करता है:
  • एसिटिलीन 3900 K से नीचे स्थिर और फ्रीज हो जाती है।
  • एथिलीन 1520 K से नीचे महत्वपूर्ण रूप से संचित होती है।
  • प्रोपीन और ब्यूटाडिएन क्रमशः 1060 K और 870 K पर उभरते हैं।
  • काइनेटिक लॉक: सभी आणविक रूपांतरण 500 K से नीचे समाप्त हो जाते हैं, जिससे अंतिम उत्पाद वितरण स्थिर हो जाता है।
• कालिख जमाव का प्रभाव: ठोस कार्बन अवक्षेपण (precipitation) को शामिल करने वाला एक तुलनात्मक बहु-चरणीय सिमुलेशन प्रकट करता है कि गैस-फेज मॉडल एसिटिलीन की प्राप्ति की काफी अधिक भविष्यवाणी करते हैं। जब 2500 K पर ठोस कार्बन का अवक्षेपण होने दिया जाता है, तो 99% कार्बन द्रव्यमान कालिख (soot) के रूप में बाहर निकल जाता है। इस कार्बन-विहीन, हाइड्रोजन-समृद्ध वातावरण में, बाद का शमन रसायन विज्ञान अंतिम उत्पादों को असंतृप्त प्रजातियों (एसिटिलीन/एथिलीन) से बदलकर स्थिर इथेन और अन-रिएक्टेड मीथेन में स्थानांतरित कर देता है, जो कि एक परिणाम है जो बहु-फेज मॉडल के विशिष्ट है और गैस-फेज बेसलाइन के विपरीत है।
• गणना दक्षता: न्यूनीकरण प्रक्रिया ने विस्तृत योजना (सैकड़ों प्रजातियां/चरण) को 26 प्रजातियों और 102 अपरिवर्तनीय प्रतिक्रियाओं के एक संक्षिप्त तंत्र में सफलतापूर्वक संकुचित कर दिया। जबकि सार (abstract) 61-गुना त्वरण का उल्लेख करता है, अनुभाग 3.5 स्पष्ट रूप से निष्पादन प्रोफाइल के लिए 14.2 गुना स्पीडअप फैक्टर रिपोर्ट करता है। यह घटा हुआ मॉडल पूर्ण तंत्र के 0.1% के भीतर प्रजाति सांद्रता को पुनरुत्पादित करता है।

महत्व और दावे
यह शोध पत्र प्रतिक्रियाशील गैस मिश्रणों द्वारा अनुभव किए जाने वाले जटिल तापमान-समय प्रोफाइल का अधिक कठोर अनुमान स्थापित करके प्लाज्मा रिएक्टर सिमुलेशन की भविष्य कहनेवाला क्षमता को उन्नत करने का दावा करता है। उच्च-तापमान कोर को तीव्र शमन क्षेत्र से सटीक रूप से अलग करके, लेखक विस्तृत रासायनिक गतिज मॉडलिंग के लिए एक भौतिक रूप से यथार्थवादी ढांचा प्रदान करते हैं।

प्राथमिक महत्व इस बात के प्रदर्शन में निहित है कि काइनेटिक फ्रीजिंग एक अनुक्रमिक प्रक्रिया है जो विशिष्ट तापीय थ्रेसहोल्ड द्वारा संचालित होती है, जो अंतिम उत्पाद विन्यास को निर्धारित करती है। इसके अलावा, अध्ययन इस बात पर प्रकाश डालता है कि विषम चरण परिवर्तन (कालिख निर्माण) मौलिक रूप से गैस-फेज प्रतिक्रिया पथों को बदल देते हैं, जिससे उत्पाद वितरण एकल-फेज मॉडल द्वारा अनुमानित ऊष्मागतिक न्यूनतम (thermodynamic minima) से हट जाता है।

अंततः, यह कार्य भविष्य के इंजीनियरिंग अनुप्रयोगों के लिए एक गणनात्मक रूप से कुशल ढांचा प्रदान करता है। पूर्ण तंत्र के विरुद्ध सत्यापित किया गया घटा हुआ तंत्र, पूर्ण-स्तरीय औद्योगिक प्लाज्मा रिएक्टरों के उच्च-निष्ठा (high-fidelity), तीन-आयामी सिमुलेशन को सक्षम बनाता है, जो फ्लूइड डायनेमिक्स कोड में इतने विस्तृत रासायनिक गतिकी के साथ जुड़े अत्यधिक गणनात्मक खर्चों की बाधा को दूर करता है।