Slug Flow Hydrodynamics Modeling for Gas–Liquid Two-Phase Flow in a Pipe

<jats:p>Gas–liquid flow in a pipeline is a very common. Slug two-phase flow is dominated in the case of slightly upward flow (+0.25°) and considered to be the comprehensive flow configuration, and can be in close contact with all the other flow patterns. The models of different flow patterns can be unified. Precise prediction of the slug flow is crucial for proper design and operation. In this paper, we develop hydrodynamics unified modeling for gas–liquid two-phase slug flow, and the bubble and droplet entrainment is optimized. For the important parameters (wall and interfacial friction factors, slug translational velocity and average slug length), the correlations of these parameters are optimized. Furthermore, the related parameters for liquid droplet and gas bubble entrainment are given. Accounting for the gas–liquid interface shape, hydrodynamics models, i.e., the flat interface model (FIM) and the double interface model (DIM), of liquid film in the slug body are applied and compared with the experimental data. The calculated results show that the predictions for the liquid holdup and pressure gradient of the DIM agree with experimental data better than those of the FIM. A comparison between the available experimental results and Zhang’s model calculations shows that the DIM model correctly describes the slug dynamics in gas–liquid pipe flow.</jats:p>

section_content: Gas-liquid flow in a pipeline is very common. Due to its remarkable economic benefit, it has been applied widely in geothermal fields, oil and gas industries, chemical reactors, nuclear reactors, aerosols, and various other heat and mass transfer processes during the past 40 years [1] [2] [3] [4]. Slug two-phase flow is dominated in the case of slightly upward flow (+0. 25 • ) and considered to be the comprehensive flow configuration [5, 6]. It can be in close contact with all other flow patterns; for example, the liquid film zone of the slug unit is similar to that of stratified flow, and the liquid slug body resembles that of bubbly or dispersed bubble flow [7, 8]. Thus, the study of the slug flow can also capture the behaviors of the other flow patterns [9, 10]. 

To optimize oil and gas production and transportation, flow behaviors such as flow pattern, liquid holdup, and pressure gradient need to be predicted accurately. Early models to predict gas-liquid hydrodynamic characteristics were based on empirical correlations [11] [12] [13]. The application of empirical correlations simplifies the complex physics and limits the opportunity to improve the models. The empirical correlations of liquid holdup and pressure drop scaled up to oilfield applications are inaccurate, because these correlations are obtained from small-scale laboratory data. 

Recently, mechanistic models were proposed to overcome this problem. Spedding et al. [14] developed a slug flow specification model to describe the slug initiation, growth, and decay, and tested 60 models on the basis of the 12 possible flow regimes. Bendiksen et al. [15] model for the simulation of two-phase oil and gas flow in pipelines, and predictions of liquid holdup and pressure gradient were found to agree with experimental data. De Henau and Raithby [16, 17] developed a one-dimensional transient two-fluid model with slug flow, which provides the necessary interphase constitutive relations for the drag coefficient and the parameters for the closure relations incorporated in the momentum conservation equations. Issa and Kempf [18] proposed a mechanistic one-dimensional model to predict hydrodynamic slug formation, growth, and decay, and the subsequent development of the slug shape into continuous slug flow, based on the numerical solution of the one-dimensional transient two-fluid model equations. The results of computations for the slug characteristics were considered to be in agreement with results of the experiments. Losi et al. [19] experimentally and theoretically investigated high viscosity oil-air slug flow, and built probability density functions for slug lengths based on superficial gas velocity. 

