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

ملخص تقني: النمذجة العددية للتحلل الحراري للميثان عند درجات حرارة عالية تحت تأثير تدرجات حرارية مستوحاة من البلازما

بيان المشكلة
تتناول الورقة البحثية تحدي تحويل احتياطيات الغاز الطبيعي البعيدة (المكونة أساساً من الميثان) إلى هيدروكربونات سائلة قابلة للنقل. وبينما توفر تقنية بلازما الموجات الدقيقة (microwave plasma) بديلاً مباشراً لمسارات التخليق غير المباشرة المستهلكة للطاقة (مثل تخليق الميثان بالبخام متبوعاً بتخليق فيشر-تروبش)، إلا أنها تواجه عقبات تشغيلية كبيرة. وتحديداً، يؤدي التكون المتزامن لجسيمات السخام الصلبة في مناطق البلازما عالية الحرارة إلى انسداد المفاعلات وتقليل الانتقائية للمواد اللقيم السائلة المستهدفة. علاوة على على ذلك، فإن المحاكاة العددية الحالية غالباً ما تبسط المناطق الحرارية المعقدة الموجودة في مفاعلات البلازما، بافتراض وجود مناطق حرارية موحدة أو تقريبات أحادية البعد بسيطة. وتفشل هذه التبسيطات في رصد المسارات الكيميائية الحقيقية للجزيئات، مما يؤدي إلى تنبؤات غير دقيقة لعوائد المنتجات ومعدلات تكون السخام. وتتمثل المشكلة الجوية في غياب إطار عمل تنبؤي صارم يمكنه بدقة فصل قلب البلازما المستقر في درجات الحرارة العالية عن منطقة التبريد السريع (quenching) عند مخرج المفاعل، مع البقاء كفؤاً حسابياً بما يكفي للربط مع ديناميكا السوائل متعددة الأبعاد.

المنهجية
تستخدم الدراسة إطار عمل حركي للتحلل الحراري لنمذجة تحول الميثان في بيئة مفاعل يعمل عند الضغط الجوي، مع استبعاد حركية الإلكترونات الصريحة لصالح شبكات الجذور الحرة المتعادلة تحت تدرجات حرارية شديدة.

• توسيع النموذج الحركي: يعتمد الميكانزم الأساسي على ميكانزم "دين" (Dean mechanism) للتحلل الحراري للميثان، والذي تم توسيعه ليشمل نطاق درجات حرارة يتراوح بين 300 و4000 كلفن، وهي الدرجات النموذجية لأقواس بلازما الموجات الدقيقة. والإضافة الحرجة هي تضمين العبور بين الأنظمة (intersystem crossing - ISC) المستحث بالاصطدام للميثيلين ذو الحالة السينجلت المثارة إلكترونياً ($^1CH_2$) إلى الميثيلين ذو الحالة تريبلت ($^3CH_2$)، والمدفوع بالاصطدامات مع غاز الأرجون الخامل. هذا الانتقال الميكانيكي الكمي يؤثر بشكل كبير على تطور مخزن الجذور الحرة، لا سيما أثناء عملية التبريد الحراري السريع.
• فك الارتباط للملف الحراري: يقوم النموذج بفصل المفاعل إلى منطقتين متميزتين:
  1. قلب البلازما النشط: منطقة تسخين مستمرة بزمن بقاء يبلغ حوالي 15.54 مللي ثانية (ذُكرت كـ 155.4 ثانية في الجدول 3، ومن المرجح أن يكون خطأ في تنسيق النص المصدر) يتم الحفاظ عليها بين 2500 و4000 كلفن، حيث يعمل امتصاص الطاقة الكهرومغناطيسية على دفع النظام إلى حالة شبه مستقرة من الأنواع الجذرية عالية الإثارة.
  2. منطقة نفث التبريد: منطقة تبريد سريع (quench) محكومة بالانتشار، حيث تخضع الغازات لتبريد سريع من 2500 كلفن إلى 300 كلفن خلال حوالي 76 مللي ثانية (تحديداً 76.9 مللي ثانية كما ورد في الجدول 3). تم اشتقاق هذا الملف عبر حل معادلة التوصيل الحراري الشعاعي أحادية البعد، والتي تم التحقق من صحتها باستخدام بيانات تشتت رايلي والتحليل الطيفي للانبعاث الضوئي (OES).
