# Dehydrogenation of Methylcyclohexane: Parametric Sensitivity of the Power Law Kinetics

Authors: Muhammad R. Usman; David L. Cresswell; Arthur A. Garforth
Journal: ISRN Chemical Engineering
 (2013)

---

## Abstract

For heterogeneous catalytic reactions, the empirical power law model is a valuable tool that explains variation in the kinetic behavior with changes in operating conditions, and therefore aids in the development of an appropriate and robust kinetic model. In the present work, experiments are performed on 1.0 wt% Pt/Al2O3 catalyst over a wide range of experimental conditions and parametric sensitivity of the power law model to the kinetics of the dehydrogenation of methylcyclohexane is studied. Power law parameters such as order of the reaction, activation energy, and kinetic rate constants are found dependent upon the operating conditions. With H2 in the feed, both apparent order of the reaction and apparent activation energy generally increase with an increase in pressure. The results suggest a kinetic model, which involves nonlinear dependence of rate on the partial pressure of hydrogen and adsorption kinetics of toluene or some intermediate.


---

## Body

## 1. Introduction

Dehydrogenation of a cycloalkane such as methylcyclohexane (MCH) is an important model reaction in reforming of naphtha [1]. A typical reformer feedstock has 20–60 vol% naphthenes [2], which principally undergo dehydrogenation to aromatics. The dehydrogenation of MCH is an essential reaction in the methylcyclohexane-toluene-hydrogen (MTH) system for the safe and economical storage and utilization of hydrogen [3]. Moreover, it can be a valuable model reaction in the refining of naphthenic-based heavy crude oils [4]. For heterogeneous catalytic reactions, the empirical power law model is an important tool in providing the variation in the kinetic behavior to elaborate the insight of the kinetics of a reaction and therefore helps in guiding towards the development of an appropriate and robust kinetic model. Kinetics of the MCH dehydrogenation over supported Pt catalysts has been studied by a number of researchers [1, 4–21]. However, the variation in kinetic behavior and kinetic parameters (power law index, rate constant, and activation energy) of the dehydrogenation reaction with variation in operating conditions is rarely studied [19, 22]. Both Jossens and Petersen [22] and Alhumaidan et al. [19] carried out experiments in the presence of hydrogen and analyzed the data with initial rate method. The present study, on the other hand, is designed to conduct a more rigorous and comprehensive kinetic investigation including experiments without hydrogen in the feed over 1.0 wt% Pt/γ-Al2O3 catalyst. Based on the power law model, kinetic analysis is carried out in which the effect of the operating conditions on the kinetic parameters of the dehydrogenation is studied.


## 2. Experimental

The experimental setup is shown in Figure1. The dehydrogenation reactor was a fixed bed reactor made of 1.02 cm I.D. stainless tube. 2.01 g of 1.0 wt% Pt/γ-Al2O3 catalyst of the size −710 + 425 μm were loaded in the reactor. The reactor tube was placed within a three zone tubular furnace in order to maintain the required temperature. γ-Al2O3 of Alfa-Aesar (Johnson-Matthey) was used as a catalyst support. The support had a pore volume of 0.58 m3/g, median pore diameter of 69 Å, and BET surface area of 208 m2/g. Chloroplatinic acid, H2PtCl6·6H2O, of Sigma-Aldrich was used as platinum source to impregnate the Pt metal over the γ-Al2O3 support. The catalyst was calcined and then reduced in situ. Calcination was carried out in air, while hydrogen was used to reduce the catalyst. The total operating pressure and reactor wall temperature were varied in the range of 1.013 to 9.0 bar and 614.2 to 653.2 K, respectively. Molal space time (feed flowrate) was studied between 3.11 × 104 and 12.44 × 104 s·g-cat/mol-MCH. Feed compositional effect was studied by introducing hydrogen and nitrogen gas. H2 to MCH and N2 to MCH molar ratios were studied between 0–8.4 and 0–1.05, respectively. Four experimental runs were carried out daily. The first experimental run was initiated with 0.25 mL/min MCH, which was followed, respectively, by 0.125 mL/min and 0.5 mL/min. Each of the previous runs lasted for 90 min. The fourth and the last run (for 45 min only) was the repeat of the first run (0.25 mL/min) and facilitated in inspecting the short-term deactivation that had taken place during the course of the dehydrogenation reaction. The products of the dehydrogenation were analyzed in Varian 3400 gas chromatograph containing 100 m nonpolar capillary column and (BP-5: 5% phenyl and 95% dimethylpolysiloxane) equipped with flame ionization detector (FID). The details of the experimental setup and experimental procedure may be found elsewhere [4].Figure 1
Experimental setup.