Arabi et al. [20] experimentally investigated the slug frequency of intermittent air-water flow in a horizontal pipe, investigated the influence of the sub-regime on slug frequency, and developed a new slug frequency correlation for intermittent flow. Gas-liquid two-phase slug flow has been investigated continuously in previous studies, partially because of its complexity. Every other flow pattern can be represented in slug flow, e. g., the film zone of slug flow resembles that of stratified or annular flow, and the slug body resembles that of bubbly or dispersed bubble flow [7]. Zhang et al. [7, 21] developed a unified model for gas-liquid pipe flow via slug dynamics. The relationship of each flow pattern can be obtained, through the transformation of the slug flow relationship, so as to achieve the unification of the calculation model for each flow pattern. However, although the model describes the phenomenon of gas-liquid entrainment, its processing method is too simple, and only considers the entrainment fraction of the two parameters of droplet and gas holdup. In fact, entrainment of droplets and bubbles occurs only under certain conditions. Furthermore, there is still a limited amount of information about the influence of the gas-liquid interface shape in slug flow hydrodynamics. A good description of interface shape is important when examining the liquid film zone of the slug unit in slug flow. Three typical interface configurations of stratified flow are presented in the literature: the flat interface model, the apparent rough surface model, and the concave down configuration. In the current study, a new mechanistic model was formulated and tested against experimental data, based on the comparison of these three interface configurations.
section_name: Introduction

section_content: The hydrodynamics model was derived for fully developed hydraulic gas-liquid slug flow, where there is no heat and mass transfer between gas and liquid phases. The continuity, momentum, and equations of slug flow were established by taking the liquid film region as the control volume. The momentum equation of slug flow can be transformed into the equation of stratified flow after the momentum exchange term is eliminated. The annular flow momentum equation can be obtained when the friction term of the gas wall is eliminated. The basic equations of each flow pattern are unified. 

Some assumptions are required because of the gas-liquid slug flow complexity: a fully developed slug flow, and negligible momentum exchange between the liquid slug and the liquid film.
section_name: Hydrodynamic Model

section_name: Slug Flow

section_content: The slug unit is defined as the control volume for the fully developed slug flow, as shown in Figure 1, including the entire liquid film zone and the gas pocket (liquid slug) zone. The liquid slug zone included in the liquid occupies the entire pipe cross-section containing the bubble. The liquid film including the liquid film and gas pocket (i. e., a Taylor bubble) is in close proximity to the liquid slug zone. U l, S l, and F l are lengths of the slug unit, liquid slug zone, and liquid film zone, respectively. The length of slug unit is defined as:

Considering the passage of a slug unit at an observation point, the following relationships hold for the liquid and gas phases, respectively:

Considering the entrainment of liquid droplets in the gas core, the gas continuity equation can be obtained by taking two cross-sections of the liquid film and liquid slug regions:

where

For a fully developed slug flow, the slug translational velocity refers to the velocity of the slug front (or the head of the long bubble). One of the key issues for hydrodynamic slug flow is the slug front velocity (or slug bubble velocity); this is one of the governing parameters for the modeling of hydrodynamic slug flow: The liquid slug zone included in the liquid occupies the entire pipe cross-section containing the bubble. The liquid film including the liquid film and gas pocket (i. e., a Taylor bubble) is in close proximity to the liquid slug zone. l U, l S, and l F are lengths of the slug unit, liquid slug zone, and liquid film zone, respectively. 

The length of slug unit is defined as:

Considering the passage of a slug unit at an observation point, the following relationships hold for the liquid and gas phases, respectively:

Considering the entrainment of liquid droplets in the gas core, the gas continuity equation can be obtained by taking two cross-sections of the liquid film and liquid slug regions:

where u T, u F, u C, and u S are velocities of the slug unit, liquid film, gas pocket, and gasliquid mixture in the liquid slug, respectively. R LC, R F and R LS are liquid holdups in the slug body, film, and liquid slug, respectively:

For a fully developed slug flow, the slug translational velocity refers to the velocity of the slug front (or the head of the long bubble). One of the key issues for hydrodynamic slug flow is the slug front velocity (or slug bubble velocity); this is one of the governing parameters for the modeling of hydrodynamic slug flow:

Energies 2022, 15, 533

where u d is the drift velocity. It is calculated according to:

where C is the empirical coefficient of speed distribution; for a fully developed pipeline cross-section speed distribution, C is approximately equal to the ratio of maximum speed to average speed. When the gas-liquid mixture in the liquid plug is sufficiently turbulent, C = 1. 2. When the gas-liquid mixture in the liquid plug is laminar flow, C = 2. 0. The liquid holdup in liquid slug, R LS, is given as:

Following the same derivation as the gas phase, the continuity equation of the liquid phase can be written as:

Adding Equation ( 4 ) to Equation ( 13 ):

where

where u SL and u SG are superficial velocities of liquid and gas:

The volume difference generated by the difference velocity between the gas pocket and liquid slug in liquid film zone is defined as

, and the volume difference generated by the difference velocity between liquid film and liquid slug is defined as ∆V 2 = R LF S P (u Fu S ). Volume differences ∆V 1 and ∆V 2 are equal, due to being inside the slug body, so:

According the Equations ( 2 ), ( 3 ), ( 4 ), (13), and ( 16 ), u S, R LF, u F, u C, and u T can be written as:
section_name: Continuity Equation

section_content: The liquid film and gas pocket in the film zone are separately used as the control volumes for the derivation of momentum equations. For fully developed slug flow, the momentum exchanges per unit time between the slug body and liquid film in the direction of the frictional and gravitational forces acting on the film are in balance. The liquid and gas phase momentum equations are written as:

Equating the pressure drop of the gas and liquid phase and subtracting Equation ( 23 ) from Equation (22) gives:

Eliminating τ i S i in both Equation (22) and Equation ( 23 ) gives:

The liquid slug momentum equation is written as:

i. e., -

where

Equation ( 27 ) can be given as:

The momentum transfer term between the liquid plug region and the liquid film region does not appear in the Equation (28). However, the momentum transfer between the two regions affects the slug flow characteristics.
section_name: Momentum Equation

section_content: The momentum equation for stratified flow is obtained by removing the momentum exchange terms from Equation ( 28 ), given as:

In annular flow, there is no contact between the gas core and the pipe wall. Equation ( 28 ) then becomes:

Energies 2022, 15, 533 6 of 16
section_name: Stratified and Annular Flows

section_content: Bubble flow includes dispersed bubble flow and bubble flow. The liquid holdup and pressure gradient of dispersed bubble flow are calculated when it is assumed that the gas-liquid mixture is uniform. For bubbly flows, the rate of bubble rise relative to the liquid must be considered; u o, is estimated by the Harmathy correlation:
section_name: Bubble Flow

section_content: The implementation of the model is dependent on both the basic equations and closure relationships. The wall friction factor, interfacial friction factor, slug length, and slug liquid holdup are the key parameters of the model.
section_name: Closure Relationships

section_content: The above continuity and momentum equations govern the characteristics of the liquid film. The components of shear stress in the combined momentum equations are calculated from:

The friction factors f L and f G of the liquid and gas phase, at the wall for the laminar flow, are calculated from:

The relationship between the above gas phase fanning friction factor, f G, and the Reynold number for the case of fully developed turbulent flow in a rough pipe, can be obtained from the Eck equation of Equation (22), as shown again here:

The liquid phase fanning friction factor, f L, can be obtained from modified Grolman model [22] :

Re sl -1 Re G 0. 25 for 2100 < Re L (37)
section_name: Wall Friction Factor

section_content: The wetted wall fraction is the fraction of the pipe perimeter wetted by the liquid phase. For a given pipe, if either the wetted wall fraction or the wetted wall perimeter are known, the other can be easily calculated. The behavior of the liquid film is similar to stratified flow. A good description of the interface shape is important when performing calculations of stratified gas-liquid two-phase flow. Figure 2 gives the schematic diagrams of three typical interface configurations of stratified flow in the literature. The flat interface model, which assumes that the interface configuration is flat and smooth, is widely used to predict gas-liquid stratified flow in a pipe, as shown in Figure 2a. The apparent rough surface model assumes a constant thickness liquid film around the pipe wall, i. e., the interface has a constant curvature circumferentially and the interface is a portion of an arc, as seen in Figure 2b. This model is not realistic, because liquid will accumulate more at the bottom of the pipe than away from the bottom due to the gravitational force. The gas-liquid interface exhibited a concave down configuration in experiments. In order to better simulate a concave down gas-liquid interface, Chen et al. [23] proposed a double-circle model, as shown in Figure 2c. 

