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

ملخص تقني: النمذجة الحركية لإعادة تشكيل الميثان/ثاني أكسيد الكربون جافاً في تفريغ الحاجز العازل

بيان المشكلة
يوفر تحويل الميثان (CH4) وثاني أكسيد الكربون (CO2) عبر إعادة التشكيل الجاف بمساعدة البلازما (DRM) في تفريغات الحاجز العازل (DBDs) مساراً غير حراري لإنتاج الهيدروجين والغاز الاصطناًعي (syngas). ومع ذلك، فإن أدوات النمذجة الحركية لهذه العملية غالباً ما تعتمد على منصات متخصصة وغير نمطية، وتستخدم تكراراً ثوابت معدل تقديرية (بترتيب المقدار) بدلاً من البيانات التي تم التحقق من صحتها بصرامة. علاوة على ذلك، هناك نقص في الأطر المفتوحة التي تربط صراحةً بين الكيمياء المتجانسة صفرية الأبعاد (0D) مع انتقال التشتت المحوري أحادي البعد (1D)، مما يجعل من الصعب تصميم بارامترات التشغيل بسرعة أو التحقق من صحة النماذج مقابل المعايير التجريبية بدقة كمية.

المنهجية
قام المؤلفون بتطوير إطار عمل مفتوح ونمطي بلغة بايثون (0D/1กัน 1D) لمحاكاة تحويل CH4/CO2 في تفريغ الحاجز العازل (DBD).

• الآلية الحركية: تم استخدام آلية مختزلة تتكون من 23 نوعاً متعادلاً و38 تفاعلاً. يتضمن ذلك التفكك الناتج عن اصطدام الإلكترونات، والتفاعلات ثنائية الجزيئات بين الجسيمات المتعادلة، والارتباطات ثلاثية الأجسام، والتفككات أحادية الجزيء.
• حركية الإلكترونات: بخلاف النسخ السابقة التي استخدمت ملاءمات تحليلية، تم اشتقاق معاملات معدل اصطدام الإلكترونات (E1–E5) من حل حقيقي لمعادلة بولتزمان للإلكترونات ثنائية الحدود (باستخدام برنامج الحل BOLOS) والمغذى بقطاعات عرض حقيقية لاصطدام الإلكترونات بالجسيمات المتعادلة من قاعدة بيانات LXCat/Morgan.
• التحفيز بالبلازما: تم نمذجة الـ DBD كقطار نبضات دوري يمثل التفريغات الشعيرية الدقيقة. تم تطبيق تقريب متوسط زمن الدورة لاستبدال المعدلات اللحظية بالقيم متوسطة الدورة، وهو ما تم التحقق منه لتكرار النتائج المعتمدة على النبضات بخطأ أقل من 1% مع تقليل التكلفة الحسابية بشكل كبير.
• تحليل الحساسية: تم إجراء تحليل حساسية "موريس" العالمي على 33 ثابت معدل للجسيمات المتعادلة. حدد هذا التحليل ستة تفاعلات محددة باعتبارها المحركات الرئيسية لتحويل CO2، ونسبة H2/CO، وكفاءة الطاقة.
• استراتيجية التحقق: تم استبدال التفاعلات الستة الحرجة المحددة (R2, R11–R14, R33) بقيم من أدبيات حركية الاحتراق والغلاف الجوي الموثقة. وتم الإبقاء على الـ 27 ثابت معدل الأخرى كمعاملات توضيحية بترتيب المقدار.
• هيكل النموذج: يربط إطار العمل نموذج "لاغرانج" صفر الأبعاد (يتتبع جزيء غاز) مع نموذج مفاعل تدفق أفقي ذو تشتت محوري "أويلر" أحادي البعد لمحاكاة الملامح المكانية.

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

