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

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.


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]. 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.

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.

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.

Newtonian Fluid Water
Non-Newtonian 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].

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 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].

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: Hubbard & Dukler [39] also predicted the actual average gas velocity as in Equation (2), which also agreed with other experimental data [41].
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: 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] Here, v 0 m 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: 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: where v m and v s l 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: 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. 6 5 0.836 + 2.7sin 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 s l / v s l + v s g , and d is the pipe diameter.
Shea et al. [46] developed a correlation describing slug frequency as a function of pipe length. In Equation (10), v S l 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 l p , 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: 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.

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.
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 = ( + ) ⁄ , and d is the pipe diameter.
Shea et al. [46] developed a correlation describing slug frequency as a function of pipe length. In Equation (10), is the superficial liquid velocity, d is the pipe diameter, and 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: Picchi et al. [47] described a slug frequency equation that considers the rheology of the shearthinning fluid. This equation is a modified version of the Gregory & Scott [41] correlation. In Equation

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 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.
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.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 7 of 31 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.
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 ( ) can be the summation of estimated coefficients , and the scale function , ( ), as shown in Equation (13). Similarly, the detail information ( ) can also be described as the summation of approximate coefficients , and the scale function ( ) in Equation (14): 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 j a (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 j d (x) can also be described as the summation of approximate coefficients d j,k and the scale function Ψ jk (x) in Equation (14): Appl. Sci. 2020, 10, 3272 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: 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: 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: where each value of a m and d m are the scaler products. U k m is the scaling signal and Ψ k m 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 2 n 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: where the two filters hk and gk are associated with the scaling function ( ) and the base wavelet function ( ) [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 = / ∑ is the percentage of coefficient energy of the original signal [55]: In this study, the norm entropy was used to analyze the pressure signal. In an orthonormal basis entropy, s is the signal, is the coefficient of s, and E is the entropy function such that (0) = 0 and ( ) = ∑ ( ) . This entropy formula was used in MATLAB to calculate norm entropy. The concentration in is a norm where 1 ≤ < 2. Now, ( ) = | | so ( ) = ∑ | | = | | 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. 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: where the two filters h k and g k are associated with the scaling function ϕ jk (x) and the base wavelet function Ψ jk (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 / n i=1 E i is the percentage of coefficient energy of the original signal [55]: In this study, the norm entropy was used to analyze the pressure signal. In an orthonormal basis entropy, s is the signal, s ij 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.

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.

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.

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. 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.

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.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 12 of 28 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 highpressure 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.

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 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.

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.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 13 of 28 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 Figures 9 and 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.
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 break down, 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.  [41]. Moreover, the slope of slug frequency vs. slug Froude number is reported in Figure 12. . 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 airwater 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. . 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 (v m = v s ) 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.
In Figure 12, the slope of slug frequency vs. slug Froude number yields an equation of = 0.0673 . . 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 airwater 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 R 2 value was found to be 88.1%, which indicates that the slug frequency and 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.  Figure 13. Relationship between slug frequency and Froude number for air-water system data (R 2 = 88.1%).
From Figure 13, the R 2 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 R 2 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 From Figure 14, it is observed that all the data points are close to the regression line, the regression has an R 2 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 R 2 value of 60.3% at a 95% confidence interval.
From Figure 14, it is observed that all the data points are close to the regression line, the regression has an R 2 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 R 2 value of 60.3% at a 95% confidence interval.  [23] model has an R 2 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. The experimental data and predictions of slug frequency of the Gregory & Scott [41] model have an R 2 value of 73.8%, and the Zabaras et al. [23] model has an R 2 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.

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.

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 0 0  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. Appl. Sci. 2020, 10, x FOR PEER REVIEW 18 of 28 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. Figure 18: Effect of gas superficial velocity on slug frequency for air/non-Newtonian flow. 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. 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 18: Effect of gas superficial velocity on slug frequency for air/non-Newtonian flow. 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 R 2 is 80.62%. In order to further validate the experimental data for an air/non-Newtonian two-phase flow system, Figures 21 and 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. 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. 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. 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.   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 R 2 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.

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].
Appl. Sci. 2020, 10, x FOR PEER REVIEW 22 of 31 modified Gregory & Scott [41] slug frequency equation with error bars. The R 2 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.

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- 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.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 23 of 31 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 ( = ). 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

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 s g v s g +v s l ). 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.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 22 of 28 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 ( = ). The void fraction ( g) is the ratio of the gas velocity and the liquid velocity. Figure 25: Change of wavelet entropy with gas void fraction for gas/Newtonian fluid. 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, bubbleelongated bubble, and slug flow regimes at different GLR combinations for the experimental setup used in this study.  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. Figures 27 and 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.

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 highresolution 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.  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.

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.