calculations of stratified gas-liquid two-phase flow. Figure 2 gives the schematic diagrams of three typical interface configurations of stratified flow in the literature. The flat interface model, which assumes that the interface configuration is flat and smooth, is widely used to predict gas-liquid stratified flow in a pipe, as shown in Figure 2a. The apparent rough surface model assumes a constant thickness liquid film around the pipe wall, i. e., the interface has a constant curvature circumferentially and the interface is a portion of an arc, as seen in Figure 2b. This model is not realistic, because liquid will accumulate more at the bottom of the pipe than away from the bottom due to the gravitational force. The gasliquid interface exhibited a concave down configuration in experiments. In order to better simulate a concave down gas-liquid interface, Chen et al. [23] proposed a double-circle model, as shown in Figure 2c.
section_name: Wetted Wall Fraction

section_content: In Figure 2a, the dashed line is the gas-liquid interface. When the gas-liquid interface is flat, the relationship between the central angle, 0 θ, and the liquid holdup, L H, can be written as:
section_name: Flat Interface Model

section_content: In Figure 2a, the dashed line is the gas-liquid interface. When the gas-liquid interface is flat, the relationship between the central angle, θ 0, and the liquid holdup, H L, can be written as:
section_name: Flat Interface Model

section_content: In Figure 2c, the liquid layer was bounded between the eccentric circle (the center O 2 ) and the pipe wall circle (the center O 1 ). The radiuses of the two circles are R 2 and R 1, respectively. The point O 2 is at a distance above the point O 1. There are two extreme cases of the double circle model. When the radius of the eccentric circle R 2 is very large, the interface approaches a flat surface, which is similar to the flat interface model. When the radius of the eccentric circle R 2 is less than the pipe radius R 1, it is a constant thickness liquid film. 

As mentioned above, the gas-liquid interface configuration (arc-shaped ADB) is a portion of circle O 2, as shown in Figure 3. In order to describe quantitatively the arc-shaped ADB, the location and size of circle O 2 must be determined, i. e., the distance, ∆, above the center of the pipe wall and the radius, R 2, of the circle O 2, or the central angle, θ 2 and the radius, R 2. In the double circle model, the relationship between the liquid holdup, H L, and the central angles, θ 1, and θ 2 is given as:

The relationship between θ 1, θ 2, R 2, and R 1 is given by:

According to the definition of the wetted wall fraction, Θ = θ 1 /π, for a given pipe, if the wetted wall fraction, Θ, is known, the other can be easily calculated. The closure equation is needed for the wetted wall fraction. 

The relationship between 1 θ, 2 θ, R2, and R1 is given by: ( )

According to the definition of the wetted wall fraction,

, for a given pipe, if the wetted wall fraction, Θ, is known, the other can be easily calculated. The closure equation is needed for the wetted wall fraction.
section_name: .2. Double Circle Model

section_content: θ π

The physical mechanism for a liquid film climbing up the pipe wall during stratified flow is still not well understood. This phenomenon is affected by many factors, such as The physical mechanism for a liquid film climbing up the pipe wall during stratified flow is still not well understood. This phenomenon is affected by many factors, such as the superficial gas velocity, gas density, and liquid phase surface tension. Under this circumstance, the wetted wall fraction can be calculated by the Fan correlation: where w SG, w SG,C = 5(1. 24/ρ G ) 0. 5 are the superficial gas velocity and the critical superficial gas velocity, m/s, respectively. σ w, σ L are the water and liquid surface tension, N/m, respectively. σ L is calculated by the Viles correlation [24]. Fr L is the liquid phase Froude number, defined as:

where u L = u SL /H L is the average velocity of the liquid phase.
section_name: -

section_content: The gas-liquid interfacial friction factor is another indispensable closure relationship which is necessary for solving the momentum equations. It is known that the interfacial shear stress greatly affects momentum transfer processes in stratified flow.
section_name: Interfacial Friction Factor

section_content: In the geometric dimensioning of oil and gas in a circular pipe, the gas-liquid interface shape is assumed to be flat. For laminar flow, the friction factor f G can be given by the Hagen-Poiseille correlation [25]. For turbulent flow, f G can be obtained by the Hall correlation [25] and the Chen explicit equation [23]. The friction factor f L can be calculated in the same way as for f G, or by the Ouyang-Aziz correlation [26] and Spedding-Hand correlation [27]. 

The Andritsos-Hanratty correlation [28] is used to obtain the interfacial friction factor f i :

where u SG,t = 5 ρ 0 G /ρ G 1/2.
section_name: Flat Interface Model

section_content: For the double circle model, the friction factor f i can be obtained by the Ottens correlation [29]. For stratified smooth flow, i. e., u SG ≤ u SG,t, f i = f G, and for stratified wavy flow, i. e., u SG > u SG,t :

where:

The liquid film velocity at the liquid-gas interface, u i, is given as:

where the average height of liquid film is given as:

(49)
section_name: Double Circle Model

section_content: The liquid holdup of the slug, H ls, can be calculated by the unified mechanistic model [7, 21] for any inclination angle. It is derived from the balance between the surface free energy of dispersed bubbles in a slug body and the liquid phase:

where:

The slug liquid holdup should be estimated before solving the combined momentum equitation. It can be given by the Gregory et al. correlation [30] : The slug length can be estimated as:

l s = 32. 0 cos 2 α + 16. 0 sin 2 α d (53)
section_name: Slug Liquid Holdup

section_content: To evaluate the hydrodynamic flat interface model (FIM) and the double interface model (DIM), they were applied to the calculation of liquid holdup and pressure gradient, and compared with the experimental data [21]. Then, the hydrodynamic model was verified with experimental measurements, including of the pressure gradient and liquid holdup. For further verification, the model was compared with that developed by Zhang et al [7]. 

Table 1 summarizes the experimental conditions. Air, kerosene, and white oil were used in the experiment for the two-phase fluids. The pipes' inside diameters were 51, 76. 3, and 40 mm. The pipe inclination angles varied from 0 to 90 deg. The flow pattern of the experiment was always that of slug flow.
section_name: Results and Discussion

section_content: Figure 4 shows the comparison for the liquid holdup (an average value of the slug unit), between the two model (i. e., FIM and DIM) predictions and the measured data of the No. 1 experiment. The liquid holdup of slug flow was calculated using the liquid film on the lengths of the slug body and the liquid holdup of the slug body. The predictions of both models and the experimental measurements show good agreement, when the liquid holdup is greater than 0. 5, as seen in Figure 4. However, the FIM resulted in greater deviation compared with the experimental data when the liquid holdup decreases to 0. 5. In addition, the deviation appears to increase with a decrease in the liquid holdup. The calculations of the DIM are consistently in good accordance with the experiment data. and compared with the experimental data [21]. Then, the hydrodynamic model was verified with experimental measurements, including of the pressure gradient and liquid holdup. For further verification, the model was compared with that developed by Zhang et al [7]. Table 1 summarizes the experimental conditions. Air, kerosene, and white oil were used in the experiment for the two-phase fluids. The pipes' inside diameters were 51, 76. 3, and 40 mm. The pipe inclination angles varied from 0 to 90 deg. The flow pattern of the experiment was always that of slug flow.
section_name: Comparison of FIM and DIM

