Chapter 2 Rate Data Analysis
Chapter 2 Rate Data Analysis
UNIVERSITY OF TECHNOLOGY
FACULTY OF CHEMICAL ENGINEERING
Chapter 2
Interpretation of Batch Reactor Data
from rate data => rate equation
• Single reaction: When a single stoichiometric equation and single rate equation
are chosen to represent the progress of the reaction.
• Multiple reactions: When more than one stoichiometric equation is chosen to
represent the observed changes, then more than one kinetic expression is
needed to follow the changing composition of all the reaction components
– Series reactions:
– Parallel reactions:
Types of reaction
Phản ứng cơ bản: phương trình tốc độ tương ứng với phương trình cân bằng hóa họcđược gọi là.
Phản ứng không cơ bản: không có sự tương ứng trực tiếp giữa cân bằng hóa học và tỷ lệ.
Order of Reaction
• For rate expressions not of the above form, it makes no sense to use the term
reaction order.
Conversion
Batch operation
A +B
Continuous operation A là limited component khi A complete
but B not complete, B: remaining
Constant volume
Integral Method of Analysis of Data
General procedure:
• Puts a particular rate equation to the test by integrating and comparing the
predicted C versus t curve with the experimental C versus t data.
• If the fit is unsatisfactory, another rate equation is guessed and tested.
• The integral method is especially useful for fitting simple reaction types
corresponding to elementary reactions.
• Linear regression is often used
The measured data can be inllustrated in a ln(CA/CA0) vs. T diagram. The slope
of the straight line leads to the reaction rate constant (Figure 2.3).
2.6 of 2.3
1. Rates of reaction
1.1. Description of reaction rates
Reaction rates depend usually in a complex manner on the concentrations, the
temperature and often on the effect introduced by catalysts:
r = f ( ci ,T, catalyst)
• First order:
A → Products: r = k(T)CA k in [s-1]
• Second order:
A + B → Products: r = k(T)CACB k in [m3 mol-1s-1]
2A → Products: r = k(T)CA2 k in [m3 mol-1s-1]
Chapter 2. Interpretation of Batch Reactor Data
Reversible Reactions:
• First order:
A R: r = k1CA – k2CR k in [s-1]
• Second order:
A + B R + S: r = k1CACB – k2CRCS k in [m3 mol-1s-1]
2A 2R: r = k1CA2 – k2CR2 k in [m3 mol-1s-1]
2A R + S: r = k1CA2 – k2CRCS k in [m3 mol-1s-1]
A + B 2R: r = k1CACB – k2CR2 k in [m3 mol-1s-1]
The dependency of the reaction rate on the temperature can be described using
Arrhenius law:
• EA is the activation energy of the reaction.
EA
k (T ) = k0 exp − • ko is the pre-exponential factor (not dependent on
RT the reaction temperature
Chapter 2. Interpretation of Batch Reactor Data
− ln A = kt
C A0
Chapter 2. Interpretation of Batch Reactor Data
C A = C A0 e − kt (2.4)
C A0 − C A
With XA = , the eq. (2.3) leads to: (2.5)
C A0
Y= − ln(1 − X A ) = kt or X A = 1 − e − kt (2.6)
Noting that the amounts of A and B that have reacted at any time t are equal and
given by CA0XA, eq. (2.7) can be written in terms of XA as
r= (2.8)
(2.9)
(2.10)
Chapter 2. Interpretation of Batch Reactor Data
After breakdown into partial fractions, integration, and rearrangement, the final
result in a number of different forms is
(2.11)
r= (2.12)
(2.13)
Chapter 2. Interpretation of Batch Reactor Data
(2.17)
(2.18)
(2.19)
Chapter 2. Interpretation of Batch Reactor Data
The k values are found using all three differential rate equations. First of all,
eq. (2.17), which is of simple first order, is integrated to give
(2.20)
(2.21)
(2.22)
Chapter 2. Interpretation of Batch Reactor Data
(2.23)
(2.24)
(2.25)
Chapter 2. Interpretation of Batch Reactor Data
(2.26)
(2.28)
Chapter 2. Interpretation of Batch Reactor Data
Noting that
(2.29)
Thus, in general, for any number of reactions in series it is the slowest step that has
the greatest influence on the overall reaction rate.
Chapter 2. Interpretation of Batch Reactor Data
(2.32)
(2.33)
Chapter 2. Interpretation of Batch Reactor Data
Eq. (2.30)
(2.34a)
(2.34b)
(2.35)
Chapter 2. Interpretation of Batch Reactor Data
(2.36)
(2.37)
Taking integration of the above equation, the following relation can be found:
(2.38)
Chapter 2. Interpretation of Batch Reactor Data
(2.39b)
(2.39c)
(2.39d)
with the restrictions that CA0 = CB0, and CR0 = CS0 = 0, the integrated rate
equations for A and B are all identical, as follows
(2.40)
Integral Method of Analysis of Data
linear requestion method
assume the rate equation
2. Determination of kinetic parameters -rA = kCA =-dCA/dt volume is
2.1. Integral method rate equation: unchange
log(rA)
Figure 2.4. Test for the particular rate Figure 2.5. Test for an nth-order rate
form rA = kf(CA) by the differential method by the differential method.
form.
Chapter 2. Interpretation of Batch Reactor Data
Problems
1. Aqueous A at a concentration CA0 = 1 mol/liter is introduced into a batch reactor
where it reacts away to form product R according to stoichiometry A → R. The
concentration of A in the reactor is monitored at various times, as shown below:
We plan to run this reaction in a batch reactor at the same catalyst concentration
as used in getting the above data. Find the time needed to lower the concentration
of A from CA0 = 10 mollliter to CAf = 2 mol/liter.
Pressure dependence
with n = r+ s+ … – a –b – …
pA N A N A0 − a. x N A0 a N − N0
CA = = = = −
RT V V V n V
a
hay p A = C A RT = p A0 − (P − P0 ) (2.3)
n
r
cho R p R = C R RT = p R0 + (P − P0 ) (2.4)
n
In which Po: Initial total pressure
P: Total pressure at t
pAo: Initial partial pressure of A
pRo: Initial partial pressure of R
34 2/24/2022 Chương 2. Xử lý dữ kiện động học
Varying-volume Batch Reactor
1 − 𝑋𝐴
𝐶𝐴 = 𝐶𝐴0
(1 + 𝜀𝐴 𝑋𝐴 )
Varying-volume Batch Reactor
1 dN i 1 d(Ci V) 1 VdCi + Ci dV
ri = = =
V dt V dt V dt
dCi Ci dV ∆𝑉
ri = + 𝑋𝐴 =
𝑉0 𝜀𝐴
dt V dt
V = V0 (1 + A X A )
VX A = 1 − VX A = 0
A = → expansion factor
VX A = 0
36 2/24/2022 Chương 2. Xử lý dữ kiện động học
Varying-volume Batch Reactor
a. 2A + B -> R
Feed mixture: 40%A, 40%B, 20% inert gas
Assume Vo = 1
2A + B -> R inert gas
t=0 0.4 0.4 0 0.2
reacted 0.4 0.2 0.2 0.2
Remain 0 0.2 0.2 0.2 (XA = 1)
0.2+0.2+0.2 −1
𝜀𝐴 = = -0.4
1
B .A + B -> 2R + S
Feed: 20%A, 60%B, 20% inert gas
0.4+0.4+0.2+0.2 −1
𝜀𝐴 = = 0.2
1
N A0d(1 − X A )
(− rA ) = − 1 dN A
=−
1
V dt V0 (1 + A X A ) dt
C A0 dX A
=
1 + A X A dt
N A0d(1 − X A )
(− rA ) = − 1 dN A
=−
1
V dt V0 (1 + A X A ) dt
C A0 dX A
=
1 + A X A dt
𝑋𝐴
(1 + 𝜀𝐴 𝑋𝐴 )𝑛−1 𝑛−1
න 𝑛
𝑑𝑋𝐴 = 𝑘𝐶𝐴0 t
0 (1 − 𝑋𝐴 )
For all rate forms other than zero-, first-, and second-order the integral
method of analysis is not useful.
• Differential method
• Integral method
• Half-lives method
• Initial rate method
• Differential reactor
• More complex kinetics
2.1.2. Differential Method of
Analysis of Data
Procedure:
1. Plot the CA vs. t data, and then by eye
carefully draw a smooth curve to
represent the data. This curve most likely will
not pass through all the
experimental points.
2. Determine the slope of this curve at
suitably selected concentration values.
These slopes dC/dt = rA are the rates of
reaction at these compositions.
Differential Method
dC j
( )
d d
Fj0 − Fj + rj V = N j → rj V = C j V → rj V = V
dt dt dt
0 0
dC j −dC A
= rj Where –rA = kCAa = kC Aa
dt dt
Average slope
a) Plot CA/t as a function of t
b) Determine dCA/dt from plot by graphical or numerical methods
c) Plot ln(-dCA/dt) vs ln CA
−dCA − dCA dt
ln = lnk + a lnCA Slope = a k=
dt C Aa
To find k, find the value of –dCA,p/dt that corresponds to a specific
concentration CA,p. Raise CA,p to the a power and divide into –dCA/dt)p
Hey, we just jumped from step a to step c. How do we get dCA/dt?
Graphical Method
C A
−
1. Plot CA/t vs t t
2. Draw rectangles on the −
dC A
dt t
(−CA/t)t=0
graph. Then draw a curve so
(−CA/t)t=1
0
that the area above the curve
that is cut off of each dC A
rectangle approximately fills −
dt t
(−CA/t)t
1
the unfilled area under the dC A
−
curve dt t
2
3. dCA/dt is read using the value 0 t1 t2
where the curve crosses a
specified time
Graphical Method Example
-
Graphical Method Example
CA t t CA −CA/t -dCA/dt
8 0 1 -4 4 4.5
4 1 1 -2 2 2.55
2 2 1 -1 1 1.35
1 3 0.5
−dC A
= kC Aa
dt
−dCA
ln = lnk + a lnCA
dt
Plot ln(-dCA/dt) vs ln CA
-
Slope = a
− dCA dt
k=
C Aa
Graphical Method Example
CA t t CA −CA/t -dCA/dt ln(-dCA/dt) ln(CA)
8 0 1 -4 4 4.5 1.5 2.1
4 1 1 -2 2 2.55 0.9 1.4
2 2 1 -1 1 1.35 0.3 0.7
1 3 0.5 -0.7 0
Plot ln(-dCA/dt) vs ln CA
−dCA
ln = lnk + a lnCA
dt
Slope = a = 1.0
− dCA dt
k=
C Aa
4 .5
→k = = 0 .6 -rA= (0.6/time)CA
81
Graphical Method Example
CA t t CA −CA/t -dCA/dt −dCA
ln = lnk + a lnCA
8 0 1 -4 4 4.5 dt
4 1 1 -2 2 2.55 Slope = a
2 2 1 -1 1 1.35 − dCA dt
k=
1 3 0.5 C Aa
Differential Method
Only for irreversible reactions
Advantages: 1 experiment
-
✓ Differential method
✓ Integral method
• Half-lives method
• Initial rate method
• Differential reactor
• More complex kinetics
Method of Half-lives
dC 1 1 1
A → products − A = kCAa t= a −1 −
dt k (a − 1) CA CA0a −1
−rA = kCAa
1
ln (t1/2) CA = CA0 at t = t1 2
2
Slope = 1- a
2a −1 − 1 1
t1 2 =
k (a − 1) CA0a −1
ln CA0
2a −1 − 1
Plot ln(t1/2) vs ln CA0. Get a straight ( )
ln t1 2 = ln
k (a − 1)
+ (1 − a ) lnCA0
line with a slope of 1-α