## 3. Data Analyses Technique and Basic Equations

### 3.1. Grouping of Experimental Data

The experimental data obtained over the 1.0 wt% Pt/γ-Al2O3 catalyst is grouped according to the feed composition, pressure, and temperature. Table 1 shows the manner in which the experimental data is grouped. The data collection in groups is carried out in order to highlight trends in the parameter values, to observe the effects of feed composition on a given fit for a specified pressure, and to signify the effects of pressure itself.Table 1
Groups formation for the experimental rate data obtained for the dehydrogenation of MCH over a 1.0 wt% Pt/γ-Al2O3catalyst.
Group
N
Feed composition
T
w (K)
p (bar)
W
/
F
A
0
×
10
-
4 (s·g-cat-/mol MCH)
y
A
0
y
C
0
y
I
0
1
11
21
0.106
0.893
0.001
614.2, 634.2, 653.2
1.013
3.11, 6.22, 12.44
21
21
0.485
0.511
0.005
614.2, 634.2, 653.2
1.013
3.11, 6.22, 12.44
31
21
0.990
0
0.010
614.2, 634.2, 653.2
1.013
3.11, 6.22, 12.44
41
21
0.485
0
0.515
614.2, 634.2, 653.2
1.013
3.11, 6.22, 12.44
2
15
21
0.106
0.893
0.001
614.2, 634.2, 653.2
5.0
3.11, 6.22, 12.44
25
21
0.485
0.511
0.005
614.2, 634.2, 653.2
5.0
3.11, 6.22, 12.44
45
21
0.485
0
0.515
614.2, 634.2, 653.2
5.0
3.11, 6.22, 12.44
3
19
21
0.106
0.893
0.001
614.2, 634.2, 653.2
9.0
3.11, 6.22, 12.44
29
21
0.485
0.511
0.005
614.2, 634.2, 653.2
9.0
3.11, 6.22, 12.44
49
21
0.485
0
0.515
614.2, 634.2, 653.2
9.0
3.11, 6.22, 12.44


### 3.2. Basic Equations

The dehydrogenation reaction was found highly selective with toluene as the only major product. The kinetic analysis was, therefore, performed for the principal reaction, (1), only:
(1)MCH⇌Toluene+3H2The power law kinetics of the following form containing the effect of short-term deactivation is employed [4]:
(2)(-r)=k·(pA-pB·pC3K)n·(1-kd·td).As mentioned in Section2, the last run was the repeat of the first run (0.25 mL/min) and upon examination shows a decrease in the final conversion of methylcyclohexane over the intervening period. To account for such reversible deactivation, time online deactivation constant, kd, is introduced in (2).The value of the equilibrium constantK was experimentally determined by Schildhauer et al. [23] and is given by
(3)K=3600·exp(-217650R(1T-1650))
with K in bar3, R in J·mol−1·K−1, and T in K.The rate constantk is assumed to follow the Arrhenius temperature dependency and is rearranged in terms of the reference temperature, Tr, as shown next in (4):
(4)k=kr·exp(B·(1-TrT)).The introduction ofTr facilitates in regression of the data and avoids correlation among the parameters [24]. The reference temperature, the central temperature of all the temperatures, is taken as 617.2 K.The termB in (4) is called the dimensionless activation energy and is given by the expression
(5)B=ER·Tr.The average temperatureT as shown in Figure 2 is calculated using the following formula assuming a parabolic temperature distribution in the radial direction [4, 19]:
(6)T=1N·∑i=1i=N(Tz,i+Tw2),
where Tz,i is the measured temperature at the ith axial position at the centerline in the catalyst bed, Tw is the reactor wall temperature, and N is the number of axial temperature measurements.Figure 2
A typical set of temperature profiles in the dehydrogenating reactor.Tw = reactor wall temperature; Tz = measured centerline temperatures in the reactor; Tave = local average temperatures between the previous two values, and T = overall average catalyst bed temperature.The following differential equation is fitted against the experimental data, and four kinetic variablesn, kr, B, and kd are collectively estimated:
(7)dXdY=(-r)×105=k×105·(pA-pB·pC3K)n·(1-kd·td),
where
(8)Y=10-5×WFA0.


