Analytical Rate Expressions Accounting for the Elementary Steps in

Publication Date (Web): April 26, 2017. Copyright ... The complexity of elementary-step reaction networks demands an adequate methodology to simulate ...
2 downloads 3 Views 725KB Size
Subscriber access provided by UB + Fachbibliothek Chemie | (FU-Bibliothekssystem)

Article

Analytical Rate Expressions Accounting for the Elementary Steps in Benzene Hydrogenation on Pt Luis Lozano, Guy B Marin, and Joris W. Thybaut Ind. Eng. Chem. Res., Just Accepted Manuscript • Publication Date (Web): 26 Apr 2017 Downloaded from http://pubs.acs.org on April 27, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Industrial & Engineering Chemistry Research is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Analytical Rate Expressions Accounting for the Elementary Steps in Benzene Hydrogenation on Pt Luis Lozano, Guy B. Marin and Joris W. Thybaut* Ghent University, Laboratory for Chemical Technology, Technologiepark 914, Ghent B-9052 Belgium KEYWORDS: Metal catalysis, microkinetic modeling, aromatics hydrogenation, Horiuti-Polanyi

mechanism

ABSTRACT: The complexity of elementary-step reaction networks demands an adequate

methodology to simulate the corresponding kinetics with a comprehensive, yet manageable rate expression. In this framework, a solution strategy has been developed for expressing the aromatics hydrogenation rate according to the Horiuti-Polanyi mechanism with a single analytical equation. Owing to the step-wise atomic hydrogen addition reaction network, a set of algebraic equations, linear in the surface coverage of free sites, is obtained. An approximate 30fold reduction in CPU time could be obtained by using the analytical solution as compared to the numerical one. The corresponding analytical solution produces simulation results that are, within the accuracy requested from the solver, identical to those obtained via the numerical solution. Finally, parameter estimates and statistical tests for the regression of the experimental data sustain the similarity between both simulation approaches. 1 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 32

INTRODUCTION Hydrogenation of aromatic compounds is an increasingly important reaction for improving diesel quality and reducing corresponding emissions.1, 2 More particularly, benzene is a compound widely produced in oil refineries, e.g. during the catalytic reforming of nafta. Benzene is frequently used as an industrial solvent; however, it is also known for its high octane number and as an important reactant for considerable (petro)chemical processes.3 Roughly 43 M ton/year of benzene are produced to supply these demands. Nonetheless, the widespread use of benzene leads to approximately 10 kg/ton of diffusive emissions, which raises environmental concerns since aromatic compounds are well-known air pollutants.4 Such concerns have been translated into environmental legislations around the globe to limit aromatic content in liquid fuels. In the last 10 years the maximum allowed concentration of benzene in gasoline has been reduced to j is assumed to be reversible, i.e. hydrogenation/dehydrogenation, and is, hence, calculated according to the general formulation in Eq.13.

z z ri → j = nei → jk (m, n)θiθH Ct − ne j→i k (m, n)−1θ jθ* Ct 2 2

(13)

ANALYTICAL SOLUTION FOR THE OVERALL REACTION RATE A methodology has been devised to calculate the overall benzene hydrogenation rate via an analytical expression while still accounting for the complex SEMK reaction mechanism. No symbolic computation was employed; hence, the arithmetic deduction provides a means to acquire more insight in the relation between the elementary reaction mechanism and its mathematical representation. In this methodology a set of equations is constructed to calculate all surfaces coverage when reactant and product adsorption and desorption are quasi equilibrated and the Pseudo SteadyState Approximation (PSSA) is applied to surface intermediates. Since hydrogen and benzene adsorption are in quasi-equilibrium (Eqs. 9 and 10), two equations are obtained to relate the surface coverage of both species to the corresponding partial pressures, see Eqs. 14 and 15.

12 ACS Paragon Plus Environment

Page 13 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

eq

θ H = K H 2 pH 2 θ*

(14)

eq

(15)

θ B = K B pBθ *

The other surface coverages are obtained by applying the PSSA and systematically elaborate the obtained equations. As result of the PSSA for all 13 “j” surface intermediates (Rj=0), Eq. 12 is reduced Eq. 16.

∑r

i→ j

g

=

∑r

