Modeling and evaluation of hydrodesulfurization and deactivation rates for partially wetted Trilobe catalyst using finite element method

• Published in: Powder technology 2019-09, Vol.354, p.779-791
• Main Authors: Limtrakul, Sunun, Bannatham, Papop, Teeraboonchaikul, Sornsawan, Vatanatham, Terdthai, Ramachandran, Palghat A.
• Format: Article
• Language: English
• Subjects: Catalyst deactivation; Catalysts; Computer applications; Deactivation; Diffusion; Diffusivity; Dimensional analysis; Finite element analysis; Finite element method; Hydrodesulfurization; Mass transfer; Mathematical analysis; Nonlinear programming; oils; Optimization; Partial wetting catalyst; Plugging; Pore size; Porosity; Reactors; Spatial variations; Trickle beds; Trilobe catalyst; Two dimensional analysis; Wettability
• ISSN: 0032-5910; 1873-328X
• DOI: 10.1016/j.powtec.2019.06.047

Hydrodesulfurization of heavy oil is usually operated in a trickle bed reactor. Analysis and design of this reactor require modeling both at the macro-reactor scale and micro-catalyst scale. The paper is focused on the catalyst-level modeling and takes into consideration complex catalyst shapes, wettability, and spatial variation in diffusivity due to catalyst pore plugging. Furthermore, the effects of mass transfer from dry and wet zones around the partially wetted catalyst surface are included. Finite element method is used for two-dimensional analysis as it is versatile and can handle complex shapes. The catalyst-level model includes an updating of effective diffusivity with time, and hence the lifetime of a complex shaped catalyst can be predicted. Long-time operation shows that thickness of plugged pore is large and also non-uniform across the two-dimensional catalyst particle with more plugging near the catalyst surface. In addition, the reaction regime changes with time, the fresh catalyst initially being in a diffusion-resistance-free regime and changing to a strong-diffusion-resistance regime at longer operating times. The computational method can be used for optimization of the operating protocol and to design optimum pore size and shape of the catalyst, although these details are not addressed in this paper. [Display omitted] •Catalyst-level models including wetting effect for a Trilobe shape are developed.•Effect of mass transfer in dry and wet zones around catalyst surface is included.•Combined reactor-level and catalyst-level models show long-term catalyst behavior.•Finite element method handles complex catalyst shape with discontinuous boundary.•Catalyst lifetime is predicted via pore plugging effect with updating diffusivity.