section_content: Figure 4 shows the comparison for the liquid holdup (an average value of the slug unit), between the two model (i. e., FIM and DIM) predictions and the measured data of the No. 1 experiment. The liquid holdup of slug flow was calculated using the liquid film on the lengths of the slug body and the liquid holdup of the slug body. The predictions of both models and the experimental measurements show good agreement, when the liquid holdup is greater than 0. 5, as seen in Figure 4. However, the FIM resulted in greater deviation compared with the experimental data when the liquid holdup decreases to 0. 5. In addition, the deviation appears to increase with a decrease in the liquid holdup. The calculations of the DIM are consistently in good accordance with the experiment data. In Figure 5, the predictions of pressure gradients for slug flow are compared against the measured data. The pressure gradients calculated in the two models include the shear stress at the pipe wall and the momentum exchange between the liquid film and the slug. In Figure 5, the predictions of pressure gradients for slug flow are compared against the measured data. The pressure gradients calculated in the two models include the shear stress at the pipe wall and the momentum exchange between the liquid film and the slug. As can be noticed, the predictions of the DIM agree with the experimental data better than those of the FIM, and a large number of the predictions of the DIM fall into the ±25% range. The liquid holdup prediction of the FIM is greater than that of measurement data. The liquid phase velocity decreases due to the dragging effect, so the liquid holdup must increase. The flat interface of the FIM generates a smaller dragging effect on the liquid phase compared with the curved interface of the DIM. 

The pressure gradient of the FIM is underestimated, probably because of the underprediction of the wetted wall fraction, due to the creep of the liquid film up the pipe wall, which tends to drag more liquid into the slug body. The level of under-prediction is more intensified when the gas superficial velocity is higher, with a constant liquid flow rate. Due to the flat interface along with lower flow resistance, in order to maintain the same flow rates, a smaller pressure gradient is necessary. With decreasing the gas flow rate, the flow resistance decreases, so a smaller pressure gradient is produced. As a result, the differences in the calculated prediction pressure gradient between the two interface assumptions are more uniform.
section_name: Comparison of FIM and DIM

section_content: The proposed DIM was validated with the measurements of the No. 2 experiment. An experimental study of gas-liquid slug flows with pipe inclinations varying from horizontal to upward vertical was carried out. Figure 6 shows the comparison between predictions of the DIM and the experimental results for liquid holdup. The predictions of liquid holdup are in good agreement with the experimental data, and all the predictions fall into the ±15% range. The calculation gives under-predictions of liquid holdup and the maximum absolute error is 14. 8%, which does not exceed 15%. The liquid holdup prediction of the FIM is greater than that of measurement data. The liquid phase velocity decreases due to the dragging effect, so the liquid holdup must increase. The flat interface of the FIM generates a smaller dragging effect on the liquid phase compared with the curved interface of the DIM. 

The pressure gradient of the FIM is underestimated, probably because of the underprediction of the wetted wall fraction, due to the creep of the liquid film up the pipe wall, which tends to drag more liquid into the slug body. The level of under-prediction is more intensified when the gas superficial velocity is higher, with a constant liquid flow rate. Due to the flat interface along with lower flow resistance, in order to maintain the same flow rates, a smaller pressure gradient is necessary. With decreasing the gas flow rate, the flow resistance decreases, so a smaller pressure gradient is produced. As a result, the differences in the calculated prediction pressure gradient between the two interface assumptions are more uniform.
section_name: Compared with Experiment

section_content: The proposed DIM was validated with the measurements of the No. 2 experiment. An experimental study of gas-liquid slug flows with pipe inclinations varying from horizontal to upward vertical was carried out. Figure 6 shows the comparison between predictions of the DIM and the experimental results for liquid holdup. The predictions of liquid holdup are in good agreement with the experimental data, and all the predictions fall into the ±15% range. The calculation gives under-predictions of liquid holdup and the maximum absolute error is 14. 8%, which does not exceed 15%. Figure 7 gives the DIM predictions of the pressure gradient compared to the experimental data. The predicted results of the pressure gradient assuming a gas-liquid concave interface exceed the measured data, and fall into the 0-25% range. Figure 7 gives the DIM predictions of the pressure gradient compared to the experimental data. The predicted results of the pressure gradient assuming a gas-liquid concave interface exceed the measured data, and fall into the 0-25% range.
section_name: Compared with Experiment

