
Yang Fluidization, Solids Handling, and Processing
.pdf
Fluidized Bed Scale-up |
35 |
at y = 0
u ′ = i/ above distribution holes
Eq. (27)
u ′ = 0 elsewhere
At all internal surfaces |
x′s |
, y′ |
|||
|
|
|
|
s |
|
Eq. (28) |
|
′=0 , |
υ |
′ =0 |
|
u |
|||||
|
|
|
|
N |
|
Eq. (29) |
|
|
|
|
|
The term Po |
(ρ f |
uo2 ) can be ignored when the fluid velocity is |
small compared to sonic velocity or the absolute pressure does not change enough to influence the thermodynamic properties of the fluid; it will be ignored in this development. Note that the fluid pressure level still influences the fluid density.
From the continuity and momentum equations for the fluid and solid phases along with the boundary conditions, the following groups of independent dimensionless parameters are found to control the hydrodynamics, noting our assumption that the particle-particle forces are only dependent on hydrodynamic parameters,
Eq. (30) |
uo2 |
, |
ρs |
, |
β L |
, |
Gs |
|
gL |
ρf |
ρs uo |
ρs uo |
|||||
|
|
|
|
as well as bed geometry.
The dimensionless fluid pressure is not included since it is a dependent parameter.
4.3Fluid-Solid Forces
The drag coefficient β can be expressed in several different limiting forms depending on the flow conditions. At low voidages typical

36 Fluidization, Solids Handling, and Processing
of bubbling beds, the Ergun relationship or similar form can be used. In that case β can be expressed as,
Eq. (31)
β Lϕ |
|
|
ε (1 - ε )2 |
|
μ |
|
|
|
(1 - ε ) |
|
- |
|
|
|
|
2 |
|
Lρ f |
s |
|
|
|
L |
|
u′ |
υ′ |
|
ε |
|||||||||
|
= 150 |
|
|
|
|
|
+ 1.75 |
ÿ |
|
|
ÿ |
|
|
|
|
|||
ρs uo |
ε 3 |
|
ρs uoϕ d p |
|
d p |
ε 3 |
|
d p |
|
ρs |
||||||||
|
|
|
|
|
|
In the limit of very high voidage, the drag coefficient can be related to the single particle drag coefficient. For the case of spherical particles,
|
β L |
|
3 |
|
|
|
|
|
f (ε ) |
ρ f |
|
L |
|
= |
|
u′-υ′ |
|||||||||
|
|
|
ρs |
|
|
|||||||
Eq. (32) |
ρs uo 4 C Dÿ |
|
|
ÿ |
|
|
d p |
where the drag coefficient, CD, in turn can be expressed in the form of
Eq. (33) C D = f ÿ ρf
uo d p
μ
ÿ
In the more general case, CD will also be a function of particle shape, sphericity, surface roughness and turbulence intensity of the fluid.
The dimensionless drag coefficient β L(ρs uo) in Eq. (30) can be expressed in terms of other fluid parameters by the use of Eqs. (31–33). For low voidages where the Ergun-like expression holds,
|
β L |
|
|
ρs uo d 2p ϕ 2 |
|
L |
ρf |
Eq. (34) |
|
= f |
ÿ |
|
, |
|
|
|
ρs uo |
|
|
μ L |
dp φ |
ρs |
ÿ
In this expression dp represents the mean diameter when a distribution of different size particles are in the bed.