(16)

j→k

h

Since θ13 can be assumed to be negligible, the CHA formation rate (rCHA) equals that of surface c-Hexyl (θ12) transformation into surface CHA (θ13) (r12->13) in Eq. 17: ~ z RCHA = rCHA = r12 →13 = ne12→13k (0,2)θ Hθ12 Ct 2

(17)

The number of single events and the rate coefficients are “given” values and θH can be obtained from Eq. 14. As a result, the calculation of the overall reaction rate, is reduced to the determination of θ12 via the equations obtained from the PSSA. Recognizing that in each of the rate expressions θH appears elevated to the first power, the PSSA results in a set of 12 linear algebraic equations in 12 unknowns that is perfectly specified thanks to hydrogen and benzene quasi-equilibrated adsorption. The set of algebraic equations resulting from the PSSA for the surface intermediates can be expressed in terms of the surface coverage of the free sites, preserving the linear character of the set of equations. First, applying the PSSA to surface BH according to Eq. 16 yields Eq. 18. r1→2 = r2→3 + r2→4 + r2→5

(18)

Expanding the reaction rate terms in Eq. 18, results in Eq. 19.

13 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 32

z z ne1→2 k (2,2)θ1θ H Ct − ne 2→1k (2,2) −1 θ 2θ* Ct 2 2 z z = ne 2→3k (1,2)θ 2θ H Ct − ne 3→2 k (1,2) −1 θ3θ* Ct 2 2 z z + ne 2→4 k (2,2)θ 2θ H Ct − ne 4→2 k (2,2) −1 θ 4θ* Ct 2 2 z z + ne 2→5 k (2,2)θ 2θ H Ct − ne 4→2 k (2,2) −1 θ5θ* Ct 2 2

(19)

Hydrogen and benzene surface coverage can be determined from Equations 14 and 15 thanks to the quasi equilibrated character of both adsorptions; through the proper substitutions and simplifications a single expression for θ2 is obtained, exhibiting a linear dependence between the surface coverage of BH and the remaining species and the coverage of the free sites, see Eq. 20. ~ ~ ~ ~ −1 −1 −1 ne1→ 2 k (2,2 )K B pB K H 2 pH 2 θ* + ne3→ 2k (1,2 ) θ3 + ne4→ 2 k (2,2 ) θ 4 + ne5→ 2 k (2,2 ) θ5 θ2 = ~ ~ ~ ~ −1 ne2→1k (2,2 ) + ne2→3k (1,2 ) + ne2→ 4 k (2,2 ) + ne2→5k (2,2 ) K H 2 pH 2

[

]

(20)

The coefficients of the surface coverage in Eq. 20 are function of the number of single events, rate coefficients and partial pressures, i.e., of “given” and/or observable quantities. These coefficients are denoted by letters to simplify the mathematical treatment of the different expressions (see Eq. 21). A2,3

~ ne3→ 2 k (1, 2 ) = ~ ~ ~ ~ −1 ne2 →1k (2,2 ) + ne2 → 3k (1,2 ) + ne2 → 4 k (2, 2 ) + ne2 → 5 k (2,2 ) K H 2 pH 2

[

(21)

]

Equation 21 can be therefore expressed in terms of these coefficients (see Eq. 22). θ 2 = A2 , 0θ * + A2 , 3θ 3 + A2 , 4θ 4 + A2 , 5θ 5

(22)

Similar operations can be performed for the surface intermediates further in the hydrogenation sequence, e.g., for 12CHD. Applying the PSSA, see Eqs. 23 and 24, provides a single expression to calculate θ3 in Eq. 25.

14 ACS Paragon Plus Environment

Page 15 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

r2 → 3 = r3→ 6 + r3→ 7

(23)

z z ne 2→3k (1,2)θ 2θ H Ct − ne3→2k (1,2) −1 θ3θ* Ct 2 2 z z = ne 3→6k (1,2)θ3θ H Ct − ne 6→3k (1,2) −1 θ 6θ* Ct 2 2 z z + ne 3→7 k (2,2)θ3θ H Ct − ne 7→3k (2,2) −1 θ7θ* Ct 2 2

(24)

