Gec 534
Gec 534
Gec 534
CHEMICAL
REACTION
ENGINEERING I
FIRST OF THREE
COURSES
Chemical Reaction
Engineering
field that studies
1. the MECHANISMS and RATES of chemical
reactions and
2. the DESIGN of the reactors in which they take
place
RATES: HOW FAST REACTION OCCURS
MECHANISMS: proposed series of steps by which
reaction is assumed to occur at molecular level
DESIGN :TYPE , SIZE and REACTION CONDITIONS
IN THE REACTOR
Chemical Reaction
Engineering
is at the heart of virtually every
chemical process.
Reactor may be small but it is crucial
to the economics of the process
It separates the chemical engineer
from other engineers.
Process Industries Draw Heavily
on Chemical Reaction Engineering
chemical treatment
step
Physical treatment
step
RAW
MAT
RECYLE
OPERATIONS
Heating
Cooling
washing
Distillation
Filtration
Drying
PRODUCTS
MANY CHEMICAL
UBIQUITOUS REACTIONS
OUR BODIES
ENVIRONMENT-decomposition of organics
OZONE LAYER DEPLETION
SMOG
POWER GENERATION
CORROSION
INDUSTRIAL PRODUCTS
DOMESTIC PRODUCTS
**NO
Chemical reaction as a
multidisciplinary subject
ommon integrals
REACTION?
How do you know a reaction has
occurred?
A chemical species is said to have
reacted when it has lost its chemical
identity-i.e. transformed into
another substance
The identity of a chemical species is
determined by
the kind,
number, and
2. COMBINATION
3. ISOMERIZATION
Rate of Reaction
is a function of the following parameters:
1. Concentration of reacting species
2. Temperature of reacting species
3. Pressure( gas phase reactants/ products)
4. the type of catalyst (if any)
Rate of Reaction is independent of the
type of reaction system (batch, plug flow
reactor, etc.)
. rj is an algebraic equation, not a
differential equation
Reaction Rate
EXAMPLE: AB
Classification of reactions
very different ways . For our purpose a broad classification is used
HOMOGENEOUS reactions
I.gas-gas
II.liquid-liquid
.HETEROGENEOUS reactions
I. Gas-solid
II.Gas-liquid
III.Gas-solid-liquid
IV. liquid-liquid
Another classification
Arrhenius equation/law
k=k0e-Ea/RT
where
Ea =activation [energy/mol]
T =absolute temperature of the reaction i.e. K;R 0
R is the universal gas constant.
The units of Ea and R must be consistent. If E a is in
kJ/mol then R must have units of kJ/mol-K
k0 =the frequency or pre-exponential" or "entropy
factor
Equation is most useful relationship between the
specific rate constant and the temperature of
reaction
T1
T2
T3
T4
k1
k2
k3
k4
Lnk
lnk1
lnk2
lnk3
lnk4
1/T
1/T1
1/T2
1/T3
1/T4
slope=Ea/R
1/T
f(Ci )=CA CB CC
where
is the order of reaction with respect to the reactant A
n= + +
n= the overall order of reaction=the sum of the orders of reaction with
respect to each of the reactants
The order of reaction with respect to A reactant is the power to which the
concentration is raised in the rate law expression.
In the above expression, A, B and C are reactants
That is f(Ci ) is the product of the concentration of the reactants each raised
to a power which is the order of reaction
Reaction Mechanisms
The mechanism of a reaction is the series of
elementary steps by which the reaction is
assumed to occur at the molecular level.
It is proposed to give the detailed explanation
at the molecular level on how a reaction
proceeds.
When searching for a mechanism:
Steps are Proposed
Rate equations are derived and
Tested with data
Elementary Reaction
A reaction is said to be
elementary when the orders of
reaction correspond to the
stoichiometric coefficients of the
reactants n the balanced
equation
It must be emphasized that the
orders of reaction are
determined experimentally
Molecularity of a reaction
NO2(g) + ClNO(g)
explains the observed rate laws for both onestep and multi-step reactions.
model assumes that the rate of any step in a
reaction depends on the frequency of
collisions between the particles involved in
that step.
This is the basis for understanding the
implications of the collision theory model for
simple, one-step reactions, such as the
following.
ClNO2(g) + NO(g)
NO2(g) + ClNO(g)
step1
(CH3)3CBr
(CH3)3C+ + Br:Slow step
step2
(CH3)3C+ + H2O
(CH3)3COH2+ :Fast step
step3
(CH3)3COH2+ + OH(CH3)3COH + H3O :Fast step
The overall rate of reaction is therefore more or less
equal to the rate of the first step.
The first step is therefore called the rate-limiting
step in this reaction because it literally limits the
rate at which the products of the reaction can be
formed.
Cl2 2Cl.
Cl. + H2
HCl + H.
H. + Cl2
HCl + Cl.
Cl. + Cl.
Cl 2
Reaction Mechanism A
eg. 2 NO(g) + 2 H2(g) ----> N2(g) + 2 H2O(g)
This appears to be a 4 particle reaction as
written. Impossible!! Hence cant be a single step
reaction
Actual reaction mechanism is :
step #1: 2 NO + H2 ----> N2 + H2O2
slow since it involves a rare three particle
collision.
This is the slowest and therefore is the step
that determines the rate of all the other
steps.
step #2: H2 + H2O2 ------> 2 H2O fast reaction
Overall rate is therefore given by
Reaction Mechanism B
-r NO =k CNO2 CH2
Termolecular reaction steps
require three molecules coming
together at the same time
They are rare because three-body
collisions in the gas phase are rare
4 particle reactions are
therefore virtually
Impossible!!
Batch reactor:
variable volume
Flow reactors
Kinetic Data
Physical Methods:1
Physical methods:2
pH and EMF methods
1. using a glass electrode sensitive to H3O+
. Reactions which occur with change in H3O+ or a
change in [OH-] to be followed with ease.
. A pH meter measures pH directly, and
. a millivoltmeter measures EMFs directly, and these are
related to H3O+
E.g.In CH3COOCH3(aq) +OH-(aq)
+CH3OH(aq)
CH3COO(aq)
Kinetic Data:
Conductance Methods
Fractional Conversion 2
NA0 =Moles of A charged into the reactor initially
(t=0)
NA =Moles of A left in the reactor after reaction
time, t
It follows that the moles of A that has been
transformed into products or consumed after
reaction time, t is
Fractional Conversion3
of the stoichiometric equation: aA + bB
cC + dD
[Moles of A reacted]/[Moles of B reacted]
[NA0 XA ]/ [NB0 XB ]=a/b
In general, the fractional conversion of a
reactant A is different from that of
another reactant B
The fractional conversion of a reactant A
increases with reaction time
At complete conversion, the fractional
conversion of the limiting reactant is 1
*
i is positive for a product
i is negative for a reactant
. If there is more than one reaction
going on, there will be an extent
term for each reaction. So for i
components and j reactions:
Stoichiometric Table
Useful in relating the relationship
between the *Concentrations of
various species in the reaction
mixture with time
*enables easy representation of the
rate law in terms of a single variable
*enables easy calculation of the
partial pressures of various species
for gas phase reactions
*enables easy calculation of the
mole fractions of various species in
charge [2]
remaining, Ni [4]
Ci [5]
NA0= NA0A
to reaction,
[3]
-NA0XA
NA0A -NA0XA
= NA0 (A -XA)
NA0 (A -XA)/ V
= CA0 (A -XA)
NB0= NA0B
-NA0XA(b/a)
NA0B-NA0XA(b/a)
= NA0 [B-XA(b/a)]
NA0 [B-XA(b/a)]/ V
= CA0 [B-XA(b/a)]
NC0= NA0C
NA0XA(c/a)
ND0= NA0D
NA0XA(d/a)
NI0= NA0I
= NA0I
=NA0I /V
Total
N0=NA0(A + B + C + D + I) N=N0 + NA0 XA [(d/a)+ (c/a)-(b/a)-1]
Stoichiometric Table3
Column 1 lists the various components of the
reaction system- reactants, products and inerts
(I)
Column 2 lists the amount of the various
components initially charged to the reactor, all
amounts expressed in terms of NA0 and i0.
The definition of i0, listed at the bottom of Table
2.1 makes this simplification possible. This
column is calculated from the composition of
the initial charge to the reactor.
i = Ni0/ NA0 ;A = NA0/ NA0 =1 ; B = NB0/ NA0
Stoichiometric Table
Column 5 gives the concentration of components at X A.
This is obtained by simply dividing the moles remaining by the
volume of the reaction mixture, which in this case is
assumed to be constant -i.e. C i = Ni / V.
The moles remaining is obtained from column 4.
Moles
remaining, Ni [4]
Concentration
Ci [5]
NA/V
=CA
[ NA0(D +d/a)-(d/a)NA]/V
CD=CA0(D +d/a)-(d/a)CA]
=2(0+1.5)-(1.5)CA]
=3-1.5CA
NI0= NA0I 0
NA0I
NA0I/V
CI= 2.0(0.05)=0.10
Total
N0=NA0(1 + 1 + 0 + 0.05)
=2.05 NA0
2.05x200=410
charge [2]
to reaction, [3]
remaining, Ni [4]
NA0= NA0A
-( ND-ND0)(a/d)
NB0= NA0B
-( ND-ND0)(b/d)
NB0-(ND-ND0)(b/d)
[NB0-( ND-ND0)(b/d)]/V
CB= CB0-(CD-CD0)(b/d)
NC0= NA0C
( ND-ND0)(c/d)
NC0+(ND-ND0)(c/d)
[NC0+( ND-ND0)(c/d)] /V
CC= CC0 + (CD-CD0)(c/d)
ND0= NA0D
( ND-ND0)
ND
ND /V
= CD
NI0= NA0I
NA0I
NA0I/V
CI= CA0I
Ci [5]
Stoichiometric Table:
last comments
It is easy to set up
May be in terms of the fractional conversion
of any reactant
In terms of the concentration of any reactant
in the reaction mixture or
In terms of the concentration of any product
in the reaction mixture
It is easy
Learn to set it up very fast
Very necessary because mathematical
manipulation generates errors and this must
be minimised
1.Integral Method
2.Differential Method
3.Method of initial rates
4.Method of half-lives
5.Method of isolation or excess
6.Least square analysis Method
Use any that is fast convenient and sufficiently
accurate
1.Integral Method
Involves the following steps:
1. Assume order(s) of reaction and hence rate
law
2. Integrate the differential form of the rate law
3. Make the appropriate linear plot based on the
integrated equation
4. If the linear plot is consistent with the
expected straight line, then the assumed
order is correct
5. Otherwise the process is repeated with
another order of reaction
6. Hence important to know linear plots for 0, 1,
2 and 3 orders of reaction
2A
products
products
Differential Method Of
Analysis
dC
dt
k f (C i )
Differential Method:
obtaining order ln(dCA/dt)= lnk + ln CA
*
ln (dCA/dt)
intercept=
lnk
slope=
ln CA
Obtaining (-dCA/dt):
Any mathematically convenient and
accurate method of obtaining the
slope
measuring the slope of the tangent to the
CA vs t. plot at various points as shown in
fit a polynomial of sufficient accuracy to
the CA vs. t plot :
CA(t)=a0 + a1t + a2t2 + a3t3+ a4t4
Microsoft Excel can be easily used for this.
Differentiating equation with respect to t
yields
dCA(t)/dt= a1 + 2a2t + 3a3t2+ 4a4t3
Differential method of
analysis
kC A of
C B time
kC A C B0 kC A
a- function
dt
kCA CB kCA0 CB k CB
dt
Method of half-life:
Half-life period or the half-life
ln (CA0/ CA )
=kt
At half-life, t= t1/2, CA=0.5CA0
ln (CA0/ 0.5CA0 ) =ln2= kt1/2
Half-life Is Independent Of
Concentration
First Order Reaction
1
1
kt
CA
C A0
kC A0
C A0
(0.5 1)
t1/ 2
(n 1)k
(0.5 1)
(1 n) ln C A0 ln
ln t1/ 2
(n 1)k
n=1-slope
(0.51 n 1)
ln
= (n 1)k
Intercept
depends on the concentration of the
reactant
and the order of the reaction.
Taking logs of both sides
ln(t1/2
)
(0.51 n 1)
int ercept
(n 1)k
Slope =1n
ln(CA
0)
lnt1/2
Orders of
reaction < 1
Orders of reaction
>1
Orders of reaction
=1
ln CA0
*
ln(-rA)
ln(CA0)
Method of excess.3
k' is a pseudo--order rate constant.
The reaction is a pseudo- -order
reaction, because it behaves as if it
were -order although it is not
really -order.
It is -order with respect to only A.
k' is not a true constant
its value varies with the initial
concentration of B, CB0.
*
By performing a series of experiments in
which CB0 varies but is always much
larger than CA0, one obtains k' as a
function of CB0
*
ln(k')
ln(CB0)
r
=
-dC
/dt
=k
C
C
-kC
A
A
A
B
C
CDb
At t=0, and with no products in
the initial charge, CC0=0;
CD0=0;
ln( rA0)
slope=
ln(k. CB0)
ln(CA0)
ln( rA0)
slope=
ln(k. CA0)
Leastsquare method
The leastsquare method can be
used to obtain the best values of the
parameters of the rate law-the
orders of reaction and the rate
constant.
For a rate law of the form:
linearizing the functional form of the
law:
and C are the variables measured
Leastsquare method2
vs. ln CA
the plot of
represents the equation of a straight
line of the form y=mx+c,
Where y=
; x= ln CA
Run
CAi (-rAj)
Xi
=lnCA
Yi
=ln(-rA)
xiyi
xi2
CA1 (-rA1)
lnCA1
ln(-rA1)
lnCA1*ln(-rA1) (lnCA1)2
CA2 (-rA2)
lnCA2
ln(-rA2)
lnCA2*ln(-rA2) (lnCA2)2
CA3 (-rA1)
lnCA3
ln(-rA3)
lnCA3*ln(-rA3) (lnCA3)2
CA4 (-rA4)
lnCA4
ln(-rA4)
lnCA4*ln(-rA4) (lnCA4)2
Can (-rAn)
lnCAn
ln(-rAn)
lnCAn*ln(-rAn) (lnCAn)2
lnCAi
ln(-rAi)
summations
Least square
The two simultaneous
equations
method
are solved for the unknowns.
These are the best values of the
parameters obtainable from the
experimental data
Linearizing
yi =
c + ax1i + ax2i
Two independent variables:
x1=lnCA and x2=lnCB
The summation of the results of all
the n runs
x2i
yi
x1i2
yi x2i
XA =(V- V0 )/ AV0
dXA/dt= [1/AV0]dV/dt
CA= NA/V =NA0(1 XA)/V0 (1 + AXA)
CA= CA0(1 -XA)/(1 + AXA)
CA
= (-1/ V0 (1 + AXA)(-NA0dXA)/dt)
A > 0
ln(V/V0)
t
A < 0
-ln[1
V/AV0]
Q & A
Nitric acid is made commercially from nitric oxide. Nitric
oxide is produced by the gas-phase oxidation of
ammonia.
4NO + 6H2O
4NH3 + 5O2
The feed consists of 15 mol% ammonia in air at 8.2 atm
and 227C.
What is the total entering concentration?
What is the entering concentration of ammonia?
Set up a stoichiometric table with ammonia as your
basis of calculation.
Express Ci for all species as functions of conversion for a
constant-volume batch reactor. Express PT as a function of X.
Express Pi and Ci for all species as functions of conversion for a
flow reactor.
initial charge
Ni0
change due to
reaction at NA
NA0
-( NA0-NA)
NA0-( NA0-NA)=NA
= NA0-(NT-N0) (a/n)
NA0B
-( NA0-NA) (b/a)
NA0B -( NA0-NA)(b/a)
=NA0B -(NT-N0)(b/n)
NA0C
( NA0-NA) (c/a)
NA0D
( NA0-NA) (d/a)
NA0I
NA0I
A
B
C
D
NA0 -(NT-N0)a/n
NA0B -( NT-N0)(b/n)
NA0C +( NT-N0)(c/n)
NA0D +( NT-N0)(d/n)
NA0I
NA=NA0-(NT-N0)a/n
PAV/RT=PA0V/RT- (PT-P0)Va/nRT
PA=PA0- (PT-P0)a/n
NB=NA0B -( NT-N0)b/n
PBV/RT=PA0BV/RT- (PT-P0)Vb/nRT
PB=PA0B- (PT-P0)b/n
NC =NA0C +( NT-N0)(c/n
PCV/RT=PA0CV/RT+ (PT-P0)Vc/nRT
PC=PA0C+ (PT-P0)c/n
ND =NA0D +( NT-N0)(d/n
PDV/RT=PA0DV/RT+ (PT-P0)Vd/nRT
PD=PA0D+ (PT-P0)d/n
NI =NA0I
PIV/RT=PA0IV/RT
PI=PA0I
Question.
The gas phase decomposition of ethylamine produces ethene and ammonia,
C2H5NH2(g)-- C2H4(g)+ NH3(g)
The following total pressures were found for the reaction at 500 0C with an initial pressure of pure ethylamine
equal to 55 mmHg.
Ptotal,mmHg: 55
:
Time,min 0 2
64
4
72 89 93
10 12
Using this and the following data, find the partial pressure of C2H5NH2 remaining at the
various times, and plot an appropriate graph showing the progress of reaction.
What will the final pressure be?
Show that the decomposition is not first order
What is the rate law of the decomposition
REVERSIBLE REACTIONS
In irreversible reactions the reactants
react until they are completely consumed
and products are formed.
In reversible reactions the products
formed do react to give back the
reactants.
Reversible reactions are thus
represented as shown, with a double
arrow:
A
D. forward reaction
The reactants, A , react to form the
product D. The reverse reaction is
represented as
A
D. reverse reaction
The products of the forward reaction D
react to form the reactants
forward and reverse reactions take
place simultaneously
Dynamic Equilibrium:
Equilibrium Constant
The larger the value for the equilibrium
constant the more the reaction goes to
completion.
The value for K is large when products
dominate the mixture.
The value for K is small when the
reactants dominate the mixture.
Irreversible reactions can be thought to
have an infinite equilibrium constant so
there are no reactants left
30
10
9
25
k1 =3/min; k1 '=1.0
k1 =3/min;
k1 '=1.0/min
20
rA =reverse reaction
CA
6
5
mole/L
4
mole/L.min 15
10
2
1
0
0
0 0.5 1 1.5 2
t,min
0.5
1.5
t, min
2.5
Ni0
Ni
Ci
NA0
-(NA0- NA) NA
CA
(NA0- NA)
CA0- CA
NA0- NA
Integrating leads to
At equilibrium,
and .
Simplifying
Hence a plot of
vs.
should give a
straight line through the origin with slope
equal to k1 .
slope=k
+k
slope=k
+k
k1
k2
k1
k2
10
C
A
10
C
C
CB
8kn;=1.0/mi
1
k1=1.0/mi
n;
k2=2.0/mi
n
k2=10.0/
min
6
mole/L
mole/L
1
0
0
1
t,min
t,min