Fluid Dynamics & Materials Processing

DOI: 10.32604/fdmp.2022.019532

ARTICLE

Systematic Investigation of Two-Phase Flow in Special Channels

Oleg R. Kuzichkin1,*, Muneam Hussein Ali2, Ayad F. Alkaim3, Lis M. Yapanto4, Gabdrakhman H. Valiev5, L. S. Abdullah6 and Mustafa M. Kadhim7

1Department of Information and Robototechnic Systems, Belgorod State University, Belgorod, 308015, Russia
2Al-Nisour University College, Baghdad, 10012, Iraq
3College of Science for Women, University of Babylon, Babil, 51002, Iraq
4Departement of Aquatic Management, Faculty of Fisheries and Marine Science Universitas Negeri Gorontalo, Gorontalo, 96128, Indonesia
5Kazan Federal University, Kazan, 420008, Russia
6Lipetsk State Technical University, Lipetsk, 398600, Russia
7Department of Dentistry, Kut University College, Kut, Wasit, 52001, Iraq
*Corresponding Author: Oleg R. Kuzichkin. Email: golubeva.julivgt@inbox.ru
Received: 28 September 2021; Accepted: 28 December 2021

1  Introduction

Gas-liquid two-phase flows have many applications in the oil and gas, petrochemical, nuclear technology, and other industries. One of the main industrial applications is two-phase currents in pipes, which are used in different slopes from horizontal to vertical. One of the most important natural and unintended phenomena of this interaction is the mixing of water passing through the structures and the surrounding airflow, which in some cases causes severe reactions between the existing climate and causes considerable damage. This type of current is called “two-phase current” [1–4]. One of the most important problems in the analysis of two-phase current in pipes is to determine the exact pattern of flow or distribution of liquid and gas phase inside the pipe [5]. Since many important design and engineering parameters such as pressure drop, mass transfer, heat transfer, etc., are closely related to the distribution of phases, so determining this distribution or in other words determining two-phase flow patterns is one of the essential requirements for two-phase flow analysis. Two-phase flow follows all the rules governing fluid mechanics. Two-phase current systems are highly complex (even one-dimensional two-phase currents in ducts) and many characteristics make these types of currents indeterminate and complex compared to single-phase currents [6]. To simplify these complexities, various researchers have conducted model experiments to use laboratory and experimental results to obtain simple equations for the analysis of practical and engineering problems, two-phase climate flows in various problems, and Various structures occur-including chemical processes, energy generation and formation, water collection and transmission systems, the oil industry, and offshore hinge lines for oil transportation [7].

It should be noted that despite many attempts to classify the types of biphasic flow regimes, all of these methods are highly qualitative and often conform to the personal views of different researchers. So far, different regimes have been defined, and a wide range of titles have been used for this purpose. These patterns are flow, layer or mould flow, layer flow, wave flow, mixing flow, annular flow, and droplet flow. Methods for predicting two-phase flow patterns within flow transmission pipelines are generally divided into two groups: using experimental and laboratory facilities and using analytical models [8]. For many years, the prediction of various two-phase air-water flow regimes in pipelines has been considered through experimental experiments. A common method for this purpose is to collect experimental data on the flow rates of the gas and liquid phases and the physical properties of each phase on a section of pipeline under test. Many studies have now been conducted in the fields of experimental, kinematic, and dynamic two-phase flow modelling, but many aspects of this phenomenon remain unknown [9–11]. Usually, the necessary research is done to identify the types of biphasic flow patterns in the channel and tubular lines of biphasic current transmission using direct current observation method. Among the methods that make these observations possible are “very fast photography” and “X-ray audiography”; So that using the received images, the type of current available can be interpreted to some extent.

The most important possible patterns in horizontal and near-horizontal ducts (tunnels and closed ducts) are floc, wave, and stratigraphic regimes. Especially when the duct length is long, other diets occur in a very small length and very quickly become the main diets mentioned with a limited length transfer, and even if they appear, they generally have many similarities to the three mentioned diets. Typically in large sloping canals, the flow velocity in operation is more than 5 to 10 m/s, and the Reynolds number is in the range of 107 to 109 and surface aeration is almost always on these Structures are observed. Due to the existence of various equations and parameters governing such high velocity and aeration flows, their numerical analysis is not possible, and in most cases to solve many problems related to hydraulics, such flows have to use hypotheses and simplifications. We are the ones who make the results differ from reality [12].