### 3.3. Regression Procedure

The regression of the kinetic data is carried out using a FORTRAN code, and the following objective function, the sum of squares of the errors (SSE) is minimized:(9)SSE=∑i=1i=N(Xi,obs-Xi,mod)2,
where Xi,obs is the ith measured or observed value of conversion, Xi,mod is the corresponding value calculated from (7), and N is the total number of data points.


## 3.1. Grouping of Experimental Data

The experimental data obtained over the 1.0 wt% Pt/γ-Al2O3 catalyst is grouped according to the feed composition, pressure, and temperature. Table 1 shows the manner in which the experimental data is grouped. The data collection in groups is carried out in order to highlight trends in the parameter values, to observe the effects of feed composition on a given fit for a specified pressure, and to signify the effects of pressure itself.Table 1
Groups formation for the experimental rate data obtained for the dehydrogenation of MCH over a 1.0 wt% Pt/γ-Al2O3catalyst.
Group
N
Feed composition
T
w (K)
p (bar)
W
/
F
A
0
×
10
-
4 (s·g-cat-/mol MCH)
y
A
0
y
C
0
y
I
0
1
11
21
0.106
0.893
0.001
614.2, 634.2, 653.2
1.013
3.11, 6.22, 12.44
21
21
0.485
0.511
0.005
614.2, 634.2, 653.2
1.013
3.11, 6.22, 12.44
31
21
0.990
0
0.010
614.2, 634.2, 653.2
1.013
3.11, 6.22, 12.44
41
21
0.485
0
0.515
614.2, 634.2, 653.2
1.013
3.11, 6.22, 12.44
2
15
21
0.106
0.893
0.001
614.2, 634.2, 653.2
5.0
3.11, 6.22, 12.44
25
21
0.485
0.511
0.005
614.2, 634.2, 653.2
5.0
3.11, 6.22, 12.44
45
21
0.485
0
0.515
614.2, 634.2, 653.2
5.0
3.11, 6.22, 12.44
3
19
21
0.106
0.893
0.001
614.2, 634.2, 653.2
9.0
3.11, 6.22, 12.44
29
21
0.485
0.511
0.005
614.2, 634.2, 653.2
9.0
3.11, 6.22, 12.44
49
21
0.485
0
0.515
614.2, 634.2, 653.2
9.0
3.11, 6.22, 12.44


## 3.2. Basic Equations

