Fluidzbed Reactor
Fluidzbed Reactor
CATALYSTBEDS
Terukatsu Miyauchi and Shintaro Furusaki
Shigeharu Morooka
and
Yoneichi lkeda
I. Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 276
A. Purpose and Outline of the Review . . . . . . . . . . . . . . . . 276
B. General Properties of Fluidized Beds . . . . . . . . . . . . . . . 277
C. Historical Development of FCB Studies . . . . . . . . . . . . . . 281
D. Quality of Fluidization in Relation to Particle Properties . . . . . . . 283
11. Flow Properties of Fluid Beds . . . . . . . . . . . . . . . . . . . . . 285
A. Particle Properties in Relation to Fluidity . . . . . . . . . . . . . 285
B. Formation, Splitting, and Coalescence of Bubbles . . . . . . . . . . 290
C. Operation of Fluid Beds . . . . . . . . . . . . . . . . . . . . . 297
D. Properties of Bulk Flow Inside Fluid Beds . . . . . . . . . . . . . 297
E. Axial Distribution to Bed Structure . . . . . . . . . . . . . . . . 305
F. Effect of Bed Internals . . . . . . . . . . . . . . . . . . . . . . 307
111. Turbulent-Flow Phenomena in Bubble Columns and Fluidized
Catalyst Beds . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 10
A. Two-Phase Bubble Flow in Recirculation . . . . . . . . . . . . . . 311
B. Comparison of Theory with Experiments on Bubble
Columns . . . . . . . . . . . . . . . . . . . . . . . . . . . . 317
275
Copyright @ 1981 by Academic Press, Inc.
ADVANCES IN CHEMICAL ENGINEERING, VOL. 11 All rights of reproduction in any form reserved.
ISBN 0-12-008511-9
276 M I Y A U C H I . F U R U S A K I . MOROOKA. A N D I K E D A
1. Introduction
A . PURPOSE
A N D OUTLINE
OF THE REVIEW
Fluidization today seems a very young and active field. even though
Charles E . Robinson originated the technique a century ago . It is being
applied today in the chemical industry to catalytic reactions. noncatalytic
F L U I D I Z E D CAT AL YST BEDS 277
B. GENERAL
PROPERTIES
OF FLUIDIZED
BEDS
In fluidization, a suspension of fine solid particles behaves like a liquid
during the upflow of a supportive gas or liquid phase. Thus the bed of
fluidized solid itself may be analyzed similarly to liquid systems. The
gas-lift effect produces internal recirculation, by providing a descending
flow of high particle concentration and an ascending flow of low particle
concentration. This effect resembles the circulation in bubble columns.
Whereas bubble columns contain dispersed gas and a continuous liquid
phase, the fluidized bed comprises the bubble phase and the emulsion
phase in which particles have gained fluid-like properties by the interstitial
gas flow.
The interrelation of pressure drop and superficial gas velocity appears
to determine the various states of fluidization, as shown in Fig. 1 by
278 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
Yerushalmi et al. (Y 11, Y 12). With increasing gas velocity, the pressure
drop increases while the bed is in a fixed state; then it stays almost
constant, equal to the fluidizing density of the bed after the minimum
fluidization velocity. Further increase of gas velocity makes the bed
fluidize vigorously; a change occurs from the bubbling bed to the turbu-
lent bed and finally to the fast fluidized bed, and the amount of particles
entrained by the gas flow likewise increases. The charge of particles in a
fluidized bed will not remain constant unless the particles are separated in
cyclones and returned to the bed.
Most industrial fluidized beds using fine powder are found to be in the
category of turbulent beds; typical of these is fluidized catalytic cracking.
On the other hand, most noncatalytic reactions and physical operations
use coarse particles and are operated as bubbling beds. Squires (S14) has
pointed out that the state of fluidization varies significantly with fluidized
particle diameter, and has identified the following two categories:
FLUIDIZED CATALYST BEDS 279
(a) Fluid bed. The fluidization of fine particles, generally all through 20
mesh; usually with a substantial part smaller than 200 mesh, often smaller
than 325 mesh. Superficial velocities generally below 60 c d s e c .
(b) Teerer bed. The fluidization of coarse particles, generally all larger
than 60 mesh; with a substantial part larger than 10 mesh. Superficial
velocities generally above 150 c d s e c .
Considering Squires’ criteria and observed industrial performance
(G15, K24, M2,010, V12,Z5), Ikeda (14) modified the criteriaas shown in
Table I, stressing that the fluidized catalyst bed should be viewed as the
typical fluid bed (11-110). The general flow features of a fluid bed operated
under relatively high gas velocity are illustrated in Fig. 2 (M25). The
dense bed consists of the emulsion phase and bubble phase; the emul-
sion phase exhibits liquid-like properties, in the case of the fluid bed.
Intense circulation of dense phase results from the central rapidly ascend-
ing bubble-rich phase and the peripheral descending emulsion phase of
low bubble content. The boundary of the dense phase and the dilute phase
is rather ambiguous; a transition region exists in which bubbles collapse
and particles are blown in from the emulsion phase. In the transition
TABLE I
TYPICAL
FLUIDBEDAND TEETERBED
Constitution Bed divides into dense phase + Density of bed is almost constant,
dilute phase, although surface of bed rather well
surface of dense phase is defined
ill-defined
Behavior of Most gas flows as uniformly small Large bubbles, which may
bubbles bubbles coalesce progressively
during rise
Behavior of Circulative flow is developed Generally relative movement of
emulsion (possibly central upward, particles; partially circula-
phase peripheral downward) tive flow
Stability Slugging unlikely, even for Frequent slugging for
large LJD, L,, > 1-1.5 m
280 M I Y A U C H I , F U R U S A K I , MOROOKA, A N D I K E D A
region, inversion of the circulating flow can occur. In the dilute phase,
particles are suspended particulately in the stream of ascending gas. The
content of suspended particles decreases with height, except in the upper
part of the dilute phase where it stays constant.
TABLE I1
SCALE-UP &TI0 AND REACTION
IN T H E DEVELOPMENT
OF FCC
Reactor
Bed diameter, DT 2 in. 15 in. 15 ft.
Relative diameter 0.133 1 12
Relative cross section 0.018 1 144
Bed height, L , [ftl 20 20 28
Relative bed height 1 1 1.4
LtI4 120 16 1.87
Regenerator
Bed diameter 4 in. 22 in. 19.5 ft.
Relative diameter 0.18 1 10.6
Relative cross section 0.03 1 112
Bed height [ft.] 20 31 37
Relative bed height 0.67 I 1.2
LtIDT 60 17 2.3
Capacity [bbVdayl 2 100 15,000
Relative capacity 0.02 1 150
Product
Gasoline [vol. %I 50.2 49.7 49.5
Gas oil [vol. %I 50.0 50.0 50.0
Dry gas [wt. %] 4.2 5.0 6.2
Carbon [wt. %] 3.2 2.8 2.8
anhydride, Betts (B11)stated that the ground catalyst easily suffers attri-
tion and produces dusty fines; these give rise to heat-exchanger fouling,
cyclone or filter blockage, and after-burning. Betts showed that the mi-
crospherical catalyst resisted attrition and eliminated the difficulties listed.
According to Graham and Way (G14), fluidized beds are ideally suited to a
moderately active catalyst because they involve contact times and large
masses of catalyst; the more highly active fixed-bed catalysts may give
reduced yield due to overoxidation and should be avoided in fluid beds.
The Fischer-Tropsch synthesis produces hydrocarbons from H2 and
CO. Fluid beds were applied to it in the Hydrocol process using ordinary
turbulent beds and in the Kellogg process using fast beds (Fig. 1). Devel-
opment of the Hydrocol process has been reported by Grekel et al. (G15)
and Hall and Taylor (Hl). The main difficulties in scale-up of these
fluidized beds were incomplete fluidization and low conversion.
On the basis of this experience with the Hydrocol process, Volk et al.
(V12) proposed providing internals in fluid beds in order to achieve good
F L U I D I Z E D CATALYST BEDS 283
OF FLUIDIZATION
D. QUALITY I N RELATION
TO PARTICLE PROPERTIES
I l l I l l 1
20 50 100 200 500 1000 2000
m i c l c size dp b m l
FIG. 3. Classification of particles by Geldart ((34, G5).
Geldart (B 1, G4, G5) has classified fluidized particles into four groups,
taking account of the fact that the behavior of fluidized beds is consid-
erably affected by properties of the fluidized particles:
Group A . Powders having a particle density less than about 1.4 g/cm3;
in particular, those with porous structure and a mean diameter in the
range of 20-100 pm.
Group B. Powders having a particle density in the range of 1.4 to about
4 g/cm3, and a mean size in the range of 40-500 pm.
TABLE 111
DESIRED PROPERTIES O F FLUIDIZED CATALYST PARTICLES
1. Effect of Fines
Particle size distribution has a great effect on most aspects of fluidiza-
tion. Zenz (24) pointed out that fluid-solid systems could be classified in
terms of a bed viscosity, which is probably controlled by the rubbing of
coarse particles against one another. According to Trawinski (T25, T26),
the fines between coarse particles act as a lubricant, reducing bed viscos-
ity. The minimum concentration of fines is the quantity required to coat
each coarse particle with a monolayer of fines. Extending this aproach to
multidispense systems, Zenz (24) calculated the optimum size distribution
of fluidized particles, shown in Fig. 4 and expressed by a skewed normal
probability distribution, which agrees with a typical size distribution of
FCC particles.
Ikedaet al. (15, IlO), investigating the effect of fines on bed expansion,
pressure fluctuation, and discharge rate of particles from a 8-cm-i.d. fluid
bed, found that bed expansion and discharge rate increase and pressure
fluctuations decrease with increasing content of fines (<44 pm). The range
of particles smaller than 44 p m is called “fines fraction”; its role has
been summarized by Ikeda (12) as shown in Table IV. The desirable prop-
erties of fluidized particles utilized as catalyst were given by Ikeda (12) as
286 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
shown in Table V. Morse and Ballow (M44), de Groot (D7) and Matsen
(M10) have also recognized that the quality of fluidization is improved by
using smaller particles with a wider size distribution. Geldart (G4, GS)
classified the fluidized beds into four categories characterized by density
difference (p, - p,) and mean size of particles, as presented in Section I.
The terms “fluid bed” and “teeter bed” (Squires, S14) seem to corre-
spond to group A and group B, respectively, although the criteria which
distinguish among the groups are not well established.
u, [cmlsrcl
FIG. 5. Expansion and pressure drop of a fluid bed (after M41).
TABLE V
DESIRABLE PROPERTIES OF FLUIDIZED
PARTICLES
~~ ~~ ~ ~~~~
Chemical Apparent contact time Stability of operation High activity: excess reaction in Dilution with inert particles
activity 3 < T < 30 sec Homogeneous temperature in dilute phase Multistage fluid bed
bed LJDT = 1-2 Low activity: decrease of reaction
yield
Particle size Weight meAn diameter Improved fluidity Severe erosion if particles are coarse Installation of internals
distribution 50 < dp < 70 pm Good fluidization Difficulty in recovering particles Multistage fluid bed with
+80 pm, 5-20 wt. % Decreased catalyst loss and (10 pm horizontal baffles
-40 pm, 10-30 wt.% attrition
Attrition Rate of attrition Avoidance of catalyst loss Very fine particles produced by Reduction of gas velocity
resistance (Forsythe Maintenance of optimum size attrition (cf. Table IV) injected through
method) distribution distributor
<O. 1-0.5
wt. %/hr
Shape of Microspherical Low tendency toward attrition High erosion rate when irregular Installation of internals
particles Decreased erosion hard particles are used Blending with MS particles
Particle Below 2-3 g/cm3 Good fluidization Unsatisfactory fluidization when Blending with low-density
density high density particles are used particles
Installation of internals
* From 12.
FLUIDIZED CATALYST BEDS 289
its volume remains constant or even decreases with increasing gas veloc-
ity; the volume of the emulsion phase is always smaller than at the bubble
point. Thus it is not correct to select the minimum-fluidization bed height
as the standard bed height for calculation of gas-bubble holdup. The
equivalent height of emulsion phase in the aggregative bed can be deter-
mined by sedimentation, extrapolating the sedimentation curve to time
zero (M41, R9). The results of Morooka et al. (M41) show that, at higher
gas velocities, the expansion ratio of the emulsion phase (L, - L,)/L, is
independent of gas velocity and bed diameter, as illustrated in Fig. 5; the
sedimentation velocity corrected by the deaeration velocity equals the
superficial gas velocity in the same manner as in a liquid-solid fluidized
bed (R7).
The expansion ratios in the emulsion phase measured by Oltrogge ( 0 5 )
and Morooka et al. (M41) at Uf >> Umfare correlated in Fig. 6 with the
parameter N,,. According to Miyauchi and Yamada (M32), the mean par-
ticle diameter in the parameter N,, is
4 = (EAn d i / x A n d j ) ” 2 (2-1)
The ratio ( L , - L,)/L, is about 0.1 when typical FCC particles are
fluidized.
loo .-
-
I I I I I 1 1 1 1 I 1 I I 1 1 1 1 J
--
Morooka rtrb(M41)
key particles
A broad PVC
-
m
4
P 8 - particles -
5 6-
-
0 -
4
* 4-
4
b
- -
Oltroggc ( 0 5 ) *\ -
-A polyprop lene
o broad FEC
o narrow F C C
B. FORMATION,
SPLITTING,
A N D COALESCENCE
OF BUBBLES
0.1
1 2 4 6 8 10 20 40
FIG. Comparison of the observed initial bubble diameter with calw-ted diameter
(after M36).
The ranges of data from which their correlation was obtained are
0.5 5 Urn,I20 cm/sec, Uf - Urn,I48 cm/sec
0.006 5 d,, 5 0.045 cm, DT 5 130 cm
The correlation of Werther (W9) gives nearly the same results as Eq. (2-5).
Both correlations predict growth of bubble diameter and decrease of bub-
ble frequency due to coalescence. However, experimental data in fluid
beds have been excluded from the above correlations.
In contrast to bubbles in teeter beds, naturally occurring bubbles in fluid
beds appear to split and recoalesce very frequently, and occasionally a
bubble frequency will increase at a height further from the distributor.
Rowe (in D5) and Toei et al. (T17) studied the mechanism of bubble
splitting. Toei et al. (T17) found that small downward cusps were gener-
ated at the two-dimensional bubble surface near the stagnation point,
which grew and thereby caused bubble deformation and, frequently, bub-
ble splitting. Figure 8 shows the effect of particle diameter on the average
total frequencies of bubble disturbance and splitting. The results of Toei et
al. (T17) indicate that the frequency of disturbances or splittings de-
creases with increasing particle diameter, and is independent of bubble
292 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
x,l
and splittings (T17).
-0.21 0.28 0.36
~ m m ~
diameter. Therefore bubbles which form in a fluid bed with fine particles
have more tendency to split than to coalesce, as discussed in Section V,C.
Morooka et al. (M43) measured properties of bubbles in fluid beds with
FCC particles. The fluid beds were 7.9 and 19.5 cm i.d., and L, was
100-200 cm. Bubble signals were detected by a hot-wire probe, originally
designed by Yamazaki and Miyauchi (Y3-3, which consists of a pair of
parallel tungsten wires 10 p m in diameter with a span of about 3 mm
and a separation of 12 mm. The electrical resistance of the hot wires
changes with the density difference between the bubble and emulsion
phase. The disturbance induced by this probe is presumably not greater
than that of the miniature capacitance probes of Werther (W6) and
Burgess and Calderbank (B17). Figure 9 shows the velocity of bubbles
ascending along the central axis of fluid beds. The effects of bubble diame-
ters, gas distributors, and axial position of the probe on the ascending
bubble velocity are rather small.
The measured bubble velocity is much greater than predicted from
equations derived by Davidson and Harrison (D3) and Taylor (TS),
respectively:
ub = uf - + 0.711(gdb)1’2 (2-7)
and
~b = 0.71 l(gd,)”’ (2-8)
FLUIDIZED CATALYST BEDS 293
I I I I I I I I I 1 1
200- L,=2 rn
I I I I
3 4 6 8 10 20 40
Uf (cmlsec)
FIG.9. Velocity of bubbles ascending along the center of a fluid bed (M43).
The reason, pointed out by Turner (T28), is that the bubble velocity in a
fluidized bed without circulation will obey Eq. (2-8), but that bubbles in
fluid beds ascend in an intense circulatory flow of the emulsion phase.
The arithmetic mean of bubble signals measured by Morooka et al.
(M43) is shown in Figs. 10 and 11. Fig. 10 illustrates the effect of superfi-
cial gas velocity on the mean length of bubble signals. Some slugs appear
in the range of U,< 30 cm/sec at the upper part of a 7.9-cm-i.d. bed.
Figure 11 shows the relationship between the mean length of bubble sig-
nals and axial probe positions for operations at higher gas velocity. If an
initial bubble diameter is larger than the equilibrium diameter, the bubble
has a tendency to split.
Actual bubble diameters can be calculated from the mean bubble sig-
nals according to the relationship given by Anderson (A6) and Ueyama
and Miyauchi (Ul). Assuming (1) the bubbles are spherical, (2) an ascend-
ing velocity is independent of a bubble diameter, and (3) the wall effect is
negligible, the following equation is obtained:
4 --e-
0
i0 FCC calalvsl
parliclrs
A
'0 10 20 30 40 50 60
Uf (cmlsecl
FIG. 10. Mean length of bubble signals in fluid beds (M43); perf. = perforated.
= (m/n) lom
P 2 M ( l )dl
(2-10)
joa
dV'(dJ d(dJ jox
F M ( 1 ) dl
where rn and n are positive integers. The actual volume-surface mean of
bubbles is calculated with Eq. (2-10) as follows:
(2- 11)
.L V
m-
1
0 100 200
Lm (cm)
FIG. 11. Mean length of bubble signals in fluid beds operated under relatively high gas
velocity (M43); perf. = perforated.
bubble growth decreased with increasing height. When the FCC particles
were fluidized at Uf = 10 cm/sec, the volume-surface mean of bubble
diameter approached the equivalent value (3.5 cm) at 40-50 cm above the
perforated plate. In a similar study, Ikeda (15) calculated' the effective
bubble diameter to be about 8 cm in a 360-cm-i.d. baffled fluid bed (equiva-
lent bed diameter = 70 cm, Uf = 45 cm/sec, U,, = 0.6 cm/sec).
Tsusui (T27) studied the effect of particle characteristics in a fluid bed of
6-cm i.d., using a hot-wire probe. The difference of flow characteristics
was most evident at high velocity of bubbles and high diameters, as shown
in Fig. 12. Fresh FCC particles showed the smoothest fluidization, ex-
pected for the high ascending velocity. Removing fine particles (d, I44
pm) from original FCC particles gave a low ascending velocity and larger
bubbles. The heavier MS catalyst showed very low ascending velocity
and a typical slugging flow. Silica microspheres, which are very light ( p ,
= 0.3 g/cm3) and average around 70 p m in diameter with a very narrow
particle size distribution, showed a fluidity similar to that of the FCC
particles without fines. The actual bubble size distributions, F(db), are
shown in Fig. 13.
All these results indicate that bubble diameters in fluid beds are strongly
influenced by the phenomenon of bubble splitting. Kehoe and Davidson
(K8, K9) observed slugging in a fluid bed with 62-pm catalyst particles,
although their results are presumably due to the narrow size distribution
of particles. An experiment conducted by Lanneau (L2) showed the pres-
ence of large bubbles in a 7.6-cm4.d. fluid bed with typical fine micro-
spherical alumina catalyst; however, his capacitance tip measured 4.8 x
6.4 x 19 mm. Since many bubbles present were probably smaller than the
296 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
--
9 -
50-
---- __-------
Eq.(2-8)
&--
- 2 -
,
- ,, r3,uz";beT' --
-
I
/
(Stewart eta/.,S21) -
/
I 1 I I
0' I
1.4 I I I I I I I I I I
1.2
1.0
n
\
5
fi 0.8
--
8 0.6
LL
0.4
0.2
0 1 2 3 4 5 6
db [em]
FIG. 13. Frequency distribution functions of actual bubble size (T27). Experimental
conditions are the same as in Fig. 12.
F L U I D I Z E D CATALYST BEDS 297
probe tip, as shown in Fig. 13, Lanneau’s measurement seems to show the
maximum size present under given conditions.
C. OPERATION
O F F L U I DBEDS
D. PROPERTIES FLUID
OF BULKFLOWINSIDE BEDS
voidage was found for fluidization of FCC particles with atmospheric air
(Z5). Bed expansion can be expressed in terms of average bubble holdup,
which also enters the calculation of contact efficiency (Section VII). Aver-
age bubble holdup is defined by
Zb = (LE - L€.)/L, (2-12)
where L, is the equivalent height of the emulsion phase measured by the
flow interruption method (R9). Experimental data by Morookaet al. (M40)
and van Swaay and Zuiderweg (V8) show (in Fig. 14) that the effect of bed
diameter on the bubble holdup is very small. At higher gas velocities,
experimental values of bubble holdup in fluid beds agree with those in
bubble columns. At U, = 50 cm/sec, z b does not exceed 0.5, and a
fluidized-bed height is at most twice that of a settled bed. In contrast with
fluid beds, the expansion of teeter beds is a function of both particle
properties and bed diameters (H17,K8, K9).
The data of Wilhelm and Kwauk (W13) plotted in Fig. 14 reveal higher
values of z b because of severe slugging in a small bed. Figure 15, given by
de Groot (D7), shows that bubble holdups in beds with narrow-range silica
are much smaller than with broad-range silica. Thus, in larger fluidized
beds of silica with a narrow size distribution, very large bubbles form with
high ascending velocities. The effect of temperature and pressure of the
I 1 I I [ I l l I I I 1
0.5 -
glass beads
c.
c. Bubble column
4
v Kato cta/.(K6)
I1
v Dr>20 cm
0.1
- 4
-
9-
0 6.6crn -
8- -
A 7.9 -
7-
0 12
6- v 19 -
0.05 - 0 /' -
1
do, I 1 1 1 1 1
I I I I &
2 4 6 8 10 20 40 60
Ur (crn/sec)
FIG.14. Holdup of gas bubbles in fluid bed (after M40).
FLUIDIZED CATALYST BEDS 299
50.':
- 1
T
I
' ' '' " ' I '
I I
1 1
10-
1 ' 1'
-
O bll 012 ' Ol, '0.'6';.;;
D, (m)
fluidizing gas on bed expansion has been studied by de Vries et al. (D8).
As Fig. 16 shows, the influence of the gas pressure is not large for wide-
range microspherical silica. These authors mention that emulsion-phase
fluidity is quite good at high temperature.
8 T = 300'C
c
'I. T= 20.C
40 r
30
1 .o
I "
0 0.5 1.5
Pabsr bar
FIG.16. Fluid bed expansion for SCP catalyst as a function of pressure (D8). Fraction
t44 p n 20%, U,= 18 c d s e c , and DT= IOcm.
J
300 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
FIG. 17. Lateral holdup distributionof gas bubbles in fluid bed with 12-cm inner diameter
(NV.
Fig. 17. The distribution is nearly independent of axial position. The hold-
up of gas bubbles at radial distancer from the center axis of the bed is
expressed as
(EbO - Eb)/EbO = [(EbO - E b ~ ) / E b l l l ( T / R ) ~ (2-13)
where ebOis the holdup of gas bubbles at the center axis of the bed, and Ebw
is the holdup of gas bubbles extrapolated to the bed wall and may be taken
as approximately zero. Parameter n in Eq. (2-13) is variable with gas
velocity and bed diameter (Figs. 18 and 19). Their experiments were
Uf Ccmhccl
FIG.20. Circulation velocity ascending at the center of fluid bed and bubble column.
Data of bubble columns are measured by: 8,Pavlov (P3), DT = 17.2 cm; €B, Pozin et al.
(P6), 30.0 cm; A,Yamakoshi (Y2), 25.0 cm; curve, Yoshitome and Shirai (Y22), 15.0 cm;
0 , Miyauchi and Shyu (M31). 10.0 cm; 0 , Hills (H9), 13.8 cm; V , Ueyama and Miyauchi
(U3), 60.0 cm; 4, Koide et al. (KIS), 504.0 cm; @, Kato and Morooka (unpublished data),
12.0 cm.
FLUIDIZED CATALYST BEDS 303
(2-16)
(2-17)
Mathesonet al. (M7), show that relatively small amounts of fines added to
coarse particles sharply decrease the bed viscosity. Based on this evi-
dence, Zenz (24) has calculated an optimum size distribution of fluidized
particles (as mentioned in Section II,A,I), The importance of the apparent
viscosity in relation to flow characteristics of the bed has been stressed by
Matheson et al. (M7), Rice and Wilhelm (R5),Finnerty et al. (F3), Grace
(G13), and others. Matheson et al. (M7) have found that the slugging
tendencies of a bed can be expressed in terms of the apparent viscosity.
Their results for FCC particles containing 20 wt.% of 46-pm-diameter fine
particles gave p b equal to 6 X g/cm . sec; this small bed viscosity
agreed with the visible high fluidity of the bed.
Meanwhile, Rice and Wilhelm (R5)have investigated the stability of
interfacial surface between the particle-free gas phase and the dense-bed
phase. According to their analysis, the lower bed interface (or upper
interface of a bubble) is inherently unstable, whereas the upper bed inter-
face (or lower interface of a bubble) is inherently stable. Thus an instabil-
ity initiated at the upper bubble surface may reduce bubble size; this
tendency is presumably stronger in a bed of lower viscosity. The apparent
viscosity which applies to the surface dynamics has been studied by Fin-
nerty et al. (F3), Takamura (T2) and Morooka and Kato (unpublished
data). They have measured wavelengths, frequencies, and attenuation
rates for waves generated at the bed surfaces. In the experiments of
Takamura (T2) and Morooka and Kato (unpublished data), rectangular
fluidized beds of 6 x 60 cm and 8 x 80 cm were utilized, respectively, and
bed heights were kept at about 30 cm. The attenuation rates of surface
waves were calibrated with glycerol-water solutions. Their results, shown
in Fig. 21, indicate that the apparent viscosity of FCC particles is the
smallest among particles employed in the experiments.
FCC particles
4- dp=60~m
$=l.Og/cm
I
-
E
2- ‘a,
&
Y
D
Y a
FIG.21. Apparent viscosity of fluidized bed (T2 and Morooka and Kato, unpublished).
FLUIDIZED CATALYST BEDS 305
E. AXIALDISTRIBUTION
TO BED STRUCTURE
1. Bed Density
With increasing gas velocity, the bed surface diffuses gradually, and a
dilute phase is formed on the dense bubbling bed owing to entrainment of
fluidizing particles. This axial distribution of bed density has been mea-
sured by several workers (B4, F2, L13, M25, M33). Fig. 22 shows the
axial distributions of bed density measured by Miyauchi et al. (M33) in
1969. Their experiment was performed in a 7.9-cm-i.d. fluid bed of FCC
particles having a mean diameter of 58 pm.
The influence of the amount of particles in the bed on the axial porosity
distribution has been reported by Bakker and Heertjes (B4) for glass
beads of 175-210-prn-diameter, fluidized in a 9.0-cm teeter bed. Their
results show that the quantity of axially suspended particles in the dilute
phase is increased noticeably by increasing the initial bed mass from 200
to 2000 g. In contrast with this behavior, the bed density distribution given
by Kaji (Kl) for the dilute phase shows less dependence on the initial bed
height Lq under constant Uf. He has measured the distribution for a 7.9-
cm-diameter bed of FCC catalyst, having a mean size of 58 pm, by vary-
ing L , at 32, 51, and 71 cm, and Uf for 27-58 cm/sec. The mean bed
density in a dense fluid bed can be calculated by the following equation
knowing the mean bubble holdup in the dense fluid bed given in Section
9 0.3
n
I
E
20.2
FIG.22. Longitudinal bulk B
density distribution (M33). 2
0.1
0
50 100 150 200
Longitudinal position [cml
306 M I Y A U C H I , F U R U S A K I , MOROOKA, A N D I K E D A
FIG.23. Entrainment rate from the bed with cracking catalyst particles.
0 Modified lnternals
I I I I
2 4 6 8 10 12 14
Equivalent Diameter, in.
cates that the equivalent diameter of a large reactor gives the same con-
version as a small reactor of the same actual diameter. The proposed
equivalent diameter is
(Deq1-I = (DeV1-l + (DeH1-l (2-20)
where Devand DeHare the component equivalent diameters based on the
vertical surfaces and the horizontal surfaces, respectively. Ikeda carried
out an acrylonitrile synthesis by the Sohio process, in fluid beds with
diameters from 8 cm to 3.8 m. Reaction yields were compared under the
same reaction conditions and catalyst, and the effective bubble diameter
was calculated by modifying the bubbling bed model originally proposed
by Kunii and Levenspiel (K22, K23). The effective bubble diameter is
given by
(db)eff = 1-9(De~)''~ (2-21)
where (db)eff andDeq are both in cm.
Several investigations have been made on the flow characteristics in
multistaged fluid beds. Nishinaka et al. (N6, N8, N9) have measured the
average bubble holdup, the lateral distribution of bubble holdup, and the
longitudinal dispersion of solid particles in four- and eight-stage fluid beds
installed with various horizontal baffles. As shown in Fig. 25 the average
bubble holdup (except for beds baffled with tube plates) is correlated by
the equation of Nishinaka et af., (N8):
Eb = 0.08 ufo.73A-O.3NO.O5
r (2-22)
where A, is the free area of baffle in %, and N is a number of horizontal
plates; Uf is in cm/sec. The longitudinal-dispersion coefficient of solid
particles has been found to decrease with decreasing free area of baffles
(N6). The intermixing mass velocity of solid particles between adjacent
FLUIDIZED CATALYST BEDS 309
FIG.25. Correlation of average gas bubble holdup in free and baffled fluid beds (N8).
FIG.26. Relationship between intermixing velocity of solid particles and gz'/V: (N8).
rate between the bubble and the emulsion phase in baffled fluid beds (M42)
is mentioned in Section V1,B.
use the same equation of motion, but different assumptions for the turbu-
lent kinematic viscosity.
The time-averaged flow pattern of liquid in a bubble column may be
visualized as shown in Fig. 27, where the column is in recirculation flow,
and the ratio LID, is large compared to unity. Under steady flow, the
basic equation of motion (H9, M31) is:
-(l/t') d(m)/dr = dP/dz + (1 - Eb)Plg (3-1)
where the left-hand side of the equation gives the force acting on a fluid
element by radial shear-stress difference, and the right-hand side gives the
forces acting by axial static-pressure gradient and by gravity. A similar
equation of motion has been applied by Wijffels and Rietema (W16) to
determine the flow properties of liquid-liquid spray columns.
The shear stress T is related to the mean axial velocity of liquid through
the sum of the molecular and turbulent kinematic viscosities:
T = -(VM + Vt)p,(dU/dr) (3-2)
where vMis usually negligible in comparison with v, except in the vicinity
of the column wall. The turbulent shear stress for the bubble column with
L / D , >> 1 is expressed by the following equation (U2):
-vt dU/dr = (1 - Eb)ULU: (3-3)
where U: and u: are the radial and axial velocity fluctuations of liquid,
respectively. Equation (3-3) is obtained by time-averaging the fluctuation
velocity components for two-phase flow.
The boundary conditions to solve the basic equations, Eqs. (3-1) and
(3-2), were given earlier (M31). As shown in Fig. 27, the liquid flow under-
goes upflow between a and b and downflow between b and c, both with a
developed turbulence, and also downflow in the laminar sublayer between
c and d near the wall; at the wall, r = R (or y = R - r = 0) and u = 0.
Consistent with the concept of the universal velocity profile (S4, SS), a
buffer layer should be present between the turbulent core and the laminar
sublayer. However, it is hard to distinguish the laminar-turbulent buffer
layer from the turbulent core, since bubbles from the turbulent core enter
irregularly into the buffer layer and agitate it. Hence, partly to simplify the
mathematical treatment and partly to satisfy a physically sound interpreta-
tion of this bubbling buffer layer, the thickness 6 used for the laminar
sublayer is taken somewhat greater than in single-phase flow and the buffer
layer is neglected.
According to the concept of the universal velocity distribution, the
velocity profile in the laminar sublayer is given by
uIv, = v.y/vM for v , y / v , 5 5 (3-4)
FLUIDIZED CATALYST BEDS 313
r -0 r =R
y =o
FIG.27. Recirculation flow pattern in a gas bubble column.
where u. = ( / ~ ~ l / pisJthe
~ ’ friction
~ velocity. The velocity distribution in
the turbulent core is
u/u, = 5.75 log(u,y/u,) + 5.5 for u,y/uM L 70 (3-5)
By equating Eqs. (3-4) and (3-9, one obtains a valuey defined as the
thickness 6 of the laminar sublayer, for which:
u,S/uM = 11.63 (3-6)
The velocity us at y = 6 is given by combining Eqs. (3-4) and (3-6):
ug = 11.63~.= 1 1 . 6 3 ( ~ ~ w ~ / ~ J 1 ’ 2 (3-7)
According to the analysis of many experimental data (M31, US), the
thickness 6 is usually much smaller than the column radius R , so that us
essentially equals the peripheral velocity luwl of the turbulent core. As a
314 M I Y A U C H I , F U R U S A K I , MOROOKA, A N D I K E D A
By introducing Eqs. (3-2) and (3-10) into Eq. (3-1) and neglecting vM in
comparison with vtrthe following basic equation is obtained for the turbu-
lent core:
-(I/r)d[v,r(du/dr)l/dr = (2/R)(TW/Pl) - (Eb - Eb)g (3-12)
This equation is solved with respect to u under the simplifying assump-
tions (M31, US) that vt is constant in the turbulent flow region, and that
bubbls gas holdup Eb is distributed radially according to the relation:
Eb/<b = 2(1 - 4') (3-13)
where 4 = r / R . The physical significance of the assumption regarding vt
will be explained later; trial calculations to find a radially variable vt giving
the best fit to the experimental velocity profile of the bubble-liquid mixed
phases have shown that a constant vt is quite workable if taken as the
radial-cross-sectionally averaged mean value of a distributed vt. This as-
sumption simplifies the mathematical treatment, and also is supported by
Figs. 17 and 18 (refer to Section II,D,2). Becausen = 2 the gas holdup QbO
at the center axis is twice Eb. The empirically adjustable viscosity term vt
has different values when n is changed, but the flow properties of the
FLUIDIZED CATALYST BEDS 315
The numerical value of vt will be given in the next section. In Eq. (3-14),
the first term of the right-hand side has been shown by Ueyama and
Miyauchi (US) to be essentially negligible for the recirculation flow re-
gime. As a consequence, one has the following simplified quantitative
relation:
@Y2 (3-15)
The undetermined term Iu,I is given by applying a mass balance for the
liquid phase:
(3- 17)
where the 2 sign has the same meaning as that in Eq. (3-16). With this
Iuwl,the radial distribution of linear liquid velocity u is given by Eq. (3-15).
Another simplified liquid velocity profile is also obtained directly from
Eq. (3-15):
(u + lu,l)/(uo + lu,l) = ( 1 - @Y (3-18)
where uo is the linear liquid velocity at the column axis, + = 0, given from
Eqs. (3-15) and (3-17) as:
(3-19)
+*
The radial position r , (or = r J R ) where the liquid velocity u relative to
the column wall is zero (point b in Fig. 27) is obtained from Eqs. (3-15) and
316 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
UG/<b f u~/(1
- Eb) = us + ( U o + 1~,1)/6(1 - iib) (3-23a)
With uo + Iu,I given by setting C#) = 0 in Eq. (3-15), Eq. (3-23a) becomes
uc/<b k uL/(l - Eb) = 6 s + (gDT/192~t)Eb/(l- Eb) (3-23b)
where the f sign has the same meaning as in Eq. (3-16). The second term
on the right-hand side of the equation is the apparent increase of slip
velocity due to recirculation flow. [Calculation shows that this term equals
U*/Eb,where U* is the net linear velocity of the recirculating liquid aver-
aged over the column cross section. '] This term is a natural consequence of
recirculation flow, and does not appear in the bubble-flow regime where
the following well-known relationship applies (N4,Y20):
uL/(l - E b ) = 6 s
uG/Zb h (3-24)
The mean gas holdup Eb in the recirculation flow regime is obtained at
UL = 0 by solving Eq. (3-23b) with respect to gb:
) "'1 (3-25)
U* =
I,"
(~rR*)-l 2nru dr.
FLUIDIZED CATALYST BEDS 317
where
(Y = U,/ii, and y = gDz,/(192 u,U,)
The reason & is approximately proportional to U$z and decreases only
slightly with column scale-up in recirculation flow (cf. Fig. 14) is ex-
plained from Eq. (3-25) by including the behavior of Us and vt, which are
discussed in the next section.
The mean flow rate o f the recirculating liquid is useful when known. For
upflow, we take the superficial value:
$ = r/R
0 0.1 02 0.3 0.4 0.5 Ok-I
profile
0.1
'0 0.l 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
fi r/R [-I
FIG. 28. Interstitial liquid velocity profile and normalized velocity profile [data by
Yamagoshi (Y2)],D, = 25 cm, I!/= , 5.2 cd s ec, air-water system at room conditions.
Adjusted u,, and I u, I are 45 and 26 c d s e c , respectively.
FCC-Fluidized bed
0
Mean gas holdup 5 [-I
FIG.29. Radial position 4%as a function of ib.
The last proportionality comes from his experimental test of the relation.
Along the column axis, u equals uo and is approximately proportional to
(U&)1’4, as is seen in Fig. 30. Hence, vt a Dq’4/V24along the column
axis.
The constant effective vt assumed in the present simplified treatment is
obtained by knowing uo and IU,~. Along the column axis ( 4 = 0), Eq.
(3-15) gives
u0 + IuwI = gPT$/(32vt) (3-30)
With uo and Iu,l obtained by trial and <b from Eq. (3-28), vt is given by Eq.
(3-30). In the case of Fig. 28, Eb is calculated to be 0.123 (the measured
value is 0.127) and vt = 33.2 cm2/sec.
In the same manner, vt has been calculated by Ueyama and Miyauchi
(US)from a group of available experimental velocity profiles. Fig. 31
shows the recalculated vt as a function of DT, matched empirically by the
power function:
vt = 0.1280$’0 (3-3 1)
150
100
’ J 7
$5;
- 3
J
2
correlation based on this dependence. The line B-B is for the bubble
columns and is expressed (in cm-sec) by
v1 = 0.160Lp,'6D3~2 (3-32)
The correlation shown in Figs, 32 and 33 as curve A-A is the turbulent
kinematic viscosity of the emulsion phase of fluidized cracking catalyst
beds estimated by an indirect method explained in Section IV,C in rela-
tion to axial dispersion of the emulsion. The correlation is given in cm-sec
units by
vl = 0.079U&'6Dk70 (3-3 la)
The exponent 1.70 on DTis different from 1S O of Eq. 3-32, but Eq. 3-3 la
gives a better correlation for axial dispersion of the emulsion.
We have seen that vt can be obtained from Iu,I or uo. In the case of Fig.
28, ul = 28.1 cm*/sec from Eq. (3-17), and vl = 36.1 cm*/sec from Eq.
(3-19). Some data are available for uo (H9, P3, p6) from which vt are
calculated and shown in Fig. 33. In general, vt as calculated from uo tends
to give a larger value than from uo + IuwI,perhaps due to lack of suffi-
0, Fml
FIG. 32. Turbulent kinematic viscosity u, as a function of DT and U,. Line A -A is for
D, [cml
FIG.33. Turbulent kinematic viscosity v, as determined by different methods: the lines
A-A and B-B are the same as those given in Fig. 32.
FLUIDIZED CATALYST BEDS 323
- 10’
<X ?
u E
u 5
4
3
I3
2
10
FIG. 34. Dependence of average slip velocity ic, on the superficial gas velocity UG;
column diameter DT has little effect on 4. A porous plate distributor seems to give
smaller &.
324 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
(3-34)
l -
D, = 19.0 cm
& = 0.054U$z
UG2 20 cm/sec
100
9, = 141 cmysec
~ = 4 6crri/sec
00
FIG. 35. Simultaneous determination of Y, and Us0 from Pb based on Eq. (3-34).
U, tcrn/secl
2 3 4 5 7 10 2 3 450 7 100
U, [cm/sec]
FIG. 36. Comparison of the calculated &, by Eq. (3-25) with the data given in Fig. 14.
Unless otherwise mentioned, U, = 49.5 c d s e c . Data are for FCC-catalyst beds of different
bed diameters. Curve E-E is an average of the data. The calculated curves BQB and BQP
are for bubble columns, FQF and FQP are for FCC-catalyst beds.
326 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
c. PHENOMENOLOGICAL BACKGROUND
OF TURBULENT
KINEMATIC
VISCOSITY
The turbulent kinematic viscosity v, has been introduced in the basic
relation, Eq. (3-2), and is correlated with the operational variables through
the experimental data (cf. Figs. 31-33) with a simplifying assumption that
vt is constant radially as an experimentally adjustable parameter.
It is desirable to establish a phenomenologically sound basis by relating
v, to well-established fluid-dynamic concepts, since we need to know the
possibility of scale-up or of applying the concept of vt to a fluidized
catalyst bed, as discussed in what follows (M27).
Thus far the experiments considered have been confined to bubble col-
umns, which provide a fair number of data over a wide range of conditions
and are therefore useful for testing the simplified theory. Equivalent data
328 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
are often not available for fluidized beds of carefully prepared particles,
especially with solids other than FCC catalyst.
The main uncertainty in applying the simplified theory to a fluidized
catalyst bed probably lies in the question as to whether the shear-stress
term in Eq. (3-10) is negligible in comparison with the buoyant-force term.
In what follows, although we still lack direct experimental proof regarding
the shear-stress term, the simplified theory will be seen to account well for
the bed performance (M27).
t o ) 4 5 7 10 2 3 & S O 7 100
FIG.37. Mean bubble velocity ubo along the column axis; experimental data are for
FCC-catalyst beds. Full curves are calculated by UbO = + u owith 4 = 49.5 c d s e c anduo
from Eq. (3-19). Umr is neglected due to Urn,<< U,.
FLUIDIZED CATALYST BEDS 329
from Eq. (3-31a), theu,, which results is shown in Fig. 37. The calculation
matches the data for columns with diameters of 5.2 and 19.5 cm, but not
for those with diameters of 7.9 cm, for reasons which are not obvious.
In Fig. 37 the bubble velocity in slugging beds [cf. Eq. (5-4b)l is also
shown. The FCC-catalyst bed of 5.2-cm diameter is fluidizing smoothly in
the recirculation-flow regime with no indication of slugging when L, is
smaller than 50-60 cm.
By reversing the above procedure of obtainingu,, by knowing iis,Eb and
vt, one can calculate vt by knowing fi,, c b and ubo. These vt are shown in
Fig. 33 as open circles. They show reasonable agreement with the correla-
tion A-A for fluidized catalyst beds, where the data for 7.9-cm diameter
with U, < 30 c d s e c are omitted.
1’1 D~ = 19.0 cm
Fluidized RI: bed
-
&b VS. UG Fig,36
u,2 7 cm/sec
P
FIG.38. Simultaneous determination of vt and &, from mean gas holdup data by using
Eq. (3-34).
FLUIDIZED CATALYST BEDS 33 1
A. THEORY
O F LONGITUDINAL
DISPERSION
FOR THE RECIRCULATION
FLOWREGIME
A gas bubble column is taken here as a model equipment undergoing
longitudinal dispersion of the continuous phase. The theory obtained is
equally applicable to a fluidized catalyst bed of good fluidity exhibiting
similar flow properties. The following procedure is from Miyauchi (M27).
1. Underlying Concept
In the bubble column the velocity profile of recirculating liquid is shown
in Fig. 27, where the momentum of the mixed gas and liquid phases
diffuses radially, controlled by the turbulent kinematic viscosity vt. When
UL = 0 (essentially no liquid feed), there is still an intense recirculation
flow inside the column. If a tracer solution is introduced at a given cross
section of the column, the solution diffuses radially with the radial diffu-
sion coefficient E , and axially with the axial diffusion coefficient E,. At the
same time the tracer solution is transported axially b y the recirculating
liquid flow. Thus, the tracer material disperses axially by virtue of both
the axial diffusivity E, and the combined effect of radial diffusion and the
radial velocity profile.
The latter mechanism assumed is the well-known Taylor dispersion
(T9, TlO, T l l ) , which has been studied extensively ( A l l , G6, L9, T14,
S2). High-speed motion pictures taken by Towell et al. (T23) in a 40-cm
bubble column (R3) have shown the presence ofturbulent eddies, on a scale
roughly equal to the column diameter, with systematic large-scale circula-
tion patterns superimposed. Their pictures showed that liquid near the
wall flowed downward, while liquid near the center of the column flowed
upward, consistent with the flow theory developed earlier and with the
Taylor dispersion mechanism.
Axial dispersion in the bubble column is usually well expressed by the
following one-dimensional diffusion model with respect to liquid concen-
tration c:
ac/ae = E~ a2C/a~2- ii ac/az - aqc) (4- 1)
2. Basic Relations
The equation of continuity for the liquid phase of the bubble column is
obtained from a mass balance on tracer material for differential gas-liquid
mixed phases. With the same procedure as for a homogeneous flow, the
following equation is obtained when @ ( c ) = 0 and the aspect ratio
LID, >> 1:
-
r ar
[r( -EreL$)] - ueL ac = cL ac as (4-3)
Here, the axial diffusion term, cLEr d2c/dz2,is tentatively omitted for
simplicity, and the radial mass flux N, (with v, = ELE~for simplicity) is
defined by
N , = -eLEr ac/ar = -u, a c / a r (4-4)
When the reference plane moves with mean net liquid velocity of
kUJ(1 - G,,) (+ indicates cocurrent and - countercurrent), Eq. (4-3)
remains valid; u , the liquid velocity relative to the moving reference
plane, is given from Eqs. (3-15) and (3-17) by setting UL= 0. This is
justified because the net flow term, which defines the moving reference
plane, cancels with the net flow term included in the linear velocity term
u , [see Eqs. (3-15) and (3-17)]. In what follows in this section we takeu as
velocity relative to the moving plane.
We now consider the case of steady back-mixing experiment to mea-
sure the total axial diffusivity EzT. A small amount of fresh liquid is con-
tinuously introduced into the column cocurrently or countercurrently with
the gas flow. A still smaller amount of tracer liquid is steadily introduced
at a downstream section. Under such circumstances the concentration
distribution of the tracer at a given cross section is steady, and the axial
concentration gradient a c / d z is constant radially. As a consequence, Eq.
(4-3) simplifies to
-- -(
i a -rv,
r ar
5) - ucL ac- oz- (4-5)
The axial dispersion coefficient E,, given by the Taylor mechanism for
the above case is given by the procedure of Tichacek et al. (T14) for
turbulent pipe flow. Integrating Eq. (4-5) twice with respect to r, one has
The net transport Q past any reference plane moving with liquid at
the mean velocity +U,/(l - Eb) is
relate vt for the columns with data for atmospheric turbulence taken by
Richardson (R6, 01). This approach gives vt approximately equal to 0.12
14/3~$3in cm-sec units; for atmospheric turbulence E,,, is assumed to be 5
cm2/sec3(01) and 1 one-half the distance between pairs of observational
points; for the bubble column 1 is assumed to be one-half the column
radius. Physically the Prandtl hypothesis may be a sounder approach.
B. LONGITUDINAL
DISPERSION
I N A BUBBLE
COLUMN
Many data are available for the longitudinal dispersion coefficient of the
bubble column in the recirculation flow regime. Most are summarized in
Fig. 39. Of these, Kato and Nishiwaki (K6, 0 and 0),Ohki and Inoue
UC [cm sec]
(02), Towell (T24, @), and Hikita and Kikukawa (H11) have utilized a
pulse response technique or step response; and the rest of the data have
been obtained by the steady back-mixing method. These are correlated as
follows (M27).
(4-12)
gas velocity. Also, the gas distributor construction has little effect in the
recirculation flow regime, provided the distributor holes are placed sym-
metric around the column axis.
Eq. (4-12) is simplified further by approximating the term in parentheses
as a function of U G .In cm-sec units this leads to
ExL = 1.35UjjzPdz (4- 13)
for U G = 5-60 cm/sec. For DT = 10-30 cm Eq. (4-13) is modified to
UGDT/E,L = 2.0( UG/(gDT)1/2)112 (4-14)
which is nearly the same as the dimensionless correlation by Kato and
Nishiwaki (K6) for UG/(gDT)*’z2 0.03. Also, Aoyama et al. (A9) report
that E Z L is proportional to Pdzfor the columns with a sieve plate as a gas
distributor. An empirical correlation by Beduraet al. (B7) gives E z L equal
to 2.4uOG.33D$4.
Some correlations are available for the apparent axial diffusivity E L E z L .
Towell and Ackerman (T24) have given ELE,L = 1.23 Uj12PJ2. Deckwer et
al. (D14) obtain the same relation with 2.0 in place of 1.23. Towell and
Ackerman’s correlation is reduced to ErL = 1.5 L13j8P2zin cm-sec units.
Equation (4-12) explains the axial dispersion coefficient reasonably well,
as shown in Fig. 39, where the range of experimental variables is DT 2.I0
cm, the perforation hole size 2 2 mm and UL 5 2.18 c d s e c , and the liquid
is water or aqueous solutions with the liquid viscosity between 0.86 and
14.5 cp. and the surface tension = 33-75.5 dyne/cm.
At low gas rate, the size of gas bubbles arid the gas hold-up are fairly
sensitive to the gas-distributor structure (A9, F7,K17, 0 2 , T15). Mixing
and bulk recirculation of liquid are induced by the interaction of liquid and
bubble motion, so that axial dispersion is influenced by the structure of
gas distributor and may show strange behavior (A9, 02). For example,
Ohki and Inoue (02, Fig. 9) show that ErLincreases rapidly with U Gto a
maximum at UG = 5 c d s e c , then decreases gradually, and approaches
typical recirculation-flow behavior regime at around UG= 16 c d s e c . In
this example a sieve plate with 91 holes of 0.4-mm diam is utilized for a
ldcm-diam column. Aoyama et al. (A9) found maximum and minimum
Ex,, in the low-gas-flow rate region of their 5.0-cm-diam bubble column
using a sintered porous plate for gas distribution. These anomalies are
closely related to the properties of bubble holdup in the transition from
bubble flow to the ultimate circulation flow with increasing UG.The prop-
erties of axial dispersion in these flow regimes may be explained by the
Taylor dispersion mechanism (T9, A1 1); thus,
(4-15)
338 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
C. LONGITUDINAL
DISPERSION BED
CATALYST
I N A FLUIDIZED
Most data available from past studies are summarized in Figs. 40 and 41
for the longitudinal dispersion coefficient of the emulsion phase in
fluidized catalyst beds of good fluidity. Such coefficients were first mea-
sured in _the pioneer work of Bart (1950, B20) for a wide range of gas
velocity ( V , = 7.5-90 c d s e c ) . His values with a 3.2-cm-diam column are
approximately one-half those by Morooka et al. (M40) in the higher flow
rate region. The basis for their correlation follows (M27).
for v, of the fluidized catalyst bed are not yet available from flow
properties.
For the FCB we adopt the notation E, instead of E,.,. Experimental data
for E, have been processed to know how they relate to DT under a given
UGand to UGunder a given DT. With these relations in hand, Eq. (4-1 1) is
utilized to find the turbulent kinematic viscosity vt of the fluid catalyst bed
as a function of DT and U,, applying the 9,. and qr utilized for the bubble
columns [Eq. (4-12)]. The vt thus obtained has been given in Fig. 32 as
curve A-A, matching Eq. (3-31a). As is evident from Fig. 32, the turbulent
kinematic viscosity for a fluidized catalyst bed with good fluidity is ap-
proximately the same as for the bubble column, and this suggests strongly
that the two systems behave similarly. In the recirculation flow, for ex-
ample, the emulsion phase should be in intense turbulence. The bubble
sizes in this flow regime should be similar for the two cases as is found in
Fig. 34.
Introducing Eq. (3-3 la) for vt into Eq. (4-1 I), again taking q 7 and q z
from Eq. (4-12), the longitudinal dispersion coefficient E , of the emulsion
phase is given in cm-sec units as follows:
E%s = (EzT)ernulsion
11
'I
w
5
4
)1
1oz
7
U, [cm/sec]
@) and those from de Groot (D7, V8, $) deviate from those for FCC
particles at smaller or larger bed diameters, respectively, the general trend
is nearly parallel.
The turbulent kinematic viscosity u, of the fluidized catalyst bed has
been determined, as Eq. (3-3la), from the use of axial dispersion coefficient
E,. This is a natural consequence of the analogy between the bubble
column and the fluidized catalyst bed of good fluidity (such as in fluidized
catalytic cracking). The mean gas holdup (Fig. 36) and the mean'bubble
velocity along the bed axis (Fig. 37) are reasonably well predicted by
applying Eq. (3-31a) for the fluidized cracking catalys! bed.
The bubble holdup Eb for a bubble column or for a fluidized bed is easily
obtained as a function of superficial gas velocity U G .Useful information
FLUIDIZED CATALYST BEDS 34 1
lo'
n 5
. de Vries et al. ( 8
lo2
Bed diameter D, [cm]
FIG. 41. Influence of particle properties on axial dispersion of the emulsion phase;
U G = 20 cm/sec with variable bed diameter D T ;keys are given in Fig. 40. Data by May
(U,= 25 cm/sec) and those by Stemerding (U,= 10 cm/sec) are extrapolated to UG = 20
cm/sec according to Fig. 40.
has been extracted from the Zb versus V , relation for the fluidized bed
(D3, N4, P8, T29). In recirculation flow an additional term appears in the
apparent mean slip velocity shown on the right-hand side of Eq. (3-23b).
This equation has been utilized to determine the turbulent kinematic vis-
cosity vt and the mean slip velocity ii, (cf. Sections 111,B,2, III,B,5, and
III,D,5). In this chapter, the equation is further examined in relation to
bed performance, since the turbulence properties of the bed result from
interaction between bubbles and the continuous phase. As shown in Fig.
34, the mean slip velocity of bubbles in a fluidized catalyst bed of good
fluidity is essentially the same as that for a bubble column when UG k 20
c d s e c . A criterion under which bubble size approaches a dynamic
equilibrium is obviously needed for predicting or evaluating the perfor-
mance of scaled-up beds.
through the emulsion phase for a fluidized catalyst bed. This one-
dimensional bubble flow has been discussed extensively (B14, D4, N4,
27). A bubble column without continuous liquid feed ( U , = 0) is taken for
simplicity.
- Free
9,
A --
gb
"b
$ Gas
1
flow
Gas flow
is stopped
(a) (b)
FIG.42. Example of the rise of swarm of bubbles: (a) for continuous bubbling; (b) for a
rising swarm of bubbles.
F L U I D I Z E D CATALYST BEDS 343
solid. The rate of downflow of solid will be a maximum when the wake is
formed at the gas inlet and released at the free surface without being
renewed during the travel; also, the rate of downflow will be essentially
zero when renewal is very fast during the travel. The actual situation lies
between these two extremes (R17); and essentially the same situation
applies in the gas bubble column.
We now consider the extreme case of negligible renewal of bubble wake
for a continuously bubbling bed. The wake fraction in the bed isfwZh.The
bubble phase is composed of a bubble and a wake, and the fraction of this
phase is (1 + fw)& In the bubbling-bed model (F12, K24) the bubble
phase is assumed to ascend at velocity ii, relative to the emulsion phase.
The rate of release of wake solid at the free surface of the bed is then
fw&iih, where iib is the rising velocity of the bubble phase relative to the
bed wall. Released solid constitutes the downward return flow through the
emulsion, its velocity being given by f&uh/[ 1 - (1 + f,).h] This return
flow retards the ascending motion of the bubble phase, and leads to the
following relations:
fib = - fwEbuh/[ 1 - (1 + fw)eb] (5-6)
or
fib = (u, - umf)/g = [l - fwEh/(l - Eh)]& (5-7)
As a consequence i i h is always smaller than the slip velocity ii, for the
extreme case of negligible renewal of the wake and equals Us whenf, = 0
or when the rate of renewal is very rapid for finitef,. In the latter example
the wake loses its identity relative to the emulsion.
In conclusion, for a homogeneously bubbling bed, the velocity of bub-
ble rise i i b relative to the bed wall is equal to or lower than the bubble slip
velocity ii, relative to the continuous phase.
for the bubble column; Us is the mean slip velocity of bubbles, equal to the
free-rise velocity Us, when the interaction between bubbles is negligible.
Towel1 et al. (T23) give p = 2 for their 40-cm-diam. bubble column under
fairly high speed aeration. In Davidson and Harrison's well-known equa-
tion (D4), p equals unity with ii, = iiso, although their equation is derived
on the basis of homogeneous bubble flow.
lo''* +
pb = ( T R z ~ , ) - ~ 21rr(u U&b dr
=I-(I--) u;/2 .1- 1.5ib
uG/cb 3(1 - Zb)
Hence, the simplification of 100% is not unreasonable.
346 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
When the wakes are renewed rapidly, they have little influence on bub-
ble motion (Section V,A,2).A relatively small influence of wakes is par-
ticularly the case when bubble motion is in the recirculation-flow regime
(Section V,B, 1).
The influence of wake motion on bulk turbulence induced in the liquid is
understood more clearly by inspecting oscillograms which show the fluc-
tuation of local liquid velocity. Figure 43 shows such oscillograms taken
by Kikuchi (K30) with a hot-wire probe. The bubble column is 8.0 cm in
diameter, water-filled to a 170-cm height, with the probe 115 cm above the
bottom gas distributor. In the column bubbles of constant volume (100 cc)
are injected successively at constant time intervals, either at 2.1 sec (case
a) or 0.50 sec (case b). Case c is an example of continuous bubbling at
UG = 6.45 c d s e c .
In case a, wake motion is obviously identified for each bubble passing
over the probe. The local velocity of liquid changes periodically in reso-
nance with the bubble frequency. The resonance, however, is consid-
erably disturbed as the bubble frequency increases. The oscillogram
shows the onset of irregular liquid motions, which have different frequen-
cies from those of bubbles. The relation between wake motion and local
liquid motion becomes quite dubious. This tendency is increased in case
b, where liquid motion exhibits fairly low frequencies. These flow proper-
ties are succeeded by continuous bubbling (case c). Here the local liquid
motion is governed by the properties of bulk liquid motion, rather than by
the bubble wakes.
Based on these discussions, the influence of bubble wakes has been
omitted in what follows, since the wake motion is difficult to distinguish
from bulk turbulence. However, the question of wake effects has not yet
been settled experimentally.
The second integral is the net liquid flow through the column and is zero
FIG.43. Hot wire oscillogram showing the fluctuation of local liquid velocity; air-water
system, DT = 8.0 cm,S = start of bubbling; data by T. Kikuchi. (a) Time interval is 2.1 sec
for bubble injection. (b) Time interval is 0.5 sec. (c) Continuous bubbling at UG = 6.45
cdsec.
FLUIDIZED CATALYST BEDS 349
The data necessary for calculating p using the above equations have al-
ready been given in Section 111.
Figure 44 shows p as calculated by Eq. (5-16) as a function of UGand
DT. Generally, p changes with these two parameters. ForD, = 40 cm, p is
UG km/sec]
FIG. 44. p as a function o f UGand DT for a bubble column, where UJZb = U s + PUG.
350 MIYAUCHI, F U R U S A K I , MOROOKA, A N D IKEDA
c. BUBBLE
SPLITTING IN FLOW
TURBULENT
The mean slip velocity Us has been shown in Fig. 34 for the bubbles
observed in bubble columns or in fluidized FCC beds. Here the mean
0' 10 20 30 40 50 0
( 0- &I 1 rcm / =I
FIG.45. (UG - Umf)/Zbas a function of (VG - Urn,)for fluidized beds. FCC-catalyst
beds show no indication of slugging.
FLUIDIZED CATALYST BEDS 35 1
&l I 11.0(cmlsac1
Z (cm)
----
- 38.0
57.0
............ 87.5
surface mean bubble size calculated from the distributions; the mean size
becomes constant withz z 3OT.
Clift et al. (C6), in their study of bubbles in fluidized beds, indicate that
instead of having a discrete “maximum stable bubble size” we can expect
bubble splitting to occur over a relatively broad and continuous range of
bubble sizes. Whether or not a particular bubble splits will depend not
only on size but also on angular position, wavelength, and amplitude of
disturbances of the bubble interface. It seems likely that measured maxi-
mum stable bubble diameters correspond to mean diameters for systems
in which dynamic equilibrium has been achieved between coalescence
and splitting.
There are three main mechanisms which determine the bubble sizes:
coalescence of bubbles, splitting of a bubble under a given disturbance,
and the occurrence of disturbances in the bubble-flow equipment. The
latter two will be discussed in what follows. (Coaledcence is discussed in
C7, M10, T16, and W10).
1. Splitting of a Bubble
In 1950 G. I. Taylor (T8) investigated whether a small disturbance
superimposed on the interface between two immiscible fluids of different
densities was amplified or damped away when the fluids were accelerated
at a constant rate. Four years later Bellman and Pennington (B 10) devel-
oped this concept further to include the influence of liquid viscosity and
surface tension. Rice and Wilhelm (RS) applied Taylor instability to the
mechanism of bubble formation in a fluidized bed. Recently, Clift and
Grace (C7) have shown that the first stage of bubble splitting in a fluidized
bed is to develop an indentation in the bubble roof, which then moves
around the bubble periphery while growing to form a curtain of particles.
Splitting occurs if the lower edge of the curtain reaches the base of the
bubble before the top passes the equator. It was suggested that this phe-
nomenon results from instability of the type discussed by Taylor, where a
heavy fluid overlies a lighter one.
For a gas bubble in liquid, Henriksen and Ostergard (H7) have observed
that the bubble is broken up by disturbances created by the downward jet
of liquid. According to their observation, a finger of liquid projects down
from the roof and eventually divides the bubble in two. The reason that
bubble breakup in methanol is easier than in water is due to lower surface
tension. Based on these observations, they support the hypothesis of Clift
and Grace that the bubble is broken up as a result of the Taylor instability.
Linearized stability analysis has been performed (B 10) by superposing a
sinusoidal disturbance of wavelength A (which equals 27r/k, with k the
FLUIDIZED CATALYST BEDS 353
WbMe phase :
density = 1.26 Xl@g/cm’
80
Continuous phase:
:U=O dyne/cm
: U=GOdyne/crr
: kinematic viscosit
0
wvelength : 1 = 2 it/ k [cm]
reduction in the growth factor and an increase in the most sensitive wave-
length, at which n is a maximum. For the gas-fluidized bed, uM is the
apparent kinematic viscosity of the particulate phase defined ((26) by uM=
pp/ppese,where pup is the effective Newtonian shear viscosity of particu-
late fluid, pp the density of single particle, and eSethe volume fraction of
particles in the emulsion phase. Clift et al. (C6) have shown that the
interstitial gas velocity has virtually no effect on the stability of a bubble
roof. In their analysis the velocity is varied in the range of 1-100 cm/sec.
When the surface tension is finite, the disturbance tends to be sup-
pressed by the surface tension, the influence of which is noticeable due to
smaller curvature of the interface as the wavelength decreases. Hence
there is a lower limit for wavelength, below which the growth factor takes
negative values. The wavelength Amin where the instability is neutral is
(H7):
Amin = 2da/ptg)”2 (5-18)
For an air-water system at room temperature Amin equals 1.70 cm, and for
air-methanol Amin is 1.04 cm. The influence of surface tension on growth
factor is shown in Fig. 47, where the behavior discussed previously is
observed clearly.
As noted previously, disturbances initiated on the roof of a bubble are
swept around the periphery, so that in practice a bubble does not split
unless the disturbance has grown sufficiently before the tip of the growing
spike reaches the side of the bubble. Clift et al. (C6) have estimated the
likelihood of splitting by comparing the time required for a disturbance to
grow by a given factor with the time available for growth.
The required growth time may be provided by T, = l / n , where T, is the
time required for a small-amplitude disturbance to grow by a factor e.
Although the T, approach is no longer accurate for disturbances which
have grown beyond the scale described by the linearizeed analysis, the
estimate is retained in the absence of a better alternative.
An estimate of the time available for growth, T,, may be obtained from
the dense-phase tangential velocity r/$ (equal to r b &/do, with $ the
angle measured from the vertical) at the bubble-emulsion interface:
(5-19)
the disturbance originates, and case B has a node located A/4 from the
bubble nose (so that the nose is an antinode in the initial disturbance).
With these the maximum time available for growth, T ~ , , is calculated from
Eq. (5-19). A bubble is liable to be split by a disturbance for which T,, >
T,; here, only disturbances with wavelength less than the arc length from
the nose to the equator, A IA,,, = m& are considered, since a distur-
bance with wavelength greater than,,,A represents a gross deformation of
the bubble rather than a perturbation on the interface.
Figure 48 shows the calculation for a bubble with a radius of 2 or 5 cm in
a gas-fluidized bed. In the bed with uM = 4 cmz/sec, the bubble withr, =
2 cm shows T ~ > , T, for A in the range 0.35-3. I cm, so that a disturbance in
this range may cause splitting. However, the same bubble in the bed with
uM = 10 cm2/sec always shows T,, < T,, and therefore should not split
whatever the wavelength of the disturbance. By contrast, the bubble of
5-cm radius in either bed is liable to be split by a broad range of A. The
most sensitive wavelength (minimum T,; point m in Fig. 48) does not
correspond to the A most likely to cause splitting; in some cases,,A is
shorter than the most sensitive wavelength. When the most sensitive
wavelength is within the range of possible disturbance, the ratio T,,/T~ is a
maximum for a wavelength less than the most sensitive.
An important conclusion of the above approximate Analysis is that the
effective kinematic viscosity of the emulsion is the dominant factor de-
termining both the initial growth of instabilities and the most sensitive
wave number. Thus, prediction of the effect of system properties on bub-
ble stability depends on prediction of the effect on the effective kinematic
FIG.48. Comparison of maximum time available for growth,ram,with time for a distur-
bance to grow by a factor e , re, for a gas-fluidized bed; disturbance is case B (after C6).
356 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
viscosity. Based on this concept, CliR et al. (C6) have suggested that
bubbles split more readly in a small-particle system and that fluidization is
more aggregative as the ratio pp/pg is increased. As for particle size dis-
tribution, Tsutsui (T27) has shown that the wide size distribution of a
standard FCC catalyst gives better fluidization than the FCC catalyst
having a very narrow size cut with the same equivalent mean diameter.
Geldart (G4, G5) indicates that fine particles of Group A (see Fig. 3) show
good fluidity. Better fluidity of given particles is primarily determined by
the volume-surface mean diameter of the particles and is not much influ-
enced by their size distribution, according to Geldart.
Basic understanding of the effective kinematic viscosity of the emulsion
is potentially important from the viewpoint of planning and operating
industrial fluidized beds. There are mechanisms which lead to the forma-
tion of very small bubbles; these are omitted here since the performance
of larger bubbles usually dominates the rate processes under high rates of
aeration.
FIG. 49. Splitting of bubble by fluid-dynamic disturbances: (a) from Henricksen and
Dstergaard (H7); (b) and (c) from T. Kikuchi. (b) Air-water system, DT = 8 cm, U , = 2.15
c d s e c , continuous bubbling. (c) Air-aqueous glycerol system, D, = 8 cm, CJ, = 4 cdsec,
continuous bubbling.
358 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
a certain high level; this is clearly observed in plate C of Fig. 49. Plate B in
the figure shows a large bubble about to be split by the turbulence induced
by continuous bubbling.
Figure 50 shows the influence of molecular kinematic viscosity v M of
liquid on the mean bubble size, when a 15-cm-diam. bubble column is
operated at the superficial air velocity of U, = 6.7-12.7 cm/sec (see data
in K10). The mean bubble size is seen to increase only slightly as uM
increases from 0.04 to 0.5-0.7 cm2/sec (aqueous glycerol solution); the
flow properties here are well into the recirculation-flow regime. As vM
increases beyond 0.7 cm2/sec the mean bubble size starts to increase
rapidly, transferring finally to the high-viscosity branch PQ of insufficient
breakup due to disturbances of low level. A similar phenomenon has also
been observed for a larger-diameter bubble column by Ueyama and
Miyauchi (U4),where a different viscosity is obtained for the flow transi-
tion effective to bubble splitting.
The mean bubble size that concerns us here is on the order of 5 cm, so
that the Eotvos number Eo (equal to dggpl/c+) is well over 40 for usual
bubble-column liquids. The bubbles are of spherical-cap type under this
condition, which is essentially equivalent to a Weber-number criterion We
(equal to dbpliii/(T) > 20, since U, = (H4, H5). The bubbles in a
fluidized catalyst bed satisfy the above criterion, since CT -
0. Conse-
quently, surface tension has relatively little effect, and instead the splitting
is closely related to disturbances induced by the bulk turbulence, the
intensity and the scale of which are mainly governed by the fluidity of the
continuous phase and the operating gas velocity.
Figure 51 shows the observed bubble diameter db in fluidized beds as a
30
20 -
___--_-
Column ----
diama
T
a ?
Id 5
0 UG: 6.7 - 12.7 cm/sec
A UG= 8.9 cm/ sec
dD1 2 34570.1 2 3 4 5 7 1 2 3 4 5 7 10 2 3
viscosity [poise]
FIG. 50. Change of mean bubble size with liquid viscosity under constant aeration.
Aqueous glycerol solution, D , = 15 cm, U c = 6.7-12.7 c d s e c . Data by Kimura (K10).
FLUIDIZED CATALYST BEDS 359
2-
FIG.51. Observed bubble diameter d b versus particle terminal velocity u,; data from
Horio and Wen (H21). Mean bubble size by Morooka et al. (M43) and that by Yamazaki
et al. (Y6)are added.
function of the particle terminal velocity ut. The data are mostly taken
from Horio and Wen (H21). Values of mean steady bubble sizes under
turbulent aeration are from Morooka et al. (M43) and Yamazaki et al.
(Y6). As the size of particles decreases while good fluidity is retained,
db tends to approach the steady bubble diameter d b s observed for bubble
columns of low-viscosity liquids (Fig. 34). In this domain wheredb = d b ,
dynamic equilibrium exists between bubbles of various sizes as a result
of turbulent motion prevailing in the emulsion; there is a steady mean
bubble size here rather than the maximum stable bubble size.
The maximum stable bubble diameter dbms defined by Davidson and
Harrison (D3) seems to apply for an emulsion of large particles, where the
turbulent motion is weak. Bubbles tend to coalesce with each other to
grow ultimately to dbms. The observed bubble diameter db is closely re-
lated to bubble splitting and coalescence as a result of turbulence. The
data in Fig. 51 are expressed by an approximation that db is the sum of dbs
and dbmsor
db = dbs + dbms
A. BUBBLEDYNAMICS
Although gas bubbles ascending through the emulsion of fine catalyst
particles are constantly splitting and coalescing (Sections 11, 111, and V),
they are largely free of particles (H14, K13,T19). Such a bubble, may be
pictured as essentially spherical, with the lower of its volume occupied
by particles (Fig. 52). As the bubble rises, it displaces solid particles
around it and cames some particles upward in its wake. Its rise velocity ii,
is proportional to the square root of its frontal diameter. If ii, is lower than
the velocity of the interstitial gas umf= Umf/cmf, the gas from the emulsion
enters the lower part of the bubble and leaves through the roof; however,
for small particles, 1, is usually larger than umf,and the circulation pattern
as shown in Fig. 52 is encountered. Gas circulates up the central core of
the bubble and, upon emerging from the roof, encounters drag from the
particles streaming past. Consequently this gas is swept downward rela-
tive to the bubble and is re-entrained at the bottom because of the prevail-
ing pressure gradient there.
It follows that as the bubble rises more quickly, the gas which emerges
from its roof will penetrate outward a smaller distance before being swept
downward; i.e., the limit of penetration will be nearer to the bubble wall.
The total volume enclosed by the limit of penetration is called the cloud,
FLUIDIZED CATALYST BEDS 36 1
yAxis 01 symmetry
FIG.52. (a) Gas and particle flow pattern near a bubble, a = 2.5 and (b) a photograph
showing cloud size re and bubble size rb by using NO,-containing-bubble technique,
a = 2.5, 230 pm ballotini (after R17).
and the region between the bubble-cavity wall and the limit of penetration
is called the cloud overlap. This region and part of the wake are the only
regions where the gas in the bubble can contact particles.
These physical pictures of bubble dynamics have been developed by
Davidson and Harrison (D3), Murray (M46, M47), Pyle and Rose (B),
Rowe et al. (R17) and others. Jackson (Jl) considered that a mantle
of bed with increased voidage exists near the roof of the bubble. Most of
the work on fluid-mechanical theories of aggregative fluidization have
been reviewed by Jackson (in DS).
The radius of the cloudr, has been given by Davidson and Harrison (D3)
for a = &/Urn, > I :
rc/rb = [(a + 2)/((r - i ) y 3 (6-1)
Alternatively, Murray (M47) proposes for (Y > 1:
(a - l)(rc/rb)4 - a(rc/rb)= cos 8 (6-2)
The experiments of Stewart (S20) show (Fig. 53) that the radius of an
observed bubble cloud at 8 = 0 falls between Eqs. (6-1) and (6-2). Figure
54 shows the thickness of the cloud-overlap region estimated by Eqs. (6-1)
and (6-2) for a fine microspherical catalyst fluidized by ambient air, where
the minimum fluidization velocity Urn,is calculated by Wen and Yu's
equation (W4). The cloud-overlap region is seen as being limited to at
most a few layers of fine particles (M30). Rieke and Pigford (R8) observed
experimentally for a two-dimensional fluidized bed that the gas in the
txprrin
El Jackror
FIG. 53. Ratio of cloud radius rc to bubble radius rb on vertical axis above three-
dimensional bubbles. Comparison of experimental values with various theoretical predic-
tions. The curve “Jackson (modified)” refers to resultsobtained by Jackson’s type of analy-
sis, but with the Davies-Taylor type rise velocity k v g r , , where k is chosen to give the best
fit to the experimental rise velocity (from S20 and D5).
2 3 4 5 678910
db (Em)
FIG.54. Thickness of gas cloud. pp = 1.0 g/cm3, emf= 0.5, fluidized by ambient air.
F L U I D I Z E D CATALYST BEDS 363
bubble wake does not return to the bubble cavity while circulating. This is
in conformity with the observation by Rowe et al. (R17); see Fig. 52. As a
consequence, the observed gas-flow patterns suggest that there is no di-
rect gas communication between the wake and the cloud-overlap region.
Progress in understanding the bubble dynamics has made it possible to
formulate transport equations for heat and mass transfer through the bub-
ble wall. which we now consider.
B. HEATA N D MASSTRANSFER
THROUGH THE BUBBLE
WALL
Most of the experimental work so far has been concerned with the
performance of an isolated bubble; few publications are available on
fluidized catalyst beds under a high aeration rate. In 1961 van Deemter
(Vl) analyzed the transient and steady-state gas mixing data by Mason
(given in V1) taken in a 7.6-cm-diam. cracking catalyst bed. Van Deemter
based his calculation of the overall mass transfer coefficient kebab on the
two-phase theory of fluidization. Further, axial mixing of gas and solid
particles has been studied for fluidized catalyst beds by de Groot (D7),
Botton (B13), van Swaay and Zuiderweg (V8), and Morookaet al. (M42).
For crushed silica fluidized at ambient temperature and pressure in beds
with diameter 0.3-1.5 m, van Swaay and Zuiderweg (V8) have reported
the apparent HTU:
(Hob)app -
- a = (1.8 - *)(3.5
k uG
ob b OFZ5
- -)LPZ5
2.5 (6-3)
where (&,)app, DT, and Lf are in meters. This relation is expected to apply
for D, and Lf up to 10 m, with 15% fines content of the solid (d, < 44 pm).
The use of microspherical particles with a larger fraction of fines showed
somewhat better mass transfer; but the dependence of HTU or of k o b a b / &
on bed dimensions was the same, as shown in Fig. 55. The coefficients are
less dependent on gas velocity and bed height for particles with good
fluidity. The data by de Groot (D7) indicate a doubtful effect of bed height,
in contrast to Eq. (6-3).
Morookaet al. (M42) studied the effect of gas adsorption on mass trans-
fer between bubbles and emulsion phase, while measuring the residence
time distributions of tracer gases for helium, carbon dioxide, Freon 12,
and Freon 22 in 12- and 19-cm-diam. free and baffled beds of cracking
catalyst. The k O & values are calculated from the distributions for differ-
ent gases according to the two-phase model, providing for the adsorption
equilibrium.
The quantity kebab increases with increasing superficial gas velocity VG;
it shows a higher value when gas is more strongly adsorbed by the parti-
364 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
-
V
Q, - @ 2-2-4.9 H2 de Groot (D7) -
A 1-3-3 H2 Botton(B13)
-
v)
\
c
2-.
0 v 1.6-3.4 H2 Botton(B13)
I
I I 1 I I l 1 1 l I I I 1 I l l 1 1 I I I l l
2 4 6 8 2 4 6 8 2 4 6 8
1 10 1 o2
m, C - 3
I 1 I
20'C'\
$5
4
\
I
i8.I
10 15 20
Fraction of fines ( < 4 4 p m ) '10
FIG.57. Effect of fines on kebab in fluid bed (data taken from D8).
cles (rn larger). The effect of increasing gas velocity or adding horizontal
baffles is mainly to increase the mean bubble holdup in the bed.
The group k & h / E h m is plotted in Fig. 56 as a function of the gas
partition ratio m,. Obviously k,,h increases with m,,owing to decreased
particle-side mass-transfer resistance. Data of Yamazaki and Miyauchi
(Y4)are plotted in Fig. 56 and match well as an extension of the Morooka
data. Using a quick-response thermoelement, Yamazaki measured the
temperature distributions of the bubble and emulsion phases, from which
the overall heat-transfer coefficient is calculated for 7.9- and 19-cm-diam.
cracking catalyst beds. Since the heat capacity of the emulsion phase is
much higher than that of the bubble phase, the overall coefficient is nearly
equal to the film coefficient h h a h for the bubble side. The heat-transfer
coefficient reduces to the mass-transfer coefficient according to the equal-
ity h h a h / C p g & = k b a h .
The effect of added fines content in increasing k&h, measured by de
Vrieset al. (D8), is shown in Fig. 57. Other work on heat and mass transfer
for teeter-type beds is reviewed by Kunii and Levenspiel (K24) and in
Davidson and Harrison's book (D5).
C. MASSTRANSFER
BETWEEN BUBBLE
A N D EMULSION
PHASE
Gas exchange during catalytic reaction between bubbles and the emul-
sion of catalyst particles is essentially a diffusional phenomenon, with
simultaneous adsorption and reaction in the emulsion. Based on this con-
cept, the overall mass transfer has been analyzed by Miyauchi and
366 M I Y A U C H I , F U R U S A K I , MOROOKA, A N D I K E D A
fqt m&c,
=bob cb
CC
4 + C.0
c3 diffusion
FIG.58. Schematic diagram of mass transfer near a bubble in fluid bed (M30).
F L U I D I Z E D CATALYST BEDS 367
for e > 0,
x= -ma . ch = ChO
x = +m: c, = ceo
The initial concentration of the cloud-overlap region c,, is obtained from
the reactant gas balance taken at the roof of the bubble, where the catalyst
particles are fed to the cloud-overlap region from the emulsion above at
volumetric flow rate q s :
CbOqf + 4qsceO = (qf +q d c C O (6-8)
When a >> 1 , the slip velocity of particles relative to the gas is negligible
in comparison with the particle velocity, so that qf/qe is nearly equal to
q s / E s , . With this equality Eq. (6-8) is modified to the form:
CcO = (EfecbO + %%eceO)/m (6-9)
The concentration distributions in each phase and the mass flux across
the bubble surface have been solved to satisfy Eqs. (6-4)-(6-6) under the
restrictions of Eqs. (6-7) and (6-9) (M26, M30). The overall mass-transfer
coefficient between the bubble cavity and the emulsion phase kob has been
given as follows:
(6- 10)
368 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
The film mass-transfer coefficients for the bubble and the emulsion phase
( k b and k,, respectively) are given by
kb = (2/7T"2)(D~/Pb)112 (6-12)
ke = ( 2 / T '/')( mDeff/Tb)' (6-13)
where Tb is the time necessary for local surface renewal; p is the Hatta
number (or reaction factor) for unsteady gas absorption with a first-order
irreversible reaction, as given by Danckwerts (Dl):
J =zerf(
4mH
$) - l*& exp[ - (4mrr:, x +E)]dx
X
(6-15)
(6-16)
083
- 1.0
-
+ 1.3
-
1.5
q
IS
L 20
2 40
35
Pe
FIG.59. Integral J as a function of Pe and m H (M26).
FLUIDIZED CATALYST BEDS 369
The correction factor r ) to Eq. (6-22) is nearly equal to unity for (Y > 10, so
that Eq. (6-22) is almost equivalent to Eq. (6-23) for smaller particles.
Equation (6-21) is independent of mH or rate of chemical reaction.
Drinkenburg and Rietema (D 17) have presented a numerical computa-
tion of kob based on the stream functions given by Davidson and Harrison
(D3) and by Murray (M47). The bubble-void resistance to mass transfer
has been neglected. Enhancement of gas transfer rate by diffusion with
simultaneous chemical reaction (Fig. 5 of D17) is reasonably well ex-
pressed by Eq. (6-11), the enhancement being expressed as the Hatta
number. Enhancement by physical adsorption (Fig. 2 of D17) is also ap-
proximated by Eqs. (6-22) or (6-23) for smaller particles.
When the reaction data by Lewis et al. (L12) are analyzed by utilizing
Eq. (6-lo), it turns out that PI.amounts to 1.45 (free-flow bubbles) for the
highest reaction rate constant. Hence it seems reasonable to take into
consideration the influence of pr on kob. Also, it has been shown expen-
mentally that k b is not necessarily negligible for fluid beds (Y4).
If the adsorption equilibrium is not attained instantaneously, a different
analysis is needed. Toei et al. (T18) studied the mechanism of heat and
mass transfer between bubbles and emulsion phase under such circum-
stances. The dependence of diffusion rate on bulk flow across the bubble
interface also becomes important when coarse particles are fluidized
(H16). For two-dimensional bubbles Chavarie and Grace (C7a) compared
various interphase mass-transfer models.
F L U I D I Z E D CATALYST BEDS 37 I
D. PARTICLE
CAPACITANCE
EFFECT
As is well known, the temperature in a fluidized catalyst bed is nearly
uniform even when considerable heat is generated in the bed (G8, L11).
This comes from the bed’s large content of solid particles, which disperse
the heat axially as a result of particle heat capacity. Similarly, the gas
concentration in the bed becomes uniform if the partition of transferring
gas component is favorable to the particle phase. The circulating particles,
holding a large amount of transferring component, iocally exchange the
component with the surrounding gas phase; the larger the partition ratio is
of the component, the higher the rate of concentration equalization in the
gas phase of the bed. Enhancement of heat and mass transport in the bed
372 M I Y A U C H I , F U R U S A K I , MOROOKA, A N D I K E D A
The particle capacitance effect in bubble mass transfer was shown ear-
lier (see Section V1,C) for the streaming emulsion outside a bubble, so
that rn is included in the emulsion-side film coefficient k,. Chiba and
Kobayashi (C3) and Drinkenburg and Rietema (D17) introduce the same
effect. Yokota et al. (Y13), in studies of mass transfer from submerged
surfaces in a gas fludized bed, find that the rate is strongly enhanced by
adsorption of the transferring gas component on fluidized particles. The
role of the particle capacitance effect in heat transfer has been discussed
by Mickley and Fairbanks (M14), Mickley et al. (M16), Drinkenburg
(D16), Yoshida et al. (Y18), and Kunii and Levenspiel (K24).
E. AXIALA N D RADIAL
MIXING
OF HEATA N D MASS
FIG. 60. Longitudinal dispersion of gas in fluid beds (after Miyauchi). The full curve
is calculated from Eqs. (6-3 I ) and (6-34) for He-air system.
I- / /
0
.0l
V
- 8
0 5 10 15 20 25 30
% FINES
Uf Ccmlsec3
particles with bed gas has been neglected, since the influence is usually
small for transient response of tracer gas at a moderate partition ratio.
The simplifications taken here are, first, that the catalyst particles are
usually several tens of pm in mean diameter, so that the interstitial gas
velocity through the emulsion phase is about umf(equal to Umf/Emf) and is
neglected in comparison with UGand emulsion-phase recirculation veloc-
ity. Also, the particles and the gas in the emulsion are in local adsorption
equilibrium the partition ratio m (equal to efe + mseSe)applying for unit
volume of the emulsion phase. Second, the central upflow and the periph-
eral downflow recirculate in the bed as ideal plug flows at approximately
the same linear velocities. Furthermore, the bubble phase ascends only
through the central upflow emulsion.
With these simplifying assumptions, the equations of continuity are
given for transient response of a tracer gas as follows (M38):
Eb &,/a6 = -uc ach/az - kohah(Cb - cfeu) (6-26a)
for the bubble phase;
mieu acfeu/a8= -mu, acfe,/dz + kohab(Ch - cfeu)
-mkexaex(cfeu - CfecJ (6-26b)
FLUIDIZED CATALYST BEDS 377
The variance of the residence times of the tracer gas with respect to
mean residence time is obtained from Eqs. (6-26) and (6-27), according
to van der Laan (V3), as follows:
0 2 - 1 mie 1
2- (ib+ mi,) + N,,
-(Nv 2 ( i b + mi,)
(6-28)
+
p = [I 4NeX/(I N~~)N,,b]’’*
Here, we note that the influence of particle motion on an impulse response
has been treated by van Deemter (V2), Yoshida and Kunii (Y15), Fryer
and Potter (F10, F1 I), and Nguyen and Fryer (N2), but their models are
all different from the one treated above.
Morookaet al. (M42) further solved the impulse response of a tracer gas
for a fluidized catalyst bed according to the one-dimensional two-phase
diffusion model (Vl), where the influence of the particle capacitance effect
was considered under the assumption of local adsorption equilibrium. The
equations of continuity for the tracer gas are:
Eb acb/ae = - u, acb/az - kobab(Cb - ce) (6-29a)
378 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
(6-3 1)
where Nob = kobabLf/UG, (PeB), = U,L,/mZ,E,, and
= [l + 4(PeB),/Nob]1’2, A = 4No,/2
For the present simplified recirculation model, the axial dispersion
coefficient is given by the following approximate equation for arbitrary m :
( m Ue)2/mkexaex
= Ee (6-32)
This relation is derived by considering a steady-state tracer exchange
between the ascending and descending emulsions and the bubbles, and
comparing the resulting equation with that of the diffusion model. Since
E, = mEeEZs,the above equation is equivalent to the following relation
from van Deemter (V2):
Glkexaex = EeEs (6-33)
Once the variance of the residence time distribution has been obtained,
the axial dispersion coefficient of the tracer gas E , is calculated by combin-
ing the following equation for the one-dimensional diffusion model (V3)
with Eq. (6-28) or Eq. (6-31):
F. REVIEWO F WALLHEATTRANSFER
I N FLUID
BEDS
Numerous measurements of heat transfer between the fluid bed and
immersed surface have been carried out; these have been reviewed by
Kunii and Levenspiel (K24), Davidson and Harrison (D5), Botterill (B12),
and Zabrodsky ef al. (Zl, 22). Wen and Leva (W3) and Wender and
Cooper (W5) correlated heat-transfer data at the walls of fluidized beds.
However, the data utilized for their correlation were obtained with a
380 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
fluidized bed less than 10.2 cm in diameter, and measurements under fluid
bed conditions were not included. Heat transfer between fluidized beds
and immersed vertical tubes was correlated by Wender and Cooper (WS).
This correlation includes data with fine porous catalyst particles by Olin
and Dean (04) and Mickley and Fairbanks (M14):
q’0L I I I I I I I I I I I 1
-
I l l l L
u
- 8-
-
yc-=---
p-
E6- ?bubble column -
- bubble column --------- -
-
- 0 dp;9pm -
/-<:- -
.I.--
--
\ -
c
f2- fluid bed
a
104 -
dp = 3 8 , 4 8 p m @ 140
I I 1 I I 1 1 1 1 1 I I I 1 1 1 1
1
renewal on the vertical heat-transfer surface in fluid beds agrees with the
data in bubble columns. The larger the particle size that is used in
fluidized beds, the smaller the rate of the packet renewal becomes com-
pared with that in bubble columns. Vreedenburg (V13) correlated data on
heat transfer between horizontal tubes and the fluidized bed, including
data for fluid beds; the data are widely scattered, and further study ap-
pears to be needed.
The bed-to-wall heat-transfer coefficient in a fluidized bed at high tem-
peratures is larger than at room temperature (Y19). Questions have been
raised about the effect of radiative heat transfer at high temperatures
(B12, Y19), and more studies are necessary on this problem.
In the transition zone and the freeboard region, heat transfer between
bed and wall is a function of bed density. Shirai et al. (S9, S11) studied
heat transfer from a sphere immersed in the fluidized bed and showed the
trend of decreasing heat-transfer coefficient with decreasing bed density.
A. EXPERIMENTAL
FACTSA N D REACTORMODELS
Early attempts to approximate gas-solid contacting in fluid catalyst
beds were based on the assumption either of isothermal plug flow of the
fluidizing gas through the bed with the catalyst uniformly distributed or of
isothermal complete mixing of the gas within the bed. The simple disper-
sion model, falling between the above two cases, was also used (G8, R4).
Evidence from both large-scale and laboratory observations (G9a, L12),
382 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
When the bed is not fixed, but is fluidized freely at the same gas velocity
as above, with UG much greater than the minimum fluidization velocity
(Imf,the bed becomes aggregative and deviates from the piston-flow reac-
tor. As a consequence, the apparent overall rate constant of reaction KO,
for the fluidized bed is smaller than k,, and the extent of reaction is
expressed by
(CO/Co)fluid bed = exp(-KoRLq/ uG) (7-3)
K O , is the same physical quantity as the specific converting power Q of
Lewis et al. (L12).
Usually, the fluidized catalyst bed constitutes an aggregative dense
phase forz between 0 and Lf, and a dilute phase forz > Lf. Hence, another
apparent overall rate constant k,, is defined on the basis of Lf. Taking Lf
instead of L,, Eq. (7-3) is rewritten:
(CdCo)fl"id bed = exp(-koRLf/ UG) (7-4)
FLUIDIZED CATALYST BEDS 383
I
- - Phase boundary
I
Err ---
I, --- Dense phase , L,
fw KJ ---
&ob ab,____ - ---Emulsion phase
k b ab. (C.,E 1
k e ab,
- -- Distributor
T Gas, UG(=EG)
FIG.65. Schematic expression of general two-phase dense bed model and the parameters
included (homogeneous return flow is by F10-12, K24, L.5).
with bubble gas v, the cloud-wake fraction fwzb, and the axial dispersion
coefficient of the emulsion Ezs.
In the case of fluid catalyst beds, the mean bubble size stays approxi-
mately constant axially, so that the above parameters are treated as re-
maining unchanged axially. Also, the interstitial gas flow through the
emulsion is neglected in comparison with the feed gas flow U G ,so as to
considerably simplify the treatment. Usually the bubble phase is assumed
to be unmixed axially.
Under the above simplifications, the two-phase consecutive model by
Shen and Johnston (S7)with vertically unmixed emulsion is written for a
differential bed height:
-uG dcb/dz - kobab(Cb - c,) =0 (7-6a)
FLUIDIZED CATALYST BEDS 385
with
kobab(Cb - ce) - gekrc, = 0 (7-6b)
where the original capacity coefficient is expressed as the product of the
overall coefficient k o b and the specific surface area of the bubbles a b in-
cluded in the bubble-emulsion mixed phases. A first-order irreversible
reaction takes place in the emulsion. Solving the equations, the overall
rate constant as defined by Eq. (7-4) is given by
koR = [(kebab)-' + (Gek)-'l-' (7-7)
When a part of the emulsion volume v is included in bubbles and con-
tacts freely with the bubble gas, the effective fraction of the emulsion
phase is Ze - Y , and Eq. (7-6) is rewritten:
-UG d C b / d Z - kobab(Cb - Ce) - VkrCb = 0 (7-8a)
kobab(Cb - ce) - (ge - v)krce = 0 (7-8b)
This case is the direct-contact model corresponding to the vertically un-
mixed emulsion (VUE) by Lewis el al. (L12), a special case of the model
by Mathis and Watson (M8). Solving the equations, the overall rate con-
stant is
{(16bab)-' + [(G - v)k]-'}-' (7-9)
The bubble flow model by Orcutt et al. (C9, 06) uses the film coefficient
kbUb between the bubble gas and the emulsion. In their vertically unmixed
emulsion k,, is given, with negligible gas flow through the emulsion, by
Eq. (7-7) with k b a b in place of kebab. When the emulsion is mixed per-
fectly, the extent of reaction is given by
c,/co = e-Nb + (1 - e+")/[l + EeNr/(l- e-"")] (7-10)
where N b = k b a b L f / u G and Nr = kLf/UG. Hence, k o R is obtained by
equating Eqs. (7-10) and (7-4).
In the bubbling bed model of Kunii and Levenspiel (K24), there are two
transfer steps for the bubble mass transfer, namely, the transfer between
bubble void and cloud-particle overlap region k b a b and that between the
cloud-particle overlap region and the emulsion phase k e a b . They further
assume that the cloud-particle overlap region and the bubble wake are
mixed perfectly, and contact freely with the cloud gas. Their basic equa-
tions in the present notation are (for their case 2):
-uG dcb/dz - kbab(Cb - c,) = 0 (7-lla)
kbab(Cb - c,) - fwgbkrc, - keab(C, - ce) = 0 (7-llb)
keab(cc - ce) - (ge - f w g b ) k r C e = (7- 1Ic)
386 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
with
kor = {(keab)-' + [(ze - fwEb)krI-'I-' + fwihkr (7-13)
Here k b is given by Davidson and Harrison (D3), and ke by the Higbie
penetration mechanism (see Section VI where k h c and k,, are utilized in
place of k b and k,, respectively).
It is interesting to observe the limiting cases where one particular
mechanism dominates the whole reactor performance. One extreme is the
case where k, << k o b U b , leading to (for all models):
koR = E,k, for k, << k o b U b (7-14)
Obviously, model identification is not possible under this condition. An-
other extreme is k , much larger than the bubble mass-transfer coefficient.
The extremes differ from one another, as follows:
koR = kebab (7-15)
for the two-phase consecutive model, Eq. (7-7);
koR = kh'b (7-16)
for the bubble flow model, Eq. (7-lo), and bubbling bed model, Eq. (7-13);
koR = vk, (7- 17)
for the direct contact model (VUE), Eq. (7-9). Here we note that koR for
the models not including directly contacting catalyst are upper-end limited
by the mass-transfer coefficient kobUb or k b U b . Only the direct-contact
model shows koRincreasing linearly with k, under a given flow condition.
These tendencies are more clearly visualized from Fig. 66, where koR is
shown as a function of k , at constant gas velocity. Numerical values given
to various parameters are those for an ethylene hydrogenation run by
Lewis et ul. (L12).
The models show considerably different k o R when k, 2 . 1 sec-'. The
bubble flow model (BFM) is for perfectly mixed emulsion, but a k,, only
slightly larger is obtained for BFM when the emulsion is vertically un-
mixed. For the bubbling bed model (BBM) koR is fairly sensitive to differ-
ent assumptions forf,; anfw of about unity is recommended ( F 1 1 , K24) to
fit the reaction data available. The value fw = 0.35 is taken from Rowe
(D5) for 75-prn-diam. spherical particles.
FLUIDIZED CATALYST BEDS 387
FIG.66. Overall reaction rate constant koR as a function of reaction rate constant k , for
various reactor models; = 25 c d s e c , ih
= 0.27, k&ab = 0.37 sec-I, k @ b = 0.94 sec-I,
k e a b = 0.61 sec-I, Y = O.lO,f, = 1.0 and 0.35; data are for a fluidized catalyst bed.
R Jsec-'J
FIG.67. Overall rate constant k,, for experimental ethylene hydrogenation runs by Ni-
or Cu-impregnated cracking catalyst beds. Dotted curves were calculated by the two-phase
consecutive model (TCM) and by the bubbling bed model (BBM). Full curves were calcu-
lated by the successive contact model.
account for the data, kbab or k&ab should be increased with increasing k;
this would imply that the bubble size decreases with increasing catalyst
activity at otherwise identical fluidization conditions, whereas the bubble
size should remain constant. A similar tendency is also observed for
catalytic ozone decomposition by a fluidized catalyst bed (see Fig. 68),
although not so decisively.
The direct contact model has some difficulties, however. In fluidized
beds, gas bubbles of very low solid content are usually considered to exist
in the dense phase (H14, K13, T19). Also, the cloud layer is negligibly
thin, due to small Umffor the usual fluid catalyst beds, according to equa-
tions of Davidson and Harrison (D3) and Murray (M47). The streamlines
of gas phase through a bubble have been observed to pass through the
cloud, but not through the bubble wake (R17). Thus there seems little
possibility of believing that the bubble gas is in direct contact with a sub-
stantial amount of catalyst in the bubble phase (see also Section V1,A).
Furthermore, the direct contact model is applied to the data by Gilliland
and Knudsen, and u in Eq. (7-9) is calculated to fit the data. Calculation
(M26) shows that the volume of catalyst, with an apparent density the
same as for the emulsion, which contacts the bubble gas freely exceeds
the volume of bubble gas itself ( u / c b = 3.3, 2.0, and 1.5, respectively, for
U G = 10, 20, and 30 cm/sec). This seems to be unsound physically.
FLUIDIZED CATALYST BEDS 389
rQ r [s~c.']
FIG.68. Summary of catalytic reaction data by fluid beds for ethylene hydrogenation and
ozone decomposition: curves were calculated by the successive contact model.
catalyst somewhere outside the dense phase. Based on the axial distribu-
tion of bed density, Miyauchi (M25) located the catalyst in the dilute
phase and developed the concept of the successive-contact mechanism. A
part of the catalyst may be placed in the jetting zone above the gas dis-
tributor, as explained in Section VIII. Catalyst erosion could be excessive
with too much jetting, however. Also, such a large amount of directly
contacting catalyst as observed for the Gilliland and Knudsen data is
difficult to locate solely in the jetting zone.
Perhaps the influence of dilute phase on the progress of reaction must
have appeared in the minds of many investigators, but it was never formu-
lated. It was stated (M33) that “Visually the particles in the dilute phase
are dispersed rather uniformly, so that under these flow conditions this
portion of the particles may take part in increasing the catalytic conver-
sion.” Interesting discussions by Riley (R12) are another example.
4. Other Reactions
Fluidized catalytic reactions have been industrially operated in the
“fluid bed” conditions, but most of the research has been carried out for
the “teeter bed.” Several studies of fluidized catalytic reaction are listed
in Table VI, which are of interest in considering transport phenomena in
fluidized catalyst beds.
B. THESUCCESSIVE
CONTACT
MECHANISM
A reactor model is developed to include reaction taking place in the
dilute phase, and to be reasonably consistent with the known flow proper-
ties of fluidized catalyst beds operated under relatively high gas velocity.
According to this model, reaction proceeds successively in the dense
phase and in the dilute phase.
4 4 UG u
rn1 Reaction rate Refer-
Author Reaction Luml [cml [cdsec] [cdsec] constant [sec-'1 ences
With (c z 10, the following equations of continuity are obtained for the
dense phase:
- UGd C b / d Z - kobab(Cb - Cf,,) - VkrCb = 0 (7-19)
kobab(Cb - Cfeu) - m k e x a e x ( c f e u - c f e d ) - ( i c e - v)krcfeu =O (7-20)
mkexaex(cfeu - cfed) - kgekrcfe, =0 (7-21)
where half of the emulsion phase is in ascending flow and the other half is
in descending flow. Equation (7-19) is for the bubble phase; Eqs. (7-20)
and (7-21) are for the ascending and descending emulsion phase, respec-
tively. The term (iGe - v) is the fraction of the ascending emulsion and 4Ge is
that of the descending emulsion.
When the exchange rate of emulsion between the upflow and the
downflow is high enough, the emulsion gas concentrations ce, and c e d
become nearly equal, so that Eqs. (7-20) and (7-21) are combined and
reduced to Lewis' VUE model, Eq. (7-8.b).
Equations (7-19)-(7-21) are combined by eliminating ce, and Ced, yield-
ing
-UG dCb/dZ - (16, + V4)Cb = 0 (7-22)
where l/kor = l/k,b& + 1/tGekr, and 6 is given in Eq. (7-28). For large-
diameter beds, cp could be smaller than 10 (see M29).
The bubble gas, with initial concentration co, leaves the dense phase at
z = 4 with concentration cLp and enters the dilute phase. For the dilute
phase the diffusion-model approach is more general, since a fair amount of
gas mixing is observed in a large-diameter bed (H19). Here, however, the
usual piston-flow assumption is made for simplicity, since in small-scale
beds the particles in the dilute phase are suspended rather uniformly.
Actually the observed flow properties of the phase, including the transi-
tion zone, seem to be more complicated.
The equation of continuity with the piston-flow assumption is
-UG dc/dZ - gedkrC = 0 (7-23)
where the physical meaning of Zed has been explained in Section VII,B, 1.
The reactant concentration coof the gas leaving the dilute phase at z =
& is obtained by integrating Eqs. (7-22) and (7-23) as follows:
co/c0 = exP{-[kor + (v + e)krlL,/UG) (7-24)
where
(7-25)
FLUIDIZED CATALYST BEDS 395
Accordingly, eL, is the total fraction of catalyst (with the same density as
the emulsion) suspended in the dilute phase. The total fraction of catalyst
in the dense phase is ZeL,, so that the following equality exists for the
successive contact model:
L, = ( e + Ze)Lf (7-26)
Comparing Eqs. (7-24) and (7-4), koR can be given by
koR = [(kebab)-' + (4Eekr)-']-' + (v + e)k, (7-27)
where 6, defined by Eq. (7-28) is the fraction of catalyst in the dense phase
effectively utilized for reaction:
4 =1 - (v/E,) + [2 + ( C A / m k e ~ d I - ' (7-28)
Equation (7-27) is reduced to the VUE model by Lewis et al., Eq. (7-9),
when e = 0 and Zek,/mkexaex<< 2; v is usually small in comparison with e
(Section VIII). In this case, with 4 = 1, Eq. (7-27) is reduced to the
simplest expression of the model:
koR = [(k,,a,)-' + (Zekr)-']-' + ek, (7-29)
3 . Experimental Evidence
Ethylene hydrogenation runs are summarized in Fig. 68, where the
data by Lewiset al. are taken from their smoothed curves (Figs. 7 and 8 of
L12). For these, approximate numerical values for cp and 4 are cp = 25-38
and4 = 0.95-0.73,for UG = lOcm/sec,andcp = 52-63and4 = 0.98-0.83
for UG = 40 cm/sec. Also, the simplest expression [Eq. (7-29)] is utilized
to calculate koR.The Lewis-Gilliland-Glass plot (see Section VII,C,2) has
been applied to obtain kebab and e from the experimental data; the other
data are similarly processed.
The data by Orcutt et al. (06) for catalytic ozone decomposition have
been processed in the same manner. Those by van Swaay and Zuiderweg
(V7, V8) for ozone decomposition are treated as follows. Their data are
presented in terms of an apparent HTU, (Hob)app, as defined by
(Hoblapp = u G / k ~ R - UG/Eekr (7-30)
where UG/ko,equals HoR,the overall height of a reaction unit by Hurt
(H18). Because (&,b)app equals UG/(kobab>app, they presume Eq. (7-7)
holds. Now koR is calculated using Eq. (7-30) by knowing the expenmen-
tally given (Hob)app; see Fig. 70.
In processing the above data, v is consistently assumed to be negligible
in comparison withe, and Eq. (7-29) is utilized in calculatingk,,, including
the full curves given in Fig. 68.
396 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
D~cm ut 1 db e 1 6b
200 - - 0 23 15 142 0.036 1 015
A 10 10 I 2.8 0.031 I 0.12
1 1 ' 'I I ' "'I
2o t
k, [sdl
FIG.70. Comparison of the calculated (Hob)app
with the data of van S w a y and Zuider-
weg (VS).
In the figure the data obtained by Gilliland and Knudsen from the use of
copper-impregnated cracking catalyst (quiescent bed density pmf = 0.33
g/cm3) show considerably higher koR than those by Lewis el al., who used
nickel-impregnated cracking catalyst (pmf = 0.52 g/cm3). Since the ex-
perimental runs have been performed under nearly isothermal conditions,
the higher koR by Gilliland seems mostly to come from their using less
dense catalyst than Lewis et al. used. In other words, at least qualita-
tively, a fair amount of catalyst is suspended in the dilute phase in con-
formity with higher e (or a ) .
Furusaki et al. use a standard cracking catalyst of wide size distribution
with smaller mean particle size, higher e, and higher UG than those of
Lewis et a f . (cf. Fig. 68).
The successive contact model seems to be sound in accounting for the
fluid-bed performance given in Fig. 68, since the model satisfies the flow
and transport properties of the beds as well. The model is also applied
successfully to catalytic oxidation of HCl (the Deacon reaction) in a fluid
bed (F14).
Recently, Yates and Rowe (YIO) have observed, on the basis of their
model for catalyst distribution in the freeboard region, that this region can
usually exert a considerable influence on the course of the reaction. Their
observation is essentially parallel with the concept of the successive con-
tact mechanism. However, they use the bubbling bed model in calculating
the reaction in the dense phase, so that the effect of directly contacting
catalyst seems to be corrected two times, first partially in the dense phase
and then in the freeboard region (see Section VII,A,3).
FLUIDIZED CATALYST BEDS 397
-
'I
T 5-
lo'
8 4-
I
3-
3= 16 cm/iec i--
Cb = 0.16 E. = 0.04
2- e = 0.074 , QOZO
d,= 5cm
10 r
7-
5-
4-
3-
FIG.71. Relative contribution of parameters to HoR = LIG/koR);P I ,P,, and P3 are the
respective locations where k,$,, e , and 0,start to be significant (M28).
398 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
Behie and Kehoe (B9) state that the effect of the jetting zone is salient,
especially for the case of very fast reaction. They use a very high mass-
transfer coefficient in their calculation (koJ = 57 cm/sec, uJ = 30 m/sec)
(B8). On the other hand, Eb could be very small in the jetting zone if the
dense phase holds only the jets of high speed and not bubbles, since Eb =
UJu,, = WJu, << 1 due to u, >> 0. Actually, the fraction of bubble phase
is high in this zone (B6, E8). In case of fluid beds, the spout above the
nozzle of the distributor is quite unsteady and splits into bubbles within a
very short distance. As for this, the observation by van Krevelen and
Hofiizer (V6) for a gas jet submerged in liquid is quite suggestive for
practical values of UG.
TABLE VII
REACTORMODELS"
KoR FOR VARIOUS
KoR= ( I - a ) / ( l / F + I / k J + ak,
koR/&R
Reactor model F U (=La/&)
~~
Direct contact
0 Data by Gilliland-Knudsen .
a1 UG= 30 cm/sec (G- 7 1
FIG.72. Overall rate constant KO, based on the quiescent bed height L , calculated by
various reactor models; U G= 30 cmlsec, Zb = 0.31, kebab = 0.17 sec-I, kbab= 0.42
sec-',k,ab = 0.27 sec-',a = 0 . 6 6 , ~= ~l,f, = 1.0 and 0.35.
sec-' for SCM and 0.057 sec-' for DCM; the latter value seems too small.
From Table VII one has
= [kobab/(&- Y)IDCM(VUE) = (kobab/ge)SCMb=O)
so that
(kobub)DCM(VUE)/(kobab)SCM(u=O) a (7-36)
= (ge - v)/ge = 1 -
In the present example a = 0.657 and 1 - a = & (see also discussion given
in Section VII,A,3).
CI '-100+140 MESH
Y
5 -
AV dp=122MICRONS
I
sr'2
Y
Y
8
1
INTERCEPT=F
1 I I I
Ob 5 10 15 20
kf (kr- K0R)
The concept of the successive contact mechanism has been given its
simplest form by dividing the fluidized catalyst bed into two parts-dense
phase and dilute phase. The concept has been found to apply to bed
performance, as shown in the preceding section. The reactor model has
been developed on the basis of several simplifying assumptions, partly to
retain mathematical simplicity as a workable design equation accounting
for the relative effects of the variables, and partly due to a relative lack of
information about bed performance. Further properties of the mechanism
are examined here, particularly as to axial distribution of reactivity inside
the bed.
A. AXIALDISTRIBUTION
OF REACTIVITY
IN A FLUIDBED
The dilute phase as defined in Section VII,B,l partly includes the transi-
tion zone. Gas-solid contacting there is much more complicated than is
FLUIDIZED CATALYST BEDS 403
0.2, is close to the observed profile. However, the axial reactivity dis-
tribution inside the bed is not always clear, although this kind of experi-
ment does give useful information. A similar experimental approach has
been utilized by other investigators (C7a, F12).
FIG.74. Axial concentration profile of bed gas, observed by taking out the gas by an
axially traveling sampler (F17).
404 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
This is the general relation to obtain the local overall rate constant koR(z)
from experimental runs. When the reaction is of first order,f(c) = c, then
koR(d/ UG = d[ln(co/cO)l/dz (8-4)
Thus, the local slope of the plot of In(co/cO)versus z for the data taken
under a constant UGwill give koR(z)/UG, and from this koR(z).
z(cm)
jetting zone near the distributor, the main dense phase, and the transition
zone between the dense and dilute phases. Besides these three zones, the
dilute phase must also be considered.
From the differential reactivity test, k,, is obtained as a measure of
reactivity. Another measure is to define the local contact efficiency 7,by
the following relation:
7, = koR(Z)/+ (8-5)
where Eed is taken instead of i,in applying 7,to the dilute phase. In the
case of the dense phase, the fraction of catalyst v is directly contacting
with the bubble gas, so that koR(z)equals k,, + vk, from Eq. (7-22); then 7,
for this case is given by
7,= t$'([k,(k,,t,Ut,)-* -k (5ce)-1]-1 + V ) (8-6)
Values of qCare given in Fig. 77 for the data of Fig. 76. The line of the
dense phase in the figure is calculated by Eq. (8-6) assuming $, = 1 and v =
0. The term ko&, is 0.4 sec-' at U, = 30.5cm/sec and is obtained from the
use of the Lewis-Gilliland-Glass plot (Section VI1,C ,2).
The contact in the dilute phase and the transition zone is different from
that in the dense phase. The former is related to mass transfer between the
gas phase and an agglomerate of solid particles, whereas the latter is
mainly related to mass transfer between bubbles and emulsion including a
certain amount of directly contacting catalyst. Upon consideration of hin-
dered settling of the swarm of particles (S16, Z7), we also find contact
efficiency to be a function of the population density of particles. Normaliz-
ing by the E, of the main dense phase, &,, E,/E,,de is chosen as a variable
' A ' O m
$0.5-
OO Ib 2b 2 ( c m30) Lb o;
FIG. 77. Dependence of T~ on the height from the distributor, catalyst 600 cm3,
U G = 30.5 cm/sec (F18). Point a indicates the upper limit of the dense phase.
FLUIDlZED CATALYST BEDS 407
FIG.78. Correlation of the contact efficiency with Ee/& for the dilute phase including
the transition zone (F18).
(Zed/Ze,de for the dilute phase) to affect vC because these are primary data
which can be obtained by measuring the solid density in the bed.
Figure 78 shows this relation. The data are rather scattered, but most of
the efficiencies are inside the shaded area of the figure. The most probable
value is given by the solid line, expressed by:
v, = 1 - 0.75(E,/E,,de)0.4 (8-7)
The equation shows that the strong interaction of particles affects the
contact efficiency. The interaction seems to increase rapidly with increas-
ing particle population. Figure 78 can be used to estimate local contact
efficiency when the density distribution of solid catalyst is known.
OF CONTACT-MECHANISM
B. AXIALDISTRIBUTION CONTRIBUTIONS
So far, the local reactivity distribution has been explained in terms of
overall rate constant koR and overall contact efficiency vc, following
Furusaki et al. (F18). The overall rate constant is further split into two
parts, the local mass-transfer term and the local fraction of directly con-
tacting catalyst, by developing a concept of the local Lewis-Gilliland-
Glass plot (M27a). The relative contribution of the variables will be ex-
plained here according to (M27a).
the figure are drawn and differentiated graphically with respect to the axial
bed height z at which one of the reactants, ethylene, has been injected.
The differentiation gives, by Eq. (8-4), the local overall rate constant koR:
k R = 2.303 UG d [ l ~ g , o ( ~ o / ~ ~ ) ] / d ~ (8-8)
Figure 79 shows the axial distribution of k,, obtained in this manner.
The local reactivity is high in the distributor region and in the transition
region (z = L,) as well. Also, the reaction proceeds in the freeboard
region, indicating the progress of reaction by freely suspended catalyst.
The reactivity is rather low in the main dense phase, in contrast to that
anticipated from the direct-contact model (L12) or the bubbling-bed
model (K24). More reaction takes place in the transition zone than in the
distributor region.
The local reactivity distribution, Fig. 79, is further processed in the
following sections to determine the relative contribution of directly and
indirectly contacting catalyst to the progress of the reaction.
FIG.79. Axial distribution of the overall rate constant koRfor different catalyst activities,
catalyst 600 cms, U G = 30.5 c d s e c (M27a).
FLUIDIZED CATALYST BEDS 409
Combining Eq. (7-37) with a and F as given in Table VII for the direct
contact model (VUE) or directly modifying Eq. (8-lo), one has the follow-
ing local relation for the dense phase section (with &Z, = Z, - u ) :
OR = vkr 4- p:(1 - kR/Z~kP)(k,bab)P,=l/(l - Y/ce) (8-15)
where is defined [see also Eq. (7-31)] by:
(8-16)
Calculation shows for the data of Fig. 79 that p: = 1.21 and 1.35 fork, =
10 and 20 sec-*, respectively.
In processing data by Eq. (8-15), first assume p: = 1 and determine
(k&,)&=l, from WhiChkb andk, are obtained (Section VI). With thek, and
k,, P;* is calculated as a function of k,. Usually, further iteration is not
needed, because pp* = 1 is a good approximation.
The plot kR/( 1 - kR/ZeQ)p: versus Q/(1 - kR/Ze,%)p: for data such as
those in Fig. 79 at a given cross section should give a straight line of slope
u and intercept (k,,ba&,=l/( 1 - v/Z,) as shown in Fig. 80a. In this way, one
can estimate the local u and (kbab)&=lfrom the reaction data for the dense
phase.
Similarly, the relation for the dilute phase is obtained by modifying Eq.
(8-14). Equation (8-15) again applies by taking kOm, ko,a,, Zed, and vd,
respectively, in place of koR,kebab, Z,, and v. Also, simply assume p: = 1
for the dilute phase, since the mass transfer term is usually negligible
here.
FIG.80. Local Lewis-Gilliland-Glass plot for the data of Fig. 7 9 (a) at z = 10-20 cm;
= 50 cm (M27a).
(b) at z
in the main dense phase, but not very different. As a consequence, mass
transfer in the dense phase can be effectively expressed by a mean value
for ( k h u h ) & = l *
In contrast to this, the fraction v distributes in a complicated manner.
Of the total directly contacting catalyst (the area enveloped by the v
412 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
FIG.81. Axial distribution of local mass transfer termk,Ga, and fraction of direct contact
catalyst v for the data of Fig. 79 (M27a).
versus z curve, Fig. 81), about 15% is located in the jetting zone (z = 0-7.5
cm), another 15% in the main dense phase (z = 7.5-30 cm), about 40% in
the transition zone (z = 30-40 cm), and the remaining 30% is in the dilute
phase (z 2.40 cm). The v distribution curve implies that direct contacting
of catalyst with the reactant gases is most noticeable in the transition zone
and in the dilute phase, less so in the jetting zone, and rather lacking in the
main dense phase. We note that the relative contribution of variables given
here is based on only one experimental example that would allow separa-
tion of v and k,,,a,. More data are needed for further progress.
The kinetic approach so far presented as the successive contact mecha-
nism, in its simplest form, has been applied to very complicated fluid-bed
behavior. Nevertheless the model has consistently paralleled the observa-
tions. A fair amount of directly contacting catalyst is shown to exist in the
transition zone and in the freeboard region, with less in the jetting zone,
under aeration with relatively high gas velocity. The relative location of
the directly contacting catalyst will certainly influence the selectivity of
complex reactions. Directly contacting catalyst in the freeboard region
may sometimes produce a drastic effect on thermal stability of the bed
(R12; also Section IX).
The successive contact mechanism seems to provide a reasonable ap-
FLUIDIZED CATALYST BEDS 413
In the dilute phase the contact efficiency is rather high but the rate of
heat transfer is low, and hence the temperature may not remain uniform in
this region. Thus, the temperature effect on conversion and selectivity
becomes important when a substantial part of the entire reaction occurs in
the dilute phase. This chapter discusses salient features of the tempera-
ture effect on fluid bed reactions.
A. EFFECT REACTION
ON STEADY
2. Overall Conversion
The effect of the dilute phase on overall conversion of the reactant is
straightforward. Conversion in the bed increases if the contact efficiency
increases due to suspended catalyst in the dilute phase. If temperature
34 - 3 1 5 0.07
5
315-400.04 3
40 - 45 0.01 1.5
45 - 50 0.001 0.1
50 - 0.005 0.5 I
50 90
z (cm)
FIG.82. Temperature profile in a fluid bed. Standard catalyst 600 cm3. Bed diameter
5.5 cm. Solid line shows calculated value using Eqs. (9-7) and (9-8). T, = temperature at
the axis.
FLUIDIZED CATALYST BEDS 415
increases in the dilute phase, conversion will further increase. Thus, from
the viewpoint of conversion of the reactant, it is recommended to use the
direct contact effect in the dilute phase and transition zone. Stability and
selectivity must now be discussed more thoroughly.
3. Effect on Selectivity
Industrially, selectivity is often as important as conversion in consider-
ing the efficiency of the reactor. In isothermal reactions, the dilute phase
and transition zone may cause better selectivity due to better contact in
that region. But in nonisothermal reactions, the effect will be different
because of the temperature effect. Mixing of gas and solids in the dilute
phase is not sufficient, and this may cause a temperature distribution for
exo- or endothermic reactions.
As an illustration, a simplified model will be considered here to show
the importance of the nonisothermal effect in the dilute phase. We assume
that cp 2 10, reaction in the bubble phase is negligible, and temperature in
the dense phase is uniform. Then the material balance equation for the
dense phase is:
uf(dcb/dz) + kbab(cb - ce) = 0 (9-1)
kbah(Ch - ce) - EehCe = (9-2)
and for the dilute phase,
U,(dc/dz) + qccdAre-E/RTc
=0 (9-3)
where the temperature is a radially averaged value. The above equations
are rearranged into dimensionless forms:
1.01
I n=o
-
1 2 3 L
Nob
1 2
FIG.83. Selectivity of isothermal consecutive reactions (A B +C).
418 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
I 2
FIG. 84. Selectivity of nonisothermal consecutive reactions (A + B +C).N,, = I,
Z, = 2, El = 83.6 kJ/mole, NArl = 2 X 108, NAa
= NArl/Io.
beds are exothermic with E2 > El and E2 >> 10 kcal/mole. In these reac-
tions it is profitable to suppress the role of the catalyst in the dilute phase.
In case of endothermic reactions, the yield of B does not vary signifi-
cantly from isothermal reactions. For E2 > 30 kcal/mole, the use of the
dilute phase even becomes profitable. Therefore, it may be desirable to
use the dilute phase in case of endothermic reactions.
The case of solids mixing with consecutive reactions is shown in Fig.
85. The magnitude of solid mixing is considerable (Npe= O.l), but the
selectivity is affected strongly by the heat of reaction, especially for the
case of highly exothermic reactions. The qualitative effect of dilute phase
on selectivity is not affected by solid mixing; therefore the effect will be
discussed for the case of negligible solid mixing, at smaller values of
effective R.
b. Denbigh's System with Optimal Temperature Distribution. Denbigh
(D11)proposed the following reaction system as a general case for organic
reactions;
3
AAX+Y
FLUIDIZED CATALYST BEDS 419
1 2
FIG. 85. Variation of selectivity of consecutive reaction (A 4 B + C ) . N , , = 0,
NM1 = 2 x lo8, NAr2= 2 x lo', n2= a,, NP. = 0.1, El = E2 = 83.6 kJ/mole.
4. Discussion of Selectivity
In discussing the yield from a reactor, the temperature distribution
inside the reactor must be investigated. In the dense phase of fluid beds,
the heat-transfer coefficient between the bed and wall has been, widely
studied (L10, M14, M15, M16, T22, V5, W3, W5). Botterill (B12) has
reviewed the recent literature; studies of heat transfer in the dilute phase
are quite limited in number. Shirai (S9, S l l ) , Furusaki (F15), and
Morooka et al. (M5O) studied the heat-transfer coefficient in the dilute
phase as well as in the dense phase. They found that the heat-transfer
coefficient between the bed and the wall decreases as the bed density
decreases, which will cause an axial distribution of temperature in bed.
0.4,
os'l
OO
, , , , . ,
5
Nob
, , , , ,
B. STABILITY BEDREACTORS
OF FLUID
phenomena in the freeboard and also in the spaces under buffle plates
inside the dense bed-the so called dilute-phase regions. Good contact
and low heat-transfer rates are the major causes of runaway. The most
significant temperature rise is considered to occur by a coupling effect of a
medium degree of temperature rise due to good contact in the dilute phase
(including the transition zone) and of noncatalytic reactions which are
initiated by the temperature rise. Thus a model reaction to consider the
possibility of runaway is:
z (-1
FIG.88. Temperature rise by noncatalytic reaction with axial mixing of solid particles.
= 2 X lo', NAr2= NArl/lO,NAr3 = 1.08 x IOI3, El = E2 = 83.6 kJlmole, E3 = 167
kJ/mole, Npe = 0.1, N,, = SO, 8, = 0.9, a, = 0, = I , O3 = SO.
large activation energy. The value NHw= 50 is that obtained for the
dilute phase with an effective diameter of about 5 cm. The values
NArl= 2 x 1O'and El = 83.6 kJ/mole give k,,L,/Uf = 5 ; 113 = 50 means
-AH3 = 2800 kJ/mole at c,! = 6.4 x lo-' mole/cm3 (-30%) at 573°K
and 1 atm. From Fig. 88 it is obvious that there is considerable temperature
rise in the dilute phase, especially due to the noncatalytic reactions.
The reaction system is very sensitive to NPe and NHw. The stabilizing
effect of the dense phase, due to good thermal conductivity, is completely
lost in the dilute phase. If Npe is larger, the temperature rise at the bottom
of the dilute phase becomes so large as to be uncontrollable.
The case of large Nre may be studied by calculation for the case of
unmixed solid particles. Calculation of this case shows that cooling of the
dilute phase is essential for steady-state operation of the reactor for such
unstable reactions. It is hard to recommend using the dilute phase, be-
cause side reactions should be avoided (Fig. 89).
The extreme case of small N P eis the case of the thermally completely
mixed dilute phase. As shown above, multiple steady states are possible.
The analyses by heat generation and rejection curves (V4)are shown in
Figs. 90 and 91. Q is the dimensionless heat being generated or rejected,
and 8 is the dimensionless temperature in the dilute phase. The value of
NHw(Lt - L f ) / L ,is 100-200 for the transition zone and 30-50 for the dilute
phase. There are two stable steady states and one unstable point for E3 of
167 and 293 kJ/mole. If E3 is reduced to 40 kJ/mole, only one stable point
424 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
z [-I
FIG.89. Runaway in the dilute phase. NArl = 2 X 108, NA, = NArl/lO,NAr3= 2.92 X
l P 4 , N , = 31, R, = sl, = 1, sl, = 50, T, = IWC, N H w= 1.17 El = E2 = 83.6 kJ/mole,
E3 = 293 kJlmole.
0 [-I
FIG. 90. Thermal effect in the dilute phase, with thermally complete mixing.
El = E2 = 83.6 kJ/mole Es = 167 kJ/mole, Tde = 3WC, R, = $2, = I , sl, = 50,
NAr1 = 2 x lon, NArz = NArl/lO, NArS = I x L,/Lf= 3.42.8, = TIT,.
e 1-1
FIG.91. Thermal effect in the dilute phase, thermally complete mixing. Parameters
are same as those in Fig. 90 except Es = 293 kJ/mole NArs = 2.9 x
FLUIDIZED CATALYST BEDS 425
A. APPLICABILITY
O F REACTOR
MODELS
The main advantages of the reactor models are:
i. Reaction performance can be expressed by the minimum number
of measurable parameters.
u. The relationship between the parameters and variables, e.g.,
catalyst properties, reactor design, and operational variables, be-
... comes clear and easy to apply.
111. The equations obtained are simple and easy to use.
Most of the models for fluidized catalytic reactions have not devoted
enough attention to flow properties of FCBs. The important features of the
FCB are as follows: (a) Bubbles grow much more slowly than in teeter
beds; (b) A circulating flow exists, centrally upward and peripherally
downward; (c) Above the dense phase are the transition region and the
dilute phase, where particle density decreases gradually with height.
These properties introduce another uncertainty in developing a large-
scale FCB, which must be operated at high superficial velocity (about 50
c d s e c ) and at a large bed height (5-10 m), with more than 30 percent of
fines fraction. The successive contact mechanism (M25, M26) is the
simplest theory dealing with FCB flow properties. The necessary pararn-
eters to calculate reactor performance are as follows: (a) bubble diameter;
(b) particle distribution through the bed height, which may be obtained
from static pressure distribution; (c) contact efficiency distribution in the
transition region and the dilute phase. If relationships between these three
parameters and the operational variables are obtained from FCB perfor-
mance in more than two laboratory steps, it is possible to calculate the
same reaction schemes in large-scale fluidized catalyst beds.
The direct contact model of Lewis et al. (L12) is a special case of the
successive contact mechanism. If the effect of circulation of the emulsion
426 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
B. DEVELOPMENT
OF INDUSTRIAL
FLUIDIZED
CATALYST
BEDS
5 . Fischer-Tropsc h Synthesis
The synthesis of hydrocarbons from H2and CO is strongly exothermic.
The Hydrocol process mentioned in Section I is a modified Fischer-
Tropsch synthesis, using powdered iron catalyst in a dense-phase
fluidized bed. In spite of elaborate scale-up based on intensive research,
many problems arose in commercial plants, and the operation was termi-
nated in 1957 for economic reasons (Z5).
The major difficulty of this process arose from the physically and chemi-
cally unstable catalyst. The iron particles broke down in size very rapidly
under the synthesis conditions. The fluffy nature of the product made
fluidization difficult, if not impossible. Also, catalyst particles were
“waxed up,” i.e., coated with carbonaceous materials, and activity de-
creased significantly (S22).
Another difficulty of this process was low conversion in the commercial
plant, which could not be predicted from the results of the pilot plants.
Grekel et al. (G15) reported the effects of particle size distribution, gas
inlet devices, and internals, on contact efficiency. Volk et al. (V12) also
emphasized the effect of bed internals. These developmental studies be-
came a very useful guide for applications of fluidized catalyst beds.
The Kellogg process has been successful for its dilute-phase transfer-
line reactor in Sasol’s South African plant. The void fraction in the reac-
tion zone is more than 95%, and the superficial gas velocity ranges from 3
to 12 d s e c (G2). Careful consideration is thus necessary in applying FCB
FLUIDIZED CATALYST BEDS 429
to processes where the catalyst can suffer physical and chemical attrition,
and is so heavy that good fluidization is hard to attain.
c. RECENT
TRENDS
IN F L U I D I Z ECATALYST
D BEDS
Recently many inventions about fluidized catalyst beds, which are use-
ful in industrial applications, have been reported (16-8), initiating trends
430 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
D. TECHNICAL
PROBLEMS
I N FCB DESIGN
of gas mixing and gas backmixing and by the ratio of length to diameter of
the bed.
2. Desired diameter: This is determined by gas velocity, which will be
in the range between the minimum fluidization velocity and the velocity
determined by excessive entrainment. The most economic velocity will be
set by balancing gas circulation cost with other costs.
3. Type of gas distributor: There are several types such as conical
bottom, packed bed, and grid.
4. Freeboard above the reactor bed: Determined by the need for dust-
removal equipment, and for free settling prior to entry to that equipment.
5 . Disengaging-space diameter: The enlarged diameter will lower the
gas velocity, aid settling of particles, and reduce carry-out.
6. Equipment required for heat exchange, if any.
Even now, the above considerations are the basis for FCB design.
Elaborate devices are installed on the gas distributor and internals for
high-performance reactions (Section I).
A panel discussion on fluid particle technology (Philadelphia, 1973) was
significant in correctly recognizing the central problems in industrial a p
plication of fluidized beds to catalytic reactions (H2). Some of the com-
ments are communicated in the following paragraphs.
Optimum size distribution is important for a fluid bed reactor (Bergoug-
nou). Models based on bubbles are not yet capable of predicting the wall
effect (Wen). Vertical baffles are most effective in breaking up large bub-
bles (Volk). The height of the bottom ends of vertical tube bundles above
the grid will set the attainable bubble size at the bottom of the bundle. The
bundles then essentially maintain the bubble size (Zenz). Horizontal per-
forated baffle plates reduce the mean residence time of elutriable fine
particles in a fluidized bed (Buckham). Observations on attrition in cy-
clones indicate that it is an exponential function of velocity (Tenney).
In general, industrial application of FCBs is possible from pilot-plant
tests (including cold runs) of appropriate scale, with the 35-year experi-
ence of industrial plants and research since the initiation of FCC. How-
ever, the following conditions are necessary. First, the catalyst should
have such properties and activity as to be suitable for FCB. Second, the
heat of reaction must be appropriate for use in FCB, considering heat
transfer and the heat capacity of fluidized beds. Third, the reaction prod-
ucts must be stable, and high selectivity must not be required if consecu-
tive reactions are involved. Generally, the catalyst is in an equilibrium
state in industrial reactors due to physical and chemical impacts, and
artificial controls from outside corresponding to those impacts.
Problems on the optimization of industrial fluidized catalyst beds are
still left unsolved. Certainly optimization in a rigorous sense is impossible.
432 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
Nomenclature
References
Al. Ahlborn, F., 2. Tech. Phys. 12, 482 (1931).
A2. Akehata, T., Shirai, T., Sugita, M., Yoshino, K., and Shirai, K., Hokkaido Meet.
Soc. Chem. Eng. Jpn. Preprinr No. B8 (1968).
A3. Akita, K., and Yoshida, F., Ind. Eng. Chem. Process Des. Dev. 12, 76 (1973).
A4. Akita, K., and Yoshida, F., Ind. Eng. Chem. Process Des. Dev. 13, 84 (1974).
AS. Alexander, B. F., and Shah, Y. T., Chem. Eng. J . 11, 153 (1976).
A6. Anderson, T. T., AlChE J . 10, 776 (1964).
A7. Anonymous, Chem. Eng., Oct., 392 (1953).
A8. Anonymous, Hydrocarbon Process. 46, 141 (1967); 48, 146 (1969); 52, 99 (1973).
A9. Aoyama, Y., Ogushi, K., Koide, K., and Kubota, H., J . Chem. Eng. Jpn. 1, 158
(1968).
AIO. Argo, W. B., and Cova, D. R . , Ind. Eng. Chem. Process Des. Dev. 4, 352 (1965).
A1 I. Ark, R., Proc. R . SOC.A235, 67 (1956).
A12. Askins, J. W., Hinds, G. P., Jr., and Kunreuther, F., Chem. Eng. Progr. 47, 401
( 1951).
438 MIYAUCHI, FURUSAKI, MOROOKA, AND IKEDA
B1. Baeyens, J., and Geldart, D., “Fluidization and Its Applications,” p. 263. Toulouse,
France, 1973.
B2. Bailie, R. C., Chung, D. S., and Fan, L. T., Ind. Eng. Chem. Fundam. 2,245 (1963).
B3. Baird, M. H. I., and Rice, R. G., Chem. Eng. J. 9, 171 (1975).
B4. Bakker, P. J., and Heertjes, P. M., Chem. Eng. Sci. 12, 260 (1%0).
B5. Bartholomew, R. N., and Casagrande, R. M., Ind. Eng. Chem. 49, 428 (1957).
B6. Basov, V. A., Makhevka, V. I., Malik-Akhnozarov, T. Kh., andOrochko, D. I., Int.
Chem. Eng. 2, 263 (1%2).
B7. Bedura, R., Deckwer, W. D., Warnecke, H. J., and Langemann, H.,Chem. Ing.
Tech. 46, 399 (1974).
B8. Behie, L. A., Doctoral dissertation, Western Ontario University, 1972.
B9. Behie, L. A., and Kehoe, P., AIChE J. 19, 1070 (1973).
B10. Bellman, R.,and Pennington, R. H., Q.Appl. Math. 12, 151 (1954).
B11. Betts, W. D., Ind. Chem. July, 370 (1963).
B12. Botterill, J. S. M., “Fluid-Bed Heat Transfer.” Academic Press, New York, 1975.
B13. Botton, R. J., Chem. Eng. Prog. Symp. Ser. 66(101), 8(1970).
B14. Bridge, A. G., Lapidus, L., and Elgin, J. C., AIChE J. 10, 819(1%4).
B15. Broadhurst, T. E., and Becker, H. A., AIChE J . 21, 238 (1975).
B16. Brotz, W., Chem. Ing. Tech. 28, 165 (1956).
B17. Burgess, J. M., and Calderbank, P. H., Chem. Eng. Sci. 30, 1511 (1975).
B18. Burkur, D., and Amundson, N. R., Chem. Eng. Sci. 30, 847, 1159 (1975).
B19. Burkur, D., Wittmann, C. V., and Amundson, N. R., Chem. Eng. Sci. 29, 173(1974).
B20. Bart, R., Sc.D. Thesis, Mass. Inst. of Tech, Cambridge, Mass., 1950, cited by
G. H. Reman (R4).
B21. Bohle, W., and Swaay. W. P. M., “Fluidization,” Proc. 2nd. Eng. Foundation
Conf. 1978, p. 167, Cambridge Univ. Press, London and New York (1978).
B22. Bukur, D. B., Ind. Eng. Chem. Fundam. 17, 120 (1978).
C1. Callahan, J. L., Grasselli, R. K., Milberger, E. C., and Strecker, H. A., Ind. Eng.
Chem. Prod. Res. Dev. 9, 134 (1970).
C2. Carlsmith, L. E., and Johnson, F. B., fnd. Eng. Chem. 37, 451 (1945).
C3. Chiba, T., and Kobayashi, H., Chem. Eng. Sci. 25, 1375 (1970).
C4. Chiba, T., Terashima, K., and Kobayashi, H., Chem. Eng. Sci. 27, 965 (1972).
C5. Chiba, T., Terashima, K., and Kobayashi, H., J . Chem. Eng. Jpn. 6, 78 (1973).
C6. Clift, R., Grace, J. R., and Weber, M. E., Ind. Eng. Chem. Fundam. 13, 45 (1974).
C7. Clift, R., and Grace, J. R., Chem. Eng. Prog. Symp. Ser. 116 67, 23 (1971).
C7a. Chavarie, C., and Grace, J. R., Ind. Eng. Chem. Fundam. 14, 75, 79, 86 (1976).
C7b. Chavarie, C., and Grace, J. R., Chem. Eng. Sci. 31, 741 (1976).
C8. Cleland, F. A., and Wilhelm, R. H., AlChE J. 2, 489 (1956).
C9. Calderbank, P. H., and Toor, F. D., in “Fluidization” (J. F. Davidson and D. Hani-
son, eds.), p. 383. Academic Press, New York, 1971.
DI. Danckwerts, P. V., Trans. Faraday SOC.46, 300 (1952).
D2. Danckwerts, P. V., Chem. Eng. Sci. 2, 1 (1953).
D3. Davidson, J. F., and Harrison, D., “Fluidized Particles.” Cambridge Univ. Press,
London and New York, 1963.
D4. Davidson, J. F., and Harrison, D., Chem. Eng. Sci. 21, 731 (1966).
D5. Davidson, J. F., and Harrison, D. (eds.), “Fluidization.” Academic Press, New
York, 1971.
D6. Davies, L., and Richardson, J. F., Trans. Inst. Chem. Eng. 44, T293 (1966).
D7. de Groot, J. H., Proc. Int. Symp. Fluidization, Eindhoven, p. 348.
FLUIDIZED CATALYST BEDS 439
D8. de Vries, R. J., van Swaaij, W. P. M., Mantovani, C., and Heijkoop, A., Proc. 5th
Eur., 2nd I n t . Symp. Chem. React. Eng. B9-59 (1972).
D9. Deckwer, W. D., Burckhart, R., and Soll, G., Chem. Eng. Sci. 29, 2177 (1974).
D10. DeMaria, F.,and Longfield, J. E., Chem. Eng. Prog. Symp. Ser. 38 58, 16 (1962).
DII. Denbigh, K. G., Chem. Eng. Sci. 8, 125 (1958).
D12. Denbigh, K. G., “Chemical Reactor Theory.” Cambridge Univ. Press, London and
New York, 1965.
D13. Dimotakis. P. E., and Brown, G . L., J . Fluid Mech. 78, 535 (1976).
D14. Deckwer, W., Graeser, U., Langemann. H., and Serpemen, Y., Chem. Eng. Sci. 28,
1223 (1973).
D14a. Donsi, G., and Massirnilla, L., “Fluidization and Its Applications,” p. 41. Capadues,
Toulouse, 1973.
D15. Donsi, G., and Massimilla, L., AIChE J. 19, 1104 (1973).
D16. Drinkenburg, A. A. H., Proc. Int. Symp. Fluidization, Eindhoven, p. 468 (1967).
D17. Drinkenburg, A. A. H., and Rietema, K., Chem. Eng. Sci. 27, 1765 (1972).
D18. Drinkenburg, A. A. H., and Rietema, K., Chem. Eng. Sci. 28, 259 (1973).
D19. Darton, R. C., Trans. Inst. Chem. Eng. 57, 134(1979).
D20. De Lasa, H. I. and Grace, J. R., AIChE J. 25, 984 (1979).
El. Eberly, P. E., J. Phys. Chem. 65, 68 (1961).
E2. Eberly, P. E., J . Phys. Chem. 66, 812 (1962).
E3. Eberly, P. E., and Kimberlin, C. N., Trans. Faraday Soc. 57, 1169 (1961).
E4. Echigoya, E.,Toyoda, K., and Morikawa, K., Kagaku Kogaku 32, 364 (1968).
E5. Echigoya, E., Iwasaki, M., Kanemoto, T., and Niiyama, H., Kagaku Kogaku 32,
571 (1%8).
E6. Eissa, S. H., El-Halwagi, M. M., and Saleh, M. A., Ind. Eng. Chem. Process Des.
Dev. 10, 31 (1971).
E7. Eissa, S. H., and Schiigerl, K . , Chem. Eng. Sci. 30, 1251 (1975).
E8. El Halwagi, M. M., and Gomezplata, A,, AIChE J. 13, 503 (1967).
E9. El Nashaie, S., Chem. Eng. Sci. 32, 295 (1977).
E10. El Nashaie, S., and Yates, J. G., Chem. Eng. Sci. 28, 515 (1973).
E l l . Endoh, I., Furusawa, T., and Matsuyama, H., Kagaku Kogaku 41, 232 (1977).
E12. Engel, W. F., and Waale, M. J., Chem. Ind. Jan. 76 (1962).
Fl. Fair, J. R., Lambright, A. J., and Anderson, J. W., Ind. Eng. Chem. Process Des.
Dev. 1, 33 (1962).
F2. Fan, L. T., Lee, C. J., and Bailie, R. C., AIChE J. 8, 239 (1962).
F3. Finnerty, R. G., Maa, J. R., Vossler, A. M., Yeh, H. S., Crouse, W. W., and Rice,
W. J., Ind. Eng. Chem. Fundam. 8, 271 (1%9).
F4. Fleurke, K. H., Chem. Eng. March, CE41 (1968).
F5. Fournol, A. B., Bergougnou, M. A., and Baker, C. G. J., Can. J. Chem. Eng. 51,401
(1973).
F6. Frantz, J. F., Chem. Eng. 69, Sept., 161 (1962).
F7. Freedman, W., and J. F. Davidson, Trans. Inst. Chem. Eng. 47, T251 (1969).
F8. Froment, G. F., Adv. Chem. Ser. 109 p. 1 (1972).
F9. Fryer, C., Ph.D. thesis, Monash University, Australia, 1974.
F10. Fryer, C., and Potter, 0. E., Ind. Eng. Chem. Fundam. I t , 338 (1972).
F11. Fryer, C., and Potter, 0. E., “Fluidization and Its Applications,” p. 440. Cepadues,
Toulouse, 1973.
F12. Fryer, C., and 0. E. Potter, AIChE J. 22, 38 (1976).
F13. Furusaki, S., Kagaku Kogakic 32, 1033 (1968).
440 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
H14. Hiraki, I., Yoshida, K., and Kunii, D., Kagaku Kogaku 29, 846 (1965).
HIS. Hirama, T., Ishida, M., and Shirai, T., Kagaku Kogaku Ronbunshu 1, 272 (1975).
H16. Hovmand, S., and Davidson, J. F., Trans. Inst. Chem. Eng. 46, T190 (1968).
H17. Hovmand, S.,and Davidson, J. F., “Fluidization” (J. F. Davidson and D. Harrison,
eds.), p. 193. Academic Press, 1971.
H18. Hurt, D. M., Znd. Eng. Chem. 35, 522 (1943).
H19. Handlos, A. E., Kunstman, R. S., and Schissler, D. O., Ind. Eng. Chem. 49, 25
(1957).
H20. Hagyard, T., and Sacerdote, A. M., Ind. Eng. Chem. Fundam. 5, 500 (1966).
H21. Horio, M., and Wen, C. Y., AIChE Symp. Ser. No. 161, 73, 9 (1977).
H22. Hoebink, J. H. B. J., and Rietema, K., “Fluidization,” Proc. 2nd. Eng. Foundation
Conf. 1978, p. 327, Cambridge Univ. Press, London and New York (1978).
11. Ikeda, Y., Kagaku Kogaku 27, 667 (1963).
12. Ikeda, Y.,Kagaku Kogaku 29, 57 (1965).
13. Ikeda, Y., Kagaku Kojo 9, (12), 63 (1965).
14. Ikeda, Y., Kagaku Kikai Gijutsu 18, 191 (1%5).
15. Ikeda, Y., Kagaku Kogaku 34, 1013 (1970).
16. Ikeda, Y., Kemikaru Enjiniyaringu 973 (1977).
17. Ikeda, Y., J. Jpn Petrol. Inst. 20, 619 (1977).
18. Ikeda, Y., and Tashiro, S . , Kagaku Kogaku 29, 956 (1965).
19. Ikeda, Y., and Tashiro, S., “Saikin no Kagaku Kogaku,” p. 195. Maruzen, Tokyo,
1%7.
110. Ikeda, Y.,Tashiro, S., Yamaguchi, T., and Kawai, M., Annu. Meet. SOC. Chem.
Eng. Jpn. 35th, Preprint D202 (1970).
I l l . Ishii, T., and Osberg, G. L., AIChE J . 11, 279 (1965).
112. Ivanov, M. E., and Bykov, V. P., Theor. Found. Chem. Eng. 4, 119 (1970).
113. Iwasaki, M.,Furuoya, I., Sueyoshi, H.,Shirasaki, H., and Echigoya, E., Kagaku
Kogaku 29, 892 (1965).
J1. Jackson, R., Trans. Inst. Chem. Eng. 41, 21 (1963).
J2. Johnstone, H. F., Batchelor, J. D., and Shen, C. Y., AIChE J . I, 318 (1955).
K1. Kaji, H., M.S. thesis, Dept. of Chem. Eng., University of Tokyo, 1%6.
K2. Kamimura, S . , J . Jpn. Petrol. Inst. 20, 627 (1977).
K3. Kamimura, S . , Kagaku Kogaku 33, 1082 (1973).
K4. Kato, K., and Wen, C. Y., Chem. Eng. Sci. 24, 1351 (1969).
K5. Kato, Y.,Kagaku Kogaku 27, 366 (1963).
K6. Kato, Y., and Nishiwaki, A., Kagaku Kogaku 35, 912 (19711, I n t . Chem. Eng. 5,
182 (1972).
K7. Kato, Y., Nishinaka, M., and Morooka, S . , Kagaku Kogaku Ronbunshu 1, 530
( 1975).
K8. Kehoe, P. W.K., and Davidson, J. F., AIChE Symp. Ser. 128 69, 34 (1973).
K9. Kehoe, P. W.K., and Davidson, J. F., AIChE Symp. Ser. 128 69, 41 (1973).
K10. Kimura, T., M.S. thesis, Dept. of Chem. Eng., University of Tokyo, March, 1972.
K11. Klinkenberg, A., Jr., Trans. I n s t . Chem. Eng. 43, TI41 (1965).
K12. Kobayashi, H., “Saikin no Kagaku Kogaku,” p. 17. SOC.Chem. Eng., Japan, 1966.
K13. Kobayashi, H., Arai, F., and Chiba, T., Kagaku Kogaku 29, 858 (1965).
K14. Kobayashi, H., and Arai, F., Kagaku Kogaku 29, 885 (1965).
K15. Koide, K.,Morooka, S., Ueyama, K., Matsuura, A., Yamashita, F., Iwamoto, S.,
Kato, Y., Inoue, H., Suzuki, S.,and Akehata, T., J. Chem. Eng. Jpn. 12, 98 (1979).
K16. Kojima, E., Akehata, T., and Shirai, T., J . Chem. Eng. Jpn. 8, 108 (1975).
K17. Kolbel, H., Beinhauer, R.,and Langemann, H., Chem. Ing. Tech. 44,697 (1972).
442 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
K18. Kolbel, H., Borchers, E., and Muller, K., Chem. Ing. Tech. 30, 729 (1958).
K19. Kolbel, H., Borcher, E., and Martins, J., Chem. Ing. Tech. 32, 84 (1960).
K20. Kumar, A., Degaleesan, T. E., Laddha, G.S.,and Hoelscher, H. E., Can. J. Chem.
Eng. 54, 503 (1976).
K21. Kumar, R.,and Kuloor, N. R., “Advances in Chemical Engineering,” Vol. 8, p. 256.
Academic Press, New York, 1970.
K22. Kunii, D., and Levenspiel, O., Ind. Eng. Chem. Fundam. 7, 446 (1968).
K23. Kunii, D., and Levenspiel, O., Ind. Eng. Chem. Process Des. Dev.7, 481 (1968).
K24. Kunii, D., and Levenspiel, O., “Fluidization Engineering.” Wiley, New York, 1969.
K25. Kunii, D., and Levenspiel, O., J. Chem. Eng. J p n . 2, 84 (1969).
K26. Kunii, D., Akiyama, T., and Takagi, K., “Shin Kagaku Kogaku Koza,” IV-5, Fluidi-
zation. Nikkan-Kogyo-Shimbun, Tokyo, 1957.
K27. Kunii, D., Yoshida, Y., and Hiraki, I., Proc. Int. Symp. Fluidization, Eindhoven,
p. 243 (1967).
K28. Kono, T., B. Eng. thesis, Dept. of Appl. Chem., Kyushu University, 1976.
K29. Khan, A. R., Richardson, J. F., and Shakiri, K. J., “Fluidization,” Proc. 2nd.
Eng. Foundation Conf. 1978, p. 351, Cambridge Univ. Press, London and New
York (1978).
, K30. Kikuchi, T., private communication (1977).
L1. LaNauze, R. D., Powder Technol. 15, 117 (1976).
L2. Lanneau, K. P., Trans. Inst. Chem. Eng. 38, 125 (1960).
L3. Larroux, G . J., Kim, Y. G.,and Bankoff, S . G.,Chem. Eng. Sci. 27, 447 (1972).
L4. Latham, R. L., and Potter, 0. E., Chem. Eng. J. 1, 152 (1970).
L5. Latham, R. L., Hamilton, C., and Potter, 0. E.,Br. Chem. Eng. 13, 666 (1%8).
L6. Leva, M., Chem. Eng. Prog. 47, 39 (1951).
L7. Leva, M., “Fluidization.” McGraw-Hill, New York, 1959.
L8. Leva, M., and Grummer, M., Chem. Eng. Prog. 48, (6), 307 (1952).
L9. Levenspiel, O., and Bischoff, K. B., “Advances in Chemical Engineering,” Vol. 4,
p. 95. Academic Press, New York, 1963.
L10. Levenspiel, O., and Walton, J. S . , Chem. Eng. Prog. Symp. Ser. 50, l(1954).
L11. Lewis, W. K., Gilliland, E. R., and Girouard, H., Chem. Eng. Prog. Symp. Ser. 38
58, 87 (1962).
L12. Lewis, W. K., Gilliland, E. R., and Glass, W., AIChE J . 5, 419 (1959).
L13. Lewis, W. K., Gilliland, E. R.,and Lang, P. M., Chem. Eng. Prog. Symp. Ser. 38 58,
65 (1962).
MI. MIT, A Bibliography of Fluidization Research at the Massachusetts Institute of
Technology, 1939- 1955.
M2. Marshall, S . , Chem. Eng. May, 219 (1953).
M3. Massimilla, L., AIChE Symp. Ser. 128 69, 11 (1973).
M4. Massimilla, L., and Johnstone, H. F., Chem. Eng. Sci. 16, 105 (1961).
M5. Massirnilla, L., and Westwater, J. W., AIChEJ. 6, 137 (1960).
M6. Massimilla, L., Donsi, G.,and Zucchini, C., Chem. Eng. Sci. 27, 2005 (1972).
M7. Matheson, G. L., Herbst, W. A., and Holt, P. H., Ind. Eng. Chem. 41, 1099 (1949).
M8. Mathis, J. F., and Watson, C. C., AIChE J. 2, 518 (1956).
M9. Matsen, J. M., Chem. Eng. Prog. Symp. Ser. I01 66, 47 (1970).
M10. Matsen, J. M., AIChE Symp. Ser. 128 69, 30 (1973).
M11. Matsen, J. M., and Tarmy, B. L., Chem. Eng. Prog. Symp. Ser. 101 66, 1 (1970).
M12. May, W. G., Chem. Eng. Prog. 55, (12), 49 (1959).
M13. Merrick, D., and Highley, J., AIChE Symp. Ser. 137 70, 336 (1974).
M14. Mickley, H. S., and Fairbanks, D. F., AlChE J. 1, 374 (1955).
FLUIDIZED CATALYST BEDS 443
M15. Mickley, H. S., and Trilling, C. A., Ind. Eng. Chem. 41, 1135 (1949).
M16. Mickley, H. S . , Fairbanks, D. F., and Hawthorn, R. D., Chem. Eng. Prog. S y m p .
Ser. 32 57, 51 (1961).
M17. Miwa, K., Mon, S., Kato, T.,and Muchi, I., Kagaku Kogaku 35, 770 (1971); I n t .
Chem. Eng. 12, 181 (1972).
M18. Miyauchi, Takeshi, and Oyamada, N., Jpn. Petrol. Inst. 20, 624 (1977).
M19. Miyauchi, T., Sato, Y.,Michiki, H., and Fujimoto, K., Japan Patent 45-39616, 1970.
M20. Miyauchi, Terukatsu, “Ryukeisosa-to-Kongotokusei.” Nikkan-Kogyo-Sinbun,
Tokyo, 1960.
M21. Miyauchi, T., Annu. Meet. SOC. Chem. Eng. Jpn. 30th. Preprint No. 6110 (1965).
M22. Miyauchi, T., Annu. Meet. SOC.Chem. Eng. Jpn. 3Ist. Preprint No. 244 (1966).
M23. Miyauchi, T., J. Chem. Eng. Jpn 4, 238 (1971).
M24. Miyauchi, T., CHISA Congr., 4th. Prague Sept. (1972).
M25. Miyauchi, T., J. Chem. Eng. Jpn. 7 , 201 (1974).
M26. Miyauchi, T., J . Chem. Eng. Jpn. 7 , 207 (1974).
M27. Miyauchi, T., Fall Meet. SOC. Chem. Eng. Jpn. 11th. Invited Lecture, Oct. (1977).
M27a. Miyauchi, T., and Kikuchi, T., Annu. Meet. SOC. Chem. Eng. Jpn., 42nd,
Hiroshima April (1977).
M28. Miyauchi, T., and Furusaki, S., AIChE J . 20, 1087 (1974).
M29. Miyauchi, T., and Morooka, S., Kagaku Kogaku 33, 369 (1969); Int. Chem. Eng. 9,
713 (1969).
M30. Miyauchi, T., and Morooka, S., Kagaku Kogaku 33, 880 (1969).
M31. Miyauchi, T., and Shyu, C. N., Kagaku Kogaku 34, 958 (1970).
M32. Miyauchi, T., and Yamada, K.,Kagaku Kogaku 35, 547 (1971).
M33. Miyauchi, T., Kaji, H., and Saito, K., J . Chem. Eng. Jpn. 1, 72 (1968).
M34. Mori, S., and Muchi, I., Kagaku Kogaku 34, 510 (1970).
M35. Mori, S., and Muchi, I., J . Chem. Eng. Jpn. 5, 251 (1972).
M36. Mori, S., and Wen, C. Y., A K h E J. 21, 109 (1975).
M37. Mori, Y.,and Nakamura, K., Kagaku Kogaku 29, 868 (1965).
M38. Morooka, S., Dr. Eng. dissertation, University of Tokyo, 1969.
M39. Morooka, S., and Miyauchi, T., Kagaku Kogaku 33, 569 (1969).
M40. Morooka, S., Kato, Y., and Miyauchi, T., J. Chem. Eng. Jpn. 5 , 161 (1972).
M41. Morooka, S., Nishinaka, M., and Kato, Y., Kagaku Kogaku 37, 485 (1973).
M42. Morooka, S., Nishinaka, M., and Kato, Y.,Kagaku Kogaku Ronbunshu 2,71(1976);
Int. Chem. Eng. 17, 254 (1977).
M43. Morooka. S., Tajima, K., and Miyauchi, T., Kagaku Kogaku 35, 680 (1971); Int.
Chem. Eng. 12, 168 (1972).
M44. Morse, R. D., and Ballow, C. O., Chem. Eng. Prog. 47, (9), 199 (1951).
M45. Murphree, E.V., Brown, C. L.,Fischer, H. G. M., Gohr, E. J., and Sweeney, W. J.,
Ind. Eng. Chem. 35, 768 (1943); Chem. Tech. 6, 523 (1976).
M46. Murray, J. D., J . Fluid Mech. 21, 465 (1965).
M47. Murray, J. D., J. Fluid Mech. 22, 57 (1965).
M48. Muchi, I., Mem. Fac. Eng. N a g o y a Univ. 17, (I), 79 (1965).
M49. Muchi, I., Mamuro, T., and Sasaki, K., Kagaku Kogaku 25, 757 (1961).
M5O. Morooka, S., Maruyama, Y., Kawazuishi, K.,Higashi, S., and Kato, Y., Kagaku
Kogaku Ronbunshu 5, 275 (1979).
M51. Morooka, S., Kawazuishi, K. and Kato, Y., Powder Technol. 26, 75 (1980).
NI. Nguyen, H. V., and Potter, 0. E.,Adv. Chem. Ser. 133, 290 (1974).
N2. Nguyen, H. V.,and Potter, 0. E., “Fluidization Technology,” Vol. 2, p. 193. Hemi-
sphere, Washington, D.C., 1976.
444 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
N3. Nguyen, X. T., Bergougnou, M. A., and Baker, C. G. J., Can. J. Chem. E n g . 51,573
(1973).
N4. Nicklin, D. J., Chem. Eng. Sci. 17, 693 (1962).
N5. Nicklin, D. J., Wilkes, J. O., and Davidson, J. F., Trans. Inst. Chem. Eng. 40, 61
( 1962).
N6. Nishinaka, M., Morooka, S., and Kato, Y., Powder Technol. 9, I (1974).
N7. Nishinaka, M., Morooka, S.,and Kato, Y., Kagaku Kogaku Ronbunshu 1,81(1975).
N8. Nishinaka, M., Morooka, S., and Kato, Y., Kagaku Kogaku Ronbunshu 1, 605
(1975).
N9. Nishinaka, M., Morooka, S., and Kato, Y., Kagaku Kogaku Ronbunshu 2, 96(1976).
N10. Nguyen, H. V., Whitehead, A. B., and Potter, 0. E., “Fluidization,” Proc. 2nd.
Eng. Foundation Conf. 1978, p. 140, Cambridge Univ. Press, London and New
York (1978).
01. Ogura, Y., “Taiki-Ranryu-Ron” (Atmospheric Turbulence), 5th. ed., pp. 2-7, 68.
Chijin-Shokan, Tokyo, 1959.
02. Ohki, Y., and Inoue, H., Chem. Eng. Sci. 25. 1 (1970).
03. Oki, K., and Shirai, T., “Fluidization Technology,” Vol. 1, p. 95. Hemisphere,
Washington, D.C., 1976.
04. O h , H. L., and Dean, 0. C., Petrol. Eng. March, C-23 (1953).
05. Oltrogge, R. D., Preliminary outline of thesis, Technical University of Eindhoven and
University of Michigan, 1971.
0 6 . Orcutt, J. C., Davidson, J. F., and Pigford, R. L., Chem. Eng. Prog. Symp. Ser. 38
58, 1 (1962).
07. Ormiston, R. M., Mitchell, F. R. G., and Davidson, J. F., Trans. Inst. Chem. Eng. 43,
T209 (1965).
08. Osberg, G. L., and Charlesworth, D. H., Chem. Eng. Prog. 47, 566 (1951).
09. Ostergaard, K., “Advances in Chemical Engineering,” Vol. 7, p. 71. Academic
Press, New York, 1968.
010. Othmer, D. F., “Fluidization,” Reinhold, New York, 1956.
011. Overcashier, R. H., Todd, D. B., and Olney, R. B., AIChE J. 5, 54 (1959).
012. Ogasawara, S., Kihara, M., Nishiyama, M., Shirai, T., and Morikawa, K., Kagaku
Kogaku 28, 59 (1964).
PI. Pansing, W. F., AIChE J . 2, 71 (1956).
P2. Partridge, B. A., and Rowe, P. N., Trans. Inst. Chem. Eng. 44, T335 (1966).
P3. Pavlov, V. P., Khim. Prom. (9), 698 (1965).
P4. Perkins, T. K., and Johnston, 0. C., SOC. Petrol. Eng. J . 70 (March 1963).
P5. Potter, 0. E., and Thiel, W. “Fluidization Technology,” Vol. 2, p. 185. Hemisphere,
Washington, D.C., 1976.
P6. Pozin, L. S.,Aerov, M. E., and Bystrova, T. A., Theor. Found. Chem. Eng. 3,714
(1969), (English trans.).
W.Prandtl, L., ZAMM 22,241 (1942); see H. Schlichting, “Boundary Layer Theory,”
Chap. 19, 4th English ed. McGraw-Hill, New York, 1960.
P8. Pyle, D. L., and Harrison, D., Chem. Eng. Sci. 22, 1199 (1967).
P9. Pyle, D. L., and Rose, P. L., Chem. Eng. Sci. 20, 25 (1965).
PIO. Pyle, D. L., and Stewart, P. S . B., Chem. Eng. Sci. 19, 842 (1964).
P11. Pyle, D. L., and Harrison, D., Chem. Eng. Sci. 22, 531 (1967).
R1. Raghuraman, J., and Verma, Y.B. G., Chem. Eng. Sci. 30, 145 (1975).
R2. Reichart, H., VDI Forschungsh. (Berlin) 414 (1942).
R3. Reith, T., Renken, S.,and Israel, B. A., Chem. Eng. Sci. 23, 619 (1968).
R4. Reman, G. H., Chem. Ind. Jan., 46 (1955).
FLUIDIZED CATALYST BEDS 445
T5. Takahashi, M., M.S. thesis, Dept. Chern. Eng., University of Tokyo, 1977.
T6. Tanaka, I., and Shinohara, H., J. Chem. Eng. Jpn. 5, 57 (1972).
“7. Tanaka, I., Shinohara, H., Hirose, H., and Tanaka, Y.,J. Chem. Eng. Jpn. 5, 5 1
( 1972).
T8. Taylor, G. I., Proc. R . SOC. (London) Ser. A201, 192 (1950).
T9. Taylor, G . I., Proc. R . SOC. (London) Ser. A219, 186 (1953).
T10. Taylor, G. I., Proc. R . SOC.(London) Ser. A223, 446 (1954).
TI1. Taylor, G. I., Proc. R . SOC.(London) Ser. A225, 473 (1954).
T12. Terahata, T., Tazawa, S., Kakizaki, T., Minoda, N., Miyazima, S., Ito, M.,Imai, H.,
Kawazu, T., Japan patent No. 49-29166, 1974.
T13. Tett, H. C., “Fluidization,” p. 1’. SOC.Chem. Ind., London, 1%4.
T14. Tichacek, L. J., Barkelew, C. H., and Baron, T., AIChE J. 3, 439 (1957).
T14. Todt, J., Liicke, J., Schiigerl, K., and Renken, A., Chem. Eng. Sci. 32, 369 (1977).
T16. Toei, R., and R. Matsuno, Proc. I n t . Symp. Fluidization, Eindhoven, p. 271 (1967).
T17. Toei, R., Matsuno, R., Oichi, M., and Yamamoto, K., J . Chem. Eng. Jpn. 7 , 447
( 1974).
T18. Toei, R., Matsuno, R., Hotta, H., Oichi, M., and Fujine, Y.,J. Chem. Eng. Jpn. 5,
273 (1972).
T19. Toei, R., Matsuno, R., Kojima, H., Nagai, Y.,Nakagawa, K., and Yu,S . , Kagaku
Kogaku 29, 851 (1965).
T20. Twi, R., Matsuno, R., Nishitani, K., Hayashi, H., and Imamoto, T., Kagaku
Kogaku 33, 668 (1969).
T21. Tone, S . , Seko, H., and Otake, T., J. Chem. Eng. Jpn. 7, 44 (1974).
T22. Toomy, R. D., and Johnstone, H. F., AIChE Symp. Ser. 5 49, 51 (1953).
T23. Towell, G. D.,Strand, C. P., and Ackerman, G. H., Paper presented at Am. Inst.
Chem. Eng. and Inst. Chem. Eng. Joint Meeting, p. 1272, London (1965).
T24. Towell, G. D., and Ackerman, G. H., Proc. 5th Eur., 2nd Int. Symp. Chem. Reaction
Eng. Amsterdam B3-1 (1972).
T25. Trawinski, H., Chem. Ing. Tech. 23, 416(1951).
T26. Trawinski, H., Chem. Ing. Tech. 25, 201 (1953).
T27. Tsutsui, T., and Miyauchi, T., Kagaksu Kogaku Ronbunshu 5, 40 (1979).
T28. Turner, J. C. R., Chem. Eng. Sci. 21, 971 (1%6).
T29. Tsutsui, T., M. S. thesis, Dept. of Chem. Eng., University of Tokyo, 1974.
T30. Tellis, C. B., and Hulburt, H. M., Proc. Pacgc Chem. Eng. Con$ 1st SOC. Chem.
Eng. Jpn. AIChE Kyoto Session 8, p. 178 (1972).
UI. Ueyama, K., and Miyauchi, T., Kagaku Kogaku Ronbunshu 2,430 (1976).
U2. Ueyama, K., and Miyauchi, T., Kagaku Kogaku Ronbunshu 2, 595 (1976).
U3. Ueyama, K., and Miyauchi, T., Kagaku Kogaku Ronbunshu 3, 19 (1977).
U4. Ueyama, K., and Miyauchi, T., Kagaku Kogaku Ronbunshu 3, I15 (1977).
U5. Ueyama, K., and Miyauchi, T., AIChE J . 25, 258 (1979).
V1. Van Deemter, J. J., Chem. Eng. Sci. 13, 143 (1961).
V2. Van Deemter, J. J., Proc. Int. Symp. Fluidization Eindhoven p. 334 (1967).
V3. Van der Laan, Chem. Eng. Sci. 7 , 187 (1958).
V4. Van Heerden, C., Ind. Eng. Chem. 45, 1242 (1953).
V5. Van Heerden, C., Nobel, P., and van Krevelen, D. W., Chem. Eng. Sci. 1,51(1951).
V6. Van Krevelen, D. W., and Hoftijzer, P. J., Chem. Eng. Prog. 46, (1). 29 (1950).
V7. Van Swaaij, W. P. M., and Zuiderweg, E J., Proc. 5th Eur. 2nd Int. Symp. Chem.
Reaction Eng. B9-25 (1972).
V8. Van Swaay, W. P. M., and Zuiderweg, F. J., “Fluidization and Its Applications,”
p. 454. Cepadues, Toulouse, 1973.
FLUIDIZED CATALYST BEDS 447
V8a. Van Swaaij, W. P. M., “Chem. Reaction Eng. Reviews,” (D. Luss and V. W.
Weekman Jr., eds.), A m . Chem. SOC. Symp, Ser., 72, 193 (1978).
V9. Veatch, F., Callaham, J. L., Idol, J. D., and Milberger, E. C., Perol. ReJner 41, ( l l ) ,
187 (1962).
VIO. Verma, Y. B. G., Powder Technol. 12, 167 (1975).
V l l . Vermeulen, T., Moon, J. S., Hennico, A., and Miyauchi T., Chem. Eng. Prog. 62,
(9), 95 (1966).
V12. Volk, W., Johnson, C. A., and Stotler, H. H., Chem. Eng. Prog. 58, (9,44 (1962).
V13. Vreedenberg, H. A., Chem. Eng. Sci. 9, 52 (1958).
V14. Vreedenberg, H. A., Chem. Eng. Sci. 11, 274 (1960).
V15. Venuto, P. B. and Habib, E. T. Jr., “Fluid Catalytic Cracking with Zeolite
Catalysts.” Marcel Dekker, 1979.
W1. Wakabayashi, T., and Kunii, D., J . Chem. Eng. Jpn. 4, 226 (1971).
W2. Wen, C. Y., and Hashinger, R. F., AIChE J . 6, 220 (1960).
W3. Wen, C. Y., and Leva, M., AIChE J. 2, 482 (1956).
W4. Wen, C. Y., and Yu, Y. H., Chem. Eng. Prog. Symp. Ser. 62 62, 100 (1966).
W5. Wender, L., and Cooper, G . T., AIChE J . 4, 15 (1958).
W6. Werther, J., Trans. Inst. Chem. Eng. 52, 149, 160 (1974).
W7. Werther, J., AIChE Symp. Ser. 141 70, 53 (1974).
W8. Werther, J., Powder Technol. 15, 155 (1976).
W9. Werther, J . , “Fluidization Technology,” Vol. 1, p. 2 IS. Hemisphere, Washington,
D.C., 1976.
W10. Werther, J., Int. J. Multiphase Flow 3, 367 (1977).
WlOa. Werther, J., Chem. Ing. Tech. 50, 850 (1978).
W I I . Wheeler, A., “Catalysis 11” (P. H. Emmett, ed.), p. 105. Reinhold, New York, 1958.
W12. Whitehead, A. B., Gartside, G., and Dent, D. C., Powder Technol. 14, 61 (1976).
W13. Wilhelm, R. H., and Kwauk, M., Chem. Eng. Prog. 44, 201 (1948).
W14. Wollard, J . M., and Potter, 0. E., AIChE J . 14, 388 (1968).
W15. Wan, C. G., and Ziegler, E. N., Ind. Eng. Chem. Fundam. 12, 55 (1973).
W16. Wijffels, Ir. J-B., and Rietema, K., Trans. Inst. Chem. Eng. 50, 224 (1972).
W17. Wijffels, Ir. J-B., and Rietema, K., Trans. Inst. Chem. Eng. 50, 233 (1972).
XI. Xavier, A. M., Lewis, D. A. and Davidson, J. F., Trans. Inst. Chem. Eng. 56,
274 (1978).
Yl. Yagi, S., and Miyauchi, T., Kagaku Kogaku 17, 382 (1953).
Y2. Yamagoshi, T., B.S. thesis, Dept. Chem. Eng., University of Tokyo, 1969.
Y3. Yamazaki, M., J. Chem. Eng. Jpn. 8, 420 (1975).
Y3a. Yamazaki, M., Annu. Meet. SOC.Chem. Eng. Jpn., 40th. Nagoya (1975).
Y4. Yamazaki, M., and Miyauchi, T., J. Chem. Eng. Jpn. 4, 324 (1971).
Y5. Yamazaki, M., and Miyauchi, T., Kagaku Kogaku Ronbunshu 3, 261 (1977).
Y6. Yamazaki, M., Ito, N., and Jimbo, G., Kagaku Kogaku Ronbunshu 3, 272 (1977).
Y7. Yang, W. C., and Keairns, D. L., AIChE Symp. Ser. 141 70,27 (1974).
Y8. Yang, W. C., and Keairns, D. L., “Fluidization Technology,” Vol. 2, p. 51. Hemi-
sphere, 1976.
Y9. Yates, J. G., and Constans, J. A. P., Chem. Eng. Sci. 28, 1341 (1973).
Y10. Yates, J. G., and Rowe, P. N., Trans. Inst. Chem. Eng. 55, 137 (1977).
Y11. Yerushalmi, J., Turner, D. H., and Squires, A. M., Ing. Eng. Chem. Process Des.
Dev. 15, 47 (1976).
Y12. Yerushalmi, J., Gluckman, M. J., Graff, R. A., Dobner, S., and Squires, A. M.,
“Fluidization Technology,” Vol. 2, p. 437. Hemisphere, 1976.
Y13. Yokota, T., Hidaka, Y., and Yasutomi, T., Kagaku Kogaku Ronbunshu 1,399( 1975).
448 MIYAUCHI, FURUSAKI, MOROOKA, A N D IKEDA
Y14. Yoshida, A., B.S. thesis, Dept. Chem. Eng., University of Tokyo, 1971.
Y 15. Yoshida, K., and Kunii, D., J . Chem. Eng. Jpn. 1, 11 (1%8).
Y16. Yoshida, K., and Kunii, D., J . Chem. Eng. J p n . 7, 34 (1974).
Y 17. Yoshida, K., Kunii, D., and Levenspiel, O., Ind. Eng. Chem. Fundam. 8,402 (1969).
Y18. Yoshida, K., Kunii, D., and Levenspiel, O., Znt. J . Heat Mass Transfer 12, 529
(1%9).
Y19. Yoshida, K., Ueno, T., and Kunii, D., Chem. Eng. Sci. 29, 77 (1974).
Y20. Yoshitome, H., Kagaku Kogaku 27, 27 (1963).
Y21. Yoshitome, H.. Dr. Eng. dissertation, Tokyo Inst. Tech., 1967.
Y22. Yoshitome, H., and Shirai, T., J . Chem. Eng. Jpn. 3, 29 (1970).
Y23. Yoshitome, H., Mannami, Y., Mukai, K., Yoshikoshi, N., and Kanazawa, T.,
Kagaku Kogaku 29, 19 (1965).
Y24. Yacono, C., Rowe, P. N. and Angelino, H., Chem. Eng. Sci. 34, 789 (1979).
Y25. Yerushalmi, J. and Cankurt, N. T., Chem. Tech. 8, 564 (1978).
Y26. Yoshida, K., Nakajima, K.,Hamatani, N. and Shimizu, F., “Fluidization,” Proc.
2nd. Eng. Foundation Conf. 1978, p. 13, Cambridge Univ. Press, London and
New York (1978).
Z1. Zabrodsky, S. S., “Hydrodynamics and Heat Transfer in Fluidized Bed.” Gosener-
goizdat, Moscow, 1963; English trans., MIT Press, Cambridge, Massachusetts, 1966.
22. Zabrodsky, S. S., Antonishin, N. V., and Parnas, A. L., Can. J . Chem. Eng. 54, 52
(1976).
23. Zalewski, W. c., and Hanesian, D., AIChE Symp. Ser. I28 69, 58 (1973).
24. Zenz, F. A., Petrol. Rejiner 36 (9,261 (1957).
Z5. Zenz, F. A., and Othmer, D. F., “Fluidization and Fluid-Particle Systems.”
Reinhold, New York, 1960.
26. Zenz, F. A., and Weil, N. A., AIChE J . 4, 472 (1958).
Z7. Zuber, N., Chem. Eng. Sci. 19, 897 (1964).
Some related literature that appeared subsequent to the preparation of this article may
be of interest to the reader and is therefore listed here:
Topic Reference