~ ~ ~ −1 −1 ne2→3k (1,2 )θ 2 K H 2 pH2 + ne6→3k (1,2 ) θ 6 + ne7→3k (2,2 ) θ 7 θ3 = ~ ~ ~ −1 ne3→2 k (1,2 ) + ne3→6 k (1,2 ) + ne3→7 k (2,2 ) K H 2 pH 2

(25)

[

]

An expression for θ3, independent of θ2, is obtained when Eq. 19 is substituted in Eq. 22 (See Eq. 23). θ3

[n =

e2→3

]

~ ~ ~ −1 −1 k (1,2 )(A2,0θ* + A2,3θ 3 + A2,4θ 4 + A2,5θ 5 ) K H2 pH2 + ne6→3 k (1,2 ) θ 6 + ne7→3 k (2,2 ) θ 7 ~ ~ ~ −1 ne3→2 k (1,2 ) + ne3→6 k (1,2 ) + ne3→7 k (2,2 ) K H2 pH2

[

]

(26)

Coefficients that only contain kinetic parameters and observable quantities are again defined, e.g., B3,0 in Eq. 27, to simplify the mathematical expression. After Eq. 27 is solved for θ3, a single expression to calculate θ3 is obtained, see Eq. 28. B3, 0

~ ne2→3k (1,2 )A2 ,0 K H 2 pH 2 = ~ ~ ~ ~ −1 ne3→ 2 k (1,2) + K H 2 pH 2 ne3→6 k (1,2 ) + K H 2 pH 2 ne3→ 7 k (2,2 ) − A2,3 K H 2 pH 2 ne2 →3k (1,2 )

θ 3 = B3,0 θ * + B3,4 θ 4 + B3,5 θ 5 + B3,6 θ 6 + B3,7 θ 7

(27)

(28)

The latter equation is independent from surface species to which a lower number of hydrogen atoms have been added. Continuing these operations up to cyclohexyl (species 12 in Figure 2) requires a total of 33 substitutions, and results in a single, yet complex expression to solve the net rate of formation for gas-phase cyclohexane in Eq. 17. The latter equation is linear in the surface coverage of free sites and, hence, allows a straightforward calculation of the cyclohexane 15 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 32

formation rate. More elaborate guidelines for the derivation of the analytical rate expression according to the proposed methodology are given in the supporting information. The applicability of the methodology presented in this work is closely related to the generic features as occurring in a Horiuti–Polanyi mechanism, i.e., the surface reactions should constitute a series of, essentially, identical steps occurring on a single model. Species adsorption should be quasi-equilibrated and the Pseudo Steady State Approximation (PSSA) should be applicable to all surface reactions. This allows the formal linearization of the rate equations in the surface coverages of the hydrocarbon species and, hence, the step wise substitution of surface coverages by those of species situated earlier in the reaction network. Finally, product desorption being assumed to be fast and irreversible, ensures that the net rate of formation of the final product is independent of any additional equation related to product re-adsorption. As a result, the extension of the presented methodology to a diversity of hydrogenation reactions such as toluene, xylene, ethylene, etc. is straightforward since the Horiuti–Polanyi mechanism is widely accepted.11-19 Thanks to the symmetry present in benzene, the corresponding hydrogenation reaction network complexity is still relatively limited. The number of surface intermediates and hydrogen addition/abstraction reactions rapidly increases with the presence of substituents. For instance, a similar SEMK model for toluene hydrogenation, comprising 40 surface intermediate species and 100 reversible surface reactions, still exhibits the same characteristics in terms of surface coverage as the benzene hydrogenation network.18 Hence, the presented analytical solution strategy remains valid for such heavier aromatic species. Other reactions that exhibit similar features as the ones described above include the total and partial oxidation of hydrocarbons.41, 42 Also for addressing sub networks of more complex ones, such as the FTS,39, 40 the presented methodology opens up interesting perspectives to derive an 16 ACS Paragon Plus Environment

Page 17 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