The dehydrogenation reaction was found highly selective with toluene as the only major product. The kinetic analysis was, therefore, performed for the principal reaction, (1), only:
(1)MCH⇌Toluene+3H2The power law kinetics of the following form containing the effect of short-term deactivation is employed [4]:
(2)(-r)=k·(pA-pB·pC3K)n·(1-kd·td).As mentioned in Section2, the last run was the repeat of the first run (0.25 mL/min) and upon examination shows a decrease in the final conversion of methylcyclohexane over the intervening period. To account for such reversible deactivation, time online deactivation constant, kd, is introduced in (2).The value of the equilibrium constantK was experimentally determined by Schildhauer et al. [23] and is given by
(3)K=3600·exp(-217650R(1T-1650))
with K in bar3, R in J·mol−1·K−1, and T in K.The rate constantk is assumed to follow the Arrhenius temperature dependency and is rearranged in terms of the reference temperature, Tr, as shown next in (4):
(4)k=kr·exp(B·(1-TrT)).The introduction ofTr facilitates in regression of the data and avoids correlation among the parameters [24]. The reference temperature, the central temperature of all the temperatures, is taken as 617.2 K.The termB in (4) is called the dimensionless activation energy and is given by the expression
(5)B=ER·Tr.The average temperatureT as shown in Figure 2 is calculated using the following formula assuming a parabolic temperature distribution in the radial direction [4, 19]:
(6)T=1N·∑i=1i=N(Tz,i+Tw2),
where Tz,i is the measured temperature at the ith axial position at the centerline in the catalyst bed, Tw is the reactor wall temperature, and N is the number of axial temperature measurements.Figure 2
A typical set of temperature profiles in the dehydrogenating reactor.Tw = reactor wall temperature; Tz = measured centerline temperatures in the reactor; Tave = local average temperatures between the previous two values, and T = overall average catalyst bed temperature.The following differential equation is fitted against the experimental data, and four kinetic variablesn, kr, B, and kd are collectively estimated:
(7)dXdY=(-r)×105=k×105·(pA-pB·pC3K)n·(1-kd·td),
where
(8)Y=10-5×WFA0.


## 3.3. Regression Procedure

The regression of the kinetic data is carried out using a FORTRAN code, and the following objective function, the sum of squares of the errors (SSE) is minimized:(9)SSE=∑i=1i=N(Xi,obs-Xi,mod)2,
where Xi,obs is the ith measured or observed value of conversion, Xi,mod is the corresponding value calculated from (7), and N is the total number of data points.


## 4. Results and Discussion

Table2 provides the parameter values and overall statistics obtained during the regression of the data in the individual groups. On the other hand, Figure 3 provides a relationship, for the selected groups, between the measured values of conversion, Xobs, and model values, Xmod, to visualize the goodness of the fit. A graphical version of the effect of pressure and composition on the kinetic parameters is shown in Figures 4, 5, and 6.Table 2
Results of the power law regression for the data in individual groups.
Group
n
E
(
-
r
)
0
×
10
5
k
d
SSE
Adj (R2)
F
—
(kJ·mol−1)
(mol MCH·g-cat−1·s−1)
day−1
11
0.99
54.55
4.14
1.35
0.00519
0.992
834.16
21
0.78
58.83
8.00
1.68
0.00989
0.983
401.26
31
0.71
54.55
8.29
1.75
0.00831
0.986
507.20
41
0.69
53.37
4.09
1.67
0.00648
0.991
810.32
15
1.51
138.7
0.30
0.69
0.00545
0.992
866.48
25
1.21
117.2
0.76
1.21
0.00682
0.989
636.61
45
0.98
83.04
1.27
1.67
0.00247
0.996
1840.9
19
1.30
93.70
0.08
0.81
0.00213
0.995
1482.0
29
1.29
119.5
0.15
0.99
0.00255
0.995
1317.2
49
1.20
85.39
0.38
2.01
0.00378
0.992
922.44Figure 3
Scatter diagrams for the power law model relating observed and model values of conversions at 1.013 bar.Figure 4
Effect of pressure on the order of the reaction for the power law kinetics.Figure 5
Effect of pressure on the activation energy of the reaction for the power law kinetics.Figure 6
Effect of pressure on the initial rate of the reaction for the power law kinetics.


### 4.1. 1.013 Bar Pressure

