Oxidative Coupling of Methane A Microkinetic Model Accounting.pdf

Dec 2, 2013 - tions in the early 80s5,6 but has still not reached industrial implementation. ... sively the oxidative transformations of light alkanes...
1 downloads 0 Views 2MB Size
Article pubs.acs.org/IECR

Oxidative Coupling of Methane: A Microkinetic Model Accounting for Intraparticle Surface-Intermediates Concentration Profiles Panagiotis N. Kechagiopoulos, Joris W. Thybaut,* and Guy B. Marin Ghent University, Laboratory for Chemical Technology, Technologiepark 914, B-9052 Ghent, Belgium S Supporting Information *

ABSTRACT: A microkinetic model for oxidative coupling of methane (OCM) has been developed that comprises a reaction network of 39 gas-phase and 26 catalytic elementary steps. It has been implemented in a heterogeneous reactor model that explicitly accounts for the interactions between gas phase and surface species. Concentration gradients arising from mass transport limitations are found to develop inside the catalyst pellet for all intermediates (i.e., surface and gas-phase) even under an intrinsic kinetics regime for the molecules and clearly affect the C2 selectivity. Special attention has been devoted to the reduction of the number of adjustable parameters in the model and the a priori determination of thermodynamic as well as kinetic parameters. A contribution analysis is conducted in order to elucidate the complex reaction pathways in OCM that lead to the desired products. Apart from the methyl radicals that couple to an extent of almost 70% in the void space between the pellets, the catalyst pellet accounts for the majority of molecules and radicals conversion, which are produced on the surface and further interact either in the catalyst pores or with other surface species. Almost 95% of CH4 consumption and more than half of the C2H6 production take place inside the catalyst pellet. A similar analysis is applied to understand the effect of various textural properties of catalysts on the performance of OCM, for example, increasing the catalyst porosity is found beneficial for the C2 yield, as long as a sufficient CH4 activation takes place, so that the coupling pathway is promoted over the heterogeneous oxidation of CH3·.

1. INTRODUCTION The oxidative coupling of methane (OCM) is a heterogeneously catalyzed gas-phase reaction.1 The catalyst produces methyl radicals, which can further react through gas-phase reactions in the catalyst pores as well as in the void space between the catalyst pellets.2,3 More specifically, methyl radicals can either couple to form ethane or can be oxidized toward carbon oxides via a branched-chain reaction mechanism.4 The reaction has been extensively investigated since the first relevant publications in the early 80s5,6 but has still not reached industrial implementation. An efficient way to extract knowledge from the large experimental data sets that have been generated over various OCM catalyst formulations is via a microkinetic analysis, which accounts for the gas-phase reactions, the fundamental catalytic surface chemistry, and the interaction between both.7 Incorporating such a microkinetic model into an advanced reactor model, which considers all relevant transport phenomena, can lead to a deeper understanding and optimization of the investigated reaction, especially in the case of OCM, given its recognized dual homogeneous−heterogeneous nature.8−10 Aparicio et al.11 were one of the first to conduct a microkinetic analysis of OCM kinetics data to gain a mechanistic insight in the occurring chemistry. Wolf12 and Buyevskaya et al.13 discriminated between various elementary step models that considered different heterogeneous reaction pathways for methane activation, surface oxidation by gas-phase oxygen and consecutive reaction steps of methyl radicals and ethane, supported by transient experiments on MgO and Sm2O3. More recently, Su et al.14 developed a thermodynamically consistent microkinetic model with which an upper bound of 28% for C2 yield was calculated. Similarly, Simon et al.15 applied various © 2013 American Chemical Society

physicochemical relations and thermodynamic constraints to calculate the parameters of a microkinetic model that was used to describe data over a La2O3 catalyst. Sinev16 has studied extensively the oxidative transformations of light alkanes. Early on, the interaction of radicals with the surface was investigated,17,18 wherein a Polanyi correlation was applied to determine the activation energies of the various proposed heterogeneous steps by postulating an analogy with equivalent homogeneous steps. The key role of surface hydroxyl groups was discussed, while a kinetic model accounting for surface interactions and diffusion of oxygen, hydrogen, and hydroxyl ions in the oxide lattice was developed to account for the participation of lattice oxygen in methane oxidation.19 More recently, the author’s findings were consolidated in an elaborate microkinetic model that was able to describe qualitatively various experimental observations.20,21 In the work of Couwenberg et al.,22 it was established that the participating gas-phase intermediates are highly reactive and lead to irreducible mass transport limitations for the latter. This effectively means that the pellet size necessary to consider the internal concentration gradients negligible is so small that it would lead to unacceptable pressure drop in the catalyst bed. As a result, during the determination of OCM kinetics these transport limitations have to be explicitly considered. The existence and importance of such gradients was also stressed, although not directly accounted for, in the work of Sinev.20 A one-dimensional heterogeneous reactor model was specifically Special Issue: Recent Advances in Natural Gas Conversion Received: Revised: Accepted: Published: 1825

September 23, 2013 November 28, 2013 December 2, 2013 December 2, 2013 dx.doi.org/10.1021/ie403160s | Ind. Eng. Chem. Res. 2014, 53, 1825−1840

Industrial & Engineering Chemistry Research

Article

developed by Couwenberg et al.22,23 for this purpose, that is, to adequately describe the interactions between the elementary catalytic reactions and the gas-phase reactions. The model was successfully applied to investigate the effects of the reactive intermediates on conversions, selectivities, and yields over Li/MgO and Sn/Li/MgO catalysts. Subsequently, Sun et al.24 expanded the surface reaction network by enhancing the detail in which deep oxidation reactions were represented and by the introduction of reversible reaction steps. Additionally, the application of catalyst descriptors was introduced in the model to facilitate knowledge extraction from high-throughput experiments. More recently, in the work of Thybaut et al.25 the surface network was further elaborated, including reaction steps for the description of the heterogeneous oxidation of ethylene. Moreover, a theoretical study was conducted so as to optimize C2 product yields varying various operating conditions and catalyst texture properties. In the present work, the aforementioned OCM model is further elaborated on the reactor as well as on the microkinetic level. Intraparticle concentration profiles for the surface intermediates are, for the first time, also accounted for explicitly, in addition to those already implemented for gas phase molecules and radicals. The implementation of the former is necessary given the reversibility of the reactions between gas phase and surface species and the importance of this interaction on products’ yield. In parallel, a considerable expansion of the surface network has been implemented to more adequately describe the complex OCM chemistry. Special attention is devoted to the a priori determination of the majority of surface network parameters, among others, by applying thermodynamic consistency within catalytic cycles. The effects of transport limitations are discussed, while an experimental data set derived from testing of a Sn/Li/MgO catalyst is used to validate the model. Via contribution analysis, the complex OCM reaction pathway is elaborated to gain a more detailed understanding of the relation between catalyst performance and properties.

The coordinate system applied in the continuity equations, in line with the work of Couwenberg et al.,22 is illustrated in more detail in Figure 1. On the reactor scale a single coordinate, z, going from 0 to L, is employed. At the pellet scale, two different coordinates are used: ξ for the intraparticle phase starting from the center of the catalyst pellet to the external surface and r for the interstitial phase starting from the center of the interstitial phase to the gas−solid interface. The respective continuity equations for all the components (i.e., molecules and radicals) in the reaction network that have to be integrated, together with the corresponding boundary conditions, are given as follows: Interstitial phase: Dm , ∂ ⎛ ∂Ci , g ⎞ Fv ∂Ci , g i = ⎜r ⎟ + R i,g εbAs ∂z r ∂r ⎝ ∂r ⎠

(1)

Intraparticle phase: −

De , i ∂ ⎛ 2 ∂Ci , c ⎞ ⎜ξ ⎟ = R i , c + εcR i , g ξ 2 ∂ξ ⎝ ∂ξ ⎠

(2)

Boundary conditions: z=0∧00∧r=

ξ=0

2. REACTOR MODEL 2.1. Continuity Equations. In order to adequately account for the irreducible transport limitations of radicals, a heterogeneous model that explicitly distinguishes between a solid phase and a fluid phase has been adopted. The solid phase contains pores in which reactants, products and intermediates diffuse toward and away from active sites and in which they can also react with each other. The two phases considered in the model (i.e., the porous solid and the fluid) will be denoted as the intraparticle and the interstitial phase respectively. Mass transfer by molecular diffusion from the interstitial phase to the gas−solid interface occurs simultaneously with gas-phase reactions (eq 1 below). Also, mass-transport by internal diffusion in the catalyst pores occurs simultaneously with chemical reactions, in this case on the catalyst pores walls as well as in the gas phase in the catalyst pores (eq 2 below). No internal heat transport or external heat transfer gradients are considered, since the present model is developed for the assessment of data derived from laboratory-scale, intrinsic kinetic measurements. This entails that during experimental design appropriate criteria should be evaluated to ensure the lack of such limitations.26 Nonetheless, Hoebink et al.27 studied the design of industrial scale fixed bed OCM reactors and found that it is specifically external heat transfer limitations that play a more important role, as compared to internal thermal gradients.

ξ=

dp 2

∂Ci , c ∂ξ

dv 2

∂r

dv 2

Ci , g = Ci ,0

=0

− ag Dm , i

=0

Ci , g = Ci , c

(3)

(4)

∂Ci , g ∂r

= acDe , i

∂Ci , c ∂ξ

(5)

(6)

(7)

where i ranges from 1 to Nmol + Nrad. The concentrations of the molecules at the inlet of the catalyst bed are determined by the boundary condition of eq 3, which expresses that no pellet scale gradient occurs at this position. The other boundary conditions are given by eqs 4−7. The gas−solid interfacial surface area, ag, and the external surface area of the catalyst, ac, as used in eq 4, are not equal since during kinetic measurements the catalyst bed is, typically, diluted with an inert material. dv, that is, the characteristic dimension of the interstitial phase, represents the average distance between two catalyst pellets. The value of dv can be calculated from the relation between the total gas-phase volume and the interfacial surface area, ag, in the reactor. Assuming that the reactor is filled with spherical particles of diameter dp and accounting for the internal porosity of the bed, the diameter dv of imaginary cylinders representing the gas phase volume in the reactor can be calculated as follows: dp ⎛ ε ⎞ dv = ⎜ b ⎟ 2 3 ⎝ 1 − εb ⎠

