Next Article in Journal
Equilibrium Study, Modeling and Optimization of Model Drug Adsorption Process by Sunflower Seed Shells
Previous Article in Journal
Method of Multi-Angle Transmission Radiowave Tomography of Dielectric Objects
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Flow Regime, Slug Frequency and Wavelet Analysis of Air/Newtonian and Air/non-Newtonian Two-Phase Flow

by
Munzarin Morshed
1,
Muhammad Saad Khan
2,
Mohammad Azizur Rahman
2,* and
Syed Imtiaz
1
1
Department of Process Engineering, Memorial University of Newfoundland, St. John’s, NL A1C 5S7, Canada
2
Department of Petroleum Engineering, Texas A&M University at Qatar, Doha 23874, Qatar
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(9), 3272; https://doi.org/10.3390/app10093272
Submission received: 7 April 2020 / Revised: 1 May 2020 / Accepted: 4 May 2020 / Published: 8 May 2020
(This article belongs to the Section Energy Science and Technology)

Abstract

:

Featured Application

Application of Newtonian and non-Newtonian fluids in horizontal flow loop conditions for directional drilling operations.

Abstract

This study focused on gas/Newtonian and gas/non-Newtonian two-phase horizontal fluid flow behavior by analyzing their flow regime identification and flow structural analysis on a horizontal flow loop apparatus. This involved the recognition of two-phase flow regimes for this flow loop and validation with existing flow maps in the literature. In addition, the study included flow pattern identification via wavelet analysis for gas/Newtonian and gas/non-Newtonian two-phase fluid flow in a horizontal flow loop apparatus. Furthermore, the study was extended to the detailed examination of slug frequency in the presence of air/Newtonian and air/non-Newtonian fluid flow, and the predicted slug frequency model was applied to the studied systems. The obtained results suggest that the flow regime maps and slug frequency analysis have a significant impact. The obtained pressure sensor results indicate that the experimental setup could not provide high-frequency and high-resolution data; nevertheless, wavelet decomposition and wavelet norm entropy were calculated. It offered recognizable flow characteristics for bubble, bubble-elongated bubble, and slug flow patterns. Therefore, this study can provide deep insight into intricate multiphase flow patterns, and the wavelet could potentially be applied for flow analysis in oil and gas pipelines.

Graphical Abstract

1. Introduction

Multiphase flows are considered as complicated flow phenomena compared to a single flow [1]. There are essential features of multiphase flows whose modeling outcome is contentious and structural explanation is still unexplored. The most common type of multiphase flow in almost all chemical, petroleum, and production industries is the two-phase gas/liquid flow [2,3]. Different forms of flow patterns may be observed when two or more phases flow simultaneously [2]. In pipe cross-sections, unpredictable turbulent flow structures generating highly asymmetric volume distribution are a challenge in experiential investigations [4,5]. This kind of unstable flow condition complicates the measurement process and may create difficulties in capturing actual flow conditions. There are also instances where the existing theoretical solution or experimental results cannot describe specific physical properties [1,6].
Flow regimes or patterns are amongst the most significant aspects of multiphase flow. The structural flow distribution of different phases in a pipe is known as the flow pattern or flow regime. The flow regime depends on the inertia force, buoyancy force, flow turbulence, and surface tension, which are altered by fluid properties, flow rates, pipe diameter, and pipe predilection [7,8]. Different forms of flow patterns may be observed when two phases of gas/Newtonian and gas/non-Newtonian flow simultaneously. Some of the conventional distributions are stratified flow, where the liquid and gas phase is separated, and the gas flows on top as it is lighter than the liquid; bubbly flow, where there is a dispersion of small-sized bubbles within the liquid; slug flow, in which each gas bubble forms a large slug shape (often a bullet shape); and annular flow where liquid flows as a film on the wall of the pipe. For gas/Newtonian and gas/non-Newtonian flow, several flow maps exist to predict the flow patterns. The Taitel and Dukler [7] and Mandhane et al. [9] flow maps for gas/Newtonian flow, and the Chhabra and Richardson [10] flow map for gas/non-Newtonian flow are the most frequently used.
Among the flow regimes, slug flow is the most frequent two-phase flow phenomena experienced in horizontal or near horizontal pipelines in practical applications [11,12,13]. The presence of slug flows in a pipeline can trigger a great range of design and operational problems, including the exertion of kinetic force on fittings and fittings, pressure cycling, control instability, and inadequate phase separation. If the flow rate is raised, much of the liquid is swept out, and the liquid handling capacity of the receiving facilities may be overwhelmed. In the case of a slug greater than the slow catcher capacity, the facilities can be flooded and destroyed. Therefore, multiple operational problems, such as pipeline network instability and damage to equipment due to high-pressure fluctuation or vibration of the system, can be caused by slug flow [14,15]. This can also be termed a water hammering effect. Therefore, slug flow and slug frequency analysis has been one of the primary research interests in multiphase flow.
To measure the hydrodynamic behavior of multiphase flow, it is essential to understand flow pattern under specific flow conditions [2,16]. Two-phase flow implies gas and liquid flow through a pipeline system, simultaneously. The gas and liquid interface is deformable, so it is hard to predict the region occupied by a gas or liquid phase [17]. When two phases flow through a pipeline, different types of interfacial distribution can form. The variety of flow patterns mostly depends upon the input flux of the two phases, the size and assembly of the pipe, the physical properties of the fluid, etc. Many experimental studies on gas/Newtonian or solid/Newtonian fluid flow have been conducted [18,19,20,21,22,23,24]. Nevertheless, a limited number of studies have observed non-Newtonian multiphase flow [25,26].
Xanthan gum is the most commonly used industrial biopolymer for thickening and stabilizing an aqueous system. Xanthan gum solution has significant thermal stability with pseudoplastic properties [27]. Due to these properties, it has a substantial application in petroleum industries. In oil industries, xanthan gum is widely used in drilling fluid [28]. Xanthan gum is also broadly used in the food industry, cosmetics, and pharmacological products [27,29].
Tutu [30] and Drahos et al. [31] characterized two-phase horizontal flow regime pressure fluctuation. Drahos et al. [31] used a probability density function (PDF), where a strain gauge pressure transducer was used in a 50 mm inner diameter (ID) Perspex pipe. Sun et al. [32] used norm entropy wavelet decomposition to analyze gas/liquid two-phase flow pressure signal data across a bluff body. Here, the inner pipe diameter was 50 mm, and, from a piezoresistive differential pressure sensor, the pressure signals were analyzed using four levels and four scales of Daubechies wavelet (db4), which provided 16 wavelet packet coefficients. This study [32] also suggested an entropy-based two-phase flow map with an identification rate of 95%. Blaney [33] used gamma rays to identify flow regimes, and continuous wavelet transforms to analyze gamma count data. Park and Kim [34] applied the wavelet packet transform to analyze pressure fluctuations in a vertical bubble column.
Furthermore, De Fang et al. [35] also used wavelet analysis to understand the gravity differential pressure fluctuation signal perpendicular to the horizontal flow of different flow patterns and the flow pattern transition of gas/liquid two-phase flow in the horizontal pipe. Here, the Haar wavelet with six levels was used to decompose the pressure signal, and the energy value was then obtained for each scale. For identifying a two-phase flow regime, Elperin and Klochko [36] also used an eight-level db4 wavelet transformation to process time series of the measured differential pressure fluctuation. In this study, using pressure transducer signal data of different flows, the pattern was analyzed using wavelet transformation to determine the pressure signal characteristics of various flow regimes.
To identify the two-phase flow regime mapping and avoid slug flow in horizontal pipelines, a clear understanding of the Newtonian and non-Newtonian fluid flow is necessary for an accurate and inherently safer design for any flow system. The uniqueness of this study is that the experiments were performed in a setup with 73.66 mm ID and approximately 19 m flow loop in the presence of gas/Newtonian and gas/non-Newtonian two-phase horizontal flow. This flow loop has a horizontal, vertical, and inclined test section connected; however, the study is only on horizontal flow aspects. Another primary focus of this study is to understand the characteristics of slug frequency analysis for both air/Newtonian and air/non-Newtonian flow conditions in the presence of Xanthan gum as a non-Newtonian fluid. This study also recognizes the characteristics of wavelet analysis for both air/Newtonian and air/non-Newtonian flow conditions, i.e., water and Xanthan gum.