It is obvious that each individual group is relatively better fitted at 1.013 bar by the power law model giving SSE values always less than 0.00989 and the corresponding values of Adj(R2) not less than 0.983.Inspecting orders of the reaction for the individual groups fitted, it is apparent that the order of the reaction decreases from group 11 to group 21 and then remains virtually the same for group 31 and group 41. This observation is crucial in explaining the fact that the concentration of hydrogen and not the MCH concentration in the feed is responsible for change in the order of the reaction. The highest value of the order of the reaction is the result for group 11 in which the concentration of hydrogen is the maximum, that is, 89.3 mol% in the feed. A comparison of the results for group 21 and group 41 allows a direct assessment of the effects of replacing H2 with N2 in the feed, while maintaining a constant partial pressure of MCH. Calculating initial rates at T=Tr=617.2 K, the initial rate of the reaction decreases from 8.0 × 10−5 to 4.09 × 10−5 mol·g-cat−1·s−1 on replacing H2 by N2. This observation should not be a result of irreversible loss of activity, since the periodic activity test showed no long-term activity loss. These observations suggest that H2, a product of the reaction, appears to act as a promoter and has a positive effect towards the kinetics of the reaction, at least at atmospheric pressure. The observation may be explained on the basis of relatively higher toluene inhibition under the conditions when no hydrogen is in the feed. This same promotion in the presence of H2 and toluene inhibition at atmospheric pressure was also observed by other researchers in the field [21, 22, 25]. A reaction order close to unity for group 11 suggests a low MCH coverage. However, a decrease in the order of the reaction with decreasing H2 in the feed suggests an increase in coverage of MCH. Combining the previous statements, it is concluded that the presence of H2 may be helpful in replacing the strongly adsorbed products species, which otherwise cover the active surface. The activation energy parameter, B, remains more or less the same, and the apparent activation energy lies within 53.4 to 58.8 kJ/mol. A wide range of apparent activation energies is reported in the literature. The values previously given tend to be towards the bottom of the range. The deactivation rate constant, kd, is always significant, confirming the importance of including the short-term deactivation.


### 4.2. 5.0 Bar Pressure

Also, atp=5.0 bar, individual groups are fitted very well. It is apparent for group 15 and group 25 in Table 2 that n>1. This seems unlikely on physical grounds. However, the manifestation of an apparent order n>1 can be reconciled with strong chemisorption of one or both of the reaction products toluene or hydrogen, or of reaction intermediates, such as methylcyclohexenes or methylcyclohexadienes. The former intermediate, methylcyclohexene, was observed in low concentrations in the condensate, corresponding to intermediate levels of conversion of methylcyclohexane. Methylcyclohexadienes, on the other hand, have never been observed. This is not to say, however, that they are not formed on the catalyst surface. It is obvious that the apparent order n is significantly greater at 5.0 bar than at 1.013 bar, when hydrogen is in the feed. Langmuir-Hinshelwood-Hougen-Watson (LHHW) postulates, on the other hand, predict a decreasing apparent order n with increasing pressure, consistent with an increasing surface coverage of methylcyclohexane.The values of the activation energy are considerably higher for group 15 and group 25. This shows that pressure has a significant effect on the activation energy of the reaction when H2 is present in the feed. Comparing group 25 and group 45, in contrast to 1.013 bar results, no promotion of hydrogen is observed. The initial rate of reaction, (-r0)×105, is reduced by up to an order of magnitude on increasing the pressure from 1.013 bar to 5.0 bar, which suggests strong product retarding effects. On average, the value of the deactivation constant, kd, is less than that for 1.013 bar and especially for group 15 (highest hydrogen feed concentration), where the 95% confidence interval involves zero, suggesting the parameter kd becomes insignificant.


### 4.3. 9.0 Bar Pressure

