MATEC Web of Conferences 49, 001 01 (2016) DOI: 10.1051/matecconf/201649001 01 C Owned by the authors, published by EDP Sciences, 2016

Modeling of with

Elena Khlebnikova 1, a, Irena Dolganova 1, Elena Ivashkina 1, Stanislav Koshkin 1, 2 1Tomsk Polytechnic University, Tomsk, Russia 2 LLC NIOST, Tomsk, Russia

Abstract. The paper considers the benzene alkylation with ethylene model development with the use of catalyst. A list of reactions occurring in the alkylation reactor was made and the thermodynamic possibility of determination of these reactions by the change of Gibbs energy was defined. The paper presents the transformation scheme, which includes the grouping of components on the basis of their reactivity and the degrees of compensation values of the corresponding reactions. Drawing on the obtained data the authors developed the kinetic model of the alkylation of benzene with ethylene.

1 Introduction Since the 70s of the last century it became possible to produce with the use of heterogeneous zeolite catalyst on a commercial scale. Heterogeneous catalytic process gradually started to replace the old The process takes place at T = 220 - 255 °Cand homogeneous catalytic technology. Nowadays most of P=3.4 MPa. This reaction is exothermic, which leads to the existing ethylbenzene production plants operate using gradual temperature increase from one catalyst bed to the the zeolite catalyst [1, 2]. According to the data of 2014 following one. This is the reason why alkylation reactors total world ethylbenzene production capacity is around are typically equipped with intermediate coolers for 37 millions tons per year [3]. keeping the total temperature rise. Technology of ethylbenzene production includes the Then ethylbenzene converts to produce small amounts stages of alkylation and . The control of of polyethylbenzenes (PEB): technological parameters such as benzene-to-ethylene ratio in the feed and temperature lets to conduct the process under the optimum conditions and monitor the output of ethylbenzene [4, 5]. Although the technology developers provide their recommendations on the values of these parameters, the development of the model of ethylbenzene synthesis process and its subsequent modeling aimed at the increase of the resource and energy efficiency of the production plant is one of the most feasible ways of customization of the process to particular cases [6]. These reactions are the main side ones and generate 2 Object of the research considerable quantities of heat being exothermic ones. PEB are separated from ethylbenzene and then The ethylbenzene technology is developed for converted to EB in the transalkyation reactor at T = 190 - ethylbenzene production with the use of the liquid phase 220 ºC and P = 3 - 3.5 MPa. and zeolite catalyst. By now the chemistry of ethylbenzene synthesis is The industrial plant operates continuously and uses 2 quite well understood [7, 8]. In addition to the main different reactors for alkylation and transalkyation. reactions leading to ethylbenzene formation it’s known a In alkylation reactor ethylene in its liquid phase goes number of side reactions. Altogether these reactions are a through almost a complete reaction accompanied by the challenging task from the perspective of mathematical excess of benzene, which is necessary to form description and its subsequent use. Representation and ethylbenzene (EB): adequacy of a mathematical model and its ability to be

a Corresponding author: [email protected]

This is an Open Access article distributed under the terms of the Creative Commons Attribution License 4.0, which permits distribution, and reproduction in any medium, provided the original work is properly cited. MATEC Web of Conferences