analytical expression to simulate the rate of the consecutive insertion of a building block in a growing chain and, consequently, reduce the extent of the reaction network that needs to be solved numerically. ANALYTICAL VERSUS NUMERICAL SIMULATION RESULTS A validation of the developed methodology for the calculation of the overall aromatic hydrogenation rate based on an analytical rate expression is performed via comparison of these simulation results with those previously reported by Bera et al.,17 using a numerical solution strategy as well as with the available experimental data. Deviations between the two solution methodologies, i.e., analytical and numerical are expected to be within tolerance of the numerical routine used. Furthermore, the analytical solution reduces the CPU time from 1.186 s to 0.043 s for the 43 experiments, when it is compared with the numerical routine. The analytical rate expression ensures that a solution is obtained while, numerically, such a solution is not guaranteed when trying to solve the set of equations describing the reaction in steady state. Moreover, even if the CPU time gain from 1.186 to 0.043 s may seem insignificant at first sight, a 30-fold reduction in required CPU time potentially makes the difference in more complex reaction networks or when many simulations need to be performed consecutively, e.g., as in a regression. A parity diagram of experimental benzene conversion against the simulated results with the analytical solution is presented in Figure 3. Simulated values are in good agreement with the experimental data. The observed deviations can entirely be attributed to experimental error.17

17 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 32

Figure 3. Experimental and simulated benzene conversion at 5.0 - 10.0 H2/Benzene inlet molar ratio, 423 – 498 K temperature range, 1 – 3 MPa total pressure and 22.5 - 69.2 kgcat s mol-1 space time. The calculated benzene conversion is obtained by computing the outlet molar flow rate of Benzene in Eq. 1. The reaction term is solved by means of the methodology described in the section Analytical Solution for the Overall Reaction Rate, using the set of kinetic parameters estimated by Bera et al.17 Apart from the overall agreement between experimental data and model simulations, a more specific performance evaluation in terms of reproduction of temperature, pressure and hydrogen/benzene inlet molar flow rate effects is studied in Figures 4 and 5. The evolution of the benzene conversion with the total pressure and temperature is simulated by the analytical solution as presented in Figure 4. The maximum in the benzene hydrogenation conversion as a function of the temperature is adequately reproduced and stems from a decrease of the benzene surface concentration with increasing temperature.7

18 ACS Paragon Plus Environment

Page 19 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Figure 4. Benzene conversion as a function of the temperature at H2/Benzene inlet molar ratio of 5 and W/Fbenzene0 of 6.38 kg s mol-1: ● 3 MPa, ▲ 2 MPa, ♦ 1MPa. Full lines are computed by calculating the benzene outlet molar flow rate in Eq.1: __ 3 MPa, _ . 2 MPa, … 1MPa. The reaction term is solved by means of the methodology described in the section Analytical Solution for the Overall Reaction Rate, using the set of kinetic parameters estimated by Bera et al.17 Figure 5 shows the effect of the hydrogen/benzene molar ratio on the overall benzene conversion at 448 K and three total pressures. According to the reaction stoichiometry, three H2 molecules are needed for each benzene molecule to yield gas-phase cyclohexane. However, it is remarkable that only moderate conversions can be achieved even with hydrogen in excess at a molar ratio amounting to 5:1. The potentially competitive character of the H2 and benzene adsorption on the catalyst surface has been extensively discussed in the literature. Our data indicate that increasing benzene inlet partial pressure has a negative effect over the activity at those working conditions.2, 7, 11-13, 17, 18, 20, 30-33

19 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 32

Figure 5. Benzene conversion as a function of the H2/Benzene molar ratio at 448 K: 3 MPa, ▲ 2 MPa, ♦ 1MPa. Full lines are computed by calculating the Benzene outlet molar flow rate in Eq.1: __ 3 MPa, _ . 2 MPa, … 1MPa. The reaction term is solved by means of the methodology described in the section Analytical Solution for the Overall Reaction Rate, using the set of kinetic parameters estimated by Bera et al.17 It is clear from Figure 5 that the model simulations based on the analytical rate expressions allowed reproducing the experimentally observed trends. A similar conclusion was obtained by Bera et al. when adopting a numerical solution strategy.17 Figure 6 compares the results of the two solution strategies directly with each other, by displaying the differences in the benzene conversion obtained with both solution strategies as a function of the cyclohexane outlet flow rate.

20 ACS Paragon Plus Environment

