Fogler - ECRE - CDROM - Book Page 830 Wednesday, September 17, 2008 5:01 PM

Download as pdf or txt
Download as pdf or txt
You are on page 1of 8

Fogler_ECRE_CDROM.

book Page 830 Wednesday, September 17, 2008 5:01 PM

830 Diffusion and Reaction Chap. 12

1.0
0.8
0.6

0.4

Sphere, zero order


0.2
Sphere, first order
Sphere, second order
0.1
1 2 4 6 8 10 20 40

Zero order φs0 = R k″0 Saρc /DeCA0 = R k 0 /DeCA0

First order φs1 = R k″1 Saρc /De = R k1/De

Second order φs2 =R k″2 SaρcCA0/De = R k 2CA0/De

(a)

Reaction
Rate
1.0
Limited
0.8
0.6 Internal
Internal Diffusion
effectiveness factor Limited
0.4
for different
reaction orders and
catalyst shapes
0.2

0.1
1 0.4 0.6 1 2 4 6 10

R
Sphere φ1 = ( R/3) k″1 Sa ρc/De = ---- k1/De
3
R
Cylinder φ1 = ( R/2) k″1 Sa ρc/De = ---- k1/De
2
Slab φ1 = L k″1 Sa ρc/De = L k 1/De

(b)

Figure 12-5 (a) Effectiveness factor plot for nth-order kinetics on spherical catalyst
particles (from Mass Transfer in Heterogeneous Catalysis, by C. N. Satterfield,
1970; reprint edition: Robert E. Krieger Publishing Co., 1981; reprinted by
permission of the author). (b) First-order reaction in different pellet geometrics
(from R. Aris, Introduction to the Analysis of Chemical Reactors, 1965, p. 131;
reprinted by permission of Prentice-Hall, Englewood Cliffs, N.J.).
Fogler_ECRE_CDROM.book Page 831 Wednesday, September 17, 2008 5:01 PM

Sec. 12.2 Internal Effectiveness Factor 831

Correction for
volume change with
reaction (i.e.,  ≠ 0 )

Figure 12-6 Effectiveness factor ratios for first-order kinetics on spherical


catalyst pellets for various values of the Thiele modulus of a sphere, s , as a
function of volume change. [From V. W. Weekman and R. L. Goring, J. Catal., 4,
260 (1965).]

How can the rate Therefore, to increase the overall rate of reaction, –r′A : (1) decrease the radius
of reaction be R (make pellets smaller); (2) increase the temperature; (3) increase the concen-
increased?
tration; and (4) increase the internal surface area. For reactions of order n, we
have, from Equation (12-20),
2 n–1 2 n–1
2 kn″SaρcR CAs knR CAs
φn = -------------------------------------
- = -------------------------
- (12-20)
De De
For large values of the Thiele modulus, the effectiveness factor is

⎛ 2 ⎞1⁄2 3 ⎛ 2 ⎞1⁄2 3 De (1  n) ⁄ 2
η = ⎜ ----------⎟ ----- = ⎜ ----------⎟ --- ------ CAs (12-35)
⎝ n + 1⎠ φn ⎝ n + 1⎠ R kn

Consequently, for reaction orders greater than 1, the effectiveness factor


decreases with increasing concentration at the external pellet surface.
The preceding discussion of effectiveness factors is valid only for iso-
thermal conditions. When a reaction is exothermic and nonisothermal, the
effectiveness factor can be significantly greater than 1 as shown in Figure 12-7.
Values of greater than 1 occur because the external surface temperature of
the pellet is less than the temperature inside the pellet where the exothermic
reaction is taking place. Therefore, the rate of reaction inside the pellet is
greater than the rate at the surface. Thus, because the effectiveness factor is the
Fogler_ECRE_CDROM.book Page 832 Wednesday, September 17, 2008 5:01 PM

832 Diffusion and Reaction Chap. 12

1000

γ = E = 30
RTs

100

β
=
0.
6
β
=
0.
4
10
β
=
0.
β 2
=
0.
1
η 1
β=
Can you find regions 0
β
where multiple solu- =
-0
tions (MSS) exist? . 8
0.1

0.01

β = CAs(-ΔHRx)De = (ΔT)max
ktTs Ts
0.001
0.01 0.1 1 10 100
φ1
Figure 12-7 Nonisothermal effectiveness factor.

ratio of the actual reaction rate to the rate at surface conditions, the effective-
ness factor can be greater than 1, depending on the magnitude of the parame-
ters  and . The parameter  is sometimes referred to as the Arrhenius
number, and the parameter  represents the maximum temperature difference
that could exist in the pellet relative to the surface temperature Ts .
E
γ = Arrhenius number = --------
RTs
ΔTmax Tmax – Ts –ΔHRx DeCAs
β = ---------------
- = --------------------- = --------------------------------
Ts Ts kt Ts
(See Problem P12-13C for the derivation of .) The Thiele modulus for a
first-order reaction, 1 , is evaluated at the external surface temperature. Typical
values of  for industrial processes range from a value of   6.5 (  0.025,
Typical parameter
values 1  0.22) for the synthesis of vinyl chloride from HCl and acetone to a value
of   29.4 (  6
105, 1  1.2) for the synthesis of ammonia. 5 The
lower the thermal conductivity kt and the higher the heat of reaction, the greater
the temperature difference (see Problems P12-13C and P12-14C). We observe
5 H. V. Hlavacek, N. Kubicek, and M. Marek, J. Catal., 15, 17 (1969).
Fogler_ECRE_CDROM.book Page 833 Wednesday, September 17, 2008 5:01 PM