approximate enough to the real data depend on reaction Formalization was performed by means of uniting scheme of the model. This scheme should be quite simple, individual substances into the pseudo-components having but at the same time accurate enough to describe the similar reactivity in respect of a particular reaction. The process [9]. determination of compensation degree, which is characterized by the partial compensation of the energy expended on the breaking of the old bonds by the energy 3 Model development of the new bonds formation on the reaction route, was Reactions that occur in alkylation and transalkyation considerably helpful in solving the problems of catalytic processes and that are based on different types of the action prediction. zeolite catalyst were investigated by certain authors [10- Therefore, further step is to study the degree of 12]. However, alkylation and transalkyation have compensation for the list of considered reactions. different conversion processes schemes, and at the same ∑ D − Eа time they have quite much in common. So the main task H = i ⋅100% (1) ∑ D is to determine the pathway of the development of i secondary reactions: transalkyation, disproportionation, where – the energy of bounds breaking, kJ / mol, dealkylation, and other ones, the result of which ∑ Di is the formation of by-products. Some steps were taken to Еа – energy of activation, kJ / mol. make up a detailed reaction scheme based on industrial The main problem in calculation of the compensation data and the existing knowledge about the process [13]. energy is to determine the activation energy of the We offer the method of designing of the math model that reaction. combines industrial data, research results and quantum- As the first approximation, the activation energy was chemical calculation. calculated according to Polanyi-Semenov equation for The list of possible reactions for benzene alkylation exothermic reactions [14]. was made with account of the quality of sources = + ⋅ Δ (ethylene and benzene), product composition and aspects Ea 11.5 0.25 H (2) of the process chemistry. Е – Δ – Thermodynamic possibility of these reactions has where а energy of activation, kcal / mol, H heat of been confirmed by determination of the isobaric- reaction, kcal / mol. Enthalpy of the reactions was calculated with the use isothermal parameter – the change of Gibbs energy ΔG at T=510 ºK and P=3.4 MPa. As for the target reaction of quantum-chemical methods of thermodynamic functions determination. The calculated enthalpy of Gibbs energy is equal to -43.4 kJ / mol and according to all defined ΔG of the considered reactions, alkylation reaction of benzene with ethylene (-117.9 kJ / thermodynamic possibility of these reactions is probable. mol) was compared with the literature data (-113 kJ / mol) The determination was performed with the use of [15] as a proof of the relevance of the method for the unknown enthalpy of some reactions in the literature. All Gaussian-98 program by PM3 semi-empirical method. The graphical description of the reactions pathway of the considered reactions are exothermic. benzene alkylation is shown in Figure 1. Among them Then the activation energy of the given reaction was there are benzene alkylation of ethylene and propylene calculated in first approximation with the use of the (ethylene feed impurity), (benzene feed impurity) empiric equation. The value of the reaction mentioned above is 78.5 kJ / mol in comparison with 63.4 kJ / mol alkylation, next step of alkylation caused by the PEB ’ formation, heavy by-product formation (diphenylethane) in [16]. It s important to note that activation energy that through the intermediate formation in the process of was found for all reaction refers to the non-catalytic incomplete benzene alkylation. process, and it was good for the first approximation on next stage of the parameter estimation. The meanings of defined activation energy according to Polanyi-Semenov equation at T=510 ºK and P=3.4 MPa are represented in Table 1. However, it might be improved with the use of literature sources of activation energy values for some reactions. After that the activation energy could be calibrated to the activation energy of the reaction. The next step in calculating the degree of compensation is the determination of the energy of bonds breaking. Let’s consider the methods based on the example of the above-mentioned reaction. Figure 1. Detailed scheme of hydrocarbons conversion in the To form ethylbenzene it’s necessary to break the two process of alkylation of benzene with ethylene bonds: C-H bond in a of benzene and C-H bond in a molecule of ethylene. The energies of these bonds The scheme needs to be formalized (simplified) in have been calculated with the use of Gaussian-98. order it could be used for future computations; at the The energy of the first bond was 433.8 kJ / mol, the same time it is to provide the highly accurate description energy of the second one was 215.5 kJ / mol. The energy of the process. of bonds breaking was equal to 649.3 kJ / mol. Similarly, the energy of bonds breaking was calculated for the rest