2. Materials and Methods

2.1. Materials

Water is used as a Newtonian fluid in this study. For non-Newtonian fluid experiments, a 0.1 wt% solution of biopolymer, i.e., Xanthan gum, was used (obtained from Merck, Germany). To conduct two-phase experiments, compressed air was used for gas/Newtonian fluid and gas/non-Newtonian experiments. The air was compressed via an in-house compression facility.

2.2. Experimental Methods

The experiments were performed in a flow loop system that has a horizontal, vertical, and inclined section. However, in this paper, we are only considering the 4-m horizontal section as our test section. The experimental setup is a 60-m-long closed cycle system for water and open cycle system for air. The liquid is pumped by a 5 HP pump that creates the required large volume water flow through DN80 or 2.9 ID PVC clear pipes. The airline of the flow loop comprises of DN15 and DN25 mild steel pipes, which supply air from the lab air supply at 670 kPa (100 psi) shut-in pressure. It also includes a DN 25 ball check valve just before the air and the liquid mixing zone to prevent any liquid from entering the air pipeline. There are two Omega PX603-100G5V pressure transducers with a range of 0 to 100 psi in the 2-m long horizontal test section. There were some specific experimental conditions used for this setup. The airflow range was approximately 85 L/min to 3300 L/min, and the water flow ranged from almost 250 L/min to 850 L/min. Within this range, the experimental setup mostly provides slug flow for two-phase flow, as well as providing bubble flow and wavy flow at some intervals. Figure 1 presents a schematic representation of the experimental setup. For this study, both gas/Newtonian and gas/non-Newtonian fluid flow cases were considered. A Universal Data Acquisition System from National Instruments was used to collect all types of data from the flowmeter and sensors, which consisted of four NI 9219 universal modules with four channels, each providing 100 samples per second. The experimental conditions are reported in Table 1, including brief details of the rheological model used for xanthan gum solutions.
Slug flow is the most frequent two-phase flow phenomenon experienced in horizontal or near horizontal pipelines in practical applications. Slug flow in pipes is encountered in different fields, such as the production and transportation of oil and gas, food industry, chemical industry, and other industrial applications. Slug frequency, in addition to water hammering, leads to various operational problems, such as pipeline network instability, equipment damage, pressure fluctuations, and vibration of the system. In the oil and gas production industries, slug flow also influences the internal corrosion rate increase of carbon steel pipelines. Slug flow creates high turbulence, which can break the pipe wall inhibitor’s protection layer [37].

2.2.1. Slug Velocity

In 1965 and 1977, Hubbard [38,39] and Otten & Fayed [40], respectively, experimentally measured the slug velocity by observing a particular slug movement in test sections. They both obtained a relationship between slug velocity and no-slip mixture velocity by plotting the experimentally measured slug velocity against the no-slip mixture velocity. The slug flow model of Hubbard [39] showed better agreement at higher slug velocity. Hubbard [39] described the relationship as:
v s = 1.25   v m
Hubbard & Dukler [39] also predicted the actual average gas velocity as in Equation (2), which also agreed with other experimental data [41].
v G = 1.19   v m
It is assumed that the Hubbard & Dukler [39] slug flow model was verified based on one dominant presupposition that the liquid slug velocity and the maximum gas phase velocity should be similar. Therefore, theoretically, the no-slip mixture velocity should be equal to the slug velocity:
v s v m = C
where C is a constant. Theoretically, C is assumed to be 1.0 for air–water two-phase flow. Hubbard [38], Rosehart et al. [42], and Gregory & Scott [41] considered C values of 1.25, 1.26, and 1.35, respectively for air–water flow. These C values may have varied because of different experimental setups and conditions [40]. For non-Newtonian/air two-phase flow, Otten & Fayed [40] compared their results with those of Rosehart et al. [42]. The obtained results suggested that the air/Carbopol 941 concentration reported by Otten & Fayed [40] increased from 0.75% to 0.2% and C values increased from 1.36 to 1.41, whereas for the same concentration, the C values of Rosehart et al. (1975) varied from 1.54 to 1.98.

2.2.2. Slug Frequency

Different correlations can predict slug frequency. However, Gregory & Scott [41] suggested a trustworthy velocity-dependent empirical equation where slug frequency was correlated with a form of Froude number, which is described as:
N f s = v l g d [ ( v m 0 ) 2 v m + v m ]
Here, v m 0 was taken as 6 m/s and from the graphs of the slug frequency vs. slug Froude number by Gregory & Scott [41], the following equation was derived:
f s = 0.0157 ( N f r ) s l u g 1.20 sec 1   .  
From Equation (5), Gregory & Scott [41] described a slug frequency correlation based on their liquid-gas two-phase flow experimental data where water and carbon dioxide are used in a 19-mm ID pipe:
f s = 0.0226 [ v l s g d ( 19.75 v m + v m ) ] 1.2
where v m   and   v l s represent the mixture velocity and superficial liquid velocity of liquid and gas, respectively. Therefore, this slug frequency can be combined with the Froude number established for superficial fluid velocity.
Greskovich & Shrier [43] reorganized the Gregory & Scott [41] correlation as follows:
f s = [ 0.0425 v l s v m ( 2.02 d + v m 2 g d ) ] 6 5
Later, Zabaras & others [44] described another correlation based on 399 data points with the smallest average absolute error and standard deviation for both horizontal and inclined pipe flow. This correlation is a modification of the Gregory & Scott [41] correlation, using the imperial unit, which is shown in Equation (8), where θ is the inclination angle. The experiment was conducted in a 2.54 cm and 10.16 cm ID pipe with air and water.
f s = [ 0.0425 v l s g d ( 1 0.0506 v m + v m ) ] 6 5 [ 0.836 + 2.7 s i n 0.25 θ ]
Heywood & Richardson [45] determined the liquid volume fraction for air–water two-phase flow utilizing the gamma-ray technique in a 4.191 cm ID horizontal pipe. To determine the liquid volume fraction, they used a power spectral density function and probability density function. These features are also helpful to decide on different slug flow characteristics, such as the value of the average film and slug volume fraction, average slug frequency, and average slug length. The slug frequency correlation was determined by curve fitting the data, resulting in Equation (9), in which λ is the liquid volume fraction, where λ = v l s / ( v l s + v g s ) , and d is the pipe diameter.
f s = [ 0.0462   λ ( 1 0.0126 d + v m 2 g d ) ] 1.02
Shea et al. [46] developed a correlation describing slug frequency as a function of pipe length. In Equation (10), v l S is the superficial liquid velocity, d is the pipe diameter, and l p is the pipe length. This correlation is based on the curve fitting of field and laboratory data, and is not based on theoretical analysis. In this equation, it is also shown that the slug frequency is inversely dependent on the pipe length lp, which does not agree with the other theoretical analysis. According to Al-Safran [14], the OLGA 2000 slug tracking model had some time delay problem between two slugs; to solve this issue, the Shea et al. [46] correlation was initially used. Moreover, pipe length can be questionable for a long-distance transmission system with hilly conditions:
f s = 0.47 [ (   v l s   ) 1.5   l p 1.1 d 2.4 ] 0.5
Picchi et al. [47] described a slug frequency equation that considers the rheology of the shear-thinning fluid. This equation is a modified version of the Gregory & Scott [41] correlation. In Equation (11), R e w = ρ w v l s d µ w is the water Reynolds number and R e n = d n v l n ρ n m 8 n 1 ( 1 + 3 n 4 n ) n 2 is the power-law fluid Reynolds number under superficial conditions, where n and m are fluid behavior indexes.
f s = 0.0448   [   v l n g d ( 32.2014 v m n + v m n ) ] 0.88 n 2.85 ( R e n R e w ) 0.07

2.2.3. Wavelet Analysis