In all the cases, at 9.0 bar, an extremely good individual group fit of the data is found with Adj(R2) > 0.99. Similar to the results at 1.013 bar and 5.0 bar, the power law model parameters appear to be group dependent. The value of n remains virtually the same in group 19, group 29, and group 49; a slight increase, however, is noticeable in the groups containing H2. The same observation is observed at 5.0 bar, however, with a greater variation. The order of the reaction n>1.0 can be described in the same way as explained previously. The activation energy is found to be quite high, greater than 85 kJ/mol in all the individual group listings. Comparing group 29 and group 49, an increased initial rate and lower activation energy are observed for group 49. Similarly, a low value of the deactivation constant (kd) is observed when hydrogen is in the feed, which confirms the beneficial effects of hydrogen in maintaining the catalyst activity. A high value of the deactivation constant (kd) in the absence of hydrogen may be explained on the basis of the formation of unsaturated intermediates which act as coke precursors.The previous discussion reveals that at each pressure, parameters are found to be group dependent. This is summarized graphically in Figures4 to 6. When H2 is in the feed, both apparent order of the reaction and apparent activation energy have increased values at high pressures. Actually, for higher concentrations of hydrogen in the feed, these pass through a maximum and then fall back, while for lower hydrogen feed concentrations, the values increase and then remain almost same. An increased value of activation energy provides the clue of some strongly adsorbed components that require higher activation energies of desorption to desorb into the gas phase. A higher order of the reaction, greater than unity, at increased pressure somewhat confirms this hypothesis. The initial rates of reaction at T=Tr=617.2 K tend to decrease with pressure for all the groups. This kind of behavior is less common, though it is compatible with an LHHW dual-site kinetic model. The effect of hydrogen at higher pressures, (say) 9.0 bar, seems to be vanishing as values of parameters somehow approach each other. Comparing group 2 at different pressures, it is observed that at low pressures, a promotion of H2 is observed as mentioned earlier, however, at increased pressure values, in fact, H2 adversely affects the rate which may be manifested that reverse reaction is important at higher pressures than at 1.013 bar or the excess H2 in the feed may be involved in the associative adsorption of some products species which otherwise are gaseous products. Another possibility that at high pressures, hydrogen may be competing for the active adsorption sites and therefore lowering the rates of the reaction may not be ruled out. This effect of hydrogen was observed to be more pronounced at 5.0 bar than at 9.0 bar.In the overall discussion, it may be summarized that the partial pressures of hydrogen and MCH and adsorption of MCH are found to be important contributors at atmospheric pressure, while partial pressures of H2 and adsorption of hydrogen and the other major product (toluene) are found important at higher pressures with adsorption of hydrogen less pronounced. This suggests a kinetic equation that includes the effects of partial pressure of MCH and hydrogen and adsorption kinetics of MCH, hydrogen, and toluene. As with pressures, the initial rates are decreased nonlinearly, so need is there for a term in the denominator that constitutes product of square or cube of the partial pressure of hydrogen (nonlinear dependence of hydrogen partial pressure) and some parameter representing adsorption kinetics of at least one of the major products other than hydrogen, that is, toluene. The best-fit kinetic rate model (based on Langmuir-Hinshelwood-Hougen-Watson single-site kinetics with loss of first hydrogen as the rate controlling step) of the overall experimental data as carried out in our previous study is shown in (10) [4]. Equation (10) clearly shows the strong nonlinear dependence of hydrogen presence and adsorption kinetics of toluene on the rate of the reaction:
(10)(-r)=k·KA·pA·(1-((pB·pC3)/(K·pA)))1+KA·pA+KB·pB+K′·pB·pC2·(1-kd·td),K′=Kr′·exp(B′·(1-TrT)),B′=Δh′R·Tr.


## 4.1. 1.013 Bar Pressure

