CHEMICAL REACTOR DESIGN

FOR PROCESS PLANTS

Volume Two:
Case Studies and Design Data
HOWARD F. RASE
W. A. Cunningham Professor of Chemical Engineering
The University of Texas at Austin
Original Illustrations by
JAMES R. HOLMES
Associate Professor of Engineering Graphics
The University of Texas at Austin
A WILEY-INTERSCIENCE PUBLICATION
JOHN WILEY & SONS, New York London Sydney Toronto
CASE STUDY 101
Styrene Polymerization
THIS STUDY is a nontrivial case of a CSTR involving the complexities of
free-radical polymerization. The weakest aspect of the model is the method
for estimating viscosity, which strongly affects heat transfer and is very
sensitive to temperature and concentration of polymer. Experimental
work with the particular polymer-styrene product mix would yield a better
correlation.
Problem Statement
Design the first stage of a two-stage thermal polymerization system for
producing 40 miIIion lb/yr of polystyrene with an overaII conversion for the
system of 95 ~ { , . The product from the first stage should have a number
average molecular weight (M
n
) of 144,000.
Feed. 99.5 wt %styrene with 10 ppm maximum polymer and maximum
weight percentages of impurities of 0.02 aldehydes as CHO, 0.01 peroxide as
H
2
0
2
, 0.0025 sulfur, and 0.01 chloride.
Chemistry and Kinetics
Details of the formation of a Diels-Alder adduct (AH) and its radical (A' )
from styrene monomer (M) have been described (l) and are partiaIIy iIIus-
trated as foIIows, together with the propagation, termination, and chain-
transfer steps.
5
6 CS-101 STYRENE POLYMERIZATION
Initiation
M+M
M +AH
M +AH
rn
~ (caIlAH)
H H Ph
M'+A'
trimer
(1)
(2)
(3)
A' + M
k
R'
~
1
M'+M
~
R'
1
Propagation
R
r
' + M
~
R
r
+ l'
Termination by Combination
R
r
' + R,'
k,
P"+s
-------->
Chain Transfer. (to AH or M)
R AH
k" P R
/.' + --------> r + l'
The rate of initiation is
(4)
(5)
(6)
(7)
(8)*
(CS-l.l )
(CS-1.2)
By using the stationary state hypothesis, expressions for concentrations
of A', M', AH, and R
r
' can be obtained as functions of monomer concentra-
tions so that Eg. CS-l.l can also be written in these terms (1)
h) = 2k1 kz[M]3 "'-
L 1 + (k
z
+ k
3
)[M] + ktr(rl/k
t
)'
where [M] is the concentration of monomer in g moles/liter. Alimiting case in
which k _1 is much greater than the remainder of the denominator appears
to fit operating data.
(CS-1.3)
AH represents a Diels-Alder adduct.
CHEMISTRY AND KINETICS
7
Other instantaneous expressions typical of free-radical polymerization with
termination by combination dominating are (l).
From Eq. CS-1.3
r
pr
= (k ~ : i 2)+[M]t = A[M]t
. t. pr
The grouping A has been found to vary with conversion.
A = A
o
exp (AtX + A
2
X
2
+ A
3
X
3
)
Also
or
- k.[M]+
(DP)-l = C + -'--
m A
(CS-1.4)
(CS-1.5)
(CS-1.6)
(CS-l.7)
(CS-1.8)
(CS-1.9)
where M
n
is the number average molecular weight. Cm, which accounts
for chain transfer, also varies with conversion.
C
m
= C
mu
+ B1X
Equations for the several constants are given in Table CS-l.l.
Table CS-l.l Kinetic Data for Thermal Polymeriza-
tion of Styrene (Ref. 1)
A
o
= 1.964 X 10
5
exp( -1O,040/Tl.(I./g mole)1 sec-
1
Al = 2.57 - 5.05 x 1O-
3
T
A
z
= 9.56 - 1.76 x lO-
z
T
A
3
= -3.03 + 7.85 x 1O-
3
T
-3 (473.12 - T)
B1 = -1.013 x 10 logl0\ 202.5
k
j
= 2.19 X 10
5
exp (-13,81O/T), (l./g-mole)Z sec-
1
e
mo
= 2.198 x 10-
1
exp (- 2820/Tl
T = OK
(CS-l.10)
8
Thermodynamics
CS-101 STYRENE POLYMERIZATION
The typical high exothermicity of polymerization, I1H = - 17,500 caljg mole
monomer converted @ 25C for this case (2,3), and the importance of
temperature in controlling molecular weight indicates a definite need for a
reactor with heat transfer. FromTable 6.1, pp. 270-271
1
, the adiabatic factor is
moderately high, and the heat generation potential is not high.
Reactor Type
A CSTR has good heat-transfer characteristics and allows continuous
operation, but it can only be used to a conversion level where heat transfer and
contacting become poor because of increasing viscosity. This problem can be
overcome by diluting with solvent, but then solvent handling and recovery
add greatly to the cost. Alternatively, a CSTR can be specified for some low
conversion range, and the remaining conversion (up to 95) can be completed
in a screw-type, tubular flow unit, such as described on p. 480
1
, or as shown
in Fig. 10.29. We will only consider the CSTR, for the tube operates at such
high viscosity that diffusion controls and modeling techniques are inadequate.
Styrene polymerizes slowly relative to many other monomers and a CSTR
is also advantageous for this reason.
Design Model (CSTR)
Equations CS-1.5-CS-1.l 0 express instantaneous values applicable to a
CSTR. The following mole and heat balances apply.
Monomer Balance. (neglecting monomer used in initiation and transfer)
Heat Balance
FMX
M
= -d[M] V = r V = A [ M ] ~ V
dt pr
(CS-l.ll)
104F
M
c
p
JI;, - To) + UAh(T. - 1j) = (-I1H
M
hJF
M
X
M
) (CS-1.l2)
where FM is the molar flow rate of monomer fed, CPM is the heat capacity of
monomer, and T., To, 1j are the exit, entrance, and cooling medium tempera-
tures, respectively.
Since the heat capacities of monomer and polymer are approximately
the same, the heat of reaction at 25C can be used in Eq. CS-1.l2.
OPERATING TEMPERATURE AND CONVERSION
Fluid Properties
C
PM
= 0.353 + (0.0014)(T - 293), caljgOC
PM = 924 - 0.918( T - 273.1), gil.
Ppr = 1084.8 - 0.605(T - 273.1), gil.
6,. = e
l
pr - :JIP ~
PM(l - X
M
)
[M] = (l 04)(l + 6,. XM)' g moles/l.
('1r - 1)(W
rr
= 109['1Jr exp {2079w
pr
(l09['1Jr + 1)/T}
9
(CS-1.13) (2,3)
(CS-1.14) (1)
(CS-1.15) (1)
(CS-1.16)
(CS-1.17)
(CS-1.18) (4)
where PM and Ppr are the densities of monomer and polymer, respectively.
'1r is the ratio of viscosity of polystyrene in styrene to that of pure styrene at
the same temperature, ['1Jr is the intrinsic viscosity of polystyrene of same
molecular weight and type in toluene at 30C in dl./g, and w
pr
is the weight
fraction as mass of polymer/mass of styrene.
The viscosity is very sensitive to temperature and concentration. Equation
CS-l.18 is corrected to agree with the original plot of experimental data given
in the reference. Ideally the value of ['1Jr should be determined for the
product in question in the laboratory. For illustrative purposes we select a
value of 0.7.
Design Calculations and Decisions
Operating Temperature and Conversion
Molecular weight in thermal polymerization is most strongly affected by
operating temperature and is not sensitive to conversion. By applying
Eqs. CS-1.6-CS-1.9 at various temperatures using a hand calculator, a
temperature of 150
a
C gave M
n
= 144,041 @ X = 0.4 and M
n
= 143,621 @
l
x = 0.45. The corresponding solution viscosities at these two conversions
are 3623 cp and 13406 cp, respectively. Styrene viscosity at 150C is 0.22 cp (5).
Thus it is seen that viscosity increases dramatically above X M = 0.4, and
heat-transfer rate will decline rapidly such that it will not be possible to
remove the heat of polymerization. Hence for safety reasons the design will
be set at X M = 0.4 as a maximum, and subsequent calculation will reveal that
this is operable.
10
Monomer Flow Rate
CS-101 STYRENE POLYMERIZATION
Basis: 95 %overall conversion for 2-stages.
operating factor
99.5 %styrene purity and no recycle
4 x 10
7
lb/yr
(0.995)(0.95)(0.95)(365)(24)(3600) = 1.413 Ib/sec or 641 g/sec
Reactor Volume
From Eq. CS-l.II
V = FMXM = (641/104)(0.4)
A[M]t (1.64243 x 10- 5)(4.9506)t
= 2752.7 liters or 727 gal
A 750 gal stainless-steel tank (60 x 60 in.) is selected from Table 8.8,
D = 60 in., D
1
= (1/3)60 = 20 in. with 160 rpm for good heat transfer
(838 ft/min tip speed, see Table 8.10). One 20-in. flat-blade turbine is selected
initially for analysis.
Heat Balance (Eq. CS-I.12)
Basis: 150C which is the reaction temperature
30C inlet, 45C average cooling water temperature.
(641)(0.451)(150 - 30) + UA
h
(l50 - 45) = (17,500)(641/104)(0.4)
Required UA
h
= (43144 - 34691)/105 = 80.5 cal/sec DC,
_ (80.5)(3600) _ 6 f 2
U - - .8 BTU/hr t
Estimate hand U. From Eq. 8.14
h = 0.73),! (C
p
fl)0.33(PND
I
2
)0.65(!!-)0.24
D),! fl flw ,
= 0.73 (0055)[(0.451 )(3623)(2.42)JO.33
20/12 . 0.055
[
(53.5)( I60)(60)(20/1 2)2JO.65( fl )0.24 _ (fl )0.24
X - 26.4
(3623)(2.42) flw flw
_ PM/1000 0.786 2 3
P = I + 6"X
M
= I _ (0.209)(0.4) = 0.858 g/cm or 53.5Ib/ft
A
r
= 0.055 BTU/(hr)ft
2
OF ft-
1
@ 150
a
C for ethylbenzene (API Data Book)
REQUIRED HORSEPOWER
11
or
Estimation of wall temperature (b
w
= t6 in., fouling factor, = 0.001,
h
o
is for water = 400 BTU/hr ft2 OF.
Assuming h = 10
T
w
- 45 = 0.044
150 - T
w
for which J1 from Eq. CS-1.l8 is 52,477 cp
corrected h = y.24 = 13.9
I . . i
V =\,9.98 + 4.367 x 10- 3)= 13.1 BTU/hr ftl OF
Thus the design is adequate compared with required V of 6.8. It should be
noted that an increase in conversion to 45 at a lower flow rate to yield the
same production would produce a h equal to 6.95 or a V of 6. 75, based on the
viscosity correlation of Eq. CS-l.18. The required V at this new condition
would be 9.9. Thus somewhere between 40 %and 45 %conversion the reac-
tion will become uncontrollable by conventional jacket cooling, but the
reaction is slow and corrections can be made. The vapor pressure of styrene
at 150
a
C is 880 mm Hg (2), but at the conversion level the actual value will be
less because of nonidealities attributable to the dissolved polymer. An
approximate operating pressure of 2.5 psi will be assumed. This pressure
or higher values at higher temperatures which could occur during upsets
can be released to cause vaporizing of the monomer and rapid cooling.
Required Horsepower
Based on manufacturer's correlation Fig. 8.8 at NRe = 163, Np = 5.5
hp = (1.1)(3.52 x + 0.5 = 4.95 hp or
specify 5 hp motor.
If one refers to Fig. 8.7, the Reynolds number is in the transition region for
which the power number is 3.5 and the corresponding hp = 3.3; but the
12 CS-101 STYRENE POLYMERIZATION
higher hp would be selected since the reactor must be periodically washed
with a low viscosity solvent which would place the operation in the turbulent
region with a higher power number.
Mixing Effectiveness
Equation 8.16 can be used to estimate the effective radius of agitation.
fhP )4.05
Rerr = c ~ --;; = c 3623= 0.033c
Using appropriate values of the coefficient, horizontal and vertical values
are 4.2 and 1.7 ft, respectively, which is adequate in the horizontal but not
the vertical direction for a 60 x 60 in. tank. Accordingly, two 18 in. impellers
should be substituted. By selecting the 18 in. impellers it is possible to main-
tain the required power at the same value, (18/20)5(0.9)(2) ;:;:: 1.0 (see instruc-
tions in caption for Fig. 8.7, p. 348
1
). The turbines should be spaced 30 in.
apart.
A better decision would be to use a single curved-blade turbine (with six
blades) since it provides superior axial mixing. The resulting hp will be lower,
but the heat-transfer coefficient will be approximately the same (see Table 8.6).
Because of the flatter power-number curve (Fig. 8.7), the power overage
required for washing will not be as great. Manufacturer's representatives
should be consulted for the final selection.
REFERENCES
I. A. W. Hui and A. E. Hamielec, J. Appl. Polym. Sci., 16, 749 (1972).
2. K. E. Coulter, H. Kehde, and B. F. Hiscock, in High Polymers, Vol. 24, Part 2, E. C.
Leonard (ed.), Wiley-Interscience, New York, 1971, p. 501.
3. R. H. Boundy and R. F. Boyer, Styrene, ACS Monograph, Reinhold, New York, 1952.
4. M. Hirose, E. O'Shima, and H. Inoue, J. Appl. Polym. Sci., 12,9 (1968).
5. K. E. Coulter, H. Kehde, and B. F. Hiscock, in High Polymers, Vol. 24, Part. 2, E. L.
Leonard (ed), Wiley-Interscience, New York, 1971, p. 499.