1. Introduction
The global interest in the study of multiphase flow has progressively grown in the last century since multiphase mixtures are widely encountered nowadays in industrial processes and they represent a crucial matter to be analyzed for the development and the optimization of many industrial processes [
1,
2].
In particular, from the 1970s onward, multiphase flow has been studied extensively in order to provide engineers with suitable design tools [
1]. A significant aspect in the characterization of a multiphase flow is represented by the flow pattern, i.e., the variety of configurations taken by the two phases during the flow. Different flow patterns strongly influence transport phenomena and therefore they have been recognized and classified and integrated, where significant, within the models in order to achieve a consistent improvement in their accuracy.
Regarding air–water adiabatic flows in horizontal ducts at atmospheric pressure, object of this paper, Mandhane et al. [
3] provided a useful map, identifying six different flow patterns as a function of the superficial velocities of the phases.
This work analyses the stratified region of the flow map, corresponding to low superficial velocities of both gas and liquid phases, which is of particular interest in many industrial areas. As remarkable examples can be mentioned gas transportation lines [
4], geothermal plants and particular operating conditions of pressurized water nuclear reactors and chemical plants [
2,
5]. From the early work by Mandhane et al. [
3], several authors provided a more detailed classification of the stratified flow pattern. In particular, Spedding and Nguyen [
6], Andritsos and Hanratty [
7] and, more recently, Tzotzi et al. [
8] identified four sub-regimes mainly differing in the shape of the interface. Specifically, at constant liquid superficial velocity, raising the gas superficial velocity, the smooth gas-liquid interface is perturbed such that initially two-dimensional “ripple” waves appear, evolving to three-dimensional “roll” waves and then growing to inertial waves. After that, liquid drop entrainment in the gas or air bubble entrainment in the liquid [
9] may take place as a beginning of transition to pseudo-slug and slug flow [
10]. Particular attention has been devoted to the mechanism of interfacial waves formation and their description [
10,
11,
12,
13], since they significantly affect interphase mass, momentum and heat transfer, and in turn the global fluid dynamic behavior characterized by pressure drop and void fraction.
Turning the attention to modelling, the typical engineering approach is based on the one-dimensional two-fluid model [
14], involving in the most general case six equations for the conservation of mass, momentum and energy of the two phases. In many of the works published in the literature and concerning the present paper [
15,
16,
17], the following assumptions greatly simplify the mathematical statement of the problem:
Steady state;
Fully developed flow conditions;
Adiabatic or isothermal flow;
Constant properties of the phases.
In particular, regarding stratified gas-liquid flows without heat transfer, phase change and with negligible entrainment of the phases, the momentum balance equations can be solved to derive the pressure gradient and the phase holdup, as reported in the early models by Govier and Aziz, [
15] Agrawal et al. [
16], Taitel and Dukler [
17], which only differ in minor details. Other authors studied the stratified flow regimes with particular attention to the interfacial stability during stratified to wavy [
18,
19,
20] and wavy to slug [
21] flow pattern transition. Over the years, particular attention was also paid to the analysis of flow conditions with a small amount of liquid, i.e., void fraction higher than 90% (e.g., [
22,
23,
24,
25,
26,
27,
28]), with the aim of modelling gas flowlines with a very small amount of liquid deposition/accumulation. Actually, as remarked above, the presence of waves at the interface significantly changes the interfacial shear stress and hence the momentum transfer between the phases, leading to very different behaviors of both phase holdups and pressure gradients compared to the case of smooth interface. In this respect, Tzotzi and Andrisitos [
29] recently provided semi-empirical models for the criteria of transition between smooth, two-dimensional and three-dimensional waves flow patterns along with empirical expressions of the interfacial friction factor
for each flow pattern. In particular, denoting
the gas to wall friction factor, it was shown that the wavy flow patterns are characterized by a ratio
, which is interpreted as a major case of deviation from the original Taitel and Dukler [
17] model. On the other hand, it is commonly assumed that for a smooth interface
. However, this statement might not be fully justified since it implies a stationary interface [
30]. On the other hand, most of the experimental data available for model validation fall near the transition between stratified and wavy flow, as shown in
Figure 1.
Accordingly, this paper focuses on operating conditions in a range of low gas superficial velocities, corresponding to stratified smooth gas-liquid interface, which is poorly studied, but significant for the application to gas flowlines where water accumulation is one of the issues affecting flow assurance [
31]. Specifically, measurements of both pressure gradient and void fraction are here reported and compared with available models. Hence, the two-fluid model is adopted to determine the interfacial shear stress and an empirical expression of the interfacial friction factor is proposed as a closure relationship.
Eventually, the same model was applied on a broader set of data, in which both the wavy and transition flow conditions and the inclination are taken under consideration, to confirm whether or not the correlation found can be applied to a nearly universal database, at least dealing with stratified flows.
2. Materials and Methods
The apparatus schematically reported in
Figure 2 was used to run the experiments. It can be divided into four sections: water feed and discharge, air feed and mixer, test section, and data acquisition system. The liquid is supplied from the bottom of a
tank by means of a centrifugal pump (CALPEDA
;
). The volume flow rate, measured by a float-type flow meter (ASAMETRO P13-2800,
, ±3% full scale), is set through a valve. The two-phase flow is obtained by a mixer (
Figure 3a) where compressed air is injected in the liquid flow.
The air volume flow rate is measured by float-type flow meters (
Table 1), the reading of which is suitably corrected to account for pressure and temperature deviation from standard conditions. In order to reduce the measurement uncertainty in both the air and water flow rates, the flow meters have been calibrated in the considered range as reported in Colombo et al. [
32].
Once air and water pass through the mixing section, the resulting two-phase flow is directed to the test section, consisting of a 24 m long horizontal pipeline composed of 13 PMMA, 60 mm bore tubes, connected by means of sealed flanges, suitably designed to minimize the misalignment within a tolerance of 0.5 mm (
Figure 3b). Furthermore, the pressure taps have been placed at a distance of at least 10 diameters from the flanges, to further dampen the disturbance caused by connections, as confirmed by the smooth behavior of the pressure along the pipeline. The loop closes on the supply tank, where the two-phase mixture is immediately separated by gravity. The test section is equipped with 4 pressure taps evenly distributed along the pipeline and connected to a differential pressure transducer.
The pressure taps are placed facing upwards as the transducer works with air only. The first tap is located 15.6 m downstream of the mixing section in order to consider fully developed flow conditions. All sensor signals are collected by means of a National Instruments acquisition board (1 kHz sampling frequency cf.
Table 2) and processed by LabVIEW®.
The operating conditions were selected in order to investigate the range shown in
Figure 1. The superficial velocities are specified in
Table 3.
Measurement repeatability of the pressure drop was checked by at least ten repetitions for each operating condition in different days: the minimum and maximum deviations turned out to be 1.3% and 11.5%, respectively, with the highest values corresponding to flow regimes with perturbed interface (highest gas superficial velocity).
Concerning the void fraction, a resistive probe was used to measure the liquid layer height. The measuring apparatus (electrical scheme in
Figure 4) is a modified version of the one developed by Kim and Kang (1992) [
33] and reported also in [
34] and [
35]. It is a single-tip, external-electrode resistive probe that consists of: (1) a needle electrode introduced radially from above into the flow, through an open pressure tap and by means of a probe support endowed with a ruler to provide the immersion depth, and (2) a second wire, submerged through a second pressure tap into the liquid phase, that acts as a reference electrode. Both the reference electrode and the probe stem itself are very thin, so their disturbance on the liquid flow is completely negligible. The needle electrode is insulated, using a non-conductive ceramic painting, apart from the final tip that is the sensing part. The reference electrode is connected to a +5 V power supply, while the needle is connected to an analogue input of the acquisition board and grounded through a resistor. The choice of connecting the reference electrode and not the probe to the +5 V power supply is to avoid any possible electrowetting phenomena, with negatively charged water attracted towards the positively charged probe and consequent distortion of the measure.
Since air is not electrically conductive, such a probe allows to sample the state density function (at a frequency of 10 kHz): the voltage signal is zero until the needle touches the liquid phase, then the circuit closes and the signal jumps to a positive value. A 214 Ω resistor was selected as a suitable compromise between the reductions of noise and of signal amplitude; in any case the precise value is not of importance as the signal is only used in a bimodal way (low: air, high: water). With such a value, the signal amplitude is about +3 V.
It is then possible to measure the difference between the pipe inner diameter and the maximum water layer depth
. Accordingly, referring to
Figure 5, void fraction values are obtained from Equations (1) and (2), where
is the angle subtended by two radii to the ends of the interface perimeter
.
The measurements were taken at two different pressure taps and at least five repetitions were done for each one, to assess repeatability. In fact, the outlined method rigorously applies with a smooth interface, which is generally confirmed by visual inspection. In any case, in the presence of weak waviness, the measurement at two different pressure taps enables the determination of the average .
Concerning the accuracy of the liquid height measurement: before the tests, the probe is fully inserted in the duct so that its tip touches the duct bottom and a reference level can be evidenced on the ruler. Then the probe is retracted by the full diameter length and progressively re-inserted until the tip touches the liquid. The tip sinking is visually evaluated. Precision is quite high also when comparing measures performed by different operators and an accuracy of 0.5 mm can be reasonably assumed. Furthermore, an uncertainty analysis on derived quantities was performed, and the results are summarized in
Table 4. The reader may refer to
Appendix A for an extended discussion.
3. Results and Discussion
3.1. Pressure Gradient
The pressure gradient was found to linearly decrease along the flow direction, which confirms the achievement of fully developed flow conditions. Accordingly,
Figure 6 reports the pressure gradient as a function of the water cut (Equation (3)):
showing a power-law dependence with an exponent equal to −1.7, regardless the liquid superficial velocity.
The experimental results were then compared to the predictions of available models suited to the range of the investigated operating conditions. In particular, Tzotzi and Andritsos [
29] developed the following criterion to mark the transition from smooth stratified flow to 2D wavy flow (Equation (4)):
being
and
the couple of superficial velocities marking the transition.
Based on such a criterion, about 25% of the data should correspond to stratified flow with smooth interface. On the other hand, visual observations show that the interface is usually smooth. The onset of 2D waviness is observed only for gas superficial velocities of about 2 m/s, i.e., about the upper limit of the investigate range, which corresponds to about 12% of the whole data.
Accordingly, the model proposed by Tzotzi and Andritsos [
29] for the prediction of the pressure gradient in the case of smooth interface is the same as developed by Taitel and Dukler [
17]. As mentioned in the introduction, it consists basically of a two-fluid formulation of the momentum equation in the assumption of one-dimensional fully developed steady flow with constant properties, which reads
Figure 7 shows a sketch of the forces accounted in the momentum balance.
Rearranging the terms for each phase, the following system of equations results
The interaction between the phases is represented by the interfacial shear stress
, which is usually given by semi-empirical expressions involving the interfacial friction factor. In particular, according to Taitel and Dukler [
17]
where for smooth interface it is assumed
with
and
for
and
and
otherwise.
However, as mentioned in the introduction, Equation (9) assumes that the interface is stationary regardless the relative velocity of the phases, which, according to Ullmann and Brauner [
30] is one of the pitfalls causing wrong predictions even in the case the interface is smooth and plane and the flow of both phases is laminar. For this reason, these authors revised the definition of the interfacial friction factor by introducing correction factors to account for the interaction between the flows in the two layers. However, the validation of the new closure equation relies upon the database of Espedal [
36], which involves a range of superficial velocities mostly corresponding to wavy flow patterns (see
Figure 1). Moreover, most of the empirical expressions of
available in the literature (see
Table 5) are only validated for wavy flows up to the transition to slugs.
No matter the model adopted for the interfacial shear stress, Equations (6) and (7) are solved simultaneously to determine both the pressure gradient and the void fraction.
Figure 8 shows the mean absolute percentage error (MAPE) as a function of the gas superficial velocity. It is seen that MAPE values lower than 30% are found for
. Below this value, the disagreement strongly increases as far as the gas superficial velocity is lowered and the liquid superficial velocity is increased. It is worth noting that the Ullman and Brauner model [
30] performs better than the Taitel and Dukler [
17] for
, but MAPE remains unacceptably high for
and
, suggesting that the cause of the disagreement still lies in the representation of the interfacial shear stress.
Accordingly, void fraction measurements were also taken in order to solve the two-fluid model, Equations (6) and (7), for the interfacial shear stress and hence provide an empirical correlation to be adopted in the considered range of the gas superficial velocities.
3.2. Void Fraction
Since both Taitel and Dukler [
17] and Ullmann and Brauner [
30] models return unsatisfactory predictions for the void fraction, especially at low air-flow rate, three correlations available in the open literature were implemented (see
Table 6), i.e., Armand [
37], being one of the simplest general correlations implemented for two-phase flow, Guzhov [
38], as it was developed specifically for plug and stratified flow in horizontal pipes, and Rouhani [
39], which is among the best predictors over a very large data set [
40].
The measured void fraction is reported against the volume quality in
Figure 9 and
Figure 10, where it is also compared with some models from the literature. However, as shown in
Figure 9 and
Figure 10, none of the implemented correlations satisfactorily predicted the void fraction in the whole range of superficial velocities, even though a similar trend can be found with the Armand correlation.
Furthermore, the agreement with the other models significantly improves as the liquid superficial velocity increases.
Hence, the collected data were processed by the drift-flux model formulated by Zuber and Findlay [
41], because this approach was proved to be the most effective in modelling void fraction [
40].
Accordingly,
Figure 11 shows the actual velocity of the gas phase, evaluated as
against the mixture velocity
, for the different liquid superficial velocities.
It is observed that each data series is linearly correlated, which is consistent with the theory. The slope of the straight line represents the so-called distribution parameter , which accounts for non-uniform phase distribution on the pipe cross-section (uniform flow: ), and ranges from 1.20 to 1.30 as increases.
The intercept represents the weighted mean drift-velocity , which accounts for the relative velocity of the gas phase with respect to the mixture; it is almost constant, with a mean value of 0.06 m/s.
Including all the data, the drift flux model
with
and
is able to predict the data with MAPE below 8%.
3.3. Interfacial Shear
Under the assumption of fully developed one-dimensional steady flow the pressure gradients in the two-phase are equal, the momentum equations of the two-fluid model Equations (6) and (7) were rearranged to express the interfacial shear stress as a function of the overall pressure gradient and void fraction (experimentally measured):
The gas-wall shear stress,
, was computed according to the Churchill correlation [
42] as suggested by Oliemans and Pots [
43]:
The interfacial shear stress is plotted against the velocity difference between gas and liquid in
Figure 12. A power law fitting with exponent 1.47 and a regression coefficient of 0.019, furthermore it shows correlation coefficient equal to 0.99.
Accordingly, the interfacial friction factor is defined as
and plotted against the superficial gas Reynolds number as shown in
Figure 13, where it is compared with some literature models suitably developed for stratified smooth/wavy flows, namely Cohen and Hanratty [
44], Agrawal [
16], Kowalski [
45] and Hart et al. [
24]. The fitting curve is:
Noticeably, the result falls in the mid-range of the literature models with a very similar trend compared to both Agrawal [
16] and Kowalski [
45] correlations, which are in the same form of a power law with the Reynolds number raised to an almost identical exponent.
3.4. Model Extension with Espedal Database
Eventually, further validation was performed based on the Espedal database reported in [
36]. A preliminary assessment was made, since all measurements were performed with inclined pipelines: for our purposes, only data concerning slightly inclined pipes were considered, i.e., inclination angle
with respect to the horizontal equal to: −0.5°, −0.25° and −0.104°. The operating conditions are reported in
Table 7.
The two-fluid model momentum equations were rearranged to account for pipe inclination. It resulted that the interfacial friction factor for small inclination angles does not vary considerably compared to the horizontal layout, so that Equation (14) would still apply. The Espedal data [
36] were processed as indicated above in the Section “Interfacial shear” in order to evaluate the interfacial friction factor. The comparison, reported in
Figure 14, shows that Equation (14) tends to overpredict the data, with MAPE = 35% and MPE = −34%. The data refer to the highest value of the inclination angle, which represents the worst case.
A correction was then introduced in Equation (14), including the Froude number, i.e., the ratio between inertia and gravity forces. For horizontal pipes () the corrective term is equal to one, i.e., no correction applies.
The modified interfacial friction factor
resulted as follows:
The two unknown parameters, i.e., and , were determined by regression as 1.3 and 0.21 respectively.
Figure 15 reports the comparison between Equation (15) and Espedal data [
36] for the highest inclination angle. The agreement is satisfactory with MAPE = 11% and MPE = 2%. However, the model is not able to explain the non-monotonic behavior derived from the experimental data.
4. Conclusions
A characterization of two-phase air–water adiabatic flows in a horizontal pipe 60 mm i.d. was experimentally performed, focusing on the pressure drop and void fraction in a range of water superficial velocities from 0.03 to 0.06 m/s and air superficial velocities from 0.41 to 2.31 m/s. To the best knowledge of the authors, similar experimental conditions were poorly investigated so far in the open literature, even though they might be typically encountered in pipelines used for industrial applications. Starting with the pressure drop measurements, the comparison with two-fluid models specifically developed for stratified and wavy flow regimes, showed a fair agreement for , whereas below this value the deviation strongly increases as far as the gas superficial velocity is lowered and the liquid superficial velocity is increased. It was argued that such a disagreement derives from the closure equations, represented by empirical expressions of the interfacial shear, which turn out to be developed for a higher range of gas superficial velocities. Accordingly, void fraction measurements were taken by means of a resistive probe and they correlated well to the gas volume fraction by means of the drift-flux model (, ) with MAPE < 8%. Hence, the momentum equations of the two-fluid model were rearranged in order to express the interfacial shear stress as a function of the overall pressure gradient and void fraction. A power-law correlation (exponent 1.47, ) between the interfacial shear stress and the gas-liquid velocity difference was found. The resulting interfacial friction factor was reported against the superficial gas Reynolds number. The comparison with some literature models shows a similar trend and values lying in the mid-range. Furthermore, the model presented was extended to slightly inclined pipelines, using the Espedal database. In particular, the model equation for the interfacial friction factor was adjusted in order to take into consideration the inclination angle through the Froude number, providing a good agreement (MAPE ~ 11%) with the data.