• إهمال حركية الإلكترونات: أكد تحليل بولتزمانما (Boltzmann analysis) صفر الأبعاد أنه عند ضغط 1 بار ودرجة حرارة 2500 كلفن، تتجاوز معدلات التفكك الحراري معدلات التفكك الناجمة عن اصطدام الإلكترونات بسبعة مراتب عشرية. وبناءً على ذلك، يعامل النموذج العملية كعملية تحلل حراري نقي محكومة بكيمياء الجذور الحرة المتعادلة.
• اختزال الميكانزم: لتمكين التكامل في ديناميكا السوائل الحسابية (CFD) متعددة الأبعاد، تم تنفيذ إجراء اختزال من خطوتين باستخدام برنامج KINALC. تضمن ذلك الاستبعاد الهيكلي للتفاعلات الزائدة، متبوعاً بتحليل المكونات الرئيسية (PCA) لمعالجة تعدد الخطية الهيكلي في معدلات التفاعل.

النتائج الرئيسية

• التحول والعائد: يتوقع النموذج صافي تحول للميثان بنسبة 98.7%. وتحت افتراض الطور الواحد المثالي (الغاز فقط)، تكون المنتجات الأساسية هي الهيدروجين والأسيتيلين (عائد 88.0%)، مع حصص أصغر من الإيثيلين (7.0%)، والبروبين (1.1%)، والبيوتاديين (1.0%).
• التجميد الحركي المتتالي: تحدد الدراسة عتبات حرارية متميزة تحكم تكوين المنتجات أثناء التبريد:
  • يستقر الأسيتيلين ويتجمد تحت 3900 كلفن.
  • يتراكم الإيثيلين بشكل ملحوظ تحت 1520 كلفن.
  • يظهر البروبين والبيوتاديين عند 1060 كلفن و870 كلفن على التوالي.
  • القفل الحركي: تتوقف جميع التحولات الجزيئية تحت 500 كلفن، مما يثبت التوزيع النهائي للمنتجات.
• تأثير ترسب السخام: كشفت محاكاة متعددة المراحل تتضمن ترسيب الكربون الصلب أن نماذج الطور الغازي تبالغ في تقدير عوائد الأسيتيلين. فعند السماح بترسيب الكربون الصلب عند 2500 كلفن، يسقط 99% من كتلة الكربون كـسخام. وفي هذه البيئة الغنية بالهيدروجين والمستنفدة من الكربون، تتحول كيمياء التبريد اللاحقة من الأنواع غير المشبعة (الأسيتيلين/الإيثيلين) إلى الإيثان المستقر والميثان غير المتفاعل، وهي نتيجة خاصة بنموذج متعدد الأطوار وتتعارض مع نموذج الطور الغازي الأساسي.
• الكفاءة الحسابية: نجح اختزال الميكانزم في ضغط المخطط التفصيلي (مئات الأنواع/الخطوات) إلى ميكانزم رشيق يتكون من 26 نوعاً و102 تفاعل غير عكسي. وبينما يشير الملخص إلى تسريع بمقدار 61 ضعفاً، ذكر القسم 3.5 صراحةً عامل تسريع قدره 14.2 مرة لملف التنفيذ. ويعيد هذا النموذج المختزل إنتاج تركيزات الأنواع بدقة تصل إلى 0.1% من الميكانزم الكامل.

الأهمية والادعاءات
تدعي الورقة أنها تساهم في تطوير القدرة التنبؤية لمحاكاة مفاعلات البلازما العددية من خلال وضع تقدير أكثر صرامة لملفات درجة الحرارة والزمن المعقدة التي تمر بها مخاليط الغاز المتفاعلة. ومن خلال الفصل الدقيق بين قلب الحرارة العالية ومنطقة التبريد السريع، يوفر المؤلف إطار عمل واقعي فيزيائياً للنمذجة الكيميائية الحركية التفصيلية.

تكمن الأهمية الرئيسية في إثبات أن التجميد الحركي هو عملية متسلسلة محكومة بعتبات حرارية محددة، وهي التي تملي التكوين النهائي للمنتجات. علاوة على ذلك، تسلط الدراسة الضوء على أن تغيرات الطور غير المتجانسة (تكون السخام) تغير جذرياً مسارات تفاعل الطور الغازي، مما يحول توزيع المنتجات بعيداً عن الحدود الدنيا الديناميكية الحرارية التي تتنبأ بها نماذج الطور الواحد.

أخيراً، يقدم العمل إطار عمل كفؤ حسابياً للتطبيقات الهندسية المستقبلية. إن الميكانزم المختزل، الذي تم التحقق من صحته مقابل الميكانزم الكامل، يتيح إجراء عمليات محاكاة ثلاثية الأبعاد عالية الدقة لمفاعلات البلازما الصناعية كاملة النطاق، متجاوزاً التكاليف الحسابية الباهظة التي كانت مرتبطة سابقاً بالكيمياء الحركية التفصيلية في أكواد ديناميكا السوائل.