(8)

Moreover, if long cylindrical pores of uniform pore diameter are assumed for the catalyst, the average pore radius can be calculated according to 1826

dx.doi.org/10.1021/ie403160s | Ind. Eng. Chem. Res. 2014, 53, 1825−1840

Industrial & Engineering Chemistry Research

Article

Figure 1. Scheme of the heterogeneous reactor model accounting for mass transport limitations.22

r̅ =

2εp ρp S

calculation of the net formation rates of molecules and radicals, Ri,c, participating in the surface network that are required for eq 2. 2.3. Solution Procedure. The above-discussed set of second-order partial differential equations is integrated applying the method of orthogonal collocation.28,29 Orthogonal collocation belongs to the methods of weighted residuals and involves an approximation of the differential equations by a trial function. In this method, it is only necessary to evaluate the residuals at the collocation points. In the orthogonal collocation method these collocation points are chosen as the roots of orthogonal polynomials. The method is very suitable for the solution of two-point boundary-value problems.30 The trial function is then a set of orthogonal polynomials which satisfy the boundary conditions, while the roots to the polynomials determine the collocation points. The values of the solution at the collocation points can then be calculated from a set of algebraic equations.28,29 More specifically, the equations of the current problem can be solved by discretization in the r and the ξ direction using orthogonal collocation, and by subsequently applying a variable-step backward differentiation formula (BDF) integration method31 for the z direction. As such, the PDE set can be converted to a differential−algebraic one (DAE), with the number of model equations depending on the number of components, molecules, radicals, and surface species considered and on the number of collocation points chosen for the discretization of the interstitial and intraparticle dimensions. For the interstitial phase 4 collocation

(9)

2.2. Surface Intermediates. The balances for the surface species in the reaction network are an explicit part of the reactor model continuity equations and, as such, belong to the state variables of the mathematical set of equations presented above. In order to calculate the concentration of the surface intermediates, the pseudo-steady-state approximation is made, meaning that the net rates of formation of the surface species are considered equal to 0: R i,s = 0

(10)

where i ranges from 1 to Nsurf. Additionally, the catalyst active site balance is accounted for via the following equation: Nsurf

θ + *

∑ θi = 1 i=1

(11)

with * representing the free site. The above equations complete the set of partial differential equations (PDE) presented in the previous paragraph and are solved simultaneously at each point of the intraparticle dimension ξ together with the rest of the equations to calculate the surface coverage by the considered species. The specific calculation way of Ri,s is presented in paragraph 3.3, while as will be described in more detail in paragraph 3.4, the obtained surface concentration is further utilized for the 1827

dx.doi.org/10.1021/ie403160s | Ind. Eng. Chem. Res. 2014, 53, 1825−1840

Industrial & Engineering Chemistry Research

Article

reduction of the number of adjustable parameters and the a priori determination of thermodynamic and kinetic parameters. 3.1. Gas-Phase Reaction Network. The reaction network in the absence of catalyst has been adopted from the work of Chen et al.34 and consists of 39 reversible reactions among 13 molecules and 10 radicals as has been presented in the aforementioned work and more recently in the work by Sun et al.24 Typical gas-phase methane oxidation networks derived from combustion chemistry comprise several hundreds of elementary steps.35 The current network was selected by reducing such comprehensive networks36,37 after conducting contribution and sensitivity analyses based on experimental data at typical OCM conditions. Moreover, the presently used gas-phase reaction network was validated with catalytic experimental data in the work of Couwenberg et al.23 over both a Li/MgO and a Sn/Li/ MgO catalyst, where the contribution of the gas-phase network is high and low, respectively. The molecules considered are dihydrogen, water, hydrogen peroxide, dioxygen, methane, methanal, carbon monoxide, carbon dioxide, ethyne, ethene, ethane, propene, and propane. The radicals are hydrogen and oxygen atoms, hydroxyl, hydrogen peroxy, formyl, methoxy radicals, and methyl, vinyl, ethyl, and propyl radicals. Preexponential factors and activation energies, as reported by Chen et al.,34 are utilized to calculate forward rate coefficients via the Arrhenius equation, while the rate coefficients for the backward directions are calculated by means of equilibrium coefficients and the equivalent forward rate coefficients. The equilibrium coefficients are calculated at a given temperature by using the thermodynamic database of CHEMKIN,38 which is interfaced to the simulation software used. All parameters of the gas-phase network, pre-exponential factors, and activation energies are fixed during regression using catalytic experimental data. Finally, in order to adequately describe the total pressure effect, the pressure falloff behavior of the rate coefficients for the unimolecular reactions in the network has been taken into account explicitly. As described in detail by Chen et al.,39 the rate of these reactions shows a stronger pressure dependency compared to the other reactions, which is ascribed to the activation by collision with any other molecule (i.e., so-called third bodies). This falloff factor is, hence, a function of the concentration of these third bodies. 3.2. Catalytic Reaction Network. The elementary steps on the catalyst surface considered in the kinetic model are shown in Table 1. This set of reversible catalytic reactions is coupled to the gas-phase branched-chain reaction network described above via the reactor model equations presented in section 2.1. The complete kinetic model thus contains 26 catalytic and 39 gas-phase reactions in which 24 gas-phase, including inert N2, and 11 surface components, including the active site, are involved. The catalytic network describes methane activation on the catalyst surface by the dissociative oxygen adsorption (reaction C1), hydrogen abstraction of methane (reaction C2), and the regeneration of the active site (reactions C5 and C6). It accounts for CO2 generation through three reaction pathways: methyl radical scavenging followed by the sequential hydrogen abstraction from methoxy species on the catalyst (reactions C7−C13); the oxidation of CO adsorbed from the gas-phase (reactions C11−C13); and the heterogeneous oxidation of ethylene (reactions C14−C16). Given the average lattice constants of metal oxides typically used as OCM catalysts and the corresponding active site densities considered, second order surface reactions, such as the

points suffice, while the intraparticle dimension is discretized using 7 points. A spherical coordinate system for the latter and a cylindrical coordinate system for the former phase are used, allowing the use of symmetrical polynomials as trial functions28,29 and, hence, automatically fulfilling the boundary conditions at r = 0 (eq 4) and ξ = 0 (eq 6). As such, the aforementioned number of collocation points represents the internal ones plus the right boundary. Twenty-four components are considered in total, fourteen of them being molecules, including inert N2, and the remaining ten being radicals. Discretization of continuity eq 1, hence, leads to 3 × 24 differential equations, while 24 algebraic ones result from eq 5 at the right boundary. Equation 2 is converted to 6 × 24 algebraic equations, with 24 additional algebraic ones defined by eq 7 that represents again the right boundary. Moreover, the balances of the 10 surface intermediates are considered for each intraparticle collocation point, meaning that eqs 10 and 11 contribute with 7 × 10 and 7 additional algebraic equations, correspondingly. In total, 341 state variables, 72 of them being differential ones, have to be integrated along the reactor axis with initial conditions as defined by eq 3. For the radicals, the initial concentration is assumed equal to zero, while for the surface species, being algebraic variables, reasonable estimates are provided to initiate the integration. The molecular binary diffusion coefficient of component i in j, Di,j, is calculated using the relation of Füller−Schettler− Giddings:32 ⎛ 1 ⎞1/2 T1.75 1 ⎜⎜ ⎟ Di , j = 3.16 × 10 + Mj ⎟⎠ Ptot((Σv)1/3 + (Σv)1/3 i j ) ⎝ Mi −8

(12)

The bulk diffusivity of component i in a gas mixture containing Nmol + Nrad components, Dm,i, can be estimated from the individual binary diffusivities, Di,j, according to the Wilke equation:32 Dm , i

⎛ Nmol + Nrad x ⎞−1 j ⎟ = (1 − xi)⎜⎜ ∑ ⎟ ⎝ j = 1(j ≠ i) Di , j ⎠

(13)

Finally, taking also into account the internal void fraction of the porous material and the tortuous nature of the pores,33 the effective diffusivity for component i, De,i, can be calculated: ε De , i = c Dm , i τc (14) Equations 13 and 14 are applied at each collocation point of the interstitial and intraparticle coordinates accordingly, given that molecular diffusivities are dependent on the local gas concentration.

3. KINETIC MODEL Oxidative coupling of methane occurs through a mechanism in which catalytic reactions interact with gas-phase chain reactions. In the present work, a kinetic model for oxidative coupling of methane which is valid over a wide range of operating conditions is implemented. The network contains 26 catalytic elementary reactions coupled with a gas-phase reaction network consisting of 39 elementary reaction steps. In what follows the kinetic network is described in more detail, while the Supporting Information elaborates on the methodology applied for the 1828

dx.doi.org/10.1021/ie403160s | Ind. Eng. Chem. Res. 2014, 53, 1825−1840

Industrial & Engineering Chemistry Research

Article

Table 1. Catalytic Elementary Reactions Considered in the Oxidative Coupling of Methane O2 + 2* ⇔ 2O*

C2H4 + O* ⇔ C2H4O*

(C1)

CH4 + O* ⇔ CH3·+ OH*

C2H4O* + O* ⇔ C2H3O* + OH*

(C2)

C2H4 + O* ⇔ C2H3·+OH*

(C3)

C2H6 + O* ⇔ C2H5·+ OH*

(C4)

2OH* ⇔ H 2O* + O* H 2O* ⇔ H 2O + *

(C7)

CH 2O* + O* ⇔ HCO* + OH*

(C9) (C10)

HCO* + O* ⇔ CO* + OH*

CO* + O* ⇔ CO2 * + *

(C17) (C18) (C19) (C20)

(C21)

H 2O2 + O* ⇔ HO2 ·+ OH* OH ·+ O* ⇔ O·+ OH*

(C22) (C23)

H 2O + O* ⇔ OH ·+ OH*

(C11)

(C12)

HO2 ·+ O* ⇔ O2 + OH* HO2 ·+* ⇔ OH ·+ O*

(C13)

(C24) (C25)

(C26)

with Ngas being the number of gas-phase reactions. This equation is used both for the interstitial phase and for the contribution of the gas phase network in the intraparticle phase, as seen in eqs 1 and 2 respectively. The rate of an elementary step in the gas phase network is calculated through the law of mass action according to Nmol + Nrad

rj = k j



Ci vj ,i

i=1

(16)

where kj is the rate coefficient of reaction j. Equivalently, the rate of an elementary step in the catalytic reactions follows from the law of mass action. For example, the rate of catalytic reaction C2 is calculated according to rC2 = k C2CCH4θO Lt *

(17)

All of catalytic reaction rates are calculated per internal surface area of the catalyst pellet and have to be multiplied with the specific surface area, Sg, and the catalyst density, ρc, so that they can be expressed per pellet volume. They are used for calculating the net production rates of surface species, Ri,s, according to

surf

quenching step 4HO2· ⎯⎯⎯→ 3O2 + 2H2O, which was previously included in the surface network.22−25 Further possible surface steps that could be considered relate, for example, to hydrocarbon surface species reacting also with gas phase molecules and radicals in an Eley−Rideal way or oxygen containing radicals directly interacting with the bare active site.21,41 Such possibilities, considering the increase in the adjustable parameters of the model they would incur and the additional complexity in relation to the mass transport coupled to microkinetics approach followed, were chosen not to be included in the current surface network. The successful description of experimental data over three different catalysts, presented elsewhere,42 indicates that the most prominent surface reactions have been accounted for in the current network. 3.3. Calculation of Reaction Rates. The net production rates of all components in the gas-phase and catalytic reactions have to be calculated in order to integrate the reactor equations presented in section 2.1. The net production rate of a component i in the gas-phase reaction network, Ri,g, is calculated through

C 26

R i,s =

∑ j = C1

vj , irj (18)

As mentioned, given the pseudosteady state approximation applied, net rates, Ri,s, are used for the determination of the surface coverages, θi, via eqs 10 and 11. However, in order to calculate the catalytic net production rate of molecules and radicals, Ri,c, used in eq 2, a separate procedure is followed, see section 3.4. 3.4. Global Catalytic Reaction Paths. Sixteen global reactions paths (a to p), built up from catalytic elementary steps, can be distinguished in this model and are presented in Table 2. These paths represent the number of independent catalytic cycles that exist in the considered surface network. Path a, for example, represents the catalytic formation of methyl radicals from methane. It can be seen that this path can be constructed from catalytic elementary steps 1, 2, 5, and 6 taking into account appropriate stoichiometry numbers, that is, 1, 4, 2, and 2. Moreover, after constructing in a similar way the combinations for all independent paths, it can be seen that

N gas

∑ vj ,irj j=1

(C16)

CH 2O + O* ⇔ CHO·+ OH* H 2 + O* ⇔ H ·+ OH*

latter, are expected to potentially significantly contribute to the reaction behavior and, hence, were considered in the model. Based on the observations by Sinev et al.,20,21 any gas phase species could potentially follow similar H-atom abstraction reactions, much as methane does in reaction C2. As such, the OCM surface network, as was presented in the work of Thybaut et al.25 and already including such reactions for ethane C3 and ethylene C4, was expanded by the inclusion of reactions C17−C25 to account for similar collisional interactions. All gasphase species were considered, apart from C3 compounds, whose concentration is negligible at OCM conditions. Reaction C25 accounts also for the quenching of hydroperoxy radicals on the catalyst surface that is considered crucial in OCM performance, since the radical is one of the most important chain carriers in the gas-phase reaction network.4,14,35,40 Similarly to the steps included in the work of Su et al.,14 reaction C26 further contributes to the quenching of the hydroperoxy radical via cleavage of the weak peroxide bond, while somehow acting also as a chain-branching step.20 Steps C25 and C26 replaced the pseudohomogeneous

R i,g =

(C15)

CH3O·+ O* ⇔ CH 2O + OH*

CHO·+ O* ⇔ CO + OH*

(C8)

CO2 + * ⇔ CO2 *

C2H5·+ O* ⇔ C2H4 + OH*

(C5)

CH3O* + O* ⇔ CH 2O* + OH*

CO + * ⇔ CO*

C2H3O* + O* ⇔ CH 2O* + HCO*

(C6)

CH3·+ O* ⇔ CH3O*

(C14)

(15) 1829

dx.doi.org/10.1021/ie403160s | Ind. Eng. Chem. Res. 2014, 53, 1825−1840

Industrial & Engineering Chemistry Research

Article

Table 2. Global Catalytic Reaction Paths Constructed Based on the Elementary Catalytic Network in Table 1 and the Respective Stoichiometric Coefficients in Table S1 of the Supporting Information a.

4CH4 + O2 ⇔ 4CH3·+ 2H 2O

i.

4CH 2O + O2 ⇔ 4CHO·+ 2H 2O

b.

4C2H4 + O2 ⇔ 4C2H3·+ 2H 2O

j.

4CHO·+O2 ⇔ 4CO + 2H 2O

c.

4C2H6 + O2 ⇔ 4C2H5·+ 2H 2O

k.

4H 2 + O2 ⇔ 4H ·+ 2H 2O

d.

2CO + O2 ⇔ 2CO2

l.

4H 2O2 + O2 ⇔ 4HO2 ·+ 2H 2O

e.

4CH3·+ 7O2 ⇔ 4CO2 + 6H 2O

m.

4HO2 · ⇔ 3O2 + 2H 2O

f.

C2H4 + 3O2 ⇔ 2CO2 + 2H 2O

n.

2OH·+ 2O2 ⇔ 2HO2 ·

g.

4C2H5·+ O2 ⇔ C2H4 + 2H 2O

o.

2H 2O + O2 ⇔ 4OH·

h.

4CH3O·+ O2 ⇔ 4CH 2O + 2H 2O

p.

4OH ·+ O2 ⇔ 4O·+ 2H 2O

The implementation of thermodynamic relationships in a microkinetic model is one approach in relating model parameters to catalyst descriptors and simultaneously reducing the number of adjustable parameters. Thermodynamic relationships can be established between surface reaction enthalpies and analogous gas-phase reactions provided that the adsorption thermodynamics of the species involved can be adequately described or can be defined as adjustable parameter.47 Applying thermodynamic consistency through such thermodynamic cycles allows to express reaction enthalpies for various elementary steps as a function of a limited number of unknown parameters, which will then be obtained, for example, by regression. Similar approaches are widely used during the development of microkinetic models,48,49 which have successfully described experimental data from a wide variety of applications.50,51 For the case of OCM, as already mentioned, the work of Su et al.14 was among the first to ensure thermodynamic consistency across their considered surface network. In the Supporting Information, this methodology is elaborated based on the OCM microkinetic network applied in this work, so as to calculate the surface reaction enthalpies and corresponding activation energies of all surface steps. On the other hand, collision theory is used to calculate the maximum value of the pre-exponential factors of all steps involving the collision of a molecule on the catalyst surface. These reactions are actually all of the adsorption steps, including the forward direction of reactions C7 and C14 that describe the capture of CH3· and C2H4 on the surface, and all of the Eley− Rideal steps. In the current work, these maximum values for the adsorption steps were corrected by a factor, denoted as sticking coefficient. These factors represent the ratio of the number of the respective molecules or radicals actually adsorbing to a clean surface to the total number of them colliding with it. Since they vary depending on catalyst, they were all defined as catalyst descriptors, meaning they are part of the adjustable parameters of the model, as also seen in section 4.1. The collision frequencies were calculated according to the following equation:52

reaction C2 participates four times in global path a, while it is not involved in any other global path. Thus, the rate of global reaction a is equal to ra = rC2/4. This way of presenting a complex reaction network follows from the work of Temkin43 and allows a straightforward calculation of the net production rates of the molecules and radicals taking part in the surface network.44 Similarly, paths b and c account for the catalytic formation of vinyl and ethyl radicals from ethylene and ethane, respectively. Path d shows the catalytic conversion of carbon monoxide into carbon dioxide, and path e describes the deep oxidation of methyl radicals into carbon dioxide. Path f represents the heterogeneous oxidation of ethylene. Likewise, the other paths describe H-atom abstraction of the remaining gas-phase species, except for path m, which accounts for the scavenging of hydroperoxy radicals on the catalyst surface. Similar analogies to path a can be found for all of the other global paths according to the coefficients presented in Table S1 of the Supporting Information, resulting in the following rates: ra = rC2/4,

rb = rC3/4,

rc = rC4 /4,

rd = rC7/4,

re = rC12/2,

rf = rC14 ,

rg = rC17/4,

rh = rC18/4,

ri = rC19/4,

rj = rC20/4,

rm = rC23/4,

rn = − 2rC24 ,

rk = rC21/4, ro = rC25/4,

rl = rC22/4, rp = rC19/4 (19)

It are these rates, of the global reaction paths, that are utilized in the reactor continuity eq 2 for the intraparticle phase according to p

R i,c =

∑ vj ,irj j=a

(20)

3.5. Calculation of Surface Reaction Network Parameters. Microkinetic model parameters are classified into kinetic and catalyst descriptors. While the former are exclusively related to the reaction kinetics, independent of the catalyst used, the latter specifically account for the physical and chemical catalyst properties and represent the essential differences between the various candidate catalysts.45,46 Special attention has been devoted to incorporating catalyst descriptors into the microkinetic model for methane oxidative coupling, as described in Sun et al.24 The main advantage of this approach is that the developed model is capable of quantifying trends between catalyst activity and selectivity on the one hand and properties on the other hand. As a result, such a model allows addressing which catalyst descriptors and, hence, corresponding properties are the most relevant for the improvement of catalytic performance.

Ai =

1 Lt n

RT 2πM

(21)

The remaining pre-exponential factor values, pertaining to reactions involving only surface species, have been left unchanged at their values as previously reported,24,25 which were calculated in those references via statistical thermodynamics and transition state theory. In order to ensure thermodynamic consistency, as seen in the Supporting Information, the surface reaction network related 1830

