Final Project 2013
Final Project 2013
Final Project 2013
Study on Determination of
Volumetric Mass Transfer
Coefficient in Stirred Vessels
Done By
Laila Moosa Al-Blushi, Maryam Said Al-Oasimy Sahar Payam Allah
Talebi, Shamsa Ahmed Al-Haddabi, and Sumaya Hamood Al- Maawali
Declaration
We hereby declare that this report is based on our original work. We also declare that it
has not been previously and concurrently submitted for any other degree or award.
Signatures
Signatures
Acknowledgment
The project group would like to thank their supervisor Dr. Ahmed AlDallal for his assistance, valuable advice, and support along the project.
They also appreciate Mr. Ibrahim Al-Ajmai and Mr. Rashid for their
never end support, and their continuous help during the design,
manufacturing and experimenting the aeration device build. Furthermore,
the group members would also like to express their gratitude to their loving parent and
friends who had helped and encouraged them. Without the mentioned parties, it was
impossible to complete this final year project.
Executive Summary
In this thesis we will report the design and fabrication of an aeration
experimental unit used to calculate the oxygen gas-liquid volumetric
mass transfer coefficient (kLa). The total cost of manufacturing the unit
was almost 500 O.R and was totally designed and arranged by the
students. The mixing tank designed to have a total volume of with
dimensions . Distilled water at atmospheric pressure and
room temperature was used to determine the value of kLa.
Table of Contents
Declaration......................................................................................................................2
Acknowledgment..................................................................................................................3
Executive Summary..............................................................................................................4
Nomenclature.......................................................................................................................6
Chapter 1 Introduction.........................................................................................................8
Chapter 2: Literature Survey............................................................................................10
2.1 Methods for measuring kLa:.....................................................................................11
2.2 Factors that affect kLa value:................................................................................13
2.2.1 Effect of salt addition and viscosity:.................................................................13
2.2.2 Effect of probe position:...................................................................................18
2.2.3 Effect of aeration rate:......................................................................................19
2.2.4 Effect of Temperature:......................................................................................22
2.2.5 Effect of medium depth:...................................................................................23
2.2.6 Effect of mixing:...............................................................................................24
2.2.6.1 Effect of impeller Type:.............................................................................24
2.2.6.2 The effect of impeller speed:.....................................................................28
2.2.6.3 Effect of impellers position:......................................................................28
Chapter 3: Theory.............................................................................................................30
3.1 Determination of kLa from experiments:.................................................................30
3.2 Determination of kLa from empirical correlations:..................................................31
REFERENCES..................................................................................................................35
Nomenclature
a
C*
kL
k La
NP
Power (W)
Qg
to
Initial time
VG
VL
Weight
Chapter 1 Introduction
are
leaching
of
metal
concentrates
and
microbiological
In this thesis, the affect of different impeller types and speed, gas flow rate, liquid
viscosity, and addition of salts on kLa will be studied and observed. At the end, an
empirical equation will be determined for possible causes of study.
Surface aeration uses a pump which is placed just above the pond
surface to circulate water from the pond surface into the air or it can
be just right at the surface. Surface aerators are preferred in shallow
ponds but it doesn't suits deep lakes.
Disadvantages
a compressor.
using the film model of mass transfer between the two phases by assuming that the
resistance at the gas phase is negligible. There are several methods that are used for
measuring kLa and are described shortly below:
low
time
requirement
and
the
results
obtained
dissolved oxygen and time and the second method requires the flow rate and the
oxygen concentration of the inlet and outlet gas. It is very important that the
reaction rate and the mass transfer rate are in right proportional to each other. If
the reaction rate is low compared to the mass transfer rate than the process is
controlled by the kinetics of the reaction, but if the reaction rate is very high than
and
experiments
have
shown
that
k La
values
are
significantly dependent on particle size of the soil and the clay content,
but are not affected noticeably with the concentration of soil organic
matter. As the soil content increases above 40% the k La value
decreases to about 60 to 70% of the value of k La in pure water. It was
also noticed that the oxygen supply is limited, diameter and the
number of bubbles also decreases if the clay content is very high [ ].
Robinson et al. (1973) and Linek et al. (1987) have studied the effect of
medium properties on the value of k La. They have found that adding
small amounts of solute to distilled water will cause reduction in bubble
size even if other properties were kept constant. They also reported in
their research paper on kLa values in electrolytic mediums that for
electrolytic solutions the decrease in bubble size is due to the electrical
effect of the concentration gradient of the ionic species between the
gas and liquid interface. This makes the medium capable of preventing
coalescence of bubbles so an increase in the value of k La is noticed
[ two refrences ].
Juarez and Orejas (2001) reported the effect of medium properties on
the value of kLa by considering two different liquids which were distilled
water and 0.5M Na2SO4 solution. The distilled water is known to act as
a coalescing system and the aqueous solution acts as a non-coalescing
system. The operating conditions used could be summarized in the
below table:
Operating variable
Liquid volume
Used Value
0.0015 m3
Temperature
0.002 m3
20C
if one impeller
if two impellers
Rotation speed
Gas flow rate
5 to 15 r.p.s
Varied between 2.5, 3.75 and 5
Experiment time
10-5 m3/s
200 to 700 s
Although water is a
coalescing system and showed low values of k La when using only one
impeller, but increasing the number of impellers to two will help to
break down the bubble size and cause an increases in kLa. On the other
hand, in sulphate solution the value of k La obtained is not affected if
one or two impellers are used because the solution prevents collection
of bubbles.
Puthli et al. (2004) experimented the effect of liquid viscosity on k La
values using three different concentrations of carboxymethyl cellulose
(CMC) which were 0.25, 0.5 and 0.75% w/v. The actual volume used
was 1.81m3. The effect of changing the viscosity was estimated at an
impeller speed of 600 rpm and gas flow rate range from (0.9-3.4)*10 -5
m3/s. The following table and graph shows the results obtained [&]:
Concentration of
Viscosity of
Value of kLa
CMC
solution
range 10
(s-1)
%w/v
0.25
mPa
4.33
-5
(m3/s)
0.9
0.0082
0.5
5.99
0.75
11
3.4
0.9
3.4
0.9
3.4
0.0161
0.0058
0.0095
0.0021
0.0037
Table 2.3: The effect of medium velocity on the kLa value obtained by
Puthli et al [].
Figure 2.2: Effect of viscosity at 600 rpm at varying gas flow rate
(triple impeller () water, () 0.25% CMC (w/v), () 0.375% CMC (w/v),
() 0.5%CMC (w/v)) obtained by Puthli et al [].
From the results obtained the writes concluded that increasing the
viscosity of the medium will reduce the k La values regardless the speed
of impeller or the gas flow rate used. This behaviour is probably
because of the decreases in the surface area of the bubbles due to the
viscous forces generated in the fluid.
Ozbek and Gayik (2000) have studied the effect of adding glycerol to
distilled water with the following composition 10,20,30,40 and 50%
glycerol solution w/w. A volume of 0.6 L at a temperature of 37C and a
pH of 7 was used .The operation condition used in the experiments and
the values obtained is described and summarized as following:
Agitatio
Air flow
Range of composition
Viscosity
kLa value
n speed
rate
of glycerol solution
(cpoise)
(min-1)
(rpm)
300
(L/min)
0.3
(% w/w)
10
50
0.935
6.948
2.65
1.4
Table 2.4: The effect of changing the viscosity of water by Ozbek and
Gayik [].
Figure 2.3: The kLa values versus the viscosity of the medium by
Ozbek and Gayik [].
The graph proves the truth that any increase in the medium viscosity
will decrease the value of kLa. The authors than investigated what
happens if medium viscosity is changed under the same operating
condition, but also the agitation speed and the flow rate of air is
varied. The following graphs were obtained when changing both
impeller speed and aeration rate:
Figure 2.4: The effect of stirrer speed on kLa in glycerol solution by Ozbek and
Gayik [O].
Figure 2.5: The kLa value verses air flow rates in glycerol solution by Ozbek and
Gayik [O].
It is noticed that although increasing the viscosity will decreases the value of k La, but
increasing the agitation rate or the gas flow rate could reduce this effect as it is noticed
from the graphs. It is also noticed from the figure that after 300 rpm the value of kLa
increased linearly with increasing the agitation rate.
conditions
Room
temperature
Atmospheric
pressure
gas 0.3 m/s
Superficial
velocity
Gas hold up
Bubble diameter
0.009
7.5 mm
2 Air diffuser. The air diffuser used was bought from a pet store and
had a diameter of 7.6 cm. However the writes suggested using
special designed air diffusers for further experiments.
3 Dissolved oxygen sensors. Three sensors were used in the
experiments, but they were not calibrated. The author
recommended that for better and accurate results the sensors
should be calibrated. A Vernier dissolved sensor was used that
contains a platinum cathode and Ag/AgCl reference anode dipped in
KCl electrolyte and are separated from the surrounding sample
solution by an oxygen permeable membrane. A fixed voltage is
applied to the platinum electrode and the following reactions take
place:
0.5 O + HO + 2e
2OH
reduction reaction
Ag + Cl
AgCl + e
oxidation reaction
As a result of the reactions, an electrical current flow is generated
and measures the concentration. The current is converted to
voltage, amplified, and recorded.
4 Other equipments such as, sensor seal, air compressor, nitrogen
cylinders, and personal computer were all used.
To determine the volumetric mass transfer coefficient first of all, any
amounts of dissolved oxygen in the water was removed by supplying
nitrogen gas until the concentration of oxygen felled below 1%. This
step usually took 30 min. After that, oxygen was introduced to the
column as compressed air through the diffuser, and the dissolved
oxygen was measured every second until water became saturated.
This process usually took less than 20 min. The experiments were
performed three times and the re-oxygenation profile for each
experiment was analyzed individually.
The author noticed that the re-oxygenation plots were almost linear
and the estimated mass transfer coefficient for the three sensors was
almost 0.31 0.01 min1. This observation proves that k La is a lumped
parameter as it didnt depend on the position of the sensors.
kLa value
Triple impellers
(rpm)
300
applied (cc/s)
8.92
33.81
8.92
33.81
8.92
33.81
8.92
33.81
(s-1)
0.0026
0.0035
0.0045
0.0081
0.0062
0.0098
0.0082
0.0161
400
500
600
Table 2.6: The experimental values of kLa obtained by Puthli et al. using
different speeds of impeller and air flow rate [].
Type of
Speed of
kLa value
(L)
impeller
impeller
Range
(s-1)
Six blades
(rpm)
500
(L/min)
0.5
6
0.02
0.12
Table 2.7: The operating conditions and the effect of air flow rate of kLa used by
Jacjson.
Figure 2.7: kLa versus Aeration Rate in distilled water by Eric Jackson [O].
Ozbek and Gayik (2000) have also studied the effect of aeration rate
on kLa values in a system containing distilled water. They have also proved that
increasing the aeration rate will cause an increase in the value of k La. The operation
conditions used in the experiments performed and the results obtained are as follows:
Volume
Temperature
(L)
(C)
pH
Agitation rate
Range of
kLa value
(rpm)
flow rate of
(min-1)
air
0.6
37
300
(L/min)
0.15
0.9
1.728
5.35
Table 2.8: The experimental values of kLa obtained using different air flow
rates for a distilled water system by Ozbek and Gayik [O].
Figure 2.8: The kLa values versus air flow rate in distilled water system obtained by
Jackson [O].
In conclusion it is noticed that the k La value is a strong function of the aeration rate with
a direct proportional relationship between them. Any increase in the gas flow rate will
cause on increases in the number of bubbles and the interfacial area, so k La value
increases.
water
temperatures
between
27.60
to
29.0
C.
Depth to
diameter ratio
in this study
H/T
shapes.
shows
(m)
Superficial gas
the
velocity (m/s)
of kLa value
of
the
0.392
(vvm)
it is noticed
depend on the type of impeller.
Geometries
condition
Vasconcelos.
turbine
had
0.013
used
several
Figure
2.8
dependence
0.125
on the shape
impeller
and
Figure 2.11: kLa results for the extreme values of constant air flow
rate QG, by Vasconcelos et al. Full lines represent the overall
correlation. Dashed lines stand for 95% confidence interval [].
Nelson et al. (1998) reported that Rushton turbine is used in most
industries but according to the limitations discussed previously which
leads to power drop due to the gas cavitations near the blades. In this
experiment, hollow blade impellers which are concave in shape, a tank
of 0.4m diameter made of clear Perspex to notice the gas dispersion,
and baffles were used. It was noticed that hollow blade impellers can
eliminate the reduction in the power demand, provide more handling
capacity and gas hold up capacity which makes them better than
Rushton turbine.
The effect of using single, double and triple impellers was studied by Puthli et al. (2005).
The impellers used in this study were single impeller contisting of a disc turbine (DT),
double impeller consisting of disc turbine-pitched blade downflow turbine (PTD), triple
impeller consisting of disc turbine-pitched blade downflow turbine- pitched blade
downflow turbine. The reactor used in this study was concentrated of glass with baffles
and it has the following dimensions:
Moreover, the impeller was located 4.5 cm from the bottom of the reactor. The ratio of
the diameter of the impeller into the diameter of the tank was 0.333 in the three cases.
Impeller
blades
Impeller
Impeller
configuration
diameter
blade
width (cm)
number
Single
(cm)
4.33
length (cm)
1.4
1.2
(Np)
4.8
Six blade DT
4.33
450 PTD (4 4.33
1.4
1.8
1.2
1.4
6.3
blades)
Six blade DT
4.33
0
45
PTD (4 4.33
1.4
1.8
1.2
1.4
7.8
blades)
450 PTD
1.8
1.4
Six blade DT
impeller
Double
impellers
Triple
impellers
(4 4.33
blades)
Table 2.10: Impellers configuration used by Puthli et al [].
During the study it was noticed that using single, double and triple impellers had an effect
on the kLa vakues at different impeller speed and gas flow rates. Distilled water was used
in this experiment at a temperature of 35 C0. With an impeller speed range of 300 to 600
rpm and the rang of the gas flow rate was 8.92-33.8 cc/s the results of kLa value at
600rpm and 33.8 cc/s are shown in the following table:
Type of impeller
Single impeller
Impeller speed
kLa value
(rpm)
(cc/s)
(s-1)
0.0073
Double impeller
Triple impeller
600
33.18
0.0098
0.0161
Table 2.11: kLa value of single, double, and triple impeller obtained by Puthli et al. [].
figure 2.10, consists of two regions first one, when the impeller speed
was lower than 300 rpm the slope was very small compared with
second region at impeller speed above 300 rpm where the k La value
increased noticeably.
Figure 2.13: the effect of impeller speed on the kLa value obtained by
Ozbek and Gayik [].
250C
0.202m
1mm
Distilled water
air
Chapter 3: Theory
= OTR OUR
(3.1)
OTR is the oxygen transfer rate (mol O2/m3.s) and is:
Where
(3.2)
= kLa (C* - C)
(3.3)
Separating the variables and integrating from C0 to C and t0 to t:
CC
ln ( CCo ) = - kLa (t t0)
(3.4)
CC
The plot between ln ( CCo ) and t will result in a line with a slope of - kLa as in
the following example graph:
Many developed correlations apply dimensionless quantities such as, Reynolds (Re). The
impeller Reynolds number represents the ratio of the inertial forces created by the
rotating impeller and the viscous forces of the liquid:
ReN =
l . N . D2
(3.5)
When the impeller Reynolds number is less than 10 4 then the liquid is in the laminar
region; while in the turbulent region, the impeller Reynolds number is more than 10 4.
Usually the mechanical stirred vessels operate in turbulent region, where the velocity
fluctuates in both time and three dimensions in space. The table below shows other
dimensionless groups used in kLa correlations:
Dimensionless quantities
Froude Number (Fr)
gas flow rate (FlG)
Sherwood (Sh)
Stanton (St)
Schmidt (Sc)
Weber (We)
Aeration (A)
Equation
2
N D
Fr =
g
FlG
Qg
= N . D3
Sh =
k La D
T
St =
k LaV L
Qg
Sc =
We =
A=
N D
ND
VG
from the bottom of the vessel, and the spacing between impellers if multi impellers are
used. The power could be calculated by, where NP is found from graph 3.2:
P = NP. l. N3. D5
(3.3)
P=KNS
(3.4)
Where
In the table below different correlations from literatures are reported. Some depends on
the gassing rate and the mixing power of the impellers, others are based on dimensionless
groups:
Researchers
Calderbank (1958) van 't Riet (1979)
Hickman (1988)
Correlation Proposed
kla =0.026(P/VL)0.4(VG)0.5
For T = 0.60 m
kLa = 0.043 (P/VL)0.4(VG)0.57
For T = 2 m
kLa = 0.027(P/VL)0.54(VG)0.68
kLa = 4.95 10-3 (P/VL)0.593(VG)0.4
kLa =3.84 10-3 (P/ VL)0.654(VG)0.4
kLa =14.6 10-3 (P/ VL)0.611(VG)0.554
kLa = (P/ VL)0.4(VG)0.35
kLa =1.25 10-4 (D/T)2.8(Fr)0.6(ReN)0.7
(FlG)0.45(D/g)-0.5
kla = 0.01(P/VL)0.475(VG)0.4
kLa =10.4 10-3 (P/ VL)0.475(VG)0.4
kLa = 5.3 10-3 (P/ VL)0.55(VG)0.32
kLa =26 10-3 (P/ VL)0.4(VG)0.5
kLa =0.57 (P/m)0.4(VG)0.55
kLa =0.031 (P/ VL)0.4(VG)0.5
kLa =1.11 10-3 (P/ VL)0.9504(VG)0.6282
kLa =10.36 10-3 (P/ VL)0.67(VG)0.56
kLa=21.2 (Re)1.11(Sc)0.5 (Vg.T/)0.45(G/a) (DL/T2)
kLa=0.06 (Sc)0.5(Re)1.5 (Vg.a/)0.6 (Fr)0.19 (A) (DL/T2)
kLa=0.368 (Re)1.38 (Sc)0.5 (Vg.a/) (Fr)0.367 (A)0.167(T/D)0.25
((P/V)/N3T5)0.75 (DL/D2)
kLa=1.41 10-3 (Re)0.67 (Sc)0.5 (We)1.29 (DL/T2)
kLa= C. [(P/VL)/ (v.g4)1/3] [(Q/VL). (v/g2)1/3]0.23 (g2/ v)
Table 3.2: Correlations proposed by different researchers to estimate the value of kLa in
stirred tank reactor [three ref].
REFERENCES
1. Gollakota, S.V. and Guln, J.A. (1984), Comparative Study of Gas-Liquid Mass
Transfer Coefficient in Stirred Autoclaves, Tubing Bomb Microreactors, and
Bubble Columns, American Chemical Society, Alabama.
2. Hashsham, A.S. (2007), Design of an Experimental Unit for the Determination of
Oxygen Gas-Liquid Volumetric Mass Transfer Coefficients using the Dynamic Reoxygenation Method, Michigan State University, East Lansing.
DETERMINATION
OF
VOLUMETRIC
MASS
TRANSFER
15. Pouliot, K. Thibault, J. Garnier, A. and Leiva, G.A. (2000), KLa Evaluation
during the course of fermentation using data reconciliation techniques,
Bioprocess Engineering 23, Canada and Chile.
16. Puthli, M.S. Rathod, V.K. and Pandit, A.B. (2005), Gas-liquid mass transfer
studies with triple impeller system on a laboratory scale bioreactor, Biochemical
Engineering Journal 23, India.
17. Robinson, C.W. and Wilke, C.R. (1973), Oxygen absorption in stirred tanks: A
correlation for ionic strength effects, Biotechnology and bioengineering.
18. Scargiali, F. Grisafi, F. and Brucato, A. Analysis of the differences in kLa values
determined by different variants of the dynamic method in stirred tanks,
University of Palermo, Italy.
19. Sorenson, K.L. (2010), COMPARATIVE STUDIES ON OXYGEN MASS
TRANSFER FOR THE DESIGN AND DEVELOPMENT OF A SINGLE-USE
FERMENTOR, Utah State University, Logan, Utah.
20.
Tervasmaki, P. (2012), Measurements of Volumetric Mass
Transfer Coefficient in a Gas Liquid Stirred Vessel, Department
of Process and environmental Engineering Bachelors Thesis 230.
21. Vasconcelos, J.M.T. Orvalho, S.C.P. Rodrigues, A.M.A.F. and Alves, S.S. (2000),
Effect of Blade Shape on the Performance of Six-Bladed Disc Turbine Impellers,
American Chemical Society, Portugal.
22. Wang, S.J. and Zhong, J.J. (1996), A Novel Centrifugal Impeller Bioreactor. II.
Oxygen Transfer and Power Consumption, Biotechnology and Bioengineering,
Volume 51, Page 520-527, China.
23. Woo, S.H. and Park, J.M. (1997), Estimation of oxygen transfer in soil slurry
bioreactor, Volume 11, Biotechnology Techniques, Korea.
24. Wu, H. (1995), AN ISSUE ON APPLICATION OF DISC TURBINE FOR GASLIQUID MASS TRANSFER, Volume 50, Elsevier Science Ltd, Italy.
25. Wu, H. Arcella, V. and Malavasi, M. (1997), A study of gas-liquid mass transfer
in reactors with two disc turbines, Volume 53, Elsevier Science Ltd, Italy.
26. Yawalkar, A.A. Heesink, A.B.M. Versteeg, G.F. and Pangarkar, V.G. (2002), GasLiquid Mass Transfer Coefficient in Stirred Tank Reactors, Volume 80, The
Canadian Journal of Chemical Engineering, India and Netherlands.