Sec. 12.3 Falsified Kinetics 833

from Figure 12-7 that multiple steady states can exist for values of the Thiele
modulus less than 1 and when  is greater than approximately 0.2. There will
Criterion for no be no multiple steady states when the criterion developed by Luss 6 is fulfilled.
MSSs in the pellet
4( 1 + β ) > βγ (12-36)

12.3 Falsified Kinetics


You may not be There are circumstances under which the measured reaction order and activa-
measuring what tion energy are not the true values. Consider the case in which we obtain reac-
you think you are.
tion rate data in a differential reactor, where precautions are taken to virtually
eliminate external mass transfer resistance (i.e., CAs = CAb). From these data
we construct a log-log plot of the measured rate of reaction –r′A as a function
of the gas-phase concentration, CAs (Figure 12-8). The slope of this plot is the
apparent reaction order n and the rate law takes the form
n′
–r′A = k′nCAs (12-37)

log

Measured rate with ´


apparent reaction
order n
T2 Slope = n´
T1

CAs log

Figure 12-8 Determining the apparent reaction order (–rA = ρb (–r′A)).

We will now proceed to relate this measured reaction order n to the true reac-
tion order n. Using the definition of the effectiveness factor, note that the
actual rate –r′A , is the product of and the rate of reaction evaluated at the
n
external surface, kn CAs , i.e.,

–r′A = η ( –r′As ) = η ( kn CAs


n
) (12-38)
For large values of the Thiele modulus n , we can use Equation (12-35) to sub-
stitute into Equation (12-38) to obtain

6 D. Luss, Chem. Eng. Sci., 23, 1249 (1968).


Fogler_ECRE_CDROM.book Page 834 Wednesday, September 17, 2008 5:01 PM

834 Diffusion and Reaction Chap. 12

3 2 n 3 D 1–n 2 n
–r′A = ----- ---------- knCAs = --- ------e CAs ---------- knCAs
φn n + 1 R kn n+1

3 2De 1/2 (n + 1)/2


= --- --------------
-k C (12-39)
R ( n + 1 ) n As
We equate the true reaction rate, Equation (12-39), to the measured reaction
rate, Equation (12-37), to get

---------- ⎛ --- ( n  1 ) ⁄ 2⎞
2 3
–r′A = De kn1 ⁄ 2CAs = k′ C n′ (12-40)
n+1 ⎝R ⎠ n As

We now compare Equations (12-39) and (12-40). Because the overall exponent
of the concentration, CAs , must be the same for both the analytical and mea-
sured rates of reaction, the apparent reaction order n is related to the true reac-
tion order n by

The true and the 1+n


apparent reaction n′ = ---------- (12-41)
order 2

In addition to an apparent reaction order, there is also an apparent activa-


tion energy, EApp . This value is the activation energy we would calculate using
the experimental data, from the slope of a plot of ln ( –r′A ) as a function of 1/T
at a fixed concentration of A. Substituting for the measured and true specific
reaction rates in terms of the activation energy gives
EApp ⁄ RT ET ⁄ RT
k′n = AAppe kn = AT e














measured true
into Equation (12-40), we find that

ntrue = 2napparent – 1 1⁄2


⎛ –E ⎞ ⎛ –EApp⎞ n′
–r′A  ⎛ --- ---------- De⎞ AT1 ⁄ 2 exp ⎜ --------T-⎟
3 2 (n1) ⁄ 2
CAs  AApp exp ⎜ -------------
-⎟ CAs
⎝R n+1 ⎠ ⎝ RT ⎠ ⎝ RT ⎠
Taking the natural log of both sides gives us

3 2 ET E App
ln --- ---------- D A 1 ⁄ 2 C (n1) ⁄ 2 – --------- n′
- = ln ( AAppCAs ) – ----------
- (12-42)
R n + 1 e T As 2RT RT
where ET is the true activation energy.
Comparing the temperature-dependent terms on the right- and left-hand
sides of Equation (12-42), we see that the true activation energy is equal to
twice the apparent activation energy.
The true activation
energy ET = 2EApp (12-43)
Fogler_ECRE_CDROM.book Page 835 Wednesday, September 17, 2008 5:01 PM

Sec. 12.4 Overall Effectiveness Factor 835

