Local Studies in Horizontal Gas - Liquid Slug Flow
Local Studies in Horizontal Gas - Liquid Slug Flow
Local Studies in Horizontal Gas - Liquid Slug Flow
175
two phase flow structure, it is essential to experimentally obtain detailed local values of fundamental
parameters.
The most sienificant and essential parameters associated with the slug flow pattern are the
distribution of gas and fiquid phases, the fiquid velocity and its fluctuating components, the gas bubble
and liquid transit frequency (or slug length), and the turbulent transport characteristics of interfacial
mass, momentum, and energy. These vriables describe the local flow conditions of the quasi-steady
slug flow, both qualitatively and quantitatively.
parameters and generalized relationships among them are necessary to understand the turbulent
transport phenomena of the two-phase flow pattern.
A comprehensive physical model describing horizontal gas-liquid slug flow was first initiated
by Dukler and Hubbard [5]. This model has been modified and extended over the years by Nicholson
et al. 6], Fabre et al. 7], Bornea and Brauner [8], Andreussi and Bendiksen 9], Andreussi et al.
I 0] and Moalem Maron et al. I I , to apply to the entire intermittent flow-pattern. The predictive
models developed by these investigators make it possible to obtain average liquid velocities both in
the liquid slug and the liquid region underneath a slug bubble, pressure drops, length of liquid slug,
and slug frequencies, if the gas and liquid mass fluxes are provided. These models seem to give
reasonable results when compared to experimental data of global measurements. However, these
models cannot give the detailed void fraction distribution due to small and large bubbles, local
velocity distnibution, and the turbulent structure throughout the liquid phase. This information is of
great importance to the eventual understanding and modeling of the basic hydrodynamics of
two-phase slug flow.
A large number of experimental investigations have been carried out to develop and verify global
slug flow models in horizontal flow configurations [10, 12 - 17]. These investigations have been
176
concerned with measurements of pressure drop, overall void fraction and of statistical characteristics
I I
g slug bubbles.
such as slug length, slug frequency and liquid film thickness underneath the lar e
However, due to the experimental difficulties associated with the intermittent nature of slug flow,
very few det"ed, local measurements have been reported in the literature. The problem of obtaining
local data is further complicated in horizontal flow configurations by the facts of axial asymmetry of
the internal structure and that the slug flows do not exhibit a quasi-fully-developed equilibrium
condition.
Kvemvold et al. [ 7 used the combination of LDV and optical two-phase probes to measure
the axial velocity distribution throughout a slug flow unit in a 24 mm ID horizontal tube of
atmospheric pressure. However, the method is limited to relatively low gas velocities since the
application of LDV induces technical difficulties in regions with high concentrations of small gas
bubbles in the liquid slug occurring at higher gas flow rates. Andreussi et al. I
and cross-sectional (conductance) probes to measure the radial void fraction distribution in the liquid
slugs, the size of the dispersed bubbles in the liquid slua, and the aeration of the liquid layer
underneath the slug bubble. Kawaji et al. [ 9 used the photochromatic dye activation technique to
visualize the instantaneous motion of the liquid and gas slugs, and to successfully measure axial and
vertical velocity profiles of the liquid phase in a horizontal slug flow. The experimental data were
obtained for both circular and rectangular channels, and the liquid flow structure was found to be
quite sirnilar between the two channels. These recent studies have provided detailed basic information
on the internal structure of the intermittent flow pattern in a horizontal configuration. Information
associated with the turbulent structure and the local void fraction contributions due to elongated large
bubbles and small bubbles are not addressed in these studies.
177
In view of the above discussion, it is evident that much experimental work is still necessary to
attain a thorough physical understanding of the internal structure of an intermittent two-phase slug
flow-pattern. In this context, an experimental investigation has been underway at the University of
Wisconsin-Nlilwaukee to clarify the local velocity behavior and the turbulence structure of this flow
pattern.
In these experiments the hot-film anemornetry technique was used to measure the
time-averaged local void fractions due to small and large slug bubbles, as well as, the local axial
velocity and turbulence in the liquid phase of an air-water ntermittent flow in a 50.3 mm IID
horizontal channel. In the folloWling, the test facility is described, and based on the data, preliminary
results are documented in terms of the local axial velocity behavior within the liquid slug and the
liquid film underneath the large gas slug bubbles.
two-phase flow, and described the probe behavior in great detail under various flow patterns. With
a careful treatment of the data, he showed that local measurements of void fraction, liquid velocity
and turbulence intensity in the liquid phase could be achieved by a proper use of the hot-film
anemometrv in air-water flows. Sce then, this technique has been used by Serizawa et al. 22, 23],
Abel and Resch 24], Wang et al. 25, 26], Liu and Bankoff 27, 28], Lance and Bataille 29] and
Grossetete 30] for describing the nternal turbulence structure and phase distributions in vertical
bubbly flow patterns. The poneering studies of Theofaneous and Sullivan 3 1] demonstrated the
utility of LDV to measure the turbulence structure in bubbly two-phase flow. However, only limited
178
efforts were made to examine two-phase flow characteristics in large scale experiments of the slug
flow-pattern.
in pnciple, the hot-film anemometry method consists of the instantaneous measurement of the
change in heat transfer from an electrically heated sensor. As the fluid flows past the constant
temperature hot-film probe, changes in the fluid velocity, including turbulent fluctuations, cools the
sensor at different rates. These changes in cooling rates result in voltage changes in the anemometer.
The voltages are digitized and recorded in a PC, where they can later be converted into fluid
velocities. In the case of an air-water two phase flow, the heat transfer rate between the two fluids
is dramatically different. This results in abrupt voltage changes as the probe encounters phase
interfaces. A typical sensor output for a two-phase slug flow is illustrated in Figure I a and b.
As seen in Figure I a, when the probe resides in the upper portion of the pipe the sensor
encounters the slug bubbles. After the sharp initial drop, caused by the probe piercing the nose of
a slug bubble, the voltage gradually continues to decrease wile the sensor is inside the gas slug. This
is due to the evaporation of a thin film of liquid that remains on the sensor. When encountering a
small gas bubble, the signal shows a sharp drop followed immediately by a sharp increase. The probe
does not encounter any slug bubbles when positioned in the lower portion of the pipe (Figure lb).
However, the voltage signal shows a quasi-penodic wave motion believed to be caused by the
passage of slug bubbles over the top of the probe.
Signal Processing
Before the two-phase voltage output data can be converted to velocities, the portion of the
signal related to the gas phase must be removed. Therefore, a phase separation technique was
developed. As has been demonstrated by Lance and Bataille 29], the peaks associated with the
sensor encountering an interface can be amplified by calculating the slope of the voltage signal,
179
(de/dt), where e is the voltage level and t is the time. By careful comparison of the voltage and slope
values at the nose and tall interfaces, threshold values were set that distinguish the interface from
normal liquid or gas turbulence signals. Once the interfaces were found, a phase identifying signal
( = I
in liquid phase) was created for use in phase separation of the voltage
signal and calculation of the local time-averaged void fraction. The voltage signal corresponding to
the liquid phase was then calibrated into velocities.
4
3
go,',
2
>
C
0 0
0. 5
1 0
T im e.
i 5
2 0
s ec.
4 0
C,
0
3. 5
(b)
2 .5
0 0
0.5
I 0
2 0
Time, sec.
1 I m/s,
As was seen in Figure lb, when the probe is positioned near the bottom of the pipe, the voltage
signal shows a wave motion induced by the ntermittent nature of slug bubble passage over the probe.
This wave causes errors in the calculation of the turbulent fluctuations. A time-domain filtering
method was introduced to remove the wave.
180
of pyrex tube, 50.3 nun ID. x 153 m in length, flanged together. Between each section, a pressure
tab is installed to measure pressure drop. The flow loop is about 15.4 m in total length, and is entirely
transparent to aow for Dow visualization, high-speed photography, and high-speed cinematography.
Water
brain
Air
INCICItaAgeJOill
ANWait' 1111-i"
C1MUCIS K Poeunwic operated Bill valeS
WIlef 110.metw$ ofapplopille sus
L Molot control
m etIto. matef contiol wl"s
m -Computer and dat aq.wllun system
Ai;::o,
site
14 250 gal. Ali 3.1,
"' 1, 1,
o 10iAjes
P 500 gal. Water unk
F M II.. (IIguiall"o "'ds
0 Aii-Wallit sertaialoi. .111iInle(n3l battles
G Al,PICSSWIltoLaio,
n Wilef still[ Oil0111
11 Al.1111ti
S 20 Ito 750 opfn Wile( pump
Wile$ P'Issu's Ithel ;INCS
press.$
laps
50.3I.G.line
25.2m 1.0.Line
h71
\K
5.64
W.-,A
1.9 M3
181
Air is injected
into the water through a cylindrical porous media of I 0 tm porosity. At the other end of the flow
loop the two phase mixture enters a separator where the air is vented to the atmosphere and the water
is returned to the storage tank. In the water st6rage tank, the water temperature is maintained by a
tap water cooling systemSeven pressure taps are mounted along the flow loop. Six diaphragm type pressure transducers
alone, with six U-tube manometers are used to measure the pressure drop. The pressure transducers
have a natural frequency of
Hz. with a range of 0 to 34.4 Pa, and an accuracy of 0.3% of the full
scale. The pressure of the air at the location of the flow meter and the two-phase system pressure
measured at the test section are both measured and used to correct for the compressibility effects of
air.
Experimental A-ocedure
The experiments were carried out usina.n hot-film anemometry and conical shaped (TSI 1264
AW) rruriature platinum hot-film probes. The base of the cone is 075 m
in diameter. Frequency
response of the probe is rated as 150 kHz in air. It is expected greater in water. The liquid and gas
volumethc superficial velocities ranged from 1 IO to 220 m s-' and 027 to 220 ni s-1, respectively,
to develop plug and slug flow regimes. For all the flow conditions, the system pressure was near
atmospheric and the temperature about 20-22 'C.
within
182
As shown in Figure 4 two hot-film anemometers were used in this experiment. One probe was
used exclusively to identify the gas and liquid phase while the other probe measured the instantaneous
axial velocity components.
a dial to facilitate traversing linearly. The two probes were mounted diametrically opposite to each
other along the vertical axis. For each probe, after the probe was centered along the axis of the tube,
the position of the probe was read on a digital linear scale. The high resolution was necessary to
evaluate probe positions in the flow stream accurately and to ensure reproducible results. To begin
with, the top probe was positioned at r/R value of 0.8 and the bottom probe was positioned at r/R
value of 0.7. The top probe was held at this location throughout the experiment. The bottom probe
was traced through the vertical axis of the pipe, stopping at I positions to take measurements. The
local, instantaneous liquid velocity was measured at each location by the bottom probe and the top
probe was used to determine whether the bottom probe was in the liquid slug or in the liquid film
underneath passing gas slug.
P-b-
FCrOUS MEDIA
5Cj'
- - - - - - - - - -- - - -- - - - - -
183
For each preset experimental condition, the pressure drop along six sections of pipe was
recorded using electronic pressure transducers and verified by U-tube manometer measurements. In
addition to the liquid and gas flow rates, the liquid temperature, gas absolute pressure, and system
absolute pressure were recorded.
Before beginning the two-phase measurements, single-phase liquid measurements were made
in the same loop to calibrate the nstrumentation, verify their consistency with known results, and to
serve as a reference for later comparison with two-phase flow measurements. The local mean axial
velocity and the turbulent fluctuations in the axial direction were measured along the radial direction.
The measured velocity profiles and turbulent fluctuations were non-dimensionalized with respect to
the characteristic velocities of ceriterline velocity and friction velocity, respectively, for the purpose
of comparing with Laufer's 32] and Liu and Bankoffs 27) data. The time averaged local axial
velocity U&), and the root-mean-square values of the turbulent fluctuations, u(r), were found to
be reasonably satisfied when compared with those results provided in these references.
These
Following ts
single-phase liquid flowing at an area-averaged mean velocity, <U.,,,>, equal to the mean superficial
liquid velocity, jf>, desired in the two-phase experiment.
The data from the probe was collected by the anemometer and stored in a PC computer. Due
to the limitations of the computer, the sampling rate was set to
meaningful sampling time. Once the data are stored in the memory of the computer, a FORTRAN
184
program was used to process the data, separating the phases, converting the voltages to velocities,
and calculating the essential parameters.
The local void fraction, x, at any point, r, can be obtained by the hot-film probe sensor. It is
defined as a time-averaged void fraction by
a(r)
lim
T 6(r, t) dt
T M fo
where , as a function of the space coordinate, r, and time, t, is equal to I if the probe sensor is in te
gas phase and
is equal to
Nib
Nb
(t 2i
T j=1
2 - 1) + E (t2i jrl
(2)
2 - I)sb
where t2,, is the time when the probe sensor enters into the small bubble, and
t2i
enters into the liquid phase. Similarly, the subscript identifies large elongated slug bubbles. N,,, and
NO,
respectively, are the number of large and small bubbles passing the probe sensor in the total
(3)
N,b
(4)
and
1,
The local, mean aal velocity and the root-mean-square values of turbulent fluctuations were
calculated using
U ,,(r)
uk(rt) IN
u'(r)
(5)
1/2
[Uk(rt - U.,(r)]2/ N
fk= 1
185
(6)
respectively. In Eq. 6), Uk(rt) is the instantaneous axial velocity for the Wth data point in the liquid
phase, and N is the total number of data points in the liquid phase of the digital sample, k = 1, . . ,
N. As described above, to remove the error caused by the intermittent wave motion, the time
averaged, mean velocityin Eq. 6), U,,(r), is replaced With a curve fit, time dependent mean velocity,
UVc c(rt) as
N
u'(r)
Y:
12
luk (r, t
-U.,,:,
CF (rt)]'/N
k=1
(7)
I
From the local values of a(r) and U,Jr) measured along the vertical axis, and with the
hypothesis that the flow variables are invariant at a given horizontal slice located at -I< r/R < , the
area- averaged liquid superficial velocity jf> was calculated as follows:
AI f [I - a r) U.
dA
(8)
be several reasons for such a consistency. First, as noted by Wang et al. 251 the void fraction
measured by the conical hot-film probe technique is underestimated due to the deformation and the
deflection of the bubbles by the probe.
pronounced in measuring the small bubble contribution toward the total void fraction expressed by
Eq. 2). Secondly, as demonstrated by Figure I b, the data shows a wavy pattern when the probe is
located underneath the slug bubble. Although a correction on the time domain was performed to
186
smooth out the behavior, this may still cause some errors in calculating the mean velocity in the liquid
phase. Finally, the small-amplitude peaks due to the ncomplete percing of small bubbles or to the
bubble sliding on the probe are difficult to detect by the present data processing. This is again
expected to be the case in the small bubble encounter of the probe. These series of experimental
errors are probably the main causes of the systematic error observed in matching the liquid superficial
velocities.
Once the statistical data is extracted from the raw data, the results are then graphed for visualinvestigation of characteristics. In the graphs the instantaneous velocity values of several slug units
are plotted and a best curve is fitted to these points. The absence of points indicate either slug units
were fewer than the minimum selected during the sampling time or there were data points which
could not be considered for analysis. The time-averaged local values of the void fraction due to small
and large bubbles, the overall mean liquid velocity, and the local turbulence fluctuations were
presented elsewhere 34]. Here we are presenting preliminary results in terms of the local axial
velocity within the liquid slug and the liquid film underneath the large gas slug bubbles. Additional
data analysis is expected to lead to a more detailed understanding of two-phase slug flow.
2)
3)
Liquid layer below the gas slug, wtch is subdMded into two regions, namely
3. 1) Near the gas slug nose, and
187
O
SILIg
(D -Gas
Sug
e
3.1
Unt
]-IqLlld
S
Sug
The
first figure shows the instantaneous velocity profile at 002 seconds ahead of incoming gas slugs. The
successive frames are presented in such a way that local liquid velocity profiles in liquid slugs can be
analyzed at every 0.02-secondsiritervals In front of incoming gas slugs. At least three moving liquid
sluas are used to construct these figures. In presenting these experimental data no other statistical
averaging techniques were used. The solid lines indicate simple fourth-order curve fittings.
The gas slug interfacial velocities which were measured by the four-sensor resistivity probe
technique 35] showed that the interfacial velocity of these gas slugs ranged between 27 to 325 m/s.
The frequency of gas slugs was about 3 to 4 whereas, the average gas slug length was about 5
60 cm. From the first few figures it is evident that there exists an acceleration region very close to
the nose of gas slugs. Ts
gas slugs.
In this region the location of the maximum velocity moves from about center line
downward below the center line of the pipe. Here the liquid moves downward due to the downward
curvature of faster moving gas slug-liquid interface near the nose resulting in the appearance of the
188
......
Z'>
. .......
.......
..
. ........
0
-1.0
-0.5
0.0
0.5
....
-1.0
-0.5
1.0
Probe Position, OR
..............
. ............... ........
0
0.0
0.5
1.0
-1.0
-0.5
2-
0.0
Probe Position.
42
E 2-
0.5
R
0. 2 sec Before
2-
1.0
I ug
..........
0
'
>
r
. ....... ................
........
0-
-1.0
................ ........
. .......
-0.5
0.0
Probe Position.
. .......
0.5
1.0
-1.0
-0.5
0.0
0.5
1.0
-1.0
........................
-0.5
0.0
0.5
1.0
2-
0
................ ........
........................
0
-1.0
-0.5
0.0
0.5
1.0
0
-1.0
-0.5
0.0
0.5
1.0
-1.0
Probe Position. rR
.......
-0.5
0.0
0.5
1.0
38
. . .........
.........
2-
....
E 2-
....
........ ........
Z'
..................7 .......
. ....... ..........................
. .......
0
-1.0
-0.5
0.0
0.5
Probe Position, fR
Figure 6
1.0
................. ........
0
-1.0
-0.5
0.0
0.5
1.0
-1.0
-0.5
189
0.5
0.0
1.0
maximum velocity below the center be. Although it is not seen here tis acceleration region is much
more pronounced when the gas superficial velocity s hgher than 0.55 M/s.
Beyond the transition region of 006 -
a sngle-phase liquid Dow. Although it is not axially symmetnc the velocity profile attains a serni-fully
developed condition. Ts
shows that the beyond the transition region the axial velocity is relatively
unaffected by the approaching gas slugs. The immediate influence of incoming slugs is confined to
a very small region in front of the gas slugs. Ts
[ 9.
However, this uniform flow situation is altered when the probe position approached to the
immediate wake re 'on of the gas slug in front of the liquid slug. Such a case can be observed around
0.24-second frame where one can see a slight deceleration. This may be attributed to the vortex
generated in front of liquid slugs at the wake of gas slugs. Detailed study of this region will be given
below.
faster moving gas slugs. Out of the acceleration zone the velocity profile exhibits a nearly fullydeveloped single-phase liquid flow character with the maximum velocity occuing slightly below the
center line.
190
....
... ...........
........
............
..
0
-0.5
0.0
0.5
1.0
-1.0
-0.5
0.0
0.5
0
.1.0
1.0
0.5
1.0
42
E
2-
...... I........
0.0
-0.5
Probe Position. rR
..........
>
0
-1.0
3-
2-
........
Z.2
................
. .. .I .. .
-1.0
-0.5
........
. .. .....
I
0.0
............
.............
0
0.5
1.0
. ......
-0.5
Probe Position, OR
0.0
0.5
1.0
-1.0
-0.5
Probe Position. /R
. .......
................
........
0.5
. ............... ........
0
0.0
1.0
......................... ........
-0.5
0.5
....... ........
2
-1.0
..
0.0
Probe Position,
.....
O
.1.0
.......
1.0
-1.0
-0.5
0.0
0.5
1.0
-1.0
-0.5
0.0
0.5
1.0
........
Z7
.2
................
........
. .......
0
.1.0
-0.5
0.0
0.5
Probe Position. rR
1.0
.. .
-1 0
.........................
. ......
. . . .. . .
-0.5
0.0
7.......................
0
0.5
1.0
-1.0
.0.5
0.0
Probe Position.
0.5
R
Figure 7 A Typical Velocity Profile Development for Liquid in the Wake Region
of the Gas Slua
191
1.0
'thin the liquid layer close to the gas slug nose whereas
Figure 9 shows the liquid velocity behavior within the liquid layer close to the gas slug tall. From
these figures it is clear that axial velocity profile never develops into a fully-developed profile. As
indicated in Figure 8, the liquid gradually decelerates fi-om the gas slug nose. The effect of interfacial
shear exerted by the faster moving gas slug seems to be confined to the very vicinity of the interface.
Such an effect causes a peculiar velocity profile with a point of inflection occurring between the wall
and gas-liquid interface. Figure 9 displays a very rapid deceleration in the liquid layer toward the
upstream of the slug tall. However, this brief deceleration is followed by a gradual acceleration
toward the slug nose.
The velocity evolution at location r/R = 0.6 is shown on Figure 10 where the time is
normalized by the liquid slug passage time T,,., and the gas slug passage time T,
slug and liquid layer, respectively. The ever developing nature of the liquid layer flow and the
relatively uniform flow behavior of the
liquid
liquid layer the liquid gradually decelerates from the slug nose. However,
acceleration is observed. On the other hand, out of the transition zones after the gas slug nose
and before the slug tall within the liquid slug, the velocity stays constant at this r/R = 0-6 location.
It Is interesting to note that there exists sgnificant differences between the liquid slug and liquid layer
velocity behaviors.
192
41
E
2-
Z
>
..............
>
I ............
0
-1.0
-0.8
-0.6
-0.4
-C-2
0.0
....
-1.0
I
06
04
02
. ............
-1.0
-0.8
2 -
E
2
2 -
06
.0.2
0.0
3 -0.10
0.0
.. .. ..
-0.8
>
... .
3 -0.12
2 .
........
......
0
>
....... .....
.
-1.0
-0.8
-0.6
-0.4
......
>
.............. .....
.......
>
-0.2
0.0
...... ............
0
-1.0
-0.8
06
04
02
00
-1.0
-0.8
Probe Position. rR
06
04
-0.2
0.0
. ...... ......
3-
2 . ......
......
2 -
......
.2
>
. ...... ......
...
-1.0
-0.6
-0.4
-0.2
0.0
-1.0
-0.
04
02
0.0
...... 7, ....
......
3 -
2 . .....................
3 -0.24
E
. ............ ......
.......
0
04
-0.2
Probe Position. OR
0.0
04
-0.2
0.0
......
.2
>
06
2-
......
-0.6
-0.8
Probe Position, OR
2 . ......
-0.8
-1.0
Probe Position, OR
-1.0
......
0
-0.8
Probe Position, UR
......
.2
......
2 . ...... ......
>
C
..
. .....
......
......
0-1.0
-0.8
06
04
-0.2
Probe Position. OR
0.0
-1.0
-0.8
06
-0.4
193
-0.2
0.0
3 -
3 -
2 . ...... ......
2 -
...
.2
>
.. ..
...... m...
0
.1.0
-0.8
-0.6
-0.4
-0.2
0.0
-1.0
-0.8
Probe Position. rR
3-
04
-0.2
0.0
-1.0
-0.8
Probe Position, rR
2 -
-0.6
......
-0.4
-0.2
0.0
3-
-0.6
..............
3-
E
21
2 -
...........
0
>
. .............
>
... ...
-1.0
-0.9
-0.6
-. 4
-0.2
0.0
-1.0
-0.8
Probe Position. OR
3-
. ...... ......
-0.6
04
-0.2
...
-1.0
Probe Position. OR
2 . ..............
2 .
. ..
-0.8
-0.6
-0.4
-0.2
0.0
3-
...........
0.18 sec
3-
2 -
........
Z'
>
....
-1.0
-0.8
-0.6
-0.4
-0.2
0.0
-1.0
-0.8
-0-6
-0.4
-0.2
0.0
-1.0
......
....... ........
.. .....
2-
30.22
........ .....
-0.4
-0.2
0.0
. ...... ......
-0.6
2--
...............
-0.8
......
.. ..... .....
2-
......
......
......
............
. ...........
0-1.0
0
-0.8
-0.6
-0.4
-0.2
Probe Position, OR
0.0
....
-1.0
0
0.8
06
-0.4
-0.2
Probe Position, OR
0.0
-1.0
-0.8
-0.6
-0.4
Probe Position, OR
194
-0.2
0.0
3
3O
2.5
..
0000
2
0
W
> 1.5
C:
co
1.0
....... . .
.......
_P
0 0 I
do
''CD
..
0 0
. .. ...
0.5
0.0
1.0
0.8
0.6
0.4
0.2
0.0
0.2
0.4
0.6
0.8
1.0
T/TFam
Figure 10. Axial Velocity Variations in Liquid Slug and Liquid Layer Under Gas Slug.
195
a 91'ven location the velocity gradually decelerates toward the gas sluc, tail. But it rapidly accelerates
towards the wake of the gas slug.
NOMENCLATURE
A
Voltage level
Subscript identifying the index number for large elongated slug bubbles
<i?
<jg>
Nib
Number of slug bubbles passing the probe in the total sampling time
N,,
Number of small bubbles passing the probe in the total sampling time
Pipe radius
Hot film sensor location measured alone, the vertical axis from the pipe center
Tfilrn
Tsjug
t2j-1
t2j
Time
U.,(r)
<U.,>
u,(rt)
Instantaneous axial velocity for the kth data point in the liquid phase
196
as bubble
U (r)
a(r)
<a(r)>
Kronecker delta
ACKNOWLEDGMENT
The work reported in this paper was performed under the auspices of the U.S. Department of
Energy, Office of Basic Energy Science. The authors would like to express their sincere appreciation
for the encouragement, support and technical comments on this program from Dr. O.P. Manley of
the U.S. DOE/BES.
197
REFERENCES
I
0. Baker, "Simultaneous Flow of Oil and Gas," Oil and Gas Journal, Vol. 53, pp. 185-190
(1954)-
2.
J. M. Mandhane, G.A. Gregory, and K. A. Aziz, "A Flow-Pattern Map for Gas-Liquid in
Horizontal Pipeline," Int. J. Multiphase Flow, Vol. 1, pp. 537-554 (194).
3.
Y. Taitel and A. E. Dukler, "A Model for Predicting Flow Regime Transitions in
Hon'zontal and Near Horizontal Gas-Liquid Flow," AICHE J. Vol. 22, pp. 47-55 1976).
4.
J. Weisman, D.Duncan, J. Gibson, and T. Crawford, "Effect of Fluid Properties and Pipe
Diameter on Two-Phase Flow Pattern in Horizontal Lines," Int. J. Multiphase Flow, Vol.
5, pp. 437-462 1979).
5.
A.E. Dukler and M.G. Hubbard, "A Model for Gas-Liquid Flow in Horizontal and Near
Horizontal Tubes," Ind. Engng. Chem. Fundam., Vol. 4, pp. 337-347 1975).
6.
M.K. Nicholson, K. Aziz and GA- Gregory, "Interrnittent Two-Phase Flow in Horizontal Pipes,
Predictive Models," Can. J. Chem. Engng., Vol. 56, pp. 653-663 1978).
7.
8.
D. Bornea. and N. Brauner, "Holdup of the Liquid Slug in Two-Phase Intermittent Flow," Int.
J. Multiphase Flow, Vol. I , pp. 43-49 1985).
9.
P. Andreussi and K.H. Bendiksen, "An Investigation of Vold Fraction in Liquid Slugs for
Horizontal and Inclined Gas-Liquid Pipe Flow," Int. J. Multiphase Flow, Vol. 15, pp. 937-946
(1989).
10.
P. Andreussi, K.H. Bendiksen and .J. Nydal, "Vold Distribution in Slug Flow," Int J.
Multiphase Flow, Vol. 19, pp. 817-828 1993).
11.
12.
13.
Z. Ruder and T. J. Hanratty, "A Definition of Gas-Liquld Plug Flow in Horizontal Pipes,"
Int J Multiphase Flow, Vol. 16, pp. 233-242 1990).
14.
G. Saether, K. Bencliksen, J. Muller and E. Froland, "The Fractal Statistics of Liquid Slug
Lengths," Int. J. Multiphase Flow, Vol. 16, pp. II 17-1126 (1990).u
198
15.
16.
Z. Fan, W. P. Jepson and T. J. Hanratty, "A Model for Stationary Slugs," Int. J. Multiphase
Flow, Vol. 8, pp. 477-494 1992).
17.
Z. Fan, Z. Ruder and T. J. Hanratty, "Pressure Profiles for Slugs in Horizontal Pipelines," Int.
J. Multiphase Flow, Vol. 19, pp. 421-437 1993).
18.
0. Kvernvold, V.
doy, T. Sontvedt, A. Saasen and S. Selmer-Olsen, "Velocity Distribution
in Horizontal Slug Flow," Int. J. Multiphase Flow, Vol. 10, pp. 441-457 1984).
19.
M. Kawaji, M. All, A. Ciastek, and C. Lorencez, "Study of Liquid Flow Structure in Horizontal
Cocurrent Gas-Liquid Slug Flow," Proc. ANS-THD, Vol. 8, pp. 79-88, 1995, Natl. Heat
Transfer Conf, August 59, 1995, Portland, Oregon.
20.
Y.Y. Hsu, F.F. Simon, and R.W. Graham, "Application of Hot-Wire Anemornetry for
Two-Phase Flow Measurements such as Void Fraction and Slip Velocity," Proc. ASN/1E Winter
Meeting, Philadelphia, PA 1963).
21.
22.
23.
24.
R. Abel and F.J. Resch, "A Method for the Analysis of Hot Film Anemometer Signals in
Two-Phase Flows," Int. J. Multiphase Flow, Vol. 4 pp. 523-533 1978).
25.
S.K. Wang, S.J. Lee, O.C. Jones and R.T. Lahey, "Local Void Fraction Measurement
Techniques in Two-Phase Bubbly Flow Using Hot-Film Anemometry," Proc. 22nd Heat
Transfer Conf., Niagara Falls, August 5-8, 1984.
26.
S.K. Wang, S.J. Lee, O.C. Jones and R.T. Lahey, "3-D Turbulence Structure and Phase
Distn'bution Measurements in Bubbly Two-Phase Flows," Int. J. Multiphase Flow, Vol. 8, pp.
327-343, 1987).
27.
T.J. Liu and S.G. Bankoff, Structure of Air-Water Bubbly Flow in a Vertical Pipe - . Liquid
Mean Velocity and Turbulence Measurements," Int. J. Heat
Mass Transfer, Vol. 36, pp.
1049-1060, 1993).
199
28.
T.J. Liu and S.G. Bankoff, "Structure of Air-Water Bubbly Flow in a Vertical Pipe - IL Void
Fraction, Bubble Velocity and Bubble Size Distributions," Int. J. Heat Mass Transfer, Vol.
36 pp. 1061-1072, 1993).
29.
30.
3 1. T.G. Theofaneous and J. Sullivan, "Turbulence in Two-Phase Dispersed Flows," J. Fluid Mech.,
Vol. 116, pp. 343-362 1982).
32.
33
34.
200