In [13], a three-dimensional Eulerian method was presented for simulating dynamic systems comprising multiple compressible solid and fluid components where internal boundaries are tracked using level-set functions. Due to the observed geometric similarities, many researchers consider these regimes as one of the three main regimes mentioned. The stratigraphic regime occurs when water or liquid phase at the bottom of a pipe or channel and close to the floor, and air with The gas-phase moves in the upper part of the liquid stream and close to the ceiling and two separate layers are formed [14]. Wave regime occurs when the airflow increases relative to the stratigraphic flow and alternating and uniform waves form and move along the tube. As the flow rate or velocity of the air increases relative to the wave current, the waves touch the roof of the tunnel under certain conditions, causing the air layer to rupture, resulting in a clotting regime. An important and significant issue in biphasic currents is that these types of currents have an oscillating and time-dependent nature. Depending on the flow regime, these oscillating characteristics change. This oscillating nature of the current causes special reactions in the biphasic current. In this paper, we try to study the kinematics of two-phase current in open channels using methods [15].

2  Methods

The model used in this research, according to Fig. 1 is a combination of the main water transmission channel with a diameter of 0.09 and a length of 10 meters. The upstream tank is a metal tank in the shape of a rectangular cube, which is installed in the middle of its longitudinal wall at the height of 5 cm from the bottom of a hole with a diameter of 1.5 m in order to connect the duct to the tank [16].

Figure 1: A general overview of parts of the laboratory model

Governing equations for above presented flow include continuity, momentum and energy [17]:

∇.U→=0; ∂u∂X+∂V∂R=0, (1)