Fluidized Bed Scale-up |
37 |
At high voidage, using Eqs. (32) and (33),
Eq. (35)
β L |
|
ρ f uo d p |
|
ρ f |
|
L |
|
|
= |
fÿ |
|
, |
|
|
|
|
μ |
ρs |
|
|
|||
ρs uo |
|
|
|
d p |
, ϕ
or shape, roughness, fluid turbulence
ÿ
Traditionally, the mean diameter is defined as the surface area averaged mean. Although this mean may be appropriate for flow resistance primarily due to surface shear forces, it is not the proper choice for drag which prevails at higher particle Reynolds numbers (nor is it the obvious choice for a mean diameter to use for bed-to-surface heat or mass transfer). It is more general to include along with the mean particle diameter, the particle size distribution, the particle size nondimensionalized with respect to the mean diameter, and the particle sphericity. To be more exact, the particle aspect ratio and surface roughness should be included at high particle Reynolds numbers. By use of isotropic material with common roughness levels, these last two parameters can be overlooked.
Using Eq. (34), the set of independent dimensionless parameters (Eq. 30) becomes,
Eq. (36)
uo2 ρ s |
|
ρs uo d 2p ϕ 2 ρf uo L |
|
Gs |
|
||||
|
, |
|
, |
|
, |
|
, |
|
, |
gL |
ρ f |
μ L |
μ |
ρs uo |
L , Bed geometry, ϕ , particle size distribution (PSD)
D
These can be rearranged by combination of parameters. It must be borne in mind that such manipulation by itself does not lead to any decrease in the number of dimensionless parameters. One such modification is,
Eq. (37)
2 |
|
ρ ρ |
uo d p |
|
ρf uo L |
|
Gs |
|
L |
|||
uo |
, |
s |
, |
s |
|
, |
|
, |
, |
|
, Bed geometry,ϕ ,PSD |
|
|
|
|
μ |
μ |
|
|
||||||
gL |
ρ f |
|
|
ρs uo |
D |