01001-p.2 ICCCP 2016 of the reactions, and then the degree of compensation was Butylene = calculated by the formula (1). The results are shown in Ethylbutylbenzene Table 1. The literature values of the energy of bonds Toluene + breaking and their comparison with the determined Ethylene = 108.4 747.6 715.0 4.6 89.9 values are also shown in Table 1. Ethyltoluene 6Benzene = 1179. Coronene + 5206.0 4980.0 4.5 90.3 Table 1. Results of determined degrees of reaction 8 compensation Butylene = cis- 183.1 952.5 1014.0 6.1 89.0 Е а ∑ D , ∑ , i (Polanyi- Di Δ , cis-Butene + cis- kJ / mol ∑Di H Reaction Semenov kJ / , Butene = (literature % 166.3 809.6 830.0 2.5 87.8 ), kJ / mol % Dicyclobutane + mol ) Hydrogen Benzene + Ethylene = 78.5 649.3 715.0 9.2 90.2 Computation results show the good comparability of Ethylbenzene the determined values of energy of bonds breaking with Ethylbenzene + the literature data. Ethylene = 86.1 633.5 715.0 11.4 92.6 Diethylbenzene Diethylbenzene + 4 Results and discussions Ethylene = 90.7 617.7 715.0 13.6 88.9 Triethylbenzene In order to formalize the transformation schemes, Triethylbenzene + reactions were divided into two groups by the criterion of Ethylene = 84.4 602.0 715.0 15.8 89.1 the similarity of their values of compensation degree Tetraethylbenzene (Table 2). The degree of compensation for the same Tetraethylbenzene + Ethylene = 84.8 586.2 715.0 18.0 88.8 group of reactions differs approximately by 1% [17-19]. Pentaethylbenzene Table 2. Alkylation reactions grouping by their degrees of Pentaethylbenzene compensation + Ethylene = 93.2 570.4 715.0 20.2 87.9 Hexaethylbenzene Reaction group Average H,% 2Benzene + for reaction Ethylene= 94.1 1083.2 1100.0 1.5 93.6 group Diphenylethane + 1. Ethylbenzene formation 90.2 Hydrogen Benzene + Ethylene = Ethylbenzene Ethylene + 2. Diethylbenzene formation 92.6 Ethylene = 77.7 882.2 955.0 7.6 92.8 Ethylbenzene + Ethylene = Butylene Diethylbenzene Benzene + 3. Triethylbenzene formation 88.9 Butylene = 105.8 762.9 685.0 11.4 90.3 Diethylbenzene + Ethylene = Butylbenzene Triethylenebenzene Benzene + 4. Heavy compounds formation 89.8 Propylene = n- 99.1 788.2 685.0 15.1 90.9 Triethylbenzene + Ethylene = Propylbenzene Polyethylenebenzene Benzene + 2Benzene + Ethylene = Diphenylethane + Propylene = 107.9 788.2 685.0 15.1 90.5 Hydrogen 5. Other monoalkylates formation 90.6 Benzene + Benzene + Butylene = Monoalkylate Ethylene = o- 63.8 1298.7 1428.0 9.1 95.6 Benzene + Propylene = Monoalkylate 95.7 Benzene + 6. Xylene formation Ethylene = m- 55.0 1298.7 1428.0 9.1 95.8 Benzene+ Ethylene=Xylene Xylene 7. Other dialkylates formation 89.4 Benzene + Ethylbenzene + Propylene = Dialkylate Ethylene = p- 55.0 1298.7 1428.0 9.1 95.8 Ethylbenzene + Butylene = Dialkylate Xylene Monoalkylate + Ethylene = Dialkylate Cumene + Toluene + Ethylene = Dialkylate Ethylene = 94.5 633.0 715.0 11.5 89.0 8. Coke formation 90.3 Ethylcumene 6Benzene=Coke+Hydrogen Ethylbenzene + 9. Cyclo and dicyclo formation 88.4 Propylene = 123.5 772.3 715.0 8.0 89.5 Butylene = cis-Butene Ethylcumene cis-Butene + cis-Butene = Dicyclobutane Butylbenzene + + Hydrogen Ethylene = 92.8 635.4 715.0 11.1 89.1 10. Butylene formation 92.8 Ethylbutylbenzene Ethylene + Ethylene = Butylene Ethylbenzene + 120.5 747.1 715.0 4.5 89.3

01001-p.3 MATEC Web of Conferences