Page 21 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Figure 6. Percentage difference on benzene conversion between the analytical solution described in section Analytical Solution for the Overall Reaction Rate, and the numerical solution achieved by Bera et al.,17 as a function of the cyclohexane outlet molar flow rate. Benzene conversion is calculated for both cases by computing the outlet molar flow rate of Benzene in Eq. 1, using the set of kinetic parameters estimated by Bera et al.17 The difference between the analytically and numerically calculated benzene conversion is randomly distributed with a maximum absolute value amounting to 8 10-2 %, which is perfectly in line with what could be expected based on the relative tolerance of the numerical solution in the DNSQE solver, which was set to 10-2.26 Kinetic parameter values were estimated making use of the analytical rate expression developed in this work. Table 1 compares the presently obtained estimates with those reported by Bera et al.17 Table 1. Kinetic parameter values (kJ mol-1) for the SEMK model for benzene hydrogenation. Kinetic parameters calculated in this work

Kinetic Parameters calculated by Bera et al.17 21

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Ea(0,2)

57.6 ± 0.3

57.5 ± 0.5

Ea(1,2)

65.3 ± 8.4

65.1 ± 8.8

∆H(0,2)

7.8 ± 0.2

7.9 ± 0.2

∆H(1,2)

1.1 ± 0.5

1.2 ± 0.5

∆HB

56.0 ± 0.5

56.0 ± 0.5

∆HH2

-59.4 ± 0.5

-59.4 ± 0.5

Page 22 of 32

Values obtained by regression are given with the corresponding 95% confidence interval, using the net rate of cyclohexane in Eq. 17 and the reactor model in Eq. 1. The reaction term is solved by means of the methodology described in the section Analytical Solution for the Overall Reaction Rate. Only small differences are encountered in the point estimates when using the analytical expression versus when relying on the numerical simulations. Moreover, the differences between the parameter estimates obtained according to both methodologies are statistically insignificant, and the already narrow confidence intervals are slightly reduced for the estimates obtained with the analytical rate expression. With the analytical rate expression a slightly higher F value for the global significance of the regression is obtained as compared to with the numerical solution, i.e., 4163 in the present work as compared to 4150 reported by Bera et al.17 The above statistical analysis, hence, indicates that the method to derive an analytical rate expression for aromatic hydrogenation leads to statistically identical parameter estimates. The differences, if any, are associated to the tolerance employed in the numerical solution and result in the latter being slightly less significant compared to the regression performed with the analytical expression. CONCLUSIONS Microkinetic models for aromatics hydrogenation are characterized by complex reaction networks involving a large number of surface intermediates. While is tempting to recur to

22 ACS Paragon Plus Environment

Page 23 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

numerical routines to solve the corresponding set of equations describing the kinetics in such a network, the particular nature of hydrogenation reactions brings an analytical solution within reach even when it doesn’t seem to be feasible at first sight due to highly tangled reaction network. According to the Horiuti – Polanyi mechanism as implemented for aromatics hydrogenation, reactant adsorption is quasi-equilibrated, hydrogen atom additions to unsaturated carbon atoms occurs step-wise, and desorption of the saturated product is fast and irreversible. Following these three assumptions, it was possible to solve a linear set of equations in the surface coverage of intermediate species and free sites allowing to calculate the cyclohexane formation rate, not having to rely on symbolic computation. The strategy proposed in our work yields statistically identical results as the numerical solution of the original set of equations and simultaneously reduces the required CPU time by a factor of 30, providing the guarantee that a solution is found in steady state simulations. The methodology has a generic character and can be directly extended towards, more complex aromatic components such as toluene and (ortho-,meta-,para-)xylene as well as to alternative reactions such as hydrocarbon oxidation or Fischer Tropsch synthesis. ASSOCIATED CONTENT

Supporting Information Further details on the methodology to derive analytical rate expressions for Benzene Hydrogenation on Pt discussed in this work. AUTHOR INFORMATION

Corresponding Author 23 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 32

*Email: [email protected]. NOMENCLATURE *

free site

1235THB

tetrahydrobenzene

1245THB

tetrahydrobenzene

123THB

trihydobenzene

124THB

trihydobenzene

135THB

trihydobenzene

12CHD

cyclohexadiene

14CHD

cyclohexadiene

13DHB

dihydrobenzene

A

