# Analytical Determination of Rational Catalyst Lifetime for Consecutive Reaction in Isothermal Plug Flow Reactor

## DOI:

https://doi.org/10.20535/1810-0546.2018.1.105885## Keywords:

Mathematical modeling, Plug flow reactor, Consecutive irreversible reaction, Deactivation of solid catalyst, Catalyst lifetime## Abstract

**Background.** The urgency of mathematical modeling of continuous chemical-technological processes in non-stationary conditions caused by the actions of various destabilizing factors is undoubted. In this case, analytical solutions have undeniable advantages over numerical ones, since they make it possible to clarify the nature of the cause-effect relationships in the properties of the modeling object under consideration and, as a practical result, to give physically grounded recommendations for increasing the efficiency of its functioning. For catalytic processes, the reason for non-stationarity is the deactivation of the solid catalyst (Kt). This leads to a decrease of the conversion degree \[{x_ \bullet } = 1 - {c_{1 \bullet }}\] of the A_{1} reagent, and hence to a negative increase of its concentration \[{c_{1 \bullet }}\] in the reaction mixture at the outlet of the plug flow reactor (PFR) of length \[{L_{0 \bullet }},\] which leads to economic losses. Therefore, a rational (maximally expedient) catalyst lifetime \[{\theta _{\max }} > > 1\] has fundamental importance and is a significant part of the individual problem of selecting Kt.