section_content: The accuracy of the model was verified by comparing the result of experiment No. 3 [31] with Zhang's model [7] and the model in this study. The comparison results are shown in Figures 8 and 9, respectively. Table 2 shows the error statistics of the calculation results of the two models. Two kinds of liquid water and white oil were used in the experi-ment, and air was used as the gas phase medium. The gas apparent velocity is 0. 06-20 m/s and the liquid apparent velocity is 0. 08-2 m/s.
section_name: Comparison with Zhang's Model

section_content: The accuracy of the model was verified by comparing the result of experiment No. 3 [31] with Zhang's model [7] and the model in this study. The comparison results are shown in Figures 8 and 9, respectively. Table 2 shows the statistics of the calculation results of the two models. Two kinds of liquid water and white oil were used in the experiment, and air was used as the gas phase medium. The gas apparent velocity is 0. 06-20 m/s and the liquid apparent velocity is 0. 08-2m/s.
section_name: Comparison with Zhang's Model

section_content: The accuracy of the model was verified by comparing the result of experiment No. 3 [31] with Zhang's model [7] and the model in this study. The comparison results are shown in Figures 8 and 9, respectively. Table 2 shows the error statistics of the calculation results of the two models. Two kinds of liquid water and white oil were used in the experiment, and air was used as the gas phase medium. The gas apparent velocity is 0. 06-20 m/s and the liquid apparent velocity is 0. 08-2m/s. The errors in the calculated result using Zhang's model are relatively large, with an average error of about 15%. In addition, the proportions of errors less than ±15% of liquid holdup and the pressure gradient calculated by Zhang's model are fewer compared to the developed model. The errors in the calculated result of the developed model are evenly distributed, and most of the errors are within the range of 15%. The proportions of errors less than ±15% of liquid holdup and the pressure gradient are 72. 5% and 82. 4% respectively. The average errors of liquid holdup and the pressure gradient are 8. 58% and 9. 95% respectively. Therefore, the calculation accuracy of the model is high, due to considering the gas-liquid entrainment condition.
section_name: Comparison with Zhang's Model

section_content: A hydrodynamic model was developed for gas-liquid two-phase slug flow. The model used the entire liquid film region as the control body, and momentum exchange between the liquid film region and the slug region was introduced into the model. The governing equation can be applied to other flow patterns after deformation, by realizing the unification of various flow patterns. 

The gas-liquid entrainment is described in more detail in the model, and entrainment conditions were added into the model. For the important parameters (wall and interfacial friction factors, slug translational velocity, and average slug length), the correlations of these parameters were optimized. Furthermore, the related parameters for the entrainment of liquid droplets and gas bubbles are given in this paper. 

Accounting for the gas-liquid interface configuration of the liquid film in the slug body is important; thus, two models (DIM assumed a concave interface, FIM assumed a flat interface) were used for the hydrodynamic predictions. The calculations of pressure gradient and liquid holdup for the DIM were in good accordance with experiment data. For further verification, the calculations of the model were compared with the experimental data and those predicted by Zhang's model. The errors in the calculated results using the developed model were evenly distributed, and most of the errors were within the range of 15%.
section_name: Conclusions

section_content: Funding: This research was funded by the Natural Science Foundation of China ( 51604282 ), Natural Science Foundation of Chongqing ( cstc2019jcyj-msxmX0268, cstc2019jcyj-msxmX0286 and cstc2019jcyj-msxmX0628 ). 

Informed Consent Statement: Not applicable.