This measurement of the apparent reaction order and activation energy results
primarily when internal diffusion limitations are present and is referred to as
disguised or falsified kinetics. Serious consequences could occur if the labora-
tory data were taken in the disguised regime and the reactor were operated in
Important industrial a different regime. For example, what if the particle size were reduced so that
consequence of
falsified kinetic
internal diffusion limitations became negligible? The higher activation energy,
ET , would cause the reaction to be much more temperature-sensitive, and there
is the possibility for runaway reaction conditions to occur.

12.4 Overall Effectiveness Factor


For first-order reactions we can use an overall effectiveness factor to help us
analyze diffusion, flow, and reaction in packed beds. We now consider a situa-
tion where external and internal resistance to mass transfer to and within the
pellet are of the same order of magnitude (Figure 12-9). At steady state, the
transport of the reactant(s) from the bulk fluid to the external surface of the cat-
alyst is equal to the net rate of reaction of the reactant within and on the pellet.

Here, both internal


and external
diffusion are
important.

Figure 12-9 Mass transfer and reaction steps.

The molar rate of mass transfer from the bulk fluid to the external surface is
Molar rate = (Molar flux) ⋅ (External surface area)
MA = WAr ⋅ (Surface area/Volume)(Reactor volume)
= WAr ⋅ ac ΔV (12-44)

where ac is the external surface area per unit reactor volume (cf. Chapter 11)
and V is the volume.
This molar rate of mass transfer to the surface, MA , is equal to the net
(total) rate of reaction on and within the pellet:
Fogler_ECRE_CDROM.book Page 836 Wednesday, September 17, 2008 5:01 PM

836 Diffusion and Reaction Chap. 12

MA = –r″A ( External area + Internal area )


External area
External area = -------------------------------------- × Reactor volume
Reactor volume
= ac Δ V

Internal area Mass of catalyst Volume of catalyst


Internal area - × --------------------------------------------- × --------------------------------------------- × Reactor volume
= --------------------------------------
Mass of catalyst Volume of catalyst Reactor volume
ρb






= Sa ρc ( 1 – φ ) Δ V

= Sa ρb Δ V

ρb = Bulk density
= ρc (1 – φ) (12-45)
φ = Porosity MA = –rA″ [ ac ΔV + Sa ρb ΔV ]
See nomenclature note
in Example 12-4.
Combining Equations (12-44) and (12-45) and canceling the volume V,
one obtains
WAr ac  –rA″ (ac  Sa b )
For most catalysts the internal surface area is much greater than the external
surface area (i.e., Sa b  ac ), in which case we have

WAr ac = –rA″ Sa ρb (12-46)

where –rA″ is the overall rate of reaction within and on the pellet per unit sur-
face area. The relationship for the rate of mass transport is
MA  WAr ac V  kc (CAb  CAs )ac V (12-47)
where kc is the external mass transfer coefficient (m/s). Because internal diffu-
sion resistance is also significant, not all of the interior surface of the pellet is
accessible to the concentration at the external surface of the pellet, CAs . We
have already learned that the effectiveness factor is a measure of this surface
accessibility [see Equation (12-38)]:

–rA″ = –rAs
″η

Assuming that the surface reaction is first order with respect to A, we can uti-
lize the internal effectiveness factor to write

–rA″  k″1 CAs (12-48)


Fogler_ECRE_CDROM.book Page 837 Wednesday, September 17, 2008 5:01 PM

Sec. 12.4 Overall Effectiveness Factor 837

We need to eliminate the surface concentration from any equation involv-


ing the rate of reaction or rate of mass transfer, because CAs cannot be mea-
sured by standard techniques. To accomplish this elimination, first substitute
Equation (12-48) into Equation (12-46):

WAr ac  k″1 SaCAs b

Then substitute for WAr ac using Equation (12-47)

kc ac (CAb  CAs )  k″1 Sa bCAs (12-49)

Solving for CAs , we obtain


Concentration at
the pellet surface as kc acCAb
a function of bulk CAs = ----------------------------------- (12-50)
gas concentration
kc ac + ηk″1Sa ρb

Substituting for CAs in Equation (12-48) gives

ηk″1 kc acCAb
–rA″ = ----------------------------------- (12-51)
kc ac + ηk″1Sa ρb

In discussing the surface accessibility, we defined the internal effective-


ness factor with respect to the concentration at the external surface of the
pellet, CAs :

Actual overall rate of reaction


η = -------------------------------------------------------------------------------------------------------------------------------------------------------- (12-28)
Rate of reaction that would result if entire interior surface were
exposed to the external pellet surface conditions, CAs, Ts

Two different We now define an overall effectiveness factor that is based on the bulk
effectiveness factors concentration:

Actual overall rate of reaction


Ω = ------------------------------------------------------------------------------------------------------------------ (12-52)
Rate that would result if the entire surface were
exposed to the bulk conditions, CAb, Tb

Dividing the numerator and denominator of Equation (12-51) by kc ac , we


obtain the net rate of reaction (total molar flow of A to the surface in terms of
the bulk fluid concentration), which is a measurable quantity:
η
–r″A = ------------------------------------------ k″1 CAb (12-53)
1 + ηk″1 Sa ρb ⁄ kc ac

You might also like