dx.doi.org/10.1021/ie403160s | Ind. Eng. Chem. Res. 2014, 53, 1825−1840

Industrial & Engineering Chemistry Research

Article

compared to the previous versions of the model, kinetic, as well as catalyst, descriptors have been estimated to describe this data set. In Table 3, the estimated catalyst descriptors are given for

calculations require the knowledge of thermodynamic properties, namely of gas-phase reaction enthalpies. Additionally, in order to calculate the equilibrium coefficients of the gas-phase network, the Gibbs energy difference is needed, which in turn requires the gas-phase species molar entropies. It is assumed that the standard-state thermodynamic properties of all species are functions of temperature only, expressed in the form of polynomial fits. In particular the format used in the NASA chemical equilibrium code53 is applied, where seven coefficients are needed for each of two temperature ranges, above and below 1000 K: α α α α α H° = α1 + 2 T + 3 T 2 + 4 T 3 + 5 T 4 + 6 RT 2 3 4 5 T

Table 3. Catalyst Descriptors Used in the Simulation of the Data Obtained on Sn/Li/MgO42 catalyst descriptor D1 D2 D3 D4 D5 D6 D7 D8 D9 D10 D11 D12 D13 D14 D15 D16

(22)

α α α S° = α1ln T + α2T + 3 T 2 + 4 T 3 + 5 T 4 + α7 R 2 3 4 (23)

Finally, the standard-state Gibb’s free energy is calculated according to G° = H ° − TS°

(24)

The actual coefficients for the above equations for all the used species were taken from the database of the CHEMKIN software from Reaction Design.38

4. RESULTS AND DISCUSSION 4.1. Model Validation. The developed model was validated against an experimental data set on a Sn promoted Li/MgO catalyst. This data set has been elaborately discussed in previous publications22−24 and will, therefore, not be further commented in the current work. In Figure 2, the simulated molar fractions

reaction enthalpy of H-atom abstraction of CH4 chemisorption enthalpy of O2 chemisorption enthalpy of CH2O chemisorption enthalpy of CHO· chemisorption enthalpy of CO chemisorption enthalpy of CO2 chemisorption enthalpy of H2O chemisorption enthalpy of C2H4O chemisorption enthalpy of C2H3O· initial sticking probability of O2 initial sticking probability of CH3· initial sticking probability of CO initial sticking probability of CO2 initial sticking probability of H2O initial sticking probability of C2H4 density of active sites

unit

value −1

(kJ mol ) (kJ (kJ (kJ (kJ (kJ (kJ (kJ (kJ

mol−1) mol−1) mol−1) mol−1) mol−1) mol−1) mol−1) mol−1)

−2

(mol m )

56.63 55.37 123.13 141.08 74.42 87.00 34.81 42.35 92.59 0.56 6.23 × 5.66 × 0.0154 0.0765 5.48 × 1.33 ×

10−05 10−05

10−05 10−06

Table 4. Reaction Families and Corresponding Parameters Used in the Polanyi Relationships Considered in Methane Oxidative Coupling42 reaction family rf1 rf2 rf3 rf4 rf5

H abstraction by Eley−Rideal reaction (steps C2−C4, C17−C25) H abstraction by surface reaction (steps C8−C10, C15) Recombination of OH (step C5) CO oxidation on the surface (step C12) C−C cleavage on the surface (step C16)

E0 (kJ mol−1)

α 0.75

58

96.8

0.50

52

141.3

0.6559 0.2660 0.9760

130.2 70.9 183.2

the Sn/Li/MgO catalyst, while Table 4 presents the values of Polanyi parameters for the considered reactions families.42 More details on these parameters in relation to the surface network parameters are found in the Supporting Information and the previous publications.24,25 Moreover, in the work of Sun et al.,24 a sensitivity analysis was presented to identify the most important descriptors in the microkinetic model. Considering the proposed methodology in studying the OCM reaction, wherein microkinetics are coupled with mass transport phenomena, the validation approach followed can be summarized as follows: The initial model developed by Couwenberg et al.22 considered a nonreversible surface network among a limited number of surface species. As such, an analytical solution of the set of equations describing the surface coverages was possible. Moreover, no intraparticle gradients were considered for the surface species in that work. The expansion of the surface reaction network as well as of the number of participating surface species and inclusion of reaction reversibility by Sun et al.24 necessitated the numerical determination of this coverage. The focus of that work was on the implementation of catalyst descriptors in the model and their validation with the Couwenberg data. Reaction reversibility, however, given the irreducible mass transport limitations of the radicals, should affect also surface coverage and, hence, OCM kinetics.

Figure 2. Parity diagram for the product outlet molar percentages on Sn/Li/MgO.

are compared to the experimental ones, where it is seen that an overall satisfactory description of the data has been achieved with the proposed model. As described in detail previously,24,25 by implementing catalyst descriptors the microkinetic model can be applied on a library of catalysts rather than a single catalyst. Such an assessment of different catalysts has also been conducted and is presented elsewhere.42 Given the significant expansion of the surface kinetic network in the current work 1831

dx.doi.org/10.1021/ie403160s | Ind. Eng. Chem. Res. 2014, 53, 1825−1840

Industrial & Engineering Chemistry Research

Article

pronounced diffusion limitations, which contribute to the (less pronounced) gradients observed for CO2. The above observations are indeed reflected by the calculated concentration profiles at the pellet-scale for some characteristic molecules at the selected operating conditions, see Figure 3. The situation is drastically different for the reactive intermediates. The diffusion lengths, λ, of all radicals are, on average, 2−3 orders of magnitude lower than the radius of the catalyst pellet; see Table S2 of the Supporting Information. The very high reactivity of these compounds is easily observed from their lifetime, τ, which is drastically smaller compared to that of the molecules and in the order of μs or lower. Hence, strong concentration profiles develop that directly affect selectivities and product distribution in the exit of the catalyst bed. CH3· radicals, being surface produced intermediates, have a much higher concentration in the catalyst pores as compared to the interstitial gas phase. However, approaching the catalyst surface, this concentration decreases significantly due to a net diffusion toward the interstitial phase. On the contrary, a surface terminated species, such as the HO2· radical, presents a much higher concentration in the interstitial phase. The radical quenching mechanism is described here by including specific elementary steps (i.e., reactions C25 and C26), in contrast to the work of Couwenberg et al.,23 where only a pseudohomogeneous reaction was considered. Nonetheless, the shape of the calculated profiles for HO2· is in line with that work confirming the existence and important role of heterogeneous termination of gas-phase reactions. The above considerations indicate the strategic advantages of applying an elaborate reactor model, such as the one used in the present work, during kinetic modeling of experimental data aimed at the measurement of intrinsic OCM kinetics 4.2.2. Surface Intermediates. Given the complex interplay between gas-phase and catalytic reactions during OCM, it can be speculated that the irreducible mass transport limitations for the reactive intermediates discussed in the previous section also affect catalyst surface coverages. This important phenomenon has not been systematically investigated till now in literature. As described by Sun et al.24 and also presented in detail in the Supporting Information of the current work, the main catalytic cycle in OCM leading to CH3· radicals consists of four steps. Initially, the dissociative adsorption of O2 leads to the production of the oxidized site O*, followed by H-atom abstraction from CH4 by this site leading to CH3· radicals and surface hydroxyl species, OH*. Two such species recombine to form O* and H2O*, the latter desorbing in the gas phase and, thus, liberating again an active site. The majority of the calculated surface coverages is determined by this catalytic cycle. At the typically high OCM operating temperatures, an equilibrium has been observed among O2 in the gas phase and the surface oxygen species.24,54,55 This holds also for the current simulations, where reaction C1 is found to be quasiequilibrated. As evident from the previous section, O2 does not develop significant concentration gradients in the catalyst pores at the investigated operating conditions and, hence, also the surface coverage by active oxygen species O*, does not exhibit significant evolutions as a function of the particle diameter, see Figure 4. Along the axial reactor coordinate there is a depletion of gaseous oxygen due to its conversion by the OCM reaction. This depletion is also reflected by the O* coverage. The OH* species possesses the highest intraparticle surface coverage together with O*. The concentration of these two surface species is closely connected via the catalytic cycle

Hence, it is aim of the current work to specifically address this latter crucial aspect. In what follows, the insight acquired by the application of the current elaborate reactor and kinetic model will be focused on and the importance of considering such pellet scale phenomena will be quantified. Moreover, a detailed reaction pathway analysis is conducted to elucidate the relative importance of the homogeneous and heterogeneous network. 4.2. Concentrations Profiles. In the next paragraph, simulation results are presented to illustrate the existence of intraparticle and interstitial concentration gradients for the reactive intermediates. The effect and importance of surface coverage gradients in the catalyst pellet pores, which represents the main reactor model improvement of the current work, is addressed in the following one. The operating conditions at which this analysis has been performed are summarized in Table 5. Table 5. Experimental Conditions Used for the Base Case Simulations 1.13 × 105 Pa 3.6 mol mol−1 970 K 5.82 kg s mol−1 10.32% 36.66%

Ptot CH4/O2,0 T W/Ft,0 XCH4 XO2

dp ρc Sg εc τc εb

2.5 × 10−4 m 2300 kg m−3 2800 m2 kg−1 0.27 m3 m−3 5 m2 m−2 0.45 m3 m−3

4.2.1. Molecules and Radicals. The reactivity of a component can be quantified by examining its lifetime, τ, defined as the reciprocal of the pseudo-first-order rate coefficient, k1, while the diffusion rate of the component can be compared to its consumption rate via the diffusion length, λ, that is calculated according to

λi =

De , i k1

(25)