• التحقق من الآلية: أدى استبدال التخمينات التوضيحية بقيم من الأدبيات إلى تصحيح أخطاء تراوحت بين 6 أضعاف إلى أكثر من 3000 ضعف. وتحديداً، تم تثبيت تفاعل CO + O (+M) → CO2 (+M) (R33)، وهو تفاعل محظور دورانياً ويصعب قياسه، بناءً على قيمة أدبية عند 200 كلفن (وهي أفضل قيمة متاحة، رغم أن النموذج يعمل عند 400 كلفن).
• اتجاهات التحويل: يحترم تحويل CO2 سقف الـ ~25% الملاحظ عبر 19 مصدراً للبلازما عند الضغط الجوي في الأدبيات لمدخلات طاقة نوعية (SEI) تصل إلى ~30,000 جول/لتر. ويرتفع تحويل CH4 بشكل أكثر حدة، مقترباً من التحويل شبه الكامل عند قيم SEI عالية.
• الانتقائية: تستقر نسبة H2/CO بين 1.1 و2.1، مما يعكس نظاماً محدوداً بـ CH4 حيث يتفكك CH4 بشكل أسرع من CO2، وهو ما يتوافق مع الملاحظات التجريبية.
• المقارنة المرجعية:
  • عند SEI مطابقة تبلغ 5.7 إلكترون فولت/جزيء، بالغ النموذج في البداية في تقدير تحويل CH4 بنحو 5.8 ضعف وتحويل CO2 بنحو 2.7 ضعف مقارنة بالمرجع التجريبي لـ Ozkan et al. [55].
  • أدى تطبيق SEI "مصحح القدرة" (الذي يأخذ في الاعتبار نقل الطاقة الإجمالي بدلاً من مجرد طاقة اصطدام الإلكترونات) إلى تقليل هذه الفجوات إلى 3.3 ضعف (CH4) و1.6 ضعف (CO2).
  • والأهم من ذلك، عند قيم SEI المنخفضة والنموذجية تجريبياً (938–2,500 جول/لتر)، أظهر النموذج توافقاً وثيقاً مع التجربة، حيث قلل من تقدير تحويل CH4 قليماً وطابق تحويل CO2 بدقة. وتزدزد الفجوة فقط عند SEI العالية، مما يشير إلى أن النموذج قوي في نطاق SEI المنخفض إلى المتوسط المعتاد في مختبرات الـ DBD.
• مصادر عدم اليقين: أكد تحليل الحساسية أنه حتى بعد التحقق، يظل R33 هو المصدر المهيمن لعدم اليقين في تحويل CO2 ونسبة H2/CO. ويُعزى التقدير المرتفع المتبقي لتحويل CH4 إلى كيمياء مجمع جذور CH4 غير المتحقق منها (R1, R3–R10) والافتراض بأن نسب التفرع لقنوات التفكك هي نسب موزعة للقنوات.

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

1. إطار عمل مفتوح: توفير كود برمجي بلغة بايثون قابل للتكرار يربط كيمياء 0D بنقل 1D، وهو متاح للنماذج السريعة والتدريس.
2. ترقية منهجية: الانتقال من الملاءمات التحليلية التعسفية إلى معاملات معدل مستمدة من حل بولتزمان الحقيقي بناءً على قطاعات عرض حقيقية.
3. التحقق المستهدف: إثبات أن استبدال ستة ثوابت معدل رئيسية فقط بقيم من الأدبيات بناءً على تحليل الحساسية يحسن بشكل كبير من دقة النموذج، مما يزيل مصدراً رئيسياً للخطأ ويقرب التنبؤات من الواقع التجريبي.
4. القيود الواقعية: يدعي المؤلفون بتواضع أنه بينما لا يزال النموذج غير دقيق كمياً عبر جميع نطاقات SEI، فإنه يوفر أساساً "نوعياً إلى شبه كمي" لنطاق SEi المنخفض إلى المتوسط. كما حددوا بوضوح خارطة الطريق المتبقية للتحقق الكامل: الحصول على بيانات درجات حرارة أقل لـ R33، والتحقق من كيمياء جذور CH4، ومعايرة بارامترات تحفيز البلازما مقابل القياسات الكهربائية المحددة للمفاعل.

لا يدعي العمل أنه قد حل مشكلة DRM كمياً، بل يقدم أساساً حسابياً شفافاً ومفتوحاً مع خارطة طريق محددة ومضيقة لتحقيق نموذج متحقق منه بالكامل.