Thus, it is possible to identify the number of pseudo- Similar work was carried out with the transalkylation components: the heavy compounds which are composed process. Kinetic models of alkylation and transalkylation of individual components such as tetra-, penta-, are presented in Table 4. hexaethylbenzene and diphenylethane; monoalkylates (butylbenzene, n-propylbenzene, cumene), dialkylates Table 4. Kinetic models of alkylation and transalkylation (ethylcumene, ethylbutylbenzene), (o-, m-, p- The kinetic model of The kinetic model of xylene), cyclic and dicyclic hydrocarbons. Formalized alkylation transalkylation scheme of hydrocarbons transformations is shown in Figure 2. ⎧dC ⎧ ⎪ benzene =−W − 2W −W −W − ⎪ 1 4.2 5.1 5.2 dCbenzene =− − − − − ⎪ dt ⎪ W1 W2 W3 W4.1 − − dt ⎪ W6 6W8 ⎪ ⎪ ⎪−W −W − 6W − 2W dC 4.2 5.1 7 8 ⎪ ethylene =− − − − − ⎪ W1 W2 W3 W4.1 dC ⎪ ethylene = dt ⎪ W8 ⎪− − − − − ⎪ dt W4.2 W6 W7..3 W7.4 2W10 ⎪ ⎪dC dC diethylbenzene =− + ⎪ butylene =− − − + ⎪ W 1 W 2 W5.1 W7.2 W9.1 W10 dt ⎪ dt ⎪ ⎪ ⎪dC dC triethylbenzene =− + ⎪ propylene =− − ⎪ W 2 W3 W5.2 W7.1 dt ⎪ dt ⎪ ⎪ dC dC ⎪ polyethylbenzene =− ⎪ cis-butene = − 2 W 3 W9.1 W9.2 ⎪ dt ⎪ dt ⎪ ⎪ dCethylbenzene dCdicyclobutane ⎪ = 2 + + + 2 + ⎪ = W W 1 W2 W3 W4.1 W5.1 dt 9.2 ⎪ dt Figure 2. Formalized scheme of hydrocarbons transformations ⎪ ⎪ dCbutylbenzene ⎪dChydrogen =− − in benzene in process of its alkylation with ethylene = + + ⎪ W4.1 W4.2 ⎪ W8 W4.2 W9.2 dt ⎪ dt ⎪ ⎪ dC dCmonoalkylate ⎨ ethylbenzene = − − − ⎨ = 2W +W According to the law of mass action the expressions W 1 W2 W7.1 W7.2 4.2 5.1 ⎪ dt ⎪ dt of the rates of alkylation and transalkylation reactions can ⎪ ⎪dC dCdiethylbenzene dioalkylate =− − ⎪ = W −W ⎪ W 5.1 W5.2 be put as follows (Table 3). ⎪ dt 2 3 ⎪ dt ⎪ dC ⎪dChydrogen triethylbenzene = − =−W +W +W +W ⎪ W 3 W4 ⎪ 5.2 6 7 8 Table 3. Reaction rate expressions in alkylation and ⎪ dt ⎪ dt ⎪ ⎪dC transalkylation dCheavy_compounds xylene = = W +W W 5.2 ⎪ 4.1 4.2 ⎪ dt ⎪ dt ⎪ dC ⎪dCmonoalkylate ⎪ = Reaction group Reaction rate expression = W +W −W W 5.2 ⎪ 5.1 5.2 7.3 ⎪ 1. Ethylbenzene −Ea(1)/RT dt dt W = k ⋅ e ⋅C ⋅C ⎪ ⎪dC 1 0(1) benzene ethylene ⎪dCdialkylate cis-butene =− formation = W +W +W +W ⎪ 2W 6 ⎪ 7.1 7.2 7.3 7.4 dt − dt ⎪ 2. Diethylbenzene Ea(2)/RT ⎪ ⎪dC = ⋅ ⋅ ⋅ dCxylene dicyclobutane = W2 k0(2) e Cethylbenzene Cethylene ⎪ = ⎪ W 6 formation W 6 ⎪ dt ⎪ dt 3. Triethylbenzene −Ea(3)/RT ⎪ ⎪dC = ⋅ ⋅ ⋅ dCtoluene diphenylethane = W k e C C ⎪ =−W W 8 3 0(3) diethylbenzene ethylene 7.4 ⎪ dt formation ⎪ dt ⎪ −Ea(4.1)/RT ⎪dC dC 4. Heavy = ⋅ ⋅ ⋅ coke = ⎪ coke = W W4.1 k0(4) e Ctriethylbenzene Cethylene ⎩⎪ W8 ⎩ 7 compounds dt dt = ⋅ −Ea(4.2)/RT ⋅ 2 ⋅ formation W4.2 k0(4) e Cbenzene Cethylene − Initial conditions are: =0, С =C ,where – 5. Other = ⋅ Ea(5.1)/RT ⋅ ⋅ t i 0i i W5.1 k0(5) e Cbenzene Cbutylene monoalkylates corresponding . = ⋅ −Ea(5.2)/RT ⋅ ⋅ formation W5.2 k0(5) e Cbenzene Cpropylene The kinetic models presented in the paper are to become the basis for the computer simulation of the 6. Xylene W = k ⋅ e−Ea(6)/RT ⋅C ⋅C formation 6 0(6) benzene ethylene ethylbenzene production system development. = ⋅ −Ea(7.1)/RT ⋅ ⋅ W7.1 k0(7) e Cethylbenzene Cpropylene 7. Other W = k ⋅ e−Ea(7.2)/RT ⋅C ⋅C 5 Conclusion dialkylates 7.2 0(7) ethylbenzene butylene = ⋅ −Ea(7.3)/RT ⋅ ⋅ formation W7.3 k0(7) e Cmonoalkylate Cethylene Thermodynamic analysis of alkylation and − transalkylation reactions were performed with the use of W = k ⋅e Ea(7.4)/RT ⋅C ⋅C 7.4 0(7) toluene ethylene quantum-chemical methods of calculating 8. Coke = ⋅ −Ea(8)/RT ⋅ 6 W8 k0(8) e Cbenzene thermodynamic functions. The method of formalization formation of the hydrocarbons transformation scheme, which 9. Cyclo and = ⋅ −Ea(9.1)/RT ⋅ W9.1 k0(9) e Cbutylene consists in the grouping of components by their reactivity dicyclo W = k ⋅e−Ea(9.2)/RT ⋅C2 drawing on the compensation degree values of the formation 9.2 0(9) cis-butene corresponding reactions, was developed. The correctness − 10. Butylene = ⋅ Ea(10 )/RT ⋅ 2 W10 k0(10) e Cethylene of the determination of the energy of bonds breaking in formation the reagent’s is confirmed by comparison of its results with the literature data. The kinetic models of where – rate of i – group reaction; – pre- Wi k0(i) alkylation of benzene with ethylene and exponential factor for i – group reaction; Ea – polialkylbenzenes transalkylation were developed on the (i, j) basis of the formalized transformations schemes. activation energy of j - subgroup for i – group reaction.