It is obvious that each individual group is relatively better fitted at 1.013 bar by the power law model giving SSE values always less than 0.00989 and the corresponding values of Adj(R2) not less than 0.983.Inspecting orders of the reaction for the individual groups fitted, it is apparent that the order of the reaction decreases from group 11 to group 21 and then remains virtually the same for group 31 and group 41. This observation is crucial in explaining the fact that the concentration of hydrogen and not the MCH concentration in the feed is responsible for change in the order of the reaction. The highest value of the order of the reaction is the result for group 11 in which the concentration of hydrogen is the maximum, that is, 89.3 mol% in the feed. A comparison of the results for group 21 and group 41 allows a direct assessment of the effects of replacing H2 with N2 in the feed, while maintaining a constant partial pressure of MCH. Calculating initial rates at T=Tr=617.2 K, the initial rate of the reaction decreases from 8.0 × 10−5 to 4.09 × 10−5 mol·g-cat−1·s−1 on replacing H2 by N2. This observation should not be a result of irreversible loss of activity, since the periodic activity test showed no long-term activity loss. These observations suggest that H2, a product of the reaction, appears to act as a promoter and has a positive effect towards the kinetics of the reaction, at least at atmospheric pressure. The observation may be explained on the basis of relatively higher toluene inhibition under the conditions when no hydrogen is in the feed. This same promotion in the presence of H2 and toluene inhibition at atmospheric pressure was also observed by other researchers in the field [21, 22, 25]. A reaction order close to unity for group 11 suggests a low MCH coverage. However, a decrease in the order of the reaction with decreasing H2 in the feed suggests an increase in coverage of MCH. Combining the previous statements, it is concluded that the presence of H2 may be helpful in replacing the strongly adsorbed products species, which otherwise cover the active surface. The activation energy parameter, B, remains more or less the same, and the apparent activation energy lies within 53.4 to 58.8 kJ/mol. A wide range of apparent activation energies is reported in the literature. The values previously given tend to be towards the bottom of the range. The deactivation rate constant, kd, is always significant, confirming the importance of including the short-term deactivation.


## 4.2. 5.0 Bar Pressure

Also, atp=5.0 bar, individual groups are fitted very well. It is apparent for group 15 and group 25 in Table 2 that n>1. This seems unlikely on physical grounds. However, the manifestation of an apparent order n>1 can be reconciled with strong chemisorption of one or both of the reaction products toluene or hydrogen, or of reaction intermediates, such as methylcyclohexenes or methylcyclohexadienes. The former intermediate, methylcyclohexene, was observed in low concentrations in the condensate, corresponding to intermediate levels of conversion of methylcyclohexane. Methylcyclohexadienes, on the other hand, have never been observed. This is not to say, however, that they are not formed on the catalyst surface. It is obvious that the apparent order n is significantly greater at 5.0 bar than at 1.013 bar, when hydrogen is in the feed. Langmuir-Hinshelwood-Hougen-Watson (LHHW) postulates, on the other hand, predict a decreasing apparent order n with increasing pressure, consistent with an increasing surface coverage of methylcyclohexane.The values of the activation energy are considerably higher for group 15 and group 25. This shows that pressure has a significant effect on the activation energy of the reaction when H2 is present in the feed. Comparing group 25 and group 45, in contrast to 1.013 bar results, no promotion of hydrogen is observed. The initial rate of reaction, (-r0)×105, is reduced by up to an order of magnitude on increasing the pressure from 1.013 bar to 5.0 bar, which suggests strong product retarding effects. On average, the value of the deactivation constant, kd, is less than that for 1.013 bar and especially for group 15 (highest hydrogen feed concentration), where the 95% confidence interval involves zero, suggesting the parameter kd becomes insignificant.


## 4.3. 9.0 Bar Pressure