All of the quantities involved in these calculations are only available after integration of the set of differential equations presented above, since they depend on the local concentration of the components and their local consumption rates. Looking at the reactants and products, it can be seen that for most of them no significant concentration gradients exist in the pellet coordinate. As shown in more detail in Table S2 of the Supporting Information, the diffusion lengths obtained for reactants and products exceed the pellet radius; for example, at the selected operating conditions, the diffusion length of CH4 is more than seven times larger than the pellet radius (8.96 × 10−4 m), while for the case of O2 the diffusion length is, logically, slightly lower (1.83 × 10−4 m) but still sufficiently larger than the pellet radius avoiding the development of significant intraparticle gradients. For the products, all C2 and C3 compounds have satisfactorily large diffusion lengths (>8.00 × 10−4 m). The same holds for all noncarbon containing compounds that are clearly not limited by mass transport inside the catalyst pellet pores. Carbon oxides, on the contrary, and formaldehyde have a diffusion length inferior to the pellet radius (1.60−6.15 × 10−5 m) and, hence, develop concentration gradients inside the catalyst pellet, that are not too pronounced, however. Via the contribution analysis presented in paragraph 4.3, it can be observed that CO2 is exclusively produced in the gas phase via pathways involving radicals HO2· and CHO·. The latter exhibit very 1832

dx.doi.org/10.1021/ie403160s | Ind. Eng. Chem. Res. 2014, 53, 1825−1840

Industrial & Engineering Chemistry Research

Article

Figure 3. Concentration profiles at pellet and reactor scale for characteristic molecules and radicals at the conditions of Table 5. The origin of the pellet-coordinate axis corresponds to the center of the catalyst pellet, while the external surface of the catalyst pellet is located at 1.25 × 10−4 m. The remaining distance till 1.95 × 10−4 m corresponds to the interstitial gas phase. (a) C2H4, (b) CH2O, (c) CH3·, (d) HO2·.

described above. As discussed by Sinev et al.,21 efficient OCM catalysts exhibit strong oxygen binding leading to a pronounced occurrence of the reverse OH* recombination reaction C5, which directly promotes the coverage by hydroxyl groups at the high operating temperatures in OCM. Across the pellet scale coordinate, a clear enhancement of the OH* is visible in Figure 4. The decreasing concentration of CH3· along this direction, as discussed above, leads to a decreasing participation of the OH* consuming reverse reaction C2 that is responsible for its calculated profile. Regarding the axial coordinate, it should be noted that the major OCM reaction product C2H6 reacts, similarly to CH4, with surface oxygen in an Eley−Rideal type of reaction, leading to C2H5· radicals. The C−H bond strength is weaker in ethane compared to methane, enhancing the generation of OH* surface intermediates as ethane gets produced. Additionally, as was shown by Chen et al.,34 the gradual increase in ethane concentration leads to higher propagation rates in the branched-chain gas-phase reaction mechanism, further increasing the total radical concentration and indirectly the OH* on the surface. As discussed in paragraph 3.2, all gas phase species are able to undergo hydrogen abstraction reactions on the same sites responsible also for CH3· radicals’ generation. Given that the concentration of these species increases in the axial direction of the reactor, the coverage in OH* is further enhanced, explaining its calculated profile in that direction. Investigating the surface coverage of species related to the heterogeneous deep oxidation of CH3· radicals, it is clear that their shape is also strongly affected by the diffusion limitations

of the radical. This is valid for all surface species that participate in the degradation of the methoxy species up to CHO*, with the much higher concentration of the CH3· radicals toward the center of the pellet clearly being reflected in the surface species coverage. Concerning the overall concentration of these species, the very low coverage by CH3O* is qualitatively in line with the IR studies of Lacombe et al.,55 who observed only very weak bands related to methoxy species. Moreover, the considerably higher percentage of the surface being covered by CH2O* compared to the rest of the species related to deep oxidation is explained by reaction C9 being much less exothermic to the rest of the relevant reactions. As described in the Supporting Information, thermodynamic consistency is applied during the calculation of the reaction enthalpies, the much higher formation enthalpy of the equivalent CHO· radical in the thermodynamic cycle leading to a mild exothermicity of reaction C9. CO2 in the gas phase is in equilibrium with its adsorbed state, explaining its steady rise along the axial coordinate of the catalyst bed. At the studied conditions, reaction C13 is clearly quasi-equilibrated. In the pellet coordinate, the aforementioned mild intraparticle gradients of CO2 are not evident in the coverage of CO2*, due to the much larger change of the variable’s value in the axial direction. Nonetheless, there are differences as high as 3 orders of magnitude from the center of the pellet to its surface in the coverage by CO2*, especially at the start of the catalyst bed. The profile of H2O* is very similar to CO2* for equivalent reasons. It was already discussed that all C2 products have large diffusion lengths and, hence, do not 1833

dx.doi.org/10.1021/ie403160s | Ind. Eng. Chem. Res. 2014, 53, 1825−1840

Industrial & Engineering Chemistry Research

Article

Figure 4. Concentration profiles at pellet and reactor scale for surface species at the conditions of Table 5. A bottom contour projection is used for O* and OH* to enhance clarity. The origin of the pellet-coordinate axis corresponds to the center of the catalyst pellet, while the external surface of the catalyst pellet is located at 1.25 × 10−4 m. (a) O*, (b) OH*, (c) CH3O*, (d) CH2O*, (e) CO2*, (f) C2H4O*.

relevant effect during kinetic modeling. The latter statement gains even more importance considering that the goal of such a microkinetic model is to facilitate information extraction from high-throughput experimentation, so as to guide catalyst design, meaning it should be able to guarantee a wide applicability versus operating conditions and catalyst properties. 4.3. Contribution Analysis. A contribution analysis, in combination with information from the surface and gas-phase species profiles, can elucidate the importance of the various reaction pathways, providing valuable insight in the occurring chemistry in relation to process and/or catalyst parameters. Such an analysis is performed in the current work at the outlet of the catalyst bed. The performed contribution analysis is differential with respect to the axial coordinate, while, due to the existence of intraparticle and interstitial gradients, it is

develop significant pellet scale concentration gradients. As such, the surface species related to C2H4 oxidation, C2H4O* and C2H3O*, are also not affected by the diffusion limitations of the radicals, while their coverage increases along the axis due to the increasing concentration of C2H4 in the gas phase. The effect of these observed profiles on calculated selectivity is significant. The production rate of CH3· is uniform along the pellet, given the steady O* and CH4 intraparticle profiles. The consumption of the radical, however, is clearly affected, not only by the decreasing CH3· concentration, which would result mainly in a decreasing rate but also from the intraparticle profile of OH*. Not considering the intraparticle profiles for the surface intermediates at the investigated operating conditions would result in a drop of C2 selectivity from 58.2% to 50.1%, clearly showing the importance of considering the 1834

dx.doi.org/10.1021/ie403160s | Ind. Eng. Chem. Res. 2014, 53, 1825−1840

Industrial & Engineering Chemistry Research

Article

Figure 5. Contribution analysis for the intraparticle (left panel) and interstitial phase (right panel) at the end of the catalyst bed for the conditions presented in Table 5. Net production rates of carbon containing compounds are considered. Color mapping is used to indicate a higher conversion percentage toward the respective compound, with blue representing the lowest considered value (0.0%) and red the highest value for conversion (100%). The thickness of the arrows is relative to the actual net rate per reactor volume.

integral with respect to the pellet coordinate. All presented rates and relative contributions over phases are based on the “gas-phase” reactor volume, that is, referring to the gaseous space available in the interstitial phase as well as in the intraparticle phase. The results of this analysis for the two considered phases are summarized graphically in Figure 5, where the net formation rates of all involved molecules and radicals are presented for the two considered phases (i.e., the intraparticle and interstitial phase). The complexity of the reaction pathways during OCM is more than obvious for both phases. Moreover, the differences in the preferred pathways between the two phases should be noted, noting that it is the intraparticle phase that accounts for the majority of the methane and oxygen consumption and production of the most crucial compounds. The actual magnitude of a reaction rate is represented in this figure via the width of the arrows, such that the relative importance of the

two phases can be assessed by comparing the two panels. For example, it is obvious that CH3· production takes place in the intraparticle phase, while their coupling proceeds in both phases. Also, the importance of direct CO2 formation routes in the intraparticle phase, either from CH3· or C2H4, is evident. A more elaborate discussion based on the actual production and consumption rates is given below. In the interstitial phase, the major methane consumption routes relate to its reaction with OH· and H· radicals at a percentage equal to 62.5 and 25.8%, respectively, both leading to CH3· radicals. In contrast, in the intraparticle phase, hydrogen abstraction via reaction C2 dominates, accounting for 97.8% of the consumption of methane, also leading to the production of CH3· radicals. In total, the intraparticle phase contributes to more than 95% of the overall methane consumption. The produced CH3· radicals are consumed via a multitude of reactions in the interstitial phase, the most important one being their coupling at a 1835

dx.doi.org/10.1021/ie403160s | Ind. Eng. Chem. Res. 2014, 53, 1825−1840

Industrial & Engineering Chemistry Research

Article

