Pergamon
Chemical En#ineerin9 Science, Vol. 52, No. 15, pp. 2471-2482, 1997
( 1997 Elsevier Science Ltd. All rights reserved
Printed in Great Britain
P I I : S0009-2509(97)00066-3
0009-2509/97 $17.00 + 0.00
Circulating fluidized-bed reactors
Daizo Kunii* and Octave Levenspiel**
"1-25-16 Nakamachi, Meguro ku, Tokyo 153, Japan; *Chemical Engineering Department,
Oregon State University, Corvallis, OR, U.S.A.
(Accepted 15 May 1996)
Abstract--This paper develops a flow and contacting model to represent a CFB. Best estimates
of contacting efficiencies are presented for the turbulent, fast fluidized, and pneumatic transport
regimes of the CFB. Material balances are presented, ending up with conversion equations for
first-order solid-catalyzed gas-phase reactions. Four examples show how to use this model and
show its predictions. 1997 Elsevier Science Ltd
Keywords: CFB reactors; conversion; circulating solids reactors.
1. INTRODUCTION
Let us first sketch the contacting regimes encountered
when a bed of solids is fluidized by gas at progressively higher velocities. Because contacting differs in
these contacting regimes, different reactor models
must be used for predicting reactor behavior; see
Fig. 1.
This paper deals with circulating solids reactors,
These contain very fine particles which are fluidized at
a rather high gas velocity, are blown out of the bed
and reactor, and have to be replaced by fresh solids.
We call these circulating fluidized beds (CFB).
We must first find what contacting regime is involved, and then deal with the reactor and its conversion equations.
2. CHARACTERIZATIONOF PARTICLES
The equivalent spherical particle diameter is defined as
diameter of a sphere
dsph =
which has the same
(~_)1/3
=
(1)
volume as the particle, V
and the particle sphericity q~s is defined as
= ( s u r f a c e o_f a_ sphere']
qS~ \ surface of particle J . . . . . . lume
relationship between ds. and dp. The best we can say
for pressure drop considerations is
for irregular particles with no
(4)
seeming longer or shorter
dimension
for irregular particles with one
somewhat longer dimension but
with length ratio not greater
(5)
than 2:1 (eggs, for example)
for irregular particles with one
shorter dimension but with
(6)
length ratio not less than 1:2
(pillows, for example)
qSsdsCr
ds~r
dp ~q52dscr
3. GAS/SOLIDCONTACTINGREGIMES
For given particles (dp and ps) and given superficial
gas velocity through the bed (Uo), we first need to find
what contacting regime is involved--packed bed,
bubbling fluidized bed (BFB), or circulating fluidized
bed (CFB), with its subregimes--turbulent fluidized,
fast fluidized, or pneumatic transport. To do this first
evaluate the dimensionless measures of particle size
and gas velocity. These are defined as
(2)
=
From the above two expressions we define the particle
size to be used for suspended solid reactors as
dp = ~b~dsph.
(3)
For fine particles we evaluate the size by screen analysis, which gives dscr. Unfortunately, there is no general
:Corresponding author.
d V p~(p~ - pg)g],/3
Y
#(p~ - pg)g
(7)
.
(8)
3.1. Minimum fluidizing velocity
The solids will be suspended when the pressure
drop exceeds the weight of solids. This happens when
the gas velocity exceeds the minimum fluidizing velocity. This velocity is given by Ergun (1952), and in
2471
2472
D. Kunii and O. Levenspiel
Freely
bubbling
bed
-~Xmt
~, 1~;,~4~:
"~'2 :::. ::.:.:-':,"~
~@I
r:!.":-~-+?-:~.%:,~1
-'"
I:'.'..:."-.::.:<f::.~ ;.-;1<
",
- - ~
]
:, :":'~.::.,~ ..".'~.. ;;
Very low
gas velocity
I.I
_ ~-_ ~,L;:"c~"~'~~: ~'"ll
_~_.,~.._._.,,-.~ ,,
",~ I~':,':,'~:~
/~
condibons
fluidizing..
"4 E ' "
with solids
present in the
freeboard
Bed at
minimum
Fixedbed
Turbulent or fast
fluidized
:- ~,~"- ~ "~I _ Some
!.
v,~. ~, r e a c t i o n
~.. .I': Z.,I occurs
~~i'~..':~,~
" " ' , --_':,':~"i~
Uo = Umt
~ t
Slugging occurs in narrower or
small beds, or in beds with
closely spaced internals
-- -.'~--~1
~
here
..
.......
"h f (
"
"-:"
so,i<
.':~,~:~.?~.~
.
~:-:.'..
5y
circulation
s,stom
A
very high
gas velocity
Fig. l. Gas/solid contacting regimes depend on the gas velocity and bed geometry.
dimensionless form is
velocity for spherical particles as
150(1 -- emS) u. s + 1.75 (U's) 2 d* = (dp*)2
f 18
0.591 ] - 1
u* = L ~
+ -(d*)
- '/2]
(9)
~mf
~s
and for irregularly shaped particles of sphericity q~s,
3.2. T e r m i n a l velocity, ut
Individual particles are blown out of the bed when
the gas velocity exceeds what is called the terminal
velocity, ut. Haider and Levenspiel (1989) give this
Bed expands, then
forms bubbles.
These split and
coalesce and stay
about 10 cm in
size. Solids mix
vigorously,
'
'
'
I18
2.335 - 1.744~b,] -1
"'* = ~
I''"I
'
'
'
I'
(-d})~
"'I
'
~:
Sandqike
(11)
"---- Bubbles grow
large. Bed
internals help
keep bubble
size down.
-I
(10)
0.5
-
Difficult to
fluidize. Solids
clump, and rise
C "////////\. . . . . ,'!
Cohesive'//////////.
Aeratabl
as slugs. Gas
channels.
,
0.1
10
iI i'','
SO
100
,,
I ,,,,I
500
,
1ooo
(~tm)
Fig. 2. Classification of solids in terms of their behavior when suspended by gas.
Very large
exploding
bubbles. Low
bed voidage.
Liable to spout.
Circulating fluidized-bed reactors
2473
10
"k
10-1
10-2
10-3
10
102
102
10
(a)
Fig. 3. The generalized map of gas/solid contacting.
3.3. Behavior of suspended solids
Based on observations of many different solids,
Geldart (1973) and Geldart and Abrahamson (1978)
came up with the A, B, C, D classification of solids, as
shown in Fig. 2. BFB reactors usually use Geldart A,
AB, and B solids. C F B reactors can go down to
Geldart C solids.
3.4. Gas-solid contacting regimes
Grace (1986) prepared a graph to show the expected behavior of gas-solid systems all the way from
BFB to CFB. Figure 3 show a somewhat modified
version of his chart.
.Q
"o
o
t~
-r4. REACTORS WITH THROUGHFLOW OF SOLIDS, uo > u,
4.1. Observed behavior
Very fine solids have a very small terminal velocity
in gas, say mm/s, so with a reasonable throughflow of
gas, say in m/s, solids are continuously carried out of
the bed and have to be replaced with make-up solids.
Let this throughflow be at a bed mass velocity of
Gs (kg/m 2 s). This leads to suspended solids reactors of
various types depending on the ratio of u, to Uo.
Roughly, the solids distribution in the vessel is as
shown in Fig. 4.
0.2
0,4
0.6
Volume fraction solids: f
Fig. 4. Vertical distribution of solids in different contacting
regimes.
2474
"' [~"~::":::.":;':
[i.,;:.,...,,:.,:;.;.
I
D. Kunii and O. Levenspiel
i,i
"
"
"
I!l, """.." ~.~rl
:
l',:~..,'.'i-.~.,,:.l
, ~:-'
k ' ,;'~
~.~' i~'' ~
":..
V:,".~.~,|'~(~'t~,~
'iI
n~,~ ~;~'.',,~.","::.";:":i:]n~v~ I~~,...'.',~. n~,rs [, ~-~ !,~q n~,~r [:i
,::s,
~."
Il
...........
. ~i~
,~
,.,.~
'''.
.,.
.....
;,
.',
'.'.-.
' "
region
Dense
.'~:~i Ha,/'C /" :"'-"" ": "'] region
;:~
v,:
" " " -
l:.
.:
.': ....
"r
Fig. 5. Contacting models of vessels with throughflow of solids.
As observed, when the gas velocity of a bubbling
fluidized bed (BFB) is progressively raised, the bubbling action becomes very violent, bubbles coalesce and
become very large and finally expand to form a core
space in the dense region of the vessel. At the same
time the cloud and emulsion merge and retreat to the
walls of the vessel. In this state we have a fast fluidized
contactor (FF). Between the BFB and the F F regimes
we have a difficult to describe turbulent bed (TB).
At even higher velocities the wall region thins, dissolves, as the vessel enters the pneumatic conveying
regime (PC).
U"
r
"
.
~,
"
"
"
bubbling bed
fa = 0.4 ~ 0.6
turbulent bed
)ca = 0.2 ~ 0.4
fast fluidized bed
fa = 0.06 ~ 0.2
pneumatic transport
fa = 0.01 ~ 0.06.
upflowing clumps of denser material of mean
density Pz rising at velocity u2,
downflowing clumps of mean density p3 and
solids sliding down the wall at u~,
lean mixture with dispersed solids of mean density pa rising at ux.
"
ul
L'.~
".
h~
"
"
..;,
i!i:;i!
~-':
. ~
.~...
~.~
,~.
.,'
~i'": , ' ~ "-,. - ~ ; ~ ' , : , ~ ; ;:~# /
~
,',:.,.
,.,..,,.:., , :.
Dense
lower
region
,:: .:
L,~'~. r,'..:.'":.7" 4,~! .," .,'...'.."
,:..~ ,-:~.:.
I.
""
:: :. ,.- ...:
".'"'"
""
,':I
"- ~-'"-~"~ ........
Fig. 6. Mechanism which yields the exponential decay of
solids with height.
Figure 6 sketches the model. The differential equations representing this model are
Gs = G~I + G~2 + Gs3 = UlPl + u2pz - u3P3
dpl
dpl
d----t = ul ~
= K~(pz + P3)
This model also includes
dp2 = (K1 + K z ) p 2
-
interchange rates between clumps and lean mixture at Kx,
direction change from up to down at rate K2.
....
r~':, ;::"-;We visualize the upper lean region to consist of three
phases:
..---
(12)
.t,C~,
4.2. Contactin9 models
In CFB, solids are found throughout the vessels, in
a dense region in the lower part of the vessel, and
having a solid fraction fa, and a lean region above,
having a solid fractionf~ which decreases with height
as shown in Fig. 5.
In the lower dense region of the vessel the fraction
of solids is found to be
:~'-'.~
u 2 dz--~
dp3
- u3 ~ z r = K2P2 - K i p 3 .
2475
Circulating fluidized-bed reactors
'-
/:::F
.-
..I
.i
. . . .
"
:'...-'.,
.- ..,
...%.-
transport saturation
capacity
~--~'~x
upper lean
region, h
n~
[\
exit value
I ~
] X
z~
exponential decay clownto f*
f~ = f* +([d-f*)e "az~ ...(13)
High u o , small a value
" " \ ~ ' ~ ~- '"~'~
~ Low u o, large a value
.i
" ~
,--.n-
constant fraction
[V
of sohds
' rag,on,
"7
f* rex
fd
Volume fraction of solids, f
Uo
Fig. 7. The exponential decay model accounts for the distribution of solids in the vessel when the solids are
in throughflow.
Solving these equations Kunii and Levenspiel
(1991, p. 184) show the following
~< 0.01 for Oeldart B solids
fd decreases with Uo increases with dp
In a very tall (infinitely high) column the solid
fraction falls exponentially from fd to a lower
limiting value f*. The value f* represents the
transport carryin9 capacity of the gas, meaning
that all the clumps of denser mixture have 'dissolved' into the flowing gas.
In a shorter column the exiting solid fraction
fex is greater than the limiting f*.
given by eq. (12)
and for the lower dense region,
fcore -~ f *
fwall
(1
(17)
--
gwall
) (1
--
6).
Lean region. Since fox ~ 1 the solid fraction in the
We show this model in Fig. 7.
This distribution of solids is calculated from the
model, and does approximate what is observed So we
retain this model.
exit stream is
4.3. Material balance
We picture the CFB (in particular the FF contactor) as shown in Figs 5(d) and 7, having a one zone
upper lean region with a solids fraction f which decreases with height The lower dense region having
a constant solids fraction J] consists of two zones,
a lean core zone having solids f~.... and a dense wall
zone having solids fw,n.
From experiments reported in the literature [15
studies for eq. (15) and 38 studies for eqs (16) and (17)],
our best estimates today (see Kunii and Levenspiel,
1991, 1995) for needed physical quantities and relationships are
Equation (14), written for the exit location, gives
2-4 s-1
auo = ~5 s -1
for Geldart AB solids
(15)
L7 s -1
for Geldart B solids
f * is independent of Uo
~< 0.02 for Geldart A solids
Gs
f~x ~ ps(uo -- u,)"
(16)
(19)
fex = f * + (fd - f * ) e x p ( - a H t ) .
(20)
So the height of the lean region is
Ht = ! ln (19c~Z-~_
f'~ - f*
f , / .~
(21)
The height of the dense region is then
Ha = Ht -- HI.
(22)
At any point in the lean region the fraction of solids
is
f = f * + (fd--f*)exp(--azt).
for Geldart A solids
(18)
(14) or (23)
This gives the mean fraction of solids in the lean
region to be
J~=~
f d z , = f * +.fdaH,--f~x
(24)
2476
D. Kunii and O. Levenspiel
and the mean solid fraction in the vessel as a whole is
But for plug flow qph,g = 1, integration gives
f v , ~ _ f H ~ + f~Hd
H~
In Cao _fp,,gk'"Hp,,g
CAn
Uo
(25)
So the weight of catalyst in these regions is
= Atp, H t f
Wd = A,p~H~f~
(26)
~=~+w~.
These equations relate u0, G~, and W. Thus, if you are
given a specific gas velocity and you want a specific
weight of solids in the reactor, these equations will tell
you what solid feed rate to use.
- r'j' = k'"CA
or
CAo
In Ca~
r-
1
lfb k'" +
1
1
P q
6BFBKbc
- r~4 = k'CA. (27)
In Fig. 8 let us sketch the contacting as viewed by our
proposed reactor models for the various contacting
regimes,
5.1. Bed at minimum fluidization
For a bed at m i n i m u m fluidization we assume plug
flow [see Fig. 8(a)]. We then write for any level in the
reactor
f~k'" +
1
(1/fiaFRK~) + 1/f~k'"
x -HBFB
u0
(30)
For a dense bed of height HnFB and containing the
same a m o u n t of solids but with solids uniformly distributed throughout, and with plug flow of gas, we
have
dCa
-- U o ~ -z =folugqr, h, gk'"Ca.
(fb + f + f)k'"HBFB
CAO
In - -
(28)
C.4ex
12o
Fixed
BFB
FF
CAex=CAd
CAex< CAd
CAex<C.Ad
~IBFB
CAex
!i :i:
Lean
--Tnd
C ~df
;..~ .ow
,...,,.r
,:': ~.'....,
%';
~
~t,I
muo
./~
Dense
;*7
. '~--';::~t.
:,
. .-.....:.
" r.~&-:
~~
(b)
i,t:.
t, :,~... "','
"."
i:~
CAO
CAO
(a)
(29)
5.2. Bubblin9 fluidized bed
For a bubbling fluidized bed with its three connected zones [-see Fig. 8(b)], we assume bubble upflow
but no flow through clouds and emulsion, only interchange Kbc and K~e. So, as shown by Kunii and
Levenspiel (1991, p. 291), the integrated conversion
expression is
5. PERFORMANCEAND CONTACTEFFICIENCY OF CFB
REACTORS
Consider a first-order solid-catalyzed reaction
A ~ R,
k'W,
uoAt
(c)
Fig. 8. Proposed reactor models for the various contacting regimes.
(d)
(31 )
Circulating fluidized-bed reactors
The contacting efficiency is then given by eq. (30)
divided by eq. (31), or
k'"
OBFBKbc q f +
1
(k'"/JBwKc~ + 1/f~
2477
ally from qd at zt = 0 to q~ = 1 very high up in the bed.
W i t h f and rh replaced by eqs (23) and (37), integration
of eq. (36) gives
--=
H~-
CAex
(1-e
+ k '" (f~ - f * ) F_I -- e - aHt
UO
q~v~, -
bn~)
UO
A + +L
(32)
5.3. Fast fluidization
For fast fluidization we have distinctly different
expressions for conversion in the upper lean and
lower dense regions; see Fig. 8(c).
For the dense region the three zones of a BFB
collapse into two. Thus,
fb ~ f ~ ....
I fore k"' +
(1/3a.vvKcw) 1/fwa,,k'"
I H d , FF
Uo
(33)
The various constants appearing in this expression
are not well known today; however, we estimate them
as follows:
fco~e ~ f *
In CAd
C A ex
k ' " f * t l ~ + k " ( f d - - f * ) ( 1 -- e ~n,). (39)
U0
U0 a
Finally, the overall conversion is given by eq. (38) or
(39) with eq. (33):
fi +fe ~fwall, Kc~ ~ m and
Kbc ~ Kcw.
So the performance expression of the BFB reduces, for
the dense region of F F reactors, to
inCa a/o:
and for the special case where r/d = 1 (very lean solids,
high gas velocity, pneumatic flow)
CAex
CAO
CAex
CAd
CAd
CAO
r/d = rh = 1
0.6 0.9
(34)
Kcw = 5-20 s- 1
ewall ~ ~mf ~--- 0.5-0.6.
The contact efficiency compared to plug flow, by an
argument similar to the BFB, is
o~ q (k'"/3e, vvK~w) + 1/fwa~
~d, FF =
fcore fwall
(35)
For the lean region of the F F reactor of height
Ht and with decreasing solid fraction with height, the
rate of any level is
dC,4
- Uo ~ z = qtfk'"CA.
(36)
The t e r m f is given by eq. (23). The efficiency term
for FF, TB, and above the BFB are not well known
today. The only studies to date (Furusaki et al., 1976;
see also Kunii and Levenspiel, 1991, p. 283) show that
the contacting efficiency in the lean phase can be
represented by an exponential expression
rh = 1 - (1
--
rld)e -bz'
(37)
where b = 6.62 m - ~, by experiment. This expression
shows that the lean region efficiency rises exponenti-
CAO"
(40)
(41)
and the regular plug flow expression, eq. (29), applies:
(faHa + f H 3 k ' "
CAe x
6d. v r =
For pneumatic conveying all the particles are evenly
dispersed in the gas. This makes contacting ideal or
close to ideal. So, if gas passes in plug flow up the
reactor the contact efficiency
Cao
In - -
~- 0.01
CAex
X A . . . . . 11
and
(42)
bl0
The reliability of the above conversion predictions all
depend on the reasonableness of the models and on
the values of the parameters chosen.
f*
b.
Finally, probably the most important value of models such as presented here is to be a guide for experimenters, and to suggest areas where research is
needed.
6. DESIGN PROBLEMS
TO illustrate the use of the material presented here,
consider a particular process requirement and let us
explore various design alternatives to meet this requirement. For this let us suppose that we want to
find the fractional conversion of A to R. We are told
that our catalyst is very active with first-order
kinetics:
A ~ R,
1 dNa
- - - k'CA,
W dt
k' = 0.01 ma/kg cats
k'" = 10 m3/m 3 cat s.
Figure 9 shows the four design alternatives that we
will consider. We refer to this figure as we proceed.
2478
D. Kunii and O. Levenspiel
Example I
Example 2
Pneumatic Transport
(low solid flow)
(high solid flow)
z,m
Example 4
5
Bubbling
fluidized
Fast fluidized Turbulent
~.?'.i
2
dense-
.z~ .'!
;::.';!
-.,.-
uo
.~.";,"':;~:
::,'-v
~" :.:;
.... "
N.
Kglm2"s
= 6 m/s
6 mls
1.5 m/s
0.375 m/s
(a)
(b)
(c)
(d)
Fig. 9. Four design alternatives for reactors treating a given flow rate of reactant gas.
Data. We want to treat ~ = 4 0 m o l / s of pure A
feed at 300 K a n d 1.3 a t m u n d e r which condition the
properties of gas a n d solid (porous catalyst particles)
are
h = 40 mol/s
Po = 1.2 kg/m 3
Gas
This high gas velocity a n d small particle size suggests t h a t solids m a y be rapidly carried out of the
reactor, so we m a y need to use solid circulation. So let
us tentatively choose a solid circulation system with
a reactor of height H, = 6 m a n d with a solid circulation rate of G, = 100 kg/m 2 s. This design is s h o w n as
Fig. 9(a).
p = 1.8 x 10 -5 k g / m s
D = 4 x 10-5 m2/s
in the bed
(1) F r o m Fig. 2 we see that we have here a G e l d a r t
A solid.
(2) F r o m eqs (7) a n d (8)
shape is spherical
Porous
catlyst
solid:
Solution:
p~ = 1000 kg/(m cat) 3
d [Po(Ps
Deft = 1 0 - 7 m3/m cat s,
a.
Pg)gl
e/
1/3
in the p o r o u s particle.
=
55
F r o m this d a t a the volumetric flow rate of feed gas is
Vo = uoA .
hRT
.
p
40(8.314)300
.
.
1.3(101 325)
0.7574 m3/s.
6.1. Example 1. Pneumatic transport reactor
Because the catalyst is said to be very active, let us
choose very fine particles, d e = 55 ~tm, a n d a vessel
diameter of dt = 0.4 m. The superficial gas velocity for
this vessel is then
Vo
0.7574
Uo - A - n(0.4) 2 = 6.0 m/s.
10-6(1'2(1
--L2)9'8 ] '/3
(1.8 10-5)2
= 1.82
u~
U0
[
= 6
--
It(ps -- Po)g
(1'2)2
11/3
1.8 X 10-5(1000 -- 1.2) 9.8
= 12.08.
2479
Circulating fluidized-bed reactors
(10) Estimate the contact efficiencies in the two
regimes:
For the dense region, from eq. (41), r/d = 1.
For the lean region, from eq. (41), q~ = 1.
(11) Since we have plug flow in both regions we can
use eq. (42) instead of eq. (33) for the dense
region and eq. (38) for the lean region. Incidentally, you would get the same answer with
these more complicated equations. Thus, we
have
(3) The terminal velocity, from eqs (10) and (8),
u* =
Ut =
I 18
0.591 ] - 1
(~p,)2 + (d,)l/~-2j = 0.1703
u, [~(p, - p.)0] '/3
t
[_
p.
= 0.1703 -
(1.2)2
= 0.084545 m/s.
CAo
ln--
This says that at 8 cm/s solids are blown out of
the bed. Note, we plan to use a velocity of 6 m/s.
(4) Figure 3 shows that we are in the Pneumatic
transport regime, so we guessed right when we
chose the vessel of Fig. 9(a).
(5) Evaluate the constants, from eqs (12), (15) and
(16):
j~ = 0.06
f * = 0.01
u o a = 3 s -1 or a = 0 . 5 m -1
(6) F r o m eq. (19)
f~x--
Gs
100
ps(uo - u~)
1000(6 -- 0.084545)
= 0.0169.
CAex
I-(0.06)(2.0390) + (0.0318)(3.9610)] 10
= 0.4138.
Therefore,
CAex = 0.6611
CAO
or
X a , total
34%.
6.2. Example 2. Pneumatic transport reactor
The design of Example 1 is unacceptable because
conversion is too low. Let us see if we can raise it
by doubling the solid throughflow rate to
Gs = 200 kg/m2s while keeping all else unchanged;
see Fig. 9(b).
(7) Locate the dense and lean regions, from eq. (21),
Solution: Following the procedure of Example 1
we find
1 t/fa - - f * "~
Ht = - In 1 7 - - - ~ l
\fo~-f /
1 l n ( 0.06--0.01 ']
=0.~
\0.~9--~1
]
pneumatic transport with r / = 1,
Ha 4.5 m, instead of 2.0 m,
Wtotal = 42.50 kg, instead of 31.20 kg.
=
3.9610 m
Ha = 6 - 3.9610 = 2.0390 m.
so
XA
. . . . . . 11 =
42%.
(8) Find f in the lean region, from eq. (24),
f =/,
+fd --f~x _ 0.01 +
aH~
O.O6 0,0169
0.5(3.9610)
= 0.0318.
(9) The weight of catalyst in the vessel, from
eq. (26),
Wa = Atp,Hafa
7~
This is better than that found in Example 1 but still
not good enough. So let us try still other alternatives.
6.3. Example 3. Turbulent or fast fluidized reactor
Maybe the chosen gas velocity was too high in
Examples 1 and 2. Let us try a lower gas velocity,
Uo = 1.5 m/s, in a squatter vessel, dt = 0.8 m and
H, = 3 m, while keeping the throughflow rates of gas
and solids unchanged from Example 1; see Fig. 9(c).
Also take
= ~ (0.4) (1000)(2.0390)(0.06)
fa = 0.16
= 15.3737 kg
Wl = Atp~Hjf
m3
of solids in the lower dense region
m 3 of dense region
m 3 of solids
f * = 0.01 m3 of vessel'
during pneumatic
conveying
= ~ (0.4) z (1000)(3.9610)(0.0318)
4
= 15.8286 kg
Therefore Wtotal = 31.2 kg.
6 = 0.7
m 3 of core
m 3 of dense region
Kcw = 5 s- x.
D. Kunii and O. Levenspiel
2480
Solution:
Following the procedures of Example 1 we
find
we have a turbulent or fast fluidized bed with
qd = 0.72,
H d = 1.5m,
Wtot.l = 165 kg,
so
XA . . . . . . 11 =
82%.
double the bed diameter and lower the gas velocity
accordingly so as to keep the gas flow unchanged; see
Fig. 9(d).
Calculate the conversion for these conditions.
Solution: We first find that we are right in the middle
of the BFB regime with negligible carryover of solids,
Uo < u,. N o throughflow of solids is needed, so let us
use a BFB with the following properties:
dt = 1.6 m,
Bubblingfluidized reactor
6.4. Example 4.
It just occurred to us that we may have chosen
particles smaller than was needed in Examples 1-3.
We should have checked for this with the Thiele
modulus, but did not. Let us do this now. What we
find, see Levenspiel (1996b), is
dp /k'"
55x10
M r = 6 ~ / D e ft
-6
6 /-~10
- / ~
M~ --
/~10
X/~
0.8 m ,
= 0.366.
This size of solid is still free from internal resistance to
diffusion. So let us adopt these solids and again
31
42
165
643
f.:.,
#"
::"
~:,:
:-
'!. .
,
%':,
=,
q = 11%
11 =
61%.
Contact Conversion
efficiency
XA
~l = 1.00
TI = 1.00
TId = 0.72
r] = 0.11
34/.
43%
82*/,
61%
(i"
Q
~. ,-
10 . 6 m
6.5. Comments about these four examples
We have explored various design alternatives and
we have shown how to determine the performance
expected for each, as summarized in Fig. 10.
Some designs give very low conversions, others
high. Some suspend a large inventory of solids (643 kg
.',,'
"....
UO = 0.375 m / s .
with
Xa . . . . . .
1
2
3
4
~.i
dp = 220
SO
~ = 0.4,
Wt=643kg
.,
spherical particles
In the fluidized bed D ---4 x 10 -5 m3/m cats.
For this system we find, by the methods outlined in
Levenspiel (1996a),
Example
number W, kg
(9
0.16 m,
fb = 0.003,
era: = 0.6,
= 0.092.
This is a much smaller Thiele modulus than needed,
so we can safely use larger particles Let us try
dp = 220 x 10 - 6 m, For these the Thiele modulus is
220 x 10
Hmf =
db =
'
~ s.'.,-
?:.';:?
.'....
:3!..' i solid
gas
\
''
solid
solid
as
gas
v
~:-t.
ttg..
~,. i
e"
(~ ...~g
Pneumatic
Turbulent
v
Solids throughflow
Bubbling
I
No solids throughflow
Fig. 10. A wide range of conversions result from the four design alternatives. Fast fluidization turns out to
be best.
2481
Circulating fluidized-bed reactors
in Example 4), others a small inventory (31 kg in
Example 1). Consequently, the pressure drop and
pumping requirement can differ greatly for these designs.
Although the BFB has a very low contacting efficiency (11%), the amount of solid used there is so
much greater than in the very much leaner solid
throughflow systems (see Examples 1-3); hence, the
intermediate conversion in the bubbling fluidized bed.
In these examples the F F turns out best, however
the conclusions found here hold only for the conditions examined in the above examples. In other
situations the advantage may be with a different contacting regime. So, in all cases it would be wise to
examine a wide range of possible operating conditions
and particle sizes,
Finally, we should point out that we have only
presented what may be called '$10 solutions' here.
Certainly, one can build on this approach.
Kcw
K1, K2
Mr
h
r~
r~'
Uo
ut
u*
V0
At
b
BFB
Ca
CFB
dp, d.... dsph
d*
Deft
fi,(i=b,c,e,
d, l, etc.)
f,
g
Gs
Ha, Ht, Htotal
k'
k'"
Kbc, Kce, Kbe
NOTATION
decay constant for solid fraction in
the lean region of a CFB; see
Fig. 6, m cross-sectional area of a CFB reactor, m 2
decay constant for gas/solid contact inefficiency in the lean region
of a CFB; see eq. (23), m bubbling fluidized bed
concentration of reactant A,
mol/m 3
circulating fluidized bed
measures of particle diameter; see
eqs (1)-(6), m
dimensionless measure of particle
diameter; see eq. (7)
diffusion coefficient of gas in the
porous catalyst, m 3 gas/m solid s
volume of solids in section i of
a slice of bed/volume of that slice
of bed, m 3 solid/m 3 reactor
solid fraction in gas stream in
pneumatic conveying conditions
acceleration due
to
gravity
(= 9.8 m/s 2)
mass velocity of solids through the
CFB, kg/m 2 of bed s
height of regions and of the
CFB, m
first-order reaction rate constant
based on unit mass of catalyst
solid, m 3 gas/kg cat s
first-order reaction rate constant
based on unit volume of catalyst
solid, m 3 gas/m 3 cat s
gas exchange coefficients between
bubble, cloud, and emulsion,
m 3 gas moving from one zone to
the other/m 3 of bubbles s
Xa
Z
gas interchange coefficient between core and wall region m 3 gas
moving from one zone to the
other/m 3 of core s
interchange coefficients; see above
eq. (13), s -1
Thiele modulus, defined in
Example 4, dimensionless
molar feed rate, mol/s
reaction rate based on mass of
catalyst, mol/kg s
reaction rate based on volume of
catalyst, mol/m 3 cat.s
superficial gas velocity in the reactor m3/m 2 vessel s
terminal velocity of a particle falling through the gas, m/s
dimensionless gas velocity; see
eq. (8)
volumetric feed rate, ma/s
weight of catalyst, kg
fraction of reactant A converted,
dimensionless
height in the lean region of a
CFB, m
Greek letters
wake volume/bubble volume, dimensionless
volume fraction of lean region, bubble or core,
in a section of reactor, dimensionless
ewan void fraction in the wall zone, dimensionless
p
viscosity of gas, kg/m s
q
efficiency of reactor when compared to plug
flow, kg/kg
p
density, kg/m 3
~bs
sphericity of solid particles; see eq. (2), dimensionless
Subscripts
b
c
core
d
f
g
l
bubble
cloud
in the core region of a CFB
dense region of a F F contactor
fluidized condition
gas phase
lean region of a CFB
mf
at minimum fluidizing conditions
s
solid
wall sliding down the wall of a CFB
REFERENCES
Ergun, S. (1952) Fluid flow through packed columns.
Chem. Engng Prog. 48, 89.
Furusaki, S., Kikuchi, T. and Miyauchi, T. (1976)
Axial distribution of reactivity inside a fluid bed
contactor. A.LCh.E.J. 22, 354-361.
Geldart, D. (1973) Types of gas fluidization. Powder
Technol. 7, 285-292.
Geldart, D. and Abrahamsen, A. R. (1978) Homogeneous fluidization of fine powders using various
gases and pressures. Powder Technol. 19, 133-136.
2482
D. Kunii and O. Levenspiel
Grace, J. R. (1986) Contacting modes and behavior
classification of gas-solid and other two-phase suspensions. Can. J. Chem. Engng 64, 353-363.
Haider, A. and Levenspiel, O. (1989) Drag coefficient
and terminal velocity of spherical and nonspherical
particles. Powder Technol. 58, 63-70.
Kunii, D. and Levenspiel, O. (1991) Fluidization
Engineering, 2nd Edn. Butterworth-Heinemann,
Boston, MA, U.S.A.
Kunii, D. and Levenspiel, O. (1995) The vertical distribution of solids in circulating fluidized beds. In
Fluidization VIII, Tours, pp. 17-24.
Levenspiel, O. (1996a) The Chemical Reactor Omnibook, Chap. 25. OSU Bookstores, Corvallis, OR,
U.S.A.
Levenspiel, O. (1996b) The Chemical Reactor Omnibook, pp. 22.3 or 23.1. OSU Bookstores, Corvallis,
OR, U.S.A.