In all the cases, at 9.0 bar, an extremely good individual group fit of the data is found with Adj(R2) > 0.99. Similar to the results at 1.013 bar and 5.0 bar, the power law model parameters appear to be group dependent. The value of n remains virtually the same in group 19, group 29, and group 49; a slight increase, however, is noticeable in the groups containing H2. The same observation is observed at 5.0 bar, however, with a greater variation. The order of the reaction n>1.0 can be described in the same way as explained previously. The activation energy is found to be quite high, greater than 85 kJ/mol in all the individual group listings. Comparing group 29 and group 49, an increased initial rate and lower activation energy are observed for group 49. Similarly, a low value of the deactivation constant (kd) is observed when hydrogen is in the feed, which confirms the beneficial effects of hydrogen in maintaining the catalyst activity. A high value of the deactivation constant (kd) in the absence of hydrogen may be explained on the basis of the formation of unsaturated intermediates which act as coke precursors.The previous discussion reveals that at each pressure, parameters are found to be group dependent. This is summarized graphically in Figures4 to 6. When H2 is in the feed, both apparent order of the reaction and apparent activation energy have increased values at high pressures. Actually, for higher concentrations of hydrogen in the feed, these pass through a maximum and then fall back, while for lower hydrogen feed concentrations, the values increase and then remain almost same. An increased value of activation energy provides the clue of some strongly adsorbed components that require higher activation energies of desorption to desorb into the gas phase. A higher order of the reaction, greater than unity, at increased pressure somewhat confirms this hypothesis. The initial rates of reaction at T=Tr=617.2 K tend to decrease with pressure for all the groups. This kind of behavior is less common, though it is compatible with an LHHW dual-site kinetic model. The effect of hydrogen at higher pressures, (say) 9.0 bar, seems to be vanishing as values of parameters somehow approach each other. Comparing group 2 at different pressures, it is observed that at low pressures, a promotion of H2 is observed as mentioned earlier, however, at increased pressure values, in fact, H2 adversely affects the rate which may be manifested that reverse reaction is important at higher pressures than at 1.013 bar or the excess H2 in the feed may be involved in the associative adsorption of some products species which otherwise are gaseous products. Another possibility that at high pressures, hydrogen may be competing for the active adsorption sites and therefore lowering the rates of the reaction may not be ruled out. This effect of hydrogen was observed to be more pronounced at 5.0 bar than at 9.0 bar.In the overall discussion, it may be summarized that the partial pressures of hydrogen and MCH and adsorption of MCH are found to be important contributors at atmospheric pressure, while partial pressures of H2 and adsorption of hydrogen and the other major product (toluene) are found important at higher pressures with adsorption of hydrogen less pronounced. This suggests a kinetic equation that includes the effects of partial pressure of MCH and hydrogen and adsorption kinetics of MCH, hydrogen, and toluene. As with pressures, the initial rates are decreased nonlinearly, so need is there for a term in the denominator that constitutes product of square or cube of the partial pressure of hydrogen (nonlinear dependence of hydrogen partial pressure) and some parameter representing adsorption kinetics of at least one of the major products other than hydrogen, that is, toluene. The best-fit kinetic rate model (based on Langmuir-Hinshelwood-Hougen-Watson single-site kinetics with loss of first hydrogen as the rate controlling step) of the overall experimental data as carried out in our previous study is shown in (10) [4]. Equation (10) clearly shows the strong nonlinear dependence of hydrogen presence and adsorption kinetics of toluene on the rate of the reaction:
(10)(-r)=k·KA·pA·(1-((pB·pC3)/(K·pA)))1+KA·pA+KB·pB+K′·pB·pC2·(1-kd·td),K′=Kr′·exp(B′·(1-TrT)),B′=Δh′R·Tr.


## 5. Conclusions

Using power law kinetics, the limited data in individual groups is fitted remarkably well, and at each pressure, parameters are found to be group dependent. With H2 in the feed both apparent order of the reaction and apparent activation energy generally increase with pressure. At atmospheric pressure with H2 in the feed, a promotion in the rate of the dehydrogenation reaction is observed. However, at increased pressures, H2 adversely affects the rate. The results suggest a kinetic model, which involves nonlinear dependence of rate on the partial pressure of hydrogen and adsorption kinetics of toluene or some intermediate.



---

*Source: 152893-2013-08-01.xml*