In the interstitial phase 62.6% of the CO is consumed with HO2· toward the production of CO2, while the rest is converted to CHO· radicals. On the contrary, in the intraparticle phase, CO reacts almost exclusively toward CO2 via adsorption and catalytic oxidation. CO2 in the interstitial phase is produced only via CO oxidation, while in the intraparticle phase, as already commented, the catalytic oxidation of C2H4 and of CH3· radicals also contributes. However, the high concentration of CO at the end of the catalyst bed as compared to that of C2H4 and even more of CH3·, results in a clear domination of the CO oxidation pathway at a percentage of 94.8%. As already discussed, the catalyst surface contributes significantly to the production of the CH3· radicals. The same holds clearly also for most of the other radicals, where it is seen that in the intraparticle phase the respective H-atom abstraction reactions account for a major part of their production, namely 87.7% for C2H3· via reaction C3, 92.3% for C2H5· via reaction C4 and 48.8% for CHO· via reaction C19. On the other hand, the surface is also responsible for radical quenching. HO2·, which, as already discussed, facilitates to a large extent gas phase deep oxidation, is consumed solely via heterogeneous reactions in the intraparticle phase and more specifically via reaction C26 at a 99.8%. Given that the experimental conditions are typical for OCM and since O2 is still far from depleted at the end of the catalyst bed, it is believed that the analysis is representative for OCM kinetics of the studied catalyst at any axial position along the reactor. A similar analysis has been conducted at the beginning of the catalyst bed to evaluate the effect of the gas mixture composition. Noted difference is the pronounced role of the homogeneous reaction of CH4 and O2 toward the production CH3· and HO2· radicals, which, given its high activation barrier compared to the reaction of CH4 with radicals, is important only at the onset of the OCM reaction.37 An equally important difference is that CH3· radicals are consumed in the interstitial phase almost exclusively by coupling (94.2%), since numerous OCM products that react with CH3·, namely H2, H2O, HO2·, CH2O, C2H6, and C2H4, do not have yet considerable concentration. Finally, at the inlet, the primary route for catalytic CO2 production is obvious with 36.5% being produced in the intraparticle phase from CH3· total oxidation, in line with former observations.56 Nonetheless, the main consumption and production routes described previously were found to be valid also at this position. 4.4. Effect of Textural Properties. Thybaut et al.25 have discussed the global effect of various textural properties on catalyst performance. In the same respect, a sensitivity analysis would provide additional understanding on their effect but goes beyond the scope of the present manuscript. We have opted to discuss the effects of two important catalyst texture properties, porosity and surface area, in detail in the following two paragraphs by means of a contribution analysis. These simulations are performed by modifying the value of the investigated texture property, while keeping the other properties reported in Table 5 constant. 4.4.1. Pellet Porosity. As reported in Table 6, increasing porosity gives rise to an increase in methane conversion and selectivity toward C2 products.25 On average, the participation of the interstitial phase related to CH4 consumption remains largely unaffected and, as also commented in the previous section, averages 5% of the total CH4 consumption. However, looking at the contribution of the various complex pathways in the intraparticle phase, it can be seen that, as porosity increases,

percentage equal to 69.6%, leading to the desired product C2H6. Other competitive CH3· consumption pathways are their reaction with HO2· (5.2%), CH2O (5.6%), and C2H6 (7.8%). In the intraparticle phase, coupling in the catalyst pores is still considerable, accounting for 21.9% of CH3· consumption. Reaction C7 leading to the formation of surface methoxy species and eventually to undesired COx consumes 12.3% of the radicals, while it is the reverse of reaction C2 that is responsible for the highest percentage of CH3· consumption equal to 64.3%. The pronounced participation of this reaction is in line with the observations of Su et al.14 and can be attributed to the low value of the reaction enthalpy of hydrogen abstraction from methane, ΔH2, for the studied catalyst. As observed by these authors and as can also be seen from eq A5 in the Supporting Information, due to thermodynamic consistency the low value of ΔH2 leads to a strong endothermicity of the hydroxyl recombination reaction C5, which in turn enhances the surface coverage of OH*. From the main catalytic cycle described above, hydroxyl recombination C5 is the most highly activated step for this catalyst. In the previous section, it was seen that multiple other reactions further contribute to the high hydroxyl surface coverage. From Figure 4, it can be seen that more than 70% of the catalyst surface is covered with OH* at these conditions, explaining the high consumption rate of CH3· radicals by the reverse reaction C2. Hence, the recognized high degree of surface hydroxylation16,21 is confirmed also in the results of the current work. The resulting high rate of CH3· radicals’ conversion back to methane was also identified in the simulations by Sinev.20 In all cases though, the net production rate of CH3· remains clearly positive in the intraparticle phase, the radical being produced by surface reactions and consumed at a lower rate in the pores by gas phase reactions, while in the interstitial phase there is clearly a net consumption of the radical. As expected, in the interstitial phase, CH3· radical coupling accounts for the majority of C2H6 production (98.2%). In the intraparticle phase, the contribution of the CH3· coupling reaction to C2H6 production is still very high and equal to 96.7%. The backward reaction C4, despite the high surface coverage by OH*, contributes only by 3.2% to the produced C2H6, since C2H5· radicals are primarily consumed via reaction C17 toward C2H4, the latter step having a much lower activation barrier. In total, the intraparticle phase accounts for about 53% of the ethane production. Moreover, 92.7% of the produced ethane is consumed via the forward reaction C4, with gas phase reaction with OH· and CH3· toward C2H5· contributing for the rest of this percentage. The production of C2H4 in the interstitial phase is mainly attributed to the dehydrogenation and oxy-dehydrogenation of C2H5· radicals with the former reaction contributing by 51.6% and the latter by 27.5%. Decomposition of C3H7· radicals is the only other reaction leading to the production of C2H4 at a relatively high percentage of 20.8%. The intraparticle production of C2H4, though, is dominated, as mentioned previously, by the surface H abstraction from C2H5· (96.9%). In total, about 75% of C2H4 production takes place in the intraparticle phase. On the contrary, a variety of gas phase reactions account for the consumption of C2H4, with its coupling with CH3· radicals (71.0%) being more pronounced, while in the intraparticle phase catalytic oxidation via reaction C14 is the primary C2H4 consumption route at 92.3%. As such, the overall contribution of this phase is very high and amounts to 83.6%. 1836

dx.doi.org/10.1021/ie403160s | Ind. Eng. Chem. Res. 2014, 53, 1825−1840

Industrial & Engineering Chemistry Research

Article

Table 6. Effect of Pellet Porosity on Conversion, Selectivity, and Most Relevant Intraparticle Rates per Reactor Volume (mol mr−3 s−1) during OCM over Sn/Li/MgOa

Table 7. Effect of Internal Surface Area on Conversion, Selectivity, and Most Relavant Intraparticle Rates per Reactor Volume (mol mr−3 s−1) during OCM over Sn/Li/MgOa

pellet porosity CH4 conversion CO selectivity CO2 selectivity C2H4 selectivity C2H6 selectivity CH4 consumption gas phase catalytic CH4 production gas phase catalytic CH3· consumption gas phase catalytic CH3· production gas phase catalytic a

surface area (m2 gc−1)

0.1

0.35

0.6

7.32 3.06 49.10 17.17 29.23

9.11 3.00 34.58 22.87 37.16

9.88 3.03 30.15 24.73 39.23

0.05 3.57

0.09 3.15

0.11 2.98

0.01 2.46

0.04 1.81

0.06 1.52

0.43 2.97

0.79 2.14

0.96 1.78

a

0.05 3.57

0.10 3.15

0.13 2.98

the increase in the rate of CH4 consumption via catalytic reactions due to the higher surface area, as clearly seen in Table 7. Moreover, it is observed that a respective rise in CH4 consumption via gas-phase reactions also occurs in the intraparticle phase with increasing surface area. This is explained via the enhancement of radicals’ production via catalytic pathways that further promotes the gas phase reactions with CH4. For example, for the change of surface area from 1 to 6 m2 g−1, the contribution of catalytic reaction C21 in the production of H· radicals increases from 21.6% to 53.2%, leading to more than an order of magnitude increase in its net production rate. For similar reasons, OH· radicals exhibit an equivalent rise, resulting in a concentration increase of H· and OH· of four and three times, respectively, the specific radicals being, as discussed, the main reactants with CH4 via homogeneous reactions. The importance of OH-mediated C−H bond activation pathways was also highlighted recently in the work of Takanabe and Iglesia,57 where it was reported that these are crucial in achieving high C2 yields. The selectivity shifts can be explained similarly via the increased participation of the heterogeneous network that directly affects the consumption of CH3· radicals. As seen in Table 7, there is an evident rise in the rate of radicals consumed via catalytic reactions with increasing surface area. The gasphase reactions are also enhanced, as commented, due to the increased concentration of radicals in the catalyst pores, leading to a small increase in the consumption of radicals via homogeneous reactions in the intraparticle phase. However, the relative contribution of the heterogeneous reactions in CH3· consumption rises steadily from 60.7% to 86.3%. Considering the analysis of this paragraph, in combination with the findings of the previous one, the textural features of an optimal catalyst can be elucidated. It has been shown via contribution analysis that the intraparticle phase accounts for the majority of the consumption and production of participating compounds, including the gas-phase intermediates. This statement holds not only for the typically considered CH3· radicals but for all gas phase species via complex surface− radicals interactions. A high catalyst porosity has been found to be very beneficial to OCM performance leading to an increase of CH4 conversion and of C2 selectivity. Due to the large pores of such a material, CH3· coupling is significantly enhanced in

CH4 conversion CO selectivity CO2 selectivity C2H4 selectivity C2H6 selectivity CH4 consumption gas phase catalytic CH3· consumption gas phase catalytic

1

3.5

6

4.59 3.59 16.83 19.71 58.47

9.93 2.68 43.53 20.74 30.77

13.46 1.61 60.67 15.80 19.63

0.05 1.38

0.08 3.94

0.09 5.97

0.50 0.73

0.76 2.95

0.79 4.89

Operating conditions as in Table 5, besides the variation of internal surface area.

Operating conditions as in Table 5, besides the variation of porosity.

a relative rise of CH4 consumption via gas phase reactions over catalytic reactions occurs. Nevertheless, the heterogeneous pathways of CH4 consumption are still considerably more important. At the same time though, a significant decrease in CH4 production via the reverse of catalytic reaction C2 can be observed, due to relative suppression of catalytic pathways, explaining the overall increase in CH4 conversion. As expected, the increase of porosity leads to a considerable increase in the consumption of CH3· radicals via gas phase reactions versus catalytic ones. For an increase of porosity from 0.1 to 0.6 a shift of 12.3−87.7% to 34.0−66.0% is observed in the respective contribution of homogeneous versus heterogeneous intraparticle consumption. The total consumption of CH3· is higher actually in the lower porosity case; however, the promotion of the gas phase pathways in the high porosity case compensates for this, resulting in an important rise of the coupling versus the deep oxidation route. In parallel, CH3· radical production occurs to a larger extent via homogeneous pathways, although catalytic reaction C2 still remains the main production route. The combined effect of these changes in contributions results in a rising net production rate for CH3· radicals in the intraparticle phase. Moreover, the increase of porosity at a given surface area essentially corresponds to a larger average pore diameter, increasing from 3.11 × 10−08 to 1.86 × 10−07 m at the investigated range. As a result, the lifetime of CH3· radicals increases four times, from 1.78 × 10−06 to 7.22 × 10−06 s. As such, their concentration in the interstitial phase is 21% higher by increasing porosity from 0.1 to 0.6, while the second order of the coupling reaction leads to a more pronounced increase of the CH 3· coupling rate compared to gas phase oxidation, namely the reaction with HO2· radicals. 4.4.2. Internal Surface Area. Increasing the internal surface area from 1 to 6 m2 g−1 leads to a sharp rise of CO2 selectivity, accompanied by a corresponding drop in the selectivity of C2 products, mainly affecting C2H6 and, to a lesser extent, also C2H4. At the same time, the CH4 conversion shows a steady rise with increasing surface area. The latter can be attributed to 1837

