DOI: https://doi.org/10.3389/fenrg.2025.1608831
تاريخ النشر: 2026-01-04
المؤلف: Mario Zuber وآخرون
الموضوع الرئيسي: المحفزات لإصلاح الميثان
نظرة عامة
تبحث هذه الدراسة في إنتاج وقود الطيران المستدام عبر مسار الشمس إلى السائل، وتحديداً من خلال التوليد الحراري الكيميائي للغاز الاصطناعي (خليط من H₂ و CO) باستخدام الطاقة الشمسية المركزة. تركز الدراسة على طريقة تُعرف بإعادة تشكيل الأكسدة الجافة، والتي تستخدم عملية دورية من خطوتين باستخدام السيريا غير المتكافئة (CeO₂). تتضمن العملية تقليل السيريا بالميثان الحراري لإنتاج الغاز الاصطناعي، تليها أكسدة السيريا المخفضة باستخدام CO₂ لإنتاج CO. يتم تقديم نمذجة حرارية ديناميكية ومحاكاة باستخدام مفاعل أنبوبي ذو سرير ثابت، مما يظهر أن ارتفاع غير المتكافئات يعزز انتقائية الغاز الاصطناعي. النموذج الديناميكي للسوائل (CFD) الذي تم تطويره في OpenFOAM يلتقط بدقة ديناميات التفاعلات غير المتجانسة العكسية وتدفق السوائل، وتم التحقق منه ضد البيانات التجريبية مع أخطاء نسبية ضمن 7%.
تشير النتائج إلى ظروف مفاعل مثالية لإنتاج الغاز الاصطناعي، موصية بمعدل تدفق قدره 10 لترات·دقيقة⁻¹ وأقطار أنابيب تتراوح بين 7-10 سم عند حوالي 990 درجة مئوية. تم تقديم استراتيجية تشغيل جديدة، تُسمى التهيئة والدورة، لاستغلال تدرجات غير المتكافئات، مما يحسن بشكل كبير انتقائية الغاز الاصطناعي ومعدلات التحويل. تسلط الدراسة أيضًا الضوء على الإمكانية للتشغيل المستمر من خلال التغذية المشتركة لـ CH₄ و CO₂، مما قد يقلل من الطلبات الطاقية ويعزز الكفاءة. بشكل عام، تسهم هذه العمل في تقدم إنتاج وقود الطيران المستدام من خلال تحسين عملية إعادة تشكيل الأكسدة الجافة وتوفير نموذج موثق للبحوث المستقبلية وجهود التوسع.
مقدمة
تناقش مقدمة ورقة البحث النهج “الشمس إلى السائل” لإنتاج وقود الطيران المستدام (SAF) من خلال العمليات الحرارية الكيميائية الشمسية، مع التركيز على استخدام الطاقة الشمسية المركزة لدفع التفاعلات الماصة للحرارة لتخليق الغاز الاصطناعي. يمكن توليد الغاز الاصطناعي، وهو خليط من الهيدروجين (H₂) وأول أكسيد الكربون (CO)، عبر إعادة تشكيل الميثان المدفوعة بالطاقة الشمسية، والتي تشمل كل من إعادة التشكيل الجاف ($\text{CH}_4 + \text{CO}_2 \rightleftharpoons 2\text{CO} + 2\text{H}_2$) وإعادة التشكيل الرطب ($\text{CH}_4 + \text{H}_2\text{O} \rightleftharpoons \text{CO} + 3\text{H}_2$). تسلط الورقة الضوء على أهمية الحصول على الميثان بيولوجياً لضمان الاستدامة وتناقش التحديات المرتبطة بتعطيل المحفز وتكوين الكربون، والتي يمكن أن تعيق كفاءة الطاقة.
بالإضافة إلى ذلك، توضح المقدمة دورة أكسدة حرارية كيميائية من خطوتين باستخدام السيريا (CeO₂) كحامل أكسجين فعال لتفكيك الماء وثاني أكسيد الكربون. لقد أظهرت دورة الأكسدة بالسيريا تقدمًا كبيرًا، بما في ذلك العروض الناجحة على الشمس لسلسلة الشمس إلى الوقود. تقدم الورقة الديناميكا الحرارية لدورة الأكسدة ومعايير تشغيلها، مثل ظروف درجة الحرارة والضغط، مع معالجة التحديات الكامنة مثل فقدان الطاقة بسبب تقلبات درجة الحرارة. تهدف الدراسة إلى تطوير نموذج ديناميكي للسوائل (CFD) لمحاكاة عملية إعادة تشكيل الأكسدة الجافة في مفاعل ذو سرير ثابت، مع دمج المبادئ الديناميكية الحرارية والتحقق التجريبي لتعزيز فهم هذه التكنولوجيا الواعدة لإنتاج الوقود المستدام.
طرق
في هذا القسم، يوضح المؤلفون التحقق التجريبي من نموذج ديناميكي حراري للتفاعلات العكسية، باستخدام مفاعل أنبوبي على نطاق المختبر في ETH زيورخ. يعمل المفاعل، المصمم للتكرار الكيميائي المعاكس، تحت ضغط جوي مع سرير ثابت من كريات السيريا، محققًا نقطة ضبط درجة حرارة حوالي 1,000 درجة مئوية ومعدل تدفق قياسي قدره 1 لتر·دقيقة⁻¹. تُظهر كريات السيريا، التي تم تصنيعها باستخدام تقنية صب رطب سريعة، معدلات تحويل ميثان محسنة—أعلى بمقدار 1.8 مرة من RPCs التقليدية—مع الحفاظ على السلامة الهيكلية طوال عملية الاختبار. يتضمن الإعداد التجريبي التحكم الدقيق في تدفق الغاز ومراقبة درجة الحرارة، مما يضمن جمع بيانات دقيقة لتكوين المنتج ومؤشرات الأداء.
تشمل الإجراءات التجريبية سلسلة من دورات الأكسدة، مع تحديد ظروف معينة لخطوات الاختزال والأكسدة. تشير النتائج إلى تحويل كامل للميثان أثناء الاختزال، تليها إنتاج الغاز الاصطناعي، مع تغييرات ملحوظة في غير المتكافئات ودرجة الحرارة طوال الدورات. تكشف النتائج عن تدرج δ داخل السرير الثابت، مما يؤثر على إنتاج CO و CO₂ أثناء الأكسدة. تدعم البيانات المجمعة الحاجة إلى نموذج حسابي لالتقاط التأثيرات المعقدة متعددة الأبعاد التي لوحظت في عملية إعادة تشكيل الأكسدة الجافة، مما يحقق صحة قدرات النموذج التنبؤية مقابل النتائج التجريبية.
نتائج
تم تقييم نتائج النموذج مقابل البيانات التجريبية لإعداد درجة حرارة قدرها \( T_{sp} = 966 \, °C \). تشير النتائج الرئيسية إلى أن النموذج يعيد إنتاج الاتجاهات التي لوحظت في دورة الأكسدة خلال كل من خطوات الاختزال والأكسدة. على وجه التحديد، خلال خطوة الاختزال، يتماشى النموذج مع البيانات التجريبية بشأن التركيزات المولية لغازات المنتج على مر الزمن، ملتقطًا نقاط الانعطاف المميزة مع زيادة تركيز الغاز الاصطناعي وانخفاض تركيز الميثان. في خطوة الأكسدة، بينما يعكس النموذج الاتجاهات العامة، يتنبأ بنقطة الانعطاف لتركيز ثاني أكسيد الكربون لتحدث قبل 149 ثانية من الملاحظات التجريبية، ويعزى ذلك إلى الاختلافات الأولية في تركيز ثاني أكسيد الكربون.
تم تحليل مقاييس الأداء، بما في ذلك تحويل المواد الأولية (\( X_e \)) وانتقائية المنتج (\( S_p \)) عبر إعدادات درجات حرارة متعددة، مما يكشف أن كلا المقياسين يزيدان مع درجة الحرارة بسبب الطبيعة المواتية للتفاعلات التي تنتج ثاني أكسيد الكربون والماء. ومع ذلك، يميل النموذج إلى المبالغة في تقدير قيم التحويل الفورية وتقليل توقيت نقطة الانعطاف للأكسدة، مما يشير إلى وجود تباين في غير المتكافئات المتوقعة عند مخرج السرير الثابت. بشكل عام، يظهر النموذج تقريبًا معقولًا للاتجاهات خلال كل من مراحل الاختزال والأكسدة، مؤكدًا ملاءمة المعلمات الديناميكية الحرارية والحركية المختارة. تم تقديم تفاصيل التحقق الإضافية، بما في ذلك ملفات تعريف درجة الحرارة والثبات العددي، في زوبر (2024).
مناقشة
يركز النموذج الحسابي المقدم في هذا البحث على محاكاة السلوك العابر والديناميات الفيزيائية للتفاعلات غير المتجانسة العكسية التي تشمل أكاسيد صلبة غير متكافئة في تكوين سرير ثابت. يتم تنفيذ النموذج في بيئة OpenFOAM من خلال الحل porousRedoxFoam، ويستخدم الديناميكا الحرارية التوازنية ويشمل نموذج حركي لأخذ المنتجات الوسيطة في الاعتبار. إنه يقوم بنمذجة المجالات المسامية وتدفقات السوائل في تكوينات ثنائية وثلاثية الأبعاد، مما يسمح بتقييمات مفصلة للتدرجات وأبعاد المفاعل. تشمل المعادلات الحاكمة، المشتقة من porousGasificationFoam، معادلات الحفظ للطور الصلب والسائل، مع معالجة الاستمرارية والكتلة والطاقة وحفظ الزخم، بالإضافة إلى حفظ الحجم لضمان توازن المسامية أثناء التفاعلات.
جانب مهم من النموذج هو معالجته للتفاعلات العكسية، التي تتميز بثابت توازن يتغير مع غير المتكافئات ($\delta$) ودرجة الحرارة ($T_s$). يلتقط النموذج السلوك الديناميكي الحراري للتفاعلات من خلال نظام من المعادلات التفاضلية العادية (ODEs) الممثلة في مصفوفة جاكوب. يؤكد التحقق ضد الآليات الكيميائية المعروفة، مثل Cantera، دقة النموذج في التنبؤ بتراكيب المنتجات تحت ظروف متغيرة من غير المتكافئات ودرجة الحرارة. يمتد تطبيق النموذج إلى تحسين تكوينات المفاعل لإعادة تشكيل الأكسدة الجافة، مما يظهر فائدته كأداة تصميم لتعزيز إنتاج الغاز الاصطناعي من خلال التعديلات الاستراتيجية في أبعاد المفاعل وظروف التشغيل.
القيود
تناقش قسم القيود الرؤى والقيود المتعلقة بالنموذج المستخدم لتحليل غير المتكافئات في مفاعل ذو سرير ثابت خلال عمليات الاختزال والأكسدة. يبرز أنه بينما يلتقط النموذج بفعالية الاتجاهات لتفاعلات الأكسدة، فإنه يعتمد على مخطط تفاعل مبسط قد يتجاهل عوامل حاسمة مثل فراغات الأكسجين وخطوات التفاعل الوسيطة. لم يتم التحقق تجريبيًا من تدرج δ، الذي يعد حيويًا لتحديد التركيزات المولية للمنتجات، بسبب التعقيدات المرتبطة بأخذ العينات وإعادة أكسدة السرير الثابت تحت ظروف مسيطر عليها.
بالإضافة إلى ذلك، تطلب اعتماد النموذج على الحركيات الظاهرة والتحديات المتمثلة في تطبيق بيانات التحليل الحراري (TG) مباشرة على ظروف مفاعل السرير الثابت ضبط معلمات الحركية. تشير النتائج إلى أن درجات الحرارة المرتفعة تعزز حركيات التفاعل وتحويلات CH₄ لكنها تؤدي إلى تدرجات δ أكثر حدة وغير متجانسة، بينما تؤدي درجات الحرارة المنخفضة إلى تدرجات أكثر تجانسًا على حساب تقليل التحويلات. تم تحديد نمط نقل الحرارة الرئيسي على أنه التوصيل، مع آثار كبيرة على دقة النموذج تحت ظروف متغيرة، خاصة فيما يتعلق بتقلبات الضغط وأوامر التفاعل، التي تم التحقق منها فقط تحت ضغوط جزئية منخفضة للمواد الأولية.
DOI: https://doi.org/10.3389/fenrg.2025.1608831
Publication Date: 2026-01-04
Author(s): Mario Zuber et al.
Primary Topic: Catalysts for Methane Reforming
Overview
This research investigates the production of sustainable aviation fuels via the sun-to-liquid pathway, specifically through the thermochemical generation of syngas (a mixture of H₂ and CO) using concentrated solar energy. The study focuses on a method known as dry redox reforming, which employs a two-step redox cyclic process utilizing non-stoichiometric ceria (CeO₂). The process involves the methanothermal reduction of ceria to produce syngas, followed by the oxidation of the reduced ceria with CO₂ to yield CO. Thermodynamic modeling and simulations using a fixed-bed tubular reactor are presented, demonstrating that higher non-stoichiometries enhance syngas selectivity. The computational fluid dynamics (CFD) model developed in OpenFOAM accurately captures the dynamics of reversible heterogeneous reactions and fluid flow, validated against experimental data with relative errors within 7%.
The findings suggest optimal reactor conditions for syngas production, recommending a flow rate of 10 L·min⁻¹ and tube diameters of 7-10 cm at approximately 990 °C. A novel operating strategy, termed priming and cycling, is introduced to exploit non-stoichiometry gradients, significantly improving syngas selectivity and conversion rates. The study also highlights the potential for continuous operation through co-feeding CH₄ and CO₂, which could lower energy demands and enhance efficiency. Overall, this work contributes to advancing sustainable aviation fuel production by optimizing the dry redox reforming process and providing a validated model for future research and scale-up efforts.
Introduction
The introduction of the research paper discusses the “sun-to-liquid” approach for producing sustainable aviation fuels (SAF) through solar thermochemical processes, emphasizing the use of concentrated solar energy to drive endothermic reactions for syngas synthesis. Syngas, a mixture of hydrogen (H₂) and carbon monoxide (CO), can be generated via solar-driven methane reforming, which includes both dry reforming ($\text{CH}_4 + \text{CO}_2 \rightleftharpoons 2\text{CO} + 2\text{H}_2$) and wet reforming ($\text{CH}_4 + \text{H}_2\text{O} \rightleftharpoons \text{CO} + 3\text{H}_2$). The paper highlights the importance of sourcing methane biogenically to ensure sustainability and discusses the challenges associated with catalyst deactivation and carbon formation, which can hinder energy efficiency.
Additionally, the introduction outlines a two-step thermochemical redox cycle using ceria (CeO₂) as an effective oxygen carrier for splitting water and carbon dioxide. The ceria redox cycle has demonstrated significant advancements, including successful on-sun demonstrations of the solar-to-fuel chain. The paper presents the thermodynamics of the redox cycle and its operational parameters, such as temperature and pressure conditions, while addressing inherent challenges like energy losses from temperature swings. The study aims to develop a computational fluid dynamics (CFD) model to simulate the dry redox reforming process in a fixed-bed reactor, integrating thermodynamic principles and experimental validation to enhance the understanding of this promising technology for sustainable fuel production.
Methods
In this section, the authors detail the experimental validation of a thermodynamic model for reversible reactions, utilizing a lab-scale tubular reactor at ETH Zürich. The reactor, designed for countercurrent chemical looping, operates under ambient pressure with a fixed-bed of ceria pellets, achieving a temperature set point around 1,000 °C and a standard flow rate of 1 L·min⁻¹. The ceria pellets, fabricated using a rapid wet casting technique, demonstrate enhanced methane conversion rates—1.8 times higher than traditional RPCs—while maintaining structural integrity throughout the testing process. The experimental setup includes precise gas flow control and temperature monitoring, ensuring accurate data collection for product composition and performance indicators.
The experimental procedure involves a series of redox cycles, with specific conditions set for reduction and oxidation steps. The results indicate a complete conversion of methane during reduction, followed by syngas production, with notable changes in non-stoichiometry and temperature throughout the cycles. The findings reveal a δ-gradient within the fixed-bed, influencing the production of CO and CO₂ during oxidation. The data collected supports the need for a computational model to capture the complex multi-dimensional effects observed in the dry redox reforming process, thereby validating the model’s predictive capabilities against experimental outcomes.
Results
The results of the model were evaluated against experimental data for a temperature setting of \( T_{sp} = 966 \, °C \). Key findings indicate that the model effectively replicates the trends observed in the redox cycle during both the reduction and oxidation steps. Specifically, during the reduction step, the model aligns with experimental data regarding the molar concentrations of product gases over time, capturing the characteristic inflection points as syngas concentration increases and methane concentration decreases. In the oxidation step, while the model reflects the overall trends, it predicts the inflection point of carbon dioxide concentration to occur 149 seconds earlier than observed experimentally, attributed to initial differences in carbon dioxide concentration.
Performance metrics, including educt conversion (\( X_e \)) and product selectivity (\( S_p \)), were analyzed across multiple temperature settings, revealing that both metrics increase with temperature due to the favorable nature of reactions producing carbon dioxide and water. However, the model tends to overestimate instantaneous conversion values and underestimate the timing of the oxidation inflection point, suggesting a discrepancy in the predicted non-stoichiometry at the outlet of the fixed-bed. Overall, the model demonstrates a reasonable approximation of the trends during both reduction and oxidation phases, affirming the appropriateness of the selected thermodynamic and kinetic parameters. Further validation details, including temperature profiles and numerical stability, are provided in Zuber (2024).
Discussion
The computational model presented in this research focuses on simulating the transient behavior and physical dynamics of reversible heterogeneous reactions involving non-stoichiometric solid oxides in a fixed-bed configuration. The model, implemented in the OpenFOAM environment through the solver porousRedoxFoam, utilizes equilibrium thermodynamics and incorporates a kinetic model to account for intermediate products. It effectively models porous domains and fluid flows in both 2-D and 3-D configurations, allowing for detailed assessments of gradients and reactor dimensions. The governing equations, derived from porousGasificationFoam, include conservation equations for solid and fluid phases, addressing continuity, mass, energy, and momentum conservation, as well as volume conservation to ensure porosity balance during reactions.
A significant aspect of the model is its treatment of reversible reactions, characterized by an equilibrium constant that varies with non-stoichiometry ($\delta$) and temperature ($T_s$). The model captures the thermodynamic behavior of the reactions through a system of ordinary differential equations (ODEs) represented in a Jacobian matrix. Validation against established chemical mechanisms, such as Cantera, confirms the accuracy of the model in predicting product compositions under varying conditions of non-stoichiometry and temperature. The model’s application extends to optimizing reactor configurations for dry redox reforming, demonstrating its utility as a design tool for enhancing syngas production through strategic adjustments in reactor dimensions and operating conditions.
Limitations
The section on limitations discusses the insights and constraints of the model used to analyze non-stoichiometry in a fixed-bed reactor during reduction and oxidation processes. It highlights that while the model effectively captures the trends of the redox reactions, it is based on a simplified reaction scheme that may overlook critical factors such as oxygen vacancies and intermediate reaction steps. The δ-gradient, which is vital for determining product molar concentrations, was not experimentally validated due to the complexities involved in sampling and reoxidizing the fixed-bed under controlled conditions.
Additionally, the model’s reliance on apparent kinetics and the challenges of applying thermogravimetric (TG) data directly to fixed-bed reactor conditions necessitated tuning of kinetic parameters. The findings indicate that higher temperatures enhance reaction kinetics and CH₄ conversions but lead to steeper and non-uniform δ-profiles, while lower temperatures yield more uniform profiles at the cost of reduced conversions. The primary mode of heat transfer identified is conduction, with significant implications for the model’s accuracy under varying conditions, particularly regarding pressure variations and reaction orders, which were only validated under low educt partial pressures.