38 Fluidization, Solids Handling, and Processing
In this form, uo2gL , the Froude number, can be viewed as a ratio of
inertial to gravity forces; |
ρs ρf |
is a ratio of particle to fluid inertial |
forces; ρs uo d p μ is the Reynolds number or ratio of particle inertial to |
||
fluid viscous forces; and ρ f |
uo L |
μ , a Reynolds number based on the bed |
dimensions and fluid density, is a ratio of fluid inertial to viscous forces. Another common form is obtained by combining the Froude and
Reynolds number to obtain the Archimedes number, which omits uo ,
Eq. (38)
ρ f ρs d 3p g |
ρ f uo dp |
|
|
=ÿ |
|
μ 2 |
μ f |
ÿ
3
gL |
μ f |
|
|
|
ÿ |
|
ÿÿ |
|
|
||
uo2 |
|
ρ f uo L |
ρ s |
|
||
ÿ |
|||
ρ |
f |
||
|
|
The list of dimensionless parameters can be rewritten as,
Eq. (39)
3 |
|
ρ |
|
|
|
|
ρ f uo L |
|
|
|
L |
|
ρ f ρs d p g |
|
|
|
2 |
|
|
Gs |
|
||||
|
, |
|
s |
, |
uo |
, |
|
, |
, |
|
, Bed geometry, ϕ , PSD |
|
μ2 |
|
|
|
μ |
|
|
||||||
|
ρ f |
gL |
|
ρs uo |
D |
Note that there isn’t anything more “fundamental” about one form compared to the others. Each has the same number of dimensionless groups which are made up of independent parameters which can be set by the bed design and operation and the choice of particles. However, when the number of dimensionless groups is simplified by omitting some phenomena, the reduction in number of groups could be influenced by the form chosen.
This set, Eqs. (36), (37) or (39), will be referred to as the full set of scaling relationships.
4.4Spouting and Slugging Beds
He, Lim and Grace (1995) have shown that for spouting beds, cohesive factors are important and this group must be augmented by including the internal friction angle and the loose packed voidage to achieve similar scale models. Since interparticle friction can occur in slugging beds, these additional parameters should be included to properly scale slugging beds.
Fluidized Bed Scale-up |
39 |
5.0SIMPLIFIED SCALING RELATIONSHIPS
In later sections, the use of the scaling relationships to design small scale models will be illustrated. For scaling to hold, all of the dimensionless parameters given in Eqs. (36), (37) or (39) must be identical in the scale model and the commercial bed under study. If the small scale model is fluidized with air at ambient conditions, then the fluid density and viscosity are fixed and it will be shown there is only one unique modeling condition which will allow complete similarity. In some cases this requires a model which is too large and unwieldy to simulate a large commercial bed.
To allow more flexibility, we will explore means of simplifying the scaling relationships by reducing the number of parameters which must be maintained constant. In most situations, one would expect that not all of the parameters are of first order importance. By reducing the number of parameters which must be maintained in the model, it may be possible to model larger commercial beds with small scale models. We will look at simplifications of the interparticle drag at the extreme of small and large Reynolds numbers based on particle diameter. This will span the range from low velocity FCC bubbling beds to high velocity fast beds or pneumatic transport using large particles. If the same simplification can be shown to hold in both of these limits, it is reasonable to consider application of the simplification over the entire range of conditions.
5.1Low Reynolds Number
At low particle Reynolds numbers for a bubbling bed, the Ergun expression can be simplified using only the first term in Eq. (31).
Thus,
Eq. (40) |
β L |
|
→ |
|
150 |
(1 - ε )2 |
|
μ L |
|
|
for |
ρ f |
uo d p |
|
→ 0 |
|||||||||
ρ |
uo |
|
ε |
2 |
|
ρ |
(φ |
dp |
)2 |
|
μ |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
s |
|
|
|
|
|
|
|
|
|
|
|
s uo |
|
|
|
|
|
|
|
|
|
|
|
At the same limit, the minimum fluidization velocity can be |
||||||||||||||||||||||||
written as, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Eq. (41) |
p |
= ( |
|
- |
ρ f |
)g (1 - |
ε mf |
) = |
150 |
(1- ε mf )2 |
|
μ umf |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
L |
|
|
ρs |
|
|
|
|
|
|
|
|
|
ε 3mf |
|
(φ dp )2 |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
40 Fluidization, Solids Handling, and Processing
For gas fluidized beds where ρs - ρf can be replaced by ρs,
umf |
= |
|
ρs g(1 - |
ε mf ) |
|
|
|||
|
2 |
|
|
|
|
|
ù |
||
Eq. (42) |
é |
|
|
|
μ |
|
|||
ê150 |
(1 - ε mf ) |
|
|
|
|
ú |
|||
|
ε 3 |
|
|
(φ |
|
)2 |
|||
|
ê |
|
|
dp |
ú |
||||
|
ë |
mf |
|
|
|
|
û |
Substituting Eq. (42) into Eq. (40),
Eq. (43) |
β L |
= |
|
|
g (1 - ε )2 |
|
|
ε mf3 L |
|||||
ρ |
umf |
|
(1 - |
) |
2 |
|
ρ |
||||||
|
s uo |
|
|
|
|
ε mf ε |
|
|
s uo |
||||
and |
|
|
|
|
|
|
|
|
|
|
|
|
|
Eq. (44) |
β L |
Fr |
= |
|
uo (1 - ε )2 ε mf3 |
||||||||
|
|
|
|
|
|
|
|
|
|||||
ρs uo |
|
umf ε 2 (1 - ε mf ) |
|||||||||||
|
|
|
|
|
|
Thus, in the low particle Reynolds number limit, maintaining uo /umf , εmf , and Fr identical between two fluidized bed guarantees that β L /ρsuo is also identical. Although φ and dp are eliminated between Eqs. (40) and (42), in general, particle sphericity and dimensionless size distribution should be held constant in the scaling since they influence εmf . The use of φ and a mean diameter in the Ergun expression only approximates the effects of these parameters. Note if the two models display identical dynamic characteristics, then ε is a dependent variable whose distribution throughout the bed should be identical for both fluidized beds. In this limit, the governing parameters given in Eq. (36) can be expressed as
Eq. (45) |
uo2 |
, |
ρs |
, |
uo |
, |
L1 |
, |
G s |
bed geometry, φ , PSD |
|
|
|
L2 |
|
||||||
|
gL ρf |
umf |
|
ρs uo |
where εmf will be a function of particle sphericity and size distribution.
Fluidized Bed Scale-up |
41 |
This will be referred to as the simplified scaling relationship. At low Reynolds numbers, this still includes the gas-to-particle density ratio.
5.2High Reynolds Numbers
Consider the limit of high particle Reynolds numbers where the inertial term in the Ergun equations dominates.
In this limit,
|
β L |
|
→ 1.75 |
(1 - ε ) ρf |
|
L | u′ - v′ | |
ρ f uo d p |
>> 1 |
|||||
Eq. (46) |
|
|
|
|
|
|
|
|
|
for |
|
||
ρ |
|
u |
|
ε |
ρ |
|
φ d |
|
μ |
||||
|
s |
o |
|
p |
|
||||||||
|
|
|
|
s |
|
|
|
|
where u´ = u/uo and v´ = v/uo and u´- v´ is a dimensionless slip velocity. The minimum fluidization velocity can be expressed as,
Eq. (47) |
p |
= (ρs - ρ f )g (1 - ε mf ) ~ |
1.75 (1 - ε mf ) |
ρ f |
umf2 |
|
|
ε 3mf |
|
φ d p |
|||
|
L |
|
rearranging and using ρs in place of ρs - ρf ,
Eq. (48) |
1.75(1 - ε mf ) |
ρ f |
= |
g (1 - ε mf ) |
||
ε mf3 |
|
ρs φ dp |
umf2 |
|||
|
|
Substituting this into Eq. (46) and multiplying by Fr,
Eq. (49)
Fr |
β L |
= |
uo2 |
|
ε mf3 L |u′ - v′| g(1 - ε ) |
= |
uo2 |
|
ε mf3 (1 - ε )|u′ - v′ | |
|
ρs uo |
gL |
|
ε umf2 |
|
umf2 |
|
ε |
|||
|
|
|
|
|
At large particle Reynolds numbers, just as at low Reynolds numbers, the dimensionless drag, βL/ρsuo, is identical when uo /umf , εmf
42 Fluidization, Solids Handling, and Processing
and Fr are identical. The variables ε, u´ and v´ are dependent dimensionless variables which are identical for two similar fluidized beds. In this limit the same set of governing dimensionless parameters applies as in the low Reynolds number limit, given by Eq. (45).
5.3Low Slip Velocity
Finally, consider the case when the magnitude of the slip velocity between the particles and the gas is close to umf /ε everywhere in the fluidized bed. With the vertical pressure drop equal to the particle weight, the following holds for any value of the particle Reynolds number,
Eq. (50) |
p |
= β | u - v | ≈ β |
umf |
= |
ρs g(1 - ε ) |
||||||||
L |
|
||||||||||||
|
|
|
|
|
|
ε |
|
|
|
|
|||
Eq. (51) |
β |
= |
ρs |
g(1 - ε ) |
ε |
|
|
|
|
|
|||
umf |
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|||||
and |
|
|
|
|
|
|
|
|
|
|
|
|
|
Eq. (52) |
β L |
Fr |
= uo / umf (1 - ε )ε |
for u′ - v′ → |
1 |
|
umf |
||||||
ρs uo |
ε uo |
||||||||||||
|
|
|
|
|
|
|
|
Again, when uo /umf and Fr are identical for two beds and the slip velocity is close to umf /ε, the dimensionless drag coefficient is also identical for two beds.
For all three limiting cases identified above, similitude can be obtained by maintaining constant values for the dimensionless parameters,
Eq. (53) |
uo2 |
, |
ρs |
, |
uo |
, |
L1 |
, |
Gs |
, bed geometry, φs , PSD |
|
|
|
L2 |
|
||||||
|
gL ρf |
umf |
|
ρs uo |
Fluidized Bed Scale-up |
43 |
5.4General Case
Since the same simplified set of dimensionless parameters holds exactly at both high and low Reynolds numbers, it is reasonable to expect that it will hold, at least approximately, over the entire range of conditions for which the drag coefficient can be determined by the Ergun equation or an equation of similar form.
The validity of the simplified parameters can be checked numerically for the intermediate range of values.
From Eq. (31),
Eq. (54)
β L |
= |
150(1 - ε )2 |
|
μ L |
|
+ |
1.75(1 - ε ) ρf |
L | u′ - v′| |
|||||||||||
ρ |
s |
uo |
|
ε 2 |
|
ρ |
(φ |
d p |
)2 |
|
ε |
|
ρ |
s |
|
φ d p |
|
||
|
|
|
|
|
|
s uo |
|
|
|
|
|
|
|
|
|
|
where u´- v´ can also be represented as uslip /uo. Using the Ergun equation to determine umf ,
Eq. (55)
|
p |
= ρs g (1 - ε mf ) |
= |
150 |
(1- ε mf )2 |
|
μ umf |
|
|
+ |
1.75(1- ε mf ) ρf |
umf2 |
|
|
|||||||||||||||||||||||
|
L |
|
ε mf3 |
|
|
(φ d p )2 |
|
|
|
ε mf3 |
|
|
|
|
|
|
φ d p |
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
Dividing Eq. (54) by 55, and rearranging, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
é |
|
1.75 |
|
φs |
Re |
|
|
|
|
|
|
|
|
ù |
||||||
|
|
|
æ |
|
|
ö |
|
Fr(1 - |
|
) |
2 |
|
|
|
ê1 + |
|
|
|
|
|
|
|
|
|
ε | u¢ - v¢|ú |
||||||||||||
|
|
|
β L |
|
ε mf |
|
|
|
150 (1 - ε ) |
|
|||||||||||||||||||||||||||
|
|
|
ç |
÷ |
|
|
|
|
ε |
|
|
ê |
|
|
|
|
|
|
|
|
|
|
ú |
||||||||||||||
|
|
|
ç |
ρ |
|
÷ |
æ |
ö |
|
|
|
|
|
= |
ê |
|
1.75 |
|
|
φ |
s |
Re |
|
|
umf |
ú |
|||||||||||
Eq. (56) |
è |
|
s uo ø |
ç |
uo ÷ |
|
2 |
3 |
|
|
|
|
1+ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
ç |
|
÷(1 - ε ) |
ε mf |
|
|
ê |
|
150 |
|
(1- ε |
|
|
|
) |
u |
|
|
ú |
||||||||||||
|
|
|
|
|
|
|
è umf ø |
|
|
|
|
|
|
ë |
|
|
|
|
|
|
|
|
|
mf |
|
|
|
o |
û |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ê |
|
|
|
|
|
|
|
|
|
|
|
|
ú |
where Re = ρf uo dp /μ .
44 Fluidization, Solids Handling, and Processing
It is easy to verify the three limits defined previously by use of Eq. (56).
For the more general case, Fig. 21 shows the value of β given by Eq. (56) relative to β at low Re over a range of conditions when uo /umf is 10 and 3, respectively, and Fr and φs remain constant. When uo /umf and the slip velocity are high, there is a larger variation of dimensionless drag coefficient with Reynolds number. Note that β does not vary with particle Reynolds number when the Reynolds number remains above about 103 or below about 10. Figure 22 illustrates the results when uo /umf is 1000, a condition approached with very fine particle bubbling beds or circulating beds. In the latter instance, the use of the Ergun relationship is questionable except for the dense lower part of the bed.
5.5Range of Validity of Simplified Scaling
To determine the validity of the simplified scaling laws over a wide range of conditions the simplified scaling laws have been used (Eq. 53) to design hypothetical models whose linear dimensions are 1/4 and 1/ 16, respectively, of the linear dimensions of a model designed using the full set of scaling laws, Eq. (37). To determine the validity of the smaller, simplified models, the dimensionless drag coefficient β L/ρsuo will be compared between the simplified models and the model using the full set of scaling laws. Figure 23 shows a comparison of the exact model and the simplified models for a pressurized fluidized combustor. Using the full set of scaling laws, the exact model, fluidized by ambient air, is approximately the same size as the combustor. The simplified models are reduced in size by their respective assumed length scale. The other parameters of the simplified model are then calculated to match the simplified parameters. For example, when the length scale is reduced to 1/4 that of the exact model, the velocity is reduced by 1/2 to keep the Froude number constant. The particle diameter is then reduced appropriately to keep the ratio of uo/umf constant. These calculations were carried out over a range of particle Reynolds numbers, RepE, based on the full scaling law, or exact, model. In the simplified scaling relationships, the Reynolds number is not maintained constant. The concern is how much the drag coefficient is impacted by the shift in Reynolds number. It was found that the particle Reynolds number for the 1/4 scale simplified model remained roughly equal to 0.34 RepE over a wide range of values for RepE, whereas the particle Reynolds number for the 1/16 scale model was roughly 0.12 RepE.