dx.doi.org/10.1021/ie403160s | Ind. Eng. Chem. Res. 2014, 53, 1825−1840

Industrial & Engineering Chemistry Research

Article

Ai: pre-exponential factor of reaction i, unit reaction dependent As: reactor cross-sectional surface area, mr2 BDEO2: bond dissociation energy of O2, J mol−1 Ci: concentration of component i, mol m−3 dv: average distance between two catalyst pellets, m dp: catalyst pellet diameter, m dt: reactor diameter, m De,i: effective diffusivity of component i, m2 s−1 Dm,i: molecular diffusivity of component i in a gas mixture, m2 s−1 Di,j: binary diffusivity of component i in j, m2 s−1 Di: catalyst descriptor, unit descriptor dependent Ea,i: activation energy of reaction i, J mol−1 E0,i: intrinsic activation barrier of reaction family i, J mol−1 Fv: total volumetric flow rate, mg3 s−1 Ft: total molar flow rate, mol s−1 G°: standard-state Gibb’s free energy, J mol−1 H°: standard-state enthalpy, J mol−1 ki: rate coefficient of reaction i, unit reaction dependent k1: pseudo-first-order rate coefficient, s−1 L: length of catalyst bed, m Lt: density of active sites, mol mc−2 Mi: molar mass of component i, kg mol−1 Ngas: total number of gas-phase reactions Nmol: total number of molecules Nrad: total number of radicals Nsurf: total number of surface species n: reaction order Ptot: total pressure, Pa QH−CH3·: C−H bond energy of methane, J mol−1 QH−X: C−H bond energy of X−H compound, J mol−1 R: universal gas constant, 8.314 J mol−1 K−1 r: interstitial coordinate, m ri: rate of reaction i, mol m−3 s−1 r:̅ average pore radius, m Ri: net production rate of component i, mol m−3 s−1 S°: standard entropy, J mol−1 K−1 Sg: specific surface area, m2 kgc−1 T: temperature, K vj,i: stoichiometric coefficient of component i in reaction j W: catalyst mass, kg xi: molar fraction of component i z: reactor axial coordinate, m

the intraparticle phase, avoiding deep oxidation routes. Of course, this trend cannot extend to extreme porosity values, since in that case the heterogeneous activation of CH4 would decrease to unacceptable limits. A high surface area on the other hand promotes the heterogeneous consumption routes at a level that becomes unselective for OCM. It can be speculated that a high porosity material with a high concentration of surface defects would represent an optimal OCM catalyst, activating sufficiently CH4 on one hand but avoiding CH3· deep oxidation on the other, such that coupling in the pores of the catalyst dominates leading to optimal C2 yields. Potentially, structured catalyst formulations providing similar textural characteristics should be considered for OCM to achieve maximal C2 yields

5. CONCLUSIONS Intraparticle gradients for reactants, products, radicals, and surface species are now accounted for in a one-dimensional heterogeneous reactor model for simulating oxidative coupling of methane (OCM). Significant gradients develop for radicals as well as for the related surface species, even at an intrinsic kinetics regime for reactants and products. These gradients directly impact conversion and selectivity and, hence, are necessary to take into account for adequate simulation of OCM behavior. The use of a combined complex microkinetic and reactor model provides unprecedented insight into the occurring phenomena. The majority of the methane is heterogeneously activated by oxidized sites on the catalyst surface in the intraparticle phase. A smaller, but non-negligible, fraction is activated homogeneously in the intraparticle phase, while only minor amounts are consumed in the interstitial phase. The coupling of the formed CH3· proceeds at significant degree both in the interstitial and in the intraparticle phase. Surface hydroxyl groups are found to affect negatively performance by consuming a considerable part of CH3· in the intraparticle phase. Optimal C2 yields can be achieved by maximizing the porosity of catalyst pellets, provided that a sufficient CH4 activation is achieved, so as to maximize the coupling route over the heterogeneous oxidation of CH3·.



ASSOCIATED CONTENT

S Supporting Information *

Detailed methodology, additional tables as described in the text. This information is available free of charge via the Internet at http://pubs.acs.org/.



Greek Symbols

AUTHOR INFORMATION

αi: transfer coefficient of Polanyi parameter for reaction family i ΔH: standard reaction enthalpy, J mol−1 ΔHf: standard formation enthalpy, J mol−1 εb: bed porosity, mg3 mr−3 εc: catalyst porosity, mg3 mc−3 θi: fraction of the surface covered by species i λi: diffusion length of component i, m ξ: pellet coordinate, m ρc: catalyst density, kg mc3 (Σv)i: diffusion volume of component i, m3 mol−1 τc: catalyst pellet tortuosity, mg2 mc−2 τ: lifetime, s

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This paper reports work undertaken in the context of the project “OCMOL, Oxidative Coupling of Methane followed by Oligomerization to Liquids”. OCMOL is a Large Scale Collaborative Project supported by the European Commission in the seventh Framework Programme (GA No. 228953).



NOMENCLATURE

Roman Symbols

ac: external surface area ag: gas−solid interfacial

Superscripts

of the catalyst, mc mr−3 surface area, mi2 mr−3 2

f: forward reaction step b: backward reaction step 1838

dx.doi.org/10.1021/ie403160s | Ind. Eng. Chem. Res. 2014, 53, 1825−1840

Industrial & Engineering Chemistry Research

Article

Subscripts

(21) Sinev, M. Y.; Fattakhova, Z. T.; Lomonosov, V. I.; Gordienko, Y. A. Kinetics of oxidative coupling of methane: Bridging the gap between comprehension and description. J. Nat. Gas Chem. 2009, 18 (3), 273−287. (22) Couwenberg, P. M.; Chen, Q.; Marin, G. B. Irreducible masstransport limitations during a heterogeneously catalyzed gas-phase chain reaction: Oxidative coupling of methane. Ind. Eng. Chem. Res. 1996, 35 (2), 415−421. (23) Couwenberg, P. M.; Chen, Q.; Marin, G. B. Kinetics of a gasphase chain reaction catalyzed by a solid: The oxidative coupling of methane over Li/MgO-based catalysts. Ind. Eng. Chem. Res. 1996, 35 (11), 3999−4011. (24) Sun, J.; Thybaut, J. W.; Marin, G. B. Microkinetics of methane oxidative coupling. Catal. Today 2008, 137 (1), 90−102. (25) Thybaut, J. W.; Sun, J.; Olivier, L.; Van Veen, A. C.; Mirodatos, C.; Marin, G. B. Catalyst design based on microkinetic models: Oxidative coupling of methane. Catal. Today 2011, 159 (1), 29−36. (26) Mears, D. E. Diagnostic criteria for heat transport limitations in fixed bed reactors. J. Catal. 1971, 20 (2), 127−131. (27) Hoebink, J. H. B. J.; Couwenberg, P. M.; Marin, G. B. Fixed bed reactor design for gas phase chain reactions catalysed by solids: The oxidative coupling of methane. Chem. Eng. Sci. 1994, 49 (24), 5453− 5463. (28) Finlayson, B. A. Packed bed reactor analysis by orthogonal collocation. Chem. Eng. Sci. 1971, 26 (7), 1081−1091. (29) Finlayson, B. A. Orthogonal collocation in chemical reaction engineering. Catal. Rev. Sci. Eng. 1974, 10 (1), 69−138. (30) Serth, R. W. Solution of stiff boundary value problems by orthogonal collocation. Int. J. Numer. Meth. Eng. 1975, 9 (3), 691−699. (31) Ascher, U. M.; Petzold, L. R. Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations; Society for Industrial Mathematics: Philadelphia, 1998. (32) Perry, R. H.; Green, D. W. Perry’s Chemical Engineers’ Handbook, 7th ed.; McGraw-Hill: New York, 1997. (33) Froment, G. F.; Bischoff, K. B. Chemical Reactor Analysis and Design; John Wiley & Sons: New York, 1990. (34) Chen, Q.; Couwenberg, P. M.; Marin, G. B. The oxidative coupling of methane with cofeeding of ethane. Catal. Today 1994, 21 (2−3), 309−319. (35) Mims, C. A.; Mauti, R.; Dean, A. M.; Rose, K. D. Radical chemistry in methane oxidative coupling: Tracing of ethylene secondary reactions with computer models and isotopes. J. Phys. Chem. 1994, 98 (50), 13357−13372. (36) Geerts, J. W. M. H.; Chen, Q.; van Kasteren, J. M. N.; van der Wiele, K. Thermodynamics and kinetic modeling of the homogeneous gas phase reactions of the oxidative coupling of methane. Catal. Today 1990, 6 (4), 519−526. (37) Zanthoff, H.; Baerns, M. Oxidative coupling of methane in the gas phase. Kinetic simulation and experimental verification. Ind. Eng. Chem. Res. 1990, 29 (1), 2−10. (38) Kee, R. J.; Rupley, F. M.; Miller, J. A.; Coltrin, M. E.; Grcar, J. F.; Meeks, E.; Moffat, H. K.; Lutz, A. E.; Dixon-Lewis, G.; Smooke, M. D.; Warnatz, J.; Evans, G. H.; Larson, R. S.; Mitchell, R. E.; Petzold, L. R.; Reynolds, W. C.; Caracotsios, M.; Stewart, W. E.; Glarborg, P.; Wang, H.; Adigun, O.; Houf, W. G.; Chou, C. P.; Miller, S. F.; Ho, P.; Young, D. J. CHEMKIN Release 4.1.1; Reaction Design: San Diego, 2007. (39) Chen, Q.; Couwenberg, P. M.; Marin, G. B. Effect of pressure on the oxidative coupling of methane in the absence of catalyst. AIChE J. 1994, 40 (3), 521−535. (40) Mackie, J. C.; Smith, J. G.; Nelson, P. F.; Tyler, R. J. Inhibition of C2 oxidation by methane under oxidative coupling conditions. Energy Fuels 1990, 4 (3), 277−285. (41) Sinev, M. Modeling of oxidative transformations of light alkanes over heterogeneous catalysts. Russ. J. Phys. Chem. B 2007, 1 (4), 412− 433. (42) Alexiadis, V.; Thybaut, J. W.; Kechagiopoulos, P. N.; Chaar, M.; Muhler, M.; Van Veen, A. C.; Mirodatos, C.; Marin, G. B. Oxidative coupling of methane: Catalytic behavior assessment via comprehensive microkinetic modelling. Appl. Catal., B Submitted.