pre-exponential factor (s-1, Pa-1 s-1)

B

benzene

BH

hydrobenzene

Ct

Total concentration of active sites (mol/kgcat)

c-Hexyl

cyclohexyl

CHA

cyclohexane

CHE

cyclohexene

DFT

Density Functional Theory

F

molar flow rate (mol s-1)

Fvalue

F value for the significance of the regression

b

model vector

h

Planck constant (6.63 x 10-34 J s)

H

enthalpy (J mol-1)

k

rate coefficient (s-1, Pa-1 s-1)

kb

Boltzmann constant (1.38 x 10-23 J K-1) 24 ACS Paragon Plus Environment

Page 25 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

m

number of unsaturated neighbor carbon atoms

n

saturation degree of the reactant carbon

ne

number of single events

nob

number of observations

npar

number of parameters

nresp

number of responses

p

pressure (bar)

S

entropy (J mol-1 K-1)

SSQ

residual sum of squares between the experimental and model calculated outlet molar flow rates (mol2 s-2)

T

temperature (K)

r

reaction rate (mol kgcat-1 s-1)

R

universal gas constant (8.314 J mol-1 K-1)

Rj

net production rate of component j (mol kgcat-1 s-1)

Wcat

mass catalyst (kg)

Greek symbols β

real parameter vector

σ

symmetry numbers

θ

surface coverage



difference

Subscripts ext

external symmetry

int

internal symmetry

gl

global symmetry

rot

rotational

t

total

Superscripts 25 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

n

number of chiral centers

r

reactant



transition state

~

single event

˄

model calculated value

o

standard state

0

inlet

Page 26 of 32

REFERENCES 1.

Metaxas, K. C.; Papayannakos, N. G., Studying the internal mass transfer phenomena

inside a Ni/Al2O3 catalyst for benzene hydrogenation. Chem Eng J 2008, 140, 352-357. 2.

Canduela-Rodriguez, G.; Sabbe, M. K.; Reyniers, M. F.; Joly, J. F.; Marin, G. B.,

Thermodynamic study of benzene and hydrogen coadsorption on Pd(111). Phys Chem Chem Phys 2014, 16, 23754-23768. 3.

Singh, U. K.; Vannice, M. A., Kinetic and thermodynamic analysis of liquid-phase

benzene hydrogenation. AIChE J 1999, 45, 1059-1071. 4.

Bahadar, H.; Mostafalou, S.; Abdollahi, M., Current understandings and perspectives on

non-cancer health effects of benzene: A global concern. Toxicol Appl Pharm 2014, 276, 83-94. 5.

Sun, H. J.; Jiang, H. B.; Li, S. H.; Dong, Y. Y.; Wang, H. X.; Pan, Y. J.; Liu, S. C.; Tang,

M. S.; Liu, Z. Y., Effect of alcohols as additives on the performance of a nano-sized RuZn(2.8%) catalyst for selective hydrogenation of benzene to cyclohexene. Chem Eng J 2013, 218, 415-424.

26 ACS Paragon Plus Environment

Page 27 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

6.

Roy, S.; Datta, S., Hydrogenation of Toluene on Zirconium-Modified Hexagonal

Molecular Sieve Supported Platinum and Palladium Catalysts. Ind Eng Chem Res 2013, 52, 17360-17368. 7.

Keane, M. A.; Patterson, P. M., The role of hydrogen partial pressure in the gas-phase

hydrogenation of aromatics over supported nickel. Ind Eng Chem Res 1999, 38, 1295-1305. 8.

Yan, X. H.; Zhang, Q.; Zhu, M. Q.; Wang, Z. B., Selective hydrogenation of benzene to

cyclohexene over Ru-Zn/ZrO2 catalysts prepared by a two-step impregnation method. J Mol Catal A-Chem 2016, 413, 85-93. 9.

Stanislaus, A.; Cooper, B. H., Aromatic Hydrogenation Catalysis - A Review. Catal Rev

