Abstract
A novel approximate solution for catalyst effectiveness factors is presented. It is based on carefully selected approximate reaction rate profiles, instead of typical assumption of composition profiles inside the catalyst. This formulation allows analytical solution of the approximate model, leading to a very simple iterative solution for effectiveness factor for general nonlinear reaction stoichiometry and arbitrary catalyst particle shape. The same model can be used with all practical Thiele modulus values, including multicomponent systems with inert compounds. Furthermore, the correct formulation of the underlying physical model equation is discussed. It is shown that an incorrect but often-used model formulation where convective mass transfer has been neglected may lead to much higher errors than the present approximation. Even with a correctly formulated physical model, rigorous discretization of the catalyst particle volume may have unexpectedly high numerical errors, even exceeding those with the present approximate solution. The proposed approximate solution was tested with a number of examples. The first was an equimolar reaction with first order kinetics, for which analytical solutions are available for the standard catalyst particle geometries (slab, long cylinder, and sphere). Then, the method was tested with a second order reaction in three cases: 1) with one pure reactant, 2) with inert present, and 3) with two reactants and non-stoichiometric surface concentrations. Finally, the method was tested with an industrially relevant catalytic toluene hydrogenation including Maxwell-Stefan formulation for the diffusion fluxes. In all the tested systems, the results were practically identical when compared to the analytical solutions or rigorous finite volume solution of the same problem.
Highlights
Calculation of effectiveness factors for catalyst particles is one of the most classical problems in Chemical Reaction Engineering [1,2,3,4]
The third test case was with elementary reaction for two reactants and non-stoichiometric surface compositions, which led to the depletion of one of the components and gradual change from second to pseudo-first order reaction as the reactants diffuse towards the catalyst particle center
An approximate method for solving reaction-diffusion-convection problems within catalyst particles was developed in this work
Summary
Calculation of effectiveness factors for catalyst particles is one of the most classical problems in Chemical Reaction Engineering [1,2,3,4]. Most of the numerical approximations found from the literature were based on physical models neglecting convective mass transfer [10,11,12,13,14,15,16,17] This leads to a situation where the resulting error is due to mathematical approximation of the underlying differential equation (which is typically minimized by elegant mathematical manipulations), but in the reaction-diffusion model itself, which does not take into account all the relevant physical phenomena. Many reactions catalyzed by a solid surface follow reaction kinetics without a clear reaction order, such as the Langmuir-Hinshelwood model In this contribution, a novel method for effectiveness factor prediction is proposed. The proposed method is suitable for any reaction rate regime, catalyst shape, and reaction stoichiometry
Published Version (Free)
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.