Chemical Reaction Engineering
Rate law is typically determined from experimental data
Elementary reaction kinetics: the reaction orders and stoichiometric coefficients are identical
Nonelementary reaction kinetics: no direct correspondence between reaction order and stoichiometry
\ce{CH3CHO -> CH4 + CO}; \qquad -r_{\ce{CH3CHO}} = k C_{\ce{CH3CHO}}^{3/2}
\ce{H2 + Br2 -> 2 HBr}; \qquad -r_{\ce{HBr}} = \frac{ k_1 C_{\ce{H2}} C_{\ce{Br2}}^{3/2} } { C_{\ce{H2}} + k_2 C_{\ce{Br2}} }
Why the rate laws are complicated? \rightarrow Reaction mechanism
Use the experimental rate law to postulate a reaction mechanism
In practice, knowledge of the reaction mechanism helps use to design better catalyst, troubleshoot, aid in troubleshooting poor reactor performance
Nonelementary kinetics are the result of multiple elementary reaction steps
One or more active intermediates are involved
A high energy species that reacts further very fast. It is consumed as fast as it is formed
It can be
It is present at extremely low concentration.
Net rate of formation of reactive intermediate \approx 0
If the active intermediate appears in n reactions, then
r_{A^*} = \sum_{i = 1}^{n} r_{i, A^*} = 0
Postulate a reaction mechanism that is a series of elementary reactions
Derive a rate equation for the postulated mechanism
Is the rate equation for the postulated mechanism consistent with the experimental results?
Species having the concentration(s) appearing in the denominator of the rate law probably collide with the active intermediate
\ce{M + A* -> [Collision products]} \qquad -r_{A1}^* = k_1 C_M C_{A \cdot}
If a constant appears in the denominator, one of the reaction steps is probably the spontaneous decomposition of the active intermediate
\ce{A* -> [Decomposition products]} \qquad -r_{A2}^* = k_2 C_{A \cdot}
Species having the concentration(s) appearing in the numerator of the rate law probably produce the active intermediate in one of the reaction steps
\ce{[Reactants] -> A* + [Other products]} \qquad r_{A3}^* = k_3 C_{Reactants}
PSSH theory helps to explain why many reactions, such as the decomposition of dimethyl ether commonly exhibit first-order kinetics
\ce{(CH3)2O -> CH4 + H2 + CO}
\ce{ A -> Products }; \qquad -r_A = k C_A
The reaction is first order but not elementary.
Mechanism
Activation:k \ce{A + M ->[k_1] A^* + M}; \qquad r_{1,\ce{A^*}} = k_1 C_A C_M
Deactivation: \ce{A^* + M ->[k_2] A + M}; \qquad r_{2,\ce{A^*}} = k_2 C_{\ce{A^*}} C_M
Decomposition: \ce{A^* ->[k_3] P}; \qquad r_{3,\ce{A^*}} = k_3 C_{\ce{A^*}}
Rate of formation of product
r_p = k_3 C_{\ce{A^*}}
Need to eliminate C_{\ce{A^*}}
PSSH: r_{A^*} = 0
\begin{aligned} r_{A^*} &= k_1 C_A C_M - k_2 C_{A^*} C_M - k_3 C_{A^*} = 0 \\ \\ k_1 C_A C_M &= C_{A^*} (k_2 C_M + k_3) \\ \\ C_{A^*} &= \frac{k_1 C_A C_M}{k_2 C_M + k_3} \\ \\ r_P &= \frac{k_3 k_1 C_A C_M}{k_2 C_M + k_3} \end{aligned}
C_M is constant
k = \frac{k_3 k_1 C_M}{k_2 C_M + k_3}
and
-r_A = k C_A
Definition
Fogler: A substance that affects the rate of a reaction but emerges from the process unchanged.
Ostwald (1895): A catalyst accelerates a chemical reaction without affecting the position of the equilibrium.
Positive catalyst increases reaction rate.
Negative catalyst decreases reaction rate.
About 75% of all chemicals are produced by catalytic process.
Catalyst can be:
Gas (\ce{NO_x}): Nitric Oxide (NO) for the Ostwald Process: In the industrial synthesis of nitric acid, nitric oxide acts as a catalyst in the oxidation of ammonia to nitrogen dioxide. This is part of a gas-phase reaction system where NO is regenerated and acts catalytically.
Liquid (\ce{H2SO4}): Sulfuric Acid in Alkylation: Used in the oil refining industry, sulfuric acid serves as a liquid phase catalyst in alkylation reactions where isobutane is alkylated with olefins to produce higher molecular weight isoparaffins. The reaction occurs in the liquid phase, with sulfuric acid enhancing the reaction by protonating the olefin to form a more reactive carbocation intermediate.
Solid (\ce{Ag}): Silver on Alumina for Ethylene Oxide Production: Silver deposited on an alumina support is a crucial catalyst for the oxidation of ethylene to ethylene oxide, a precursor in the production of ethylene glycol. This is an example of a solid catalyst used in a gas-phase reaction, highlighting the versatility of catalyst forms and applications.
In this lecture only solid catalysts will be considered.
Chemisorption is a necessary part of catalytic processes
Adsorption data are often presented using adsorption isotherms, which show how much gas is adsorbed by a solid at various pressures at a constant temperature.
We first propose an adsorption mechanism and derive an isotherm from it. This theoretical isotherm is then compared with experimental data.
If the theoretical isotherm aligns well with the experimental data, it suggests that the model accurately reflects the physical processes occurring within the system. Discrepancies between the predicted and actual data indicate that the model may not capture one or more crucial aspects of the physical system.
\ce{CO + S <=> CO*S}
\ce{CO + 2S <=> C*S + O*S}
The adsorption of A (gas phase) on an active site S is represented by:
\ce{A(g) + S <=> A*S }
S is a vacant surface site. \ce{A*S} is A bound to a surface site.
Rate of adsorption = rate of attachment - rate of detachment
r_{AD} = k_A P_A C_v - k_{-A} C_{\ce{A*S}}
The rate is directly proportional to the number of collisions between the reactant molecules and the catalyst surface, which is influenced by the partial pressure of the reactant (P_A).
Additionally, the reaction rate is dependent on the availability of vacant active sites on the catalyst surface.
Adsorption equilibrium constant K_A = k_A/k_{-A}
r_{\text{AD}} = k_A P_A C_v - k_{-A} C_{\ce{A*S}} \qquad r_{\text{AD}} = k_A \left( P_A C_v - \frac{k_{-A}}{k_A} C_{\ce{A*S}} \right)
\Rightarrow r_{\text{AD}} = k_A \left( P_A C_v - \frac{C_{\ce{A*S}}}{K_A} \right)
C_t: Total number of active sites per unit mass of catalyst divided by Avogadro’s number (mol/g cat)
C_v: Number of vacant sites per unit mass of catalyst divided by Avogadro’s number
C_v is not measurable, but the total number of sites C_t can be measured
Site balance
In the absence of catalyst deactivation, assume the total number of active sites remains constant:
C_t = C_v + C_{\ce{A*S}} + C_{\ce{B*S}}
We will use the site balance equation to put C_v in terms of measurable species
Site balance
C_t = C_v + C_{\ce{A*S}}
At equilibrium: The net rate of adsorption is zero. Setting r_{AD}=0 we can obtain equation for concentration of A adsorbed on the surface
C_{\ce{A*S}} = K_A C_v P_A \qquad C_{\ce{A*S}} = K_A P_A (C_t - C_{\ce{A*S}})
Langmuir isotherm: By rearranging, we obtain the Langmuir isotherm
C_{\ce{A*S}} = \frac{K_A P_A C_t} {1 + K_A P_A}
This equation relates the equilibrium concentration of A adsorbed on the surface, C_{\ce{A*S}}, as a function of the partial pressure of A
Linearized form
\frac{P_A}{C_{\ce{A*S}}} = \frac{1}{K_A C_t} + \frac{P_A}{C_t}
Consider the adsorption of CO on iron where CO dissociates into separate atoms
\ce{CO + 2S <=> C*S + O*S}
Dissociative adsorption requires two adjacent vacant sites, unlike molecular adsorption that needs only one.
The likelihood of finding two adjacent vacant sites is proportional to the square of the concentration of vacant sites, in line with the law of mass action.
Consequently, the rate of adsorption is proportional to the product of the carbon monoxide partial pressure and the square of the vacant-site concentration.
For desorption, two adjacent occupied sites are needed, making the rate of desorption proportional to the product of the occupied-site concentrations C_\ce{C*S} and C_{\ce{O*S}}.
The net rate of adsorption
r_{AD} = k_A P_{\ce{CO}} C_v^2 - k_{-A} C_{\ce{C*S}} C_{\ce{O*S}}
r_{AD} = k_A \left( P_{\ce{CO}} C_v^2 - \frac{ C_{\ce{C*S}} C_{\ce{O*S}} } { K_A }\right) \qquad K_A = \frac{k_A}{k_{-A}}
At equilibrium r_{AD} = 0 \qquad \Rightarrow \qquad k_A P_{\ce{CO}} C_v^2 = k_{-A} C_{\ce{C*S}} C_{\ce{O*S}}
For C_{\ce{C*S}} = C_{\ce{O*S}}: \qquad (K_A P_{\ce{CO}})^{1/2} C_v = C_{\ce{O*S}}
Substituting for C_{\ce{C*S}}, and C_{\ce{O*S}} in site balance equation
C_t = C_v + C_{\ce{C*S}} + C_{\ce{O*S}}
C_t = C_v + (K_A P_{\ce{CO}})^{1/2} C_v + (K_A P_{\ce{CO}})^{1/2} C_v = C_v (1 + 2(K_A P_{\ce{CO}})^{1/2})
Solving for C_v C_v = \frac{C_t}{1 + 2(K_A P_{\ce{CO}})^{1/2}}
Langmuir isotherm:
C_{\ce{O*S}} = \frac{(K_A P_{\ce{CO}})^{1/2} C_t}{1 + 2(K_A P_{\ce{CO}})^{1/2}} \qquad; \qquad \frac{(P_{\ce{CO}})^{1/2}} {C_{\ce{O*S}}} = \frac{1}{C_t (K_A)^{1/2}} + \frac{2 (P_{\ce{CO}})^{1/2}}{C_t}
After a reactant has been adsorbed onto the surface, \ce{A*S}, is capable of reacting in a number of ways to form the reaction product
Single site mechanism
The surface reaction may be a single-site mechanism in which only the site on which the reactant is adsorbed is involved in the reaction.
An adsorbed molecule of A may isomerize (or perhaps decompose) directly on the site to which it is attached, for example, such as pentene isomerization
\ce{A*S <=> B*S}
r_S = k_S C_{\ce{A*S}} - k_{-S} C_{\ce{B*S}}
r_S = k_S \left( C_{\ce{A*S}} -\frac{C_{\ce{B*S}}}{K_S}\right)
K_S = k_S/k_{-S}
\ce{A*S + S <=> B*S + S}
r_S = k_S \left( C_{\ce{A*S}} C_v -\frac{C_{\ce{B*S}} C_v}{K_S}\right)
\ce{A*S + B*S <=> C*S + D*S}
r_S = k_S \left( C_{\ce{A*S}} C_{\ce{B*S}} -\frac{C_{\ce{C*S}} C_{\ce{D*S}}}{K_S}\right)
\ce{A*S + B*S' <=> C*S' + D*S}
r_S = k_S \left( C_{\ce{A*S}} C_{\ce{B*S'}} -\frac{C_{\ce{C*S'}} C_{\ce{D*S}}}{K_S}\right)
Reactions involving either single- or dual-site mechanisms, are referred to as following Langmuir-Hinshelwood kinetics.
Reaction between adsorbed molecule and a molecule in the gas phase
Unlike the Langmuir-Hinshelwood mechanism it postulates reactions between adsorbed species and molecules directly from the gas phase.
Examples:
\ce{A*S + B (g) <=> C*S}
r_S = k_S \left( C_{\ce{A*S}} P_B -\frac{C_{\ce{C*S}}}{K_S}\right)
Products are desorbed into the gas phase
\ce{C*S <=> C + S}
Rate of desorption
r_{DC} = k_D \left( C_{\ce{C*S}} -\frac{P_C C_v}{K_{DC}}\right)
The desorption of C is the reverse of the adsorption of C
r_{DC} = - r_{AD,C}
The desorption equilibrium constant K_{DC} is the reciprocal of the adsorption equilibrium constant K_C
K_{DC} = \frac{1}{K_C}
r_{DC} = k_D \left[ C_{\ce{C*S}} - K_C P_C C_v \right ]
When heterogeneous reactions are carried out at steady state, the rates of each of the three reaction steps in series (adsorption, surface reaction, and desorption) are equal to one another
-r'_A = r_{AD} = r_S = r_D
One step is usually rate-limiting or rate-controlling
Postulate catalytic mechanism, and then derive the rate law for that mechanism
Assume pseudo-steady state hypothesis (-r'_A = r_{AD} = r_S = r_D)
Identify rate limiting step
Steps to derive the rate law
The approach in determining catalytic and heterogeneous mechanisms is usually termed the Langmuir-Hinshelwood approach, since it is derived from ideas proposed by Hinshelwood based on Langmuir’s principles for adsorption. The Langmuir-Hinshelwood approach was popularized by Hougen and Watson and occasionally includes their names (LHHW kinetics).
Solid catalyst activity changes with operation time.
Catalyst life time is from few minutes to 10 years.
Regeneration of solid catalyst is important
Effects of catalyst decay on solid catalysed reaction
It is very rare to see catalyst activity increases with the operation time
The nickel-based catalysts used in hydroprocessing is often susceptible to sulfur poisoning. However, certain conditioning processes involving controlled exposure to sulfur compounds can rearrange the surface structure, enhancing the catalyst’s activity or selectivity over time.
Amount of sintering is usually measured in terms of active surface area of the catalyst S_a
S_a(t) = \frac{S_{a0}}{1 + k_d t}
Loss of catalytic activity due to a loss of active surface area
Decay law is second order with respect to the present activity
r_d = k_d a^2 = - \frac{da}{dt}
Integrating with a = 1 at time t=0 for constant k_d
a(t) = \frac{1}{1 + k_d t}
The sintering decay constant follows Arrhenius equation
k_d = k_d (T_0) \exp \left[ \frac{E_d}{R} \left( \frac{1}{T_0} - \frac{1}{T} \right) \right ]
When possible, coking can be reduced by running at elevated pressures (2000–3000 kPa) and hydrogen-rich streams.
Catalysts deactivated by coking can usually be regenerated by burning off the carbon.
Results of carbonaceous (coke) material deposited on the surface of a catalyst
This mechanism is common to reactions involving hydrocarbons.
The amount of coke on the surface after a time t
C_C = A t^n
For cracking of crude oil in fixed bed
C_C = 0.52 t^{0.38}
Empirical relationships between a and C_C
\begin{aligned} a &=& \frac{1}{k_{Ck} C_C^p + 1} \\ a &=& e^{-\alpha_1 C_C} \\ a &=& \frac{1}{1 + \alpha_2 C_C} \end{aligned}
\text{Main reaction:} \left\{ \begin{array}{ll} \ce{A + S <=> A*S} \\ \ce{A*S <=> B*S + C(g)} \\ \ce{B*S <=> B + S} \end{array} \right\} \quad -r'_{A} = a(t) \frac{k_{A} C_{A}}{1 + K_{A} C_{A} + K_{B} C_{B}}
\text{Poisoning reaction:} \qquad \ce{P + S -> P*S} \qquad -r'_{d} = -\frac{da}{dt} = k'_d C_p^m a^q
Rate of depletion of the poison from the reactant gas stream onto the catalyst sites, r_{\ce{P*S}}, is proportional to the number of sites that are unpoisoned (C_{t0} - C_{\ce{P*S}}) and the concentration of poison in the gas phase is C_P
r_{\ce{P*S}} = k_d (C_{t0} - C_{\ce{P*S}}) C_P
This rate is also equal to the rate of removal of total active sites (C_t) from the surface
- \frac{dC_t}{dt} = \frac{dC_{\ce{P*S}}}{dt} = r_{\ce{P*S}} = k_d (C_{t0} - C_{\ce{P*S}}) C_P
Dividing through by C_{t0} \frac{df}{dt} = k_d (1 - f) C_P \qquad f: \text{fraction of sites that have been poisoned}
The fraction of sites available for adsorption (1 - f) is the activity a(t)
-\frac{da}{dt} = a(t) k_d C_P
Functional Form of Activity | Decay Reaction Order | Differential Form | Integral Form | Examples |
---|---|---|---|---|
Linear | 0 | -da/dt = \beta_0 | a = 1 - \beta_0 t | Conversion of para-hydrogen on tungsten when poisoned with O2 |
Exponential | 1 | -da/dt = \beta_1 a | a = e^{-\beta_1 t} | Ethylene hydrogenation on Cu poisoned with CO; Cracking of gas oil; Vinyl chloride monomer formation |
Hyperbolic | 2 | -da/dt = \beta_2 a^2 | 1/a = 1 + \beta_2 t | Cyclohexane dehydrogenation on Pt/Al2O3; Isobutylene hydrogenation on Ni |
Reciprocal power | - | -da/dt = \beta_3 a^{1/5} | a = A_0 t^{-\beta_3} | Cracking of gas oil and gasoline on clay |
Chemical Reaction Engineering