1994, 36, 75-123. 10. Kumar, S. A. K.; John, M.; Pai, S. M.; Niwate, Y.; Newalkar, B. L., Low temperature hydrogenation of aromatics over Pt-Pd/SiO2-Al2O3 catalyst. Fuel Process Technol 2014, 128, 303-309. 11. Vanmeerten, R. Z. C.; Verhaak, A. C. M.; Coenen, J. W. E., Gas-Phase Benzene Hydrogenation on a Nickel-Silica Catalyst .2. Gravimetric Experiments of Benzene, Cyclohexene, and Cyclohexane Adsorption and Benzene Hydrogenation. J Catal 1976, 44, 217225. 12. Vanmeerten, R. Z. C.; Degraaf, T. F. M.; Coenen, J. W. E., Gas-Phase Benzene Hydrogenation on a Nickel-Silica Catalyst .3. Low-Field Magnetization Measurements on Hydrogen, Benzene, Cyclohexene and Cyclohexane Adsorption, and Benzene Hydrogenation. J Catal 1977, 46, 1-12. 27 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 32

13. Vanmeerten, R. Z. C.; Coenen, J. W. E., Gas-Phase Benzene Hydrogenation on a NickelSilica Catalyst .4. Rate Equations and Curve Fitting. J Catal 1977, 46, 13-24. 14. Lin, S. D.; Vannice, M. A., Hydrogenation of Aromatic-Hydrocarbons over Supported Pt Catalysts .1. Benzene Hydrogenation. J Catal 1993, 143, 539-553. 15. Lin, S. D.; Vannice, M. A., Hydrogenation of Aromatic-Hydrocarbons over Supported Pt Catalysts .2. Toluene Hydrogenation. J Catal 1993, 143, 554-562. 16. Lin, S. D.; Vannice, M. A., Hydrogenation of Aromatic-Hydrocarbons over Supported Pt Catalysts .3. Reaction Models for Metal-Surfaces and Acidic Sites on Oxide Supports. J Catal

1993, 143, 563-572. 17. Bera, T.; Thybaut, J. W.; Marin, G. B., Single-Event Microkinetics of Aromatics Hydrogenation on Pt/H-ZSM22. Ind Eng Chem Res 2011, 50, 12933-12945. 18. Bera, T.; Thybaut, J. W.; Marin, G. B., Extension of the Single-Event Microkinetic Model to Alkyl Substituted Monoaromatics Hydrogenation on a Pt Catalyst. ACS Catal 2012, 2, 1305-1318. 19. Mattson, B.; Foster, W.; Greimann, J.; Hoette, T.; Le, N.; Mirich, A.; Wankum, S.; Cabri, A.; Reichenbacher, C.; Schwanke, E., Heterogeneous. Catalysis: The Horiuti-Polanyi Mechanism and Alkene Hydrogenation. J Chem Educ 2013, 90, 613-619. 20. Saeys, M.; Reyniers, M. F.; Neurock, M.; Marin, G. B., Ab initio reaction path analysis of benzene hydrogenation to cyclohexane on Pt(111). J Phys Chem B 2005, 109, 2064-2073.

28 ACS Paragon Plus Environment

Page 29 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

21. Saeys, M.; Reyniers, M. F.; Thybaut, J. W.; Neurock, M.; Marin, G. B., First-principles based kinetic model for the hydrogenation of toluene. J Catal 2005, 236, 129-138. 22. Thybaut, J. W.; Saeys, M.; Marin, G. B., Hydrogenation kinetics of toluene on Pt/ZSM22. Chem Eng J 2002, 90, 117-129. 23. Bos, A. N. R.; Lefferts, L.; Marin, G. B.; Steijns, M. H. G. M., Kinetic research on heterogeneously catalysed processes: A questionnaire on the state-of-the-art in industry. Appl Catal A-Gen 1997, 160, 185-190. 24. Steijns, M.; Froment, G. F., Hydroisomerization and Hydrocracking .3. Kinetic-Analysis of Rate Data for Normal-Decane and Normal-Dodecane. Ind Eng Chem Prod Rd 1981, 20, 660668. 25. Berty, J. M., Reactor for Vapor-Phase Catalytic Studies. Chem Eng Prog 1974, 70, 7885. 26. Hiebert, K. L. DNSQE. http://www.netlib.org/ 27. Rogers, J. E. ODRPACK Software for Weighted Orthogonal Distance Regression. http://www.netlib.org 28. Martens, J. A.; Vanbutsele, G.; Jacobs, P. A.; Denayer, J.; Ocakoglu, R.; Baron, G.; Arroyo, J. A. M.; Thybaut, J.; Marin, G. B., Evidences for pore mouth and key-lock catalysis in hydroisomerization of long n-alkanes over 10-ring tubular pore bifunctional zeolites. Catal Today 2001, 65, 111-116.

