Single event kinetic modeling of solid acid alkylation of isobutane with butenes over proton-exchanged Y-zeolites

Complex reaction kinetics of the solid acid alkylation of isobutane with butenes over a proton-exchanged Y-zeolite has been modeled at the elementary step level. Starting with a computer algorithm that generated the reaction network based on the fundamentals of the carbenium ion chemistry, the formation of over 100+ product species has been modeled in order to gain understanding of the underlying phenomena leading to rapid catalyst deactivation and product selectivity shifts observed in experimental runs. An experimental investigation of the solid acid alkylation process was carried out in a fixed bed catalytic reactor operating with an excess of isobutane under isothermal conditions at moderate temperatures (353-393 K) in liquid phase. Experimental data varying with run-time for a set of butene space-times and reaction temperatures were collected for parameter estimation purposes. A kinetic model was formulated in terms of rate expressions at the elementary step level including a rigorous modeling of deactivation through site coverage. The single event concept was applied to each rate coefficient at the elementary step level to achieve a significant reduction in the number of model parameters. Based on the identification of structural changes leading to the creation or destruction of symmetry axes and chiral centers in an elementary step, formulae have been developed for the calculation of the number of single events. The Evans-Polanyi relationship and the concept of stabilization energy were introduced to account for energy levels in surface-bonded carbenium ions. A novel functional dependency of the stabilization energy with the nature of the carbenium ion and the carbon number was proposed to account for energy effects from the acid sites on the catalyst. Further reductions in the number of parameters and simplification of the equations for the transient pseudohomogeneous one-dimensional plug-flow model of the reactor were achieved by means of thermodynamic constraints. Altogether, the single event concept, the Evans-Polanyi relationship, the stabilization energy approach and the thermodynamic constraints led to a set of 14 parameters necessary for a complete description of solid acid alkylation at the elementary step level.