Thesis 15
Thesis 15
Thesis 15
A Thesis
Submitted to the College of Engineering
of Nahrain University in Partial Fulfillment
of the Requirements for the Degree of
Master of Science
in
Chemical Engineering
by
Zahra,a Hamid Mohammed
(B. Sc. In Chemical Engineering 2004)
December 2007
CERTIFICATION
Signature:
Name:
Dr. Kamal Shakir
Abdulmasih
( Supervisor )
Date:
9/ 1 /2008
Signature:
Name:
Date: 9 / 1 / 2008
Summary
Contents Page
Abstract I
List of Contents 5
Notations V
List of Tables X
List of Figures XI
References 65
Appendixes
Appendix A Minimum Fluidization.
IV
Notations
Symbols Notations Units
as = Proportionally factor m2/ s
bs = Packing factor
Cb = Bulk concentration. mole / m3
Cf = Specific heat of gas J / kg. K
Cm , C m = Relative and mean relative mass capacity
respectively
Cmf ,Cms = Specific mass capacity of gas and particles kg / kg
respectively
Cs = Concentration at the surface. mole /m3
e = Void fraction. -
V
Symbols Notations Units
fο = Hydrodynamic parameter -
Fr = Froude number. u 2 / dp g -
VI
M = Molecular weight. g/mole
Rmp, Rmw = Mass transfer packet and contact resistance m2. s /kg
respectively
Rep = Reynolds number based on the diameter of the -
inert particles. ρud p / µ
Sc = Schmidt number. µ / ρ Dv -
Sh = Sherwood number. kg ds / Dv -
t = Temperature. ْC
Ts = Surface temperature. K
Greek Letters
∆Pb = Bed pressure drop. mm Hg
ψ = Sphericity. π d p2 / aext -
VIII
ρf , ρs = Gas and particle density respectively kg/m3
δ = Carbon transfer coefficient m/s
α = Parameter of the simplified packet model -
λmp = Effective mass conductivity for a packet kg / m.s
θ = Mass transfer potential (X/Xw) Potential
unit
θο = Potential in bulk medium of a bed Potential
unit
θw = Potential at an immersed contact time and its mean
value respectively Potential
unit
τ = Time s
τ b ,τ b = Bubble residence contact time and its mean value s
respectively
τ p ,τ p = Packet residence contact time and its mean value s
respectively
εp = Packet porosity -
Subscripts
b Bubble
e Equilibrium
f Gas
m Mass (minimum)
0 Bulk medium of a fluidized bed
s Solid(particle)
w Immersed surface
IX
List of Tables
X
List of Figures
XII
Chapter One
Introduction
2
Chapter Two
Literature Survey
2.1 Fluidization
5
which is a sufficient fraction of the bed pressure drop. Many operating
problems can be traced back to poor distributor design. Some distributor
designs in common use are shown in Fig. (2.2)[14].
(a) drilled plate (b) cap design (c) continuous horizontal slots
3
e mf ρs − ρ 2
U mf = 0.0055 d g --- (2.1)
(1 − e ) µ
7
Davidson et al., 1966[17], has shown that the Umf is a function of the
square of particle diameter; as a result, the quantity of air required for the Umf
changes as the products particle size change.
Wen and Yu, 1966[18], produced an empirical correlation for Umf for
gas fluidization the Wen and Yu correlation is often taken as being most
suitable for particles larger than 100 µm, where as the correlation of Baeyens
in 1974 [19], shown below in equation (2.2), is best for particles less than
100µm.
p (ρ p − ρ )
d 1.8 0.934 0.934
g
U mf = --- (2.2)
110µ 0.87 ρ 0.066
ρu mf d p
= (C 2 + 0.0408Ar ) − C
1/2
--- (2.3)
µ
using the value 33.4 and 33.7. It would seem reasonable to use a mean value
of 30.
8
Gupta and Sathiyamoorthy, 1999[22] and Kunii and Levenspiel,
1991[13], studied the pressure drop vs. velocity relationships during
fluidization and defluidization. They plotted curves for different sand sizes for
increasing flow rates (fluidization) and decreasing flow rates (defluidization).
In a typical curve showing the variation of bed pressure drop with superficial
air velocity for the sand size of 375 microns is shown in Fig. (2.4), they
showed that the pressure drop was higher on fluidization than on
defluidization. It was evident because initially the bed was fewer perms.
( Re p )
C3
Sh = 2.0 + C1 ( Sc )
C2
--- (2.4)
12
2.4 Mass Transfer in Fixed Bed.
If the gas is pumped upwards through a granular bed at low flow rate
the fluid percolates through the pores with no perturbation of the bed height,
it's called a fixed bed.
For fluid-solid reactions, the pressure drop for flowing through the
fixed bed is higher than for a fluidized bed with a same flow rate. A schematic
of the pressure drop versus flow rate is shown is shown in Fig. (2.5). At low
fluid velocities the pressure drop is approximately liner with flow rate, this is
expected behavior for packed beds. After achieving incipient fluidization
increasing the fluid flow velocity doesn't result in any significant increase in
the pressure drop as the bed expands to reduce resistance to flow [1].
For
d pu ρ
Re p = > 80
µ
Wakao and Funazkri, 1978[31], correlated the published mass transfer
data, of particle to fluid mass transfer coefficient in packed bed, for the axial
fluid effective dispersion coefficient.
The corrected Sherwood numbers in the range of Reynolds number
from about 3 to 10000 are correlated by:
14
that in the proposed correlations not only the influence of gas velocity, u, on
the fluid to particle mass transfer coefficient, kg, is hidden but also that of the
void fraction, e, on kg [33].
Several correlations for the calculation of the Nusselt and Sherwood
dimensionless numbers are reported. Since these correlations are mainly
based on experimental investigations performed under laboratory conditions,
they may be different to the situation in large scale reactors.
Hurt, 1934[34], measured the height of transfer unit for gas film
controlled system in fluidized beds with different sizes and shapes of packing.
No spherical particles were used, but the cylindrical particle data showed
close agreement between heat and mass transfer factors as Chilton and
Colburn 1934 had suggested. No such agreement between heat and mass
transfer data was obtained for other packing shapes. Hurt did not report
fraction voids or surface area of the beds he used.
McCune and Wilhelm, 1949[35], attempted to relate the mass transfer
group ejd with the modified Reynolds [dpG/eu], but found that the group
obtained from their data for fluidized beds did not correlate with the
corresponding groups for packed beds.
Resnick and White, 1949 [36], calculated the mass transfer coefficient
of naphthalene crystals of five different size ranges (250-1000microns) in air,
hydrogen, and carbon dioxide at a temperature of 298K and rates between (
0.01 and 1.5 kg/m2.s ). They expressed (J-factor) and plotted against Reynolds
number as shown in Fig. (2.6).
Gamson, 1951[37], utilizing the available mass transfer data for packed
and fluidized beds related the mass transfer modulus [ j d /(1 − e )0.2 ] to the
15
Figure (2.6): (jd) Factor for The Transfer of Naphthalene
Vapor to Air in Fixed and Fluidized Beds
0.863
ej d = 0.01 + 0.58
--- (2.7)
Re p − 0.483
Re p > 1 For
16
In fluidized bed, the solid particles are sufficiently separated so that in effect
there is mass transfer between a gas and single particles. The most widely
used correlation for this purpose is the equation of Froessling for mass
transfer to single sphere:
17
between particles and fluid. They concluded that the mechanism of heat is not
analogous to that for mass transfer, because heat may be transported by both
gas and the fluidizing particles but mass can be transferred only by the gas
phase, when the particles are non-adsorptive.
Ziegler and Holmes, 1966 [40], studied mass transfer from fixed
surface to gas fluidized beds, mass transfer coefficients were measured for the
diffusion of water vapor from a saturated porous sphere into various air-
fluidized beds of solid particles. Naphthalene diffusion from coated flat plate
into fluidized beds was also studied. They found that for the case in which
particle adsorption is negligible, the Sherwood number could be correlated
by:
Sh = She + f ( y ) --- (2.8)
--- (2.9) f ( y ) = C1 y C 2
From experimental data, the value of C1 and C2 found are 33.7 and 0.4
respectively with an average standard deviation of 16.5%.
For Naphthalene diffusion, unfortunately adsorption on the particles
increased the transfer rate, for which reason data are inconsistent [33].
Ciborowski and Kopec, 1985[41], reported previous works, for mass
transfer from an immersed object to a gas fluidized bed, which can be
18
summarized as follows:
1- The ratio of surface-to-inert bed mass transfer coefficient to that between
the surface and a gas; (kg / kge), varies from 1.1 to 10.
2- ( kg ) is a slightly increasing function of the superficial gas velocity and for
large velocities is independent of u.
3- ( kg )rises with an increase in fluidizing particle size, while h diminishes as
dp grows.
−0.4
⎛ Re ⎞
jd ε = 0.64 ⎜⎜ p ⎟⎟
⎝ 1− ε ⎠ --- (2.11)
With a mean absolute deviation of (16.88 %) and a mean relative deviation of
(- 3.71 %).
1
jd
0.1
Rep
Weimin GAO et al., 2004 [50], studied the mass transfer in fluidized
bed. A steel work piece covered with carbon was used in their investigation.
The carbon transfer coefficient was determined from the carbon distribution
21
within the diffusion layer of the sample. An empirical relationship of the
carbon potential as a function of carburizing atmosphere, bed temperature (-
90ْ,-55ْ,-30ْ, 0ْ, 30ْ, 60ْ) and fluidizing velocity (1.32Umf, 2.12Umf) was
( )
δο = 2.54 ×10−5 + 3.3 ×10−7C p6.22 exp ⎜ −
⎝ RT
⎛ 40.365 ⎞
⎟
⎠
--- (2.13)
⎛ ⎞
⎜ 8.805 ⎟
γ = exp ⎜ − 2⎟
--- (2.14)
⎝ (
⎜ ρ p U −U
mf
⎟
⎠ )
22
Chapter Three
3.1 Introduction
k y ′ = (1 − f ο )k yp ′ + f ο k yb ′ --- (3.1)
Where
τb
fο = --- (3.2)
τb +τ p
In contrast to heat transfer theory [54], where the heat within a packet
is transferred through gas and particles and the accumulation of heat within
particles plays a dominant role, these workers assumed that mass within a
packet is transferred only by gas between particles occurs. Thus the mass
transfer coefficient to a packet was found [51, 52] to be
23
1/2
⎛D ⎞
k yp ′ = 2 ρf ⎜ p ⎟ --- (3.3)
⎜ πτ p ⎟
⎝ ⎠
1/2
⎛D ⎞
k yb ′ = 2 ρf ⎜⎜ b ⎟ --- (3.4)
⎝ πτ b
⎟
⎠
1/2
⎛ ε D M ρ (1 − ε p ) ⎞
k yp = ρf ⎜ p f s s ⎟ --- (3.5)
⎜
⎝ τ p
⎟
⎠
k yp L
Where Sh p = transformed into the dimensionless form:
D f ρf
1/2
⎛ ε p (1 − ε )M s ρs L 2 ⎞
Sh p = ⎜ ⎟ --- (3.6)
⎜
⎝
D f τ p ⎟
⎠
For particles with large mass capacities, the mass transfer coefficients
become higher and greater similarity between surface-to-bed mass and heat
transfer mechanism.
In this chapter, the mass capacity process was investigated and
24
described on the basis of the modified packet model including the mass
contact resistance. For the contact resistance control region the alternative
simplified packet model is developed.
25
The different equation which describes non-steady-state mass transfer into the
packet is:
∂Y ∂X ∂ ⎛ ∂θ ⎞
ε p ρf + (1 − ε p ) ρs = ⎜ λmp --- (3.9)
∂τ ∂τ ∂z ⎝ ∂z ⎟⎠
θ (τ = 0, z ) = θο , θ (τ , z → ∞ ) = θο , θ (τ , z = 0 ) = θw ---(3.11)
For adsorbing particles the first term in curly brackets in equation (3.10) can
be neglected, since ( ρsC ms >> ρf C mf ).
Where
(1/ ε p − 1)(1 − 1/ B ) 1
β= + --- (3.15)
1
+ 0.28ε p0.63B
0.18 εp
B
D s ρs C m
and B =
D f ρf
The solution of equation (3.10) for constant λmp , Cms and Cm and for boundary
⎡ ⎛ (1 − ε ) ρ C 1/ 2 ⎤
θ − θw ⎢z ⎜
⎞
⎥
= erf p s ms
⎟ --- (3.16)
θο − θw ⎢2⎜ λmpτ ⎟ ⎥
⎣⎢ ⎝ ⎠ ⎦⎥
1/ 2
(
⎛ 1− ε p
=⎜
) ρsC ms λmp ⎞⎟ (θ − θ ) --- (3.17)
ο
⎜ πτ ⎟
w
⎝ ⎠
In order to get the time-averaged mass flux Np information about the packet
residence time distribution is necessary. In this work an exponential
distribution is assumed, so that
1 ∞
Np = N (τ )exp(−τ / τ p )d τ --- (3.18)
τ p 0∫ pτ
After substitution of equation (3.17) into equation (3.18) and integration, the
packet mass transfer resistance is calculated as:
1/ 2
Np ⎛ τp ⎞
R mp = = 2⎜ ⎟ --- (3.19)
θw − θο ⎜ (1 − ε p ) ρs λmpC ms ⎟
⎝ ⎠
Cms = Cm Cmf = 1
Cms = Xw Cmf = Xw / Cm
λmpC ms = βε p D f ρf C ms
2
/C m for θ = X / X w
28
In assumption that there is an additional contact mass resistance Rmw in
the packet zone adjacent to the surface which is proportional to the particle
diameter ds. The formula for mass transfer coefficient will then be:
From equations (3.19) and (3.22) it follows that the mass transfer coefficient
for packets ( k θ p ) is a function of the driving force (θw − θο ) if
θ =Y orθ = X / X w . Hence the driving force should be selected
experimentally to assure the smallest variation of k θ p .
From equations (3.19) and (3.22) found that for a negligible contact
resistance, i.e. for sufficiently small values of ( d s2C ms / τ p ), the mass transfer
⎛ ∂X r ⎞
D s ρs ⎜ ⎟
⎝ ∂r ⎠r =d s / 2
(
= k ys Y w −Y r =ds / 2 )
--- (3.23)
The well known solution of equation (3.22) leads to the following expression
describing the time dependence of the mean concentration X.
Xw −X ∞
2
= 1.5Bi ms ∑ ηi exp(− µi2 Fo ms ) --- (3.25)
Xw −X ο i =1
Where
sin 2 µi
ηi = and µi is the root of (1-Bims)
µi3 ( µi − 0.5sin 2µi′ )
30
After deriving the result of different of equation (3.25) by equation (3.25):
∞
ηi µi2 exp(− µi2 Fo ms )
∂X 4D s i∑
= 2 =1∞ (X w − X ) --- (3.26)
∂τ ds ∑ η exp(− µ Fo ms )
2
i i
i =1
When (Bims < 10, and Foms < 10-3) this expression can be approximated with
an accuracy of + 10 % by the form resulting from Foms = 0.
∂X hD s --- (3.27)
= 2 (X w − X )
∂τ ds
∞ ∞
Where h = 4 ∑ ηi µi2 / ∑ ηi is a constant while a packet is in contact with the
i =1 i =1
surface. Note that this approximation should be confused with that well
known for large Fourier numbers.
Equation (3.27) may be extended to the more general case of non-
spherical particles with an unknown or undetermined diffusivity as follows:
∂X as ∂as
= (X − X ) with =0 ---(3.28)
∂τ d s2 w ∂τ
or ,after integration
∂X as aτ
= 2 exp(− s 2 )(X w − X ο ) --- (3.29)
∂τ d s ds
The instantaneous mass flux from the surface to a packet is expressed as:
bs ∂X
N pτ = m
2 s ∂τ
--- (3.30)
ds
Where ms is the particle mass and ( bs / ds2 ) is the number of particles per unit
area of the surface.
For non-spherical particles only an approximate value of the packing factor (
bs ≈ 1 ) can be found. While for spherical particles ( bs = 1 → 4 / 3 ).
The time-averaged mass flux results from equations (3.18), (3.29) and
31
(3.30) as follows:
π as bs ρs / d s
Np = (X w − X ο ) --- (3.31)
6(1 + asτ p / d s2 )
It follows that both coefficients (as and bs) cannot be determine directly.
Hence, the parameter of the simplified packet model ( α ) cannot be calculated
from equation (3.33) and must be found on the basis of experimental data.
The total mass flux from an immersed surface to a mass capacity
fluidized bed should be calculated from the expression analogous to heat
transfer:
N = (1 − f ο )k θ p (θw − θο ) + N ′ --- (3.34)
Where N ′ = k y′ (Y w −Y ο )
transfer coefficient for an inert fluidized bed and can be predicted from
equation (3.1).
32
Chapter Four
Experimental Work
Table (4.1)
Properties of Sand Particles
33
4.1 General Description of Flow System
A Q.V.F glass column of (8 cm) internal diameter, and (70 cm) in high
was used as fluidization column. A porous material was placed at the bottom
of the column to support the packing material.
The pressure drop across the bed was measured by the use of (U-tube)
manometer made of glass with total length of (0.75 m). The manometer was
placed on a wide sheet of wood with a measuring tape for the measurement of
the level difference of the liquid (water) inside the tube.
The heater has been used was an electrical heater placed inside 2"
Q.V.F. glass tube. The variation in the heat supplied from the heater was
controlled by the use of a variac connected directly to the heater.
Table 4.2
Conditions of Mass Transfer
36
37
4
1
5
6
3
Figure (4.2)
1. Rotameter
2. Variac
3. Heating Equipment
4. Temperature Reader
5. Fluidization Column
6. Manometer
38
4.2 Experimental Procedure
4.2.1 Determination of Minimum Fluidizing Velocity
The pressure drop of the bed ( ∆pb ) can be determined with the range of
gas velocities by subtracting the pressure drop of distributor ( ∆pd ) that find
out for a range of superficial gas velocities, from total pressure drop ( ∆pt )
that find out for a range of superficial gas velocities after loading known
weight of sand particles having known diameter into the bed to a static level
of 30 cm.
A curve of pressure drop across the bed ( ∆pb ) versus superficial gas
40
7-When conditions reach to steady state (constant flow rate and constant
temperature); the coated sphere was lowered inside the column 15 cm
above the distributor surface.
8-Every 5 minutes, the sphere was taken out of the bed and the change of
weight was measured by digital balance.
9-The previous steps 1-8, were repeated for arrange of air superficial
velocities and a range bed temperatures.
N = k g (C s − C b ) --- (4.2)
3729.3 ---
log Ps = − + 11.450
Ts (4.4)
Where Ps in mHg and Ts in K.
5- Values of Sherwood number for mass transfer from the sphere to a bed of
fluidized particles were calculated by the equation:
c c3
She = 2 + C1 Re p2 Sc --- (4.6)
43
Chapter Five
The above equation shows good agreement with the previous work by
Ranz and Marshal [28]. Fig.(5.2) illustrates the experimental results and
correlation.
For mass transfer in fluidized bed, a set of experiments were performed
to determine the value of mass transfer coefficient from the sphere to the
fluidized bed. The variables in this investigation which effect mass transfer
were: air velocity, sphere surface temperature, size of fluidizing particles and
sphere size. The range of sphere surface temperature varied from ambient
temperature to a temperature below the melting point of naphthalene. The air
44
velocity is chosen to be within the range (1-1.4) Umf, because this range of
flow is usually used in industrial practice. The particle size of sand was
selected to be as fine particles in order to get a smooth fluidization. The
experimental conditions and results are listed in Tables (A-3 and A-4)
respectively.
• Samples of data sheet for experiments are listed in Tables (5.5, 5.6 and
5.7) respectively.
• Samples of the experimental result showing effect of air temperature
on Sherwood number are illustrated in Figs. (5.3 and 5.4).
• Samples of the experimental result showing effect of air flow rate on
Sherwood number are illustrated in Figs. (5.5).
• The effects of both air temperature and particle size on Sherwood
number are showed in Figs. (5.8, 5.9 and 5.10).
• The effects of both air flow rate and particle size on Sherwood number
are showed in Figs. (5.11, 5.12, 5.13).
• The effects of both particle size and different temperature on
Sherwood number are showed in Figs. (5.14, 5.15, 5.16, and 5.17).
45
Sh = She + f ( y ) --- (5.2)
is chosen as a power function of Y ,that is
C2
f ( y ) = C1 y --- (5.3)
Two attempts have been made to correlate the experimental results:
1- The first attempt was made by choosing the dimensionless function ,Y,
to be equal to [(G − G mf ) µ /[(ψ d p )2 ( ρ p − ρ ) ρ g ]] as given by Ziegler and
Holmes[40],i.e.:
(G − G mf ) µ C2
Sh = She + C1[ 2
] --- (5.4)
(ψ d p ) ( ρ p − ρ ) ρ g
( ρ p − ρ ) ρ g (ψ d p ) C (ψ d p )(G − G mf ) C3
Sh = She + C1[ ] 2[ ] --- (5.5)
(G − G mf )2 µ
46
agreement with experiments, in which 97% of the points have an error less
than 25%, consequently this correlation obtained from the present work.
5.2 Comparison of Experimental Results With previous Works
and Model
47
Table (5.1) Comparison of the Orders of Magnitude of
The Experimental Parameters
Silica Sand-air-water Sand-air-
gel-air- Rmp/Rmw~ ∞ naphthalene
water Rmp/Rmw~ ∞
Rmp/Rmw
48
]Table (5.2) Selected Data Sheet for Experiment
of Mass Transfer in Empty Bed
Note:
T1 = Temperature below the sphere.
T2 = Temperature above the sphere.
49
Table (5.3) Selected Data Sheet for Experiment
of Mass Transfer in Empty Bed
Note:
T1 = Temperature below the sphere.
T2 = Temperature above the sphere.
50
Table (5.4) Selected Data Sheet for Experiment
of Mass Transfer in Empty Bed
Note:
T1 = Temperature below the sphere.
T2 = Temperature above the sphere.
51
Table (5.5) Selected Data Sheet for Experiment
of Mass Transfer in Fluidized Bed
Note:
T1 = Temperature below the sphere.
T2 = Temperature above the sphere.
52
Table (5.6) Selected Data Sheet for Experiment
of Mass Transfer in Fluidized Bed
Note:
T1 = Temperature below the sphere.
T2 = Temperature above the sphere.
53
Table (5.7) Selected Data Sheet for Experiment
of Mass Transfer in Fluidized Bed
Note:
T1 = Temperature below the sphere.
T2 = Temperature above the sphere.
54
25.00
20.00
15.00
Sh (exp. )
10.00
5.00
0.00
200.00
100.00
300.00
Sh (exp.)
200.00
100.00
300.00
Sh (exp. )
100.00
Sh (calc.) 300.00
200.00
100.00
300.00
Sh (calc.)
200.00
100.00
Sh (calc.)
200.00
150.00
100.00
1.4 Umf
1.3 Umf
1.2 Umf
300.00
Sh (calc.)
200.00
100.00
1.2 Umf
360.00
Sh (calc.)
320.00
280.00
39 C
250.00
51 C
Sh (calc.)
55 C
200.00
66 C
150.00
100.00
51 C
300.00
55 C
250.00 66 C
Sh (calc.)
200.00
150.00
100.00
300.00
Sh (calc.)
200.00
100.00
1.4 Umf
300.00
1.3 Umf
1.2 Umf
Sh (calc.)
200.00
100.00
Temperature = 51 C
300.00
Sh (calc.)
1.4 Umf
1.3 Umf
1.2 Umf
200.00
100.00
Temperature = 55 C
300.00
Sh (calc.)
1.4 Umf
1.3 Umf
1.2 Umf
200.00
100.00
Temperature = 66 C
300.00
Sh (calc.)
1.4 Umf
1.3 Umf
200.00 1.2 Umf
100.00
Conclusion and
Recommendations for Future work
6.1 Conclusions
63
6.2 Recommendations for Future work
From the present study it was noticed that further studies in the
following areas would be desired:
1- Studying the effect of geometry on transfer phenomenon by using
different shapes of immersed objects which have important applications
in industry.
2- Using different systems to study the common uses for mass transfer in
fluidized bed.
3- Study of the mass transfer of the solid materials applied in Iraqi
industries, by fluidized mass transfer to improve products quality.
64
References
http://www.informaworld.com/smpp/content~content=a777773479~db=all
47- Guedes de Carvalho, J. R. F. and J. M. Delegado, "Mass Transfer from a
Large Sphere Buried in Packed Bed along Which Liquid Flows", (1999).
48- Philipp Schlichthaerle, Joachim Werther," Influnce of the Particle Size
and Superficial Gas Velocity on the Sublimation of Pure Substances in
Fluidized Beds of Different Sizes ", (2000).
http://www.informaworld.com/smpp/content~content=a7777911662~db=all
49- Delgado JMPQ, Guedes de Carvalho JRF," Measurement of the
Coefficient of Transverse Dispersion in Packed Beds Over a Range of
68
Values of Schmidt Number (50-1000)",(2001).
50- Weimin Gao, John M. Long, Lingxue Kong and Peter D. Hodgson,
"Influence of the Geometry of an Immersed Steel Work piece on Mass
Transfer Coefficient in Fluidized Bed", School of Eng. and Tech.,
Deakin University, Australia,(2004).
51- A.P.Baskakov and W.M. Suprun, Khim. Promst., (1970).
52- M. N. Markova , Autoreferat Dissert.,(1972).
53- E. N. Prozorov , Izuest. VUZ-ov Khim. Technol.,(1976).
54- L. A. Kirk and F. L. Hudson, Trans. Inst. Chem. Eng.,(1966).
55-Teiji Yokota, Y. Hidaka and T. Yasutomi, Kagaku Kogaku
Ronbunshu,(1975).
56- J. Kopec, Ph. D. Thesis, Warsaw Technical University,(1981).
57- Perry, R.H. and Chiton, C.H., Chemical Engineers Handbook,(1973).
58- Howard W. Newton "Fluidized Bed Reactor ".
http://www.engin.umich.edu./~cre/12chap/html/fluidizedbed.pdf.
59- H.Scott Fogler, "Element of Chemical Reaction Engineering",(2006).
http://www.eng.buffalo.edu/Courses/ce561/Day31.pdf
60- J.M.P.Q. Delgado "A Critical Review of Dispersion in Packed Beds ",
Heat and Mass Transfer,(2006).
http://paginas.fe.up.pt/ceft/pdfs/2006_HTM1_JMPQD.pdf
61- Tatemoto, Y., Mawatari, Y., Yasukawa, T., and Noda, K., "Numerical
Simulation of Particle Motion in Vibrated Fluidized Bed ", (2004).
62- Krell, L., Kunne, H. J. and Morl, L., "Flow Regimes and Heat and Mass
Transport in Gas-Fluidized Beds of Solids ", Inter. Chem. Eng. Vol.30,
(1990).
69
Appendix - A
10.00
1.00
Pbed (CmH2O)
0.10
0.01
10.00
^ P bed (CmH2O)
1.00
Air velocity (Cm/s)
1.00
1.00
Air Velocity (Cm/s)
10.00
^P Bed (CmH2O)
1.00
Air Velocity (Cm/s)
A-2
***Calibration of Air Flow meter:
The variable area rotameter has been calibrated for air at 101.3 KN /m2
and 39 0C ( ρ = 0.6013 kg / m3 and ν = 1.326 * 10-5 m3. s-1).
The flow rate at other conditions may be calculated from:
1.204
Actual flow rate = indicated flow rate ×
ρm
Where ρm is the density of the air in the meter in kg/m3 .
P
Note: from the ideal gas equation, ρ =
RT
0.6013
ρm
A-3
Dimensional Analysis:
Y = [ψ d p ]a [( ρ p − ρ ) ρ ]b [G − G mf ]c g d µe ---- (A-2)
The common groups for mass transfer are the Sherwood number, the
Schmidt number, and the Reynolds number. The Sherwood number contains
the mass transfer coefficient and the diffusivity. The Schmidt number
contains only the physical properties of the fluid and its active component.
The Reynolds number is, of course, a measure of flow rate.
The number of dimensionless group obtained from the dimensional
analysis is equal to the number of variables, n=5, minus the number of
fundamental dimensions, r=3, and hence two dimensionless groups will be
obtained.
In term of fundamental dimensions:
m2 b m c L d m e
1 = [L ]a [ ] [ ][ ] [ ]
L6 L2θ θ 2 Lθ --- (A-3)
For L 0 = a – 6b – 2c + d – e ---(2)
These three equations can be solved to give all the exponents in terms of
two:
From (1) c = - 2b - e
From (3) d = - 0.5c - 0.5e
= -0.5 (-2b-e) - 0.5 e
=b
From (2) 0 = a - 6b - 2 (-2b-e) + b- e
=a+e-b
a =- e + b
( ρ p − ρ ) ρ g (ψ d p ) (ψ d p )(G − G mf )
Y =[ ]b [ ]−e --- (A-5)
(G − G mf ) 2 µ
From equation above, one can notice that the first term is the invert of Froude
number (Fr) and the second is the modified Reynolds number (Re).
A-5
Table (A-1) Experimental Conditions for mass Transfer in
Empty Bed
A-6
Table (A-2) Experimental Results for Mass Transfer in
Empty Bed
A-7
Table (A-3) Experimental Conditions for mass Transfer
in Fluidized Bed
A-11
ﺍﳋﻼﺻﺔ
ﻋﻨﻮﺍﻥ ﺍﻟﺒﺤﺚ :ﺍﻧﺘﻘﺎﻝ ﺍﳌﺎﺩﺓ ﻣﻦ ﺳﻄﺢ ﺍﻟﻨﻔﺜﺎﻟﲔ ﺍﻟﻐﺎﻃﺲ ﺍﱃ ﻃﺒﻘﺔ )ﺭﻣﻞ_ﻫﻮﺍﺀ( ﳑﻴﻌﺔ.
ﻭﺗﺴﺘﻌﻤﻞ ﺣﺎﻟﻴﺎ ﰲ ﻋﻤﻠﻴﺎﺕ :ﺍﻟﻔﺼﻞ ,ﺍﻟﺘﺼﻨﻴﻒ ,ﺍﻟﺘﺠﻔﻴﻒ ﻭﺧﻠﻂ ﺍﻟﺪﻗﺎﺋﻖ ,ﺍﻟﺘﻔﺎﻋﻼﺕ ﺍﻟﻜﻴﻤﻴﺎﻭﻳﺔ ﻭﻋﻤﻠﻴـﺎﺕ ﺍﻋـﺎﺩﺓ
ﺍﻟﺘﺼﻨﻴﻊ .ﻭﺍﺣﺪﺓ ﻣﻦ ﻫﺬﻩ ﺍﻟﻌﻤﻠﻴﺎﺕ ﻫﻲ ﺍﻧﺘﻘﺎﻝ ﺍﳌﺎﺩﺓ ﻣﻦ ﺳﻄﺢ ﻏﺎﻃﺲ ﺍﱃ ﻃﺒﻘﺔ ﻏﺎﺯﻳﺔ ﳑﻴﻌﺔ.
ﻣﻌﺎﻣﻞ ﺍﻧﺘﻘﺎﻝ ﺍﳌﺎﺩﺓ ﻫﻮ ﺛﺎﺑﺖ ﺍﻻﻧﺘﺸﺎﺭ ﻭﺍﻟﺬﻱ ﻳﺘﻌﻠﻖ ﺑﻨﺴﺒﺔ ﺍﻧﺘﻘﺎﻝ ﺍﳌﺎﺩﺓ ,ﻣﺴﺎﺣﺔ ﺍﻧﺘﻘـﺎﻝ ﺍﳌـﺎﺩﺓ ,ﻭﳏﺼـﻠﺔ
ﺍﻟﺘﺮﻛﻴﺰ ﻭﺍﻟﱵ ﲤﺜﻞ ﻓﺎﻋﻠﻴﺔ ﺍﻻﻧﺘﻘﺎﻝ .ﻭﳑﻜﻦ ﺍﻥ ﲢﺴﺐ ﻣﻦ ﻋﺪﺓ ﻣﻌﺎﺩﻻﺕ ﻭﻋﻼﻗﺎﺕ ﻭﻧﻈﺮﻳﺎﺕ ﳐﺘﻠﻔﺔ ﺩﺍﻟﺔ ﻋﻠﻰ ﺧﺼﺎﺋﺺ
ﺍﳌﺎﺩﺓ.
ﰎ ﺍﺳﺘﺨﺪﺍﻡ ﻋﻤﻮﺩ ﺯﺟﺎﺟﻲ ﺑﻘﻄﺮ ﺩﺍﺧﻠﻲ ) 8ﺳﻢ( ﻭﺑﻄﻮﻝ ) 70ﺳﻢ(,ﻭﻛﺎﻥ ﺍﻟﻐﺮﺽ ﻣﻦ ﻫﺬﺍ ﺍﻟﺒﺤﺚ ﻫﻮ ﺍﳚﺎﺩ
ﻋﻼﻗﺔ ﺭﻳﺎﺿﻴﺔ ﺗﺮﺑﻂ ﻣﻌﺎﻣﻞ ﺍﻧﺘﻘﺎﻝ ﺍﳌﺎﺩﺓ ﺍﱃ ﻃﺒﻘﺔ ﻏﺎﺯﻳﺔ ﳑﻴﻌﺔ ﺑﻌﺪﺓ ﻣـﺘﻐﲑﺍﺕ ﰎ ﺩﺭﺍﺳـﺘﻬﺎ ﻭﻫـﻲ :ﺩﺭﺟـﺔ ﺣـﺮﺍﺭﺓ
( ρ p − ρ ) ρ g (ψ d p ) C2 (ψ d p ) )(G − G m f C3
[S h = S he + C 1 ] [ ]
(G − G m f ) 2 µ
ﺍﻟﺜﻮﺍﺑﺖ ﻫﻲ C3, C2, C1:ﺗﺴﺎﻭﻱ ﻋﻠﻰ ﺍﻟﺘﻮﺍﱄ 0.1284, 0.07497 , 16.8574ﻭﺑﻨﺴـﺒﺔ ﺧﻄـﺄ ﺗﺴـﺎﻭﻱ
.1.54%
ﰲ ﻫﺬﺍ ﺍﻟﻌﻤﻞ ﻭﺟﺪ ﺍﻥ *: Sherwoodﻟﻪ ﻋﻼﻗﺔ ﻋﻜﺴﻴﺔ ﻣﻊ ﺣﺠﻢ ﺍﻟﺪﻗﺎﺋﻖ
* ﻟﻪ ﻋﻼﻗﺔ ﻃﺮﺩﻳﺔ ﻣﻊ ﺳﺮﻋﺔ ﺟﺮﻳﺎﻥ ﺍﳌﺎﺋﻊ
* ﻟﻪ ﻋﻼﻗﺔ ﻋﻜﺴﻴﺔ ﻣﻊ ﺩﺭﺟﺔ ﺣﺮﺍﺭﺓ ﺍﻟﺴﻄﺢ
ﰎ ﺍﺳﺘﻌﺮﺍﺽ ﻣﻮﺩﻳﻞ ﺭﻳﺎﺿﻲ ﻻﻧﺘﻘﺎﻝ ﺍﳌﺎﺩﺓ ﻣﻦ ﺳﻄﺢ ﻏﺎﻃﺲ ﺍﱃ ﻃﺒﻘﺔ ﳑﻴﻌﺔ ﳍﺎ ﺍﻟﻘﺎﺑﻠﻴﺔ ﻋﻠﻰ ﺍﻣﺘﺼﺎﺹ ﺍﳌﺎﺩﺓ .ﻭﰎ ﻣﻘﺎﺭﻧﺔ ﻧﺘﺎﺋﺞ
ﺍﳌﻮﺩﻳﻞ ﻣﻊ ﺍﻟﻨﺘﺎﺋﺞ ﺍﻟﻌﻤﻠﻴﺔ ﳍﺬﺍ ﺍﻟﺒﺤﺚ ﻭ ﺍﻟﺒﺤﻮﺙ ﺍﻟﺴﺎﺑﻘﺔ.
ﺍﻨﺘﻘﺎل ﺍﻟﻤﺎﺩﺓ ﻤﻥ ﺴﻁﺢ ﺍﻟﻨﻔﺜﺎﻟﻴﻥ ﺍﻟﻐﺎﻁﺱ ﺍﻟﻰ ﻁﺒﻘﺔ )ﺭﻤل-ﻫﻭﺍﺀ(ﻤﻤﻴﻌﺔ
ﺭﺴﺎﻟﺔ
ﻤﻘﺩﻤﺔ ﺇﻟﻰ ﻜﻠﻴﺔ ﺍﻟﻬﻨﺩﺴﺔ ﻓﻲ ﺠﺎﻤﻌﺔ ﻨﻬﺭﻴﻥ
ﻭﻫﻲ ﺠﺯﺀ ﻤﻥ ﻤﺘﻁﻠﺒﺎﺕ ﻨﻴل ﺩﺭﺠﺔ ﻤﺎﺠﺴﺘﻴﺭ ﻋﻠﻭﻡ
ﻓﻲ ﺍﻟﻬﻨﺩﺴﺔ ﺍﻟﻜﻴﻤﻴﺎﻭﻴﺔ
ﻤﻥ ﻗﺒل