Abstract
This study introduces a multiscale model for ammonia synthesis catalysts that integrates intrinsic kinetics, intraparticle diffusion, and deactivation mechanisms, specifically water vapor self-poisoning and long-term aging. Extending the Temkin-Pyzhev kinetic framework, the model incorporates a size-dependent self-poisoning coefficient (γ(dp)), a time-dependent aging factor (af(t)), and a Thiele modulus-based effectiveness factor (η). Calibrated with experimental data, it accurately predicts nitrogen consumption rates (rN2) for catalyst particle sizes ranging from 0.6 to 9.0 mm and operational lifetimes of 2 to 5 years, with errors as low as 0.8% for larger particles. Unlike traditional models, this approach quantifies reduction-induced deactivation, which significantly impacts larger particles by markedly reducing activity. Implemented in MATLAB, the model provides a predictive tool for optimizing catalyst design and reactor performance under industrial conditions. By linking microkinetic, transport, and deactivation phenomena, this work enhances the efficiency and longevity of ammonia synthesis processes, both traditional and novel.




Similar content being viewed by others
Data availability
The author declares that the data supporting the findings of this study are available within the paper and its supplementary information files. Should any raw data files be needed in another format they are available from the corresponding author upon reasonable request. Source data are provided with this paper.
Abbreviations
- CFD:
-
Computational fluid dynamic
- A :
-
Pre-exponential factor, h−1
- a i :
-
Activity, i = N2, H2, NH3
- af(t) :
-
Time-dependent aging factor, dimensionless
- C SM :
-
Stefan-Maxwell correction, dimensionless
- D e :
-
Effective diffusivity, m2/h
- dp :
-
Particle diameter, mm
- E a :
-
Activation energy, J/mol
- f i :
-
Fugacity, i = N2, H2, NH3
- f i * :
-
Fugacity at reference state, i = N2, H2, NH3
- k :
-
Rate constant, h−1
- K a :
-
Equilibrium constant, dimensionless
- P :
-
Pressure, atm
- R :
-
Universal gas constant, J/mol/K
- R p :
-
Pellet radius, m
- r NH3 :
-
Intrinsic rate of ammonia formation, mol/m3s
- r N2 :
-
Intrinsic rate of nitrogen consumption, mol/m3s
- t :
-
Time, years
- T :
-
Thermodynamic temperature, Kelvins
- x :
-
Fractional conversion of N2,
- α :
-
Constant
- ε:
-
Catalyst porosity, dimensionless
- γ(dp) :
-
Size-dependent self-poisoning coefficient, dimensionless
- ϕ :
-
Thiele modulus, dimensionless
- \(\Phi\) :
-
Fugacity coefficient, i = N2, H2, NH3
- η :
-
Effectiveness factor, dimensionless
- τ :
-
Tortuosity (labyrinth) factor, dimensionless
References
IEA (2021) Ammonia technology roadmap. IEA, Abu Dhabi
IRENA (2022) Innovation outlook renewable ammonia. IRENA, Abu Dhabi
Skubic L, Gyergyek S, Husans M, Likozar B (2024) A review of multiscale modelling approaches for understanding catalytic ammonia synthesis and decomposition. J Catal 429:115217. https://doi.org/10.1016/j.jcat.2023.115217
Ertl G (1991) Mechanism of ammonia synthesis. In: Jennings JR (ed) Catalytic ammonia synthesis: fundamentals and practice. Plenum Press, New York
Nielsen A, Kjaer J, Hansen B (1964) Rate equation and mechanism of ammonia synthesis at industrial conditions. J Catal 5(1):68–79. https://doi.org/10.1016/0021-9517(64)90094-6
Ponikvar Z, Likozar B, Gyergyek B (2022) Electrification of catalytic ammonia production and decomposition reactions: from resistance, induction, and dielectric reactor heating to electrolysis. ACS Appl Energy Mater 5:5457–5472
Belluci Sessa A, Abdoos B, Turco H, Di Benedetto A, Russo D (2025) Ammonia-based processes integrated with fuel cells for electrical power generation: a review. Renew Sustain Energy Rev 224:116134. https://doi.org/10.1016/j.rser.2025.116134
Rakic E, Grilc M, Likozar B (2023) Liquid organic hydrogen carrier hydrogenation–dehydrogenation: from ab initio catalysis to reaction micro-kinetics modelling. Chem Eng J 472:144836. https://doi.org/10.1016/j.cej.2023.144836
Pavlisic A, Hus M, Prasnikar A, Likozar B (2020) Multiscale modelling of CO2 reduction to methanol over industrial Cu/ZnO/Al2O3 heterogeneous catalyst: linking ab initio surface reaction kinetics with reactor fluid dynamics. J Clean Prod 275:122958. https://doi.org/10.1016/j.jclepro.2020.122958
Salciccioli M (2011) Multiscale modeling of heterogeneous catalysis: From quantum mechanics to reactor design. Chem Eng Sci 66(20):4319–4355
Temkin M, Pyzhev V (1940) Kinetics of ammonia synthesis on promoted iron catalyst. Acta Phys-Chim 12:327–356
Aika K, Christiansen LJ, Dybkjaer I, Hansen JB, Hφjlund Nielsen PE, Nielsen A, Stolze P, Tamaru K (1995) Ammonia catalysis and manufacture. Springer, Berlin Heidelberg. https://doi.org/10.1007/978-3-642-79197-0
Twigg MV (1989) Catalyst handbook. Wolfe Publishing Ltd., London
Nielsen A (1968) An Investigation on Promoted Iron Catalysts for the Synthesis of Ammonia. Doctoral Thesis, Haldor Topsøe A/S
Levenspiel O (1999) Chemical reaction engineering. Wiley, New York
Bartholomew CH (2001) Mechanism of catalyst deactivation. Appl Catal A 212(1–2):17–60
Argyle MD, Bartholomew CH (2015) Heterogeneous catalyst deactivation and regeneration: a review. Catalysts 5(1):145–169. https://doi.org/10.3390/catal5010145
Moghaddam AA, Krewer U (2020) Catalysts 10(11):1225. https://doi.org/10.3390/catal10111225
Jennings JR (1991) Catalytic ammonia synthesis: fundamentals and practice. Plenum Press, New York
Cholewa T, Steinbach B, Heim C, Nestler F, Nanba T, Güttel R, Salem O (2024) Reaction kinetics for ammonia synthesis using ruthenium and iron based catalysts under low temperature and pressure conditions. Sustain Energy Fuels 8:2245–2255. https://doi.org/10.1039/d4se00254g
Dyson DC, Simon JM (1989) Kinetic modeling of ammonia synthesis. Ind Eng Chem Fundam 8(3):362–370
Froment GF, Bischoff KB, De Wilde J (2010) Chemical reactor analysis and design. Wiley, New York
Raróg-Pilecka W (2003) Particle size effects in ammonia synthesis over ruthenium catalysts. J Catal 218(2):465–469
Liu H (2014) Ammonia synthesis catalyst 100 years: practice, enlightenment and challenge. Chin J Catal 35(10):1619–1640
Zybert M (2022) Structure sensitivity of ammonia synthesis on cobalt: effect of the cobalt particle size on the activity of promoted cobalt catalysts supported on carbon. Catalysts 12(10):1285. https://doi.org/10.3390/catal12101285
Peng X, Chen X, Zhou Y, Sun F, Zhang T, Zheng L, Jiang L, Wang X (2022) Size-dependent activity of supported Ru catalysts for ammonia synthesis at mild conditions. J Catal 408:98. https://doi.org/10.1016/j.jcat.2022.02.024
Gillespie LJ, Beattie JA (1930) The thermodynamic treatment of chemical equilibria in systems composed of real gases. I. An approximate equation for the mass action function applied to the existing data on the Haber equilibrium. Phys Rev 36:743. https://doi.org/10.1103/PhysRev.36.743
Funding
The author has no relevant financial or non-financial interests to disclose. The author certifies that he has no affiliations with or involvement in any organization or entity with any financial interest or non-financial interest in the subject matter or materials discussed in this manuscript. The author has no financial or proprietary interests in any material discussed in this article.
Author information
Authors and Affiliations
Contributions
Nenad Zecevic conceptualized the framework, performed the analysis, and wrote the manuscript.
Corresponding author
Ethics declarations
Conflict of interest
The authors have no competing interests to declare that are relevant to the content of this article.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Zecevic, N. Multiscale catalyst model for ammonia synthesis: coupling kinetics, diffusion and deactivation. Reac Kinet Mech Cat 138, 3645–3664 (2025). https://doi.org/10.1007/s11144-025-02947-4
Received:
Accepted:
Published:
Version of record:
Issue date:
DOI: https://doi.org/10.1007/s11144-025-02947-4