29 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 30 of 32

29. Saeys, M.; Reyniers, M. F.; Marin, G. B.; Neurock, M., Density functional study of benzene adsorption on Pt(111). J Phys Chem B 2002, 106, 7489-7498. 30. Saeys, M.; Reyniers, M. F.; Neurock, M.; Marin, G. B., Density functional theory analysis of benzene (De)hydrogenation on Pt(111): Addition and removal of the first two Hatoms. J Phys Chem B 2003, 107, 3844-3855. 31. Taimoor, A. A.; Pitault, I., Kinetics of toluene hydrogenation-integrating a dynamic approach regarding catalyst activity. React Kinet Mech Cat 2011, 102, 263-282. 32. Bratlie, K. M.; Flores, L. D.; Somorjai, G. A., In situ sum frequency generation vibrational spectroscopy observation of a reactive surface intermediate during high-pressure benzene hydrogenation. J Phys Chem B 2006, 110, 10051-10057. 33. Fergusson, G. A.; Vorotnikov, V.; Wunder, N.; Clark, J.; Gruchalla, K.; Bartholomew, T.; Robichaud, D. J.; Beckham, G. T., Ab Initio Surface Phase Diagrams for Coadsorption of Aromatics and Hydrogen on the Pt(111) Surface. J Phys Chem C 2016, 120, 26249-26258. 34. Toppinen, S.; Rantakyla, T. K.; Salmi, T.; Aittamaa, J., Kinetics of the liquid-phase hydrogenation of benzene and some monosubstituted alkylbenzenes over a nickel catalyst. Ind Eng Chem Res 1996, 35, 1824-1833. 35. Dewachtere, N. V.; Santaella, F.; Froment, G. F., Application of a single-event kinetic model in the simulation of an industrial riser reactor for the catalytic cracking of vacuum gas oil. Chem Eng Sci 1999, 54, 3653-3660.

30 ACS Paragon Plus Environment

Page 31 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

36. Baltanas, M. A.; Vanraemdonck, K. K.; Froment, G. F.; Mohedas, S. R., Fundamental Kinetic Modeling of Hydroisomerization and Hydrocracking on Noble-Metal-Loaded Faujasites .1. Rate Parameters for Hydroisomerization. Ind Eng Chem Res 1989, 28, 899-910. 37. Baltanas, M. A.; Froment, G. F., Computer-Generation of Reaction Networks and Calculation of Product Distributions in the Hydroisomerization and Hydrocracking of Paraffins on Pt-Containing Bifunctional Catalysts. Comput Chem Eng 1985, 9, 71-81. 38. Froment, G. F., Single event kinetic modeling of complex catalytic processes. Catal Rev

2005, 47, 83-124. 39. Lozano-Blanco, G.; Thybaut, J. W.; Surla, K.; Galtier, P.; Marin, G. B., Single-event microkinetic model for Fischer-Tropsch synthesis on iron-based catalysts. Ind Eng Chem Res

2008, 47, 5879-5891. 40. Lozano-Blanco, G.; Thybaut, J. W.; Surla, K.; Galtier, P.; Marin, G. B., Simulation of a Slurry-Bubble Column Reactor for Fischer-Tropsch Synthesis Using Single-Event Microkinetics. AIChE J 2009, 55, 2159-2170. 41. Auer, R.; Thyrion, F. C., Kinetics of the total oxidation of methane over a La(0.9)Ce(0.1)CoO(3) perovskite catalyst. Ind Eng Chem Res 2002, 41, 680-690. 42. Heracleous, E.; Machli, M.; Lemonidou, A. A.; Vasalos, L. A., Oxidative dehydrogenation of ethane and propane over vanadia and molybdena supported catalysts. J Mol Catal A-Chem 2005, 232, , 29-39.

31 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 32

FOR TABLE OF CONTENTS ONLY

32 ACS Paragon Plus Environment