Issa PDF
Issa PDF
Issa PDF
&OWVFEFMPCUFOUJPOEV
%0$503"5%&-6/*7&34*5%&506-064&
%MJWSQBS
Institut National Polytechnique de Toulouse (INP Toulouse)
1STFOUFFUTPVUFOVFQBS
Hayder Mohammed ISSA
24 Octobre 2013
le
5JUSF
6OJUEFSFDIFSDIF
Laboratoire de Gnie Chimique (LGC), Toulouse
%JSFDUFVS T
EFTF
Martine POUX, Ingnieur de Recherches (HDR), INP-ENSIACET, Toulouse
Catherine XUEREB, Directrice de Recherches CNRS-LGC/INPT, Toulouse
Jury :
Cathy CASTELAIN, Directrice de Recherches CNRS-LTN/INP, Nantes, Examinateur
Michel SARDIN, Professeur, ENSIC/INPL, Nancy, Rapporteur et Examinateur
Denis BOUYER, Professeur, Universit de Montpellier II, Rapporteur et Examinateur
Jean-Pierre GRASA, Prsident-Directeur Gnral , Biotrade, Toulouse, Examinateur
Remerciements et Ddicace
Je voudrais tout dabord remercier mes directrices de thse Catherine Xuereb et
Martine Poux, pour la confiance, le temps et la patience quelles mont accord au
long de ces trois annes, ainsi que pour leurs conseils aviss qu'elles ont port au
projet.
Jadresse ma reconnaissance M. Michel Sardin, Professeur INPL Nancy et M.
Denis Bouyer, Professeur lUniversit de Montpelliers II, qui ont accept de juger
ce travail et den tre les rapporteurs.
Je remercie galement Mme. Cathy Castelain, professeur lUniversit de Nantes
davoir accept de prsider le jury de thse.
Je remercie M. Jean-Pierre Grasa, PDG Biotrade pour son intrt ce travail, pour
son aide et ses conseils aviss et pour accepter de juger ce travail.
Mes remerciements sadressent Jolle Aubin et Karine Loubiere pour leur soutien
quelles mont fourni tout au long de ce travail.
Je souhaite remercier lquipe technique du LGC pour leur aide, tout particulirement,
Jacques Labadie, Lahcen Farhi et Alain Muller, pour leur disponibilit, et pour
leur aide dans la mise au point exprimentale. Mes remerciements s'adressent
galement aussi au service administratif du LGC particulirement Danile
Bouscary pour sa disponibilit. Je remercie tous mes collgues du laboratoire LGC.
Je ddie ce travail la mmoire de ma chre mre. Je ddie ce travail aussi ma
chre pouse Aseel, merci pour ton soutien et mes enfants Abdullah, Danya
et Mohammed, qui mattendent avec impatience.
Enfin, je tiens remercier tous ceux qui, de prs ou de loin, ont contribu la
ralisation de ce projet.
II
Abstract:
A new surface aeration system for water and wastewater treatment has been studied.
Its uniqueness lies in its ability to operate in two modes: aeration or simply blending
(mixing) by just reversing the direction of rotation. An experimental plant has enabled
to focus on mass transfer performance and hydrodynamics. The flow pattern and the
velocity field measurements inside the agitated tank were performed by both the Laser
Doppler Velocimetry (LDV) and the Particle Image Velocimetry (PIV) techniques for
the single phase (Mixing) mode and for the two phases (Aeration) mode. The oxygen
mass transfer occurs both in the water bulk and in the spray above water surface and
has been independently investigated. Different configurations and operational
conditions were tested during the experimental part in order to interpret phenomenon
effect of the draft tube and RTP propeller, rotational speed, turbine blades
submergence and else on the flow field and the oxygen mass transfer in the agitated
system that produced mainly by a cone shape turbine. The experimental part dealing
with hydrodynamics and flow field shows that the down-pumping operation mode
with the draft tube has the most convenient results in the mixing mode with respect to
turbulent flow field and mixing time. Whilst for the up-pumping aeration mode the
hydrodynamics experimental results show the whole system configuration is the most
convenient with regarded to mean velocities, turbulent flow intensity and mixing
time. For the oxygen mass transfer experimental part, it is found that the highest
standard liquid bulk aeration efficiency is achieved (SAEb = 2.65 kgO2 kw-1h-1) when
the whole system configuration is used. The highest standard aeration efficiency at 20
o
C for the water spray zone is accomplished ((Esp)20 = 51.3 %) with the whole system
configuration. Several correlations models have been derived for the oxygen mass
transfer in water bulk and spray zones, power consumption and mixing time, on the
basis of experimental results. They can be used as tools to estimate these parameters
for geometrical and dynamical similar systems at industrial scales.
Keywords: Surface aeration, Agitated tank, Mass Transfer, Oxygen mass transfer,
Hydrodynamics, Multiphase Flow, LDV, PIV, Dimensional analysis, Modelling
III
Rsum :
Un nouveau systme daration de surface pour le traitement des eaux uses a t
tudi. Sa spcificit rside dans sa capacit fonctionner selon deux modes : aration
ou simple brassage, en modifiant uniquement le sens de rotation du systme. Un
pilote a permis de cibler le travail sur ltude exprimentale du transfert de matire et
de lhydrodynamique.
Les champs d'coulement et les mesures de vitesse l'intrieur de la cuve agite ont
t raliss par vlocimtrie laser effet Doppler (LDV) et par vlocimtrie par
images des particules (PIV) pour le mode monophasique (brassage) et pour le mode
diphasique (aration). Le transfert d'oxygne se produit la fois dans la cuve et dans
le spray au-dessus de la surface de l'eau. Il a t tudi dans les deux zones.
Diffrentes configurations et conditions opratoires ont t testes
afin de
comprendre les phnomnes dinteraction : tube de guidage, hlice complmentaire
RTP, vitesse de rotation, niveau de submersion des pales de la turbine. La partie
exprimentale sur lhydrodynamique et les champs d'coulement montre que le mode
de fonctionnement en pompage vers le bas (brassage) avec tube de guidage procure
les meilleurs rsultats en termes de mlange si on se rfre aux champs d'coulement
et la mesure du temps de mlange. Pour le mode de fonctionnement en pompage
vers le haut (aration), les rsultats exprimentaux montrent que la configuration du
systme complet est la plus efficace si on considre le transfert doxygne, les vitesses
moyennes, l'intensit de l'coulement turbulent et le temps de mlange. Il est constat
que la meilleure efficacit d'aration standard est atteinte (SAEb = 2.65 kgO2kw-1h-1)
lorsque le systme complet est utilis. L'efficacit d'aration standard 20C la plus
leve au niveau du spray d'eau est obtenue ((ESP)20 = 51,3%) avec la configuration du
systme complet.
Plusieurs modles sont proposs pour calculer le transfert d'oxygne dans la cuve et
dans le spray, la consommation nergique et le temps de mlange. Ces relations
permettent dvaluer linfluence des diffrents paramtres gomtriques et de
fonctionnement dans des systmes similaires une chelle industrielle.
IV
Table of Contents
10
I. Porous Diffusers
10
10
11
11
11
11
11
11
11
11
12
13
14
14
15
16
16
16
16
16
18
V
18
19
I. Principals
19
23
24
27
28
28
29
29
30
31
A. Tank Geometry
31
B. Baffles Effect
32
33
34
36
36
37
V. Hydrodynamics
38
38
38
B. Vortex Formation
40
40
41
41
42
VI
43
43
44
44
45
46
1.3. Conclusions
46
51
2.1. Introduction
51
51
54
54
54
54
55
55
58
58
I. Theory
58
59
59
60
V. Laser Source
61
61
62
62
VII
63
63
64
65
66
67
67
69
69
I.1. Introduction
69
71
71
72
73
73
I.7. Determination Model of the Bulk Zone Oxygen Mass Transfer Coefficient
75
I.8. Measurement Procedure for the Bulk Zone Oxygen Mass Transfer Coefficient
76
77
I.10. Oxygen Transfer Rate for the Bulk Mass Transfer Zone (OTRb)
77
I.11. Standard Oxygen Transfer Rate for the Bulk Mass Transfer Zone (SOTRb)
77
I.12. Standard Aeration Efficiency for the Bulk Mass Transfer Zone (SAEb)
78
78
II.1. Introduction
78
79
II.3. Determination Model for the Spray Zone Oxygen Mass Transfer Coefficient
79
81
II.5. Oxygen Transfer Rate in the Spray Mass Transfer Zone (OTRsp)
82
VIII
82
83
84
2.7. Conclusions
87
91
91
3.1.1. Introduction
91
92
92
93
95
97
98
VI. Standard Aeration Efficiency (SAEb) and Standard Oxygen Transfer Rate (SOTRb) for the
Water Bulk Zone
101
106
I. Mass Transfer
106
108
110
3.2.1. Introduction
110
110
110
111
115
I.3. Surface Aeration Water Spray Discharge Velocity and Volumetric Flow Rate
117
118
IX
I.5. Contribution Percentage of the Spray and Bulk Zones in the Overall Mass Transfer
Operation
119
121
121
122
126
127
II.5. Contribution Percentage of the Spray and Bulk Zones in the Overall Mass Transfer
Operation
129
III. Effect of Propeller and the Draft Tube
130
130
133
134
135
IV. Comparing the OTRsp for the Whole System and Turbine Alone Configurations
136
137
3.2.4. Conclusions
142
147
147
150
150
150
158
162
162
165
X
168
170
4.5. Agitation Index and Liquid Volume Quantification for the Down-Pumping
Mode with Draft Tube Configuration
171
172
166
4.7. Conclusions
178
183
5.2. Flow Pattern and Mean Velocity Field in the Aerated Tank
183
185
5.3.1. Flow Pattern and Mean Velocity Field in the Aeration Tank
185
186
5.4.1. Flow Pattern and Mean Velocity Field in the Aeration Tank
186
199
5.6. Agitation Index and Liquid Quantification for the Whole System
Configuration
199
200
203
204
205
206
207
209
5.9. Conclusions
211
XI
215
223
225
225
226
226
226
228
A. Nitrogen Removal
228
B. Phosphorus Removal
229
231
232
232
233
237
List of Symbols
241
List of Figures
247
List of Tables
257
References
259
XII
-1-
-2-
Introduction
Clean water is growingly demanded in the different fields of human activities, for
example water is more and more being used in the industry (Roubaty and Boeglin,
2007). On the other hand the continuous diminution of existing water resources made
water and wastewater treatment to become a truly developing and problematic
question. One of the ways used to maintain clean water resources for the diverse
industrial or urban demands, is wastewater treatment.
The implementation of the surface aeration process in the water and wastewater
treatment is established as an effective treatment for various wastewater types
especially in activated sludge biological and aerobic water treatment processes. This
technology has an important capacity of delivering the needed oxygen to the aerobic
micro-organisms for respiration and ensures efficient mixed condition for the entire
treatment tank through maintaining the microbial flocs in continuous state of agitated
suspension by accompanied mixing in order to achieve maximum contact surface area
between the flocs and wastewater (Gary, 2004). Surface aeration has various desirable
characteristics such minimum sludge residual is produced for the used activated
sludge process as a continuous operation of recycling the used sludge is implemented
for wastewater treatment plant (Nair et al., 2008; Ramalho, 1977).
Taking in to consideration the capacity of now used surface aerators in the water
treatment field with respect to the accomplished aeration efficiency, energy
consumption and the complicated maintenance as described in the related works.
It is expected from this work to provide the necessary investigations to prove the
flexibility and capability of the purposed innovative surface aerator to work in two
ways, aeration and mixing by simply reversing the sense of the rotation and acting on
a clutching system. The novel investigated surface aerator (FR Patent Demand, 2012)
is found useful and promising after the new method for combination between the
delivering the necessary oxygen into the water treatment tank by the up-pumping
aeration mode and to achieve an efficient mixing for the treatment tank constituents
by the down-pumping mixing mode.
Implementing the new surface aeration technology enables easy maintenance and
more energy saving. In addition the operation may be performed with minimum cost.
This research can be regarded as an approach that may open new opportunities to
enhance the aeration efficiency with optimized operation condition.
During the last decades, the surface aeration is considered as an effective oxidation
mean among the existed wastewater treatments tools. It consists in the dispersion of
the waste water into droplets and there projection through the atmospheric air, where
-3-
-4-
Outlines
The main objective of this work is to characterize the performance of the surface
aeration system in aeration mode and its blending capacity in the mixing mode. This
will be fulfilled by identifying the affecting parameters on the oxygen mass transfer
process that developed within the liquid bulk inside the tank and in the spray at the
water surface. These parameters are the operation conditions and the geometrical
configurations such as; the rotational speed, mixing time, impellers configurations and
the power consumption.
The other main objective of this work is to acquire the flow behavior for the two
phase (gas-liquid) condition (Aerated mode) and for the single phase condition
(Mixing mode) with related power consumption and impellers configurations.
-6-
-7-
-8-
Chapter One
Surface Aeration Processes for Water Treatment
-8-
I. Porous Diffusers
They also in turn have many types as follows:
I.1. Plate Diffusers: These are usually having a form of (30 cm) square surface and
(25-38 mm) thick; most are constructed of ceramic media or made of porous plastic
media of (30 cm x 61 cm) surface area. Air is introduced below the plates through a
plenum (See Fig. 1.2).
I.2. Panel Diffuser: These usually employ a plastic membrane which is stretched over
(122 cm) wide placed on base material of reinforced cement compound or fiber
reinforced plastic.
I.3. Tube Diffuser: These are constructed from stainless steel or a durable plastics,
they have generally (51-61 cm) long with (6.4 7.7 cm) diameter (Figure 1.2,a).
I.4. Dome Diffusers: They are usually used with the dimensions (18 cm) in diameter
and (38mm) high, the used medium is usually is ceramic materials.
I.5. Disc Diffuser: they are relatively flat but they differ in size, shape, method of
attachment and kinds of diffuser materials. Generally they have configuration of (1851 cm) diameter (Figure 2.13, b).
(a)
(b)
Figure (1. 2) : (a) Tube diffuser aerator, (b) Disc diffuser aerator, (Gemgate GmbH)
- 00 -
For the diffused aeration system with deep tanks between (4-6m), a combination of
aeration systems is used to the water treatment plants in order to improve the transfer
of oxygen. Usually one of applied means is using turbo-compressors to increase the
flowrate of injected air.
- 02 -
- 03 -
Figure (1.6): Flow diagram for uncovered pure oxygen aeration (Mueller et al., 2002).
- 04 -
Figure (1.8): Low speed vertical flow aerator (Up-ward flow) (Praxair Technology).
(a)
(b)
Figure (1.9): Low speed vertical flow aerator (Down-ward flow), (a) Turboxal (Aire
Liquide), (b) Praxair (Praxair Technology).
- 06 -
Figure (1.10): Low speed horizontal flow aerator (Twin mini rotor aeration,
(Botjheng Water Ltd.).
- 07 -
Figure (1.11): High speed surface aerator, Aqua turbo (AER-AS), (AQUATURBO
SYTEMS inc.) .
- 11 -
Figure (1.12): The surface aeration regimes applying air entrainment from free liquid
surface (A) Direct entraining of the atmospheric air, (B) Spray formation and
entraining the air with droplets impingement at the surface,(Patwardhan and Joshi,
1998).
Various geometries are developed for surface aeration turbines, where the design is
emphasized to project highest quantity of water droplets into the atmospheric air by
achieving the largest contact area between the two phases. The limiting factor for the
surface aeration turbines is always taking in account reducing the radial discharge
flow of the droplets toward water surface level. The circulation of the water and
comprised air bubbles in the treatment tanks is generally consist of a main loop in
entire the tank, many secondary loops can be generated depending on many system
specifications like the number and type of impellers, tank geometry and impeller
position in addition to many system characterizations such as air bubbles hold-up and
retention time and else. For instance when the impeller is placed in the water bulk, the
secondary loops may appear around the impellers, or they can be developed in the
upper or lower part of the tank. When surface aeration turbine is positioned at the
water surface, the water bulk is usually engaged with one main circulation loop.
Deeper treatment tanks are generally preferred to ensure the needed residence time of
the air bubbles but on the other hand these tanks need sophisticated tools to let air
bubble reach the bottom of the tanks (Jakobson, 2008; McCabe et al., 1985; Nagata,
1975; Roustan, 2003; Tatterson, 1994; Xuereb et al., 2006) .
The circulation or mixing time is considered as a measurement indicator of the
average water bulk motion that generated by the impeller in tank. Measurements of
circulation and mixing times are considered as an indicator to understand the scalar
transport in the tank (Edwards and Baker, 2001 ). The circulation time in the surface
aeration is generally associated with; tank overall flow rate, air entrainment flow and
impeller pumping capacity. For more complicated configurations with multiple
impellers, the mixing or circulation time behavior depends on the created circulation
resultant of these impellers, where its not evident always when the number of
impellers increased that leads to shorten the mixing time (Wang et al., 2010).
For continuous flow surface aeration systems that is the case with open channels, the
circulation time is related with; the impeller speed, bulk motion and convective
- 10 -
Table (1.1): The standard aeration efficiency (SAE) for various aerators types
SAE
(KgO2/kWh)
Aerator Type
High Oxygen Purity Aerator
(3.5-5.5)b
(2.1-2.55)h
Submerged Turbines(Axial)
(1.0-1.6)b
Submerged Turbines(Radial)
(1.1-1.5)b
2.50a
1.00a, 0.9e
(0.4-0.9)b, 0.42c, 1.6e
0.60a
b
c
d
(Duchene and Cotteux, 2002); (Mueller et al., 2002); (Kumar et al., 2010a); (Moulick and Mal, 2009);
e
f
g
h
(Cancino, 2004a); (Boyd, 1998); (Thakre et al., 2009);
(Taricska et al., 2009).
- 11 -
Table (1.2): The standard aeration efficiency of the fine bubble diffused aeration in
different water treatment basins (Duchene and Cotteux, 2002)
SAE
(KgO2/kWh)
Basin Type
Large Open Channels
3.41
1.95
3.11
2.12
- 13 -
(1.3)
Where, MLSS is the mixed liquor suspended solids (the water with containing
microorganisms); r is the specific respiration rate; z is the length of each stage and us
is the liquid velocity. The modified model takes in account the correction from water
conditions to mixed liquor. This general model proposes a relative improvement in
the oxygen mass transfer coefficient (kla), but it still didnt reach the required true
value depending on standard operation conditions. (Ju and Sundarajan, 1992) found in
their studies on the oxygen transfer in surface aerated bioreactors with containing
microorganisms that the presence of the biomass showed no effect on the oxygen
transfer rate to water because of their formed film adjacent to the gas-liquid film is
- 14 -
(1.4)
Where, A = Q Emd /VL; B = klsas; Cd*and CLS* are the saturated concentrations in the
droplet and water respectively and Emd is Murphree efficiency for the operation.
(Oliveira and Franca, 1998) modified the previously developed model for oxygen
transfer by (McWhirter et al., 1995) to represent the surface aeration for the turbines
that placed in water sub-surface position. They tried to simplify the model by
applying boundary initial condition and other considerations and assumptions to
eliminate the low effecting parameters, so they found the following model for oxygen
mass transfer coefficient with certain conditions:
dy/ dz = - At/G kla (CL* - CL) K2
(1.5)
With applying Henry's low and rearranging the equation the found following model
was:
C*L = C*st (1+yo/yo) [[Pb-Pv+L g (Zs-Z) / 1-Pv ] y/y+1
(1.6)
(1.7)
Where C*st is the standard DO level; Pb is the barometric pressure; Pv is the vapor
pressure; G is the gas flow rate; K2 conversion factor; At is the tank cross-sectional
area; C*L represents the true bulk liquid DO level; yo is the oxygen concentration in
the bubbles at (z = 0); and H is Henry's law constant.(Oliveira and Franca, 1998)
tested these models with the previous models experimentally, where they found that
acceptable fitting between the experimental and theoretical results, where the
dissolved oxygen at equilibrium state is decreased with increasing the temperature,
the highest oxygen transfer is noticed with low temperature values.
(Stukenberg et al., 1977) have studied the probability of the calculation errors with
several types of aeration equipment including the surface aeration. They studied the
- 15 -
Figure (1.13): The histogram distribution of measured aeration efficiency for (111)
low speed surface aerators in field, the average is 1.49 kgO2/kWh, (Heduit and
Racault, 1983b).
- 16 -
(1.8)
(1.10)
(1.11)
Where: is the blade tip angle. This model is applied within these ranges of Re
*103(50-84) and S/D (0.17-0.25).
- 18 -
OTR, (kgO2h-1)
(Backhurst et al., 1988) examined the effect of blade submergence on the oxygen
transfer rate for different impeller types. Their results showed that there is an
optimum submergence (starting with submergence an efficient aeration is noticed) for
all tested impellers as illustrated in the Fig. 1.14.
Figure (1.14): The relation between surface aeration impeller blades submergence and
oxygen transfer rate for different blades number, (where H is the liquid level in the
tank) (Backhurst et al., 1988).
The clearance of surface aerator impeller in the treatment tank is commonly defined
as the distance between the lowest point of surface aerator and the tank bottom. While
the submergence is defined as the distance the water surface level and the specific
point on surface aerator blade. (Patwardhan and Joshi, 1998) concluded that with
increasing the submergence of surface aerator impeller the intensity of surface
aeration, oxygen transfer rate and gas hold up are decreased. They explained that the
amount of energy reaching the liquid surface is decreased. That also agrees with the
results obtained by (Backhurst et al., 1988). It is better always to set the impellers in
closer position to the liquid surface to enhance the reached energy to the liquid
surface and increase the gas holdup. This persuasion was found by(Deshmukh and
Joshi, 2006) by testing three types of surface aerators impeller of PBTU, PBTD and
- 21 -
A. Tank Geometry
The surface aeration for water treatment is usually performed in cylindrical shape
tanks or basin, which are the most usual among the used tanks, but the geometry of
these cylindrical tanks may vary between plate bottom shape to curved and conical
shape (with 150 - 300 degree angle depending of the existed activated sludge
properties).The volume of the tank is commonly related with the height of the liquid
in the tank (Jakobson, 2008). Square or rectangular shape are also used but in very
limited way and for especial uses. On the whole, the water treatment tank volume that
equals height of the liquid is considered as standard geometric ratio for design
considerations.
(Rao and Kumar, 2007b) have derived a model for circular shape aeration tank, where
the derived model relates the mass transfer coefficient of oxygen, kla, with circular
- 20 -
(1.12)
The experimental results of (Rao and Kumar, 2007a) showed that the square tank
aeration was also effective in surface aeration process, where higher values of mass
transfer coefficient k are achieved in shorter duration (that doesn't agree with previous
woks where the square shape not preferred because the formation of dead angles
exist) but in the power requirement point of view the circular aeration tanks were
more effective for less amount of power was required to reach the value of mass
transfer coefficient with keeping the other conditions constant during the experiment.
For square tanks, (Rao and Kumar, 2007a) verified a correlation that was developed
earlier in the previous work (Rao, 1999), which represents the mass transfer and
power measurement on electrical measured basis. The general correlation found by
(Rao, 1999) was for the mass transfer parameter k (where k = kla20(2/g)1/3), as a
function of geometric and physical properties that referred as, X, the ratio of the
Froude Number Fr, to Reynolds No., Re, for a baffled tank:
k = [17.32 exp(-0.3/X 1.05)+3.68 -0.925 exp(-750 )X-0.057)2)]10-6 X
(1.13)
Where; X is (Fr4/3/Re1/3). This model is applicable within X range of (0.01 8.0) and
it is developed for the geometric and operational ranges; h/D =1.0, S/W=1.26 and
W/D = 0.24
(Fuchs et al., 1971) have studied the surface aeration performance by examining the
effect the volume of the aeration tank according to volumetric power provided to the
operation, where they tried to keep the mass transfer coefficient constant during the
tests. (Fuchs et al., 1971) founded that the oxygen mass transfer coefficient was
generally increased as the volume of aeration tank is decreased for large volume
tanks, a satisfactory results found for high levels of provided power per volume ratios
B. Baffles Effect
Baffles are commonly used in the water treatment tank, where they are fixed near the
walls of the tank to reduce or prevent the formation of vortex that are generated
because of the centrifugal force created by impeller rotation especially when
cylindrical vessels used and when the impellers are centrally positioned in the vessel
- 21 -
(1.14)
(1.15)
Where, PV = actual or measured power/ volume and k = kla20 (2/g)1/3. These models
are developed for the geometric and operational ranges; h/D =1.0, S/W=1.26 and W/D
= 0.24.
(1.16)
AE is the aeration efficiency (it is defined as the transferred oxygen mass rate to the
liquid per the power consumed at actual condition).
For pitched blade impellers
AE= (Q/P)-o.o2678 (D N)-5.7148(Re)-0.3388(Fr)4.8695(2)-0.3676(D/S)-0.1256(1)0.2
(1.17)
Where Q is the water flow splashed by the aerator; 1 is 1/Inlet angle of paddle (rad1);
2 is 1/Outlet angle of paddle (rad1). The model is developed for impeller blade
submergence (58% - 228%).
For paddle wheel surface aerators performance, (Moulick et al., 2002) established a
general model by finding a modified standard aeration efficiency termed SAE'. They
verified the model by performing several set of experiments with varying the
geometric configuration for each one with SAE' to reach the optimal geometric
parameter, where they varied the paddle width / impeller diameter ratio, liquid
volume/ impeller diameter, impeller pitch/ impeller diameter and horizontal projection
of bent length / impeller diameter and bent angle. The correlated general equation
was:
(1.18)
Where (SAE') = (SAE) (v/g2)1/3 (c)-1 . .N3. D2and X= Fr4/3/Re1/3. The ranges of this
model are; Re (2105 - 8105), Fr (0.05 0.25) and (S/D)paddle wheel (0.025-0.225).
(Zlokarnik, 1979) compared between different surface aerator types with respect to
the oxygen mass transfer. He found that the surface aeration is highly influenced by
several geometric ratios. The investigation he made was as a function of aerator type
and their geometrical specification with other operational parameters like Froude
- 24 -
(Patil et al., 2004) found out the surface aerator with the PBTD type was the best in
the general in comparison between the impellers for surface aerator; PBTD, PBTU
and DT impellers within the tested range of system configuration with respect to the
values of mass transfer coefficient achieved by these aerators.
(Deshmukh and Joshi, 2006) have tested the impeller geometry effects on the liquid
flow profile for the surface aeration system with impeller placed in water sub-surface
position. They tested different impeller types PBTU, PBTD and DT. The results
showed that the performance of PBTU was the worst between the impellers according
flow patterns in the tank with increasing the rotation speed, where the created flow
was unsatisfactory and didnt cover all parts of the tank in contrary to the rest two
impellers PBTD and DT.
Figure (1.15): Relation between sorption number (Y) with the Froude number (Fr) for
the conical shape turbine surface aerator (Zlokarnik, 1979).
(Backhurst et al., 1988) built a correlation for the OTR (The oxygen mass transfer rate
to the liquid during the aeration (kgO2/h)) with affecting factors like the geometric
parameters, the flow and water surface conditions in the surface aeration system:
OTR=10-3(DCs Dm) Re1.90 Fr0.15 (D/h)0.2 (S/h) (n/n8)0.20 (Tv/D)0.05
(1.19)
- 25 -
(1.21)
V. Hydrodynamics
V.1. Flow Patterns
A. Flow Patterns Characterization
The flow pattern and circulation inside of surface aeration in the wastewater tank has
a very crucial importance to ensure distribution and efficient mixing for the dissolved
oxygen, the homogenization and suspension of the mixed liquor that contains the
activated sludge (Roustan et al., 1984).
Most of the agitation and mixing means that used for the surface aeration are
operating with including impellers or turbines that are placed on a turning shaft. There
are many of these impellers that can be characterized according to the generated
pattern flow. These impellers accomplish two main flow types. The axial flow
- 27 -
B. Vortex Formation
The depth and the form of formed vortex has the predominate effectiveness on the
process of un-baffled surface aerator. (Rao et al., 2009) built a general correlation for
specific configuration with emphasizing on the vortex forms. In order to correlate the
relevant parameters of the process with depending on previous work of (McWhirter et
al., 1995) to describe the formation and performance of vortex especially, the role of
critical speed of the impeller and other effective parameters in the system.
(Rao et al., 2009) performed experimental runs to deduce the mass transfer
coefficient related with vortex formation and they found that the mass transfer
coefficient at standard condition increase sharply near or above the critical speed. The
general modified correlation for scale-up purposes was:
hv / D = 43.2 Fr(0.1 Ga0.18)(h Cu W/L)0.16
(1.22)
Where hv is the depth of the vortex; Cu is the turbine blade clearance (in this case it
represents the distance between the horizontal bottom of the tank and the top of the
blades); L represents blades length; Ga = Re2/Fr. The model was developed within the
ranges hv / D (0.003-0.2) and Fr (0.005-0.02).
Figure (1.16): Bubble distribution of surface aerator system of rotation speed N=110
rpm, liquid level h=0.66 m. (Lee et al., 2001).
(Deshmukh and Joshi, 2006) pointed out that the gas holdup plays an important role
on the velocity field in the impeller region for the surface aeration system, which in
turn is affected by the impeller design, submergence and rotational speed of the
impeller.
These observations were noticed by (Sun et al., 2006). They found that the air bubbles
holdup profile was non-uniform in the surface aerated agitated tank, the bubbles
holdup was high in two regions, first close the liquid surface and second in the
impeller region. The weakest gas holdup region was under the impeller region. The
air bubbles holdup was measured at different operation conditions; the impeller was
Rushton disc turbine with surface baffle used.
- 30 -
Figure (1.17): The relation between the rotation speed and mixing time for surface
aeration dual impeller system, (Kang et al., 2001).
(Guillard and Trgrdh, 2003) found out the general approach is difficultly applicable
in actual sized tanks and they deduced that the mixing time under aerated condition is
longer than non-aerated and depends on the tracer injection position also. They used
tracers by pouring pulses of concentrated acid; the injection was in three different
positions in the tank then the concentration gradient determined by PH electrode.
(1.23)
Where, S is the upper impeller submergence, and Fg, is the modified aeration number,
which equals G/(Nd3nim). The model was developed within the ranges Sp/D =2.12, Re
(2.0103-14103), D/Tv= 0.16 and h/Tv=3.5.
For micro- mixing in surface aeration process, (Rao and Kumar, 2009) derived a
model for mixing time by depending on the basis of theoretical power per unit volume
to simulate the process time with X = Fr4/3Re1/3,the ranges for the model are same in
equations 1.14 and 1.15.
(1.24)
- 31 -
Figure (1.18): The relation between the Reynolds number and power number for three
different liquid condition ns, clean water and two types of activated sludge mixed
liquid (Takase et al., 1982).
(1.25)
(1.26)
( )
Where de is the geometric equivalent dimension (
), c is dissolved
(1.27)
(1.28)
K2 is the equation constant. The models 1.25 to 1.28 were developed within the ranges
D / Tv =0.125, h/D= 3.0, Np (1.0-3.0) and Fr (0.05-4.0).
(1.29)
(Cancino et al., 2004) deduced that it is required to increase the contact time of water
droplets in the air to enhance the mass transfer operation and this can be done in the
way of increasing the height of water spray that projected or by increasing the
dispersion of the water droplets (decrease the size of the droplets). (Cancino et al.,
2004) demonstrated that to achieve these objectives it is required to increase also the
projected water flow rate to power consumption ratio (Q/P). (Cancino, 2004a)
- 34 -
(2.31)
Where VOCER and VOCERs represent VOC emission rates; Cvoc is the dissolved
VOC concentration in the bulk liquid; Hc and Hcs are Henrys law constant of VOC at
water and air wet-bulb temperatures respectively; and CG represents VOC
concentration in the air.
Water and wastewater treatment plants that implement surface aeration are well
known as sources for noise pollution, so it is recommended to build these plants away
from urban zones.
1.3. Conclusions
In this chapter a brief presentation was made for the present day aeration systems for
the water and wastewater treatment with introducing the types and efficiencies for
each system. The surface aeration for water treatment was described somehow in
detail.
It is found that; many modification essays are made on the surface aerators to improve
the performance. There are many trials have performed to enhance the operation
efficiency by either of increasing the mass transfer rate or by reducing the power
- 35 -
- 36 -
- 37 -
- 49 -
- 50 -
Chapter Two
Experimental Setup and Calculation Methods
2.1. Introduction
The experimental design and setup for surface aeration runs depend on many
important influencing factors, that generally can be classified as the geometric
configuration of the system, the operation conditions, the theoretical assumption that
describes the occurred process during the experimentation such as the models of the
transferred oxygen mass in water, the flow pattern of the air-water flow inside the
tank and at water surface, the physical and chemical properties for both water and air
used. The economic factor is an important factor that should be taken into account.
The experimental work was performed to interpret the important aims: the energy or
power consumed by the surface aeration system, the oxygen mass transfer coefficient
in water bulk kla and in the water spray (droplets) klad, the impeller configurations and
operational influences on the oxygen mass transfer and power consumption. Where,
the experimental runs can be classified into two main categories, the mass transfer and
hydrodynamics.
- 15 -
Vessel
Turbine
Draft tube
Propeller
Cone
Baffles
Draft tube
baffles
Impellers
spacing
Water
height
Diameter
(m)
0.8 (Tv)
0.19 (D)
0.15 (df)
0.12 (dpr)
0.15
Clearance
Tv/3.13
Tv/10
Tv/5.13
Height
(m)
0.6
Width
(m)
Blade width
(m)
No .of
blades
0.024
12
0.018
Numbers
0.1
0.06
0.6
Tv/10
0.06
df /8
0.076
0.28
- 15 -
(2.1)
Where, To, and Toe, are the measured torques in filled and empty vessel respectively
in (Nm). The torque meter used with a torque capture transducer all are made by
(HBM). A motor (LEROY SOMER), (LS90SL, rated power 1.1 kW, max. rated speed
1420 rpm) utilized to achieve the desired rotation speed.
(4)
Figure (2.4): The Laser Doppler Velocimetry (LDV) testing apparatus, (1) Laser
source. (2) Traverse system, (3) Tested tank, (4) Flow velocity analyser.
(2.2)
The quality of the LDV results depends on the size of the ellipsoidal measuring
volume formed at the beam intersection point and the number of the interference
fringes in this volume. The size is refined as much as possible in order to avoid the
presence of too many particles with different velocities and thus a random out of
phase signal. At the same time the number of interference fringes must kept
maximized in order to maximize the precision of the frequency of the signal (Aubin,
2001).
The dimensions of the measuring volume can be determined as (See Fig. 2.6):
- 15 -
(2.3)
2bv = dp / cos(/2)
(2.4)
2cv = dp
(2.5)
The two velocity components vr and vz, are measured by changing the position of the
laser beams in the vessel. To measure the velocity component, the laser beams was
oriented such that they are in the same plane as the component and the bisector of the
angle formed by the two beams is perpendicular to it. Aligning the beams in a
horizontal plane allows the radial velocity component, vr to be measured and when in
a vertical plane the axial velocity component, vz can be determined. So, the measured
velocity components by the LDV laser beams are the axial and radial components
only. The other component of the flow velocity that the tangential is not measured,
especially the measurement concerned area applied by the LDV is the propeller
vicinity and the draft tube inlets, where the majority of occurred velocities are the
axial and the radial, where the laser beams should move in a direction perpendicular
to the optical axis (See Fig. 2.7).
IV. Signal Post-Processing
A Dantec burst spectrum analyzer (BSA) for time-resolved measurements is used as
the LDV signal processing, where the signal from the photomultiplier that contained
in the apparatus was filtered to remove the undesired low and high frequencies then
amplified by a variable gain. The number of samples required to achieve statistically
independent results depends on the flow field. It was found that at least 300 samples
- 15 -
- 15 -
- 15 -
Figure (2.9): The Particle Image Velocimetry (PIV) testing apparatus, (1) Laser
source, (2) Recording camera, (3) Tested tank.
The light is not blocked by the small particles but it spread in all directions, so
massive multi-scattering occurs when sufficient number of tracer particles exists
inside the light sheet, where the light scattered by more than one particle is imaged so
the recorded light by the lens is not only because of the direct illumination but also
because of the fraction of the light (See Figure 2.10).
Figure (2.10): Light scattering by the (10 m) glass particles in the water (Raffel et
al., 2007).
V. Laser Source
The laser is used in the PIV process is related to their ability to emit monochromatic
light with high energy density, where it easily can be bundled into light sheet for
illuminating and recording the tracer particles.
- 56 -
(2.7)
- 55 -
(2.8)
Where M: is the magnification factor. The image intensity field for the time of the
second exposure:
N
(2.9)
Where, Vo: is the volume of the interrogation during second exposure. So the crosscorrelation function of the two interrogation areas will be:
RII (s, , Ddsp)= 1/aI Vo (Xi) Vo(Xj+ Ddsp) aI (X-Xi) (X- Xj + S - ddsp)dX
(2.10)
The comparison between the PIV and LDV techniques showed good agreement
between them for the quantitative comparison of the radial and axial velocities (Myers
et al., 1997).
- 55 -
Common types of propellers have the pumping numbers around 0.5. For ordinary flat
blade turbines and blade turbines are 1.3 and 0.87 respectively (Roustan, 2005).
The pumping number, NQp, for the propeller is obtained in depending upon the radial
and axial flow balance; the total inflow should equal the total outflow. The discharge
flow, Qp, is generated from propeller blades rotation is calculated from the measured
values of axial and radial velocities at known distance from propeller blades within a
determined volume calculated by equations 2.12a, 2.12b (See Figure 2.11). At
propeller vicinity the axial flow, Qpz, is calculated by measuring the axial velocities at
a distances of 1 mm above and below from the propeller edges, where these distances
are usually depend on the rotation position of the propeller. To calculate the radial
flow, Qpr the outer radial borders was taken at 1 mm from propeller edges side. The
formula of pumping rate calculation is depending on the velocity component occurred
within the controlled area as it explained with following equations (Sardeing et al.,
2003):
- 55 -
(2.12b)
dr
VZ+
Vr
dz
-
VZ
Figure (2.11): The pumping number and pumping flowrate measuring volume.
- 51 -
(2.13)
The maximum occurred axial flowrate in the vessel is assumed to be equal to the
maximum radial flowrate (Jaworski et al., 1996). The circulation flowrate is
normalized by dividing by ND3 and known as circulation number.
(2.14)
Since the circulation number presents the overall flowrate entrained by the impeller in
the vessel, therefore the NQc is commonly greater than NQp (Jaworski et al., 1996).
(2.15)
The 3-D flow inside the stirred vessel can represented by the corresponding 2-D
velocities for each volume grid that related to the vessel dimensions (See Figure
2.12), where the volume weighted average 2-D velocities can be used to compare the
agitation indices. These velocities are determined in assuming that each composite
mean local 2-D velocity (vij) corresponds to a 3-D cell liquid volume and related to
the grid coordinates (Garcia-Cortes et al., 2006; Mavros and Baudou, 1997). The
volume weighted velocity for the vessel is calculated as:
the volume swept by the impellers and the volume occupied by the cone and baffle
edges, the liquid volume calculations details were carried out as determined in
(Mavros and Baudou, 1997).
- 55 -
The 2-D composite mean local velocities was calculated for each grid point in
depending on the measured mean velocities in the r-z plane performed by the PIV, as
following (Garcia-Cortes et al., 2006; Mavros and Baudou, 1997; Mavros et al.,
1997):
(2.17)
Figure (2.12): The 3-D liquid volume cell and the vessel volume grids, (Garcia-Cortes
et al., 2006).
Protective Cap
Cap Sleeve
the Membrane
the Cathode
the Anode
The optical probe (HACH, LDO type) and the dissolved oxygen meter (HACH, type
HQ 10 (OENODEN)) were used. With this probe, the recalibration is not necessary
and the operation principle depends on the fluorescence properties of the electrode.
The amount of dissolved oxygen is determined from zero concentration to the
saturation. The probe includes blue LED, red LED, filter, lens, photo detector and
oxygen preamble foil. The concentration of the dissolved oxygen depends on the
delay by the oxygen molecules held back on the gas permeable foil in the probe cover
cap for the emitted light by blue LED inside the probe. The returning signals are
detected by a photo detector. This shifting or delay between the returning red light
and the blue excitation is measured to determine the dissolved oxygen quantity in
water (see fig. 2.14). The average percentage error of probes readings was calculated
by least squares best fit between experimental and theoretical readings. Average
percentages error were (3.49 %) and (4.8 %) for polarographic and optical probes
respectively.
LED
source
photodetector
cover
cap
foil
- 55 -
- 56 -
Because the measured kla for seven different positions is close and the difference
noticed is very small, the assumption of well mixed tank can be verified for the used
probes accuracy and for applied experimental conditions (Philichi and Stenstrom,
1989).
- 55 -
Figure (2.16): (a) Experimental oxygen probes positions in the liquid volume,
(b) Experimental oxygen probes positions in the liquid volume, N=2.5 rps, (h/Tv) =
0.35.
- 55 -
polargraphic probe 1
optical probe 2
ln(Cs-Ct/Cs-Co)
2.5
2
y = 0.0024x + 0.1455
1.5
y = 0.0024x + 0.044
0.5
0
200
400
600
Time, (s)
800
1000
Figure (2.17): The repeatability of the mass transfer (kla) experimental results for both
oxygen probes, (N= 2.08 rps).
(2.18)
While the experimental DO values are determined directly during the experimental
measurements of the probes response time through the direct changing from the zero
DO level solution (prepared by passing the Nitrogen gas through the water), then into
- 55 -
Ct , Experimental
10
8
R = 0.9973
6
4
2
0
4
6
8
Ct , Theoretical
10
Ct , (mg/l), Experimental
10
8
R = 0.9979
4
2
0
10
- 55 -
(2.20a)
By rearranging equation (2.20a):
(2.20b)
On the other hand, the effect of probe response time on the dissolved oxygen
concentration measurement for aeration processes is identified according to the
following equation (Sardeing et al., 2005):
(2.21)
To evaluate the probe response time effect on the experimental results and the
calculated mass transfer coefficient kla, a comparison was made for same
experimental and calculated conditions. The dissolved oxygen concentrations in both
experimental and theoretical that calculated by equation (2.21) for the polarographic
and optical probes respectively are shown in the figure (2.20) and (2.21). The oxygen
mass transfer relation with time with considering the probe response time was implied
in the applied mass transfer equation that represents the experimental runs to take in
account the experimental errors that are produced due to that. The oxygen mass
transfer coefficient is calculated theoretically with applying equation 2.21, where the
effect of time lag constant for the oxygen probes are implied in our calculations of
oxygen mass transfer in water bulk zone kla.
- 51 -
8
R = 0.978
6
4
2
0
2
4
6
8
Ct, (mg/l), (Experimental)
10
Figure (2.20): The probe response time relationship verification for both theoretical
and experimental DO values for the optical probe.
Ct , (mg/l) , (Theoretical)
10
8
R = 0.9902
6
4
2
0
2
4
6
8
Ct , (mg/l), (Experimental)
10
Figure (2.21): The probe response time relationship verification for both theoretical
and experimental DO values for polarographic probe.
I.8. Measurement Procedure for the Bulk Zone Oxygen Mass Transfer
Coefficient
At the beginning of the experimental runs, the water volume was de-oxygenated by
the bubbling of Nitrogen gas through it for at least 15 minutes until DO concentration
becomes below (0.5 mg/l). The DO measurements were conducted directly by the two
probes that situated in their positions as shown in the Figure (2.16), where the
dissolved oxygen was measured each 30 seconds till reaching the saturated level,
where run was terminated. The operation parameters and geometric configuration
- 55 -
(2.22)
The volumetric coefficient of oxygen mass transfer is derived by depending on the all
the resistances in the water to oxygen transfer. The volumetric coefficient of oxygen
transfer, (kla)T is measured at operation temperature T, and corrected to standard
temperature here usually T=10 oC or 20 oC by the equation 2.22 (Roustan, 2003).
I.10. Oxygen Transfer Rate for the Bulk Mass Transfer Zone (OTRb)
The oxygen transfer rate (OTR), which is the transferred mass of oxygen rate at
experimental condition in the total volume of the filled vessel with water.
(2.23)
To have complete view for oxygen mass transfer for all surface aeration process and
to realize the contribution of droplets spray mass transfer zone in the oxygen mass
transfer operation, oxygen mass transfer rate in other side of water bulk re-aeration
mass transfer zone OTRb is needed to be calculated. OTRb was determined by
applying the following relation, which is modified to be agreeable for implemented
system depending upon mass transfer expression derived by (Chern and Yang, 2004):
(2.24)
I.11. Standard Oxygen Transfer Rate for the Bulk Mass Transfer Zone (SOTRb)
The standard oxygen transfer rate (SOTR)b is the transferred mass of oxygen rate at
standard condition in the total volume of tank filled with water. The SOTRb is the
standard criterions that employed to determine the surface aeration performance
- 55 -
(2.25)
I.12. Standard Aeration Efficiency for the Bulk Mass Transfer Zone (SAEb)
The standard aeration efficiency (SAE)b sometimes it called the overall transfer
efficiency which represents the mass of oxygen transferred to the water per the
consumed power (Cancino et al., 2004; Sardeing et al., 2005)
- 55 -
II.3. Determination Model for the Spray Zone Oxygen Mass Transfer Coefficient
(klad)
The representative model for the spray zone that explains the oxygen mass transfer
during the travel of the droplets in the atmospheric air from turbine blades tips until
their impingement points at the water surface. The droplets trajectory are determined
by the following unsteady state oxygen mass balance equation (McWhirter et al.,
1995):
(2.27)
In equation (2.27) several assumptions are made to facilitate the explanation of the
oxygen transfer process for the intended zone. Each water droplet is assumed to be
completely mixed; the dissolved oxygen distribution is uniform within the droplets at
any time when they traverse the atmospheric air until their impingement points at
water surface with no back mixing might take place during this flight. The dissolved
oxygen concentrations at the limits of the droplet zone are related with other
presented zones, where at the zone inlet the concentration at any time is the same of
the well mixed liquid bulk for the dissolved oxygen concentration level inside the
- 55 -
(2.28)
Equation (2.29) shows that the contacting efficiency can be varied from its minimum
value (0.0) (no mass transfer occurred) to its maximum value of (1.0), which is an
ideal condition that can be reached at infinite time when equilibrium is achieved
between the water droplets and the atmospheric air (McWhirter et al., 1995). The
equilibrium level for the dissolved oxygen at the droplets outer surface depends on the
- 56 -
In the equation (2.30), the definition of Murphree contacting efficiency relates the
boundary levels for the actual and theoretical oxygen concentration gradients that take
place from water droplets propelling position to the impingement with water surface
location. In order to facilitate the calculation of the oxygen mass transfer results, the
previous equation (2.30) is rearranged in more appropriate correlation form (Baylar
and Bagatur, 2000a; McWhirter et al., 1995):
(2.31)
From equation (2.31) it can be estimated that the plot between the resulted oxygen
concentrations at the droplets impingement point with water surface at the
corresponding recorded water bulk dissolved oxygen gives a straight line with a slope
of (1-Esp), and the intercept with y-axis will be (Esp Cds).
Same as the surface aeration the similar efficiency calculation technique is applied in
other aeration processes like the aeration by weirs (Baylar et al., 2006; Kim and
Walters, 2001), where it is required to find a uniform basis for comparison purposes
between the obtained results, which can be achieved by normalizing the calculated
efficiencies to standard condition.
(2.32)
)
- 55 -
Where the (Esp)20 is the spray aeration efficiency at 20 oC, and Esp is spray transfer
aeration efficiency at the water bulk temperature, f is the exponent determined by the
following equation:
(2.33)
II.5. Oxygen Transfer Rate in the Spray Mass Transfer Zone (OTRsp)
It is necessary to take in account the relative humidity of the air when this approach of
OTRsp is applied if a remarkable change is occurred with the relative humidity during
the experimentation because the saturation level that used at the air wet-bulb
temperature is also depends on the air relative humidity.
Another important factor in the oxygen transfer in the water droplets spray zone is the
oxygen transfer rate OTRsp, which beside the standard aeration efficiency depends on
the water volumetric flow rate discharged by the turbine blades (Chern and Yang,
2004; McWhirter et al., 1995):
(2.34)
Figure (2.22): The droplets zone mass transfer coefficient (klad) measurement.
(2.35)
Equation (2.35) is derived principally from Newtons second law of motion, where
the initial vertical velocity of the water droplets has a changeable value during
droplets throwing operation (as it will be seen in next section). The droplets flight
time is affected by its trajectory, which in turn is affected by many parameters beside
rotation speed that act as lift force on droplets, like gravity, buoyancy and drag force
(Matsuura et al., 2003).
- 55 -
Figure (2.23): Schematic diagram illustrates the surface aeration water droplets spray
from turbine blades till the impingement point with some important relevant
dimensions.
(2.36)
- 55 -
(2.37)
The used straight tip blade turbine propels water droplets toward radial and tangential
directions in a horizontal plane, is parallel to the water surface, as it has been
described in Figure 2.24. To calculate these velocities for the given applied geometry
the following relationships (See equations 2.38 and 2.39) are employed (McWhirter et
al., 1995) with equation (2.35) for droplets flight time:
(2.38)
(2.39)
Figure (2.24): Droplets radial and tangential velocities propelled in horizontal plane
by turbine blades.
- 51 -
(2.40)
(2.41)
(2.42)
When the droplets spray velocity is known, the volumetric flow rate of the droplets
that are propelled at turbine blades edges can be determined by applying the overall
conservation of energy to surface aeration turbine and the overall power consumed
through blades to project these droplets into atmospheric air. It is assumed the power
consumption is delivered to the free water surface as a kinetic energy, which
contributes to elevate water velocity up to discharge velocity leaving the blades as
droplets thrown in the air with power conversion is 100% effectively achieved (Baylar
et al., 2001). The energy balance relation gives an appropriate relation between the
spray volumetric flow rate and the power consumption (Baylar et al., 2001;
McWhirter et al., 1995). (Baylar et al., 2001) used in their approach the pumping
power relation in weir aeration process by water droplets creation. The relation was
applied within the boundaries between jet point and plunging point of droplets. In
surface aeration process the previous relation can be considered as a convenient
- 55 -
Where, the constant 1.38910-7 depends on the applied geometry of the surface aerated
equipment. For eliminating errors that may occur when whole system is used during
the energy balance calculation, the power consumption was measured only with
presence of the turbine to assure all the energy was converted to propel water droplets
into atmospheric air (i. e. the power consumption by the propeller with draft tube do
not have any relation with droplets generation outside the water surface but they
consume power only for water bulk mixing purposes inside the tank).
2.8. Conclusions
The procedures for measurements and calculation methods for the oxygen mass
transfer coefficients in both water bulk and water droplets spray zones, power
consumption, mixing time, water spray discharge flowrate, water droplets flight time
and other related parameters are represented in this chapter. These methods are
applied in the experimentation to propose the results that are showed in the next
chapters. The measurements of the flow field and mean velocity profiles in all parts of
the tank that achieved by Laser Doppler Velocimetry LDV and Particle Image
Velocity PIV are also have been described in this chapter.
- 55 -
- 55 -
- 89 -
- 90 -
Chapter Three
The Oxygen Mass Transfer in the Aeration Mode
The surface aeration system has the capacity to operate in two different operation
modes;(i) Firstly, The aeration mode, where two phases (water and air bubble)
condition occurs, this mode is always exist when the turbine is employed. Mode,
where the single phase condition occurs, the mixing is accomplished by the RTP
propeller (no oxygen mass transfer takes place in the mixing mode).
The oxygen mass transfer in the surface aeration system is divided into two main
mass transfer operation zones of water bulk and spray zones according to the way of
contact between the water and the air (See Figure 2.14). The two occurring oxygen
transfer zones have two phases; the dispersed phase of the water droplets in the
continuous phase of atmospheric air. In water bulk mass transfer zone these two
phases are reversed, where the continuous phase is the water bulk and the dispersed
phase is the air bubbles. For these two zones there will be two different oxygen mass
transfer coefficients (McWhirter and Hutter, 1989). The water spray zone represents
the oxygen mass transfer that takes between the atmospheric air and water droplets in
the spray. While the water bulk zone represents the oxygen mass transfer that takes
place inside aeration tank between air bubbles and water bulk. It is so hard to identify
clearly the limits between the two zones. Theoretically the water bulk zone doesnt
begin directly when the water droplets impinge water surface where air bubbles are
created and entrained inside water bulk. Due to this impingement and the direct
contact with the atmospheric air a separate mass transfer zone occurs, where a reaeration takes place between the water surface and the atmospheric air in a manner
close to that happened in the water spray zone. This water surface aeration is actually
not a part of water bulk zone. But it is highly combined with water bulk transfer zone
due to flow profile and mixing effects that takes place in the tank zone, so practically
the effect of surface re-aeration is included in this water bulk zone.
Figure (3.1): Different tested geometric configurations, (a) Whole system, (b) Turbine
alone, (c) Turbine + Propeller, (d) Turbine + Draft tube.
- 19 -
7.5
6.5
5.5
4.5
3.5
2.5
1.5
0.5
1.5
2.5
3.5
N, (s-1)
Figure (3.2): Impellers speed effect on the oxygen transfer coefficient with four
different geometrical configurations; (Whole System, without draft tube, turbine alone
and without propeller), D/Tv=0.24, C/Tv=0.31, h/D=1.47, Cpr/Tv=0.2, S/W=1.
10
Saturated DO line
9
8
7
6
1
0
200
400
600
800
Time, (s)
1000
1200
1400
1600
Figure (3.3): The DO profile for different rotation speeds, h/Tv= 0.35, C/Tv=0.313,
Cpr/Tv=0.2, temperature =15 oC.
For the different geometrical configurations tested, the oxygen transfer coefficient kla
always shows a dependence on the impellers rotational speed, as it can be seen in the
- 19 -
(a)
- 19 -
N=2.08
rps
S/W=1.42
S/W=1.00
S/W=0.58
kla*10-3 , (s-1)
1
0
1.30
6
5
kla*10-3 , (s-1)
1 (S/W), (-) 2
1.40
1.50
(h/D), (-)
1.60
1.70
1.5
2.5
N, (s-1)
(b)
(c)
Figure (3.4): (a) The limits of the tested turbine blades submergence, (b) Mass
transfer Coefficient kla relation with turbine blades submergence and water height in
the liquid bulk for three levels of rotational speed, (b) Mass transfer Coefficient kla
relation with the rotation speed for three levels of turbine submergence. D/Tv= 0.24,
C/Tv = 0.31, Cpr/Tv=0.2.
Figure 3.4 shows that higher S/W ratio leads to higher kla until the ratio S/W=1.75
after this point the kla values are decreased. This relation can be explained that with
higher (S/W) more water droplets are thrown into the air, consequently larger
interfacial area achieved between air and water and also the flight trip of these
droplets goes farther, so as a result the kla mount up, but at the lower impellers speed
1.67 rps, the kla stays relatively constant despite elevation the S/W ratio.
At higher impellers speeds, further increasing with S/W higher than 1.58 the k la goes
down, as illustrated in (Figure 3.4). The dropping down of kla with further increasing
of water level ratio is happened, when it is noticed with higher water levels the spray
begins to be deformed; the interfacial contact area is constricted. For the water levels
lower that 100% submergence (h/D = 1.47) of the aerator blades, the values of kla are
moderate. This relation is depending on the desired shape of droplets spray as well as
preserving interfacial area between air and water.
Figure 3.4, displays the influence of water level on the kla, as the higher (h/D) ratio
leads to higher kla produced, This relation can be explained that with higher (h/D),
more efficient fluid circulation loop is accomplished in the tank, so as a result the kla
is increased, but at the lower impellers speed 1.67 rps, the kla stays relatively constant
despite elevation the h/D ratio.
- 19 -
1.2
1.7
2.2
2.7
N,
3.2
3.7
(s-1)
Figure (3.5): Spacing effect on the oxygen transfer coefficient, D/Tv= 0.24,
C/Tv = 0.31, Sp (Reference) =7.6 mm.
V. Power Consumption Measurements
In the aeration mode the power consumption was studied for various turbine blades
submergences (S/W), for the spacing variation between two impellers and for the
different geometrical configuration modifications (Figs. 3.6, 3.7 and 3.8).
It is observed that the power consumed is changing due to all these modifications; the
power consumed is highly relevant to each replacement or modification performed in
different manners. The difference in power consumption was very slight for the
spacing alteration and for the draft tube elimination, while it was remarkably affected
by the increasing of turbine blade submergence levels.
Figures 3.6 a and b illustrate that the power consumption in aeration mode is
dependent on both the Reynolds number (the diameter of the turbine D, was taken for
calculating Re) and water levels or turbine submergence. The power number
decreases with increasing the Reynolds number, as higher power is consumed with
higher impellers rotation speed (See Figure 3.6b). The power number increases when
the turbine blade submergence increases, wherein higher power consumed with
turbine blade submergence and higher water level. In the case of the water level is
higher than upper tip of turbine blades h/D =1.53 with blade submergence ratio
S/W=1.42, the turbine blades need more power to propel water droplets through the
air. The resulted power num er s for the three tur ine su m ergence levels are
- 19 -
6
(S/W)=1.42
(S/W)=1.0
(S/W)=0.58
Np, (-)
0
20000
60000
100000
140000
180000
220000
Re, (-)
(a)
40.0
35.0
S/W=1.42
S/W= 1.00
30.0
S/W= 0.58
P, (watt)
25.0
20.0
15.0
10.0
5.0
0.0
1.0
2.0
3.0
N, (s-1)
4.0
5.0
6.0
(b)
Figure (3.6): (a) The relationship between power number and Reynolds No., for
different turbine blades submergence levels C/Tv = 0.31, (b) The relation between
power consumption and rotation speed, C/Tv = 0.31.
- 11 -
5
Position A (Spacing=7.6 mm)
Position B (Spacing =8.3 mm)
Np, (-)
1
40000
60000
80000
100000
120000
140000
Re, (-)
Figure (3.7): The spacing between two impellers effect on the power consumption.
Figure 3.8, shows that the power number Np relation with Re for geometrical
configuration effect on whole system, without draft tube and propeller and without
draft tube configurations, where the relation appears similar behavior with Figure 3.6.
It is important to explain that power consumption is less relevant to geometry
modification than the water level h elevation.
Figure 3.8 illustrates the effect of the draft tube and the RTP propeller on the power
number for the condition, where for all different configurations the water level is just
at the upper edge of turbine blades. It points to the power consumption was slightly
changed from whole system configuration for Reynolds number range (6104 to
18104). The power number profile was very similar for the draft tube elimination and
the turbine alone aeration mode. This refers that the presence of the draft tube and the
propeller doesnt contri ute to an increase of the power consumption.
Comparing the power consumption relation with mixing mode (down-pumping), the
power number in that case was relatively constant and it was around 0.6. In mixing
mode the RTP propeller operates only in liquid phase with low gas presence, therefore
the power number has close values for different Reynolds numbers, is for that power
number can be referred by one value.
- 911 -
Turbine+Propeller(Aeration)
Turbine alone(Aeration)
3.5
Propeller+DT(Mixing)
Np, (-)
3
2.5
2
1.5
1
0.5
0
50000
100000
150000
200000
Re, (-)
VI. Standard Aeration Efficiency (SAEb) and Standard Oxygen Transfer Rate
(SOTRb) for the Water Bulk Zone
The standard oxygen efficiency of the water bulk mass transfer zone SAEb refers to
the relationship between the power consumption with achieved oxygen transfer rate
and calculated by equation 2.26, the SOTRb was calculated by equation 2.25. It is
noticed that SAEb was increased according to the elevation of rotation speed (See
Figure 3.9). The SAEb was affected by the increasing of impellers rotation speed for
the different tested geometrical configurations. For the whole system configuration,
when N was elevated from 1.67 rps to 3.33 rps the SAEb value was approximately
doubled this behavior was repeated for all tested geometrical configuration. Starting
from N=2.08 the whole system configuration shows competitive behavior with
turbine alone configuration till higher N >= 2.92 rps, where it looks clearly the whole
system achieves highest SAEb this behavior can explained that the contribution of
draft tube in redirecting the flow up-ward. As illustrated in Figure 3.9 the SAEb
relation with N wasnt similar for all tested configurations. At lower N levels it seems
that the propeller and draft tube have no contribution in the SAEb. With higher N
levels, the turbine and propeller performance become closer to turbine alone
configuration. But in general it can be notice that the draft tube presence looks
somehow necessary to assist the propeller in its work of pushing more gas liquid
dispersion to turbine intake region. The other advantage of draft tube is its power
consumption was quite low. On the other hand it was found that the SAEb didnt
changed effectively and proportionally with the turbine blades submergence or water
level variations. From Table 3.1 it was noticed that higher turbine blades submergence
or higher water level leads to higher accomplished SAEb till reaching the turbine
- 919 -
3
Whole system
Turbine + Propeller
Turbine alone
2.5
1.5
1.5
2.5
N, (s-1)
3.5
Figure (3.9): The relationship between the standard aeration efficiency SAEb for three
aeration mode configurations (Whole System, turbine + propeller and turbine alone)
with the impeller rotation speed, D/Tv=0.24, C/Tv=0.31, h/D=1.47, S/W=1.
- 919 -
Table (3.1): The SAEb and SOTRb for different water level ratios (h/D) and turbine
blade width to blades submergence ratio (S/W), D/T=0.24, C/T=0.31, Cpr/T=0.2,
N=2.803 rps, (Whole system Configuration).
Power
(watt)
2.932
5.418
7.552
9.148
9.698
10.025
h/D
(-)
1.37
1.42
1.47
1.53
1.58
1.63
S/W
(-)
0.17
0.58
1.00
1.42
1.83
2.25
(kla)T
(s-1)
9.80E-04
2.09E-03
2.84E-03
3.02E-03
3.10E-03
2.79E-03
(kla)10 C
(s-1)
8.38E-04
1.72E-03
2.52E-03
2.64E-03
2.74E-03
2.51E-03
SOTRb
(kgO2.h-1)
0.00448
0.00953
0.01451
0.01572
0.01686
0.01596
SAEb
(kgO2.kW-1h-1)
1.5265
1.7591
1.9208
1.7185
1.7386
1.5923
The SOTRb and SAEb were decreased when turbine blades submergence increased
S/W from 1.83 to 2.25 (water level ratio h/D was increased from 1.58 to 1.63), where
the lowest SOTRb took place (0.01596 kgO2/h). As turbine blades submergence was
increased to S/W= 2.25 (water level ratio was increased to h/D=1.63) that will cause
deformation of water spray shape created by turbine blades, so variation in contact
interfacial area was available.
The standard oxygen transfer rate SOTRb was increased with increasing impellers
rotation speed or water level for the three configurations examined that can be
explained as same manner for kla results discussed before (see Tables from 3.1 to
3.4). From Table 3.2, it can be observed that SOTRb is higher when higher water
height level used. Table 3.4 elucidate the relation between the experimental results of
SOTRb and with SAEb; where high SOTRb was reached with high surface aeration
efficiency SAEb was achieved, but this increasing is applicable for water level ratio up
to 1.58(h = 0.3 m).
Table (3.2) demonstrates that standard oxygen transfer rate for water bulk mass
transfer zone SOTRb is acted same as the standard aeration efficiency SAEb. SOTRb
was increased according to increasing of impellers rotation speed N, when N was
increased that will leads to more air bubbles entrapped into well mixed water bulk, so
larger contact interfacial area between the water bulk and air bubbles is provided. In
other hand water droplets propelled by turbine blades take longer travel time before
reaching impingement at water surface, this meaning larger contact area between
droplets and atmospheric air; as a consequence higher dissolved oxygen concentration
was accomplished.
- 919 -
N
(rps)
(kla)T
(s-1)
(kla)10 C
(s-1)
SOTRb
(kgO2.h-1)
SAEb
(kgO2.kW-1h-1)
4.503
1.67
1.28E-03
1.1E-03
0.0063
1.4078
7.552
2.08
2.84E-03
2.52E-03
0.0145
1.9205
10.603
15.184
2.5
3.33
4.25E-03
7.89E-03
3.79E-03
7.01E-03
0.0218
0.0403
2.0570
2.6540
Tables (3.3, 3.4 and 3.5) are representing the SOTRb with different impellers rotation
speeds for three different configurations with the achieved oxygen mass transfer
coefficient at actual and standard conditions. From Figure 3.10 and from tables 3.3,
3.4 and 3.5 it can be noticed that the SOTRb is lowered when both the draft tube and
the propeller were removed or when the draft tube alone was removed from the
system. The explanation for this is similar to that mentioned earlier in this chapter for
the oxygen mass transfer. The same behavior was observed for spacing changing
between the two impellers and for the turbine alone implementation cases
respectively. The SOTRb was increased with impellers rotational speed (N)
increasing; this is logic as with higher impellers rotation speed higher oxygen mass
transfer coefficient is obtained.
Table (3.3): The SAEb and SOTRb for different impellers speed, D/Tv= 0.24, C/Tv=
0.31, Cpr/Tv= 0.2, h/D = 1.47 (Turbine and propeller configuration).
Power
(watt)
4.409
6.991
10.053
14.472
N
(rps)
1.67
2.08
2.5
3.33
(kla)T
(s-1)
1.23E-03
2.38E-03
3.77E-03
7.23E-03
(kla)10 C
(s-1)
1.10E-03
2.13E-03
3.35E-03
6.42E-03
- 919 -
SOTRb
(kgO2.h-1)
0.00631
0.01224
0.01926
0.03693
SAEb
(kgO2.kW-1h-1)
1.4318
1.7515
1.9154
2.5516
0.04
Turbine + Propeller
Turbine alone
SOTRb, ( kg O2 h-1)
0.035
0.03
0.025
0.02
0.015
0.01
0.005
0
1.5
2.5
N,
3.5
(s-1)
Figure (3.10): The relation between the standard oxygen transfer rate SOTRb for three
aeration mode configurations (Whole System, turbine + propeller and turbine alone)
with the impeller rotation speed, D/Tv=0.24, C/Tv=0.31, h/D=1.47, S/W=1.
Table (3.4): The SAEb and SOTRb for different impellers speed, D/Tv= 0.24, C/Tv=
0.31, h/D = 1.47, (Turbine alone configuration).
Power
(watt)
4.472
6.82
9.943
14.263
N
(rps)
1.67
2.08
2.5
3.33
(kla)T
(s-1)
1.28E-03
2.37E-03
4.07E-03
6.93E-03
(kla)10 C
(s-1)
1.21E-03
2.22E-03
3.75E-03
6.38E-03
SOTRb
(kgO2.h-1)
0.00697
0.01275
0.02154
0.03668
SAEb
(kgO2.kW-1h-1)
1.5588
1.8700
2.1664
2.5715
Table (3.5): The effect the spacing between impellers on SOTRb and SAEb, D/Tv=
0.24, C/Tv= 0.31, h/D = 1.47, Sp=8.3 mm.
Power
(watt)
4.932
7.448
10.352
15.038
N
(rps)
1.67
2.08
2.5
3.33
(kla)T
(s-1)
1.40E-03
2.50E-03
4.00E-03
7.20E-03
(kla)10 C
(s-1)
1.25E-03
2.22E-03
3.52E-03
6.33E-03
- 919 -
SOTRb
(kgO2.h-1)
0.00720
0.01277
0.02024
0.03643
SAEb
(kgO2.kW-1h-1)
1.46
1.71
1.96
2.42
I. Mass Transfer
From the previous studies that dealt with surface aeration dimensional analysis, it was
found that there are various trends to correlate different relevant variables. (Fuchs et
al., 1971) have correlated the mass transfer coefficient kla with the volumetric power
consumption (P/V) for tank volumes of (V > 200 liter). (Zlokarnik, 1979) proposed
two correlation models for both mass transfer and power consumption for various
types of surface aerator. The main distinction in the correlations made by (Zlokarnik,
1979) is the mass transfer parameter was represented as a dimensionless parameter
that called the sorption number Y, which represents both the transferred oxygen and
turbulence intensity (Y=klaV(/g2)1/3D3). (Zlokarnik, 1979) has related in his model
the oxygen mass transfer with the water surface flow characteristics. (Zeybek et al.,
1997) have applied the Box-Wilson method to optimize the values of mass transfer
correlation for both linear and non-linear models of kla in general application of
aerated agitated tanks. (Patil et al., 2004) performed various tests for different types of
surface aerators and they emphasized on the geometrical configuration that affecting
parameters such as liquid height and volume, tank diameter and impeller clearance.
The model was for optimum values of kla with respect to impeller submergence.
(Cancino, 2004; Chandrasekharan and Calderbank, 1981; Chisti and Jauregui-Haza,
2002; Deshmukh and Joshi, 2006; Kumar et al., 2010; Moulick et al., 2002; Rao,
1999) have attempted to correlate the oxygen transfer and / or the power consumption
for diverse affecting variables (for the aeration process), each one of them depending
on its implemented geometrical configuration and proposed relevant geometric,
material and process parameters.
- 919 -
(3.2)
Since the surface aeration process predominantly within the turbulent flow state, so
the Reynolds No. can be ignored, where for air-water submergence aerator the value
of Re is always ( > 104), so the Re, is considered irrelevant to the process objective
(Zlokarnik, 1979).
So equation (4.5) was reformed as:
(
(3.3)
Equation (3.3) was solved to determine the values of the constants by applying
multiple non-linear regressions. The following correlation was developed:
- 919 -
For surface aeration cases the constant (K1) which have the value 0.0322 in equation
3.4 it normally represents the effect of the system geometry, but it in our case in all
implemented geometrical configurations have the same effect because actually we
have kept the essential geometry for all the tested configurations, where the turbine
and tank geometries didnt changed during experimental runs.
From equation (3.4) it is noticeable that the water level has higher effect on the
dimensionless mass transfer parameter (kla/N) than the wave condition that
represented by Fr or the power consumption that represented by Np. This elevated
value for h/D parameters shows that the oxygen mass transfer is highly related with
the liquid volume and flow circulation condition in the tank
Thee mass transfer correlation (equation 3.4) is applicable for the three ranges, (P/V,
watt/m3) = (22 -100), (Fr) = (0.054 0.214), (h/D) = (1.37 - 1.58) with respect to the
applied geometry in the experimentation.
The comparison between the predicted values from applying above correlation
(equation 3.4) with the experimental results that obtained from 74 experimental runs
showed that the coefficient of determination of (0.985), which is accepted in the
frame of the error as shown in Figure 3.11a.
(3.5)
- 919 -
(3.6)
Since the flow regime is always kept in turbulent region (Re > 104), the Re is
considered irrelevant to the process objective same as mass transfer model. Equation
3.6 was solved to determine the values of the constants by applying multiple nonlinear regressions, the following correlations developed:
(
The power correlation applicable for the ranges (22< P/V, watt/m3 <100), (1.37< h/D
<1.63) and (0.054<Fr<0.214) according to the applied geometry in the
experimentation. The comparison of predicted values from power consumption
correlation with experimental results values obtained from 58 experimental runs
found that the coefficient of determination of (0.972), (See Fig .3.11b).
2.5
2
Predicted (kla/N)*10-3
R = 0.985
1.5
0.5
0.5
1
1.5
2
-3
Experimental (kla/N)*10
(a)
- 911 -
2.5
Predicted Np
R = 0.972
4
3
2
1
0
Experimental Np
(b)
Figure (3.11): (a) The comparison between the kla/N, predicted by the correlation
model (Eq. 3.4) with the experimentally resulted kla/N, (b) The Comparison between
the (Np) values predicted by the Eq. 3.7 with the experimentally resulted (Np) values.
- 999 -
DO (mg l-1)
8
C
1
CLt
Run1
Lt Run
Cdt Run1
CLt
Run2
C
2
Lt Run
Cds Run1
Cdt Run2
Cds Run2
200
400
600
800
Time (s)
1000
1200
1400
1600
Figure (3.12): Dissolved oxygen concentration profile with time of experiment in the
water bulk zone (Cdt), and the water spray zone (CLt) at N = 1.67 rps for two
experimental runs.
10
DO (mg l-1)
CLt
CLt Run1
Cdt Run1
CLt Run2
CLt
Run2
Cdt Run2
Cds Run1
Cds Run2
200
400
600
800
1000
1200
1400
Time (s)
Figure (3.13): Dissolved oxygen concentration profile with time of experiment in the
water bulk zone (Cdt), and the water spray zone (CLt) at N = 2.08 rps for two
experimental runs.
- 999 -
10
DO (mg l-1)
8
CLt
Run 1
C
Lt Run1
Cdt Run1
CLt
Run2
C
Lt Run2
Cdt Run2
Cds Run1
Cds Run2
100
200
300
400
500
Time (s)
600
700
800
900
Figure (3.14): Dissolved oxygen concentration profile with time of experiment in the
water bulk zone (Cdt), and the water spray zone (CLt) at N = 2.5 rps for two
experimental runs.
12
10
DO (mg l-1)
8
CLt
Run1
CLt Run1
Cdt Run1
Cds Run1
CLt
Run2
CLt Run1
Cdt Run2
2
0
Cds Run2
100
200
300
400
Time (s)
500
600
700
Figure (3.15): Dissolved oxygen concentration profile with time of experiment in the
water bulk zone (Cdt), and the water spray zone (CLt) at N = 3.33 rps for two
experimental runs.
- 999 -
10
Run 1
Run 1
Run2
y = 0.776x + 2.4913
R = 0.9884
Run 2
y = 0.791x + 2.3126
R = 0.9928
y = 0.772x + 2.2727
R = 0.996
y = 0.745x + 2.5914
R = 0.9909
Esp
Run 1= 0.224
Run 2= 0.209
2
0
Esp
2
0
10
Run 1= 0.228
Run 2= 0.255
10
12
Run 1
Run 2
Run 1
10
8
y = 0.6223x + 4.0444
R = 0.9912
y = 0.6689x + 3.0932
R = 0.9925
4
Esp
y = 0.518x + 5.0191
R = 0.9956
y = 0.529x + 4.4901
R = 0.9968
6
4
Esp
Run 1= 0.471
Run 2= 0.482
Run 1= 0.331
Run 2= 0.379
Run 2
10
10
10
Figure (3.16): The linear regression correlation of water spray zone aeration
efficiency by plotting Cdt versus CLt for various rotation speeds.
It is noticed from Figures 3.12 to 3.15, the impellers rotation speed gradually
increasing from 1.67 rps to 2.08, 2.5 and 3.33 rps causes the increasing of oxygen
mass transfer coefficient klad in the spray zone. The increased oxygen mass transfer
coefficient klad means that the DO inside the droplets reaches its higher concentration
level (nearest to air wet-bulb saturated oxygen condition) in shorter time, but in the
same time the flight time for the water droplets from turbine blades to the water
- 999 -
Tbulk
(oC)
(mg L-1)
(mm)
(mm)
Rm
Esp
(Esp)20
1.67
1.67
2.08
2.08
2.50
2.50
3.33
3.33
15.8
15.2
16.1
15.9
16.5
16.0
16.3
15.9
10.7
10.3
9.6
9.85
9.5
10.07
9.78
10.3
15
15
20
20
35
35
45
45
110
110
130
130
170
170
280
280
0.224
0.209
0.228
0.255
0.331
0.379
0.471
0.482
0.242
0.229
0.245
0.275
0.352
0.405
0.498
0.513
(rps)
Cds
Ym
tf
klad
0.11
0.11
0.13
0.13
0.17
0.17
0.19
0.19
2.305
2.131
2.000
2.264
2.365
2.802
3.351
3.462
(s)
(s-1)
4
3.5
klad, (s-1)
3
2.5
2
1.5
1
1.5
2.5
N,
3.5
(s-1)
Figure (3.17): The relation between the oxygen mass transfer coefficient klad and
impellers rotation speed and in the spray zone, S/W = 1.
- 999 -
(mm)
Ym
(mm)
Rm
(mm)
Rsp
Esp
1.67
1.67
2.08
2.08
2.50
2.50
3.33
3.33
15
15
20
20
35
35
45
45
85
85
115
115
175
175
270
270
180
180
210
210
270
270
365
365
0.224
0.209
0.228
0.255
0.331
0.379
0.471
0.482
tf
Vsp
(s)
(Watt)
(m s-1)
(L h-1)
0.111
0.111
0.128
0.128
0.169
0.169
0.19
0.19
4.472
4.472
6.820
6.820
9.943
9.943
14.263
14.263
1.71
1.71
1.76
1.76
1.80
1.80
2.14
2.14
11010
11010
15851
15851
22094
22094
22422
22422
From the first glance on the results illustrated in Table 3.7, it is obvious that the
droplets spray flow rate was increased when the turbine blades rotation speed was
increased, this agrees with (Takase et al., 1982), this relation can be easily explained
as the higher rotation speed leads to higher water droplets were projected through the
air by the turbine blades. The droplets velocities in the spray at the impingement
position at water surface is approximately considered equal along the water spray
from turbine blades tip to plunging point and the energy consumption within water
spray path was ignored since the droplets flight time is extremely short. With
elevating the turbine rotation speeds it was noticed that the spray velocity was
increased, this relation was achieved with higher rotation speed or in other word with
higher power consumed, so higher mechanical energy was converted to a kinetic
- 999 -
Esp
1.67
1.67
2.08
2.08
2.5
2.5
3.33
3.33
0.22
0.21
0.23
0.26
0.33
0.38
0.47
0.48
(rps)
tf
Vsp
OTRsp
(s)
(m s-1)
(L h-1)
(gO2 h-1)
0.111
0.111
0.128
0.128
0.169
0.169
0.19
0.19
1.71
1.71
1.76
1.76
1.80
1.80
2.14
2.14
11010
11010
15851
15851
22094
22094
22422
22422
25.649
23.011
33.972
39.005
68.012
80.973
99.050
107.642
- 999 -
120
100
80
60
40
20
0
1.5
N, (s-1)
2.5
3.5
Figure (3.18): The relation between the oxygen transfer rate OTRsp and impellers
rotation speed and, in the spray zone, S/W = 1.
I.5. Contribution Percentage of the Spray and Bulk Zones in the Overall Mass
Transfer Operation
To have complete view for the oxygen mass transfer for all surface aeration process
and to realize the contribution of spray mass transfer zone in the overall oxygen mass
transfer operation, the oxygen mass transfer rate in the water bulk mass transfer zone
OTRb is needed to be calculated. OTRb was determined by applying the equation
(2.24) (See Appendix II).
- 991 -
(s)
tf
(gO2 h-1)
OTRsp
(gO2 h-1)
OTRr
OTRsp%
OTRb%
1.67
1.67
2.08
2.08
2.50
2.50
3.33
3.33
0.11
0.11
0.13
0.13
0.17
0.17
0.19
0.19
25.649
23.011
33.972
39.005
68.012
80.973
99.050
107.642
6.753
6.810
10.183
10.214
20.621
19.935
35.307
35.541
79.16
77.16
76.94
79.24
76.73
80.02
73.72
75.18
20.84
22.84
23.06
20.76
23.27
19.98
26.28
24.82
Results from Table 3.9 and Figure 3.19 show that for different impellers rotation
speed, N levels (1.67-3.33 rps); the spray zone oxygen mass transfer OTRsp has a
contribution range between (80.02% - 73.72%) of the overall oxygen mass transferred
during the operation, which is higher than the contribution of bulk mass transfer zone
in the overall oxygen transfer rate. At higher N, the percentage of contribution for the
oxygen mass transfer rate in the bulk zone OTRb in the overall oxygen transfer rate
stays the same although the OTRb was increased with N elevation, where more
turbulent state of water surface was occurred by water droplets at the impingement
and plunging position due to more mechanical energy was converted to kinetic
energy. As a consequence of this condition more air bubbles were generated and
entrained into the water bulk, beside that more fine air bubbles were diffused into the
tank (Ozkan et al., 2006). So the range of the oxygen mass transfer rate of the bulk
zone is improved. The improvement of the OTRb assists the bulk zone to keep its
contribution somehow similar within the tested range of N. It is also important to
mention that both the OTRb and OTRsp were calculated at the same water bulk
temperature and ambient temperature for each run.
Percentage in the overall OTR, (%)
100%
80%
60%
OTRb%
40%
OTRsp%
20%
0%
1.67
2.08
2.5
3.33
N, (s-1)
Figure (3.19): The influence of the impellers rotational speed on the spray zone
oxygen transfer rate OTRsp and in the bulk zone oxygen transfer rate OTRsp
contributions the overall transfer rate OTR, S/W = 1.
- 991 -
- 999 -
h/D
(-)
1.37
1.37
1.42
1,42
1.47
1.47
1.53
1.53
1.58
1.58
Ym
Rm
Rsp
(mm)
(mm)
(mm)
23
23
26
26
33
33
35
35
38
38
55
55
85
85
170
170
180
180
187
187
150
155
180
185
265
265
275
275
282
282
Esp
(-)
0.182
0.185
0.217
0.223
0.353
0.340
0.362
0.360
0.392
0.384
tf
Vsp
(s)
(Watt)
(m s-1)
(L h-1)
0.093
0.093
0.109
0.109
0.164
0.164
0.169
0.169
0.176
0.176
4.214
4.214
6.754
6.754
9.943
9.943
11.624
11.624
12.872
12.872
1.74
1.74
1.79
1.79
1.81
1.81
1.82
1.82
1.80
1.80
10021
10021
15176
15176
21850
21850
25264
25264
28602
28602
- 999 -
10
8
S/W=0.17
S/W=0.58
S/W=1.00
S/W=1.42
S/W=1.83
4
2
0
200
400
600
800
Time, (s)
1000
1200
1400
Figure (3.20) Effect of the turbine blade submergence on spray mass transfer zone
dissolved oxygen concentration at rotation speed (N=2.5 rps) and turbine clearance
(C/Tv = 0.35).
From Figure 3.20 it was found also that elevating S/W ratio to 1.42 has no
enhancement on DO behavior, on the contrary higher S/W levels causes deformation
of water spray, in consequence the achieved DO is lowered. In the other hand at lower
S/W ratio (1.37 and 1.42) it seems that DO profiles with time are close together.
- 999 -
Run 2
y = 0.8177x + 2.8226
R = 0.9892
y = 0.8155x + 2.8716
R = 0.9804
Esp
Run 1= 0.182
Run 2= 0.185
2
0
Run 1
Run 2
10
Run 1
4
6
CLt (mg l-1)
y = 0.7823x + 2.2402
R = 0.9928
y = 0.7776x + 2.2324
R = 0.9936
Esp
Run1= 0.217
Run2= 0.223
10
5
CLt (mg l-1)
Run 1
Run 2
10
y = 0.6468x + 3.4357
R = 0.9938
12
Run 1
y = 0.6599x + 3.4707
R = 0.9915
4
Esp
Run1= 0.353
Run2= 0.34
2
0
4
6
CLt (mg l-1)
Run 2
y = 0.6377x + 3.9956
R = 0.979
y = 0.6398x + 3.8801
R = 0.9877
Esp
Run1= 0.389
Run2= 0.36
2
0
10
10
5
CLt (mg l-1)
10
y = 0.6076x + 3.8876
R = 0.9907
y = 0.6156x + 3.8405
R = 0.9946
Esp
2
0
Run 1= 0.41
Run 2= 0.42
0
10
Figure (3.21): The linear regression correlation of water spray zone aeration
efficiency by plotting (Cdt) versus (CLt) for various turbine blades submergence, N =
2.5 rps.
- 999 -
0.45
0.4
Esp, (-)
0.35
0.3
0.25
0.2
0.15
0.1
0.5
1
S/W, (-)
1.5
Figure (3.22): The relation between the spray zone aeration efficiency Esp and turbine
blade submergence ratio, N = 2.5 rps.
- 999 -
Table (3.11): The calculated water spray mass transfer zone (klad), droplets flight
time(tf) and some related measured parameters and turbine blades submergence and
water level, at constant rotation speed (N =2.5 rps).
S/W
(-)
0.17
0.17
0.58
0.58
1.00
1.00
1.42
1.42
1.83
1.83
h/D
(-)
1.38
1.38
1.42
1.42
1.47
1.47
1.53
1.53
1.58
1.58
Tbulk
Cds
Ym
Rm
(oC)
(mg L-1)
(mm)
(mm)
17.7
17.8
18.4
18.3
16.8
16.7
16.7
16.5
15.2
15.1
10.73
10.9
9.45
9.94
9.5
10.07
10.5
10.2
9.94
10.0
23
23
26
26
33
33
35
35
38
38
60
60
90
90
170
170
180
180
187
187
- 999 -
Esp
(-)
0.182
0.185
0.217
0.223
0.353
0.340
0.362
0.360
0.392
0.384
(Esp)20
(-)
0.190
0.193
0.224
0.230
0.373
0.360
0.383
0.382
0.424
0.417
tf
klad
0.093
0.093
0.109
0.109
0.164
0.164
0.169
0.169
0.176
0.176
2.160
2.199
2.244
2.449
2.655
2.534
2.659
2.641
2.827
2.785
(s)
(s-1)
0.5
1.5
2.5
2
1.35
1.4
1.45
1.5
1.55
1.6
klad , (s-1)
2.75
2.5
2.25
(h/D), (-)
Figure (3.23): The relation between the spray zone oxygen mass transfer coefficient,
klad with the turbine submergence and liquid level impellers rotation speed, N = 2.5
rps.
- 999 -
h/D
(-)
1.37
1.37
1.42
1.42
1.47
1.47
1.53
1.53
1.58
1.58
Esp
(-)
0.182
0.185
0.217
0.223
0.353
0.340
0.362
0.360
0.392
0.384
(s)
tf
(m s-1)
Vsp
(L h-1)
OTRsp
0.093
0.093
0.109
0.103
0.164
0.164
0.169
0.169
0.176
0.174
1.74
1.74
1.79
1.79
1.81
1.81
1.82
1.82
1.80
1.80
10021
10021
15176
15176
21850
21850
25264
25264
28602
28602
16.232
16.407
29.211
30.424
71.346
67.975
85.511
85.221
106.514
104.340
0.5
(g h-1)
(S/W), (-)
120
1.5
2.5
100
80
60
40
20
0
1.3
1.4
1.5
1.6
(h/D), (-)
Figure (3.24): The relation between spray zone oxygen transfer rate OTRsp with the
turbine blades submergence and liquid level, N = 2.5 rps.
- 999 -
100%
80%
60%
OTRb%
40%
OTRsp%
20%
0%
0.17
0.58
1.42
1.83
S/W, (-)
Figure (3.25): The influence of the turbine submergence on the spray zone oxygen
transfer rate, OTRsp and in the bulk zone oxygen transfer rate OTRsp contributions the
overall transfer rate OTR, N = 2.5 rps.
- 991 -
Table (3.13): Oxygen mass transfer rate and percentage of contribution for both water
spray and bulk (re-aeration) zones, (N= 2.5 rps).
S/W
(-)
0.17
0.17
0.58
0.58
1.00
1.00
1.42
1.42
1.83
1.83
h/D
(-)
1.37
1.37
1.42
1.42
1.47
1.47
1.53
1.53
1.58
1.58
tf
OTRsp
OTRb
OTRsp%
OTRb%
0.093
0.093
0.109
0.103
0.164
0.164
0.169
0.169
0.176
0.174
16.232
16.407
29.211
30.424
71.346
67.975
85.511
85.221
106.514
104.340
10.778
10.437
14.673
16.607
19.319
20.578
23.540
22.841
29.969
30.157
60.10
61.12
66.56
64.69
78.69
76.76
78.41
78.86
78.04
77.58
39.90
38.88
33.44
35.31
21.31
23.24
21.59
21.14
21.96
22.42
(s)
(g h-1)
(g h-1)
(%)
(%)
DO (mg l-1)
8
6
N=1.67 rps
4
N=2.08 rps
N=2.5 rps
N=3.33 rps
250
500
750
1000
1250
1500
1750
2000
Time (s)
Figure (3.26): The effect of the turbine rotation speed on the spray mass transfer zone
dissolved oxygen concentration at h/D=1.47and C/T=0.35, (Turbine alone
configuration).
The spray zone aeration efficiency Esp for the turbine alone configuration was
calculated by the linear regression for the Cdt versus CLt plot. Aeration efficiency Esp
was determined for the four different rotation levels by applying best fit with the least
squares method. The results are shown in Figure 3.27.
As in the whole system configuration, with higher N, higher water spray zone aeration
efficiency Esp is achieved. Figure 3.26 shows that when N was increased from 1.67
rps to 2.08 rps; the Esp was increased from 0.183 to 0.206. At moderate N=2.5 rps, the
aeration efficiency Esp was noticeably elevated to 0.292. At relatively elevated
N=3.33 rps, the Esp was found to be 0.356. From Figures 3.27 and 3.28 it is found that
the Esp is remarkably lowered when the propeller and the draft tube are removed.
The Esp values are relatively lower than the whole system configuration aeration
efficiencies for all tested N levels, for N= 1.67rps, the Esp was lowered from 0.21 to
0.183, while for N= 3.33rps, the Esp was lowered from 0.471 to 0.356. From these
comparisons it can be deduced that the presence of the propeller and draft tube has a
relative influence on the water spray efficiency Esp especially at higher rotation
speeds.
- 999 -
10
y = 0.817x + 1.5911
R = 0.998
y = 0.7949x + 1.8793
R = 0.9959
4
Esp= 0.206
Esp= 0.183
2
0
CLt (mg
l-1)
10
10
10
y = 0.7076x + 2.6796
R = 0.9919
6
4
y = 0.6449x + 3.1303
R = 0.9762
Esp = 0.292
2
0
10
6
4
Esp= 0.356
10
10
Figure (3.27): The linear regression of the spray zone aeration efficiency of the plot
(Cdt) versus (CLt) for various rotation speeds, (Turbine alone configuration).
The performances of both the draft tube and the propeller are to enable the system to
obtain an efficient mixing condition of the DO in the water bulk mass transfer zone
inside the aeration vessel, which has also a relative effect in the enhancing of the
aeration efficiency Esp in the spray mass transfer zone through increasing the
dissolved oxygen DO in the liquid bulk. The water bulk here is considered in the same
time as the intake of water spray droplets.
- 999 -
0.5
Whole system
Esp, (-)
0.4
0.3
0.2
0.1
0
1.5
2.5
N, (s-1)
3.5
Figure (3.28): The relation between the impellers speed N and the spray zone aeration
efficiency Esp for both turbine alone and whole system configurations, S/W = 1.0.
3.5
Whole system
klad, (s-1)
3
2.5
2
1.5
1
1.5
2.5
N, (s-1)
3.5
Figure (3.29): The relation between N and the oxygen mass transfer coefficient klad
for both turbine alone and whole system configurations, S/W = 1.0.
- 999 -
Tbulk
(oC)
(mg L-1)
(mm)
(mm)
1.67
2.08
2.5
3.33
19.6
20.8
21.0
19.7
8.85
8.74
9.28
9.01
15
20
35
45
110
130
170
280
(rps)
Cds
Ym
Rm
Esp
(-)
0.183
0.206
0.292
0.356
(Esp)20
(-)
0.184
0.203
0.287
0.358
tf
klad
0.11
0.13
0.17
0.19
1.837
1.774
2.031
2.316
(s)
(s-1)
Ym
Rm
Rsp
(rps)
(mm)
(mm)
(mm)
1.67
2.08
2.5
3.33
15
20
35
45
85
115
175
270
180
210
270
365
Esp
(-)
0.183
0.206
0.292
0.356
- 999 -
tf
Vsp
(s)
(Watt)
(m s-1)
(L h-1)
0.111
0.128
0.169
0.19
4.472
6.820
9.943
14.263
1.71
1.76
1.80
2.14
11010
15851
22094
22422
120
Turbine alone
Whole system
100
80
60
40
20
0
1.5
N, (s-1)
2.5
3.5
Figure (3.30): The relation between spray zone oxygen transfer rate with the
rotational speed for the whole system and turbine alone configurations, S/W = 1.0.
- 999 -
(rps)
1.67
2.08
2.5
3.33
Esp
(-)
0.183
0.206
0.292
0.356
(s)
tf
(m s-1)
Vsp
(L h-1)
OTRsp
0.111
0.128
0.169
0.19
1.71
1.76
1.80
2.14
11010
15851
22094
22422
16.622
26.253
55.353
67.131
(g h-1)
IV. Comparing the OTRsp for the Whole System and Turbine Alone
Configurations
To find out the effect of the RTP propeller and draft tube presence on oxygen mass
transfer rate in spray mass transfer zone OTRsp, a comparison was made between the
calculated OTRsp in both whole system and turbine alone configurations. It is found
that OTRsp in whole system configuration has moderately higher values than the
turbine alone configuration as illustrated in table 3.17 and Figure 3.31. The
improvement in OTRsp when propeller and draft tube placed under the turbine is
generally has an excess percentage around (30 %). The experimental results showed
that OTRsp is always greatly depending on the aeration efficiency for various
conditions tested in this work as seen in Table 3.16. The excess in the oxygen mass
transfer rate for the spray mass transfer zone, OTRsp, when the RTP propeller and the
draft tube are positioned in their places is due to the improvement of the mass transfer
rate in the tank as a result of flow pattern improvement which is highly related to the
oxygen mass transfer in the spray zone since the initial concentration of the DO in the
spray zone is in the same time the final concentration of DO in the bulk zone.
100%
80%
60%
OTRsp% (Whole system)
OTRsp% (Turbine alone)
40%
20%
0%
0.17
0.58
N, (s-1)
1.42
Figure (3.31): The influence of the tank internals (Propeller and draft tube) on the
spray zone oxygen transfer rate, OTRsp contribution the overall transfer rate OTR,
S/W = 1.0.
- 999 -
Table (3.17): The Comparison for the spray zone oxygen transfer rates OTRsp
between the whole system and turbine alone configurations, S/W =1.00, h/D=1.47.
N
OTRsp
OTRsp
(rps)
(L h-1)
Turbine alone
Whole system
1.67
2.08
2.5
3.33
11010
15851
22094
22422
(g h-1)
16.622
26.253
55.353
67.131
OTRsp
(g h-1)
24.330
36.488
74.493
103.346
31.68
28.10
25.69
35.04
(Whole system)
- 999 -
Geometrical Parameter:
Turbine diameter
Turbine blades width
D
W
Material Parameters:
Density
Viscosity
Process Parameters:
Gravitational constant
Turbine speed
Turbine blades submergence
g
N
S
The spray zone standard efficiency varies noticeably with the turbine rotation speed.
The standard spray efficiency depends also on the submergence level of the turbine
blades, where the turbine blades submergence S/W showed also a remarkable effect
on the standard spray zone efficiency.
As a primary essay to correlate the standard efficiency of the spray mass transfer zone
with the influencing parameters in the same way that applied for the water bulk zone
as represented in the equation 3.1.
(Esp)20 = f (N, , , D, g, S, W)
(3.8)
(3.9)
( )
Fr=0.54-0.215
(3.10)
S/W=0.17-1.83
Where Reynolds number Re was ignored because it is considered irrelevant to the
process objective (Esp )20 according to the multiple non-linear regressions, as the spray
flow is predominantly turbulent the value of Re is always > 104 (Zlokarnik, 1979).
The plot of the predicted values of spray zone standard efficiency calculated from
equation 3.9 with the experimental values showed a coefficient of determination of
(0.9095) as shown in Figure 3.32.
Equation 3.10 means the correspondences (Esp)20 (N)1.23 or (Esp)20 (S/W)0.355 are
achieved when the turbulent regime was prevailing (Re > 104) and within indicated
ranges of Fr and S/W.
- 999 -
R = 0.9
Predicted (Esp)20
0.5
0.4
0.3
0.2
0.1
0
0.1
0.2
0.3
0.4
Experimental (Esp)20
0.5
0.6
Figure (3.32): The Comparison between the (Esp)20, Predicted by the Correlation
model (Eq. 3.10) with the experimentally resulted (Esp)20.
From the earlier works (Baylar and Bagatur, 2006; Baylar et al., 2010; Nakasone,
1987; Tarshish et al., 2000; Wormleaton and Soufiani, 1998), they have found that the
form of the water spray standard efficiency correlation is more complicated than the
power low correlation applied in equation (3.10). For the same influencing parameters
as shown in equation (3.8), this was solved by applying multiple non-linear
regressions and least square method. The following model was derived:
( )
(3.11)
- 991 -
R = 0.9
Predicted (Esp)20
0.5
0.4
0.3
0.2
0.1
0
0.1
0.2
0.3
0.4
Experimental (Esp)20
0.5
0.6
Figure (3.33): The Comparison between the (Esp)20, Predicted by the Correlation
model (Eq. 3.11) with the experimentally resulted (Esp)20.
(Baylar and Bagatur, 2006; Nakasone, 1987; Wormleaton and Soufiani, 1998)
developed a model of weir type aeration, where they showed that spray discharge
flow rate and droplets height involves the effects on the standard spray zone
efficiency.
(Tarshish et al., 2000) have built a model for the achieved oxygen mass flow rate by
the surface aeration in the wastewater treatment tank, they showed the oxygen mass
flow parameter in the model is affected by power consumption and jet flow rates of
waste water and activated sludge. (Wormleaton and Tsang, 2000)studied the aeration
behavior in weirs with. For the spray zone, they derived the model of standard spray
aeration in terms of a spray Froude number and Reynolds number.
The effect of spray discharge flow rate and droplets height included in the previous
models are not considered separately here, but they are implicated in more general
dimensionless parameters.
In dependence on experimental results and earlier studies, it is found that standard
spray zone efficiency for the surface aeration, it is considered that the target quantity
of standard spray zone efficiency (Esp)20 for the given process conditions (the
discharge spray flow rate/spray radius, qsp and the maximum height of the spray Ym)
depends on the turbine blades submergence S, turbine blades width W, gravitational
constant g, and kinematic viscosity of the water . The following relevance list is:
[(Esp)20 ;qsp, Ym ;W, S; g, ]
(3.12)
- 991 -
(3.13)
(3.14)
(3.15)
The spray Reynolds number is calculated by equation (3.14), the spray Froude
number is calculated by the equation (3.15). Equation (3.13) was solved to determine
the values of the constants by applying multiple non-linear regressions and least
square method. The following correlation was developed:
(
( )
(3.16)
The mass transfer model (equation 3.16) is applicable within the ranges, (Resp*10-3) =
(22.1 - 46.2), (S/W) = (0.17 - 1.83), (Frsp) = (0.044 1.80) and (water spray maximum
height Ym/D) = (0.08 0.24).
Equation 3.16 shows that the determined powers for the three parameters are
somehow close and they are proportionally related the spray zone standard efficiency.
This means the effect of these parameters have close effects on the objective
parameter (Esp)20. For example the equation 3.16 refers that the standard spray zone
efficiency (Esp)20 increases by the factor 1.3 when the turbine blade submergence ratio
(S/W) was altered from 0.17 to 0.58.
The standard error of estimation for the equation is 0.0215, which indicates a 98%
probability that the equation 3.16 will predict (Esp)20 within + 0.01 of its true value.
The correlation coefficient for this equation is R2 = 0.9568. The measured (Esp)20
values are compared with those predicted by the equation 3.16 are given in the Figure
3.34. An acceptable agreement between the measured values and the values computed
from the predictive equation is obtained.
- 999 -
Predicted (Esp)20
0.5
0.4
0.3
0.2
0.1
0
0.1
0.2
0.3
0.4
0.5
0.6
Experimental (Esp)20
Figure (3.34): The Comparison between the (Esp)20, Predicted by the Correlation
model (Eq.3.16) with the experimentally resulted (Esp)20.
3.2.4. Conclusions
With the totality of experimental runs performed for the power consumption in water
bulk zone, it was found that the power consumed is related with rotation speed and
water height level (turbine blades submergence) variations more than other
parameters occurred in aeration mode.
For the given impellers rotational speed range, it was observed that the desirable
performance of the aeration mode and flow patterns in the bulk zone is achieved
starting with impellers rotation speed of 2.08 rps, the formation of air bubbles appears
in the vessel due to plunging water droplets with water surface, so the impellers
rotational speed 2.08 is considered as the critical speed to accomplish aeration
process.
At the given water level ratios and turbine blades submergence, the water bulk zone
oxygen mass transfer coefficient kla was increased with the rotation speed N. It was
observed that the kla was highly dependent on the impellers rotation speed, however
this relation is limited: the outer edge of the spray should not colloid the vessel wall
what can be observed with rotational speed N lower up to 3.33 rps.
An optimum water height level came across from the water bulk zone surface aeration
experimental testing and it depends on a satisfactory formation created by the aerator
turbine blades. Where, neither too low nor too high blade submergence is preferred.
- 999 -
- 999 -
- 999 -
- 145 -
- 146 -
Chapter Four
Hydrodynamics in the Single Phase Non-Aerated Agitated
Tank for Up and Down Pumping Directions Modes
The investigated surface aeration system consists in two agitation tools; the turbine
positioned at liquid surface which propels droplets into atmospheric air, thus forming
the water spray. The second tool is the mixing assembly (RTP propeller and draft
tube) placed below the turbine. The proposed function of the mixing assembly is to
improve mixing capability of the surface aeration system by enhancing air bubbles
distribution inside the vessel and redirecting the flow toward the upper turbine. The
mode of mixing assembly employment is either individual or conjugated with turbine,
where mixing assembly operates separately from the turbine when its needed to
achieve well-mixed condition inside the vessel, when sufficient dissolved oxygen is
accomplished but without satisfying distribution of dissolved oxygen throughout the
liquid bulk.
In this chapter the turbulent flow regime in liquid phase generated by the reversible
twisted pithed, axial 4-bladed and 45o pitched blade propeller (RTP) (Milton Roy
Mixing HPM204D) is investigated to determine velocity field for different areas of
propeller vicinity, inside the draft tube and in the vessel, whereas the upper turbine is
not combined. Thus there is only the liquid phase and very few gas bubbles are
entrapped inside the liquid bulk can be neglected. Comparison between flow patterns
resulting from up-pumping and down-pumping operation modes of RTP propeller is
made. The up-pumping condition is performed ordinarily with surface aeration
system; so the propeller turns in same manner than the upper turbine and pushes the
fluid upward. The down-flow is applied mainly when propeller and draft tube are used
separately aiming that mixing objectives are to be fulfilled for the liquid bulk.
Figure (4.1): LDV implementation for measuring velocity in a (r-z) plane, (Plane B).
The LDV is a Dantec Fiberflow system with two green and blue laser beams,
wavelengths of 514.4 nm and 488nm respectively. The laser source in the LDV is 4W Argon with a focal length of 600 mm (Spectra Physics). The LDV measurements
were performed through an r-z plane (Plane B) at (90o) to posterior baffle and (30o) to
two lateral baffles (See Figure 4.1), where for each measurement, about (4000 to
5000) samples are applied with a maximum measurement time of 300 seconds. The
tracer particles that were used to seed the liquid are (Iriodin 111) (dp =15m).
The PIV were performed in order that flow patterns can be determined for all vessel
areas. the vessel dimensions are considered somehow large in comparison with same
kinds of experimentations, so it is preferred to conduct flow field measurements by
- 148 -
Figure (4.2) PIV implementation for measuring velocity in a (r-z) plane, (Plane A)
For single (liquid) phase flow the used propeller rotation speed was kept at (5 rps) to
ensure reasonable vector displacement between two frames and since upper limit of
rotation speed due to turbine presence and consequence droplets impingement
position does not exist. The applied software for acquisition system is (DaVis 8.0.5
flow master), while image processing is (PIV imager pro x2M system).
- 149 -
Figure (4.3): The axial 4-bladed, reversible 45o pitched and twisted blade propeller
RTP propeller (HPM204D, Milton Roy Mixing)
- 152 -
- 153 -
Figure (4.4): The (r - z) flow field map for the RTP-D propeller with draft tube configuration in the entire vessel.
- 154 -
- 155 -
Figure (4.5): The (r - z) velocities contour map for the RTP-D propeller with draft tube configuration in the entire vessel.
- 156 -
0.05
0
Vz/Vtip, (-)
-0.05
-0.1
-0.15
-0.2
-0.25
-0.3
0.2
0.4
0.6
0.8
r/R, (-)
Figure (4.6): Dimensionless axial velocity profile at different liquid levels (below
impeller z/h=0.55, above impeller z/h=0.62) for RTP-D with draft tube, Cpr /Tv =0.2,
N= 5 rps in two different measuring planes; plane A at 00 and plane B at 900 from the
posterior baffle respectively.
The axial velocity in the discharge region of the impeller looks higher in the
measurement plane A than those in plane B and also it has a different profile more
than in inward region above the impeller. This difference is due to the effect of the
- 157 -
(a)
(b)
Figure (4.7): (a) The (r - z) flow field map for the RTP-D propeller without draft tube
configuration in the vessel core, (b) The (r - z) velocities contour map for the RTP-D
propeller with draft tube configuration in the vessel core.
As a result of the presence of the cone at the base of the vessel, a radial flow in the
discharge stream begins to appear relatively and gradually, where the majority of the
discharge flow is radial near the bottom. The other secondary circulation below the
impeller exists as in the propeller and draft tube configuration but with less intensity.
(Mavros et al. 2002) reported a close observation for pitch blade impeller (Mixel TT
- 159 -
- 160 -
0.05
0
-0.05
Vz/Vtip, (-)
-0.1
-0.15
-0.2
-0.25
-0.3
-0.35
0.2
0.4
0.6
0.8
r/R, (-)
Figure (4.8): The dimensionless axial velocity profile at different liquid heights for
RTP-D propeller without draft tube (below impeller z/h=0.55, above impeller
z/h=0.62), Cpr /Tv =0.2, N= 5 rps in two different measuring planes; plane A at 00 and
plane B at 900 from the posterior baffle respectively.
It has to be mentionned also that for the propeller alone configuration, sometimes with
the rotational speed N=4.6 rps, a swilling appears from the liquid surface and extends
to the impeller baldes tips (See Fig. 4.9). This phenomineon becomes more
continuous starting from 5.8rps, where the vortex is formed for the implemented
geometry. Despite of using the baffled tank, this vertex occurs because of low ratio of
dpr/Tv = 0.15 is implemented, where the baffles (Lbf/Tv=0.1) are somehow far away
- 161 -
Figure (4.9): The creation of vortex in the down-pumping mode with draft tube.
- 162 -
(a)
(b)
Figure (4.10): (a) The (r - z) flow field map for the RTP-U propeller with draft tube
configuration in the vessel core, (b) The (r - z) velocities contour map for the RTP-U
propeller with draft tube configuration in the vessel core.
- 163 -
- 164 -
0.3
0.25
Vz/Vtip, (-)
0.2
0.15
0.1
0.05
0
-0.05
-0.1
0.2
0.4
0.6
r/R, (-)
0.8
Figure (4.11): The dimensionless axial velocity profile at different liquid heights for
RTP-U propeller with draft tube (below impeller z/h=0.55, above impeller z/h=0.62),
Cpr /Tv =0.2, N= 5 rps in two different measuring planes; plane A at 00 and plane B at
900 from the posterior baffle respectively.
Similar to the propeller and draft tube configuration in down-pumping direction, the
effect of the draft tube baffles appears and alters the velocities in the inward zone
below the impeller and inside the draft tube, but on contrary of down-pumping
condition the alteration occurs in the discharge zone. Beside the fact that the
divergence is due to the difference of the measurement planes A and B cannot be
neglected.
- 166 -
(a)
(b)
Figure (4.12): (a) The (r - z) flow field map for the RTP-U propeller without draft
tube configuration in the vessel core, (b) The (r - z) velocities contour map for the
RTP-U propeller with draft tube configuration in the vessel core.
The up-pumping RTP-U propeller alone generates high radial velocities in the
discharge flow especially of the upper part of the vessel core near liquid surface. The
maximum radial velocity in this part is 0.3Vtip, whilst in the lower part of the vessel
core the radial velocities are of the order of 0.156Vtip. These values are relatively
close to those reported by and (Mishra et al., 1998) for the APV B-2 of 0.3Vtip in the
impeller discharge and 0.15Vtip in the lower zone below the impeller, while the
reported radial velocities are higher than of those reported by (Gabriele et al., 2011)
for PBT of 0.15 Vtip in the discharge flow and 0.12Vtip in the lower zone below the
impeller.
- 167 -
0.3
Vz/Vtip, (-)
0.25
0.2
0.15
0.1
0.05
0
0.2
0.4
0.6
0.8
r/R, (-)
Figure (4.13): The dimensionless axial velocity profile at different liquid heights for
RTP-U propeller without draft tube (below impeller z/h=0.55, above impeller
z/h=0.62), Cpr /Tv =0.2, N= 5 rps in two different measuring planes; plane A at 00 and
plane B at 900 from the posterior baffle respectively.
Figure 4.13 illustrates the quantitative analysis of the normalized axial velocities
above and below the up-pumping RTP-U propeller without the draft tube in the two rz measurement planes A and B at the vertical heights of z/h= 0.62 and z/h=0.55
respectively. As it is noticed the axial velocities profile above and below the propeller
in the planes A and B are generally similar, with highest axial velocities above the
impeller of the order 0.31Vtip and 0.3Vtip in the planes A and B respectively. These
values are close to those of propeller and draft tube configuration. Below the RTP-U
propeller the highest axial velocities are of the order of 0.21 Vtip and 0.2 Vtip
respectively for the planes A and B. With these values of the axial velocities it is quite
clear in this case where the draft tube effect no longer exists, that the produced
velocities especially in planes below the impeller are relatively similar. The axial
velocities below the RTP-U propeller are significantly different from those in the
propeller and draft tube configuration case, where maximum axial velocities of
0.27Vtip and 0.18Vtip are achieved for the planes A and B respectively.
Cim/Tv
dim/Tv
0.2
Np
Np
(With DT*)
(Without DT)
0.15
0.53
0.54
0.2
0.15
0.61
0.61
HFa -D
0.50
0.33
0.20
0.27
--
PBTb -U
0.3
--
PBTc -D
PBTc -U
0.25
0.33
0.5
0.33
0.5
--
0.25
0.43
0.45
0.45
---
0.33
0.5
--
0.25
0.45
--
1.16
0.18
1.4
--
0.33
0.33
--
0.33
0.33
0.3
0.35
--
APV-B2d -D
APV-B2d -U
MTTc -D
MTTc -U
EEe -D
EEe -U
PBTf U
(Lf =0.54 m)
HF4g D
(Lf = 0.125)
(Lf = 0.167)
Propellerd
*
0.40
0.41
0.46
1.3
1.2
1.1
1.93
2.58
0.95
0.91
0.74
0.67
2.1
1.7
0.89
c
DT : Draft tube ; Kumaresan and Joshi (2006); Chapple et al. (2002) ; Aubin et al. (2001)
d
Mishra et al. (1998); e Zhu et al. (2009); f Broz et al. (2004); g Kumaresan et al. (2005)
Table 4.1 illustrates clearly that the power characteristics for all types of pitched blade
impellers are very dependent on the impeller to tank diameter ratio or and on the offbottom clearance to tank diameter. As observed, the power number changes for the
same pitched blade impeller when the geometry is changed, so it is difficult to
compare the power characteristic between the pitched blade impellers with different
geometrical ratios. This agrees with (Chapple et al., 2002) who found that the power
number for pitched blade impellers is highly related to impeller position in the tank.
(Chapple et al., 2002) explained this changing as the flow at the impeller interacts
strongly with the proximity of the tank walls, so changes in the position of the
- 169 -
NQp
dim/Tv
0.2
0.15
0.2
0.15
0.50
0.33
0.20
0.27
0.33
0.5
0.25
0.43
0.45
0.45
0.33
0.5
0.25
0.45
1.16
0.18
--
0.5
--
0.33
0.33
0.33
0.41
0.65
--
--
0.33
0.67
0.89
1.3
NQc
NQp
Cim/Tv
1.186
(With DT)
1.22
(With DT)
-0.96
-0.91
1.82
1.19
0.975
1.13
1.35
---
(With DT*)
(Without DT)
0.56
0.51
0.55
0.495
-0.59
-0.75
0.68
0.83
0.58
0.67
0.61
0.73
0.72
------
0.73
DT : Draft tube ; Kumaresan and Joshi (2006); Chapple et al. (2002) ; Aubin et al. (2001)
d
Mishra et al. (1998); e Zhu et al. (2009); f Bro et al. (2004); g Kumaresan et al. (2005)
For the RTP propeller the pumping numbers in the up-pumping mode are slightly
higher than those in the down-pumping mode for both configurations with and
- 170 -
Table (4.3): The agitation Index for RTP-D and various pithed blade impellers, (this
work and literature).
Impeller
Type
RTP-D
(Lf = 0.1m)
PBTa -D
PBTa -U
MTTa -D
MTTa -U
MTTb -D
b
Lightnin A310 -D
*
DT : Draft tube ;
Cim/Tv
dim/Tv
0.2
0.15
0.33
0.5
0.33
0.5
0.33
0.5
6.7
16.0%
22.5%
14.1%
13.7%
12.0%
0.5
6.7
8.0%
0.33
Ig
(rps)
3.47%
(With DT*)
The liquid volume quantification or distribution for the down-pumping mode with the
draft tube is illustrated in the Figure 4.14, where the liquid volumes distribution in the
vessel is related to determine ranges of the composite mean local 2-D velocities in the
vessel. The composite mean velocities are illustrated as a 3-D surface chart as seen in
Figure 4.14; the liquid volume movement is elucidated in the vessel except in the
swept volume of the impeller, in the volume of the cone at the vessel base, and in the
one of the shaft and baffles edges. From figure 4.14 it is found that in the majority of
the liquid volume, the mean velocity is low, smaller than of the order 0.11Vtip except
in three zones. The most elevated is the discharge and intake zone of the RTP-D
where mean velocity is less than 0.34Vtip excluding a very small part in this with
mean local velocities of the order between 0.34Vtip and 0.45Vtip. The other two zones
near the vessel bottom and walls have lower local velocities less than 0.23Vtip; these
two zones are generated due to the interaction of the flow induced by the RTP-D
propeller with both the vessel bottom and walls.
- 173 -
- 174 -
0.5-0.6
0.4-0.5
0.3-0.4
0.2-0.3
0.6
0.1-0.2
0.5
0-0.1
- 175 -
vm, (m S-1)
0.4
390
360
330
300
270
240
210
0.3
0.2
0.1
0
180
280
260
150
240
220
120
200
180
160
90
140
120
60
100
80
60
30
40
20
Figure (4.14): The distribution of the liquid volume in the vessel related with mean local composite velocity ranges for RTP-D propeller with
draft tube configuration.
- 176 -
120
tm, (s)
100
80
60
40
20
0
1.50
2.50
3.50
N, (rps)
4.50
5.50
Figure (4.15): The effect of RTP propeller rotational speed on the mixing time for
different flow and geometry conditions, h/dpr=1.47, Cpr/Tv = 0.2.
- 177 -
240
Down-pumping without DT
Up-pumping without DT
Ntm, (-)
220
200
180
160
140
120
20
30
40
50
60
-3
Re*10 , (-)
70
80
Figure (4.16): The relation of Reynolds number with the dimensionless mixing time
for different pumping modes and geometry configurations , h/dpr = 1.47, Cpr /Tv = 0.2.
4.7. Conclusions
In this chapter the flow patterns and r-z velocity vector fields and contour maps in a
single phase for the RTP propeller agitated vessel were measured by applying both
PIV and LDV for up-pumping and down-pumping operation modes with and without
draft tube configuration. The LDV measurements were restricted inside the draft tube
and in the measurement plane at 90o from the posterior tank baffle.
For the down-pumping mode, mainly one recirculation loop is generated with RTP-D
propeller with two notable secondary loops; one below the propeller and the second
- 178 -
- 179 -
- 180 -
- 181 -
- 182 -
Chapter Five
Hydrodynamics and Flow Pattern in Aerated Agitated Tank
5.1. Introduction
In this chapter the flow patterns and velocity fields for the aerated vessel have been
investigated. Two phases air-water are occurring since the water droplets are sprayed
in the air and then with their impacting at the water surface a sufficient amount of air
bubbles are entrapped and entrained inside the water bulk. To evaluate the aeration
performance of the system and to estimate the contribution of each element (Turbine,
RTP propeller and draft tube) on the process different experimental sets were
performed.
The PIV measurements for the flow pattern and mean velocity were carried out in the
r-z plane (Plane A) same as described in chapter four, where r-z measurement plane
was positioned directly on posterior baffle (0o) and at (60o) from the two lateral
baffles from front face (See Figure 4.1). The number of image pairs was (300) and
exposure time delay (t) between each two sequent images was (2500 s) depending
on operation condition such as upper turbine tip velocity and image scale factor
adapted during the acquisition. A pixel interrogation of (32%x32%) with (50%
overlap) was applied for the velocity measurements.
The conical shape turbine was placed at the water surface level (C/Tv=0.313, D/Tv=
0.238, h=0.28m). The used vessel is the 3-baffled cylindrical vessel as described
earlier (Tv =0.8 m). The four 45o pitched and twisted blade propeller RTP effect was
studied at clearance (Cpr/Tv = 0.2 m) with and without the draft tube (df =0.15 m, Lf =
0.1 m) that was positioned in the vessel with clearance (Cf = 0.1m). The RTP
propeller and the draft tube were removed in a part of the measurements. The
measurements were performed with turbine rotation speed of (N= 2.08 rps) to ensure
staying always in turbulent regime for created flow field (Re = 75000) and prevent
excessive air bubbles presence with higher rotation speed, which leads to
miscalculation of r-z mean velocities vector plot by PIV due to the high reflection of
laser bean occurred by these bubbles. The PIV measurements were performed in the
entire tank to identify the effect for the both the draft tube and/or RTP propeller
configurations on the produced mean velocity field and liquid flow patterns.
5.2. Flow Pattern and Mean Velocity Field in the Aerated Tank
The flow pattern and mean velocity fields in the aerated tank for the whole system is
illustrated in the Figures 5.1 and 5.2. From these Figures the time average velocities
vector plot and contour inside the aerated tank at the r-z plane at 0o of the posterior
- 183 -
- 187 -
Figure (5.1): The (r - z) flow field map for the turbine and the RTP-U propeller with draft tube (Whole system) configuration in the entire vessel.
- 188 -
- 189 -
Figure (5.2): The (r - z) velocities contour map for the turbine and the RTP-U propeller with draft tube (Whole system) configuration in the
entire vessel.
- 190 -
- 191 -
Figure (5.3): The (r - z) flow field map for the turbine alone configuration in the entire vessel.
- 192 -
- 193 -
Figure (5.4): The (r - z) velocities contour map for the turbine alone configuration in the entire vessel.
- 194 -
- 195 -
Figure (5.5): The (r - z) flow field map for the turbine and the RTP-U propeller configuration in the entire vessel.
- 196 -
- 197 -
Figure (5.6): The (r - z) velocities contour map for the turbine and the RTP-U propeller configuration in the entire vessel.
- 198 -
NQc
NQp
Whole system
2.42
0.49
Turbine + RTP-U
2.33
0.44
Turbine alone
2.24
0.40
Table 5.1 shows that for the whole system configuration, the turbine pumping number
and the system circulation number are higher than those for both with and without
draft tube and turbine alone configurations. The turbine pumping number is slightly
increased with draft tube presence than the turbine alone configuration. From these
values it is found that the presence of the draft tube assists the RTP-U propeller
performance by redirecting part of the RTP-U propeller radial discharge flow
discharge flow upward as it was seen in the up-pumping in the Mixing mode.
The RTP propeller function seems to have an obvious effect on the turbine pumping
number and the circulation number of the system, where these parameters are
relatively elevated, when the RTP-U propeller was positioned below the turbine, this
elevation occurs due to the improvement of the intake flow of the turbine that
achieved by the propeller action.
5.6. Agitation Index and Liquid Quantification for the Whole System
Configuration
The agitation index (the volume-weighted average velocity) for the whole system
configuration was calculated according to the equations (2.15 to 2.17). The agitation
index is equal to 5.26%, this value is higher than for the non-aerated vessel with
down-pumping mixing mode agitation index of 3.47% (N=5 rps) despite that for the
aerated mode the rotation speed is higher 2.08 rps and the presence of the air bubbles
may hinder the mixing to be efficient easily. The agitation index in the aerated mode
- 199 -
- 200 -
0.5-0.6
0.4-0.5
0.3-0.4
0.2-0.3
0.6
0.1-0.2
0.5
- 201 -
vm, (m s-1)
(vm) levels
0-0.1
0.4
390
360
330
300
270
240
210
0.3
0.2
0.1
0
180
280
260
150
240
220
120
200
180
160
90
140
120
60
100
80
60
30
40
20
Figure (5.7): Distribution of the liquid volume in the vessel related with mean local composite velocity ranges for whole system configuration.
- 202 -
Ntm, (-)
80
60
40
20
1.60
2.10
2.60
N, (rps)
3.10
3.60
Figure (5.8): The relation between impellers rotation speed and the dimensionless
mixing number Ntm for whole system configuration in the aerated vessel.
The relationship of mixing time tm with the impellers rotational speed N is illustrated
in Figure 5.9. For the whole system configuration, the mixing time appears to be very
dependent on the elevation of impellers rotation speed like the single phase condition
(Mixing mode), where it decreases with increasing of the impellers speed. As a
consequence this indicates that the mixing time in the aeration mode is highly effected
by the flow pattern and the bubbles presence in the liquid phase.
For the turbine alone configuration, both the draft tube and the propeller are removed.
Figure 5.9 shows the mixing time tm is highly affected by the elevation of N. The
mixing time tm is lessened remarkably from 33 sec to 17 sec when the rotation speed
was doubled (1.67-3.33).
- 203 -
30
tm, (s)
25
20
15
10
1.50
2.00
2.50
N, (rps)
3.00
3.50
Figure (5.9): The effect of impellers rotation speed on the mixing time for two
different geometrical configurations, h/D=1.47, Cpr/T= 0.2, C/T = 0.313 (Aerated
mode).
- 204 -
Ntm, (-)
80
60
40
20
20
40
60
80
Figure (5.10): The relation between the Reynolds number and the dimensionless
mixing time of the aerated vessel for whole system and turbine alone configurations,
h/D= 1.47, Cpr/Tv = 0.2, C/Tv = 0.313.
- 205 -
31
tm, (s)
26
21
16
11
1.60
2.10
2.60
3.10
3.60
N, (s)
Figure (5.11): The relation between impellers rotation speed and the spacing between
impellers; h/D= 1.47, C/Tv = 0.313, (Aerated mode).
70
S/W=1.42
h/D=1.53
Ntm, (-)
60
50
40
30
20
20
25
30
35
40
-3
Re *10 , (-)
45
50
Figure (5.12): The relation between the Reynolds number and the dimensionless
mixing time of the aerated vessel for three turbine blade submergence levels, Cpr /T=
0.2, C/T = 0.313.
(5.4)
From equation 5.4 it can be found that the decreasing of the turbine blade
submergence lengthens the mixing time within the range tm (S)-0.16. This means
when the submergence of the turbine blades decreased from 24 to 14 mm, the mixing
time would lengthen by the factor 1.1.
The comparison between the model and the experimental results is shown in Figure
5.13. The agreement is quite satisfying. It is somehow hard to compare the constants
in the equation (5.4) with the previous models because of the difference in the
geometry and mode of operation applied, where the already developed aerated
condition mixing time models didnt applied the water droplets projection principles
to achieve the aimed aeration. But generally for the negative exponent of power
number it seem to be logical and refers to that mixing time number is decreased with
higher power consumption (Nienow, 1997). The (S/W) ratio exponent refers to the
mixing time number is decreased when higher turbine blades submergence (higher
liquid level) which doesnt agree with the model that developed by (Kang et al., 2001)
with considering the difference in geometry with their system.
- 208 -
Ntm, Experimental
60
50
40
30
20
10
0
10
20
30
40
Ntm, Predicted
50
60
70
- 209 -
4.0
Turbine alone
Np, (-)
3.5
3.0
2.5
2.0
1.5
1.0
20
30
40
Re*10-3,
50
(-)
Figure (5.14): The effect of Reynolds number on power number for the two different
configurations (Whole system and turbine alone); h/D= 1.47, Cpr/Tv = 0.2, C/T =
0.313, (Aerated mode).
250
Whole system
Turbine alone
P/V, (Watt/m3)
200
150
100
50
4
N,
(s-1)
Figure (5.15): The relation between impellers rotation speed and the consumed power
per water volume for the aerated mode, h/D= 1.47, Cpr/Tv = 0.2, C/Tv = 0.313.
- 210 -
- 211 -
- 212 -
- 213 -
- 214 -
The general global purpose of this thesis is to evaluate and interpret the capacity for
the novel technology applied with the investigated surface aerator related to the
possibility to perform the aeration or the mixing by the same equipment.
The first main objective of this study consists in the characterization of the aeration
and the mixing modes capacities related to the operational and system configuration
parameters. The second major objective involved the investigation of the flow pattern
and velocity field takes place in the tank using advanced measuring techniques for
both up-pumping aeration mode and down-pumping mixing mode.
The implemented surface aerator has a distinctive feature of easy putting in place
between the aeration and the mixing modes by reversing the sense of the rotation
where the oversaturation condition is frequently occur in the liquid bulk with
existence of poor aerated zones. This interchanging is performed by applying the
decoupling (clutching) system that enables the lower propeller to carry out the mixing
with setting aside the turbine from the rotation.
A brief description was presented for the nowadays aeration systems for the water
and wastewater treatment with exposing the types and efficiencies for each system,
with emphasizing on the surface aeration.
From the literature review it is found that many modifications and developments are
made on the surface aerators to improve the oxygen mass transfer rate and to reduce
the power consumed. Many parameters have been studied like hydrodynamic, power
consumption, geometric configurations, oxygen mass transfer rate in the water,
modeling and scale up process. It can be presumed that in general the limiting factors
for successful surface aeration can be resumed as follows:
- The turbulent regime should be ensured in the entire water treatment tank that
contains sufficient dissolved oxygen entrained from atmospheric air.
- The generated water flow by surface aerators must be sufficient to reach all parts of
the treatment tank.
- The aeration process accomplished with the aerator turbine must be accompanied by
mixing and agitation performance to ensure an efficient distribution of the dissolved
oxygen.
- These constrains can be overcame by using an appropriate surface aerator with
suitable pumping capacity to handle efficiently and effectively the large quantities of
water.
- 215 -
- 218 -
The changing between the aeration mode and the mixing mode is applied using a
coupling (clutching) system. Other ways of exchanging between the mixing and
- 219 -
- 220 -
- 221 -
- 222 -
Appendix I
A Review for Biological Wastewater Treatment with
Activated Sludge
The clean water availability question becomes more and more as a developing
important issue especially during the last decades. The necessity for wastewater
treatment regarded as a very important compromise among the available solutions to
preserve the water resources and to reduce the water pollution. The aerobic biological
purification that implied in the wastewater plants consists of many steps; generally it
can be divided into two main treatments: urban or domestic wastewater treatment and
industrial wastewater treatment.
- 222 -
- 222 -
Figure (I.2): Domestic wastewater treatments steps (Try and Price, 1995).
Figure (I.3): Development and control of biomass in the activated sludge process
(Gary, 2004).
The activated sludge treatment depends on several process factors, especially the
flocculation as described earlier. The successful agglomeration of the various organic
materials is related to their flocculent growth when mixed with activated sludge, the
settled biological materials from this operation are removed later by the separation
from the wastewater. The flocs can be determined as a biomass cluster of several
million heterotrophic bacteria bounded together with some inert organic and inorganic
material (Gary, 2004). The process of activated sludge includes several
implementations. Among of them there are nutrient removal processes:
- 222 -
A. Nitrogen Removal
Nitrogen presents in different forms in the wastewater due to various oxidation states,
which can be changed from one state to state to another a according to the different
physical and biochemical conditions involving wastewater. The source of organic
nitrogen may come from domestic wastes, farming activities and food industry. Most
of bio-nitrogen in wastewater is present as ammonia of unstable organic compounds,
but in urban wastewater, it contains excess nitrogen of the microbial requirement to
oxidize the amount of carbon present, so only a part of nitrogen is being combined
into the biomass (Crites et al., 2006).
The removal of the ammonia prevents its toxic effects on the aquatic life systems and
the deterioration of water quality (Gupta and Sharma, 1996). In some cases the N2O
gas could be emitted during nitrates removal from waste water as a secondary
product. The emitted gas is controlled by additional techniques like intermittent
aeration (Park et al., 2000).
The nitrogen is used by the biomass effectively, where ammonia is oxidized by
autotrophic nitrifying bacteria to nitrite. The biological treatment of ammonia by
oxidation is occurred in two stages each achieved by different type of nitrifying
bacteria; they utilize the ammonia as an energy source in the first step, or the nitrite
for same purpose in the second step as following:
NH4+ + 1.5 O2
0.5 O2 + NO3-
The nitrification implemented with slow growing bacteria therefore, it proceeds at low
rate then, in order to achieve an acceptable treatment it needs maintaining more longer
retention time for the activated sludge in the tank or the basin by keeping maximum
contact with activated sludge to ensure maximum nitrification rate and to avoid washout of the nitrifiers. It is preferred in the nitrification process with activated sludge
maintaining the dissolved oxygen concentration about (2.0 mg/l), where it is remarked
that the nitrification doesnt occur below (0.5 mg/l), where with high dissolved
- 222 -
B. Phosphorus Removal
The phosphorus exists in the wastewater in form of orthophosphates (PO4 3-, HPO4 2-,
H2PO4 -2, H3PO4), polyphosphates and organic phosphates, where its source is
generally comes from fertilizing implementation. The total phosphorus concentration
in wastewater is about (5-20 mg P/l), where only (1-2 mg P/l) is removed by
biological treatment (Sincero and Sincero, 2003). The phosphorus in wastewater
stream is removed by both chemical and biological methods (Figures I.5 and I.6);
where the main biological phosphorus removal is anaerobic-oxic system. The bacteria
- 222 -
2CaCO3 + 2H2O
Then the calcium ions combined with orthophosphate to form insoluble and
gelatinous calcium hydroxyapatite.
5Ca 2+ + 4OH - + 3HPO4 2-
Ca5(OH)(PO4)3 + 2H2O
Figure (I.5): The Two stage biological phosphorus treatment (Gary, 2004).
- 222 -
- 222 -
Figure (I.7): Covered aerated tank using pure oxygen (Mueller et al., 2002).
Usually the aeration is performed with active sludge in tanks or basins, their depth
vary from (3.5 6.0 m) to about (20-30 m ) equipped with air diffuser below the
liquid and inject the air or oxygen , where these types are preferred to treat industrial
wastewater especially with proper type of aeration tools was chosen. An effective
process that preferred to treat synthetic wastewater can remove more than 95% of
COD from the effluent (Polprasert and Raghunandana, 1985; Sincero and Sincero,
2003).
- 222 -
- 222 -
- 222 -
- 235 -
- 236 -
Appendix II: Derivation of Spray Flowrate for the Surface Aeration System
Appendix II
Derivation of Spray Flowrate for the Surface Aeration System
Fluid mechanics analysis of the surface aerator system is performed to find a logical way
to calculate the water spray flowrate Q, where direct experimental measurement of Q
would be extremely difficult and expensive to accomplish this objective (McWhirter et
al., 1995). This calculation depends on the basic physical principles that occurred in the
surface aeration system related to the volumetric flow rate Q of the water spray
discharged by the surface turbine blades into atmospheric air.
The spray flowrate can be de determined by consideration of the overall conservation of
energy applied to a surface aeration system and the measured overall power input to the
aerator (Baylar et al., 2001; McWhirter et al., 1995). (McWhirter et al., 1995) envisioned
the surface aeration system as impeller blades rotating in the free liquid surface of a
relatively large body of water, they considered the surface aeration turbine impeller
blades accelerate a flow rate of liquid Q from a relatively low velocity at the inlet to the
impeller blades up to a relatively high velocity at the discharge from the tip of the aerator
impeller blades.
The general power input to the surface aeration system is a mechanical energy transfer to
the water bulk and then it is converted to kinetic energy in the water spray as explained in
chapter two, this power input can be presented as the following equation (Baylar et al.,
2001).
The fluid friction losses for the water flow across the surface aerator turbine blades is
considered negligible, the efficiency is assumed 100% as a result all the mechanical
energy is converted with 100% to a kinetic energy, so the efficiency, , equals to 1. Ym
represents the height between the upper point of the water spray maximum point.
By applying the general energy balance of Bernoulli equation between two points 1 and
2 as illustrated in Figure (II.1), where the total energy of a fluid in motion consists of the
following components; potential, pressure and kinetic energies (Holland and Bragg,
1995).
- 237 -
Appendix II: Derivation of Spray Flowrate for the Surface Aeration System
Figure (II.1): The energy balance points for the spray discharge flowrate.
Equation II.3 is reduced to equation II.4 as the pressures for the two points are same (the
atmospheric pressure).
Since the spray velocity v2 at the upper point (point 2) of the spray is equal to zero
The water is propelled from the tip of the turbine blades will be subject to the same
kinematic laws of physics as any free falling projectile (McWhirter et al., 1995), so the
velocity at the impingement point equals the velocity at the turbine blades tip.
Equation II.8 is written as follows, for power input in kW and flowrate in liter/h
Where the constant 1.38910-7 includes the unit conversion factor and it depends on
system geometry as the discharge flowrate depends on the turbine geometry.
- 238 -
Appendix II: Derivation of Spray Flowrate for the Surface Aeration System
The flow rate Q is assumed to be constant along the spray, where the flow rate Q into the
atmospheric air at a relatively high total discharge velocity Vsp from the tip of the turbine
blades can be applied at the impingement point at the water surface. With negligible
friction loss all of the aerator shaft power input goes into accelerating the flow rate of
liquid spray Q up to the discharge velocity Vsp (McWhirter et al., 1995), so equation II.8
can be applied to relate the total power input with the discharged spray flowrate and
spray velocity Vsp at impingement point.
- 239 -
- 240 -
List of Symbols
List of Symbols
Latin symbols
a
a term in kla
ad term in klad interfacial area per unit volume in water droplets, (m2 m-3)
Ad
Cd
Cdt
Ce
Cl
CLS
CLt
Cpr
Cs
Csd
Cs10oC
Ct
Ctr
Co
df
dpr
Emd
List of Symbols
Esp
(Esp)20
Fr
Frsp
Ig
kla
klad
kla10oC
klaT
Np
NQp
NQc
OTRb
oxygen transfer rate in water bulk mass transfer zone, (kgO2 h-1)
OTRsp
oxygen transfer rate in water spray mass transfer zone, (kgO2 h-1)
Pin
the input power to the turbine (that assumed 100% converted to kinetic
energy to propel water spray into the air), (k watt)
(kla20(2/g)1/3),
in
equations
List of Symbols
Qc
Qp
Qpr
Qpz
qsp
rcent
Re
Resp
Rm
Rsp
radius of the water spray calculated from the shaft center, (m)
SOTRb
SAEb
time, (s)
tf
tm
To
Toe
Tv
Va
the mean velocity of the water spray radial and tangential velocities
components, (m s-1)
Vd
VL
List of Symbols
vm
Vsp
Vr
Vy
Vtip
vr
vz
Vy
Voy
initial axial component of the water spray velocity at time tf =0, (m s-1)
Vr'
Vy'
WT
Ym
zcent
Greek Symbols
- 244 -
List of Symbols
EE D, U
HF-D
HF4-D
MTT D, U
PBT D, U
RTP D, U
- 245 -
- 246 -
List of Figures
List of Figures
Chapter One
Figure (1.1): Plate diffuser aerator,(Permox H ceramic plate diffuser, (Supratec Co.
Ltd.)
10
Figure (1.2): (a)
Tube diffuser aerator,
(b)
Disc diffuser aerator,
(Gemgate GmbH)
11
Figure (1.3): Static tube diffusers (Process Engineering s.r.l)
12
14
Figure (1.6): Flow diagram for uncovered pure oxygen aeration (Mueller et al.,
15
2002)
Figure (1.7): Horizontal flow aspirating aerator, Aspirator, (AIRE-O2 Aeration
Industries International)
16
Figure (1.8): Low speed vertical flow aerator Up-ward flow,Praxair Technology.17
Figure (1.9): Low speed vertical flow aerator (Down-ward flow), (a) Turboxal (Aire
Liquide), (b) Praxair (Praxair Technology)
17
Figure (1.10): Low speed horizontal flow aerator (Twin mini rotor aeration,
(Botjheng Water Ltd.)
18
Figure (1.11): High speed surface aerator, Aqua turbo (AER-AS), (AQUATURBO
SYTEMS inc.)
19
Figure (1.12): The surface aeration regimes applying air entrainment from free liquid
surface (A) Direct entraining of the atmospheric air, (B) Spray
formation and entraining the air with droplets impingement at the
surface,(Patwardhan and Joshi, 1998)
21
Figure (1.13): The histogram distribution of measured aeration efficiency for (111)
low speed surface aerators in field, the average is 1.49 kgO2/kWh,
(Heduit and Racault, 1983b)
27
Figure (1.14): The relation between surface aeration impeller blades submergence
and oxygen transfer rate for different blades number, (where H is the
liquid level in the tank) (Backhurst et al., 1988)
30
Figure (1.15): Relation between sorption number (Y) with the Froude number (Fr) for
the conical shape turbine surface aerator (Zlokarnik, 1979)
36
- 247 -
List of Figures
Figure (1.16): Bubble distribution of surface aerator system of rotation speed N=110
rpm, liquid level h=0.66 m. (Lee et al., 2001)
41
Figure (1.17): The relation between the rotation speed and mixing time for surface
aeration dual impeller system, (Kang et al., 2001)
42
Figure (1.18): The relation between the Reynolds number and power number for
three different liquid condition ns, clean water and two types of
44
activated sludge mixed liquid (Takase et al., 1982)
Chapter Two
Figure (2.1): The schematic diagram of the experimental apparatus
52
52
56
57
58
59
Figure (2.9): The Particle Image Velocimetry (PIV) testing apparatus, (1) Laser
source, (2) Recording camera, (3) Tested tank
60
Figure (2.10): Light scattering by the (10 m) glass particles in the water (Raffel et
al., 2007)
60
Figure (2.11): The pumping number and pumping flowrate measuring volume
65
Figure (2.12): The 3-D liquid volume cell and the vessel volume grids, (GarciaCortes et al., 2006)
67
Figure (2.13): The schematic diagram of the polarographic probe
68
68
69
- 248 -
List of Figures
Figure (2.16): (a) Experimental oxygen probes positions in the liquid volume,
(b) Experimental oxygen probes positions in the liquid volume, N=2.5
rps, (h/Tv) = 0.35
72
Figure (2.17): The repeatability of the mass transfer (kla) experimental results for
both oxygen probes, (N= 2.08 rps)
73
Figure (2.18): Optical probe response time verification
74
74
Figure (2.20): The probe response time relationship verification for both theoretical
and experimental DO values for the optical probe
76
Figure (2.21): The probe response time relationship verification for both theoretical
and experimental DO values for polarographic probe
76
Figure (2.22): The droplets zone mass transfer coefficient (klad) measurement
83
Figure (2.23): Schematic diagram illustrates the surface aeration water droplets spray
from turbine blades till the impingement point with some important
relevant dimensions
84
Figure (2.24): Droplets radial and tangential velocities propelled in horizontal plane
by turbine blades
85
Chapter Three
Figure (3.1): Different tested geometric configurations, (a) Whole system, (b)
Turbine alone, (c) Turbine + Propeller, (d) Turbine + Draft tube 92
Figure (3.2): Impellers speed effect on the oxygen transfer coefficient with four
different geometrical configurations; (Whole System, without draft
tube, turbine alone and without propeller), D/Tv=0.24, C/Tv=0.31,
h/D=1.47, Cpr/Tv=0.2, S/W=1
93
Figure (3.3): The DO profile for different rotation speeds, h/Tv= 0.35, C/Tv=0.313,
Cpr/Tv=0.2, temperature =15 oC
94
Figure (3.4): (a) The limits of the tested turbine blades submergence, (b) Mass
transfer Coefficient kla relation with turbine blades submergence and
water height in the liquid bulk for three levels of rotational speed, (b)
Mass transfer Coefficient kla relation with the rotation speed for three
levels of turbine submergence. D/Tv= 0.24, C/Tv = 0.31,
Cpr/Tv=0.2.
96
- 249 -
List of Figures
Figure (3.5): Spacing effect on the oxygen transfer coefficient, D/Tv= 0.24,
C/Tv = 0.31, Sp (Reference) =7.6 mm Figure (3.6): radial flow rate by
propeller, (m3/s)
98
Figure (3.6): (a) The relationship between power number and Reynolds No., for
different turbine blades submergence levels C/Tv = 0.31, (b) The
relation between power consumption and rotation speed,
C/Tv = 0.31
99
Figure (3.7): The spacing
consumption
between
two
impellers
effect
on
the
power
100
List of Figures
Figure (3.17): The relation between the oxygen mass transfer coefficient k lad and
impellers rotation speed and in the spray zone, S/W = 1
116
Figure (3.18): The relation between the oxygen transfer rate OTRsp and impellers
rotation speed and, in the spray zone, S/W = 1
119
Figure (3.19): The influence of the impellers rotational speed on the spray zone
oxygen transfer rate OTRsp and in the bulk zone oxygen transfer rate
OTRsp contributions the overall transfer rate OTR, S/W = 1
120
Figure (3.20): Effect of the turbine blade submergence on spray mass transfer zone
dissolved oxygen concentration at rotation speed (N=2.5 rps) and
123
turbine clearance (C/Tv = 0.35)
Figure (3.21): The linear regression correlation of water spray zone aeration
efficiency by plotting (Cdt) versus (CLt) for various turbine blades
submergence, N = 2.5 rps
124
Figure (3.22): The relation between the spray zone aeration efficiency Esp and turbine
blade submergence ratio, N = 2.5 rps
125
Figure (3.23): The relation between the spray zone oxygen mass transfer coefficient,
klad with the turbine submergence and liquid level impellers rotation
speed, N = 2.5 rps
127
Figure (3.24): The relation between spray zone oxygen transfer rate OTRsp with the
turbine blades submergence and liquid level, N = 2.5 rps
128
Figure (3.25): The influence of the turbine submergence on the spray zone oxygen
transfer rate, OTRsp and in the bulk zone oxygen transfer rate OTRsp
contributions the overall transfer rate OTR, N = 2.5 rps
129
Figure (3.26): The effect of the turbine rotation speed on the spray mass transfer zone
dissolved oxygen concentration at h/D=1.47and C/T=0.35, (Turbine
alone configuration)
131
Figure (3.27): The linear regression of the spray zone aeration efficiency of the plot
(Cdt) versus (CLt) for various rotation speeds, (Turbine alone
configuration)
132
Figure (3.28): The relation between the impellers speed N and the spray zone aeration
efficiency Esp for both turbine alone and whole system configurations,
S/W = 1.0
133
Figure (3.29): The relation between N and the oxygen mass transfer coefficient klad
for both turbine alone and whole system configurations,
S/W = 1.0
133
- 251 -
List of Figures
Figure (3.30): The relation between spray zone oxygen transfer rate with the
rotational speed for the whole system and turbine alone configurations,
S/W = 1.0
135
Figure (3.31): The influence of the tank internals (Propeller and draft tube) on the
spray zone oxygen transfer rate, OTRsp contribution the overall transfer
rate OTR, S/W = 1.0
136
Figure (3.32): The Comparison between the (Esp)20, Predicted by the Correlation
model (Eq. 3.10) with the experimentally resulted (Esp)20
139
Figure (3.33): The Comparison between the (Esp)20, Predicted by the Correlation
model (Eq. 3.11) with the experimentally resulted (Esp)20
140
Figure (3.34): The Comparison between the (Esp)20, Predicted by the Correlation
model (Eq.3.16) with the experimentally resulted (Esp)20
142
Chapter Four
Figure (4.1): LDV implementation for measuring velocity in a (r-z) plane,
(Plane B)
148
Figure (4.2): PIV implementation for measuring velocity in a (r-z) plane,
(Plane A)
149
Figure (4.3): The axial 4-bladed, reversible 45o pitched and twisted blade propeller
RTP propeller (HPM204D, Milton Roy Mixing)
150
Figure (4.4): The (r-z) flow field map for the RTP-D propeller with draft tube
configuration in the entire vessel
153
Figure (4.5): The (r-z) velocities contour map for the RTP-D propeller with draft
tube configuration in the entire vessel
155
Figure (4.6): Dimensionless axial velocity profile at different liquid levels (below
impeller z/h=0.55, above impeller z/h=0.62) for RTP-D with draft
tube, Cpr /Tv =0.2, N= 5 rps in two different measuring planes; plane A
at 00 and plane B at 900 from the posterior baffle respectively
157
Figure (4.7): (a) The (r - z) flow field map for the RTP-D propeller without draft
tube configuration in the vessel core, (b) The (r - z) velocities contour
map for the RTP-D propeller with draft tube configuration in the vessel
core
159
- 252 -
List of Figures
Figure (4.8): The dimensionless axial velocity profile at different liquid heights for
RTP-D propeller without draft tube (below impeller z/h=0.55, above
impeller z/h=0.62), Cpr /Tv =0.2, N= 5 rps in two different measuring
planes; plane A at 00 and plane B at 900 from the posterior baffle
respectively
161
Figure (4.9): The creation of vortex in the down-pumping mode with draft tube 162
Figure (4.10): (a) The (r - z) flow field map for the RTP-U propeller with draft tube
configuration in the vessel core, (b) The (r - z) velocities contour map
for the RTP-U propeller with draft tube configuration in the vessel
core
163
Figure (4.11): The dimensionless axial velocity profile at different liquid heights for
RTP-U propeller with draft tube (below impeller z/h=0.55, above
impeller z/h=0.62), Cpr /Tv =0.2, N= 5 rps in two different measuring
planes; plane A at 00 and plane B at 900 from the posterior baffle
respectively
165
Figure (4.12): (a) The (r - z) flow field map for the RTP-U propeller without draft
tube configuration in the vessel core, (b) The (r - z) velocities contour
map for the RTP-U propeller with draft tube configuration in the vessel
core
167
Figure (4.13): The dimensionless axial velocity profile at different liquid heights for
RTP-U propeller without draft tube (below impeller z/h=0.55, above
impeller z/h=0.62), Cpr /Tv =0.2, N= 5 rps in two different measuring
planes; plane A at 00 and plane B at 900 from the posterior baffle
respectively
168
Figure (4.14): The distribution of the liquid volume in the vessel related with mean
local composite velocity ranges for RTP-D propeller with draft tube
configuration
175
Figure (4.15): The effect of RTP propeller rotational speed on the mixing time for
different flow and geometry conditions, h/dpr=1.47, Cpr/Tv = 0.2
177
Figure (4.16): The relation of Reynolds number with the dimensionless mixing time
for different pumping modes and geometry configurations , h/dpr =
1.47, Cpr /Tv = 0.2
178
Chapter Five
Figure (5.1): The (r-z) flow field map for the turbine and the RTP-U propeller with
draft tube (Whole system) configuration in the entire vessel
187
- 253 -
List of Figures
Figure (5.2): The (r-z) velocities contour map for the turbine and the RTP-U
propeller with draft tube (Whole system) configuration in the entire
vessel
189
Figure (5.3): The (r-z) flow field map for the turbine alone configuration in the
entire vessel
191
Figure (5.4): The (r-z) velocities contour map for the turbine alone configuration in
the entire vessel
193
Figure (5.5): The (r-z) flow field map for the turbine and the RTP-U propeller
configuration in the entire vessel
195
Figure (5.6): The (r-z) velocities contour map for the turbine and the RTP-U
propeller configuration in the entire vessel
197
Figure (5.7): Distribution of the liquid volume in the vessel related with mean local
composite velocity ranges for whole system configuration
201
Figure (5.8): The relation between impellers rotation speed and the dimensionless
mixing number Ntm for whole system configuration in the aerated
vessel
203
Figure (5.9): The effect of impellers rotation speed on the mixing time for two
different geometrical configurations, h/D=1.47, Cpr/T= 0.2, C/T =
0.313 (Aerated mode)
204
Figure (5.10): The relation between the Reynolds number and the dimensionless
mixing time of the aerated vessel for whole system and turbine alone
configurations, h/D= 1.47, Cpr/Tv = 0.2, C/Tv = 0.313
205
Figure (5.11): The relation between impellers rotation speed and the spacing between
impellers; h/D= 1.47, C/Tv = 0.313, (Aerated mode)
206
Figure (5.12): The relation between the Reynolds number and the dimensionless
mixing time of the aerated vessel for three turbine blade submergence
levels, Cpr /T= 0.2, C/T = 0.313
207
Figure (5.13): The comparison between experimental and predicted values of
dimensionless mixing time number (tmN) by equation (5.4)
209
Figure (5.14): The effect of Reynolds number on power number for the two different
configurations (Whole system and turbine alone); h/D= 1.47, Cpr/Tv =
0.2, C/T = 0.313, (Aerated mode)
210
Figure (5.15): The relation between impellers rotation speed and the consumed power
per water volume for the aerated mode, h/D= 1.47, Cpr/Tv = 0.2, C/Tv =
0.313
210
- 254 -
List of Figures
Appendix I
Figure (I.1):
Figure (I.2):
Figure (I.3):
Figure (I.4):
Figure (I.5):
Figure (I.6):
Figure (I.7):
225
230
232
Appendix II
Figure (II.1): The energy balance points for the spray discharge flowrate
- 255 -
238
- 256 -
List of Tables
List of Tables
Chapter One
Table (1.1):
Table (1.2):
Chapter Two
Table (2.1):
53
Chapter Three
Table (3.1):
The SAEb and SOTRb for different water level ratios (h/D) and turbine
blade width to blades submergence ratio (S/W), D/T=0.24, C/T=0.31,
Cpr/T=0.2, N=2.803 rps, (Whole system Configuration)
103
Table (3.2):
The SAEb and SOTRb for different impellers speed, D/T= 0.24, C/T=
0.31, Cpr/T= 0.2, h/D = 1.47 (Whole system Configuration, Aeration
mode)
104
Table (3.3):
The SAEb and SOTRb for different impellers speed, D/Tv= 0.24, C/Tv=
0.31, Cpr/Tv= 0.2, h/D = 1.47 (Turbine and propeller
configuration)
104
Table (3.4):
The SAEb and SOTRb for different impellers speed, D/Tv= 0.24, C/Tv=
0.31, h/D = 1.47, (Turbine alone configuration)
105
Table (3.5):
The effect the spacing between impellers on SOTRb and SAE b, D/Tv=
0.24, C/Tv= 0.31, h/D = 1.47, Sp=8.3 mm
105
Table (3.6):
The spray droplets mass transfer zone of klad, flight time t f and the
related measured parameters, (of 100% submerged turbine blades, H =
zero), h/D = 1.47, S/W = 1
116
Table (3.7):
Table (3.8):
- 257 -
List of Tables
Table (3.9):
The oxygen mass transfer rates for spray and bulk zones OTR b and
OTRsp with the percentage of contribution
120
Table (3.10): The water spray velocity and volumetric flow rate for the whole
system geometrical configuration with various turbine blades
submergence and water levels, N= 2.5 rps
122
Table (3.11): The calculated water spray mass transfer zone (klad), droplets flight
time(tf) and some related measured parameters and turbine blades
submergence and water level, at constant rotation speed (N =
2.5 rps)
126
Table (3.12): Water spray zone oxygen transfer rate (OTRsp) variation with different
turbine blades submergence and water levels, (N = 2.5 rps)
128
Table (3.13): Oxygen mass transfer rate and percentage of contribution for both
water spray and bulk (re-aeration) zones, (N= 2.5 rps)
130
Table (3.14): The spray mass transfer zone (klad), flight time (t f) and some related
measured parameters, S/W=1.00, h/D =1.47, (Turbine alone
configuration)
134
Table (3.15): The spray velocity and volumetric flow rate for various turbine rotation
speeds (Turbine alone configuration)
134
Table (3.16): The Spray zone oxygen transfer rate OTRsp for different rotation speed
levels, (Turbine alone configuration)
136
Table (3.17): The Comparison for the spray zone oxygen transfer rates OTR sp
between the whole system and turbine alone configurations, S/W
=1.00, h/D=1.47
137
Chapter Four
Table (4.1):
The RTP propeller and several pithed blade impellers power numbers
with different configurations and pumping modes
169
Table (4.2):
The RTP propeller and several pithed blade impellers pumping and
circulation numbers with different configurations and pumping
modes
170
Table (4.3):
The agitation Index for RTP-D and various pithed blade impellers,
(this work and literature)
172
Chapter Five
Table (5.1):
199
References
References
-AAdrian, R., 1991. Particle-imaging techniques for experimental fluid mechanics.
Annual Review of Fluid Mechanics 23, 261-304.
Aeschbach, S., Bourne, J.R., 1972. The attainment of homogeneous suspension in a
continuous stirred tank. The Chemical Engineering Journal 4, 234-242.
Ag, J.C., Jimenez, J., 1987. On the performance of particle tracking. Journal of
Fluid Mechanics 185, 447-468.
Aeration
AIRE-O2
www.aireo2.com.
Industries
International,
Commercial
Document,
References
Aubin, J., 2001. Mixing Capabilities of Down- and Up-Pumping Axial Flow
Impellers in Single Phase and Gas-Liquid Systems: Experimental and CFD Studies.
PhD Thesis, University of Sydney, Institut National Polytechnique de Toulouse,
France.
Aubin, J., Le Sauze, N., Bertrand, J., Fletcher, D., Xuereb, C., 2004. PIV
Measurements of Flow In An Aerated Tank Stirred By A Down-and An Up-Pumping
Axial Flow Impeller. Exp. Therm. & Fluid Sci. 28, 447-456.
Aubin, J., Mavros, P., Fletcher, D.F., Bertrand, J., Xuereb, C., 2001. Effect of Axial
Agitator Configuration (Up-Pumping, Down-Pumping, Reverse Rotation) on Flow
Patterns Generated in Stirred Vessels. Chemical Engineering Research and Design 79,
845-856.
-BBackhurst, J.R., Harker, J.H., Kaul, S.N., 1988. The performance of pilot and fullscale vertical shaft aerators. Water Research 22, 1239-1243.
Baeza, J.A., Gabriel, D., Lafuente, J., 2002. Improving the nitrogen removal
efficiency of an A2/O based WWTP by using an on-line Knowledge Based Expert
System. Water Research 36, 2109-2123.
Bandaiphet, C., Prasertsan, P., 2006. Effect of aeration and agitation rates and scaleup on oxygen transfer, kla in exopolysaccharide production from enterobacter
cloacae. Carbohydrate Polymers 66, 216-228.
Barigou, M., Greaves, M., 1992. Bubble-size distributions in a mechanically agitated
gasliquid contactor. Chemical Engineering Science 47, 2009-2025.
Barth, E.F., Stensel, H.D., 1981. International nutrient control technology for
municipal effluents. Journal (Water Pollution Control Federation), 1691-1701.
Baylar, A., Bagatur, T., 2000a. Aeration performance of weirs. WATER SAPRETORIA- 26, 521-526.
Baylar, A., Bagatur, T., 2000b. Study of aeration efficiency at weirs. Turk. J. Eng.
Environ. Sci. 24, 255-264.
Baylar, A., Bagatur, T., 2006. Experimental studies on air entrainment and oxygen
content downstream of sharpcrested weirs. Water and Environment Journal 20, 210216.
- 260 -
References
Baylar, A., Bagatur, T., Tuna, A., 2001. Aeration performance of trangular notch
weirs at recirculating system. Water Qual. Res. J. Can. 36, 121-132.
Baylar, A., Emiroglu, M.E., Bagatur, T., 2006. An experimental investigation of
aeration performance in stepped spillways. Water and Environment Journal 20, 35-42.
Baylar, A., Unsal, M., Ozkan, F., 2010. Hydraulic structures in water aeration
processes. Water, Air, & Soil Pollution 210, 87-100.
Bhattacharya, S., Hebert, D., Kresta, S.M., 2007. Air Entrainment in Baffled Stirred
Tanks. Chemical Engineering Research and Design 85, 654-664.
Bi, A.K., 1993. Gas entrainment by plunging liquid jets. Chemical Engineering
Science 48, 3585-3630.
Biotrade Co., Commercial Document, www.biotrade.fr.
Botjheng Water Ltd., Commercial Document, www.botjheng.com.
Bouaifi, M., Roustan, M., 2001. Power consumption, mixing time and
homogenisation energy in dual-impeller agitated gasliquid reactors. Chemical
Engineering and Processing: Process Intensification 40, 87-95.
Boyd, C.E., 1998. Pond water aeration systems. Aquacultural Engineering 18, 9-40.
Boyd, C.E., Martinson, D.J., 1984. Evaluation of propeller-aspirator-pump aerators.
Aquaculture 36, 283-292.
Brdjanovic, D., Slamet, A., Van Loosdrecht, M.C.M., Hooijmans, C.M., Alaerts, G.J.,
Heijnen, J.J., 1998. Impact of Excessive Aeration on Biological Phosphorus Removal
from Wastewater. Water Research 32, 200-208.
Brown, D., Jones, P., Middleton, J., Papadopoulos, G., Arik, E., 2004. Experimental
methods Handbook of Industrial Mixing Science and Practice. John Weily & Son Inc.
Bro, J., Fot, I., Sperling, R., Jambere, S., Heiser, M., Rieger, F., 2004. Pumping
Capacity of Pitched Blade Impellers in a Tall Vessel with a Draught Tube. Acta
Polytechnica 44, 48-53.
-CCancino, B., 2004a. Design of high efficiency surface aerators: Part 2. Rating of
surface aerator rotors. Aquacultural Engineering 31, 99-115.
- 261 -
References
Cancino, B., 2004b. Design of high efficiency surface aerators: Part 3. Dimensional
analysis of rotor performance. Aquacultural Engineering 31, 117-121.
Cancino, B., Roth, P., Reu, M., 2004. Design of high efficiency surface aerators:
Part 1. Development of new rotors for surface aerators. Aquacultural Engineering 31,
83-98.
Capela, S., Gillot, S., Hduit, A., 2004. Comparison of oxygen-transfer measurement
methods under process condition. Water Env. Res. 76, 183-188.
Caprari S. p. A., Commercial Document, www.caprari.com.
Carrera, J., Baeza, J., Vicent, T., Lafuente, J., 2003. Biological nitrogen removal of
high-strength ammonium industrial wastewater with two-sludge system. Water
Research 37, 4211-4221.
Chandrasekharan, K., Calderbank, P.H., 1981. Further observations on the scale-up of
aerated mixing vessels. Chemical Engineering Science 36, 818-823.
Chanson, H., 2003. Interaction of Strong Turbulence with Free Surfaces, in: Liu, P.
(Ed.), Advances in Coastal and Ocean Engineering Interaction of Strong Turbulence
with Free Surfaces. World Scientific Publishing, Singapore.
Chapple, D., Kresta, S.M., Wall, A., Afacan, A., 2002. The Effect of Impeller and
Tank Geometry on Power Number for a Pitched Blade Turbine. Chemical
Engineering Research and Design 80, 364-372.
Chern, J.-M., Chou, S.-R., 1999. Volatile Organic Compound Emission Rates from
Mechanical Surface Aerators: Mass-Transfer Modeling. Ind. Eng. Chem. Res. 38,
3176-3185.
Chern, J.-M., Yang, S.-P., 2003. Oxygen transfer rate in a coarse-bubble diffused
aeration system. Ind. Eng. Chem. Res. 42, 6653-6660.
Chern, J.-M., Yang, S.-P., 2004. Measuring and Modeling of Oxygen Transfer Rate in
a Drop Structure. Ind. Eng. Chem. Res. 43, 7657-7663.
Chisti, Y., Jauregui-Haza, U.J., 2002. Oxygen transfer and mixing in mechanically
agitated airlift bioreactors. Biochemical Engineering Journal 10, 143-153.
Chung, K., Barigou, M., Simmons, J., 2007. Reconstruction of 3-D flow field inside
miniature stirred vessel using a 2-D PIV technique. Trans IChemE, Chem. Eng,
Res.&Des. 85, 560-567.
- 262 -
References
Clesceri, N., 2009. Ch. 3 Treatment by subsurface application, Biological treatment
process, 1st ed. Humana press, NY- USA.
Costes, J., Couderc, J.P., 1988. Study by laser Doppler anemometry of the turbulent
flow induced by a Rushton turbine in a stirred tank: Influence of the size of the
unitsI. Mean flow and turbulence. Chemical Engineering Science 43, 2751-2764.
Crites, R., Middlebrooks, J., Reed, S., 2006. Natural wastewater treatment. CRC press
Taylor& Francis Group, FL-USA.
Cui, Y.Q., van der Lans, R.G.J.M., Luyben, K.C.A.M., 1996. Local power uptake in
gas-liquid systems with single and multiple rushton turbines. Chemical Engineering
Science 51, 2631-2636.
Cumby, T.R., 1987a. A review of slurry aeration 1. Factors affecting oxygen transfer.
Journal of Agricultural Engineering Research 36, 141-156.
Cumby, T.R., 1987b. A review of slurry aeration 3. Performance of aerators. Journal
of Agricultural Engineering Research 36, 175-206.
-DDehkordi, A., Savari, C., 2011. Determination of interfacial area and overall
volumetric mass-transfer coefficient in a novel type of two impinging streams reactor
by chemical method. Ind. Eng. Chem. Res. 50, 6426-6435.
Demoyer, C., Schierholz, E., Gulliver, J., Wilhelms, S., 2003. Impact of bubble and
free surface oxygen transfer on diffused aeration systems. Water Research 37, 18901904.
Deshmukh, N.A., Joshi, J.B., 2006. Surface Aerators: Power Number, Mass Transfer
Coefficient, Gas Hold up Profiles and Flow Patterns. Chemical Engineering Research
and Design 84, 977-992.
Duchene, P., Cotteux, ., 2002. Insufflation d'air fines bulles: application aux stations
d'puration en boues actives des petites collectivits (Document technique FNDAE
N 26).
Dudley, J., 1995. Process testing of aerators in oxidation ditches. Water Research 29,
2217-2219.
- 263 -
References
-EEckenfelder, W., Davis, J., Ford, L., Englande, A., 2009. Industrial water quality, 4th
ed. McGraw-Hill USA.
Edwards, M., Baker, M., 2001 A Review of Liquid Mixing Equipment, in: Nienow,
A.W., Harnby, N., Edwards, M. (Eds.), Mixing in the Process Industries, Second ed.
Butterworth-Heinemann.
El-Temtamy, S., Khalil, S., Nour-El-Din, A., Gaber A., 1984. Oxygen mass transfer
in a bubble column bioreactor containing lysed yeast suspension. Appl. Microbiol.
Biotechnol. 19, 376-381.
Engineers, A.S.o.C., Federation, W.P.C., 1988. Aeration: A Wastewater Treatment
Process. Water Pollution Control Federation.
Escudie, R., Line, A., 2003. Experimental analysis of hydrodynamics in a radially
agitated tank. AIChE Journal 49, 585-603.
-FFan, L., Xu, N., Wang, Z., Shi, H., 2010. PDA experiments and CFD simulation of a
lab-scale oxidation ditch with surface aerators. Chemical Engineering Research and
Design 88, 23-33.
FR Patent Demand, 2012. FR 2971951, Apparaeil De Traitement Mechanique D'une
Composion Liquide er Procede de Tratement Mechanique D'une Telle Composiotion
Liquide, , France.
Fuchs, R., Ryu, D.D.Y., Humphrey, A.E., 1971. Effect of Surface Aeration on ScaleUp Procedures for Fermentation Processes. Industrial & Engineering Chemistry
Process Design and Development 10, 190-196.
-GGabriele, A., Tsoligkas, A.N., Kings, I.N., Simmons, M.J.H., 2011. Use of PIV to
measure turbulence modulation in a high throughput stirred vessel with the addition of
high Stokes number particles for both up- and down-pumping configurations.
Chemical Engineering Science 66, 5862-5874.
Garcia-Cortes, D., Xuereb, C., Taillandier, P., Jauregui-Haza, U.J., Bertrand, J., 2006.
Flow Induced by Dual-Turbine of Different Diameters in a Gas-Liquid Agitation
- 264 -
References
System: the Agitation and Turbulence Indices. International Journal of Chemical
Reactor Engineering 4.
Gary, N.F., 2004. Biology of Wastewater Treatment, 1st ed. ICP press, London.
Gemgate GmbH, Commercial Document, www.gemgate.ie.
Gresch, M., Armbruster, M., Braun, D., Gujer, W., 2011. Effects of aeration patterns
on the flow field in wastewater aeration tanks. Water Research 45, 810-818.
Guillard, F., Trgrdh, C., 2003. Mixing in industrial Rushton turbine-agitated
reactors under aerated conditions. Chemical Engineering and Processing: Process
Intensification 42, 373-386.
Gulliver, J.S., Wilhelms, S.C., Parkhill, K.L., 1998. Predictive capabilities in oxygen
transfer at hydraulic structures. Journal of hydraulic engineering 124, 664-671.
Gupta, K., Sharma, R., 1996. Biological oxidation of high strength nitrogenous
wastewater. Water Research 30, 593-600.
-HHadjiev, D., Sabiri, N.E., Zanati, A., 2006. Mixing time in bioreactors under aerated
conditions. Biochemical Engineering Journal 27, 323-330.
Hecht, E., Zajac, A., 2001. Optics. Addison-Wesley Pub, USA.
Heduit, A., Racault, Y., 1983a. Essais d'aerateurs: Enseignements tires de 500 essais
en eau claire effectues dans 200 stations d'epuration differentesI. Methodologie.
Water Research 17, 97-103.
Heduit, A., Racault, Y., 1983b. Essais d'aerateurs: Enseignements tires de 500 essais
en eau claire effectues dans 200 stations d'epuration differentesII Resultats. Water
Research 17, 289-297.
Heim, A., Krasa wski, A., Rzyski, E., Stelmach, J., 1995. Aeration of bioreactors by
self-aspirating impellers. The Chemical Engineering Journal and the Biochemical
Engineering Journal 58, 59-63.
Holland, F., Bragg, R., 1995. Fluid Flow for Chemical Engineers, 2nd ed. Edward
Arnold, Great Britain.
- 265 -
References
Houcine, I., Plasari, E., David, R., 2000. Effects of the Stirred Tank's Design on
Power Consumption and Mixing Time in Liquid Phase. Chemical Engineering &
Technology 23, 605-613.
Hsu, Y., Huang, K., 1997. Effects of geometrical factors on liquid mixing in a gasinduced agitated tank. J. Chem. Tech. Biotechnol. 68, 222-228.
Huang, W., Wu, C., Xia, W., 2009. Oxygen transfer in high-speed surface aeration
tank for wastewater treatment: Full-scale test and numerical modeling. Journal of
Environmental Engineering 135, 684-691.
-JJakobson, H., 2008. Agitation and Fluid Mixing Technology, Chemical Reactor
Modeling Multiphase Reactive Flows, 1 ed. Springer-Verlag Heidelberg.
Jaworski, Z., Nienow, A., Dyster, K., 1996. An LDA study of the turbulent flow field
in a baffled vessel agitated by an axial, downpumping hydrofoil impeller. The
Canadian Journal of Chemical Engineering 74, 3-15.
Ju, L.-K., Sundarajan, A., 1992. Model analysis of biological oxygen transfer
enhancement in surface-aerated bioreactors. Biotechnology and Bioengineering 40,
1343-1352.
-KKang, J., Lee, C.H., Haam, S., Koo, K.K., Kim, W.S., 2001. Studies on the Overall
Oxygen Transfer Rate and Mixing Time in Pilot-Scale Surface Aeration Vessel.
Environmental Technology 22, 1055-1068.
Kawase, Y., Moo-Young, M., 1989. Mixing time in bioreactors. Journal of Chemical
Technology & Biotechnology 44, 63-75.
Kawecki, W., Reith, T., Van Heuven, J., Beek, W., 1967. Bubble size distribution in
the impeller region of a stirred vessel. Chemical Engineering Science 22, 1519-1523.
Keane, R., Adriane, R., 1992. Theory of Cross-Correlation Analysis of PIV Images.
App. Sci. Res. 49, 191-215.
Khopkar, A.R., Mavros, P., Ranade, V.V., Bertrand, J., 2004. Simulation of Flow
Generated by an Axial-Flow Impeller: Batch and Continuous Operation. Chemical
Engineering Research and Design 82, 737-751.
- 266 -
References
Khudenko, B.M., 1983. Development of self-propelled aerators. Journal of
Environmental Engineering 109, 868-885.
Kim, J., Walters, R.W., 2001. Oxygen transfer at low drop weirs. Journal of
Environmental Engineering 127, 604-610.
Kirke, B., El Gezawy, A., 1997. Design and model tests for an efficient mechanical
circulator/aerator for lakes and reservoirs. Water Research 31, 1283-1290.
Kleinheinz, G., Wright, P., 2009. Biological odor and VOC control process, in: Wang,
L., Pereira, N., Hung, Y.-T. (Eds.), Biological Treatment Processes. Humana Press,
pp. 733-757.
Kucukali, S., Cokgor, S., 2009. Energy concept for predicting hydraulic jump aeration
efficiency. Journal of Environmental Engineering 135, 105-107.
Kumar, A., Moulick, S., Mal, B.C., 2010a. Performance evaluation of propelleraspirator-pump aerator. Aquacultural Engineering 42, 70-74.
Kumar, B., Patel, A., Rao, A., 2010b. Shape effect on optimal geometric conditions in
surface aeration systems. Korean J. Chem. Eng. 27, 159-162.
Kumaresan, T., Joshi, J.B., 2006. Effect of impeller design on the flow pattern and
mixing in stirred tanks. Chemical Engineering Journal 115, 173-193.
Kumaresan, T., Nere, N.K., Joshi, J.B., 2005. Effect of Internals on the Flow Pattern
and Mixing in Stirred Tanks. Ind. Eng. Chem. Res. 44, 9951-9961.
Kusabiraki, D., Murota, M., Ohno, S., Yamagiwa, K., Yasuda, M., Ohkawa, A., 1990.
Gas Entrainment Rate and Flow Pattern in a Plunging Liquid Jet Aeration System
Using Inclined Nozzles. Journal of Chemical Engineering of Japan 23, 704-710.
-LLaakkonen, M., Honkanen, M., Saarenrinne, P., Aittamaa, J., 2005. Local bubble size
distributions, gasliquid interfacial areas and gas holdups in a stirred vessel with
particle image velocimetry. Chemical Engineering Journal 109, 37-47.
Lee, K., Yianneskis, M., 1998. Turbulence properties of the impeller stream of a
Rushton turbine. AIChE Journal 44, 13-24.
- 267 -
References
Lee, M., Kang, J., Lee, C.H., Haam, S., Hun-Hwee, P., Kim, W.S., 2001. Oxygen
Transfer Characteristics in a Pilot Scale Surface Aeration Vessel with Simcar Aerator.
Environmental Technology 22, 57-68.
Leng, D., Katti, S., Obeng, V., 2008. Industrial Mixing Technology, in: Albright, L.
(Ed.), Albrights Chemical Engineering Handbook, 1 ed. CRC press, NW.
Li, X., Yu, G., Yang, C., Mao, Z.-S., 2009. Experimental Study on Surface Aerators
Stirred by Triple Impellers. Ind. Eng. Chem. Res. 48, 8752-8756.
Li, Z., Bao, Y., Gao, Z., 2011. PIV experiments and large eddy simulation of singleloop flow fields in Rushton turbine stirred tanks. Chem. Eng. Sci. 66, 1219-1231.
Lines, P.C., 2000. Gas-Liquid Mass Transfer Using Surface-Aeration in Stirred
Vessels, with Dual Impellers. Chemical Engineering Research and Design 78, 342347.
Lu, W.-M., Ju, S.-J., 1989. Cavity configuration, flooding and pumping capacity of
disc-type turbines in aerated stirred tanks. Chemical Engineering Science 44, 333342.
-MMatsuura, K., Kataoka, I., Serizawa, A., 2003. Prediction of Turbulent Diffusion
Coefficient of Liquid Droplet Using Lagrangian Simulation. GENES4/ANP2003,
Kyoto, Japan.
Mavros, P., 2001. Flow visualization in stirred vessels A review of experimental
techniques. Trans IChemE 97, 113- 127.
Mavros, P., Baudou, C., 1997. Quantification of the Performance of Agitator in
Stirred Vessele: Definition and Use of an Agitation Index. Trans IChemE 75.
Mavros, P., Naude, I., Xuereb, C., Bertrand, J., 1997. Laser Doppler Velocimetry in
agitated vessels: effect of continuous liquid stream on flow patterns. IChemE 75, 763776.
Mavros, P., Xuereb, C., Bertrand, J., 1998. Determination of 3-D flow field in
agitated vessels by Laser Doppler Velocimetry: use and interpretation of RMS
velocities. Trans. Industrial Chemistry 76, 223-233.
Mavros, P., Xuereb, C., Fot, I., Bertrand, J., 2002. Investigation by laser Doppler
velocimetry of the effects of liquid flow rates and feed positions on the flow patterns
- 268 -
References
induced in a stirred tank by an axial-flow impeller. Chemical Engineering Science 57,
3939-3952.
McCabe, W., Smith, J., Harriott, P., 1985. Agitation and Mixing of Liquid, Unit
Operation of Chemical Engineering, 4 ed. McGraw-Hill, NY.
McWhirter, J., Hutter, J., 1989. Improved oxygen mass transfer modeling for
diffused/subsurface aeration systems. AIChE Journal 35, 1527-1534.
McWhirter, J.R., Chern, J.-M., Hutter, J.C., 1995. Oxygen Mass Transfer
Fundamentals of Surface Aerators. Ind. Eng. Chem. Res. 34, 2644-2654.
Merchuk, J.C., Yona, S., Siegel, M.H., Zvi, A.B., 1990. On the first-order
approximation to the response of dissolved oxygen electrodes for dynamic KLa
estimation. Biotechnology and Bioengineering 35, 1161-1163.
Mishra, V.P., Dyster, K.N., Nienow, A.W., McKemmie, J., Jaworski, Z., 1998. A
study of an up- and a down-pumping wide blade hydrofoil impeller: Part I. LDA
measurements. The Canadian Journal of Chemical Engineering 76, 577-588.
Moulick, S., Mal, B., 2009. Performance evaluation of double-hub paddle wheel
aerator. Journal of Environmental Engineering 135, 562-566.
Moulick, S., Mal, B.C., Bandyopadhyay, S., 2002. Prediction of aeration performance
of paddle wheel aerators. Aquacultural Engineering 25, 217-237.
Mueller, J., Boyle, W., Lightfoot, E., 1967. Effect of the response time of a dissolved
oxygen probe on the oxygen uptake rate. Applied microbiology 15, 674.
Mueller, J., Boyle, W.C., Popel, H.J., 2002. Aeration: principles and practice. CRC.
Myers, K., Ward, R., Bakker, A., 1997. A digital particle image velocimetry
investigation of flow field instabilities of axial-flow impellers. Journal of fluids
engineering 119, 623.
-NNagata, S., 1975. Mixing Principles and Applications. John Weily & Sons, NY.
Nair, R., Dhamole, P., Lele, S.S., DSouza, S., 2008. Biotreatment of High Strength
Nitrate Waste Using Immobilized Preadapted Sludge. Appl Biochem Biotechnol 151,
193-200.
- 269 -
References
Nakanoh, M., Yoshida, F., 1980. Gas Absorption by Newtonian and Non-Newtonian
Liquids in a Bubble Column. Industrial & Engineering Chemistry Process Design and
Development 19, 190-195.
Nakasone, H., 1987. Study of aeration at weirs and cascades. Journal of
Environmental Engineering 113, 64-81.
Nere, N.K., Patwardhan, A.W., Joshi, J.B., 2003. Liquid-Phase Mixing in Stirred
Vessels: Turbulent Flow Regime. Ind. Eng. Chem. Res. 42, 2661-2698.
Nienow, A.W., 1997. On impeller circulation and mixing effectiveness in the
turbulent flow regime. Chemical Engineering Science 52, 2557-2565.
-OOgnean, T., 1993a. Aspects concerning scale-up criteria for surface aerators. Water
Research 27, 477-484.
Ognean, T., 1993b. Dimensionless criteria for estimating oxygen transfer in aeration
systems. Biotechnology and Bioengineering 41, 1014-1020.
Ognean, T., 1997. Relationship between oxygen mass transfer rate and power
consumption by vertical shaft aerators. Water Research 31, 1325-1332.
Ohkawa, A., Kusabiraki, D., Kawai, Y., Sakai, N., Endoh, K., 1986. Some flow
characteristics of a vertical liquid jet system having downcomers. Chemical
Engineering Science 41, 2347-2361.
Okawa, T., Shiraishi, T., Mori, T., 2008. Effect of impingement angle on the outcome
of single water drop impact onto a plane water surface. Exp Fluids 44, 331-339.
Oldshue, J., 1983. Fluid Mixing Technology. McGraw-Hill, NY.
Oliveira, M.E.C., Franca, A.S., 1998. Simulation of oxygen mass transfer in aeration
systems. International Communications in Heat and Mass Transfer 25, 853-862.
Ozkan, F., Baylar, A., Tugal, M., 2006. The performance of two phase flow systems
in pond aeration. Int. J. Sci. Technol 1, 65-74.
- 270 -
References
-PPark, K., Inamori, Y., Mizouchi, M., Ahn, K., 2000. Emission and control of nitrous
oxide from a biological wastewater treatment system with intermittent aeration. Water
Research 90, 247-252.
Patil, S.S., Deshmukh, N.A., Joshi, J.B., 2004. Mass-Transfer Characteristics of
Surface Aerators and Gas-Inducing Impellers. Ind. Eng. Chem. Res. 43, 2765-2774.
Patwardhan, A.W., Joshi, J.B., 1998. Design of stirred vessels with gas entrained from
free liquid surface. The Canadian Journal of Chemical Engineering 76, 339-364.
Patwardhan, A.W., Joshi, J.B., 1999. Relation between Flow Pattern and Blending in
Stirred Tanks. Ind. Eng. Chem. Res. 38, 3131-3143.
Philichi, T.L., Stenstrom, M.K., 1989. Effects of dissolved oxygen probe lag on
oxygen transfer parameter estimation. Journal (Water Pollution Control Federation),
83-86.
Polprasert, C., Raghunandana, H., 1985. Wastewater treatment in Deep Aeration
Tank. Water Research 19, 257-264.
Praxair Technology, Commercial Document, www.praxair.com.
Process Engineering s.r.l, Commercial Document, www.pe-process.it.
Puthli, M.S., Rathod, V.K., Pandit, A.B., 2005. Gasliquid mass transfer studies with
triple impeller system on a laboratory scale bioreactor. Biochemical Engineering
Journal 23, 25-30.
-QQu, X.L., Khezzar, L., Danciu, D., Labois, M., Lakehal, D., 2011. Characterization of
plunging liquid jets: A combined experimental and numerical investigation.
International Journal of Multiphase Flow 37, 722-731.
-RRaffel, M., Willert, C., Wereley, S., Kompenhans, J., 2007. Particle Image
Velocimetry A Practical Guide, 2nd ed. Springer-Verlag.
- 271 -
References
Ramalho, R.S., 1977. Introduction to wastewater treatment processes. Academic
Press.
Ranade, V., Joshi, J., 1989. Flow generated by pitched blade turbines I: measurements
using laser Doppler anemometer. Chemical Engineering Communications 81, 197224.
Rao, A.R., 1999. Prediction of reaeration rates in square, stirred tanks. Journal of
Environmental Engineering 125, 215-223.
Rao, A.R., Kumar, B., 2007a. Scale-up criteria of square tank surface aerator.
Biotechnology and Bioengineering 96, 464-470.
Rao, A.R., Kumar, B., 2007b. The use of circular surface aerators in wastewater
treatment tanks. Journal of Chemical Technology & Biotechnology 82, 101-107.
Rao, A.R., Kumar, B., 2009. Simulating Surface Aeration Systems at Different Scale
of Mixing Time. Chinese Journal of Chemical Engineering 17, 355-358.
Rao, A.R., Kumar, B., Patel, A.K., 2009. Vortex behaviour of an unbaffled surface
aerator. Science Asia 35, 183-188.
Rodgers, T.L., Gangolf, L., Vannier, C., Parriaud, M., Cooke, M., 2011. Mixing times
for process vessels with aspect ratios greater than one. Chemical Engineering Science
66, 2935-2944.
Roubaty, J.-L., Boeglin, J.-C., 2007. Pollution Industrielle de leau Caractrisation,
Classification, Mesure, Techniques de Lingnieur.
Roustan, M., 2003. Transferts gaz-liquide dans les procds de traitement des eaux et
des effluents gazeux. Paris : d. Tec & doc Paris.
Roustan, M., 2005. Agitation Mlange Caratristiques Des Mobiles DAgitation,
Techniques De Lingnieur
Roustan, M., Ghauri, B., Roques, H., 1984. Etude de la circulation dans les bassins
d'aeration a boues activees. Water Research 18, 495-499.
Roustan, M., Pharamond, J.C., Roques, H., 1975. Etude comparative de divers
systemes d'aeration. Water Research 9, 1065-1068.
- 272 -
References
-SSano, Y., Usui, H., 1985. Interrelations among mixing time, power number and
discharge flow rate number in baffled mixing vessels. Journal of Chemical
Engineering of Japan 18, 47-52.
Sardeing, R., 2002. Hydrodynamique induite par un oxygnateur de surface :
influence du design des agitateurs sur la qualit de la dispersion, INP Toulouse.
Sardeing, R., Ferrand, F., Poux, M., Avrillier, P., Xuereb, C., 2003. Hydrodynamics
and Gas Dispersion Characterization in a System Equipped with a New Gas-Inducing
Impeller. Engineering in Life Sciences 3, 31-37.
Sardeing, R., Poux, M., Melen, S., Avrillier, P., Xuereb, C., 2005a. Aeration of Large
Size Tanks by a Surface Agitator. Chemical Engineering & Technology 28, 587-595.
Sardeing, R., Poux, M., Xuereb, C., 2005b. Procd doxygnation et de brassage
pour le traitement biologique des eaux uses, Techniques de Lingnieur, pp. IN 1-10
Schultz, T.E., 2005. Biological Wastewater Treatment, Chemical Engineering
Magazine, pp. 1-5.
Shammas, N., Chang, C., Wang, L., 2006. Tertiary micro-screening, Advanced
physico-chemical treatment process, 1st ed. Humana press Inc., NJ-USA.
Sharma, K., 2007. Principles of mass transfer. Prentice-Hall India.
Shluter, V., Decker, W., 1992. Gas/liquid mass transfer in stirred vessels. Chem. Eng.
Sci. 47, 2357-2362.
Sincero, G., Sincero, A., 2003. Physical-chemical treatment of water & wastewater.
CRC press UK.
Sirivasan, V., Aiken, R., 1988. Mass transfer to droplets formed by the controlled
breakup of cylindrical jet physical absorption. Chem. Eng. Sci. 43, 3141-3150.
Stenstorm, M., Rosso, D., 2008. Aeration and Mixing, in: Henze, M. (Ed.), Biological
Wastewater Treatment Principles Modeling and Design. IWA Publishing, UK
Stukenberg, J.R., 1984. Physical aspects of surface aeration design. Journal (Water
Pollution Control Federation), 1014-1021.
Stukenberg, J.R., Wahbeh, V.N., McKinney, R.E., 1977. Experiences in evaluating
and specifying aeration equipment. Journal (Water Pollution Control Federation), 6682.
- 273 -
References
Sun, H., Mao, Z.-S., Yu, G., 2006. Experimental and numerical study of gas hold-up
in surface aerated stirred tanks. Chemical Engineering Science 61, 4098-4110.
Supratec Co. Ltd., Commercial Document, www.supratec.cc.
-TTakase, H., Unno, H., Akehata, T., 1982a. Power Consumption of Surface Aerator in
a Square Tank. Kagaku Kogaku Ronbunshu 8, 560-565.
Takase, H., Unno, H., Akehata, T., 1984. Oxygen transfer in a surface aeration tank
with square cross section. Int. Chem. Eng 21, 128-134.
Takase, H., Unno, H., Takashi, A., 1982b. Discharge flow rate of impeller in a surface
aeration tank with square cross section. Journal of Chemical Engineering of Japan 8,
211-213.
Taricska, J., Chen, J.P., Hung, Y.-T., Wang, L., Zou, S.-W., 2009. Surface and Spray
Aeration, in: Wang, L., Pereira, N., Hung, Y.-T. (Eds.), Biological Treatment
Processes. Humana Press, pp. 151-206.
Tarshish, M., Arviv, R., Aharoni, A., 2000. Mechanical aerators' power and aeration
control in bioreactors. Journal of Environmental Engineering 126, 382-385.
Tatterson, G., 1994. Scale-Up and Design of Industrial Mixing Process McGraw-Hill
Inc, USA.
Tatterson, G.B., 1982. The effect of draft tubes on circulation and mixing times.
Chemical Engineering Communications 19, 141-147.
Tatterson, G.B., 1991. Fluid mixing and gas dispersion in agitated tanks. McGrawHill New York.
Tatterson, G.B., Yuan, H.-H.S., Brodkey, R.S., 1980. Stereoscopic visualization of the
flows for pitched blade turbines. Chemical Engineering Science 35, 1369-1375.
Thakre, S., Bhuyar, L., Deshmukh, S., 2009. Oxidation ditch process using curved
blade rotor as aerator. Int. J. Environ. Sci. Tech 6, 113-122.
Toombes, L., Chanson, H., 2005. Airwater mass transfer on a stepped waterway.
Journal of Environmental Engineering 131, 1377-1386.
Treybal, R., 1980. Mass Transfer Operations, 3 ed. McGraw-Hill Inc.
- 274 -
References
Try, P., Price, G., 1995. Sewage and Industrial Effluents, Wastewater Treatment and
Disposal. RS.C press, Manchester-UK.
Tsui, Y.-Y., Hu, Y.-C., 2008. Mixing Flow Characteristics in a Vessel Agitated by the
Screw Impeller With a Draught Tube. Journal of fluids engineering 130.
-UUdrea, D., Bryanston-Cross, P., Driver, C., Calvert, G., 1997. The Application of PIV
(Particle Image Velocimetry) and Flow Visualisation to the Coolant Flow Through an
Automotive Engine, SPIE.
-VVan't Riet, K., 1979. Review of Measuring Methods and Results in Nonviscous GasLiquid Mass Transfer in Stirred Vessels. Industrial & Engineering Chemistry Process
Design and Development 18, 357-364.
Veljkovi, V., Skala, D., 1989. Effect of number of turbine impellers on surface
aeration in laboratory fermentor. Biotechnology and Bioengineering 34, 207-213.
Vouk, D., Gjetvaj, G., Malus, D., 2005. Wastewater Aeration Using by Plunging
Water Jet Aerators, Symposium on Water Management and Hydraulic Engineering,
Ottenstein/Austria, pp. 313-322.
-WWang, T., Yu, G., Yong, Y., Yang, C., Mao, Z.-S., 2010. Hydrodynamic
characteristics of dual-impeller configurations in a multiple-phase stirred tank.
Industrial and Engineering Chemistry Research 49, 1001.
Watson, C.C., Walters, R.W., Hogan, S.A., 1998. Aeration performance of low drop
weirs. Journal of hydraulic engineering 124, 65-71.
White, D., De Villiers, J., 1977. Rates of Induced Aeration in Agitated Vessels.
Chem. Eng. 121, 113 -118.
Wichterle, K., 1994. Heat transfer in agitated vessels. Chemical Engineering Science
49, 1480-1483.
- 275 -
References
Wiesmann, U., Choi, L., Domdrowski, E., 2007. Fundamentals of biological
wastewater treatment, 1st ed. WILEY-VCH verlag, Berline-Germany.
Wormleaton, P.R., Soufiani, E., 1998. Aeration performance of triangular planform
labyrinth weirs. Journal of Environmental Engineering 124, 709-719.
Wormleaton, P.R., Tsang, C.C., 2000. Aeration performance of rectangular planform
labyrinth weirs. Journal of Environmental Engineering 126, 456-465.
Wu, H., 1995. An issue on applications of a disk turbine for gas-liquid mass transfer.
Chemical Engineering Science 50, 2801-2811.
-XXuereb, C., Poux, M., Bertrand, J., 2006. Agitation et Mlange, Aspects
Fondamentaux et Application Industrielles. Dound, Paris.
-YYatomi, R., Takenaka, K., Takahashi, K., Tanguy, P., 2008. Mass-transfer
characteristics by surface aeration of large paddle impeller: application to a
polymerization reactor with liquid level change. Chem. Eng. Res. Des. 86, 1345-1349.
Yeh, N., Rochelle, G., 2003. Liquid-phase mass transfer in spray contactors. AIChE
Journal 49, 2363-2373.
Yoshida, F., Ikeda, A., Imakawa, S., Y., M., 1960. Oxygen absorption rates in stirred
gas-liquid contactors. Ind. Eng. Chem. 52, 435-438.
-ZZeybek, Z., Abilov, A.G., Alpbaz, M., 1997. Optimization Strategies of Gas/Liquid
Mass Transfer in Mechanically Agitated Tanks. Chemical Engineering Research and
Design 75, 480-486.
Zhu, H., Nienow, A.W., Bujalski, W., Simmons, M.J.H., 2009. Mixing studies in a
model aerated bioreactor equipped with an up- or a down-pumping Elephant Ear
agitator: Power, hold-up and aerated flow field measurements. Chemical Engineering
Research and Design 87, 307-317.
Zlokarnik, M., 1979. Scale-up of surface aerators for waste water treatment, Advances
in Biochemical Engineering, Volume 11. Springer Berlin Heidelberg, pp. 157-180.
- 276 -