Wavelet analysis is an effective means of signal processing. Wavelets are uneven and asymmetrical waveforms of adequately limited duration with an average value of zero. Wavelet analysis breaks up the mother wavelet signal into a shifted and scaled version. The typical wavelet transformation of a sine wave is presented in Figure 2.
In Fourier analysis, the signals are decomposed into different sine waves. Therefore, an irregular wavelet performs better than the steady sine wave for rapidly changing signals as it can provide better information about specific and appropriate locations. Wavelet analysis can also show any discontinuity, breakdown, trend, noise, coefficient, etc., of signals [48]. In this study, wavelet analysis was conducted using the MATLAB toolbox. There are two types of wavelet analyses, namely, discrete wavelet transform and continuous wavelet transform. There are various subgroups of these two types of wavelet transforms.

Contentious Wavelet Transform (CWT)

The continuous wavelet transforms (CWT) is a function of the shifted and scaled version of the wavelet function Ψ multiplied by the summation of the overall time of the signal. However, scaling means compressing or stretching the wavelets. A scale factor is used to represent the scaling, where the wavelet is more compressed when the scale factor is smaller. Wavelet sifting means hastening or detaining its onset.
C ( p o s i t i o n ,   s c a l e ) = f ( x ) Ψ ( p o s i t i o n , s c a l e , x ) d t
where C is the wavelet coefficient of the CWT as a function of position and scale [48].

Discrete Wavelet Transform (DWT)

