A Review of The WGS Reaction Kinetics
A Review of The WGS Reaction Kinetics
A Review of The WGS Reaction Kinetics
R EACTOR E NGINEERING
Volume 8
2010
Review R4
Muruganandam Loganathan
Abstract
The worlds progression towards the Hydrogen economy is facilitating the
production of hydrogen from various resources. In the carbon based hydrogen
production, Water gas shift reaction is the intermediate step used for hydrogen enrichment and CO reduction in the synthesis gas. This paper makes a critical review
of the developments in the modeling approaches of the reaction for use in designing and simulating the water gas shift reactor. Considering the fact that the rate of
the reaction is dependent on various parameters including the composition of the
catalyst, the active surface and structure of the catalyst, the size of the catalyst, age
of the catalyst, its operating temperature and pressure and the composition of the
gases, it is difficult to narrow down the expression for the shift reaction. With different authors conducting experiments still to validate the kinetic expressions for
the shift reaction, continuous research on different composition and new catalysts
are also reported periodically. Moreover the commercial catalyst manufacturers
seldom provide information on the catalyst. This makes the task of designers difficult to model the shift reaction. This review provides a consolidated listing of the
various important kinetic expressions published for both the high temperature and
the low temperature water gas shift reaction along with the details of the catalysts
and the operating conditions at which they have been validated.
KEYWORDS: water gas shift reaction, reaction kinetics
1. INTRODUCTION
The Water Gas Shift Reaction (WGSR) is a reaction traditionally used for the
production of Hydrogen from synthesis gas which is further used for ammonia
production in the fertilizer industry, petroleum refineries for a variety of
operations and recently as fuel for power generation and transportation. The use
of gasification for power generation has also increased the use of water gas shift
reactors multifold. The earliest recording of the reaction dates back to 1888
[Rhodes et al., 1995], and its prominence came with the Haber ammonia synthesis
process and development of catalyst by Bosch and Wilde in 1912 [Twigg, 1989].
The catalyst developed containing iron and chromium was capable of catalyzing
the reaction at 4000C to 5000C and reduced the exit carbon monoxide content to
around 2%.
Water gas shift reaction is a moderately exothermic reversible reaction and
is expressed by
CO + H2O CO2 + H2
(1)
the catalyst. Pasel et al. (2005) has suggested using isothermal reactors based on
his experiments.
Both the high temperature and low temperature catalysts are available
commercially and every vendor has a specific composition and characteristics for
the catalysts. Active research is also carried out across the globe to develop better
catalysts for the reaction. With the commercialization of fuel cell powered
vehicles and fuel cell based power generation, the need for active water gas shift
catalyst is felt as this reaction consumes the bulk space requirements in the power
generation system. Moreover designing the reactor for this reaction requires the
prediction of the rate of the reaction which in turn is predicted by the kinetics of
the reaction. Hence it is imperative that the kinetics of the reaction be known to
designers. Keeping in mind this requirement, a literature search was carried out to
list the various kinetic expressions available in the open literature that could help
in determining the rate of the reaction. Even though the focus was on the kinetic
parameters, the details of the catalysts and the operating conditions on which
these models were developed were also consolidated to help understand the
reaction conditions.
2. HIGH TEMPERATURE SHIFT CATALYSTS
The high temperature catalyst can operate in the temperature range of 3100C to
4500C and are called ferrochrome catalysts because of their composition [Rhodes
et al., 1995]. The reaction is operated adiabatically in the industrial scale where
the temperature increases along the length of the reactor. Hence to prevent the
catalyst bed temperature from going higher that damages the catalyst, the inlet
temperature is normally maintained at 3500C. This inlet temperature can give a
maximum temperature of approximately 5500C at the exit. Newsome (1980) has
reported the typical composition of high temperature shift catalyst as 74.2% Fe2O3,
10.0% Cr2O3, 0.2% MgO and the remaining being volatiles. The Cr2O3 acts as a
stabilizer and prevents the sintering of the iron oxide and its optimal content in the
catalyst is reported by Newsome (1980) to be 14%. To prevent the compromise
on surface area, 8% is used industrially [Rhodes et al., 2002]. The high
temperature reactor is capable of reducing the CO concentration to 3%, which is
the equilibrium concentration at 4500C. The industrial reactors can operate from
atmospheric pressure to 8375 kPa and varying compositions of gases with CO
composition ranging from 3% to 80% [Newsome, 1980]. Rase (1977) has listed
inorganic salts, Boron, Oils, Phosphorous compounds, liquid water (temporary
poison) and sulphur compounds greater than 50 ppm as poisons for the iron
chromia catalyst.
The pretreatment of HT shift catalyst is carried out by partially reducing
the Hematite (Fe2O3) to Magnetite (Fe3O4) using the addition of the process gas
http://www.bepress.com/ijcre/vol8/R4
mixtures to activate the catalyst [Rhodes et al., 1995]. This also converts any
CrO3 present in the catalyst to Cr2O3. The reactions are
3Fe2O3 + H2 2Fe3O4 + H2O
3Fe2O3 + CO 2Fe3O4 + CO2
H = -16.3 kJ/mol
H = +24.8 kJ/mol
(2)
(3)
This partial reduction process should be controlled by controlling the gas feed
mixture concentrations of the constituent gases to achieve the required reduction.
Over reduction of the catalyst can cause damage to the catalyst pellets. A
convenient tool to control the reduction process is the reduction factor proposed
by Rhodes et al. (2002). The typical value of the reduction factor (R) should be
less than 1.2 and values more than 1.6 can cause over reduction.
R = ({[CO] + [H2]} / {[CO2] + [H2O]}]
(4)
H = - 80.8 kJ/mol
(5)
Since the reduction reaction is exothermic, the catalyst should be exposed to the
process stream at temperatures of 2300C and higher temperatures will lead to
sintering of the catalyst [Callaghan et al., 2003]. Moreover the process stream
should contain Hydrogen and steam is not preferable as condensation could affect
the catalyst. Tanaka et al. (2003) found CuMn2O4 and CuAl2O4 mixed oxide
catalyst having conversion of CO more than the commercial catalyst. Henrik
Kusar et al. (2006) have found the copper ceria catalyst to be non-pyrophoric and
stable. Studies on Mn promoted Cu/Al2O3 by Dinesh et al. (2006) showed that
with 8.55 wt% of Mn, 513 K and space time of 5.33 hour, CO conversion of 90%
could be achieved. Many more studies on various modifications in the catalysts
http://www.bepress.com/ijcre/vol8/R4
are being proposed routinely in the literature. But the mainstay of all the catalysts
still is copper.
4. NOBLE METAL CATALYSTS
For compact fuel cell applications, several authors have listed the inherent
disadvantages with the existing commercial catalysts and hence have tried the use
of noble metals as catalysts. In the transportation sector, efforts are on to integrate
reforming of fossil fuel and water gas shift reaction with the fuel cell technologies
to improve modularity, reduce the emission and increase the efficiency of fuel
utilization. The commercial Iron based catalysts are prone to coke formation in
the presence of excess fuel from the reformer [Wheeler et al., 2004]. The use of
the commercial copper based catalysts, due to kinetic limitations will occupy
more volume. Moreover the Cu catalyst is pyrophoric in reduced state and gets
deactivated in the presence of condensed water due to leaching of active
component or formation of surface carbonates [Henrik Kusar et al.(2006),
Wheeler et al.(2004), Mhadeshwar and Vlachos (2005)]. The catalysts should also
be capable of withstanding thermal cycling and faster response during startup and
shutdown. Hence the search for a better noble metal or promoted catalyst is still
being carried out.
Hilaire et al. (2001) found the ceria supported transition metals acting as
catalysts for the water gas shift reaction. Apanee et al. (2003) found the Pt/CeO2
catalyst more active than Au/CeO2 and Au/Fe2O3. Wheeler et al. (2004) studied
the possibility of the WGSR using noble metals and metals with Ceria in the
temperature range of 3000C to 10000C and found the activity of the metals in the
order, Ni>Ru>Rh>Pt>Pd. The measured frequency factor and activation energy
are given in Table 4. Phatak et al. (2007) carried out experiments on Pt supported
on alumina and ceria at varying compositions and reported their order of reactions
and activation energy. Pt exhibited lower turnover rate than Cu.
Gonzaleza et al. (2010) studied the effect of Pt catalysts supported on TiO2,
CeO2 and CeTiO2 and found the Pt supported on Ce-modified TiO2 support
exhibiting better activity than those corresponding to individual ceria and titania
supported catalysts. Most of the recent studies have been directed in using any of
the precious metals like Pt, Rh, Pd and Au deposited on Ceria, Zirconia, Alumina,
Titania, Thoria or Magnesia supports [Andrea et al., 2010]. Latest research works
have started using mixed oxides of the above supports with the precious metals
dispersed in them to effect the low temperature water gas shift reaction.
Ratnasamy and Wagner (2009) have recently published a comprehensive review
on the latest developments on the catalysts for the water gas shift reaction.
5. REACTION KINETICS
The kinetic models provide the easiest way to represent the reaction and help the
designers in determining the rate of the reaction and thus design the reactors.
Basically the kinetic expressions can be classified as microkinetic approach and
the empirical method. The microkinetic method is based on the knowledge about
the elementary steps that are involved in the reaction and its energetics. This
method explores the detailed chemistry of the reaction. Using this method it is
possible to estimate the surface coverage, reaction order and activation enthalpy.
This method provides the accurate pathway and prediction of the reaction, but is
computationally intensive. On the other hand, the empirical models are based on
the experimental results and are typically expressed in the Arrhenius model and
provide an easy and computationally lighter way to predict the rate of reaction.
Most of the design works use the empirical models. The WGSR has been
explained through both the microkinetic approach and empirical approach.
5.1 Micro Kinetic Models
The WGSR mechanism over the metal oxide catalyst has been broadly classified
as
1. Regenerative Mechanism and
2. Associative Mechanism
In the regenerative mechanism also known as the redox mechanism, the oxidation
reduction cycle occurring on the catalyst surface is responsible for the reaction. In
the redox mechanism it is proposed that the catalyst surface is oxidized by H2O,
producing H2 as by product followed by reduction of surface to convert CO to
CO2 as represented in equations 6 and 7.
H2O + red H2 + ox
(6)
CO + ox CO2 +red
(7)
The redox mechanism is used to explain the high temperature water gas shift
reaction. The low temperature shift reaction has been explained by several authors
using both the regenerative mechanism and the associative mechanism.
The mechanism for the water gas shift reaction has been elucidated from
single copper crystal studies. Ovesen et al. (1992) proposed an eight step model
based on surface redox mechanism using single crystal studies on copper. Citing
his studies, he found the reaction structure sensitive with the reaction more on Cu
(1 1 0) than Cu (1 1 1). Wang et al. (2003) investigated the kinetic properties for
http://www.bepress.com/ijcre/vol8/R4
both the forward and reverse reactions and supported the surface redox
mechanism at low temperatures. In their energetic analysis they found the
dissociation of water as the rate controlling step. They too supported the
sensitivity of the reaction to Cu structure. Rhodes et al. (1995) has questioned the
redox mechanism by showing that the reduction reaction as possible and raised
doubt about the possibility of the oxidative step.
The associative mechanism is an adsorption - desorption model where the
adsorbed species interact to form an adsorbed intermediate which then
decomposes to form H2 and CO2 as explained in reaction step 8.
CO + H2 (intermediate) CO2 + H2
(8)
Even though many authors support this mechanism, there is still no convergence
of idea on the nature of the intermediate. The conflict basically in the associative
mechanism is based on the nature of the intermediate formed namely Formate
intermediate mechanism and Carboxyl intermediate mechanism.
Rhodes et al. (1995) explained the associative mechanism of Cu
Chromite as a Langmuir Hinshelwood process with the reaction steps 9 to 13.
CO(g) CO(ads)
(9)
H2O(g) H2O(ads)
(10)
(11)
CO2(ads) CO2(g)
(12)
H2(ads) H2(g)
(13)
They found the intermediate species to be formate. Grenoble and Estadt (1981)
proposed formic acid as the intermediate based on the fact that the principal
decomposition pathways for formic acid involved both the reactants and products
of the water gas shift reaction. Ovesen et al., (1996) further modified the 8 step
microkinetic model he proposed in 1992 based on redox mechanism, by adding 3
more steps to account for the formate mechanism and suggested that the formate
coverage was negligible at atmospheric conditions and was significant at higher
pressures. He found the water dissociation and carbon monoxide oxidation as the
rate limiting steps. Gideon Botes (2007) listed out three different mechanisms and
applied them to the data generated from SASOL (South African Coal, Oil and Gas
Corporation) and found the associative formate mechanism providing better fit for
the data. Waugh (1999) using the microkinetic approach has questioned the
formate mechanism and supported the redox mechanism. Sun et al. (2005) studied
the mechanism of precious metals for the water gas shift reaction and found the
Langmuir Hinshelwood model giving better predictions.
Salmi et al. (1989) carried out transient studies on the water gas shift
reaction and observed that H2 is liberated slowly which contradicted with the
associative mechanism and put forth a mechanism with both adsorptive and
regenerative reaction steps. Hakkarainen et al. (1994) found the reaction
mechanisms different for the ferrochrome catalyst and the Co Mo catalysts
using transient response studies. Keiski et al. (1996) carried out stationary and
transient studies on the high temperature water gas shift reaction and found that
the CO adsorption, CO2 desorption and H2 formation were rate determining steps
in the transient studies. Fishtik and Dutta (2002) found out that the formate and
associative mechanisms were dominant in the low temperature region whereas the
redox mechanism was dominant at higher temperatures.
Gokhale et al. (2008) has proposed the carboxyl mechanism as the
intermediate route for the water gas shift reaction based on studies with Cu(1 1 1).
They found the formate species non reactive and at high pressure contribute to
blocking of active sites which explained the loss in activity with increase in CO2
concentration. Mao et al. (2008) using the Density Functional Theory
investigation has reported that the formate mechanism had low probability
whereas the carboxyl and the redox mechanisms had higher probability on copper
catalyst. Grabow et al. (2008) carried out studies on low temperature water gas
shift reaction on Pt (1 1 1) and found the carboxyl mechanism better suited than
the formate mechanism. Tang et al. (2009) using density functional method and
slab models studied the mechanism of the water gas shift reaction on Cu (1 1 1).
They found the carboxyl and redox mechanisms feasible rather than the formate
and associative mechanisms. Thus an undisputed conclusive mechanism is not
available for the water gas shift reaction so far with many researchers taking sides
for these mechanisms based on their studies. So to carry out the microkinetic
simulation of the water gas shift reaction, the kinetics of the elementary steps can
be used. The 18 step microkinetic model to explain the water gas shift reaction
based on Cu (1 1 1) is given by Callaghan et al. (2003), Fishtik and Dutta (2002)
is listed out in Table 1.
http://www.bepress.com/ijcre/vol8/R4
Table 1 Microkinetic model for forward water gas shift reaction on Cu(1 1 1)
[Callaghan et al. (2003), Fishtik and Dutta (2002)]
Ko
Elementary Step
Adsorption/desorption(1/atm s)
Ea
(kcal/mol)
0
0
25.4
10.7
0
15.5
0
1.4
4
29
26.3
1.3
0.9
14.6
5.3
15.3
5.5
15.3
10
Kinetic Expression
High Temperature Shift [Ref. Twigg (1989)]
[CO 2][ H 2]
k [CO][ H 2O]
r=
(1 + KCO[CO] + KH 2O[ H 2O] + KCO 2[CO 2] + KH 2[ H 2])
Hulburt - Vasan
r=
k [ H 2O ]
1 + K [ H 2O ] / [ H 2 ]
Langmuir
Hinshelwood Model
[CO 2][ H 2]
kKCOKH 2O [CO][ H 2O ]
r=
2
(1 + KCO[CO] + KH 2O[ H 2O] + KCO 2[CO 2] + KH 2[ H 2])
Oxidation
Model
Reduction
Bohlboro et al.
r=
The equation used by Twigg (1989) to calculate the equilibrium constants was
Keq = exp (Z (Z (0.63508-0.29353Z) +4.1778) +0.31688)
(14)
(15)
http://www.bepress.com/ijcre/vol8/R4
11
Keq
200
250
300
350
400
450
500
550
210.82 83.956 38.833 20.303 11.723 7.3369 4.9035 3.4586
The equilibrium constant for the reaction derived from thermodynamics is given
in equation 16.
ln( Keq ) =
5693.5
49170
+ 1.077 ln(T ) + 5.44 104 T 1.125 107 T 2
13.148
T
T2
(16)
Moe (1962) has derived a simple to use empirical model to represent the
equilibrium constant which is sufficiently good for design computations and is
given by
Keq = exp (
4577.8
4.33)
T
(17)
The detailed kinetic expressions for the high temperature, low temperature
and noble metal catalysts were compiled and consolidated in Tables 4 and 5. In
table 4, the kinetic parameters corresponding to the power law type kinetic
expression are listed. These parameters are segregated based on the type of
catalysts and similar catalysts are grouped together. Rhodes et al. (2003) has
carried out a series of experiments on the high temperature catalysts and have
reported their arrhenius parameters which are listed in Table 4. The table also lists
the parameters put forth by Keiski et al. (1996) and San et al. (2009). The power
law type reaction rates for the low temperature catalysts for different conditions
are also listed in Table 4. Choi and Stenger (2003) and Koryabkina et al. (2003)
have provided a comprehensive analysis of the previously published parameters.
A good review of the noble metal catalysts was available in the work of Grenoble
and Estadt (1981), Wheeler et al. (2004) and Phatak et al. (2004). By substituting
the pre exponential factor, activation energy and the powers of the component
gases, the rate of reaction can be computed for various catalysts.
Apart from the power law type of kinetic expression where the orders of
the components are listed, complex kinetic expressions have also been developed
to represent the water gas shift reaction. Podolski and Kim (1974) made a detailed
study of the available models and concluded that their experimental data was
reproduced by both the Langmuir Hinshelwood type model and the power law
type model. The same is supported by Newsome (1980) where he has listed the
12
power law model by Bohlboro. This model was capable of predicting the rate in
the presence of H2S. Rase (1977) has provided a kinetic expression for both the
high temperature and low temperature shift reactions in his book which has been
used by Elnashaie and Elshishini (1993) in their modeling works on the shift
reaction. Singh and Saraf (1977, 1980) extended the kinetics of laboratory
catalysts for both low temperature and high temperature to industrial scale by
using factors for diffusion limitation, age of catalyst, pressure correction and
effect of exposure to H2S. San et al. (2009) carried out the high temperature
reaction for varying catalyst composition and derived two rate equations. Keiski
et al. (1996) too has provided an expression for the high temperature shift reaction.
Thomas and Barton (2009) have recently used the models of San et al. and Keiski
et al. for modeling the heterogeneous reactor by introducing correction factors for
porosity.
Earlier, the kinetic expressions for High temperature were applied to the
low temperature shift also. Newsome (1980) has reported the Temkin model for
the low temperature range. Van Herwijnen and De Jong (1980) have listed the
various kinetic models for the low temperature WGSR. Among them are the
models by Campbell and Moe based on industrial data, Shchibrya, Kulkova and
Temkin based on the redox mechanism and the power law model by Goodridge
and Quazi. Ovesen et al. (1992) introduced a factor to account for the pressure
correction in a power law expression. A pressure correction formulae is also
provided by Rase (1977) to extrapolate the experimental data normally carried out
at atmospheric pressure to industrial data. Amadeo and Laborde (1995) made an
analysis of five model equation, 2 representing redox mechanism and 3
representing Langmuir Hinshelwood model. They found the LH model which
included the adsorption of the four components and surface reaction controlling
the reaction provided best fit to their experimental results.Based on their
experimental results, they provided a modified LH model for the low temperature
shift reaction.
Zhao et al. (1999) using the reduced rate method analyzed the pressure
effects and found the LH model to be better than the other models. Choi and
Stenger (2003) tested the validity of two adsorptive, two regenerative and one
empirical model using Cu/ZnO/Al2O3 catalyst and found that all the models
giving good results with the best fit provided by the single step regenerative
model. They also developed their own empirical equation which is listed in Table
5. Wei et al. (2008) has included the steam to CO ratio in the kinetic expression.
The noble metal catalysts are still evolving and have not reached the
commercial stage. Hence a large number of newer catalyst combinations are being
reported in the literature. The nole metal that shows more promise is platinum and
the kinetic expressions derived for platinum has been reported by Ding and Chan
(2008). Another expression is by Sun et al. (2005) for noble metals.
http://www.bepress.com/ijcre/vol8/R4
13
Table 4 Power Law parameters for various WGSR catalysts and operating conditions
Catalyst
Fe3O4/Cr2O3
8 wt% Cr2O3
Fe3O4/Cr2O3
180 -250m
Fe3O4/Cr2O3
180 -250m
Fe3O4/Cr2O3
180 -250m
Fe3O4/Cr2O3
89%Fe2O3,
9% Cr2O3
80-90%Fe2O3,
8-13% Cr2O3,
1-2%CuO,
80-95%Fe2O3,
5-10% Cr2O3,
1-5%CuO,
CuO/Fe3O4/Cr2O3
180 -250m
CuO/Fe3O4/Cr2O3
180 -250m
CuO/Fe3O4/Cr2O3
180 -250m
Operating
Conditions
Ko
Arrhenius Parameters
Ea
Order of reaction
m
n
q
Reference
0.53
0.47
-0.18
Rhodes et al
( 2003)
Rhodes et al
( 2003)
Rhodes et al
( 2003)
Rhodes et al
( 2003)
Keiski et al.
(1996)
Keiski et al.
(1992)
San et al.
(2009)
1 atm, 4500C
6mm x 6mm
102.845
(mol/ gcat s)
111
(kJ/mol)
-0.36
-0.09
1 atm, 4500C
6mm x 6mm
100.659
(mol/ gcat s)
88
(kJ/mol)
0.9
0.31
-0.156
-0.05
San et al.
(2009)
751
(kJ/mol)
852
(kJ/mol)
851
(kJ/mol)
Rhodes et al
( 2003)
Keiski et al.
(1992)
Keiski et al.
(1992)
14
Rase (1977)
Rase (1977)
Rase (1977)
1 atm, 2000C
1 atm, 2500C
0.45
1.07
0.07
0.55
41.8
(kJ/mol)
47.4
(kJ/mol)
Choi
and
Stenger
(2003)
Salmi and
Hakkarainen
(1989)
Choi
and
Stenger
(2003)
Choi
and
Stenger
(2003)
71
(kJ/mol)
Henrik et al.
(2006)
67
(kJ/mol)
79
(kJ/mol)
0.2
0.6
0.8
0.8
-0.9
-0.9
Koryabkina
et al. (2003)
Koryabkina
et al. (2003)
40%CuO-ZnO
Al2O3(United
Catalysts)
1 atm, 1900C
21.4
(kCal/g mol)
23.8
(kCal/g mol)
27.0
(kCal/g mol)
http://www.bepress.com/ijcre/vol8/R4
Cu-ZnO- Al2O3
Cu/ Al2O3
1 atm,
1800C-2000C
-
10%Cu- Al2O3
1 atm, 1300C
8%CuO-Al2O3
1 atm, 2000C
CuO/MnO2
8%CuO-CeO2
1 atm, 2400C
8%CuO15%CeO2- Al2O3
Cu(111)
1 atm, 2000C
1 atm, 3400C
Cu(110)
1 atm, 3400C
Ru
Ru/ceria
Ni
Ni/Ceria
15
86
(kJ/mol)
69.3
(kJ/mol)
1.4
-0.7
-0.9
1.9
55
(kJ/mol)
62
(kJ/mol)
55
(kJ/mol)
0.3
0.38
0.9
0.8
-0.7
-0.8
0.9
0.4
-0.6
-0.6
0.7
0.6
-0.6
-0.6
0.5-1
(2003)
56
(kJ/mol)
32
(kJ/mol)
71
(kJ/mol)
42
(kJ/mol)
Koryabkina
et al. (2003)
Choi and
Stenger
(2003)
Koryabkina
et al. (2003)
Koryabkina
et al. (2003)
Choi and
Stenger
Koryabkina
et al. (2003)
Koryabkina
et al. (2003)
Koryabkina
et al. (2003)
Koryabkina
et al. (2003)
Wheeler
al. (2004)
Wheeler
al. (2004)
Wheeler
al. (2004)
Wheeler
et
et
et
et
16
Rh
alumina
support
Rh/Ceria
5wt% loading
Rh/ Al2O3
3300C
Rh/ SiO2
3500C
Pd
1.66%Pt/ Al2O3
3000C
10000C
0.008 to 0.05 s
contact time
Coated
on
alumina
support
5wt% loading
1atm, 2250C 2850C
1atm,2850C
3450C
1atm, 2850C
1.66%Pt/ Al2O3
1atm, 3000C
2%Pt/ Al2O3
1atm, 2700C
0.9%Pt/ Al2O3
1atm, 1000C
Pd/Ceria
Pt
Pt/Ceria
1%Pt/ Al2O3
1%Pt/ Al2O3
5.10 x 106
(molecules/s/ site)
3.23 x 105
(molecules/s/ site)
- 4.0 x 106 (s-1)
4.0 x 107 (s-1)
1.0 x 106 (s-1)
2.5 x 107 (s-1)
-
(kJ/mol)
130
(kJ/mol)
130
(kJ/mol)
23.01.3
-0.10
0.44
-0.24
0.53
0.1
1.1
-0.07
-0.44
0.06
-0.09
-0.44
0.11
0.82
-0.06
-0.49
0.1
0.77
-0.08
-0.46
-0.21
0.75
0.02
0.55
-0.22
(kcal/mole)
22.82.5
(kcal/mole)
100
(kJ/mol)
100
(kJ/mol)
80
(kJ/mol)
80
(kJ/mol)
68
(kJ/mol)
84
(kJ/mol)
81
(kJ/mol)
81
(kJ/mol)
82
(kJ/mol)
http://www.bepress.com/ijcre/vol8/R4
al. (2004)
Wheeler
al. (2004)
Wheeler
al. (2004)
Grenoble
al. (1981)
Grenoble
al. (1981)
Wheeler
al. (2004)
Wheeler
al. (2004)
Wheeler
al. (2004)
Wheeler
al. (2004)
Phatak et
(2007)
Phatak et
(2007)
Phatak et
(2007)
Phatak et
(2007)
Phatak et
(2007)
Phatak et
(2007)
et
et
et
et
et
et
et
et
al.
al.
al.
al.
al.
al.
0.4%Pt/ Al2O3
1atm, 5440C
Pt/ Al2O3
2700C
1%Pt/ CeO2
1atm, 2000C
1.9 x 106
(molecules/s/ site)
-
1%Pt/ CeO2
1atm, 2400C
1% Pt/CeO2
Pt/SiO2
3400C
Pt/C
3400C
1% Pt/TiO2
1.18 x 105
(molecules/s/ site)
3.84 x 106
(molecules/s/ site)
-
1.4%Pt8.3%
CeO2/Al2O3
2%Pt1%
Re/CeO2-ZrO2
CuO.2CeO.8O2-y
(Cu ceria)
1atm, 2600C
1atm,
2100C-2600C
473 K 623K 1.8 x 103 (s-1)
CO/H2O =1/3
CuO.1CeO.9O2-y
(Cu ceria)
Ce(La)Ox
5
at.%
Cu- 175 300 0C
Ce(10%La)Ox
CO/H2O=1.5
5
at.%
Ni- 275 300 0C
Ce(10%La)Ox
CO/H2O=1.5
17
39
(kJ/mol)
19.61.3
0.45
0.37
-0.73
-0.21
0.75
-0.03
0.44
-0.09
-0.38
0.14
0.66
-0.54
-0.08
-0.08
0.69
0.13
0.35
0.3
0.85
-0.67
0.13
0.49
-0.12
-0.45
-0.05
-0.32
0.85
-0.05
(kcal/mole)
75
(kJ/mol)
46
(kJ/mol)
91
(kJ/mol)
19.10.8
(kcal/mole)
25.51.4
(kcal/mole)
59
(kJ/mol)
86
(kJ/mol)
71
(kJ/mol)
61
(kJ/mol)
78
(kJ/mol)
58.5
(kJ/mol)
30.4
(kJ/mol)
38.2
(kJ/mol)
Phatak et al.
(2007)
Grenoble et
al. (1981)
Phatak et al.
(2007)
Phatak et al.
(2007)
Thinon et al.
(2009)
Grenoble et
al. (1981)
Grenoble et
al. (1981)
Thinon et al.
(2009)
Phatak et al.
(2007)
Phatak et al.
(2007)
Henrik et al.
(2006)
Henrik et
(2006)
Li et
(2000)
Li et
(2000)
Li et
(2000)
al.
al.
al.
al.
18
For Table 4, the rate is given by r = kPlCOPmH2OPnCO2Pq H2 (1-), rate constant k =Ko exp (-Ea/RT), Ko is pre
exponential or frequency factor, Ea is activation energy and is the approach to equilibrium and is given by
= PCO2PH2/KeqPCOPH2O
Table 5 Macro kinetic Rate expressions for WGSR
Catalyst
information
Operating
Conditions
Rate Expression
Girdler(G3-b)
http://www.bepress.com/ijcre/vol8/R4
72% Fe2O3
8% Cr2O3
Density
g/cm3
19
Sud Chemie
SHT-4
Density =
7633.65 kg/m3
r is rate of reaction (cm3/gcat h), XCO is mole fraction of CO, XCO* is mole
fraction of CO in equilibrium condition P is pressure (atm), Keq is equilibrium
constant, R is gas constant(cal/gmol K), T is temperature (K), Eff is effectiveness
factor, Ra is relative activity, Agf is aging factor, Pf is pressure factor, fs is rate
reduction factor due to H2S concentration, [H2S] is H2S concentration in ppm and
is age of catalyst in days
Reference: Wei et al. (2008)
r = k (PCOPH2O PCO2PH2 /keq)
k = 1.78 x 1022 (1+ 0.0097 1.1364 2) T-8 exp(-70/RT)
r is rate of reaction, Pi is partial pressure of species ( kPa), Keq is equilibrium
constant, R is gas constant, T is temperature (K) and is steam to CO ratio
20
ferrochrome
catalysts
Bohlboro
Power Law
Model
0.25
-0.6
25
0.75 0.8
0.50
-0.10
-0.15
0.40
-0.10
-0.10
2000
For small
Catalysts
(crushed/experi
mental)
For Large
Catalysts
(Commercial)
0.93
0.85
0.24
r = k [CO] [H2O]
Ea=105.9 kJ/mol
-0.31
[CO2]
[H2] (1-):
http://www.bepress.com/ijcre/vol8/R4
80-90wt%
Fe2O3,
8-13 % Cr2O3,
1-2% CuO
375 0C to 475
0
C
Temkin Model
21
PCO
R
T 678
r is rate of reaction (mol/kg s), Pi is partial pressure of species (kPa) and T is temperature (K)
= PCO2PH2 / KeqPCOPH2O
Low Temperature Catalysts
Reference: Newsome (1980)
kPH 2OPCO (1 )
4577.8
26800
Rate Constant, k = 6 1011 exp
1/atm s
1.987 T
21500
Constant, A = 2.5 109 exp
1.987 T
r is rate of reaction (1/s), Pi is partial pressure of species (Pa), T is
temperature(K)
r=
22
kKCOKH 2O PCOPH 2O
cat
Keq
r=
2
(1 + KCOPCO + KH 2OPH 2O + KCO 2 PCO 2 ) 60
4577.8
40.32
29364
Rate Constant k = exp
+
mol/g cat. Min
1.987 T 1.987
6.74
6216
12.77
3064
KCO = exp
+
, KH 2O = exp
1.987 T 1.987
1.987 T 1.987
18.45
12542
KCO 2 = exp
1.987 T 1.987
r is rate of reaction (mol/cm3s), Pi is partial pressure of species (Pa), cat is
density of catalyst(g cat/cm3), T is temperature(K)
1855.5
http://www.bepress.com/ijcre/vol8/R4
Girdler(G3-b)
Catalyst
composition,
(%p/p) CuO
(32.7), ZnO
(47) and Al2O3
(11)
particle size
0.62 cm
equivalent
diameter
(1/4 x 1/4)
1 atm,
453 K - 503 K
23
454.3
0.92 e
PCOPH 2O(1 )
T
r=
2
101.5
158.3
2737.9
1596.1
1
2.2
0.4
e
P
e
P
P
P
+
+
CO
H
CO
H
2 O + 0.0047e
2 + 0.05e
2
T
T
T
T
24
1 atm,
Sud Chemie
120 0C - 2500C,
(EX-2248)
Cu/ZnO/ Al2O3 particle size 200
250 m,
CO:H2 1:2
RT
Keq
Girdler/Sud
Chemie (G66
b). unreduced
wt% - CuO
(32.2),
ZnO(61.8) and
Fe2O3(1.6)
CuO (33%)
ZnO (66%)
5.09g/cm3
http://www.bepress.com/ijcre/vol8/R4
25
Sud Chemie
MDC 7
Cu-Zn based
catalyst
Cu/ZnO/Al2O3
40%Cu,
22%Zn, 5%Al
r = k (PCOPH2O keq-1PCO2PH2)
k = 1.74 x 1017 (1- 0.1540 + 0.008 2) T-8.5exp(-35/RT)
r is rate of reaction, Pi is partial pressure (Pa), Keq is equilibrium constant, R is
gas constant, T is temperature (K) and is steam to CO ratio
Copper surface
area =10m2/g
r = k [CO]0.8 [H2O]0.5[CO2]-0.15[H2]0(1-)
[ ] represents the species concentration
Reference: Ovesen et al. (1992)
r = A exp(-Ea/RT)PlCOPmH2OPnCO2PqH2Ptot(1-)
is fudge factor correcting total pressure dependence, = PCO2PH2 / KeqPCOPH2O
Catalyst
P(bar)
Ea(kJ/mol)
l
m
n
q
Cu/ZnO/Al2O3
5
86.5
1
1.4
-0.7
-0.9
Cu/ZnO/Al2O3
20
78.2
1
1.5
-0.7
-0.7
Cu/Al2O3
20
59.3
1.9
-1.4
-0.9
26
Platinum
Precious
Metals
Particle Size =
3.3 mm
Surface area =
80 m2/g
k PCOPH 2O
Keq
r=
2
(1 + KCOPCO + KH 2OPH 2O + KCO 2 PCO 2 + KH 2 PH 2 )
Ko (mol/kgcat s) = 748.824/bar2
Ea (J/mol) = 53,821
KCO2=0.036 (1/bar), KCO=2.222(1/bar)
KH2=2.197 x 10-5(1/bar), KH2O=2.006 x 10-5(1/bar)
Reference: Sun et al. (2005)
PCOPH 2O (1 )
Ea
r = A exp
H 1
H 2
RT
1 + A1 exp RT PCO 1 + A2 exp RT PH 2
7
3
2
-3
A=9 x 10 mol/m atm s, Ea = 44 kJ/mol, A1= 1 x 10 1/s, H1=-46 kJ/mol,
A2=43 1/s, H2=11 kJ/mol
http://www.bepress.com/ijcre/vol8/R4
27
Almost similar kinetic expressions have been proposed for both the high
temperature and low temperature catalysts. The kinetic models that have been
endorsed by many authors from their experiments with various catalysts have
been the Langmuir Hinshelwood model and the power law model. The recent
literature publications use the kinetic expressions of Keiski et al. (1996), San et al.
(2009) for high temperature and Choi and Stenger (2003) for low temperature
water gas shift reaction.
The difference in opinion on the nature of kinetics for the water gas shift
reaction has been attributed to the presence of impurities, mass transfer limitations,
experiments carried out at atmospheric pressure and the use of integral reactor for
kinetic studies rather than the differential reactor [Levent, 2001]. Since all the
experiments were reported at atmospheric pressures, the kinetic models can be
corrected for pressure using the modification recommended by Rase (1977) which
also takes care of the diffusional effects of the catalysts or the pressure correction
factor of Singh and Saraf (1977).
6. CONCLUSION
The use of water gas shift reaction is growing due to the shift towards the
Hydrogen economy. This necessitated the need for designers to have kinetic
expressions for the reaction to help them in designing the water gas shift reactors.
A detailed view of the various catalysts used in the high temperature shift, low
temperature shift and the noble metal catalysts was made. The microkinetic
approach towards the WGSR was explained. A comprehensive listing of the
various kinetic models for the high temperature shift, low temperature shift and
noble metal shift reaction was provided in Table 4 and 5 along with the relevant
numerical values. This consolidated table could provide a ready reference for
computing the rate of the water gas shift reaction. Even though different authors
have obtained different ranges of rates by these models for their experimental data,
most of the works conform to the power law type of kinetic expression and hence
can be conveniently used for design calculations.
28
REFERENCES
http://www.bepress.com/ijcre/vol8/R4
29
30
http://www.bepress.com/ijcre/vol8/R4
31
Rhodes C and Graham J. Hutchings, Studies of the role of the copper promoter
in the iron oxide/ chromia high temperature water gas shift catalyst, Phys.
Chem. Chem. Phys., 5, 2003, 27192723
Salmi T, S. Bostrom and L E Lindfors, A Dynamic Study of the Water-Gas Shift
Reaction over an Industrial Ferrochrome Catalyst Journal of Catalysis,
112, 1988, 345-356
Salmi T and R Hakkarainen, Kinetic Study of the Low-Temperature Water-Gas
Shift Reaction over a Cu-ZnO Catalyst, Applied Catalysis, 49, 1989, 285306
San Shwe Hla, D. Park, G.J. Duffy, J.H. Edwards, D.G. Roberts, A. Ilyushechkin,
L.D. Morpeth, T. Nguyen, Kinetics of high-temperature water-gas shift
reaction over two iron-based commercial catalysts using simulated coalderived syngases, Chemical Engineering Journal 146, 2009, 148154
Seo Y S, Dong-Joo Seo, Yu-Taek Seo and Wang-Lai Yoon, , Investigation of the
characteristics of a compact steam reformer integrated with a water-gas
shift reactor, Journal of Power Sources, 161, 2006, 12081216
Singh C. P. P. and Deoki N. Saraf, Simulation of High-Temperature Water-Gas
Shift Reactors, Ind. Eng. Chem. Proc. Des. Dev., 16 (3), 1977, 313319
Singh C P P and Deoki N. Saraf, "Simulation of Low-Temperature Water-Gas
Shift Reactor", Ind. Eng. Chem. Process Des. Dev., 19 (3), 1980, 393396
Sun J, Joel DesJardins, John Buglass and Ke Liu, Noble metal water gas shift
catalysis: kinetics study and reactor design, International Journal of
Hydrogen Energy, 30, 2005, 1259-1264
Tang Q L, Zhao-Xu Chen and Xiang He, A theoretical study of the water gas
shift reaction mechanism on Cu (1 1 1) model system, Surface Science,
603, 2009, 21382144
Thinon Olivier, Karim Rachedi, Fabrice Diehl,Priscilla Avenier and Yves
Schuurman, Kinetics and Mechanism of the WaterGas Shift Reaction
Over Platinum Supported Catalysts, Top Catal, 52, 2009, 19401945
Thomas A. Adams II, Paul I. Barton, A dynamic two-dimensional heterogeneous
model for water gas shift reactors, International Journal of Hydrogen
Energy, 34, 2009, 8877- 8891
Twigg M V, Catalyst Handbook, second edition, Wolfe Publishing Ltd., 1989
Twigg M V and M S Spencer, Deactivation of supported copper metal catalysts
for hydrogenation reactions, Applied Catalysis A: General, 212, 2001,
161174
Van Herwijnen T and W A De Jong, Kinetics and Mechanism of the CO shift on
Cu/ZnO 1. Kinetics of the Forward and Reverse CO Shift Reactions,
Journal of Catalysis, 63, 1980, 83-93
32
Wang G, Ling Jiang, Yuhua Zho, Zunsheng Cai, Yinming Pan, Xuezhuang Zhao,
Yongwang Li, Yuhan Sun, Bing Zhong, Xianyong Pang, Wei Huang,
Kechang Xie, Investigation of the kinetic properties for the forward and
reverse WGS reaction by energetic analysis, Journal of Molecular
Structure (Theochem), 634, 2003, 2330
Waugh K C, Prediction of global reaction kinetics by solution of the Arrhenius
parameterised component elementary reactions: microkinetic analysis,
Catalysis Today, 53, 1999, 161176
Wei-Hsin Chen, Mu-Rong Lin, Tsung Leo Jiang and Ming-Hong Chen,
Modeling and simulation of hydrogen generation from high temperature
and low-temperature water gas shift reactions, International Journal of
Hydrogen Energy, 33, 2008, 6644 6656
Wei-Hsin Chen, Tai-Ching Hsieh and Tsung Leo Jiang An experimental study
on carbon monoxide conversion and hydrogen generation from water gas
shift reaction, Energy Conversion and Management, 49, 2008, 2801
2808
Wheeler C, A. Jhalani, E.J. Klein, S. Tummala, and L.D. Schmidt, The water
gas-shift reaction at short contact times Journal of Catalysis, 223, 2004,
191199
Tanaka Yohei, Toshimasa Utaka, Ryuji Kikuchi, Kazunari Sasaki and Koichi
Eguchi, Water gas shift reaction over Cu-based mixed oxides for CO
removal from the reformed fuels, Applied Catalysis A: General, 242,
2003, 287295
Zhao Huanqi, Yongqi Hu and Jianjun Li, Reduced rate method for
discrimination of the kinetic models for the watergas shift reaction,
Journal of Molecular Catalysis A: Chemical, 149, 1999, 141146.
http://www.bepress.com/ijcre/vol8/R4