**Objective****.** The aim of the paper is analytical calculation of the maximally expedient lifetime \[{\theta _{\max }} = {\tau _{\max }}/{\tau _{L \bullet }}\] of industrial Kt at the point \[{x_{0 \bullet }}\] of maximum of the nominal yield \[{\eta _{02 \bullet }} \equiv \eta _{02}^{\max }\] of the product A_{2} for the isothermal system “PFR at the optimum residence time \[{\tau _{L \bullet }} = {L_{0 \bullet }}/{u_0}\] of the reactants + consecutive catalytic reaction \[\mathrm{A}_{1}\xrightarrow[\mathrm{Kt},k_{\mathrm{d}1}]{k_{01},n_{1}=1}\alpha _{2}\mathrm{A}_{2}\xrightarrow[\mathrm{Kt},k_{\mathrm{d}2}]{k_{02},n_{2}=1}\alpha _{3}\mathrm{A}_{3}"\] under the influence of the Kt deactivation destabilizing factor.

**Methods****. **A linearized mathematical model in the form of a system of ordinary differential equations of characteristics for calculating the relatively small Kt deactivation influence on the system operating mode stationarity has been used.

**Results****. **In the case of the first-order reaction for the conditions of the industrial Kt deactivation, the relative deviations \[\left | \varepsilon _{x\bullet } \right |=\left | (x_{\bullet }/x_{0\bullet })-1 \right |\sim k_{\mathrm{d}1}\tau _{\mathrm{max}}< < 1\] of the degree of A_{1} conversion, the relative deviations \[|{\varepsilon _{\eta 2\bullet }}|\, \sim {k_{{\text{d1}}}}\tau_{max}\] of the A_{2} yield and the relative deviations \[\varepsilon _{s2\bullet }\sim {k_{{\text{d1}}}}\tau_{max}\] of selectivity \[{s_{2 \bullet }} = {\eta _{2 \bullet }}/{x_ \bullet }\] from the nominal values are analytically calculated in the linear approximation. It is established that the magnitudes \[{\varepsilon _{x \bullet }} < 0,\] \[|{\varepsilon _{\eta 2 \bullet }}|\geqslant 0\] and \[{\varepsilon _{s2 \bullet }} > 0\] are determined by the simplex \[\gamma _{0k}=k_{01}/k_{02}\] of the nominal rate constants and by the simplex \[\gamma _{\mathrm{d}}=k_{\mathrm{d}2}/k_{\mathrm{d}1}\] of the Kt deactivation rate constants of the reaction stages.

**Conclusions****. **It is proved that with respect to the yield of A_{2} product, there is a self-regulation effect \[({\varepsilon _{\eta 2 \bullet }} \approx 0)\] of the stationary mode at the condition of the equality \[\gamma _{\mathrm{d}}=1\] of the Kt deactivation rate constants. A nomogram for determining \[1 < <\theta _{max}< < (k_{\mathrm{d}1}\tau _{L\bullet })^{-1}\] on the maximum admissible value of \[\left | \varepsilon _{x\bullet } \right |_{max}^{\mathrm{adm}}< < 1.\] is calculated. For example, at a nominal degree of conversion \[{x_{0 \bullet }} = 75\% \Leftrightarrow {\gamma _{0k}} = 2\] of A_{1} reagent and \[\left | \varepsilon _{x\bullet } \right |_{max}^{\mathrm{adm}}= 5\%\Rightarrow \theta _{max}\approx 1,1\cdot 10^{3}(k_{\mathrm{d}1}\tau _{L\bullet }=10^{-4}).\] It is shown that the rational catalyst lifetime \[{\theta _{\max }}\] is inversely proportional to the complex \[k_{\mathrm{d}1}\tau _{L\bullet }\] of the Kt deactivation rate constant of the first stage.

## References

A.Yu. Zakgeym, *Introduction** **to** **Modeling** **of** **Chemical**-**Technological** **Processes*. Moscow, Russia: Khimiya, 1982 (in Russian).

K.R. Westerterp *et al.*, *Chemical Reactor Design and Operation*.* *New York: Wiley, 1987. doi: 10.1002/aic.690320124

Yu.M. Zhorov, *Kinetics** **of** **Industrial Organic Reactions*. Moscow, Russia: Khimiya, 1989 (in Russian).

M.E. Davis and R.J. Davis, *Fundamentals of Chemical Reaction Engineering*. Boston: McGraw-Hill, 2003.

D.V. Evdokimov *et al.*, “Analysis of development tendencies of modern mathematical and numerical modelling”, *Visnyk DNU,* *Ser. Modeliuvannia*, no. 8, pp. 3–17, 2009 (in Russian). doi: 10.15421/140901

U. Mann, *Principles of Chemical Reactor Analysis and Design*. Hoboken, New Jersey: John Wiley & Sons, 2009. doi: 10.1002/ 9780470385821

*Catalysis: From Principles to Applications. *M. Beller *et al.*, eds. New York: Wiley-VCH, 2012. doi: 10.1002/anie.201210089

J.B. Rawlings and J.G. Ekerdt, *Chemical Reactor Analysis and Design Fundamentals*. Madison, Wisconsin: Nob Hill Publishing, LLC, 2012.

A. Rasmuson* et** **al**.*, *Mathematical Modeling in Chemical Engineering*. Cambridge, UK: MPG Print-group Ltd., 2014. doi: 10.1002/cite.201590043

N.N. Kulov and L.S. Gordeev, “Mathematical modeling in chemical engineering and biotechnology”, *Theor**. **Found**. **Chem**. **Eng*., vol. 48, no. 3, pp. 243–248, 2014. doi: 10.1134/S0040579514030099

A.K. Verma, *Process **Modelling and Simulation in Chemical, Biochemical and Environmental Engineering*. Boca Raton: CRC Press, 2014. doi: 10.1002/cite.201690094

L.K Doraiswamy and D. Uner, *Chemical Reaction Engineering. Beyond the Fundamentals*. Boca Raton, London, New York: CRC Press, 2014.

S.R. Upreti,* **Process Modeling and Simulation for Chemical Engineers*: *Theory and Practice*. New York: Wiley, 2017. doi: 10.1002/9781118914670.ch8

I.D. Lucheyko, “Deactivation of catalyst in system “consecutive reaction A_{1 }® a_{2}A_{2 }® a_{3}A_{3 }+ catalyst + mixing continuous reactor”, *Naukovi Visti NTUU KPI*, no. 6, pp. 145–151, 2012 (in Ukrainian).

I.D. Lucheyko, “Mathematical modeling system “mixing continuous reactor + reaction A_{1 }®a_{2}A_{2}”at conditions of catalyst deactivation”, *Izvestiya Vuzov. Khimiya i Khim. Tekhnologiya*,* *vol. 57, no. 12, pp. 88–92, 2014 (in Russian).

I.D.* *Lucheyko,* *“Criterion Damkohler as parameter defining mathematical description of non-stationary operating mode of CSTR”, in *Proc**. XX Mendeleev** **Congress** **on** **General** **and** **Applied** **Chemistry*, vol. 3, Ekaterinburg, Russia, 2016.

I.D.* *Lucheyko,* *“Mathematical modeling of isothermal plug flow reactor with consecutive reaction taking into account the catalyst deactivation”, *Naukovi Visti NTUU KPI*, no. 5, pp. 106–115, 2016 (in Ukrainian). doi: 10.20535/1810-0546.2016.5.71920

## Downloads

## Published

## Issue

## Section

## License

Copyright (c) 2018 Igor Sikorsky Kyiv Polytechnic Institute

This work is licensed under a Creative Commons Attribution 4.0 International License.

Authors who publish with this journal agree to the following terms:- Authors retain copyright and grant the journal right of first publication with the work simultaneously licensed under CC BY 4.0 that allows others to share the work with an acknowledgement of the work's authorship and initial publication in this journal.
- Authors are able to enter into separate, additional contractual arrangements for the non-exclusive distribution of the journal's published version of the work (e.g., post it to an institutional repository or publish it in a book), with an acknowledgement of its initial publication in this journal.
- Authors are permitted and encouraged to post their work online (e.g., in institutional repositories or on their website) prior to and during the submission process, as it can lead to productive exchanges, as well as earlier and greater citation of published work