Effect of Viscosity, Surface Tension and Nozzle Size On Atomization in Two-Phase Nozzles
Effect of Viscosity, Surface Tension and Nozzle Size On Atomization in Two-Phase Nozzles
Effect of Viscosity, Surface Tension and Nozzle Size On Atomization in Two-Phase Nozzles
Abstract
This study compares the mean drop size (SMD) produced from sprays of different viscosities, surface tensions
and sizes of small-scale two-phase nozzles. A drop size correlation was established as a function of nozzle size,
viscosity, surface tension and flow conditions. The small-scale nozzles are geometrically similar to commercial
scale nozzles used in fluid coker (FC) nozzles. The experiments were performed using mixtures of air with various
liquids - water, canola oil and a glycerine-water mixture. The liquid viscosities varied from 1mPa-s to 67mPa-s, and
the surface tensions varied from 25mN/m to 70mN/m. All fluid properties were at 21°C. The nozzles used were a
one-quarter and one third scale nozzle, which were geometrically similar to a full-size (or commercial) FC nozzle.
The liquid flow rates varied from 95g/s to 196 g/s, and the GLR was fixed at 1%, similar to the commercial FC
nozzles. . The mean drop size (SMD) within the spray was measured using a 2-D Phase Doppler Particle Analyzer
(PDPA). Measurements were performed at an axial distance of 100 mm from the nozzle exit and across the spray
within the range -50 mm to 50 mm in the horizontal plane. The results show that at the representative spray radius,
the SMD increased with viscosity by a maximum of 40%. A decrease in surface tension increased the SMD by a
maximum of 9%, which is not significant. Negligible difference in SMD with increase in nozzle size was observed
across the spray for the glycerine-water spray. However, for the water spray, the SMD produced at representative
spray radius increased by up to 32%. Finally, the correlation obtained using an area-averaged SMD across the spray
showed a strong relationship with liquid viscosity and nozzle exit diameter, but a weak relationship with surface
tension. The correlation gave a maximum deviation of about 17% from measured drop sizes. Results from this study
provide a comprehensive means of improving the design of two-phase nozzles, which can be used in the FC.
2
fluid atomizer studies, where the SMD increases with nozzle exit orifice. For a nozzle discharging into
nozzle exit orifice diameter. ambient surroundings of given temperature and
From the atomization studies mentioned above, it is pressure, the following can suffice:
evident that the different designs and classes of two-
phase nozzles exhibit different atomization SMD = f (D, L, ρ G , ρ L , QG , Q L , μ G , μ L , γ, g) (1)
characteristics. As such, the atomization characteristics
of FC nozzles using various liquid properties are likely In the above equation, L, g and Q denote nozzle
to behave differently from the other classes of gas-liquid length, acceleration due to gravity and volume flow rate,
nozzles. To the best of the author’s knowledge there has respectively. Subscripts “L” and “G” denote liquid and
been no study on the effects of liquid properties and gas properties, respectively. Applying the Buckingham
nozzle size on the atomization behaviour of FC nozzles. PI principle, and re-writing the gas density using the
This is the main motivation of the current study. ideal gas equation, one can obtain the following
dimensionless groups:
Specific Objectives
This study is aimed towards comparing the
SMD ⎛ L μG mG PG ρ Q Q2 ρ Q2 ⎞
differences in radial profiles of SMD in sprays produced = f ⎜⎜ , , , , L L , L5 , L 3L ⎟⎟ (2)
from two-phase nozzles. All comparisons are made at D ⎝ D μL mL RGTabsρL μL D gD γD ⎠
an axial position near the nozzle exit i.e. ~ 100 mm. The
flow conditions studied have a gas-to-liquid ratio by In Eqn. (2), PG is the absolute pressure of the gas
mass (GLR) of 1%, which is similar to those in the FC. phase at the nozzle just before the exit, which is also the
The first and second objectives of this study are to pressure of the liquid phase at the nozzle exit. The
quantify the variations in SMD for different liquid symbols RG and Tabs denote the characteristic gas
viscosities and surface tensions, respectively, in sprays constant and absolute temperature of the gas,
produced from a typical one-quarter scale FC nozzle. respectively. Note that the third term in the brackets is
The third objective is to quantify the variation of SMD the GLR, whereas the last three terms are the Reynolds
for two nozzle exit orifice sizes. The final objective is to number (Re), Froude number (Fr), and Weber number
establish an empirical relationship for an area-averaged (We) based on the liquid superficial conditions. Ideally,
SMD in the spray as a function of nozzle size, liquid to study the effect of surface tension on SMD, all terms
viscosity, surface tension and flow conditions. It is within the brackets in Eqn. (2) must be matched, but the
intended that the results from this small-scale nozzle We term must be significantly different. Similarly, to
study could be applied to the full-scale FC nozzles to study liquid viscosity effects on SMD, all terms within
better understand their operating performance and aid in the brackets in Eqn. (2) must be matched, but the Re
the design of future of nozzles. term must be significantly different. In most cases
matching the required terms may not be practically
Methods and Apparatus feasible [12]. Hence, similar to standard dimensionless
The smallest nozzle assembly used during the test analysis, restriction on some of the quantities may be
had conduit length and internal diameter of 368 and 5.2 relaxed depending on knowledge of the dominant forces
mm, respectively. The nozzle exit diameter (D) was (in this case forces affecting atomization) in the system
3.1 mm. This nozzle was geometrically one-quarter as stated by Douglas et al. [14]. In this study, the
scale of a patented full-scale design (US Patent #: nozzles used were geometrically similar hence the first
6003789) employed in a fluidized bed coker for heavy condition within the brackets in Eqn. (2) is satisfied.
oil upgrading. The second nozzle used in the test was Secondly, since for all fluids used μG<<μL, the second
manufactured to be geometrically-similar and was one- condition is also satisfied. The GLR (third term in the
third scale of the full-scale nozzle used in commercial brackets) is an important atomization parameter in gas-
applications in the fluid coker. During the test, the liquid flows (as stated earlier), hence it was kept
nozzles were mounted on a traversing rig, which was constant (GLR = 1%) in our analysis. Note that GLR is
capable of three-dimensional motion. fixed at 1% because it is the flow condition in the
To study the effects of surface tension, viscosity or commercial FC. The mixture pressure, PG is an
nozzle geometry on SMD the flow conditions were important parameter in the gas-liquid atomization of a
determined using the dynamic similarity procedure for given liquid density. Therefore the gas-liquid density
two-phase (gas-liquid) flow stated by Chesters [11]. ratio (DR), i.e. the fourth term in the brackets in Eqn.
This procedure was applied experimentally by Geraets (2) was also conserved in this study.
[12] to estimate the pressure drop along two-phase (gas- The test liquids were water, canola oil and aqueous
liquid) horizontal flows with an uncertainty of ±5%. A glycerine mixture. Liquid density, ρL, was obtained by
similar method was also applied by Hurlbert et al. [13] measuring the volume and mass of liquid samples, using
in predicting pressure drops along horizontal tubes. a graduated cylinder and electronic mass balance,
Note that the above technique is applicable to respectively. The liquid viscosity was measured using
isothermal gas-liquid flows with without mass transfer. Cannon-Fenske Routine viscometers, whereas surface
In this study SMD is the dependent variable (instead of tension, γ, was measured using the pendant-drop
pressure drop), and flow conditions are based at the technique in conjunction with the computer code
3
ADSA-P (Axisymmetric Drop Shape Analysis-Profile). circulated back to the nozzle. The mist produced during
From the image of the drop and known drop density, the the spraying process was removed from the top of the
software computes the liquid surface tension. The liquid collector tank/reservoir using a mist
properties of the test liquid are presented in Table 1. extraction/ventilation system. The system consists of a
blower and 2.2 kW (3 horsepower) motor assembly with
Table 1 – Properties of the test liquids at 21°C. Data the capacity to extract about 56.6 m3/min (2000 CFM)
with superscripts ‘∗’ were adopted from White [14]. of free air at a static pressure of 0.14 m (5.5”) water
Absolute column.
Density, Surface The SMD within the spray was measured using a
Liquid viscosity,
ρL tension, γ 2-D Phase Doppler Particle Analyzer (PDPA). The focal
μL
(kg/m3) (mN/m) lengths of the PDPA transmitter and receiver lenses
(mPa-s)
were 400 and 310 mm, respectively. The PDPA unit
Water 998* 70 ± 2 1* consists of an Nd-YAG and He-Ne laser with
Glycerine wavelengths of 532 and 632.8 nm, respectively. During
1200 ± 30 61 ± 4 67 ± 3
solution data collection, the PDPA was operated in forward-
Canola Oil 905 ± 1 25 ± 3 66 ± 1 scatter and refraction mode, and the receiver was set to
a scattering (or off-axis) angle (φ) of 30° for the air-
water tests. Forward scattering is chosen in this study
Figure 1 presents the schematic of the test set-up used in
since from laser theory scattered light from particles in
this study. Compressed dry air and filtered liquid were
this mode is about 102 orders of magnitude higher
mixed at a tee junction and fed through the nozzle
compared to the backward scattering mode [16].
assembly. Air and liquid flow rates were measured with
Furthermore, first order refraction is the most dominant
a thermal mass flow meter and variable area piston-type
flowmeter, respectively. Static pressures close to the scattering mode at φ = 30° for water drops in air [16].
mixing point of the fluids and at the nozzle exit were However, a specific scattering angle was not stated for
measured using a pressure gauge and Validyne canola oil drops or drops from the glycerine-water
transducers, respectively. The room and liquid mixture. For consistency and also based on the high
temperatures were measured using thermocouples. The confidence level in drop size measurement at this
temperatures near the mixing point of the fluids, at the scattering angle [16], the scattering angle was also set to
nozzle exit and within the spray were also measured 30° for these tests. Each flow condition consisted of
using thermocouples. The liquid flow rates varied from about 2 to 3 Runs. The typical sample size and sampling
95g/s to 195 g/s, and the GLR was fixed at 1%, similar time during data collection were set to 12000 or 360
seconds, respectively. This sample size and time were
Mixing tee PDA measurement To mist
sufficient to ensure data was collected for steady state
Duct
Nozzle positions extractor
conditions, and independent of velocity samples [16].
Spray
Profiles of SMD (and corresponding mean axial
Filtered air from velocities) were obtained within radial (or horizontal)
compressor
positions, y in the range -50mm≤ y ≤ 50mm about the
spray axis, and axial positions, x = 100 mm downstream
of the nozzle exit. Data at this axial distance was the
main focus of study since this is the region of interest
downstream of the full-scale (or commercial) nozzle in
the fluid coker.
Liquid by-pass line Liquid
tank
Results and Discussion
Liquid supply line The SMD profiles are presented in subsequent Figures
below. Comparisons of the profiles in each plot are
Legend
Flowmeter; made at the radius corresponding to the highest liquid
Liquid filter;
volume per unit radius for each spray. These radii are
Pump; Flow valve; termed representative radii in this study.
Pressure gauge; Check valve;
Effect of viscosity
The SMD profiles across the spray in the air-water
(μL = 1 mPa-s) and glycerine-water systems (μL = 67
Figure 1 - Schematic of the experimental set-up. mPa-s) at an axial distance of 100 mm downstream of
to the commercial FC nozzles. The corresponding the nozzle exit are presented in Figure 2. The
pressures at the fluid mixing section varied from 221 to corresponding liquid mass flow rates for same GLR
327 kPa, whereas the nozzle exit pressure varied from (=1%) and density ratio (~0.0032) were 95 g/s and 137
508 to 990 kPa. Downstream of the gas and liquid g/s for the water and glycerine solution systems,
mixing point the resulting spray was discharged into a respectively. The representative radii are 20 mm for the
liquid collector tank/reservoir, with the liquid re- 67 mPa-s liquid and 15 mm for the 1 mPa-s liquid.
4
250 250
200 200
SMD (microns)
SMD (microns)
150
150
100
100
50
67 mPa-s
50 61 mN/m
1 mPa-s
25 mN/m
0
-60 -40 -20 0 20 40 60 0
Horizontal distance, y (mm) -60 -40 -20 0 20 40 60
Figure 2: Profiles of SMD in two sprays with different Horizontal distance, y (mm)
liquid viscosities. Figure 3: Profiles of SMD in two sprays with different
surface tensions.
Comparing the SMD data at the representative radii for 250
both systems the glycerine solution produces an SMD
about 40% greater than that of the water system. A 5%
test of significance suggests that this difference in SMD 200
is significant. The presence of the large drops in the 67
S M D (m ic ro n s )
5
Effect of Nozzle size geometry open squares denoted ‘correlation-fit data’ are the data
Figure 4 presents the SMD profiles within the sprays used to obtain the correlation in Eqn. (3). The Figure
produced by the two nozzle sizes using air and the shows that the maximum deviation between estimated
aqueous glycerine solution. The liquid mass flow rates and actual SMDgb data is about 17%.
at 1% GLR and density ratio of 0.0022 were 95 g/s and
196 g/s for the small (D = 3.1 mm) and big (D = 4.1 225
6
Acknowledgements 14. White, F.M., Fluid Mechanics, McGraw-Hill
The authors acknowledge the Natural Sciences and Series, U.S.A, p. 771, 1999.
Engineering Research Council (NSERC) of Canada for
funding. CEE also acknowledges NSERC for an
Industrial Partnership (IP) scholarship. B. Knapper and
E. Chan are also acknowledged.
References
1. Lefebvre, A.H., Atomization and Sprays,
Hemisphere Publishing Corporation, New York,
pp. 27-260, 1989.
2. Elkotb, M.M, Mahdy, M.A., and Montaser, M.E.
Investigation of external-mixing airblast atomizers,
Proceedings of the 2nd International Conference
on Liquid Atomization and Sprays, Madison,
Wisconsin, U.S.A, pp. 107-115, 1982.
3. El-Shanawany, M.S.M.R. and Lefebvre, A.H.,
Airblast atomization: The effect of linear scale on
mean drop size, Journal of Energy, vol. 4, No. 4.
pp. 184-189, 1980.
4. Buckner, H.N. and Sojka, P.E., Effervescent
atomization of high-viscosity fluids: Part I,
Newtonian liquids, Atomization and Sprays, vol. 1,
pp. 239-252, 1991.
5. Lund, M.T., Sojka, P.E., Lefebvre, A.H., and
Gosselin, P.G., Effervescent atomization at low
mass flow rates. Part 1: the influence of surface
tension, Atomization and Sprays, vol. 3, pp. 77-89,
1993.
6. Santangelo, P.J., and Sojka, P.E. A holographic
investigation of the near nozzle structure of an
effervescent atomizer produced spray, Atomization
and Sprays, vol. 5, pp. 137-155, 1995.
7. Sovani, S.D., Sokja, P.E. and Lefebvre, A.H.,
Effervescent atomization, Progress in Energy and
Combustion Science, vol. 27, pp. 483-521, 2001.
8. Lefebvre, A.H., Wang, X.F., and Martin, C.A.,
Spray characteristics of aerated-liquid pressure
atomizers, AIAA Journal of Propulsion and Power,
vol. 4, No. 6, pp. 293-298, 1988.
9. Roesler, T.C. and Lefebvre, A.H., Studies on
aerated-liquid atomization, International Journal of
Turbo Jet Engines, vol. 6, pp.221-230, 1989.
10. Chesters, A.K., The applicability of dynamic-
similarity criteria to isothermal liquid-gas two-
phase flows without mass transfer, International
Journal of Multiphase Flow, vol. 2, pp. 191-212,
1975.
11. Geraets, J.J.M., Centrifugal scaling of isothermal
gas-liquid flow in horizontal tubes, International
Journal of Multiphase Flow, vol. 14, No. 3, pp.
287-303, 1988.
12. Hurlbert, K.M., Witte, L.C., Best, F.R, and
Kurwitz, C., Scaling two-phase flows to Mars and
Moon gravity conditions, International Journal of
Multiphase Flow, vol. 30, pp. 351-368, 2004.
13. Douglas, J.F., Gasiorek, J.M, and Swaffield, J.A.,
Fluid Mechanics, Singapore, Longman Scientific
and Technical, p. 268, 1995.