01001-p.4 ICCCP 2016

References 11. F. Canavi, MCorazzari, E. Bencini, G.Goffredi, App. Cat. A:General J. 226, 31-40 (2004) 1. Thomas F. Degnan Jr, C Morris Smith, Chaya R, 12. Tseng-Chang Tsai, Shang-Bin Liu, Ikai Wang, App. Venkat, Appl. Catal. A: Gen. J. 221, 283-294 (2001) Cat. A:General J. 181, 355-398 (1999) 2. C.Perego, P. Ingallina, Catal. Today, 73, 3-22 (2002) 13. N. Hamendi, D. Iranshahi, M.R. Rahimpour, S. 3. Jennifer Carvajal-Diaz, David Byrne CEH (2015) Raeissi, H. Rajaei, Journal of the Taiwan Inst. Of 4. Robert A. Meyers, Handbook of Heterogeneous Chem. Eng J. 48, 56-72 (2015) , McGRAW-HILL, USA, 665 (2005) 14. Anthony G. Dixon, Advances in Chemical 5. N. Ebrahimi Ali, Z. Sharak Ashkan, A. Mousavi Engineering, USA, 301 (2014) Seyyed, Faeze Aghazadeh, Ahmad Soltani, Chemical 15. Kgutso Mokaema -Transfer (Transalkylation) Engineering and Processing, 50, 31-36 (2011) reactions of alkylaromatics on solid catalyst, 6. Xiaoming Chen, Shiping Huang, Dapeng Cao Johannesburg, 68 (2005) Fluid.Ph.Equilibria, 260, 146-152 (2007) 16. Hamid Ganji, Jafar S, Ahari et all, & J. 7. V.A. Fetisova, E.D. Ivashkina, E.N., Ivanchina, A.V. 46, 1, 55-63 (2004) Kravtsov, Catal. In Industry 2-1, 55-61 (2010) 17. Howard F. Rase, Handbook of Commercial Catalyst. 8. T.Odedairo, S. Al-Khattaf, Chemical Engineering J. Heterogeneous Catalyst, CRC Press, USA, 551 167, 240-254 (2011) (2000) 9. Evgeniya Frantsina, Elena Ivashkina, Emiliya 18. Nurmakanova A. E., Salishcheva A. A., Chudinova Ivanchina, Rostislav Romanovskii, Chemical A. A. , Syskina A. A., Ivashkina Procedia Chemistry, Engineering Journal (available online 28 February 10, 284-288 (2014) 2015) 19. Khlebnikova E. S., Bekker A. V., Ivashkina E. N., 10. L. Forni, G. Cremona, F. Missineo, G. Bellusi, C. Procedia Chemistry, 10, 297-304 (2014) Perego, G. Pazzuconi, App. Cat. A:General J. 121, 261-272 (1995)

01001-p.5