CRE10 Non Isothermal NonAdiabatic Reactors
CRE10 Non Isothermal NonAdiabatic Reactors
CRE10 Non Isothermal NonAdiabatic Reactors
Topic 4
Non-Isothermal, Non-Adiabatic
Reactor Design and Multiple
Steady States
(Chapter 12- Foglers Essentials)
1
cCpc (Ta1 Ta 2 )
Q m
UA(Ta1 Ta 2 )
ln[ (T Ta1 ) /(T Ta 2 )]
8.46
UA
c Cpc (Ta1 T ) 1 exp
Qm
c Cpc )]
(m
8.49
For high coolant flow rates , Ta1 Ta2=Ta and the 2nd order terms
are neglected in the Taylor expansion of exp function. Therefore,
Q UA(Ta T )
8.50
Q U (Ta T ) dA
U a (Ta T ) dV
dQ
Ta
T
dE
T i Cpi dT ] FA0 X H Rxn (T ) dt
0
Now, that we have seen how Q can be described, let us reexamine the general energy balance equation
At Steady State, the energy balance equation for
a single rxn system is
8.48
Cp
C p
dT
Tref
T Tref
8.28
~
Cpi
Cp
dT
T0
8.29
T T0
C
p
(
T
T
)
0
i i
i 1
X
H Rxn (Tref ) C p(T Tref )
8.49
XEB
Q 0
W 0
n
UA(Ta T )
~
H Rxn (Tref ) C p(T Tref ) X i Cpi (T T0 )
FA0
i 1
F
i 1
i0
Hi0
FH
i 1
FA0
T0
dQ
Ua (T T )
a
dV
Ta
FAe
T
Te
V V+V
0
H i Fi
0
dV dV
dV
i 1 dV
i 1
10
Let us evaluate each term to obtain an equation in terms of dT/dV
dFi
(ri )
dV
i (rA )
H i (T ) H io (Tref )
dH i
dT
Cpi
dV
dV
C pidT
T Tref
T
i 1
i 1
dT
0
dV
Rearranging,
n
dT
dV
U a (Ta T ) i (rA ) H i (T )
i 1
F Cp (T )
i 1
11
dT
dV
U a (Ta T ) i (rA ) H i (T )
i 1
F Cp (T )
i 1
dFi
i (rA )
dV
or
g( X ,T )
dX
FA0
(rA )
dV
f ( X ,T )
dT
dV
U a (Ta T ) i (rA ) H i (T )
i 1
F Cp (T )
i 1
n
dV
Fi Cpi (T )
i 1
n
dV
FA0 [ i Cpi (T ) X Cp]
i 1
1 Ae E / RT
Q Ws FA0 i C pi (T Ti 0 ) FA0 X H RX
(TR ) C p (T TR ) 0
i 1
n
UA(Ta T )
~
X H RX (TR ) C p (T TR ) i C pi (T Ti 0 )
FA0
i 1
~
iC pi (T Ti 0 ) UA(T Ta ) / FA0
i 1
H RX
(TR ) C p (T TR )
15
Ae E / RT
X
1 Ae E / RT
H RX
(TR ) H C (68 F ) H B (68 F ) H A (68 F )
pC
pB
pA
~
~
~
~
C
i pi A pA B pB M pM
i 1
FB 0 ~
FM 0 ~
~
C pA
C pB
C pM
FA0
FA0
403.3Btu / lb mol propylene oxide ( A) F
C
i pi (T Ti 0 ) UA(T Ta ) / FA0
i 1
H RX
(TR ) C p (T TR )
V
V
0.1299h
v0 v A0 vB 0 vC 0
Ti 0 75 F
TR 68 F
Ta 85 F
UA 100 40
F
i 1
i0
Hi0
FH
i 1
FA0
T0
dQ
Ua (T T )
a
dV
Ta
FAe
T
Te
V V+V
0
H i Fi
0
dV dV
dV
i 1 dV
i 1
17
Let us evaluate each term to obtain an equation in terms of dT/dV
dFi
(ri )
dV
i (rA )
H i (T ) H io (Tref )
dH i
dT
Cpi
dV
dV
C pidT
T Tref
T
i 1
i 1
dT
0
dV
Rearranging,
n
dT
dV
U a (Ta T ) i (rA ) H i (T )
i 1
F Cp (T )
i 1
18
dT
dV
U a (Ta T ) i (rA ) H i (T )
i 1
F Cp (T )
i 1
dFi
i (rA )
dV
or
g( X ,T )
dX
FA0
(rA )
dV
f ( X ,T )
dT
dV
U a (Ta T ) i (rA ) H i (T )
i 1
F Cp (T )
i 1
n
dV
Fi Cpi (T )
i 1
n
dV
FA0 [ i Cpi (T ) X Cp]
i 1
d /a
De
b/a
Be
C C
KC [
C C
Equilibrium
]
X
~
i C pi (T T0 )
H Rxn (T )
T0
21
Equilibrium
Increasing the inlet
temperature results in
shifting of the EB
equation to the
RIGHT
Xe
Adiabatic Temperature
T01
T02
22
23
Xe
T01
T02
i 1
Ti 0
X EB
C
i pi (T T0 )
i 1
H RX (T )
24
AB
Make a plot of equilibrium conversion as a function of time.
Determine the adiabatic equilibrium temperature and conversion
when pure A is fed to the reactor at a temperature of 300 K.
25
AB
Rate law : rA k C A B
K e (T )
equilibrium
-rA = 0
C A0 X e
C A0 (1 X e )
C Ae
CBe
Ke
K e (T )
Xe
1 K e (T )
d (ln K e ) H
H RX (T ) H RX
(TR ) C p dT
TR
C p C pB C p A 0
H 1 1 H RX H B H A 20000
K e K e (T1 ) exp
R T1 T T1 298K
Only thermodynamics!
Xe = f (T) nothing to do with the energy
balance
RX
26
X EB
~
iC pi (T T0 )
i 1
C
i 1
1 C p A
pi
H RX (T )
X EB
X EB
XEB
From thermodynamics
CPA (T T0 )
H RX (T )
50(T 300)
20000
27
T0 = 600
T0 = 350
0.15
350
500
600
28
29
XEB
final conversion
cooling process
(with interstage cooling)
30
X e 0.41T0 300K
Temperature leaving
the first reactor
X
T1 460K
(T 300)
EB
400
cooling to 350 K
Cooling duty ?
95 % of equilibrium conversion
X 0.4
31
n
Q W s Fi 0 H i 0 Fi H i 0
i 1
i 1
Q Fi H i Fi 0 H i 0
i 1
i 1
Fi Fi 0
F
i 1
i0
(Hi Hi0 )
Q Fi 0C pi (T2 T1 )
i 1
Q
220000
m C
94 mol / s
Q m C C pC (Tout Tin )
being absorbed by
Q FA0C p A (T2 T1 )
the coolant stream
40 50 (350 460)
220 kcal / s
32
Th2= 460 K
Th1= 350 K
Reaction mixture
Heat Exchanger
Tc2= 400 K
Tc1= 270 K
Coolant
A = 31.6 m2
X EB
(T 300)
400
energy balance
X 0.6
95 % of equilibrium
conversion
X EB 0.63
33
Endothermic reactions
XEB
final conversion
heating process
(with interstage heating)
Xe
Adiabatic Temperature
T01
T02
34
Xe
X3
X EB
X2
X1
~
i C pi (T T0 )
H Rxn (T )
35
Non-adiabatic reactor
operation :
36
S O2 SO2
1
2O5
SO2 O2 V
SO3
2
2 k
rSO
PSO2
PSO3
PO
2
PSO3
K p PSO
2
The feed to an SO2 converter is 7900 lb mol/h and consists of 11% SO2, 10%
O2, and 79% inerts. The converter consists of 4631 tubes packed with
catalyst,
each 20 ft long.
The tubes are 3 in. o.d. and 2.782 in i.d. The tubes will be cooled by a boiling
liquid at 805 F, so the coolant temperature is constant over this value. The
entering pressure is 2 atm.
For inlet temperatures of 740 and 940 F, plot the conversion, temperature,
equilbrium conversion, and reaction rate profile down the reactor.
37
42311
K p exp
11.24
RT
T in R, Kp in atm-1/2
176008
k exp
(110.1ln T ) 912.8
T
T in R
1
A B C
2
1
2O5
SO2 O2 V
SO3
2
Mole balance
(PFR) :
dX
FA0
rA
dW
dX
k
dW FA0
Rate law :
PSO2
PO
2 k
rSO
PSO3 2
38
Stoichiometry :
(i i X )( RT ) P
Pi Ci RT C A0
(1 X )(T / T0 ) P0
2
PSO3
K p PSO
2
PA0
( i i X ) P
(1 X ) P0
X
O2
1 X P
2 SO3 X 1
P
A0
1 X K 2
SO3 X P0
1 X
per tube
dX
1 X P 0.2 0.11X
X
dX f ( X , T , P)
5.32k
dW
X P0 1 0.055 X (1 X ) K p dW
39
n
dW
FA0 ( i C pi XC p )
i 1
Ua(Ta T ) (rA ) H RX
(TR ) T T 2 dT
TR
dT
n
dW
FA0 ( i C pi XC p )
i 1
dT 5.11(Ta T ) (rA ) H RX (T )
n
dW
0.188( i C pi XC p )
dT
f ( X , T , P)
dW
i 1
40
Ergun equation:
dP
T P0 FT
dW
2 T0 P / P0 FT 0
2 0
AC C (1 ) P0
G
0
0 gc Dp
1 150(1 )
3
1.75G
Dp
GTP0 (1 X )
dP
dW
b Ac 0T0 PD p g c
1 150(1 )
3
1.75G
Dp
Viscosity is also a function of
temperature,
but we ignore here.
dP
1.12 108 (1 0.55 X )T
5500 2288
dW
P
dP
f ( X , T , P)
dW
41
Mole balance
dX
1 X P 0.2 0.11X
X
5.32k
dW
X P0 1 0.055 X (1 X ) K p
Energy balance
dT 5.11(Ta T ) (rA ) H RX (T )
n
dW
0.188( i C pi XC p )
i 1
dP
1.12 108 (1 0.55 X )T
5500 2288
Momentum balance
dW
P
T
T
T
TR
i 1 i 0
Set Cp = 0
UA(T Ta )
~
FA0 i C pi (T T0 )
FA0 XH RX
(T )
FA0
i 1
n
= Cp0
CSTR
FA0 X rAV
42
i 1
T0
W
F
Q
s
A0
dT
dt
43
and UA
FA0C p 0
and Tc T0 FA0 C p 0 U A Ta
U A FA0 C p 0
as
C p 0 (1 ) (T Tc ) (rA V ) H Rxn (T ) Ni Cpi
i 1
dT
dt
i.e.
44
T0
R(T)
R(T ) C p 0 (1 ) (T Tc )
Tc
T0 FA0 C p 0 U A Ta
U A FA0 C p 0
R(T)
Ta
T0
T
UA
FA0C p 0
45
G(T)
Low Ea
High Ea
46
G(T)
&
R(T)
T
47
48
UA
C p 0 FA0
Tc
T0 FA0C p 0 UATa
UA C p 0 FA0
R(T)
Ta T0
Increase T0
T
C p 0 (1 )(T Tc )
rAV
(H RX
)
FA0
R(T)
=0
Increase
For Ta < T0
Ta
T0
G(T ) (H RX
)X
function of T ...
i.e. mole balance, energy balance .
49
ignition temperature
extinction temperature
50