DOI: https://doi.org/10.1016/j.cej.2025.164593
تاريخ النشر: 2025-06-06
المؤلف: Hilbert Keestra وآخرون
الموضوع الرئيسي: تركيب الأمونيا وتقليل النيتروجين
نظرة عامة
تستعرض هذه القسم التحسينات التي تم إدخالها على مقال بحثي قبل نشره النهائي، مع التركيز على تطوير طريقة جديدة لتناسب عامل تصحيح القطبية في معادلة الحالة المعدلة Soave-Redlich-Kwong. إن عامل تصحيح القطبية -0.2368 يأخذ في الاعتبار بفعالية تأثيرات قطبية الأمونيا (NH₃) عند الضغوط العالية، مما يمكّن من وصف التوازن بدقة تصل إلى 1000 بار. يستخدم البحث نموذجًا من نوع Arrhenius لتوازن الأمونيا، والذي يتوافق مع برامج تخطيط التدفقات، ويستخدم خوارزمية استراتيجية تطور مصفوفة التغاير لضمان التقارب الموثوق في النماذج المعقدة.
بالإضافة إلى ذلك، يقدم الدراسة نموذجًا حركيًا من نوع Langmuir-Hinshelwood يصف بنجاح جميع مجموعات البيانات الحركية باستخدام تعبير موحد لمعدل التفاعل وعامل نشاط تحفيزي نسبي واحد لكل محفز. تسهل هذه التبسيط تقييمات النشاط التحفيزي وتدعم فرضية وجود آلية تفاعل مشابهة عبر مختلف المحفزات القائمة على الحديد، بغض النظر عن الاختلافات في المواد الأولية أو طرق التحضير. يظهر النموذج الحركي فعاليته عبر مجموعة واسعة من ظروف التشغيل، ويعزز دمج أكثر من 100 عام من بيانات تخليق الأمونيا بشكل كبير من قوة كل من النماذج الديناميكية الحرارية والحركية، كما هو مطبق في Aspen Plus.
مقدمة
تسلط المقدمة الضوء على أهمية الأمونيا (NH₃)، التي كان الطلب على إنتاجها حوالي 185 مليون طن في عام 2020. تُستخدم بشكل أساسي في الأسمدة وكمواد خام كيميائية، كما تكتسب الأمونيا أيضًا اهتمامًا كوقود خالٍ من الكربون وناقل للهيدروجين. تشمل تطبيقاتها المحتملة الاستخدام المباشر كوقود في المحركات البحرية وتوليد الطاقة الثابتة في محطات الطاقة الحرارية، والتوربينات الغازية، وخلايا الوقود. بالإضافة إلى ذلك، يمكن تكسير الأمونيا حراريًا إلى هيدروجين ونيتروجين، مما يسمح باستخدام الهيدروجين لاحقًا في قطاعات الطاقة والصناعة.
يتم تخليق الأمونيا من خلال عملية هابر-بوش، التي تجمع بين الهيدروجين (H₂) والنيتروجين (N₂) كما هو موضح في المعادلة ذات الصلة. تعتبر هذه الطريقة في التخليق ضرورية لتلبية الطلب المتزايد على الأمونيا في تطبيقات متنوعة، بما في ذلك دورها في الانتقال إلى حلول الطاقة المستدامة.
النتائج
تقدم نتائج الدراسة تحليلًا شاملاً لنموذج Langmuir-Hinshelwood الذي تم ملاءمته عبر 11 مجموعة بيانات، مع تسليط الضوء على الخصائص الإحصائية لنماذج مختلفة تأخذ في الاعتبار أنواعًا مختلفة من المواد الممتصة. تشير النتائج إلى أن النماذج التي تتضمن فقط النيتروجين ($N^*$) والهيدروجين ($H^*$) كأنواع سطحية تحقق أفضل أداء تنبؤي، مع متوسط خطأ تربيعي (MSE) قدره $4.03 \times 10^{-5}$. بالمقابل، فإن تضمين أنواع إضافية مثل الأمونيا ($NH^*$) لا يعزز دقة النموذج بشكل كبير، مما يشير إلى تغطيتها الضئيلة تحت الظروف المختبرة. كما تؤكد الدراسة على أهمية أخذ عدم المثالية في التفاعلات الغازية في الاعتبار، كما يتضح من مقارنة الفوغاسيتي والضغط الجزئي في النماذج.
يظهر النموذج الحركي المستمد من الدراسة حاجز تنشيط واضح قدره 100.7 كيلوجول·مول$^{-1}$، متوافقًا مع القيم الموجودة في الأدبيات، بينما تم الإبلاغ عن طاقات الربط لـ $N^*$ و $H^*$ كـ -17.95 كيلوجول·مول$^{-1}$ و -4.948 كيلوجول·مول$^{-1}$، على التوالي. يتم دعم قوة النموذج بشكل أكبر من خلال وجود ارتباط قوي بين كسور مول الأمونيا التجريبية والمتوقعة، كما هو موضح في رسم التوزيع. بالإضافة إلى ذلك، يوفر معامل النشاط التحفيزي النسبي ($b_{cat}$) وسيلة لمقارنة أداء مختلف المحفزات القائمة على الحديد، مما يعكس التباينات في كثافة المواقع النشطة وتأثيرات المحفزات دون تغيير النموذج الحركي الأساسي. بشكل عام، تؤكد الدراسة أن النموذج الحركي المقترح كافٍ لنطاق درجات الحرارة من 251-550 °م، بينما تعالج أيضًا القيود والافتراضات المتأصلة في نهج النمذجة.
المناقشة
تستعرض قسم المناقشة في الورقة التقدمات الكبيرة في نمذجة تخليق الأمونيا باستخدام المحفزات القائمة على الحديد. يقترح المؤلفون نموذجًا حركيًا جديدًا يدمج إطار عمل Langmuir-Hinshelwood مع عامل نشاط تحفيزي نسبي واحد ($b_{cat}$)، مما يبسط تقييم أداء المحفز عبر ظروف متنوعة. يبني هذا النموذج على العمل الأساسي لـ Temkin، الذي حدد تفكك النيتروجين كخطوة تحدد المعدل وقدم مفهوم الأسطح التحفيزية غير المتجانسة. من خلال توسيع النموذج لاستيعاب مجموعات بيانات متعددة وأنواع محفزات متنوعة، يظهر المؤلفون دقة محسّنة في التنبؤ بحركيات تخليق الأمونيا عبر نطاق واسع من درجات الحرارة (251-550 °م) والضغوط (1-324 بار).
بالإضافة إلى ذلك، تناقش الورقة تطوير معادلة الحالة المعدلة Soave-Redlich-Kwong (SRK EOS) التي تتضمن عامل تصحيح القطبية للأمونيا، مما يعزز قدرة النموذج على وصف سلوك الغاز غير المثالي عند الضغوط العالية. تم تصميم تعبير التوازن من نوع Arrhenius لتسهيل التنفيذ في برامج تخطيط التدفقات التجارية، مما يسهل التطبيقات العملية في البيئات الصناعية. تشير النتائج إلى أن آليات التفاعل لمختلف المحفزات القائمة على الحديد متشابهة من الناحية الأساسية، على الرغم من التباينات في تركيبها وتحضيرها، مما يدعم فائدة النموذج الحركي المبسط لاختيار المحفزات وتقييمها.
DOI: https://doi.org/10.1016/j.cej.2025.164593
Publication Date: 2025-06-06
Author(s): Hilbert Keestra et al.
Primary Topic: Ammonia Synthesis and Nitrogen Reduction
Overview
This section outlines the enhancements made to a research article prior to its final publication, emphasizing the development of a novel fitting method for the polarity correction factor in the modified Soave-Redlich-Kwong Equation-of-State. The polarity correction factor of -0.2368 effectively accounts for ammonia (NH₃) polarity effects at high pressures, enabling accurate equilibrium descriptions up to 1000 bar. The research employs an Arrhenius-type model for ammonia equilibrium, which is compatible with flowsheeting software, and utilizes the Covariance Matrix Adaptation Evolution Strategy algorithm to ensure reliable convergence in complex models.
Additionally, the study introduces a Langmuir-Hinshelwood kinetic model that successfully describes all kinetic data sets using a unified rate expression and a single Relative Catalytic Activity factor per catalyst. This simplification facilitates catalytic activity evaluations and supports the hypothesis of a similar reaction mechanism across various iron-based catalysts, regardless of differences in precursor or preparation methods. The kinetic model demonstrates its efficacy across a broad range of operational conditions, and the integration of over 100 years of ammonia synthesis data significantly enhances the robustness of both the thermodynamic and kinetic models, as implemented in Aspen Plus.
Introduction
The introduction highlights the significance of ammonia (NH₃), which had a production demand of approximately 185 million tons in 2020. Primarily utilized in fertilizers and as a chemical feedstock, ammonia is also gaining attention as a zero-carbon fuel and hydrogen carrier. Its potential applications include direct use as fuel in maritime engines and stationary power generation in thermal power plants, gas turbines, and fuel cells. Additionally, ammonia can be thermally cracked into hydrogen and nitrogen, allowing for hydrogen’s subsequent use in energy and industrial sectors.
Ammonia is synthesized through the Haber-Bosch process, which combines hydrogen (H₂) and nitrogen (N₂) as described in the relevant equation. This synthesis method is crucial for meeting the growing demand for ammonia in various applications, including its role in the transition to sustainable energy solutions.
Results
The results of the study present a comprehensive analysis of the Langmuir-Hinshelwood model fitted across 11 datasets, highlighting the statistical properties of various models that account for different adsorbed species. The findings indicate that models incorporating only nitrogen ($N^*$) and hydrogen ($H^*$) as surface species yield the best predictive performance, with a mean squared error (MSE) of $4.03 \times 10^{-5}$. In contrast, including additional species such as ammonia ($NH^*$) does not significantly enhance model accuracy, suggesting their negligible coverage under the tested conditions. The study also emphasizes the importance of accounting for non-ideality in gas-phase interactions, as evidenced by the comparison of fugacity and partial pressure in the models.
The kinetic model derived from the study shows an apparent activation barrier of 100.7 kJ·mol$^{-1}$, consistent with literature values, while the binding energies for $N^*$ and $H^*$ are reported as -17.95 kJ·mol$^{-1}$ and -4.948 kJ·mol$^{-1}$, respectively. The model’s robustness is further supported by a strong correlation between experimental and predicted ammonia mole fractions, as illustrated in the parity plot. Additionally, the relative catalytic activity ($b_{cat}$) parameter provides a means to compare the performance of different iron-based catalysts, reflecting variations in active site density and the effects of promoters without altering the fundamental kinetic model. Overall, the study confirms that the proposed kinetic model is sufficient for the temperature range of 251-550 °C, while also addressing the limitations and assumptions inherent in the modeling approach.
Discussion
The discussion section of the paper outlines significant advancements in the modeling of ammonia synthesis using iron-based catalysts. The authors propose a novel kinetic model that integrates the Langmuir-Hinshelwood framework with a single Relative Catalytic Activity ($b_{cat}$) factor, simplifying the evaluation of catalyst performance across various conditions. This model builds on the foundational work of Temkin, which identified nitrogen dissociation as the rate-determining step and introduced the concept of non-uniform catalytic surfaces. By expanding the model to accommodate multiple datasets and varying catalyst types, the authors demonstrate improved accuracy in predicting ammonia synthesis kinetics over a wide range of temperatures (251-550 °C) and pressures (1-324 bar).
Additionally, the paper discusses the development of a modified Soave-Redlich-Kwong Equation-of-State (SRK EOS) that incorporates a polarity correction factor for ammonia, enhancing the model’s ability to describe non-ideal gas behavior at high pressures. The proposed Arrhenius-type equilibrium expression is designed for ease of implementation in commercial flowsheeting software, facilitating practical applications in industrial settings. The findings indicate that the reaction mechanisms for different iron-based catalysts are fundamentally similar, despite variations in their composition and preparation, thereby supporting the utility of the simplified kinetic model for catalyst selection and evaluation.