c: catalyst pellet g: gas phase r: reactor s: refers to Ri calculated via the global reaction network presented in section 3.4 0: initial conditions



REFERENCES

(1) Lunsford, J. H. The catalytic conversion of methane to higher hydrocarbons. Catal. Today 1990, 6 (3), 235−259. (2) Lunsford, J. H. The role of surface-generated gas-phase radicals in catalysis. Langmuir 1989, 5 (1), 12−16. (3) Driscoll, D. J.; Campbell, K. D.; Lunsford, J. H.; D.D. Eley, H. P.; Paul, B. W. Surface-generated gas-phase radicals: Formation, detection, and role in catalysis. Adv. Catal. 1987, 35, 139−186. (4) Chen, Q.; Hoebink, J. H. B. J.; Marin, G. B. Kinetics of the oxidative coupling of methane at atmospheric pressure in the absence of catalyst. Ind. Eng. Chem. Res. 1991, 30 (9), 2088−2097. (5) Fang, T.; Yeh, C. t. Interactions of methane with ThO2 SiO2 surface at 1073 K. J. Catal. 1981, 69 (1), 227−229. (6) Keller, G. E.; Bhasin, M. M. Synthesis of ethylene via oxidative coupling of methane. I. Determination of active catalysts. J. Catal. 1982, 73 (1), 9−19. (7) Boudart, M. From the century of the rate equation to the century of the rate constants: A revolution in catalytic kinetics and assisted catalyst design. Catal. Lett. 2000, 65 (1−3), 1−3. (8) Reyes, S.; Kelkar, C. P.; Iglesia, E. Kinetic-transport models and the design of catalysts and reactors for the oxidative coupling of methane. Catal. Lett. 1993, 19 (2−3), 167−180. (9) Reyes, S. C.; Iglesia, E.; Kelkar, C. P. Kinetic-transport models of bimodal reaction sequences. I. Homogeneous and heterogeneous pathways in oxidative coupling of methane. Chem. Eng. Sci. 1993, 48 (14), 2643−2661. (10) Follmer, G.; Lehmann, L.; Baerns, M. Effect of transport limitations on C2+ selectivity in the oxidative methane coupling reaction using a NaOH/CaO catalyst. Catal. Today 1989, 4 (3−4), 323−332. (11) Aparicio, L. M.; Rossini, S. A.; Sanfilippo, D. G.; Rekoske, J. E.; Treviño, A. A.; Dumesic, J. A. Microkinetic analysis of methane dimerization reaction. Ind. Eng. Chem. Res. 1991, 30 (9), 2114−2123. (12) Wolf, D. Microkinetic analysis of the oxidative conversion of methane. Dependence of rate constants on the electrical properties of (CaO)x(CeO2)1−x catalysts. Catal. Lett. 1994, 27 (1), 207−220. (13) Buyevskaya, O.; Wolf, D.; Baerns, M. Surface processes in the catalytic oxidative coupling of methane to ethane. Recl. Trav. Chim. Pays-Bas 1994, 113 (10), 459−464. (14) Su, Y. S.; Ying, J. Y.; Green, W. H., Jr Upper bound on the yield for oxidative coupling of methane. J. Catal. 2003, 218 (2), 321−333. (15) Simon, Y.; Baronnet, F.; Marquaire, P. M. Kinetic modeling of the oxidative coupling of methane. Ind. Eng. Chem. Res. 2007, 46 (7), 1914−1922. (16) Sinev, M.; Arutyunov, V.; Romanets, A. Kinetic models of C1− C4 alkane oxidation as applied to processing of hydrocarbon gases: Principles, approaches, and developments. Adv. Chem. Eng. 2007, 32, 167−258. (17) Sinev, M. Y. Elementary steps of radical-surface interactions in oxidative coupling of methane. Catal. Today 1992, 13 (4), 561−564. (18) Sinev, M. Y. Kinetic modeling of heterogeneous−homogeneous radical processes of the partial oxidation of low paraffins. Catal. Today 1995, 24 (3), 389−393. (19) Sinev, M. Y.; Bychkov, V. Y.; Korchak, V. N.; Tulenin, Y. P.; Fattakhova, Z. T.; Kalashnikova, O. V. Redox properties and catalytic performance of complex oxides in oxidative coupling of methane. Catal. Today 1994, 21 (2−3), 371−376. (20) Sinev, M. Y. Modeling of oxidative transformations of light alkanes over heterogeneous catalysts. Russ. J. Phys. Chem. B 2007, 1 (4), 412−433. 1839

dx.doi.org/10.1021/ie403160s | Ind. Eng. Chem. Res. 2014, 53, 1825−1840

Industrial & Engineering Chemistry Research

Article

(43) Temkin, M. I. The kinetics of steady-state complex reactions. Int. Chem. Eng. 1971, 11, 709−717. (44) Marin, G. B.; Yablonsky, G. S. Kinetics of Chemical Reactions: Decoding Complexity; Wiley-VCH: Weinheim, 2011. (45) Thybaut, J. W.; Choudhury, I. R.; Denayer, J. F.; Baron, G. V.; Jacobs, P. A.; Martens, J. A.; Marin, G. B. Design of optimum zeolite pore system for central hydrocracking of long-chain n-alkanes based on a single-event microkinetic model. Top. Catal. 2009, 52 (9), 1251− 1260. (46) Katare, S.; Caruthers, J. M.; Delgass, W. N.; Venkatasubramanian, V. An intelligent system for reaction kinetic modeling and catalyst design. Ind. Eng. Chem. Res. 2004, 43 (14), 3484−3512. (47) Mhadeshwar, A. B.; Wang, H.; Vlachos, D. G. Thermodynamic consistency in microkinetic development of surface reaction mechanisms. J. Phys. Chem. B 2003, 107 (46), 12721−12733. (48) Madon, R. J.; Braden, D.; Kandoi, S.; Nagel, P.; Mavrikakis, M.; Dumesic, J. A. Microkinetic analysis and mechanism of the water gas shift reaction over copper catalysts. J. Catal. 2011, 281 (1), 1−11. (49) Maestri, M.; Vlachos, D. G.; Beretta, A.; Groppi, G.; Tronconi, E. A C1 microkinetic model for methane conversion to syngas on Rh/ Al203. AIChE J. 2009, 55 (4), 993−1008. (50) Maestri, M.; Vlachos, D. G.; Beretta, A.; Groppi, G.; Tronconi, E. Steam and dry reforming of methane on Rh: Microkinetic analysis and hierarchy of kinetic models. J. Catal. 2008, 259 (2), 211−222. (51) Donazzi, A.; Maestri, M.; Beretta, A.; Groppi, G.; Tronconi, E.; Forzatti, P. Microkinetic analysis of CH4 CPO tests with CO2-diluted feed streams. Appl. Catal., A 2011, 391 (1−2), 350−359. (52) Dumesic, J. A.; Rudd, D. F.; Aparicio, L. M.; Rekoske, J. E.; Treviño, A. A. The Microkinetics of Heterogeneous Catalysis; American Chemical Society: Washington, DC, 1993. (53) Gordon, S.; McBride, B. J. Computer Program for Calculation of Complex Chemical Equilibrium Compositions, Rocket Performance, Incident and Reflected Shocks and Chapman−Jouguet Detonations; NASA Report SP-273; NASA: Cleveland, OH, 1971. (54) Ito, T.; Wang, J. X.; Lin, C. H.; Lunsford, J. H. Oxidative dimerization of methane over a lithium-promoted magnesium oxide catalyst. J. Am. Chem. Soc. 1985, 107 (18), 5062−5068. (55) Lacombe, S.; Geantet, C.; Mirodatos, C. Oxidative coupling of methane over lanthana catalysts: I. Identification and role of specific active-sites. J. Catal. 1995, 151 (2), 439−452. (56) Buyevskaya, O. V.; Wolf, D.; Baerns, M. Rhodium-catalyzed partial oxidation of methane to CO and H2. Transient studies on its mechanism. Catal. Lett. 1994, 29 (1−2), 249−260. (57) Takanabe, K.; Iglesia, E. Rate and selectivity enhancements mediated by OH radicals in the oxidative coupling of methane catalyzed by Mn/Na2WO4/SiO2. Angew. Chem., Int. Ed. 2008, 47 (40), 7689−7693. (58) Krylov, O. V. Catalytic reactions of partial methane oxidation. Catal. Today 1993, 18 (3), 209−302. (59) Schumacher, N.; Boisen, A.; Dahl, S.; Gokhale, A. A.; Kandoi, S.; Grabow, L. C.; Dumesic, J. A.; Mavrikakis, M.; Chorkendorff, I. Trends in low-temperature water-gas shift reactivity on transition metals. J. Catal. 2005, 229 (2), 265. (60) Michaelides, A.; Liu, Z. P.; Zhang, C. J.; Alavi, A.; King, D. A.; Hu, P. Identification of general linear relationships between activation energies and enthalpy changes for dissociation reactions at surfaces. J. Am. Chem. Soc. 2003, 125 (13), 3704−3705.

1840

dx.doi.org/10.1021/ie403160s | Ind. Eng. Chem. Res. 2014, 53, 1825−1840