Extra Exercises For Chemical Reactor Analysis and Design Fundamentals
Extra Exercises For Chemical Reactor Analysis and Design Fundamentals
for
Chemical Reactor
Analysis and Design
Fundamentals
James B. Rawlings
John G. Ekerdt
D
Department of Chemical Engineering
Austin, Texas
b Hi
No ll Publishing
orders@nobhillpublishing.com
http://www.nobhillpublishing.com
The extra exercises are intended for use by students and course in-
structors.
1 4
P 2 Cl 2 )* 3
- - PCl 5 (2.73)
(b) What is the rank of A ? How many linearly independent reactions can satisfy
Reaction 2.73?
(c) Write all of the possible stoichiometries for these species using the smallest pos-
1 3
H PO 4 2 SrCO 3 )* 3
- - Sr 3 PO
4
4 2 CO 2 5 2
H O
(b) How many linearly independent reactions are possible with these species?
(c) Write all of the possible stoichiometries for these species using the smallest pos-
1
2 The Stoichiometry of Reactions
(b) How many linearly independent reactions are possible with these species?
AT
Hint: you will not obtain integer coef®cients if you use standard software to ®nd
null space of because the vectors of the null space are made into an or-
thonormal set. Try the Octave program nullint which gives an integer-valued
null space.
(a) What is the difference between a set of reactions that is linearly independent and
2NO 2 )*
O - - 2NO
2
is it likely that this overall reaction would occur also as a molecular event? Why
or why not?
A B )*
2C - - D E
)*
A - - B E
2B )*
2C - - D
1. Write out the species list and stoichiometric matrix. For ease of grading,
A student in your class asks for your help with the following dif®cult she is having.
A )*
B - - C
I assume the gas phase is an ideal-gas mixture and the liquid phase is an
ideal-liquid solution. The phase rule says that I have two free intensive vari-
ables
f n 2 r
f 3 2 2 1
f 2
one. I equate the chemical potentials for each component in the two phases
jliq jgas j A; B; C
for three equations, and I know the condition for reaction equilibrium is
n
Y
K aj j
j 1
for the fourth equation. But if I express the reaction equilibrium equation
for the gas phase, I get a different equation than for the liquid phase
aj PyyCy
Y j
gas-phase
j A B
Y j x
aj x xC liquid-phase
j A B
3
4 Review of Chemical Equilibrium
in which I have ignored the Poynting correction for the liquid phase because
the pressure is low. These cannot give the same solution because the phase
Pyj Pj0 xj j A; B; C
in which Pj0 is the vapor pressure of component j at the system T . Using
K PyyCy
A B
PyPy C
A
PyB
0 PC x C 0
0
PA xA
PB xB
0
0 P C x C
PA PB xA xB
K xxxC
A B
which does not agree with the choice of using the liquid phase for the reac-
tion equilibrium. Can I conclude that the component vapor pressures must
be related by
PC0
PA0 PB
1
If not, how do I get consistent results for the equilibrium composition? Na-
tion!
The following reaction takes place at low pressure in the gas phase
A )*
B - - C
You charge a batch reactor with equal amounts of A and B and no C, run the reactor
for a long time at a given temperature and pressure, and then open it up and ®nd that
But you need to reach 75 percent conversion for this process to be pro®table. Your
boss asks you to troubleshoot this process and make suggestions to save the project,
which is about to be canceled. What do you tell your boss? Here are some options to
consider. Feel free to choose one or more of these, or suggest your own options. The
(a) Just run the reactor for a longer time. Eventually it will reach 75 percent conver-
sion, guaranteed. You just have to be patient. State what other thermochemical
(b) Change the operating temperature. If you choose this option, what new tem-
perature do you select? State what other thermochemical data you require, if
any, to select the new temperature, and calculate the new temperature assuming
database.
(c) Change the operating pressure. If you choose this option, what new pressure do
you select? Again state what other thermochemical data, if any, you require to
select the new pressure, and calculate the new pressure assuming whatever new
Consider the volume element V of total volume VR with bounding surface S and outward
normal n depicted in Figure 4.39. The volume element has one entrance and one exit
vj relative to some ®xed laboratory coordinates. The velocity of the element's outer
bounding surface is vs . Decide whether the following statements are true or false, and
Z
@cj d
d Z c d
(a)
Vt
@t dt Vt
j if the volume element is well-mixed
Z
@cj d
dcj V
(b)
Vt
@t dt R if the volume element is well-mixed
d Z c d
Z @cj d
Z c v n
d
j s
(c)
dt Vt
j Vt
@t St
Rj
Q0 V
Q1
cj0 cj1
n
S
vs
Figure 4.39: Volume element V with one entrance stream and one exit stream. The
7
8 The Material Balance for Chemical Reactors
QA
cAf
Q
QB cA
VR
cBf cB
Figure 4.40: Organic acid and base streams fed into a CSTR.
Z
(d) cj vj n
d Q0 cj0 Q1 cj if the volume element is well-mixed
St
Z
(e) cj vj n
d Q0 cj0 Q1 cj cj dV
dt
R if the volume element is well-mixed
St
d Z c d
Z
cj vj vs
nd
Z
Rj d
(f)
dt Vt
j St
Vt
Consider the liquid-phase organic esteri®cation reaction taking place in a CSTR depicted
in Figure 4.40. Two streams, an acid stream containing no base, and a base stream
containing no acid, are fed into the CSTR. The esteri®cation reaction and its rate are
given by
A ! B - ester H O 2 r kcA cB
in which A is the organic acid and B is the organic base. The acid and base are dissolved
in an organic solvent and the acid and base feed streams have feed concentrations cAf
and cBf , respectively. You may assume that the density of the ¯uid is independent of
concentration over the concentration range of interest here. The reactor's volume is
(b) What is the volumetric ¯owrate Q of the ef¯uent stream in terms of the feed
(c) Write out the transient material balances for components A and B. You should
(d) Now consider the steady-state problem. Write the steady-state balances for the
NA0 NA ; NB ; NC
VR
T0 Tout
P0 P0
Figure 4.41: PFR with nonconstant number of moles and nonconstant temperature.
(e) Can you solve these two equations for cAs cBs
, in the general form for all values
(f) Whether or not you were able to solve the equations in the general form, given the
cAs ; cBs
10 min QA QB cAf 8 mol/L cBf 4 mol/L
(g) What are the steady-state molar conversions of A and B for these speci®c param-
eter values?
Consider the gas-phase PFR depicted in Figure 4.41 in which the following reaction
takes place
A - ! B C r kcA
The feed is pure A with molar ¯ow NA0 at temperature T0 , pressure P0 and volumetric
¯owrate Q0 . The outlet of the reactor is at temperature Tout . Neglect pressure drop in
the tube P P0 . The gas may be assumed to be an ideal gas over this temperatures
range at this pressure. The rate constant may be assumed independent of temperature
(a) List the ideal gas equation of state relating temperature T , pressure P and total
molar concentration c ?
(b) What is the relationship between total molar ¯ow N , volumetric ¯owrate Q , tem-
(c) For this reaction stoichiometry, what is the total molar molar ¯ow N expressed
involving NB and NC ?
(d) Assume ®rst that the volumetric ¯owrate does not depend on the change in the
number of moles with reaction, nor on the temperature, i.e. QV
Q0 . Find the
10 The Material Balance for Chemical Reactors
outlet).
(e) Assume next that the tube is isothermal, but account for the change in volumetric
¯owrate due to change in the number of moles with reaction. What reactor vol-
ume is required to achieve 90% conversion? Call this volume VR2 . Find the ratio
VR2 =VR1 . You can obtain a numerical value for this ratio with the given informa-
(f) Finally, assume the change in temperature is well approximated by the linear
relationship
TV
T0 Tout T0
VV
R
Solve the problem again accounting for both the change in temperature and the
change in the number of moles with reaction. Assume the reactor temperature at
the outlet is one and a half times the temperature at the inlet
Tout : T0
1 5 for T in absolute temperature units
What reactor volume is required to achieve 90% conversion? Call this volume VR 3 .
Find the ratio VR3 =VR2 . You can obtain a numerical value for this ratio with the
(g) For this stoichiometry, conversion and temperature change, which effect has more
impact on the reactor size, changing the number of moles by reaction or changing
Revisit Exercise 4.8, but solve the transient reactor starting from a water-®lled reactor of
1000 L assuming an ideal mixture. Be sure to plot the molar numbers of all components,
the reactor volume and the ef¯uent ¯owrate for each case below.
(a) First assume the reactor is operated at constant volume. Plot the solution to
the model showing the transient approach to steady state. Which set of reactor
balances from Tables 4.1 and 4.2 do you choose and why? You may want to read
Exercise 4.19 to get some ideas on what is going on with these different model
forms.
(b) Next assume the reactor has a ®xed and constant outlet ¯owrate
Q Qf
The Material Balance for Chemical Reactors 11
A B
cA0 nB 0
cB 0 nA0
Figure 4.42: Species A and B in a well-mixed volume element. Continuum and molec-
ular settings.
For what range of does the model admit a steady-state solution? Provide a
over this range. Plot this steady-state reactor volume versus over this range.
You may want to display this result using a semi-log plot for the y -axis.
(c) Choose : 0 987 and simulate the transient reactor behavior. Which set of
reactor balances from Tables 4.1 and 4.2 work for this situation?
Let the concentration of A and B molecules in the volume of interest be denoted cA0 ,
cB0 . Consider the three possible irreversible reactions between these species using the
k1
A A - ! C r1 k1 cA2
k2
A B !
- D r2 k2 cA cB
k3
B B !
- E r3 k3 cB2
Consider also the total rate of reaction
r r1 r2 r3
(a) If the A and B species are chemically similar so the different reactions' rate con-
r kcA2 0
3
12 The Material Balance for Chemical Reactors
But if we erase the distinctions between A and B completely and relabel the B
and B as cA cA0 cB
2 , 0 and the total rate is then
r r1 r2 r3
r k1 cA2 k2 cA cB k3 cB2
r k cA0
2 k cA0
k
2
2 2 0 0
r kcA2 0
4
Why are these two total rates different and which one is correct?
(b) Repeat your analysis of the reaction rates if we reduce the length scale and con-
sider the molecular kinetic setting in which we have integer-valued nA0 nB0
,
A chemist colleague comes to you with the following reaction, kinetic data, and pro-
cessing requirements
O OH
2
O
k 2
2 H O
cyclohexanone phenol
or
A 2
O !k
- B 2H O2 r kcA cO2
in which A is cyclohexanone and B is phenol. It is required for downstream processes
available so the reactor is required to achieve the entire 90% conversion. The reactor
also needs to be operated at 0.8 bar so that if there is a leak the hydrocarbons do not
You have several laboratory PFRs available in stock as shown in the table below, and
Reactor Volume
Small A 100 mL
Small B 200 mL
Medium A 1 L
Medium B 2 L
Large 50 L
Other data: T 200 C, k 1500 L/mol hr, R : 0 08314 L bar/mol K; all components
(a) Write out the mass balances for the four components.
(c) Assuming pure reactants in the feed with NAf NO2 f solve the PFR mass balance
for component A.
(d) What molar feed of A and O 2 is required to get 90% conversion in each of the
(e) You need to produce 105 mol/yr of phenol (100 mol/yr for downstream processes
and 5 mol/yr to cover losses). What single reactor do you recommend to minimize
(f) If you can use multiple reactors, what reactors do you recommend to minimize
both the number of reactors and the reactor down time? (i.e., how close can you
get to 105 mol/yr production rate using the minimal number of reactors). You
A - !k B r kcA
is taking place in the single CSTR system shown in Figure 4.43(a). You need to increase
the steady-state conversion without lowering the overall production rate. A creative
but slightly unstable colleague has proposed the modi®cation using a second reactor
VR2
cA2 cA
Q2 Q2 r Q 1
(a) (b)
Figure 4.43: (a) Original single CSTR system. (b) Modi®ed two-CSTR system with recy-
cle.
(a) Compute the steady-state exit dimensionless concentration cA =cAf for the single
(b) Write a total mass balance at the splitting point and determine Q1 and Q2 as a
(c) Write component A mass balances over both reactors and ®nd the concentrations
(d) Evaluate the steady-state overall conversion of A using the two CSTRs for recycle
ratio, r = 1 2. Which conversion is higher, the single CSTR or the two-CSTR with
(e) What do you think is the best recycle ratio to use to maximize the conversion?
What are the advantages and disadvantages of your colleague's proposal to use
the second reactor? Can you suggest a better use of the second reactor?
A k!1
B - C r1 k1 cA2 cB
A k!2
B - D r2 k2 cA cB
There is a large excess of reactant B initially, and therefore we can assume that cB is ap-
proximately constant throughout the batch time. Note that, despite the stoichiometry,
(a) Write out the component balance for species A in the reactor. How does it simplify
(b) Solve the component balance for cA t . Check that your solution satis®es the
initial condition.
(d) How would you de®ne the instantaneous selectivity of product C, sC (production
is measuring.
(e) How would you de®ne the overall selectivity of product C, SC as a function of the
Are the two selectivities you have de®ned the same for this application? Why or
why not?
(f) How would you de®ne overall conversion of A for this reactor, xA ? Explain in
Z
dx
x a
x b
and, since WolframAlpha is unreachable, instead notice that you can express
x a
x b
b a x a x b
1 1 1 1
and you know how to integrate both terms on the right-hand side. This is called a
partial fraction expansion, and you will use it later in CBE 470.
(a) Solve the problem numerically for cA t ; cC t ; cD t using the following parame-
ter values
(b) De®ne
(c) Let's say that you have done an economic analysis of this process and decided that
you need greater than 80 percent conversion of A and greater than 50 percent yield
of C. Are these two requirements feasible? If so, what batch time do you choose
(d) Make a plot of xA versus SC . What is the largest yield of C possible for 80 percent
2A - ! B
The reactor is fed with a mixture of of reactant A and inert diluent I at feedrates NAf
and NIfat total molar concentration cf . The feed conditions and rate constant are
cf :
0 2 mol = L NAf : 2 0 mol = min NIf =NAf 3 k : = 20 0 L mol min
The pressure and temperature are constant in the PFR, and the gas may be assumed to
(a) Write a balance for the steady-state molar ¯ow of A, NA . Make sure that your
2. Double the feed ¯ows NAf and NIf , but keep feed concentration constant?
3. Double the feed concentration, but keep the feed ¯ows constant?
A - ! B
complex as
r kcKcA
2
1A
The Material Balance for Chemical Reactors 17
(a) What are the units of k and K ? Draw a sketch of rcA as a function of cA . Show
(b) Draw a sketch of 1 =rcA as a function of cA . Find the minimum of 1 =rcA . Draw
(c) Assume we run this reaction in liquid phase in a CSTR and also in a PFR.
values, which reactor has more volume, the CSTR or the PFR?
Justify your answer. You might ®nd your sketch of 1 =rcA in the previous part
useful.
(d) To achieve the same 50 percent conversion of A and with the following parameter
values, which reactor has more volume, the CSTR or the PFR?
Justify your answer. You might ®nd your sketch of 1 =rcA in the previous part
useful.
Let's revisit Exercise 4.33 and see if we can make the PFR and CSTR reactors the same
size.
What conversion can you choose so that the sizes of the PFR and CSTR are iden-
tical?
what inlet feed concentration can you choose so that the sizes of the PFR and
(c) Why do you not need to know the reactor rate constant k in Exercise 4.33 and this
exercise?
5
Chemical Kinetics
Consider the reaction mechanism depicted in Figure 5.26 consisting of the following
reactions [2]
k1
A X )*
- - A X fast
k 1
k2
A
X H 2 )*
- - F X
k 2
k3
A
X H 2 )*
- - B X
k 3
k4
G H 2 )*
- - F fast
k 4
k5
H 2 F - ! C 2
H O
k6
B )*
- - D H 2
k 6
k7
D - ! E CO
k8
H 2 B - ! C 2
H O
The rates of reactions 1 and 4 are known to be fast, so these two reactions may be
assumed to be at equilibrium with respect to the other six reactions. A gas-phase CSTR
Typical steady-state concentrations of the gas-phase species in the ef¯uent are listed
in the table.
(a) From the table, what species are good candidates for making the quasi-steady-
choice.
(b) Apply the QSSA to these species and the equilibrium assumption to reactions 1
and 4 and express the production rate of butane (E) in terms of reactants and
products present in large concentration in the gas phase (greater than 5 mole
percent).
19
20 Chemical Kinetics
The H2 are not shown in the diagram for clarity but are listed in the mech-
anism below.
species A
A X G F B D E C 2
H O H 2 CO
mole percent 20 0 23 2 0 0 10 14 14 7 10
Table 5.9: Steady-state gas-phase mole percentages at 545 K, 1.3 bar in the ef¯uent of
a ¯ow reactor. Adsorbed methyl THF (AX) is adsorbed on the catalyst and
amount.
Chemical Kinetics 21
(c) Find the production rate of n-pentane (C) in terms of reactants and products
present in large concentration in the gas phase (greater than 5 mole percent).
Consider the following simple reaction mechanism taking place in a well-mixed, constant-
k1 k2
A -! B - ! C
the QSSA.
cAs cA0 e k1 t
The concentration of C is always available if desired from the total species balance
(b) The B species has two-time-scale behavior. On the fast time scale, it changes
dcB k c c k 1
1 A B
d 2
cBi Y0 Y1 2 Y2
The initial condition, CBi CB0 , must be valid for all , which gives for the initial
conditions of the Yn
Y0
cB0
0 Yn
; n ; ; : : :
0 0 1 2
22 Chemical Kinetics
Substitute the series expansion into B's material balance, collect like powers of
and show the following differential equations govern the Yn
0 dY0 Y
:
d 0
1 dY1 k c Y
1 A
:
d 1
n dYn Y n
n
:
d 2
Y0 c B 0 e
Y1 cA0 k =kk1
e e k1 =k2
1 2 1
Yn 0 n 2
Because Yn vanishes for n 2, show you obtain the exact solution for the B
concentration for all by using the ®rst two terms. Compare your result to Equa-
tion 5.59.
(d) Next we analyze B's large-time-scale behavior, also called the outer solution. Di-
vide B's material balance by k2 again but do not rescale time and obtain
dcB
dt k1 cA cB
Expand cB again in a power series of
cBo B0 B1 2 B2
Substitute the power series into the material balance and collect like powers of
to obtain the following equations
0 : B0 0
1 dB0 k c B
1 A
:
dt 1
dBn B
n n 1 n
:
dt 1
B0 0
B1 k 1 c A
Bn kn1 cA n 2
(e) Show that the classic QSSA analysis is the ®rst-order outer solution.
(f) To obtain a uniform solution valid for both short and long times, we add the inner
and outer solution and subtract any common terms. Plot the uniform zeroth-order
cA0 1 cB 0 = 1 2 k1 1 k2 10
Compare to the exact solution and the ®rst-order outer solution (QSSA solution).
(g) Show that the in®nite order uniform solution is also the exact solution.
Consider again Exercise 5.21 with a slightly more complex reaction mechanism
k1 k2
)*
A - - B - ! C
k 1
species for which we wish to examine the QSSA. Notice that either k 1 or k2 may be
large with respect to the other without invalidating the QSSA assumption for B. Only if
k1 K2
cAs cA0 cB0 K e 1K2 t 1
2 1
k
cBs k 1K cA0 cB0 K e
1 1
k1 K2
1K2 t
1 2 1 2 1
in which K2 k2 =k 1 . Notice that these are in agreement with but more general
than Equations 5.70 and 5.71, which have assumed in addition that cB0 0.
(b) With this mechanism, both the A and B species have two-time-scale behavior, so
we use a series expansion for both cA and cB . Let the inner solution be given by
cAi X0 X1 2 X2
cBi Y0 Y1 2 Y2
in which the small parameter is the inverse of the largest rate constant in the
of and show
0 dX0 Y dY0 K
Y
:
d 0
d 2 0 1
1 dX1 k X Y dY1 k X K
Y
:
d 1 0 1
d 1 0 2 1 1
dXn k X Y dYn k X
n 1 n 1 n 1 n 1 K2
Yn n
:
d d 1 1
(c) Solve these for the zero-order inner solution and show
X0 cA0 cB0 K 1
1 e 1K2
Y0 cB0 e 1K2
2
1
(d) Next we construct the outer solution valid for large times. Postulate a series
cAo A0 A1 2 A2
cBo B0 B1 2 B2
0 : B0 0 K2
B0
1 0
1 dA0 k A B dB0 k A K
B
:
dt 1 0 1
dt 1 0 2 1 1
dAn 1 k A B dBn 1 k A
n 1 n n 1 n 1 K2
Bn n
:
dt dt 1 1
k1 K2
t
A0 A0
e 0
1K2 B0 0
Again we see that to zero order, the B concentration is zero after a short time.
Note also that, unlike in Exercise 5.21, we require an initial condition for the
!1
lim X0
t!0 A0 t
lim
Chemical Kinetics 25
In other words, the long-time solution (steady state) on the fast-time scale is the
short-time solution (initial condition) on the slow-time scale. Using this matching
condition show
A0
cA0 cB0 K
0
1
2 1
(f) Find also the ®rst-order solution, B1 , and show that the QSSA solution corresponds
to the zero-order outer solution for cA and the ®rst-order outer solution for cB .
Consider again the reaction equilibrium assumption and the following reactions
k1 k2
)* ;
A - - B B -! C
k 1
in which rate constants k1 ; k 1 are much larger than the rate constants k2 , so the ®rst
(a) Apply the method discussed in Chapter 5 to obtain differential equations for cA
and cB valid on the slow time scale. What are the correct initial conditions for
(b) Now consider the problem using the singular perturbation approach of Exer-
cises 5.21 and 5.22. Assume k 1 is the largest rate constant and de®ne the fast
kt K1 kk1
time scale as
1 1
Expand cA and cB on the fast time scale as
cA X0 X1 2 X2 cB Y0 Y1 2 Y2
Show the fast time-scale model for zero order is
dX0 K X Y X0
cA0
d 1 0 0 0
dY0 K X Y Y0
c B 0
d 1 0 0 0
X0 cA0 e 1K1
c 1K1
A0 cB 0
e
1
K1 1
1
cA A0 A1 2 A2 cB B0 B1 2 B2
Show the zero order terms give the single algebraic equation
0 K1 A 0 B 0
and the ®rst order terms give the two differential equations
dA0 K A B
dt 1 0 1
dB0 K A B k B
dt 1 0 1 2 0
Use the matching conditions with the fast time scale solution to derive the initial
A0
cA0 cB0
K
0
1
B0
cA0 cB0
K1K
0
1 1 1 1
(e) We now wish to eliminate the ®rst-order term B1 appearing in the zero-order
B1 1
. Use the algebraic equation and show the two zero-order differential equations
can be rewritten as
dA0 k2 K1 A
dt 1 K1 0
dB0 k2 K1 B
dt 1 K1 0
(f) Solve these differential equations and show the zero-order outer solution is
k2 K1
t
cAo cA0 cB0
K e 1
1K1
1 1
K
cBo cA0 cB0
K e
1
k2 K1
1K1 t
1 1
Compare this solution to the reaction equilibrium solution of the ®rst part. What
1 Notice that this addition of the two equations to remove the unknown B1 is analogous to the step that
k1
E S )*
- - ES
k 1
k2
ES - ! P E
in which the free enzyme E binds with substrate S to form bound substrate ES in the
®rst reaction, and the bound substrate is converted to product P and releases free en-
zyme in the second reaction. This mechanism has become known as Michaelis-Menten
kinetics [6], but it was proposed earlier by Henri [4]. If the rates of these two reactions
are such that either the free or bound enzyme is present in small concentration, the
of S to P
S - ! P r
(a) For a well-stirred batch reactor, show the total enzyme concentration satis®es
(b) Find an expression for the QSS concentration of E. What is the corresponding
concentration of ES?
(c) Show the rate expression for the reduced model's single reaction is
r kcKc
S k k2 KE0 K k k1 k E0 cE
cES
0 0 (5.163)
1 S 1 2
which depends solely on the substrate concentration. The inverse of the constant
RS r RP r
Notice we have reduced the number of reactions from two to one; we have reduced
(d) Plot the concentrations versus time for the full model and QSSA model for the
k1 5 k 1 k2 1 10
cE
0 1 cES
0 cS
0 0 50 cP
0 0
28 Chemical Kinetics
k1
E S )*
- - ES
k 1
k2
ES - ! P E
Now assume the rate constants satisfy k1 ; k 1 k2 so that the ®rst reaction is at
RP kcKS c
e
ke k2 K1 E0 K1 k1 =k 1 (5.164)
1 1 S
in which K1 is the equilibrium constant for the ®rst reaction. Notice this form is
identical to the production rate of P given in the QSSA approach. For this reason,
these two assumptions for reducing enzyme kinetics are often mistakenly labeled
considered a special case of the QSSA [3, p.97]. That viewpoint is erroneous be-
It is interesting to note that in their original work in 1913, Michaelis and Menten
in which the second step is slow compared to the ®rst step [6]. Michaelis and
Menten credit Henri with proposing this mechanism to explain the experimental
observations that (i) production rate of P increases linearly with substrate at low
The QSSA analysis of enzyme kinetics was introduced by Briggs and Haldane in
1925, in which the enzyme concentration is assumed small compared to the sub-
strate [1]. Since that time, the QSSA approach has become the more popular
of Henri, Michaelis and Menten, and Briggs and Haldane, and includes English
translations of Henri's 1901 paper in French[5] and Michaelis and Menten's 1913
paper in German [7]. Boyde concludes that Henri's pioneering contribution has
been largely overlooked in favor of Michaelis and Menten, despite Michaelis and
The reader should be aware that either approximation may be appropriate de-
pending on the values of the rate constants and initial conditions. Although both
reduced models give the same form for the production rate of P, they are often
quite different in other respects. Finally, for some values of rate constants, in
assumption apply.
(c) Show that the slow-time-scale reduced model for the reaction equilibrium assump-
ES - ! E S re1
!
S - P re2
with the following rate expressions
!
K 1 cE
r K c c
e1
kc
e S
ke k2 K1 E0
1 E
1 S 1 K1 c S
re2 kcKS c
e
K1 k1 =k 1
1 S
1
Notice here we have not reduced the number of reactions; we still have two re-
actions, but as before we have reduced the number of rate constants from three
( k1 ; k 1 ; k2 ) to two ( k;
e K1 ). The ®rst rate expression here depends on cS and cE
rather than only cS as in the previous QSSA reduction. Therefore the production
(d) Plot the concentrations versus time for the full model and reaction equilibrium
model for the following values of the rate constants and initial conditions.
k1 : 0 5k 1 k2 : 1 0 5
cE
0 20 cES
0 cS
10 cP
0 50 0 0
Recall that you must modify the initial conditions for the slow-time-scale model
k1 k2 k3
A - ! B B - ! C 2B - ! D
(a) Write the material balances for all four species: A, B, C, and D.
(b) Assume species B is a highly reactive intermediate species. Write down the alge-
(c) Solve this equation for intermediate B's concentration in terms of concentrations
(d) Remove the concentration of B from the material balances of reactants and prod-
ucts and write the reduced model for the time evolution of the concentrations of
A, C and D.
You have been continuing studies on the following reaction and have noticed limitations
O OH
2
O
k 2
2 H O
cyclohexanone phenol
or
A 2
O !k
- C 2H O 2 r kcA cO2
Your chemist colleague does some experiments and comes up with the following mech-
anism consisting of six reactions involving two different types of catalyst sites S and S
k1
)-*
AS OS - BS 2
H OS
O O OH
k1
BS OS
-k!2 CS 2
H OS
K3
A S )*
- - AS
K
O2 2 S
)*4
- - 2 OS
K5 cyclohexanone cyclohexenone phenol
C S )*
- - CS A B C
K6
H O 2 S
)*
- -
2
H OS
Chemical Kinetics 31
1. Reaction 2 is irreversible.
5. O 2 2
and H O adsorb on surface sites S
due to their small molecular sizes.
You have a measurement device that is calibrated to detect the total concentration
of phenyl groups in the system, both in the gas phase and on the surface. What is
the production rate of phenyl groups in the above mechanism in terms of gas-phase
sites cm cm
(for S) and
(for S )? No surface species concentration should appear in
k1
CO S )*
- - CO ads
k 1
k2
O 2 2S )*
- - 2O ads
k 2
k3
CO ads O ads !
- CO 2 2S
(a) Express the full kinetic model for these reactions taking place in a well-mixed
S, the vacant surface sites. Solve the full model using the initial condition of an
equimolar concentration of CO and O in the gas phase and a bare catalyst surface.
Use the following parameter values (expressed in some consistent set of units)
cCO
; cO
; c CO c O ; c v c m
0 1 0 1 0
c m ; k1 k2 : ; k 1 k 2 : ; k3
1 0 01 0 05 1
Solve the resulting ODE model with Octave or MATLAB . Plot all species concentra-
(b) Next increase the rate constants of the ®rst two reactions by two orders of mag-
nitude
k1 k2 ; k 1 k 2 ; k3
1 5 1
32 Chemical Kinetics
and solve the model again. Does the solution indicate that it is reasonable to
(c) Derive the reduced model that eliminates the ®rst two fast reactions' rate con-
K1 k1 =k 1 ; K2 k2 =k 2
Note that your reduced model will be a mixture of differential equations and
algebraic equations. Use an implicit ODE solver to solve this model. See the com-
or MATLAB .
Hint: note that before solving the reduced model, you need to adjust the initial
conditions by equilibrating the two fast reactions. This step requires that you
tions, and use these equilibrated concentrations as the initial conditions for the
ceq c
T "
0
We are keeping the slow third reaction extent "3 ®xed at zero, so we have two
unknowns, "1 ; "2 . We have two constraints, the equilibrium conditions for the
fast reactions,
c CO K1 cCO c v ;0 c O K2 cO c 2v 0
Solve the reduced model using ceq as the initial condition. Does your reduced
model agree with the full model? Explain why or why not.
Revisit Exercise 5.28, but this time let the reactions takes place in a constant volume
CSTR. Assume the reactor is initially charged with inert N 2 and a bare catalyst surface,
the feed is equimolar concentrations of CO and O , 2 cCOf cO2 f 1, the residence time
of the CSTR is 10. Assume that the inert N 2 feed concentration is large enough that
we may assume Q Qf without signi®cant error (even though moles are not conserved
in these reactions).
Consider the process depicted in Figure 6.43 in which the reactor contents and the
streams entering and leaving the reactor consist of multiple, well-mixed phases. With-
out loss of generality, we shall consider two phases: and . Consider the state of
extensive variables
T ; P ; m ; cj j ; : : : ns
1
T ; P ; m ; cj j ; : : : ns
1
Let us assume the phases equilibrate with each other even though the reactor is
conditions of phase equilibrium, Equations 3.50, imply that both phases are at the
same temperature, T , and both phases are at the same pressure, P , which allows us
to describe the state of the reactor with 2 ns 2 intensive variables and 2 extensive
1 JBR would like to thank Vasilios Manousiouthakis of UCLA for helpful discussion of this and the next
two exercises.
QÇ
T P VR
cj m Q
Qf
cjf
Tf cj m Q
33
34 The Energy Balance for Chemical Reactors
variables
T; P; m ; m ; cj ; cj j ; : : : ns1
Since the process in Figure 6.43 is an open system, we also have the two ef¯uent
Q ; Q ; Q Ç
We assume we have speci®ed the feed conditions, Qf ; cjf ; Tf . Let us also assume that
the the system pressure and volume are known constants and the rate of heat trans-
ferred to the system is speci®ed. So we have speci®ed two of our unknown variables
P; QÇ
which brings us back to 2 ns 5 unknowns, and we have one constraint on the total
system volume.
(a) Is the process now fully speci®ed? If so, write the 2 ns 5 equations that fully
(b) Let's apply this result to a simple system consisting of a single species that does
not undergo chemical reaction. An example would be a boiler fed with a pure
water stream in which liquid water and steam are the two different phase ef¯uent
Can you write seven independent equations to specify this system? If so, write
If not, provide a physical argument for the boiler that shows the system is under-
Exercise 6.25: Writing DAE models for systems with multiple phases
(a) Describe the system of Exercise 6.24 with a set of DAEs of the form
dx f x; y
dt
g x ; y
0
in which x are the differential states and y are the algebraic states. Be sure that
have access to complete thermochemical data, which means you can evaluate any
Notice that neither the heat of reaction nor the heat of phase change are required in
the energy balance when one has complete thermochemical data available. These two
enthalpy changes are only useful for simplifying the required thermochemical data.
When full thermochemical data are available on the other hand, using heats of reaction
and phase change merely complicates the energy balance, and their use should there-
fore be avoided. Experience indicates that always using heat of reaction and heat of
phase change can be a dif®cult habit for chemical engineers to break. It seems new
training on how to express energy balances would be useful given the current wide
state. Say a given database has de®ned a particular reference pressure, temperature,
composition and phase of matter, which determines all intensive thermochemical prop-
erties. Denote the speci®c internal energy and enthalpy of the reference state by
Uref Href
à Ã
ref
Consider these common material and energy balances for a two-phase system
dnj Q c Q c Q c X r r j ; : : : ; n
f jf j j ij i i s
dt i
1
dU P dVR Q H Q H Q H Q
dt dt f f f à à à Ç
Ue U Uref
à à Ue U Uref
à à He H Href
à à He H Href
à Ã
Calling the thermochemical database functions returns these shifted quantities. De®ne
Ue m U Uref
m U Uref
à à à Ã
m U m U m m
Uref
à à Ã
Ue U mUref Ã
(a) We would like to use a new energy balance expressed in terms of these shifted
dUe P dVR Q He Q He Q He Q
dt dt f f f Ç (6.95)
(b) Consider a second approach in which the same reference energy, Ã Eref is subtracted
Ue U Eref
à à Ue U Eref He H Eref
à à à à He H Eref
à Ã
Ue U mEref Ã
Is the shifted energy balance, Equation 6.95, correct with these de®ned quantities? If
In deriving the energy balances in Chapter 6, we often found enthalpy rather than inter-
nal energy to be the more convenient energy function. For situations in which internal
energy is more convenient, provide derivations for the following useful relations
@U TV PV
(a)
T
@P T;nj
(b)
!
@U H j TV j
@nj T;V;nkj T
The reaction
A )*
B - - C D
shown in Figure 6.44. The feed is equimolar in A and B and enters each reactor at 27 C.
The heat removed between the reactors is :
87 5 kcal/min.
T0
27 C T0 27 C
Q Ç :
87 5 kcal/min
State any assumptions that you make while solving the problem.
Data:
HR 30 kcal/mol
Cp 25 cal/mol K
K
at 50 C :
5 0 10
5
NA0 NB0
, 10 mol/min
Q
Ç :
87 5 kcal/min
R :
1 987 cal/mol K
Exercise 6.29: A colleague's batch reactor energy balance
A colleague presents you with the following derivation of the energy balance for a
You know, the people writing our textbooks sure make things complicated. I
found a much simpler energy balance for the constant-volume batch reactor.
dU Q
dt
Ç (6.96)
in T V
, and nj by
!
dU @U dT @U dV X @U dnj
@T V;nj @V T;nj j @nj T;V;nkj
(6.97)
The de®nitions of constant volume heat capacity and partial molar internal
energy are
CV @U@T V;nj (6.98)
!
@U
U j @n (6.99)
j T;V;nkj
38 The Energy Balance for Chemical Reactors
Since the reactor volume is constant, dV 0, and forming the time derivative
dU C dT X U dnj
V
dt dt j j dt (6.100)
H j U j PV j (6.101)
in which Hj is the partial molar enthalpy and Vj is the partial molar volume.
dnj R V X r V ; j ; : : : ; ns
j R ij i R
dt i
1 (6.102)
in which ri is the reaction rate per reactor volume and ij is the stoichio-
ns
X ns
X
HRi ij H j VRi ij V j (6.103)
j 1 j 1
I put the material balance and these de®nitions in Equation 6.100 and I obtain
CV dT
X
dt HRi P VRi
ri VR Q Ç (6.104)
Notice that my result doesn't agree with Equation 6.66 in our text!
CV dT HRi T VRi ri VR Q
X
dt
Ç
i T
(a) Is Equation 6.96 correct for the constant-volume batch reactor given the other
assumptions?
(c) Do Equations 6.98 and 6.99 agree with the usual de®nitions of constant-volume
(f) Is Equation 6.102 a correct statement of the material balance for the batch reactor?
(g) Does Equation 6.103 agree with the usual de®nitions of enthalpy change and vol-
If your answers to all questions above were yes, do you think the two energy
balances are equivalent? If so, how would you show this equivalence?
If your answers to some of the questions above were no, how would you repair
Start with the plug-¯ow energy balance of Section 6.5, rewritten here,
@ U
@ QH
q
@t @V
à à Ç
(a) Derive the dynamic energy balance in the variables T; P; cj for a single-phase sys-
tem.
2
(b) Does your dynamic energy balance reduce at steady state to Equation 6.50? Dis-
ple 6.6.
(a) Solve the resulting boundary-value problem for the parameters given in Table 6.6
with collocation and compare your three steady-state pro®les to the results given
in Figures 6.39 and 6.40. Be sure to note that the ¯owrate should be Qf
: 3
0 05713 m /s in Table 6.6, and that the the middle steady-state pro®le (B) in
Figures 6.39 and 6.40 is shifted somewhat from the correct answer. See the er-
rect pro®les.
CP @T T @P Q C @T @P
T
@V
X
HRi ri q
@t
Ã
@t P @V Ã 1 Ç
i
40 The Energy Balance for Chemical Reactors
(b) What initial guesses for the collocation equations did you use to converge to the
(c) Vary the number of collocation points and demonstrate convergence for the upper
steady-state solution. How many collocation points are required for a reasonably
accurate solution?
(d) Describe the advantages and disadvantages of the shooting method versus collo-
k1
)*
A - - B
k 1
takes place in two adiabatic PFRs with interstage cooling as depicted in Figure 6.33. In
Example 6.4 the PFRs were assumed to be long and the reaction reached equilibrium at
Q
Ç ;
200 000 BTU/hr
Vt : 2 22 ft
3
k1 :
3 47 s
1 at T Tf Ea =R 9000
R
(a) Solve for the NA and T pro®les in both reactors. Plot NA and T versus reactor
volume for both reactors. Compute the conversion at the exit of each reactor
and compare to the answers given in Figure 6.34. You should see that the overall
(b) Next let's see if we can achieve xA2 : 0 7 by changing the heat duty. Vary Q
Ç
between ;
125 000 and ;
250 000 BTU/hr, and make a plot of xA2 Q
versus Ç over
this range. Can you reach xA2 : 0 7 by changing the heat duty? If so, what value
Q
of Ç do you use? If not, why not?
The Energy Balance for Chemical Reactors 41
k1 k2
A - ! B B -! C
starting with no B or C. It was shown previously that for an isothermal, constant volume
cB cA0 k k1 k e
batch reactor
c A c A0 e k1 t k1 t e k2 t
2 1
We wish to maintain an isothermal reactor, but the reactions are quite exothermic.
To size the heat exchange equipment appropriately, we must know how much heat to
Q
(a) Give an expression for Ç , the heat removed by the heat exchanger. You may use
(b) At what time is the maximum heat generated in the reactor? Do not assume
How much heat should the heat exchanger be able to remove? When does the
maximum occur? If you were not able to answer part (b) analytically, estimate the
maximum heat removal rate and time of the maximum for the parameter values
given.
(d) After the reaction has gone for 3 min, how many moles of A have reacted? Moles
A !kB - C
takes place in a CSTR equipped with a heat exchanger. The feed consists of A and B
dissolved in an inert solvent. An inhibitor is stripped from the feed just before addition
to the reactor.
(a) You plan to start up the CSTR with a solvent ®lled reactor at the feed temperature.
Simulate the reactor behavior ( cA t
; cB t
; Tt
) for six hours after start up.
42 The Energy Balance for Chemical Reactors
(b) Simulate the reactor again for the following situation. After running for one hour,
the coolant pump fails for 2 1/2 hours (the reactor runs adiabatically), and then
the pump is restored. Simulate the response of the reactor for the six-hour start
up period. Does the system return to the steady state observed in part (a)? Explain
(c) Calculate the adiabatic temperature rise for the feed. Start the reactor ®lled with
solvent with initial temperature equal to the feed temperature plus the adiabatic
temperature rise of the feed. Simulate for six hours. Compare the steady state
The reaction parameters are listed in the following table in which kT
km e E1=T 1=Tm
and U A
= C
à V
.
p R
km 5 10
4 min
1
Tm 298 K
E 8000 K
HR 3 5
10 kJ =
kmol
1000 kg = 3
m
CP
à 4 0 : kJ =
kg K
:
0 015 min
1
Ta 298 K
Tf 298 K
cAf 3 0 : kmol = m
3
cBf 2 0 : kmol = m
3
30 min
7
Fixed-Bed Catalytic Reactors
pellet is DA : 2
0 009 cm /s at 525 K and that diffusion is in the regime of Knudsen
¯ow. The bulk ¯uid and the external surface concentrations can be assumed the same.
Find the production rate if the catalyst is changed to a cylindrical pellet that is 0.5 cm
®xed-bed reactor:
A - ! B r kcA2
The intrinsic rate constant is 16 L = mol s. The feed is 0 8 L : = s of a solution of A with
bed density is 0 5 g : = 3
cm . The effective diffusivity of A inside the catalyst pellet is
:
3 14 10
3 cm2 =
s. You wish to achieve 99.9% conversion of A.
(a) First neglect both diffusional resistance and mass transfer resistance. What is the
(b) Next consider diffusional resistance but neglect mass transfer resistance. What
are the Thiele modulus and effectiveness factor at the entrance of the reactor?
What are the Thiele modulus and effectiveness factor at the exit of the reactor?
Given these two effectiveness factors, compute an upper bound and a lower bound
(c) Next consider both diffusional resistance and mass transfer resistance. Assume
43
44 Fixed-Bed Catalytic Reactors
the effectiveness factor at the entrance of the reactor. Find the effectiveness factor
Given these two effectiveness factors, compute an upper bound and a lower bound
(d) Compare your answers to parts (a)±(c). Discuss whether diffusion and/or mass
A - ! ; B r kcA1=2
Assume this reaction takes place in the the PFR in Figure 8.44, which is operating at
steady state with residence time and feed A concentration cAf . The total overall
reaction rate is r1 , which is the moles of A converted per time in the entire reactor
volume, divided by the volume of the reactor. At time zero, we turn on a stirrer and
perfectly mix the reactor contents. Immediately upon mixing the reactor achieves an
overall reaction rate r2 . Also consider the steady-state overall reaction rate achieved at
your chemistry colleague comes to you with a catalytic breakthrough. By using a ho-
mogeneous catalyst the reaction can now tolerate various substituents on the cyclohex-
anone ring
PFR
Figure 8.44: Reactor in three conditions: PRF at steady state, just after a perfect stirrer
is started, and stirred reactor at steady state. The overall reaction rates
45
46 Mixing in Chemical Reactors
OH O
R5 R1 R5 R1
O2 2
2 H O
R4 R2 R4 R2
R3 R3
This new reaction has a downside though. The substituted cyclohexanones and
phenols can no longer tolerate high temperatures and must remain in the liquid phase,
but the O 2 is still in the gas phase. These new issues now need to be understood.
Assume the reaction rate is given by r kcA when excess O 2 is present. To provide
the required O 2 the reaction takes place in an isothermal, isobaric tubular membrane
reactor inside a large cylinder with a continuous supply of O 2 as shown in the following
schematic.
2
O supply
VR ; Qf
66 mL = 3 mL min ; cAf : = ; cA VR
:
0 1 mol L 0 01 mol =;
L
(a) Ideal PFR. What is the residence time for this reactor? Solve the steady-state PFR
(b) Nonideal ¯ow. To perform an RTD step test on the reactor, two inert solvents
are used, toluene and cyclohexane. The reactor (PFR) begins with toluene ¯owing
ane ¯ow and cutting off the toluene ¯ow. The measured RTD is shown in Figures
(c) Dispersion number. From Figure 8.47, which is similar to Figure 8.10 in the
text, and the experimental data shown in Figures 8.45 and 8.46, estimate the
(d) Dispersed PFR model. What differential equations would you use to model the
dispersed PFR? Include all required boundary conditions. Make sure that all pa-
How would you solve this model? Describe what you would do; you do not need
to actually solve it. Do you need any additional information to solve the problem?
If so, what?
(e) Engineering intuition. Would you expect the outlet A concentration to be higher,
lower, or the same as the ideal PFR case? Assume both cases have the same mean
residence times, ¯ow rates, rate constants, temperatures, pressures, etc. and base
your argument on only the nonideal ¯ow and dispersion. Justify your answer.
48 Mixing in Chemical Reactors
0.9
0.8
0.7
0.6
0.5
xc
0.4
0.3
0.2
0.1
-0.1
time (s)
Figure 8.45: Mole fraction of cyclohexane in the ef¯uent versus time after step change
at t 0.
0.9
0.8
0.7
0.6
0.5
xc
0.4
0.3
0.2
0.1
-0.1
time (s)
Figure 8.46: Mole fraction of cyclohexane in the ef¯uent versus time after step change
:
0 9
:
0 8
:
0 7
:
0 6
P
:
0 5
:
0 4
:
0 3
:
0 2 D :
0 0001
D :
0 0005
:
0 1 D :
0 001
D :
0 005
0
:
0 8 :
0 85 :
0 9 :
0 95 1 :
1 05 :
1 1 :
1 15 :
1 2
=
Figure 8.47: P
versus = for dispersed plug ¯ow with different dispersion number
D.
9
Parameter Estimation for Reactor Models
Exercise 9.16: Least squares estimation and estimating measurement variance from
data
H O2 )*C0 - - CO 2 H 2 (9.59)
H O2 )* H - - H 2 OH (9.60)
2 3
RH
6 R H2
6 7
" # 7
6 7
r rr1 R 66 RROH
6 7
7
7
2 6 H2 O 7
6 7
4 RCO 5
RCO2
We have replicate measurements of the six species production rates. We wish to es-
timate the two reaction rates. Twenty replicate measurements are given on the class
website in the Exams folder. I also list the actual variance of the measurement error,
Qac . Download the ®le confestvar_measurements.dat and load these into Octave
using the load command; you should obtain two matrices, Qac and Rmeas. Perform the
following analysis.
(a) Assume we do not know Qac (the usual case) so we assume the measurement
all 20 estimates, and the 95% con®dence interval with center located at the mean
estimate.
Notice you should use Qac to compute the con®dence interval, not QI . We
want to analyze how the estimates are actually distributed, not how they would
(b) Next calculate the covariance matrix for the twenty measurement samples. Call
51
52 Parameter Estimation for Reactor Models
(c) Next using Q Qest , reestimate the values of r using weighted least squares. Plot
all twenty estimates and the 95% con®dence intervals with center located at the
mean estimate. How different are the mean estimates using QI and Q Qest ?
(d) Finally, using zero as the center, plot the 95% con®dence intervals for estimates
using Q I Q Qest
, , and Q Qac . How do these ellipses compare in size? Do
you obtain tighter con®dence interval (smaller variance) when you use Qest instead
of I in the least-squares estimate? How does the con®dence interval using Qest
compare to the one using Qac ? Why?
Bibliography
[1] G. E. Briggs and J. B. S. Haldane. A note on the kinetics of enzyme action. Biochem. J., 19:
338±339, 1925.
[3] G. G. Hammes. Thermodynamics and Kinetics for the Biological Sciences. Wiley-Interscience,
919, 1901.
[5] M. V. Henri. General theory of the action of certain enzymes. In T. R. C. Boyde, editor,
Foundation Stones of Biochemistry. Voile et Aviron, Hong Kong, 1980. English translation.
[6] L. Michaelis and M. L. Menten. Die Kinetik der Invertinwirkung. Biochem. Z., 49:333±369,
1913.
[7] L. Michaelis and M. L. Menten. The kinetics of invertin action. In T. R. C. Boyde, editor,
Foundation Stones of Biochemistry. Voile et Aviron, Hong Kong, 1980. English translation.
[8] D. L. Nelson and M. M. Cox. Lehninger Principles of Biochemistry. Worth Publishers, New
53