u∂u∂X+V∂u∂R=−∂P∂X+1Renf[∂∂X(∂u∂X)n+∂∂R(∂u∂R)]nu∂V∂X+V∂V∂R=−∂P∂X+1Renf[[∂∂X(∂V∂X)n+1R∂∂RR(∂V∂R)]n−V2R (2)

U→.∇→T=(knf/ρnfCpnf) ∇2T;

u∂T∂X+V∂T∂R=1Renf.Prnf[∂∂X(∂T∂X)n+1R∂∂RR(∂T∂R)]n (3)

where, P is the pressure, U→ is the velocity vector, T is temperature, and ρnf, knf, Cpnf are the density, thermal conductivity, and specific heat of nano-fluid, respectively. For the pressure corrector can be used the PISO algorithm which described in the [18]. The energy conservation equation consists of several terms, which are written separately here. Relationships governing two-phase currents Equations governing two-phase currents, like single-phase currents, include the equations of mass, momentum and energy [19]. Energy changes over time and the rate of energy input and output of the two-phase control volume is the internal energy of each phase. Rate of work on two-phase control volume:

∂∂t[αkρk(ek+uk22)Aδz]+mk(ek+uk22)δz−[mk(ek+uk22)−δz∂∂zmk(ek+uk22)] (4)

where e is the internal energy of each phase. Also, the rate of work on two-phase control volume is calculated by:

[mkpρk−(mkpρk+δz∂∂z(mkpρk))Aδz]−mkg.sinθδz−pAδz∂αk∂z+pkδzpρk+uk∑1n⁡τknpknδz (5)

Energy rate added to each phase by mass transfer:

Pkδz=(ek+uk22) (6)

Finally, the energy conservation equation for each phase of the steady stream without the presence of a heat source is written in the form of Eqs. (7) and (8):

d[mg(ig+ug22)]+mgg.sinθ.δz=qwgPwg+qgIPgIδz+ugτgIPgI+pgδz(ig+ug22) (7)

d[mI(iI+uI22)]+mIg.sinθ.δz=qwIPwI+qlgPlgδz+uIτIgPIg+plδz(il+ul22)

ik=uk+pρk (8)

It should be noted that the phrase of Eq. (9) indicates the input energy rate to each phase:

qkwPkwδz+∑1n⁡qknPknδz+qkAαkδz (9)

Laboratory flow conditions In order to evaluate the turbulence characteristics of the mixed airflow on a physical model of fast water, the range of discharges is between 13 to 38 litres per second. The existing pumping system, the volume of water in the laboratory tanks, the geometric dimensions and the discharge-Ashle curve of the rapid water model were considered [18]. The design and construction of the physical model were based on the law of landing similarity. It is necessary to mention that in order to ensure that the “scale effect” is insignificant, an attempt was made to provide conditions in the model so that forces that are not large and effective in principle are not effective in the model as well [19]. These conditions were created by selecting the appropriate scale and geometry, controlling the minimum flow depth in the model, controlling the minimum surface roughness, and so on. In this research, the properties of the liquid used in the experiment are as follows:

The flow simulation was performed for the geometry similar to the device made in three dimensions so that the results could be compared well [20]. The effect of gravity was also considered, and it was also assumed that there is no mass transfer between the two phases and that both phases are incompressible; Because the length of the pipe is short and the pressure changes at this length are not large, according to the two-phase flow analysis in the transient state, a time step of 0.001 s was used, and the calculations continued until reaching the steady-state [21].

Apparent weather speeds are defined as follows:

VSL=QLAp (10)

VSG=QGAp (11)

The fluid which used in all cases enters the pipe vertically 998 kg/m3 and the water density was set at 225.1 kg/m3. Speed input limit condition was used for input, and the output current limit condition was used for output [22].

3  Results

3.1 Mesh Independence Analysis

The first step in examining the independence of the solution from the network is to generate a large network and simulate it. Under these conditions and after convergence solves, the desired result as shown in Fig. 2 for pressure in two-phase compensation. To do this, first shred the grid and increase the number of elements and cells to repeat the simulation using the new grid. Then, to check the independence of the solution from the network, the network is slightly smaller and the pressure drop is recorded. This process of refinement and simulation continues until the difference of the desired parameter in the two final simulations is very small and their difference is less than the desired error. In this case, the process of checking the independence of the solution from the network is stopped and the last one network that has a similar answer to the last network is considered as the optimal solution network.

Figure 2: Mesh independent analysis

3.2 Validation of the Results

There is currently a severe lack of information regarding appropriate standards and design methods for predicting the hydraulic characteristics of biphasic flows, especially in hydraulic structures due to lack of information and insufficient research, whether theoretical or laboratory. In such cases, it is usually based on previous experience or speculation [23]. Carefully in this form, the closeness of the information obtained in this study and the results of previous research can be seen. Therefore, at first glance, the flow pattern map obtained in this research can be a suitable tool for predicting two-phase flow patterns.

As shown in Fig. 3, the results of measurements related to the velocity and turbulence characteristics of two-phase airflow were performed completely for a number of dimensionless flow parameters [24]. The maximum water load on the inlet opening of the canal is between 25.89 to 95.93 cm and is equivalent to the flow rate of 0.0167 to 0.02 cubic meters per second per unit width. The landing number of the stream approaching the overflow is in the range of 0.663 to 1/007. The results of distribution and dispersion of variable values of climate flow velocity and parameter without turbulence intensity dimension indicate some of the characteristic features of this type of two-phase flow that the conditions of occurrence of two-phase flow patterns were obtained from the flow map according to the

Figure 3: Comparison of current test results with previous research

Reynolds number, Re=UDhv , where v is the kinetic viscosity, Dh is the hydraulic diameter, U is the mean single velocity of the single-phase liquid. Motion viscosity was calculated using the temperature of the fluid in the inlet manifold. The test conditions are given in Table 3 [26].

3.3 Numerical Results

Fig. 4 shows the two-phase flow pressure drop profiles. In general, channels are when the channels are hollow, which is intended to be compared to a porous media channel. The common feature of these three curves is the decrease in pressure drop, followed by a sudden increase. The reason is the transfer of the two-phase flow pattern from the stratified flow in one channel and the stagnant fluid in the other channel to the stratified currents in both channels.

Figure 4: Pressure drop profiles

Fig. 4 shows typical images of two-phase currents observed in a circular microchannel with different gas flow velocities and low liquid flow rates (Fig. 5a), medium liquid flow rates (Fig. 5b), and high liquid flow rates (Fig. 5c).

Figure 5: Typical images of two-phase currents

In a given flow condition, different flow patterns are observed which, depending on the surface configuration, can be classified into five distinct patterns: “liquid alone (or liquid snail) “,” gas core with the smooth thin liquid film “,” Gas core with the smooth liquid film “,” Gas core with the ring-shaped liquid film “and” Gas core with deformed interface “Liquid alone pattern “No bubbles or liquefied gas interface in the field Visibility did not show 1 mm and appeared under any circumstances. The smooth liquid films around the gas core can be classified as thin or thick, and the wavy liquid films can be classified as a ring or deformed [27,28].

When the biphasic current from the header is divided into parallel T channels, the surface tension effect becomes more apparent. In the two-phase air/water flow test, we find that all the liquids are divided into partially branched channels, but the rest of the branch channels are blocked by alternating air coils and liquid coils. Therefore, it is necessary to first determine the effect of surface tension and dredge all channels.

In this study, the flow behaviour of two phases of air and water in parallel channels of specified width with experimental porous inserts was investigated using a self-made and produced transparent assembly. In this work, we created a two-phase model for channel flow. Theoretical analysis was performed to calculate the liquid water saturation, the starting point of the two-phase flow and the concentration of species along the channel. Due to the expansion of industrial applications of two-phase downstream flow, its modelling and simulation are of great practical importance for predicting its various conditions in pipes. The dependence of these physical characteristics on cell performance of the main working parameters, such as air stoichiometry and relative humidity, is examined. The results show that liquid water accumulates rapidly in the inlet area, followed by a gentle increase downstream, at full moisture inlet conditions. The predicted liquid saturation rate reaches 30%. Subsequently, biphasic flows in heterogeneous channels were theoretically analyzed and numerically simulated. In this study, the flow pattern of downstream, bubble, Coal and foamy in the downstream flow of the pipe was clearly identified and photographed and using the flow map obtained from the laboratory study, the conditions for the occurrence of two-phase flow patterns were determined. The results of current research can be summarized as follows. First, in general, the experimental results agree well with the numerical predictions; however, some deviations are obvious when the fluid pressure drop is large. Second, porous channels consisting of spheres with smaller diameters have a higher heat transfer coefficient due to their larger heat transfer area.

4  Conclusion

The observed biphasic flow patterns were mainly alternating currents, But a closer look at the structure of the liquid film showed that the gas core flows with a smooth or ring-shaped film and a helical gas core flow surrounded by a deformed liquid. The probability of the emergence of different flow patterns indicates a gradual change in the dominant flow pattern with increasing gas and liquid flow velocities. For low liquid flow velocities, gas core currents with a thin, smooth liquid film often occur at low gas flow velocities, but gas kernel flows with a ring-shaped liquid film increase at high gas flow velocities. For the inlet flow, the gas and liquid distribution assume a transitional property. As the inlet gas surface velocity increases, more gas tends to the three middle channels, while more fluid is occupied in the first and end channels. Increasing the surface velocity of the inlet fluid does not change the inlet flow pattern in the header much, so the flow distribution changes slightly. When the input current changes to a ring, the inhomogeneous distribution of the two phases reaches its maximum. Liquid flows in large numbers to the first and last channels, while gas is mostly divided into three intermediate channels. However, in general, the relatively good agreement between the simulation results, especially the flow map and the laboratory flow map, shows the good ability of computational fluid dynamics in predicting low two-phase flow patterns.

Acknowledgement: The article was prepared as part of the state task “Research and development of complex energy-saving and thermoelectric regenerative systems” Application No. 2019–1497, Subject No. FZWG-2020-0034

Funding Statement: The authors received no specific funding for this study.

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study

This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.