The discrete wavelet transform is a wavelet transform where the wavelets are separately sampled. In this analysis, the original signal is divided into two parts, namely, approximations and details. The estimate a is the low pass filter where the low-frequency components of the original signal are separated, and the detail d is the high pass filter where high-frequency components pass. Moreover, the original signal x is not only isolated in one level, but the approximation a is also decomposed into much lower level (k = 3) components, known as multiple level decomposition, which is shown in Figure 3.
The significant difference between CWT and DWT is that CWT operates at every scale up to the maximum value, whereas, in DWT, the range and positions can be preselected and, in that way, the size of the analysis reduces its size and is thereby faster and more precise.
Mathematically, for j scales and k levels, the approximate information f a j ( x ) can be the summation of estimated coefficients a j , k and the scale function φ j , k ( x ) , as shown in Equation (13). Similarly, the detail information f d j ( x ) can also be described as the summation of approximate coefficients d j , k and the scale function Ψ j k ( x ) in Equation (14):
f a j ( x ) = k a j , k φ j , k ( x )
f d j ( x ) = k d j , k Ψ j k ( x )
A standard means of applying this is via logarithmic discretization of the scale s and then connecting it to the step size. The step size is the values between the translation parameter τ. The equation is adapted from Gao & Yan (2010) [49] and shown below, where Z is an integer parameter:
{ τ = k τ 0 s 0 j s = s 0 j   τ 0 0 ;   s 0 < 1   ( j ϵ Z , k ϵ Z )
Ψ j k ( x ) = s 0 0.5 j Ψ ( x s 0 j k τ 0 )  
Ψ j k ( x ) = 2 0.5 j Ψ ( x 2 j k )  
where j is the scale and k is the level of the wavelet. Equation (18) is the base wavelet equation. Addison (2017) [50] assumed s 0 = 2 and τ 0 = 1 ; therefore, Equation (19) can be derived and, finally, the discrete wavelet transform is obtained:
W ( j , k ) = [ f ( x ) , Ψ j k ( x ) ] = 2 0.5 j f ( x )   Ψ ( x 2 j k ) d x
f ( x ) = j , k C j , k Ψ j k ( x )
In Equation (18), f ( x ) is the original signal and, in Equation (19), C j , k is the wavelet coefficient. For multilevel wavelet analysis, there are many types of orthogonal wavelet transformations, which determine the shape of the wavelet. Among them, the Daubechies wavelet is one of the most common orthogonal wavelet transformations.

Daubechies Wavelet

The Daubechies wavelet uses scalar products with scaling wavelets and signals to calculate moving average and difference. This method allows for obtaining a good range of signal data to compute the median and deviation. Daub4 is the most accepted and straightforward way of analyzing wavelets. If we consider a signal x constituting n number of values, then the daub4 transformation creates the mapping x D k ( a k | d k ) to its approximation a k and details d k sub-signals for k levels:
a m = x . U m k
d m = x . Ψ m k
where each value of a m and d m are the scaler products. U m k is the scaling signal and Ψ m k is the wavelet at level k [51].

Wavelet Packet Analysis

In DWT, the primary signal is decomposed into an approximation and details, and the approximation is divided into a second-level approximation and features; thus, n levels of decomposition can be undertaken. In wavelet packet analysis, both the details and the estimate can be decomposed, which means the signal can be encoded in 2n ways. The wavelet packet decomposition tree is shown in Figure 4.
In the MATLAB toolbox, the entropy-based criterion is used to find the most desirable wavelet decomposition. Wavelet packet transformation provides many bases, and the best tree base can be found by the entropy criterion [48]. Wavelet packets are the general form of orthogonal wavelets. This splits detailed spaces to yield more prime decomposition. Coifman & Wickerhauser (1992) [52] explained the wavelet packet transformation equation as the following:
{ v 2 i j ( x ) = 2 k h k v i j ( 2 x k ) v 2 i + 1 j ( x ) = 2 k g k v i j ( 2 x k ) ;   i = 0 , 1 , 2 , a n d   k = 0 , 1 , m
where the two filters hk and gk are associated with the scaling function φ j k ( x ) and the base wavelet function Ψ j k ( x ) [49].

Wavelet Entropy

Wavelet entropy represents the non-uniformity of states, which is an ideal parameter for measuring the ordering of unsteady signals [53]. It can also provide information about the dynamic process and the signal potential. When the coefficient matrix of the wavelet transformation is represented by a probability distribution, the calculated wavelet entropy represents the randomness of the pattern [54]. The wavelet packet decomposition is an orthogonal function, which means the total energy entropy of the original signal should be the summation of the coefficient energy entropy [32]. The wavelet entropy energy can be defined as in Equation (23), where P i = E i / i = 1 n E i is the percentage of coefficient energy of the original signal [55]:
E N = i = 1 n P i l o g P i
In this study, the norm entropy was used to analyze the pressure signal. In an orthonormal basis entropy, s is the signal, s i j is the coefficient of s, and E is the entropy function such that E ( 0 ) = 0 and E ( s ) = i E ( s i ) . This entropy formula was used in MATLAB to calculate norm entropy. The concentration in l P is a norm where 1 P < 2 . Now, E ( s ) = | s i | P so E ( s ) = i | s | P = | | s | | P P for the norm entropy method [48]. The wavelet entropy can find small or abnormal frequencies, and wavelet entropy can thus find different characteristics of multiphase flow. Therefore, this study aims to characterize a two-phase flow pattern using norm entropy based on wavelet packet decomposition of the pressure signal. This method followed the steps shown in Figure 5.

Wavelet Packet Analysis of the Experimental Data

In this study, a pressure transducer provided time-domain pressure fluctuations, which were analyzed using wavelet packet analysis. As mentioned previously, the existing experimental setup only provides slug flow, and the dispersed bubble flow regime and pressure signal also show specific characteristics for each kind of flow regime. The data acquisition system collected pressure transducer signals with a sampling frequency of 100 Hz. Overall, 10,000 data points were used in 1-D wavelet analysis using MATLAB (details can be found in the above sections). The MATLAB processed data were taken into account for further analysis.

3. Results and Discussions

3.1. Flow Regime Mapping for Horizontal Flow

3.1.1. Air/Newtonian Flow Map

The experimental values were used to derive a flow regime map for horizontal pipe flow. This flow regime map was compared with that of Taitel & Dukler [7] and Mandhane et al. [9], where water and superficial air velocity was used. In the Taitel & Dukler [7] flow map for a horizontal pipe in Figure 6, most of the experimental data points fall in the respective flow regime area. However, the Taitel & Dukler [7] flow map better predicted the dispersed bubble flow for high gas/water velocity than the Mandhane et al. [9] flow map for this experimental setup.
In Figure 7, the Mandhane et al. [9] flow map is provided where the data for the slug and dispersed bubble flow data were fitted in the graph accordingly. The map can predict the slug and bubble flow regimes. Nevertheless, for high gas and water flow rates, this map cannot predict the dispersed bubble flow regime due to experimental conditions and the pipe diameter.

3.1.2. Air/non-Newtonian flow map

Researchers have also developed different flow pattern maps for horizontal, vertical, and inclined gas/non-Newtonian flow. For horizontal gas/non-Newtonian fluid, Chhabra & Richardson [56] developed a flow pattern map by slightly modifying the Mandhane et al. [9] horizontal flow pattern map (Figure 8). This map was developed for evaluating the literature and verified using 3700 data points of gas/non-Newtonian shear-thinning liquid mixture flow with 70% certainty. However, there was not enough data to confirm the Chhabra & Richardson [56] flow map for annular and slug flows.
In Figure 8, the experimental flow regime almost matches with the Chhabra & Richardson [10] flow map, but the slug to dispersed bubble flow transition started slightly earlier for this experiment. Chhabra & Richardson [10] used the particulate suspension of china clay, aqueous polymer solutions, limestone, and coal, which is a significantly more viscous shear-thinning non-Newtonian fluid compared to the 0.1% solution of xanthan gum that was used in this experiment. This is why the dispersed bubble flow regime started earlier. It is well-known that flow patterns of gas/non-Newtonian fluid are not significantly different from gas/Newtonian fluid for horizontal flow. However, due to high viscosity, the bubbles and slug could not break easily and collided together to form more prominent and well-defined bubbles.
The experimental data were used to reconstruct the flow maps and validated with the existing literature for the identification of the two-phase flow regimes for this experimental setup. Due to the targeted velocity condition pre-requisite for hole cleaning operations and the limitation of the existing flow loop system (that is, not designed for the extreme conditions of higher flow rates and high-pressure environments), the experiment was performed under moderate pressure and velocity conditions. Therefore, the results are effectively limited by the bubble flow and plug flow range, which can provide stratified, wavy, or annular flow. The Taitel & Dukler [7] and Mandhane et al. [9] flow maps for air/water two-phase horizontal flow and the Chhabra & Richardson [56] flow map for air/xanthan gum solution horizontal two-phase flow represented the flow regimes of the experimental setup relatively accurately. Nonetheless, the boundary of the flow regime varied due to the unpredictable characteristics of the transition of the flow pattern.

3.1.3. Air/Newtonian Two-phase Flow

The flow behavior of the air/Newtonian two-phase flow is shown in Figure 9. Figure 10 demonstrates that slug frequency increases with the increase of liquid superficial velocity for all studied combinations, while the superficial gas velocity was held constant for each set of data. Similar behavior was reported by Abed and Ghoben [22] for a Newtonian fluid (water). This happened due to the increase in liquid volume fraction. The liquid occupies more space in the liquid film region as the elongated bubble unit becomes smaller, which is why the slug unit increased in number.
Figure 10 shows the effect of a superficial gas ratio on slug frequency. Results from Figure 10 reveal that, for a constant liquid flow rate, slug frequency decreases with increasing gas velocity. The slug frequency initially decreases with the gas velocity up to around 5 m/s gas velocity, and then increases again, indicating additional slug formation. Moreover, to provide further insight, Figure 11 represents the slug frequency vs. the mixture velocity for air/water flow.
Comparing Figure 9 and Figure 11, it is evident that the slug frequency curves mainly depend on superficial gas velocity. Two of these graphs also show that in the range of 5 m/s to 6.5 m/s, the slug frequency reaches a minimum, and the slug frequency increases with increasing mixture velocity or gas superficial velocity. The sudden increase in slug frequency perhaps indicates the transition of the flow regime from slug to dispersed bubble flow. The higher gas flow rates induce more turbulence in the fluid flow, initiating the slug breakdown, which is attributed to the increase in the overall number of slugs in the system. It has also been observed that the number of dispersed bubbles increases in the slug pocket and wet film area. This indicates the beginning of the transition of the flow pattern. These findings also agree with the experimental data of Otten & Fayed [40] and Gregory & Scott [41]. Moreover, the slope of slug frequency vs. slug Froude number is reported in Figure 12.
In Figure 12, the slope of slug frequency vs. slug Froude number yields an equation of f s = 0.0673   (   N f r ) 0.098 . This equation shows a slight deviation from Gregory & Scott [41], which is shown in Equation (5), perhaps due to the experimental conditions and the assumption (vm = vs) for the air–water flow of this experiment. To evaluate the standard deviation and experimental accuracy of this study, Figure 13 reports the regression analysis of slug frequency and Froude number.
From Figure 13, the R2 value was found to be 88.1%, which indicates that the slug frequency and Froude number have a good relationship. Since the experimental and modeled data are close to the regression line, the equation f s = 0.0673   (   N f r ) 0.098 can explain the variability of the data around its mean.
In order to further validate and test the reliability of the obtained experimental data, the results are reported with correlations with different model results reported in various previous studies. Figure 14 shows the experimental and modeled validation with the Gregory & Scott [41] correlation model for air–water two-phase flow.
From Figure 14, it is observed that all the data points are close to the regression line, the regression has an R2 value of 73.8%, and all the experimental data fit well within in the 95% confidence interval. The experimental slug frequency results are also compared with the Zabaras et al. [23] predicted model for the air–water system in Figure 15. The results reveal that the data points of the Zabaras et al. [23] predicted model have comparatively deviated from the regression line, with an R2 value of 60.3% at a 95% confidence interval.
The experimental data and predictions of slug frequency of the Gregory & Scott [41] model have an R2 value of 73.8%, and the Zabaras et al. [23] model has an R2 value of 60%. Therefore, the Gregory & Scott [41] model is found to be closer to the experimental data. In the above graphs, the differences between experimental and predicted data varied mainly due to the variation in experimental conditions and setups [22]. Since the experimental data are more similar to the Gregory & Scott [41] model, Figure 16 shows the Gregory & Scott [41] correlation with the standard deviation of three samples under the same experimental condition. The error bar indicates a 95% confidence interval of the data, and none of the confidence intervals includes zero, which means the data are statistically significant and repeatable for air/water two-phase flow.

3.1.4. Air/non-Newtonian Two-Phase Flow

In this air/non-Newtonian fluid experiment, a 0.1% solution of xanthan gum was used as the non-Newtonian fluid. Non-Newtonian fluids are generally superior drilling fluids compared to Newtonian fluids due to their viscosity and the holding/uplifting capacity of the cuttings. Figure 17 illustrates the impact of liquid superficial velocity on slug frequency for the air/non-Newtonian flow system.
In Figure 17, the slug frequency increases with the increment of liquid non-Newtonian superficial velocity when the superficial gas velocity is held constant. Therefore, at lower superficial liquid velocity, the slug frequency increases sharply, and at higher liquid velocity, the slug frequency decreases. As the liquid velocity increases, the air requires more energy to drive the viscous fluid; nevertheless, the airflow rate is constant for each set. This explains why the number of slugs decreases as the liquid velocity increases at continuous airflow (see Figure 18). Furthermore, Figure 18 highlights the relationship between gas superficial velocity and slug frequency for air/non-Newtonian flow conditions.
In Figure 18, slug frequency change is shown with superficial gas velocity for a constant liquid superficial velocity. Here, 0.1% xanthan gum solution was used as a non-Newtonian fluid where the power-law index, n, is equal to 0.81 and k = 0.009344. From Figure 18 it is apparent that the slug frequency decreases as the gas velocity rises to 6.5 m/s. Because the gas flow rate increases in a constant liquid velocity, Taylor bubble formation becomes relatively more significant; therefore, the length of the slug unit increases, and slug frequency decreases. The evaluation of the slug frequency vs. mixture velocity for air/non-Newtonian fluid flow is reported in Figure 19.
Figure 19 represents the change of slug frequency with the mixture velocity of air–xanthan gum flow. In addition, it is seen that until 6.5 m/s, the mixture velocity slug frequency is minimal. Otten & Fayed [40] also obtained similar patterns for air/non-Newtonian flow. The same phenomena also occurred in Figure 11 for gas/Newtonian flow. However, the minimum slug frequency was around 5 m/s mixture velocity, which happened significantly earlier than the gas/non-Newtonian two-phase flow. Here, we can observe the inevitable effect of viscosity. Water viscosity at 20 °C room temperature is around 1 cP, and the experimental viscosity of 0.1% xanthan gum is 2.4 cP, which is slightly more viscous than water.
The flow mechanism of slug flow is that the gas bubble is trapped in the water and drives the water forward almost at the same velocity as the gas velocity. Nevertheless, when the liquid becomes viscous, the gas requires more energy to drive the liquid forward. At a constant airflow rate, it is hard to achieve extra energy, so the whole process becomes slow, and the slug velocity and number of slugs decrease [42]. If further experiments were conducted for gas/non-Newtonian fluids, slug frequency would increase again, with an increased gas flow rate in the slug to bubbly flow transition zone as the gas/water two-phase flow results in the transition from slug to dispersed bubble flow. The turbulence of the flow structure begins to increase, the unit slug starts to break down, and the number of slugs increases at higher gas flow rates. It has also been observed that the number of dispersed bubbles increases in the slug pocket and wet film area. This indicates the beginning of the transition of the flow pattern.
To understand this result further and demonstrate its wider applicability, the present study was extended to a non-dimensional number, i.e., the Froude number, presented in Figure 20 as the relationship between slug frequencies and Froude number for an air/xanthan gum solution.
In Figure 20, it is shown that the slope of slug frequency vs. slug Froude number for air/xanthan gum solution can be modeled using the equation v s = 0.0082   (   N f r ) 1.5624 , where the model strength R2 is 80.62%. In order to further validate the experimental data for an air/non-Newtonian two-phase flow system, Figure 21 and Figure 22 represent the relationship between the experimental and predicted data of the Gregory & Scott [41] proposed model with standard error bars and 95% confidence level, respectively.
The error bar in Figure 21 shows the 95% confidence interval of the data; none of the confidence intervals includes zero, which means the data are statistically significant and repeatable for air/water two-phase flow.
Picchi et al. [47] modified the slug frequency equation of Gregory & Scott [41] for a shear-thinning non-Newtonian fluid. Figure 21 represents the comparison of experimental results with the modified Gregory & Scott [41] slug frequency equation with error bars. The R2 value of 75.3% also represents the reliability and repeatability of the experimental data of this study. Figure 22 also shows the 95% confidence interval of the experimental data.

3.2. Wavelet Spectrum Analysis

The wavelet packet analysis decomposed the pressure signals into four levels. Among the wavelet decomposition methods, Daubechies’ four-scale base wavelet (db4) has been used most frequently in multiphase flow time series decompaction [57]. In this study, Daubechies’ four-scale base wavelet (db4) and norm entropy analysis method generated 16 wavelet packet coefficients. The pressure fluctuation signal achieved from the experimental data only yields a frequency of 100Hz. Thus, only four-level decomposition is sufficient because the pressure signals do not have sufficiently high-frequency and high-resolution data to obtain more detailed frequency analysis. The spectrum of the packet wavelet analysis represents the time-frequency plot, which provides a coefficient of decomposed frequencies at a different level. This spectrum represents the time and location of the fluctuation of the signal.
From Figure 23, it can be observed that the time-frequency plot divides the time-frequency plane into concentrated rectangles as a two-dimensional representation of signals. The pink color intensity of each box depends on the coefficient of the wavelet packet [34].
A similar plot is shown in Figure 24, which is the time-frequency plot for the bubbly flow regime for various flow conditions. The intensity of the pink shade represents the energy level of a time-frequency cell; the lower the energy content, the lighter the shade. For the bubbly flow regime, the bubbles are smaller, so the pressure fluctuation intensity is less, which means the low-frequency response has less energy content and the pink shade is light and scattered. In the slug flow, the Taylor bubble size is more prominent. Therefore, the low-frequency cells have more energy and a darker shade [34]. In addition, the repetition of the intense pink shade after a specific time interval can be evidence of the picks of the pressure signal. Thus, with a higher-resolution and better-quality sensor, where the pressure signal picks are more precise, this map can be a helpful way to understand the flow phenomena inside the pipe. When comparing the wavelet spectrum analysis of bubbly and slug flows for the same water flow rate, it has been observed that, for bubbly flow, the color intensity is comparatively less in the low-frequency response area. However, the use of a higher-resolution pressure transducer may enhance the wavelet spectrum quality with more precise fluctuation characteristics.

Wavelet Entropy Analysis

The wavelet entropy analysis of the pressure fluctuation data represents the nonlinearity of the gas/liquid two-phase flow. After calculating the wavelet entropy of 10,000 pressure signal data points of gas/liquid two-phase horizontal flow, it was seen that the wavelet entropy increased with the increase of the pressure signal fluctuation.
The entropy values of the pressure fluctuation data were compared with the gas to liquid volume flowrate ratio (GLR) and the void fraction ( α g = v g s v g s + v l s ) . The void fraction (αg) is the ratio of the gas velocity and the liquid velocity.
In Figure 25, it can also be observed that wavelet entropy increased with the increase of gas void fraction for gas/non-Newtonian fluid. This means the fluctuation of the pressure increases with the rise in the GLR for gas/water flow. Another observation is that the wavelet entropy value for gas/non-Newtonian flow is less than that for the gas/Newtonian fluid flow. This phenomenon occurred due to the viscosity effect of the non-Newtonian fluid, which means the disturbance of the pressure signals is less than that of the Newtonian fluid.
Figure 26 shows that with the growth of the GLR, the norm wavelet entropy increased, which means the fluctuation of the pressure increases with the increase of the GLR. However, the norm entropy change at the low GLR is not consistent. Fan et al. (2013) also saw similar behavior for low GLR values, mostly in the bubble flow or bubble-slug transition flow region. As bubble flow and motion are random and complicated, it is hard for low-resolution sensors or wavelet norm entropy to detect the change. Figure 27 and Figure 28 show the wavelet norm entropy of bubble, bubble-elongated bubble, and slug flow regimes at different GLR combinations for the experimental setup used in this study.
The wavelet norm entropy value may vary with higher resolution sensors. Wavelet norm entropy depends on the pressure fluctuation signal; therefore, the more sensitively and precisely the sensor can detect the flow condition, the more accurate the wavelet norm entropy will be. However, the wavelet entropy change pattern with a different type of flow regime should remain similar. Fan et al. (2013) and Sun et al. (2013) [32] both observed similar wavelet entropy change patterns. However, their sensors, experimental setups, and experimental conditions were different.

4. Conclusions

In this experimental work, air/Newtonian and air/non-Newtonian fluids were evaluated in horizontal flow loop conditions. The flow maps of both systems were established, and the type of flow regimes were determined for both air/Newtonian and air/non-Newtonian fluid flow conditions. The flow regime analysis revealed that the slug flow regime condition was observed for most of the studied conditions. The slug frequency analysis illustrated that the gas/Newtonian and gas/non-Newtonian fluids have a significant difference in slug properties, mainly due to the difference in viscosity of the systems. The air/water slug frequency decreased to an air velocity of approximately 5 m/s and then increased with increased air velocity. On the contrary, the air/xanthan gum solution did not demonstrate a similar effect, with the slug frequency slowly decreasing with increased air velocity within the experimental data range. Moreover, both the air/Newtonian and air/non-Newtonian slug frequency was validated with the slug frequency correlations available in the literature, and experimental data were found to be a satisfactory analog of model predictions. Although the pressure sensor used in this study was not able to provide high-frequency and high-resolution data, the wavelet decomposition and wavelet norm entropy were nonetheless calculated, and these offered recognizable flow characteristics for bubble, bubble-elongated bubble, and slug flow patterns. In addition, 1-D wavelet packet decomposition was shown to be a useful method to identify different features of multiphase flow and for recognizing different flow patterns. Therefore, the present study can pave the way for the development of safer designs in pipe flow, especially in oil and gas pipelines, where undesired slug flow can lead to corrosion issues.

Author Contributions

M.M., M.A.R. and S.I. conceived the presented idea. M.A.R. and S.I. developed the theory, and M.M. performed the experimentations. M.A.R. and S.I. supervised the findings of this work. M.S.K. and M.M. wrote the manuscript, and M.A.R. and S.I. reviewed the article. Overall, all the authors contributed to the betterment of the final manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The APC was partially funded by Qatar National Research Fund via the grant NPRP10-0101-170091.

Acknowledgments

The authors would like to acknowledge the Faculty of Engineering and Applied Science at the Memorial University of Newfoundland for providing laboratory facilities for this research. In addition, this publication was made possible by the grant NPRP10-0101-170091 from Qatar National Research Fund (a member of the Qatar Foundation). Statements made herein are solely the responsibility of the authors.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

v l s   Superficial Liquid Velocity
v l   Liquid Inlet Velocity
ρwWater Density = 996 Kg/M3
l p Pipe Length
v m No-Slip Mixture Velocity
fsSlug Frequency
fsnSlug Frequency for Non-Newtonian Fluid
N f s Froude Number
v m n No-Slip Mixture Velocity for Non-Newtonian Fluid
N f s n Froude Number for Non-Newtonian Fluid
RewWater Reynolds Number
ρnNon-Newtonian Density = 1002 Kg/M3
σShear Stress
xSignal
Ψ j k ( x ) Wavelet Base
kWavelet Level
nPower Law Index
v g s Superficial Gas Velocity
v g Gas Inlet Velocity
μwWater Viscosity 998 At 20 °C
dPipe Diameter
λLiquid Volume Fraction
v G True Average Gas Velocity in Multiphase Flow
gAcceleration due to gravity
v m 0 Minimum Slug Frequency in The Graph = 6 m/s
v l n   Non-Newtonian Liquid Inlet Velocity
RenNon-Newtonian Reynolds Number
γ y x . Shear Rate
μ n Apparent Viscosity
W ( j , k ) Wavelet Transform
φ j , k ( x ) Scale Function
jWavelet Scales
αgVoid fraction

References

  1. Epelle, E.I.; Gerogiorgis, D.I. A review of technological advances and open challenges for oil and gas drilling systems engineering. AIChE J. 2020, 66, e16842. [Google Scholar] [CrossRef]
  2. Zahid, A.A.; Rehman, S.R.; Rushd, S.; Hasan, A.; Rahman, M.A. Experimental investigation of multiphase flow behavior in drilling annuli using high speed visualization technique. Front. Energy 2018. [Google Scholar] [CrossRef]
  3. Khan, M.S.; Bavoh, C.B.; Partoon, B.; Nashed, O.; Lal, B.; Mellon, N.B. Impacts of ammonium based ionic liquids alkyl chain on thermodynamic hydrate inhibition for carbon dioxide rich binary gas. J. Mol. Liq. 2018, 261, 283–290. [Google Scholar] [CrossRef]
  4. Kamyab, M.; Rasouli, V.; Cavanough, G.; Mandal, S. Challenges of cuttings transport in micro- borehole coiled tubing drilling for mineral exploration. WIT Trans. Eng. Sci. 2013, 81, 109–120. [Google Scholar]
  5. Kelin, W.; Tie, Y.; Xiaofeng, S.; Shuai, S.; Shizhu, L. Review and analysis of cuttings transport in complex structural wells. Open Fuels Energy Sci. J. 2013, 6, 9–17. [Google Scholar] [CrossRef] [Green Version]
  6. Qasim, A.; Khan, M.S.; Lal, B.; Shariff, A.M. A perspective on dual purpose gas hydrate and corrosion inhibitors for flow assurance. J. Pet. Sci. Eng. 2019, 183, 106418. [Google Scholar] [CrossRef]
  7. Taitel, Y.; Dukler, A.E. A model for predicting flow regime transitions in horizontal and near horizontal gas-liquid flow. AIChE J. 1976, 22, 47–55. [Google Scholar] [CrossRef]
  8. Shojaei, M.J.; Rodríguez de Castro, A.; Méheust, Y.; Shokri, N. Dynamics of foam flow in a rock fracture: Effects of aperture variation on apparent shear viscosity and bubble morphology. J. Colloid Interface Sci. 2019, 552, 464–475. [Google Scholar] [CrossRef]
  9. Mandhane, J.M.; Gregory, G.A.; Aziz, K. A flow pattern map for gas-liquid flow in horizontal pipes. Int. J. Multiph. Flow 1974, 1, 537–553. [Google Scholar] [CrossRef]
  10. Chhabra, R.P.; Richardson, J.F. Prediction of flow pattern for the co-current flow of gas and non-newtonian liquid in horizontal pipes. Can. J. Chem. Eng. 1984, 62, 449–454. [Google Scholar] [CrossRef]
  11. Azarinezhad-Mohammadi, R. A Chemical Based Wet Cold Flow Approach for Addressing Hydrate Flow Assurance Problems. Ph.D. Thesis, Heriot-Watt University , Edinburgh, UK, 2010. [Google Scholar]
  12. Joshi, S.V.; Grasso, G.A.; Lafond, P.G.; Rao, I.; Webb, E.; Zerpa, L.E.; Sloan, E.D.; Koh, C.A.; Sum, A.K. Experimental flowloop investigations of gas hydrate formation in high water cut systems. Chem. Eng. Sci. 2013, 97, 198–209. [Google Scholar] [CrossRef]
  13. Khan, M.S.; Lal, B.; Bavoh, C.B.; Keong, L.K.; Bustam, A. Influence of ammonium based compounds for gas hydrate mitigation: A short review. Indian J. Sci. Technol. 2017, 10, 1–6. [Google Scholar] [CrossRef]
  14. Al-Safran, E. Investigation and prediction of slug frequency in gas/liquid horizontal pipe flow. J. Pet. Sci. Eng. 2009, 69, 143–155. [Google Scholar] [CrossRef]
  15. Sultan, R.A.; Alfarek, S.; Rahman, M.A.; Zendehboudi, S. CFD and experimental approach on three phase gas-liquid-solid Newtonian fluid flow in horizontal pipes. Int. J. Comput. Methods Exp. Meas. 2018, 7, 33–44. [Google Scholar] [CrossRef]
  16. Zadrazil, I.; Markides, C.N. An experimental characterization of liquid films in downwards co-current gas-liquid annular flow by particle image and tracking velocimetry. Int. J. Multiph. Flow 2014, 67, 42–53. [Google Scholar] [CrossRef] [Green Version]
  17. Cerqueira, R.F.L.; Paladino, E.E.; Ynumaru, B.K.; Maliska, C.R. Image processing techniques for the measurement of two-phase bubbly pipe flows using particle image and tracking velocimetry (PIV/PTV). Chem. Eng. Sci. 2018, 189, 1–23. [Google Scholar] [CrossRef]
  18. Talley, J.D.; Worosz, T.; Kim, S.; Buchanan, J.R. Characterization of horizontal air–water two-phase flow in a round pipe part I: Flow visualization. Int. J. Multiph. Flow 2015, 76, 212–222. [Google Scholar] [CrossRef] [Green Version]
  19. Majumder, S.K.; Ghosh, S.; Kundu, G.; Mitra, A.K. Frictional pressure drop of gas-Newtonian and gas-non Newtonian slug flow in vertical pipe. Int. J. Chem. React. Eng. 2011. [Google Scholar] [CrossRef]
  20. Ehinmowo, A.B.; Cao, Y.; Yeung, H.C. Slug flow in large diameter pipeline—Riser systems: Prediction and mitigation. Covenant J. Eng. Technol. 2018, 2, 1–21. [Google Scholar]
  21. Majumder, S.K.; Ghosh, S.; Mitra, A.K.; Kundu, G. Gas-Newtonian and gas-non-Newtonian slug flow in vertical pipe, part I: Gas holdup characteristics. Int. J. Chem. React. Eng. 2010. [Google Scholar] [CrossRef]
  22. Abed, E.M.; Ghoben, K. Gas-liquid slug frequency and slug unit length in horizontal pipes. Iraqi J. Mech. Mater. Eng. 2015, 15, 166–180. [Google Scholar]
  23. Zabaras, G.J.; Shell, E.; Co, P.T. Prediction of Slug Frequency for Gas/Liquid Flows; Society of Petroleum Engineers: Richardson, TX, USA, 1999. [Google Scholar]
  24. Sultan, R.A.; Rahman, M.A.; Rushd, S.; Zendehboudi, S.; Kelessidis, V.C. CFD analysis of pressure losses and deposition velocities in horizontal annuli. Int. J. Chem. Eng. 2019, 2019, 1–17. [Google Scholar] [CrossRef] [Green Version]
  25. Daza Gómez, M.A.; Ratkovich, N.; Pereyra, E.; Poesio, P. CFD simulations of two-phase gas/non-Newtonian fluid flows in pipes. In Proceedings of the 3rd World Congress on Momentum, Heat and Mass Transfer (MHMT’18), Budapest, Hungary, 12–14 April 2018. [Google Scholar]
  26. Dziubinski, M.; Fidos, H.; Sosno, M. The flow pattern map of a two-phase non-Newtonian liquid-gas flow in the vertical pipe. Int. J. Multiph. Flow 2004, 30, 551–563. [Google Scholar] [CrossRef]
  27. Katzbauer, B. Properties and applications of xanthan gum. Polym. Degrad. Stab. 1998, 59, 81–84. [Google Scholar] [CrossRef]
  28. Bavoh, C.B.; Md Yuha, Y.B.; Tay, W.H.; Ofei, T.N.; Lal, B.; Mukhtar, H. Experimental and modelling of the impact of quaternary ammonium salts/ionic liquid on the rheological and hydrate inhibition properties of xanthan gum water-based muds for drilling gas hydrate-bearing rocks. J. Pet. Sci. Eng. 2019, 183, 106468. [Google Scholar] [CrossRef]
  29. Owusu, E.B.; Tsegab, H.; Sum, C.W.; Padmanabhan, E.; Qasim, A.; Khan, M.S.; Lal, B.; Ismail, M.C.; Rostani, K.; Shojaei, M.J.; et al. Synthesis, characterization, stability and thermal conductivity of multi-walled carbon nanotubes (MWCNTs) and eco-friendly jatropha seed oil based nanofluid: An experimental investigation and modeling approach. J. Mol. Liq. 2019, 293, 6–14. [Google Scholar]
  30. Tutu, N.K. Pressure fluctuations and flow pattern recognition in vertical two phase gas-liquid flows. Int. J. Multiph. Flow 1982, 8, 443–447. [Google Scholar] [CrossRef]
  31. Drahoš, J.; Čermák, J.; Selucky, K.; Ebner, L. Characterization of hydrodynamic regimes in horizontal two-phase flow: Part II: Analysis of wall pressure fluctuations. Chem. Eng. Process. Process Intensif. 1987, 22, 45–52. [Google Scholar] [CrossRef]
  32. Sun, Z.; Shao, S.; Gong, H. Gas-liquid flow pattern recognition based on wavelet packet energy entropy of vortex-induced pressure fluctuation. Meas. Sci. Rev. 2013, 13, 83–88. [Google Scholar] [CrossRef] [Green Version]
  33. Blaney, S. Gamma Radiation Methods for Clamp-On Multiphase Flow Metering; Cranfield University: Bedford, UK, 2008. [Google Scholar]
  34. Park, S.H.; Kim, S.D. Characterization of pressure signals in a bubble column by wavelet packet transform. Korean J. Chem. Eng. 2003, 20, 128–132. [Google Scholar] [CrossRef]
  35. De Fang, L.; Liu, R.; Lu, Q.H.; Wang, X.J.; Liang, Y.J. The flow pattern transition identification and interphases force detection of gas-liquid two-phase flow. AASRI Procedia 2012, 3, 534–539. [Google Scholar] [CrossRef]
  36. Elperin, T.; Klochko, M. Flow regime identification in a two-phase flow using wavelet transform. Exp. Fluids 2002, 32, 674–682. [Google Scholar] [CrossRef]
  37. Kouba, G.E.; Jepson, W.P. The flow of slugs in horizontal, two-phase pipelines. J. Energy Resour. Technol. 1990, 112, 20–24. [Google Scholar] [CrossRef]
  38. Hubbard, M.G. An Analysis of Horizontal Gas-Liquid Slug Flow; University of Houston: Houston, TX, USA, 1965. [Google Scholar]
  39. Hubbard, M.G.; Dukler, A.E. The characterization of flow regimes for horizontal two-phase flow. Proc. Heat Transf. Fluid 1966, 1996, 100–121. [Google Scholar]
  40. Otten, L.; Fayed, A.S. Slug velocity and slug frequency measurements in concurrent air-non-Newtonian slug flow. Trans. Inst. Chem. Eng. 1977, 55, 64–67. [Google Scholar]
  41. Gregory, G.A.; Scott, D.S. Correlation of liquid slug velocity and frequency in horizontal concurrent gas-liquid slug flow. AIChE J. 1969, 15, 933–935. [Google Scholar] [CrossRef]
  42. Rosehart, R.G.; Rhodes, E.; Scott, D.S. Studies of gas liquid (non-Newtonian) slug flow: Void fraction meter, void fraction and slug characteristics. Chem. Eng. J. 1975, 10, 57–64. [Google Scholar] [CrossRef]
  43. Greskovich, E.J.; Shrier, A.L. Slug frequency in horizontal gas-liquid slug flow. Ind. Eng. Chem. Process Des. Dev. 1972, 11, 317–318. [Google Scholar] [CrossRef]
  44. Zabaras, G.J. Prediction of slug frequency for gas-liquid flows. In Proceedings of the SPE Annual Technical Conference and Exhibition, Houston, TX, USA, 3–6 October 1999. [Google Scholar]
  45. Heywood, N.I.; Richardson, J.F. Slug flow of air water mixtures in a horizontal pipe: Determination of liquid holdup by γ-ray absorption. Chem. Eng. Sci. 1979, 34, 17–30. [Google Scholar] [CrossRef]
  46. Shea, R.; Eidsmoen, H.; Nordsveen, M.; Rasmussen, J.; Xu, Z.; Nossen, J. Slug frequency prediction method comparison. In Proceedings of the 4th North American Conference on Multiphase Technology, Banff, AB, Canada, 3–4 June 2004; pp. 227–237. [Google Scholar]
  47. Picchi, D.; Manerba, Y.; Correra, S.; Margarone, M.; Poesio, P. Gas/shear-thinning liquid flows through pipes: Modeling and experiments. Int. J. Multiph. Flow 2015, 73, 217–226. [Google Scholar] [CrossRef]
  48. Misiti, M.; Misiti, Y.; Oppenheim, G.; Poggi, J.-M. Wavelet Toolbox; MathWorks Inc.: Natick, MA, USA, 1996; Volume 15, p. 21. [Google Scholar]
  49. Gao, R.X.; Yan, R. Wavelets: Theory and Applications for Manufacturing; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  50. Addison, P.S. The Illustrated Wavelet Transform Handbook: Introductory Theory and Applications in Science, Engineering, Medicine and Finance; CRC Press: Boca Raton, FL, USA, 2017. [Google Scholar]
  51. Walker, J.S. A Primer on Wavelets and Their Scientific Applications; CRC Press: Boca Raton, FL, USA, 2008. [Google Scholar]
  52. Coifman, R.R.; Wickerhauser, M.V. Entropy-based algorithms for best basis selection. IEEE Trans. Inf. Theory 1992, 38, 713–718. [Google Scholar] [CrossRef] [Green Version]
  53. Uyar, M.; Yildirim, S.; Gencoglu, M.T. An effective wavelet-based feature extraction method for classification of power quality disturbance signals. Electr. Power Syst. Res. 2008, 78, 1747–1755. [Google Scholar] [CrossRef]
  54. Fan, C.; Ding, Y.; Ren, X. Wavelet entropy applied in gas-liquid two-phase flow. In Proceedings of the 32nd Chinese Control Conference (CCC), Xi’an, China, 26–28 July 2013; pp. 8623–8627. [Google Scholar]
  55. Yu, Y.; Junsheng, C. A roller bearing fault diagnosis method based on EMD energy entropy and ANN. J. Sound Vib. 2006, 294, 269–277. [Google Scholar] [CrossRef]
  56. Chhabra, R.P.; Richardson, J.F. Non-Newtonian Flow: Fundamentals and Engineering Applications; Butterworth-Heinemann: Oxford, UK, 1999. [Google Scholar]
  57. Shaikh, A.; Al-Dahhan, M.H. A review on flow regime transition in bubble columns. Int. J. Chem. React. Eng. 2007. [Google Scholar] [CrossRef]
Figure 1. Schematic of the experimental setup (horizontal test section) with a data acquisition system.
Figure 1. Schematic of the experimental setup (horizontal test section) with a data acquisition system.
Applsci 10 03272 g001
Figure 2. Wavelet transformation of a sine wave.
Figure 2. Wavelet transformation of a sine wave.
Applsci 10 03272 g002
Figure 3. Multiple level discrete wavelet analysis.
Figure 3. Multiple level discrete wavelet analysis.
Applsci 10 03272 g003
Figure 4. Wavelet packet analysis decomposition tree.
Figure 4. Wavelet packet analysis decomposition tree.
Applsci 10 03272 g004
Figure 5. The steps of wavelet decomposition for different flow pattern identification.
Figure 5. The steps of wavelet decomposition for different flow pattern identification.
Applsci 10 03272 g005
Figure 6. Comparison of the Taitel & Dukler [7] (adapted) flow map with experimental data for horizontal gas/Newtonian flow. Copyright AIChE Journal, 2020.
Figure 6. Comparison of the Taitel & Dukler [7] (adapted) flow map with experimental data for horizontal gas/Newtonian flow. Copyright AIChE Journal, 2020.
Applsci 10 03272 g006
Figure 7. Comparison of the Mandhane et al. [9] (adapted) flow regime map with experimental data obtained for horizontal gas/Newtonian flow. Copyright Elsevier, 1974.
Figure 7. Comparison of the Mandhane et al. [9] (adapted) flow regime map with experimental data obtained for horizontal gas/Newtonian flow. Copyright Elsevier, 1974.
Applsci 10 03272 g007
Figure 8. Comparison of the Chhabra & Richardson [10] (adapted) flow regime map with experimental data obtained for horizontal gas/non-Newtonian flow. Copyright Canadian Journal of Chemical Engineering, 2020.
Figure 8. Comparison of the Chhabra & Richardson [10] (adapted) flow regime map with experimental data obtained for horizontal gas/non-Newtonian flow. Copyright Canadian Journal of Chemical Engineering, 2020.
Applsci 10 03272 g008
Figure 9. Effect of liquid superficial velocity on slug frequency for air/water flow.
Figure 9. Effect of liquid superficial velocity on slug frequency for air/water flow.
Applsci 10 03272 g009
Figure 10. Effect of gas superficial velocity on slug frequency of air/water two-phase flow.
Figure 10. Effect of gas superficial velocity on slug frequency of air/water two-phase flow.
Applsci 10 03272 g010
Figure 11. Slug frequency vs. mixture velocity for air/water flow.
Figure 11. Slug frequency vs. mixture velocity for air/water flow.
Applsci 10 03272 g011
Figure 12. Slug frequency vs. Froude number for air/water flow.
Figure 12. Slug frequency vs. Froude number for air/water flow.
Applsci 10 03272 g012
Figure 13. Relationship between slug frequency and Froude number for air–water system data (R2 = 88.1%).
Figure 13. Relationship between slug frequency and Froude number for air–water system data (R2 = 88.1%).
Applsci 10 03272 g013
Figure 14. Experimental slug frequency for the air–water system compared with the prediction model of Gregory & Scott [41] correlation (R2 = 73.8%).
Figure 14. Experimental slug frequency for the air–water system compared with the prediction model of Gregory & Scott [41] correlation (R2 = 73.8%).
Applsci 10 03272 g014
Figure 15. Experimental slug frequency for the air–water system compared with the prediction model of Zabaras et al. [23] correlation (R2 = 60%).
Figure 15. Experimental slug frequency for the air–water system compared with the prediction model of Zabaras et al. [23] correlation (R2 = 60%).
Applsci 10 03272 g015
Figure 16. Experimental slug frequency for the air–water system compared to the predictions of Gregory & Scott [41] correlation with standard deviations.
Figure 16. Experimental slug frequency for the air–water system compared to the predictions of Gregory & Scott [41] correlation with standard deviations.
Applsci 10 03272 g016
Figure 17. Effect of liquid superficial velocity on slug frequency for air/non-Newtonian flow.
Figure 17. Effect of liquid superficial velocity on slug frequency for air/non-Newtonian flow.
Applsci 10 03272 g017
Figure 18. Effect of gas superficial velocity on slug frequency for air/non-Newtonian flow.
Figure 18. Effect of gas superficial velocity on slug frequency for air/non-Newtonian flow.
Applsci 10 03272 g018
Figure 19. Slug frequency vs. mixture velocity for air/non-Newtonian fluid flow.
Figure 19. Slug frequency vs. mixture velocity for air/non-Newtonian fluid flow.
Applsci 10 03272 g019
Figure 20. Slug frequency vs. Froude number for an air/xanthan gum solution.
Figure 20. Slug frequency vs. Froude number for an air/xanthan gum solution.
Applsci 10 03272 g020
Figure 21. Experimental vs. predicted data [41] with standard error bars.
Figure 21. Experimental vs. predicted data [41] with standard error bars.
Applsci 10 03272 g021
Figure 22. Experimental slug frequency for the air–xanthan gum system compared to predictions by the Gregory & Scott [41] correlation, where R2 = 74.6%.
Figure 22. Experimental slug frequency for the air–xanthan gum system compared to predictions by the Gregory & Scott [41] correlation, where R2 = 74.6%.
Applsci 10 03272 g022
Figure 23. Spectrum for slug flow under different flow conditions.
Figure 23. Spectrum for slug flow under different flow conditions.
Applsci 10 03272 g023
Figure 24. Spectrum for bubbly flow under different flow conditions.
Figure 24. Spectrum for bubbly flow under different flow conditions.
Applsci 10 03272 g024
Figure 25. Change of wavelet entropy with gas void fraction for gas/Newtonian fluid.
Figure 25. Change of wavelet entropy with gas void fraction for gas/Newtonian fluid.
Applsci 10 03272 g025
Figure 26. Change of wavelet entropy with gas to liquid ratio for gas/Newtonian flow.
Figure 26. Change of wavelet entropy with gas to liquid ratio for gas/Newtonian flow.
Applsci 10 03272 g026
Figure 27. Wavelet entropy flow map for gas/Newtonian flow.
Figure 27. Wavelet entropy flow map for gas/Newtonian flow.
Applsci 10 03272 g027
Figure 28. Wavelet entropy flow map for gas/non-Newtonian flow.
Figure 28. Wavelet entropy flow map for gas/non-Newtonian flow.
Applsci 10 03272 g028
Table 1. The chemicals and experimental parameters used in this study.
Table 1. The chemicals and experimental parameters used in this study.
Newtonian FluidWater
Non-Newtonian Fluid0.1% xanthan gum solution (Power Law Index, n = 0.81) and (Power Law Index, m or k = 0.009344)
Liquid velocity range1.5 to 2.5 m/s
Air velocity2.9 to 6.4 m/s

Share and Cite

MDPI and ACS Style

Morshed, M.; Khan, M.S.; Rahman, M.A.; Imtiaz, S. Flow Regime, Slug Frequency and Wavelet Analysis of Air/Newtonian and Air/non-Newtonian Two-Phase Flow. Appl. Sci. 2020, 10, 3272. https://doi.org/10.3390/app10093272

AMA Style

Morshed M, Khan MS, Rahman MA, Imtiaz S. Flow Regime, Slug Frequency and Wavelet Analysis of Air/Newtonian and Air/non-Newtonian Two-Phase Flow. Applied Sciences. 2020; 10(9):3272. https://doi.org/10.3390/app10093272

Chicago/Turabian Style

Morshed, Munzarin, Muhammad Saad Khan, Mohammad Azizur Rahman, and Syed Imtiaz. 2020. "Flow Regime, Slug Frequency and Wavelet Analysis of Air/Newtonian and Air/non-Newtonian Two-Phase Flow" Applied Sciences 10, no. 9: 3272. https://doi.org/10.3390/app10093272

APA Style

Morshed, M., Khan, M. S., Rahman, M. A., & Imtiaz, S. (2020). Flow Regime, Slug Frequency and Wavelet Analysis of Air/Newtonian and Air/non-Newtonian Two-Phase Flow. Applied Sciences, 10(9), 3272. https://doi.org/10.3390/app10093272

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop