# Mathematical Model of Air Pocket Evolution during Water Filling a Long-Slope Pipeline and Its Application to Air Removal Prediction

Authors: Yuanyuan Chen; Yanyu Cui; Wenhui Chen; Tao Deng; Jing Gong
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

During water filling a long-slope pipeline, air pocket is very likely to be entrapped at peak point. In order to track air movement and predict air removal conditions, a mathematical model of air pocket evolution, including its formation, compression, and entrainment, is proposed in this paper. The simulation results were compared with the engineering field data and the two are basically consistent. Furthermore, the two most important factors which play a great role in the removal of air pocket, i.e., the terrain category and the inlet flow rate, are analyzed in detail. It is concluded that the removal conditions reach three outcomes: air pocket compressed and partly removed, compressed and completely removed, and compressed without any removal. In this paper, the terrain which leads to the last outcome is called the “Dangerous Terrain.” And for the “Dangerous Terrain,” it is of great importance that the inlet flow rate should be strictly confined within a certain level. While for the other two categories of terrains, an increased flow rate is in any respect beneficial to the removal of air pocket.


---

## Body

## 1. Introduction

Air pocket is very likely to be entrapped at high points during water filling a long-slope pipe. The entrapped air, if cannot be removed properly and timely, may cause a series of probable damages, such as water column separation, pipeline transport efficiency reduction, power consumption increase, etc. [1–3]. In some severe cases, it may even cause a pipeline burst accident due to the insufficient awareness of the air–water interaction flow characteristics in the pipeline.Pipeline systems with entrapped air pocket have been widely discussed both experimentally and theoretically. Leon et al. [4] used the lumped approach to calculate the pressure of air pocket entrapped between two pressurized columns. Pothof and Clemens [5] established the mathematical models for the calculation of air pocket motion which neglected the effect of surface tension and viscosity. They pointed out that air pockets can be avoided provided that pipe diameter exceeds 191 mm. Deng et al. [6, 7] conducted a series of experiments both in laboratory and field trials to simulate the water pressure testing for long pipeline to ensure the safety and efficient construction thereof. They found that, during the last stage of water filling, an instantaneous pressure impulse will be induced due to the existence of air pocket in the sloping reach of pipeline.Different from the traditional two-phase flow, air entrapment flow is characterized by the inlet boundary of only liquid velocity rather than gas–liquid two phases. The air entrapment herein refers to a “gaseous section” where the gas is entrapped by the liquid slug with a constant mass at a stagnant position. Under certain circumstances, the liquid slug front downstream is likely to entrain the pocket tail upstream, resulting in the entrainment of air pocket as well as aeration of liquid. This phenomenon is often referred to as “air entrainment” in literatures. It is suggested by Fan et al. [8] and Bonetto and Lahey [9] that the air entrainment process is similar to a plunging of vertical jet into a static pool. Jepson [10] observed from the experiments that the “aeration” of liquid caused by air entrainment is quite similar to the process where air is entrained by water hydraulic jump [11–14]. By using clean water, surfactant-laden water, and wastewater as test media, respectively, Pothof et al. [15] conducted a systematic study on the influence of surface tension on air entrainment. It was found in the study that the smaller the surface tension, the easier the bubble to be entrained in respect of declined pipes.However, it should be noted that most of the current research is still based on horizontal or near-horizontal pipes or individual downward-sloping pipes, rather than pipelines laid along mountainous areas with continuous long-slope and typical V-shaped pipes. For the latter, air-water flow will be much more complex. Not only because it is highly dependent on the terrain category, but also because it is related to the complicated evolution of the air pocket, i.e., formation, compression, and entrainment thereof. Although many considerable works on air entrapment have been carried out to date, due to the weak understanding of the air–water interaction mechanism in the long-slope pipeline, the flow characteristics associated with air pocket formation, compression, and entrainment have been rarely studied.Therefore, based on a comprehensive consideration of formation, compression, and entrainment of air pocket, a mathematical model of air pocket evolution is proposed in this paper to track air movement and predict air removal conditions in long-slope pipelines. With the model proposed, a computational support can be provided for the determination of pressure and location and length of air pocket, and the remaining gas within the pipeline can also be estimated and evaluated beforehand without having to stop the transmission. Based on all the above, the model is considered to have a significant meaning in laying the theoretical and scientific basis on the optimal design and safe operation for water filling pipelines with long slopes.


## 2. Mathematical Model

Before establishing the model, the following is assumed:(1)
The fluid system is thermostatic(2)
The pipeline outlet keeps open to the atmosphere, i.e., the outlet pressure remains constant as atmospheric(3)
The pipeline inlet flow rate remains constant throughout the process(4)
The liquid is incompressible, and the gas is ideal(5)
The flow is unidimensional, and the uneven distribution of pressure, velocity, and liquid holdup within the pipeline cross section is neglected


### 2.1. Air Pocket Formation

Before water filling, the pipeline is emptied with static atmospheric air. After the water is pumped and reaches the peak point, it will flow towards the valley point in a form of liquid film similar to the open-channel flow (Figure1).Figure 1
Schematic diagram of open channel flow.The average liquid velocity is given by the open-channel flow model as follows:(1)vlf=1.49nAlSl2/3sinθ11/2.The liquid cross-sectional area and wetted perimeter can be expressed as follows:(2)Al=D24δ−12sin2δ,(3)Sl=Dδ.In addition, it is important to note that(4)Al=Qlvlf.When the water reaches the valley point, it will start accumulating there and gradually grow into a liquid slug which blocks and entraps the air upstream. In this way, the air pocket begins to take shape along the sloping reach. Based on (1)–(4), the initial liquid holdup by the time the air pocket is formed can be calculated.


### 2.2. Air Pocket Compression

#### 2.2.1. Governing Equations

After the air pocket is formed, the liquid slug at the valley point will start growing towards both upstream and downstream sloping reaches, resulting in the compression of air pocket. In this stage, the entrainment of air pocket is not considered. By specifying the V-shaped segment of gas–liquid two-phase flow as the object of study, the mathematical model for liquid slug growth and air pocket compression is established as shown in Figure2.Figure 2
Schematic diagram of liquid slug growth and air pocket compression.(1) Air Continuity Equation. Since throughout the water filling process, only liquid phase is pumped in from the inlet without any addition of air to the flow, thus, there is no mass exchange between the air pocket and the external environment. By specifying the air pocket segment as the controlled volume, the air continuity equation can be expressed as follows:(5)dmgdt=0.According to the assumptions, the air pocket meets the ideal gas equation. Thus, the mass of air pocket in the controlled volume can be expressed as(6)mg=MgARTPgLgHg.By substituting the above into (5), it can be deduced that(7)dPgLgHgdt=0.(2) Liquid Continuity Equation. By specifying the liquid in the V-shaped segment (the blue section in Figure 2) as the controlled volume, the liquid continuity equation is then given by(8)dmldt=ρlQl.The mass of liquid within the controlled volume is given by (9) based on the assumptions of D/h1≪1 and D/h2≪1:(9)ml=ρlAL1−LgHg+h2sinθ2.By substituting (9) into (8), it can be deduced that(10)1sinθ2dh2dt−dLgHgdt=vsl.(3) Momentum Equation for the Stratified Flow. The most probable flow pattern that could occur along the downward slope is the stratified flow. The quasisteady state force balance equation is used as the momentum equation:(11)ρgAggsinθ1−τgSg+τiSi−AgdPgdx=0,(12)ρlAlgsinθ1−τlSl−τiSi−AldPgdx=0.Combining the two equations above by eliminating the pressure gradient term, the combined momentum equation is then deduced as(13)ρl−ρggsinθ1−τlSlA1−Hg+τgSgAHg−τiSi1A1−Hg+1AHg=0.Among which,(14)τg=fgρgvg22,τl=flρlvlf22,τi=fiρgvg−vlf22.During the filling, the air pocket is compressed under the effect of the incoming flow upstream and liquid slug downstream and after all entrapped at a stagnant position, i.e.,vg=0. In this way, (13) and (14) can be further written as(15)ρl−ρggsinθ1−τlSlA1−Hg−τiSi1A1−Hg+1AHg=0,(16)τg=0,τl=flρlvlf22,τi=fiρgvlf22.(4) Momentum Equation for the Liquid Slug. Since the compression process of air pocket is directly related to the motion of liquid slug, the slug is hereby specified as the controlled volume for which a momentum equation is established. Based on the fact that the resultant force that acts on the controlled volume is equal to the difference of liquid momentum per unit time by subtracting those flowing in the controlled volume from those flowing out, the following equation can be established:(17)Pg−PaA−ΔPfh1sinθ1+h2sinθ2−ρlAh2−h1g=ddtρlAh1sinθ1+h2sinθ2u−ρlQlvlf,where ΔPf is the friction force per unit length and the expression is specified as(18)ΔPf=18πDρλu2,And u is the speed at which the liquid slug advances upward in the upward slope, and it can be calculated using(19)u=1sinθ2dh2dt.During the growth of liquid slug toward the upstream in the downward sloping reach, the liquid rises at a rate equal to the shortening of air pocket:(20)1sinθ1dh1dt+dLgdt=0.


#### 2.2.2. Closure Models

(1) Geometric Relationship. Refer to Figure 1; for the air pocket, the following hold: (21)Ag=D24π−δ+12sin2δ,(22)Sg=Dπ−δ,(23)Si=Dsinδ.Based on the above, the void fraction can be calculated as follows:(24)Hg=1−1πδ−12sin2δ.(2) Frictional Coefficient. Since the air pocket is nearly stagnant throughout the filling process, the air frictional coefficient fg is therefore neglected.Spedding and Hand [16] used (25) to calculate the the liquid frictional coefficient between the liquid and pipe wall:(25)fl=24Relf−1,Relf<2100,0.0262HlResl−0.139,Relf>2100.Among which,(26)Relf=Dvlfρlμl,Resl=Dvslρlμl,Hl=1−Hg.Xiao et al. [17] commented on 16 computational methods for interphase frictional coefficient fi calculation after testing them with the multiphase flow database of Tulsa University. The Andritsos–Hanratty correlation is after all recommended when D≤0.127m. When D>0.127m, the corrected Duns and Ros correlation by Baker is preferred. In the situation where the flow regime is to be determined, the recommended interphase frictional coefficient by Shoham and Taitel [18] is adopted, i.e., fi=0.0142, which is the preferable average for interphase wave of small amplitude and is applicable to a diversity of pipe diameters. According to the field pressure drop data of AGA large-diameter pipeline, the computational method for fi recommended by Xiao is well supported.


### 2.3. Air Pocket Entrainment

As the liquid slug grows, the liquid level heights both in the downward and in the upward sloping reaches continue to increase accordingly. This paper assumes that the highest liquid level of upward sloping reach, i.e., the waterfront arrives at the upward summit, marks the end of the slug growth. At this moment, if the liquid level height of the downward sloping reach is lower than that of the upward one, the difference of water levels in upstream and downstream reaches induces backpressure. Furthermore, assuming that only the excess of the backpressure above the air pocket pressure is responsible for the air pocket entrainment, the criterion for the onset of entrainment is thus given by(27)Pb=ρlgL2sinθ2−h1′>Pg.Once entrainment occurs, due to the intaking of large numbers of entrained small bubbles from the pocket tail, the fluid between pocket tail and downstream liquid slug front starts to occur in a form of aerated liquid slug as shown in Figure3.Figure 3
Schematic diagram of air entrainment at the pocket tail.In the figure, “T” represents the tail of air pocket. Bubbles are entrained from the pocket tail at the rate of Φen. Meanwhile, part of the entrained bubbles will flow back and recoalesce with the pocket tail. However, considering the fact that the elongated air pockets in long-distance pipeline are always with great length (normally several kilometers), the amount of entrained small bubbles is usually much more than that of the flow back. Therefore, the effect of bubbles flow back on the entrainment is not considered in this study. All of these bubbles entrained at the rate of Φen will enter the liquid slug downstream to form a bubbly flow.Brauner and Ullmann [19] assume that the entrainment rate Φen is dependent on the balance between the rate of turbulent kinetic energy production and the rate of bubble surface energy production and subsequently deduced that(28)Φenvsl=1400CJd˜critWe−Wec,where(29)We=ρlDvlf−vsl2σ,(30)Wec=100C′d˜crit.In fact,CJ and C′ are tunable constants, where CJ can be tuned based on turbulence intensity and critical bubble size. It is suggested by Brauner and Ullmann that CJ=1,C′=2/3. And this value is found to be in agreement with the experimental data of Nydal and Andreussi [20] in a nearly horizontal 5 cm pipe.The critical size of deformable bubbles is suggested by Brodkey [21]:(31)d˜crit=0.4σρl−ρggcosθ1′D21/2=0.224cosθ1′EoD,EoD≥0.2,0.25,EoD<0.2,Among which,(32)EoD=ΔρgD28σ,θ1′=θ1,θ1≤45∘,90−θ1,θ1>45∘.Regardless of the possible coalescence and buildup of entrained small bubbles in the slug under the effect of buoyancy, it is considered that the entrained small bubbles will simply mix with the liquid downstream to form a bubbly flow. The bubbly flow in the slug region can be described by the drift-flux model. The bubble velocity in liquid slug at pocket tailvblsT is given by(33)vblsT=C0Tvm−vdT,where C0T is the velocity distribution parameter in the near bubble tail region. The rigorous computation of C0T requires the bubble velocity distributions in these regions, which are quite difficult to obtain. However, wall peaking in void fraction distributions generally results in C0T<1, whereas C0T>1 characterizes center peaking. In the case of horizontal or inclined pipeline especially when the liquid velocity is relatively low, C0T<1 is generally considered. This is because the liquid carrying capacity is relatively weak in this situation and the entrained bubbles will normally gather around the upper wall under the buoyancy effect. According to the experimental results by Hernandez and Fabre [22], the recommended value of C0T as 0.5 is adopted in this paper.vm is the velocity of mixture flux which is calculated as follows:(34)vm=vsl+Φen.vdT is the drift velocity. By using the Harmathy model corrected by Wallis [23] for the effect of pipe inclination, the drift velocity is given as follows:(35)vdT=1.53σgρl−ρgρl20.25sinθ11−HglsT1.5.According to the mass conservation, the gas void fraction in liquid slug at the pocket tailHglsT is(36)HglsT=ΦenvblsT.Due to the air entrainment, the pressure, length, and mass of the air pocket decreases constantly. The pressure during the entrainment is given by(37)Pgk+1=Pgk−ΔPg,where k is the current calculation time step and k+1 is the next one. ΔPg is obtained by integrating (11).The change rate of mass of entrained air can be calculated by(38)dmgendt=−ρgTAHglsTvblsT,where ρgT can be deduced from the pressure Pgk+1. And then the remaining air mass and air length after the air pocket is entrained can thus be expressed as(39)mgk+1=mgk−Δmgen,(40)Lgk+1=1−ΔmgenmgkLgk.


### 2.4. Computational Procedure

By specifying the time when the liquid phase is pumped into the pipeline from the inlet as the simulation start time, the initial conditions at the moment of air pocket taking shape are given by(41)Hg0=1−1πδ0−12sin2δ0,(42)Pg0=Pa,(43)Lg0=L1,(44)ρg0=ρa,(45)mg0=ρg0AHg0Lg0,(46)h10=h20=D,(47)u0=Ql2A,where δ0 is iteratively solved by the open-channel flow model in Section 2.1. Then, the air pocket compression model in Section 2.2 is initiated to conduct the calculation. The five ordinary differential equations (ODE set) of (7), (10), (17), (18), and (19) are solved using the forward difference method to obtain Pg,Lg,h1,h2,u and thenHg by (15). The waterfront arriving at the summit of the upward sloping reach marks the end of the slug growth. At this moment, whether the criterion for the onset of entrainment, i.e., (27), is satisfied is judged. If it does not meet the criterion, the simulation stops. Otherwise, the entrainment model in Section 2.3 is activated to continue the calculation as follows. At each calculation time step, (28)–(32) are firstly replaced by the previously solved variables such as Pg,Lg,Hg and other required information as input to obtain Φen and then vm by (34). Secondly, based on the calculated Φen and vm, the remaining three implicit algebraic equations, i.e., (33), (35), and (36), are iteratively solved to obtain HglsT,vblsT,vdT. Finally, the pressure during the entrainment Pgk+1, the remaining air mass mgk+1, and the remaining air length Lgk+1 are deduced by (37), (39), and (40), respectively. After the current time step is finished, the above entrainment process is repeated again to calculate the next time step until the end of the simulation. The flowchart for the entire computational procedure is shown in Figure 4. The time steps were taken as 10 s and 300 s for compression process and entrainment process, respectively.Figure 4
Flow chart of the entire computational procedure.


## 2.1. Air Pocket Formation

Before water filling, the pipeline is emptied with static atmospheric air. After the water is pumped and reaches the peak point, it will flow towards the valley point in a form of liquid film similar to the open-channel flow (Figure1).Figure 1
Schematic diagram of open channel flow.The average liquid velocity is given by the open-channel flow model as follows:(1)vlf=1.49nAlSl2/3sinθ11/2.The liquid cross-sectional area and wetted perimeter can be expressed as follows:(2)Al=D24δ−12sin2δ,(3)Sl=Dδ.In addition, it is important to note that(4)Al=Qlvlf.When the water reaches the valley point, it will start accumulating there and gradually grow into a liquid slug which blocks and entraps the air upstream. In this way, the air pocket begins to take shape along the sloping reach. Based on (1)–(4), the initial liquid holdup by the time the air pocket is formed can be calculated.


## 2.2. Air Pocket Compression

### 2.2.1. Governing Equations

After the air pocket is formed, the liquid slug at the valley point will start growing towards both upstream and downstream sloping reaches, resulting in the compression of air pocket. In this stage, the entrainment of air pocket is not considered. By specifying the V-shaped segment of gas–liquid two-phase flow as the object of study, the mathematical model for liquid slug growth and air pocket compression is established as shown in Figure2.Figure 2
Schematic diagram of liquid slug growth and air pocket compression.(1) Air Continuity Equation. Since throughout the water filling process, only liquid phase is pumped in from the inlet without any addition of air to the flow, thus, there is no mass exchange between the air pocket and the external environment. By specifying the air pocket segment as the controlled volume, the air continuity equation can be expressed as follows:(5)dmgdt=0.According to the assumptions, the air pocket meets the ideal gas equation. Thus, the mass of air pocket in the controlled volume can be expressed as(6)mg=MgARTPgLgHg.By substituting the above into (5), it can be deduced that(7)dPgLgHgdt=0.(2) Liquid Continuity Equation. By specifying the liquid in the V-shaped segment (the blue section in Figure 2) as the controlled volume, the liquid continuity equation is then given by(8)dmldt=ρlQl.The mass of liquid within the controlled volume is given by (9) based on the assumptions of D/h1≪1 and D/h2≪1:(9)ml=ρlAL1−LgHg+h2sinθ2.By substituting (9) into (8), it can be deduced that(10)1sinθ2dh2dt−dLgHgdt=vsl.(3) Momentum Equation for the Stratified Flow. The most probable flow pattern that could occur along the downward slope is the stratified flow. The quasisteady state force balance equation is used as the momentum equation:(11)ρgAggsinθ1−τgSg+τiSi−AgdPgdx=0,(12)ρlAlgsinθ1−τlSl−τiSi−AldPgdx=0.Combining the two equations above by eliminating the pressure gradient term, the combined momentum equation is then deduced as(13)ρl−ρggsinθ1−τlSlA1−Hg+τgSgAHg−τiSi1A1−Hg+1AHg=0.Among which,(14)τg=fgρgvg22,τl=flρlvlf22,τi=fiρgvg−vlf22.During the filling, the air pocket is compressed under the effect of the incoming flow upstream and liquid slug downstream and after all entrapped at a stagnant position, i.e.,vg=0. In this way, (13) and (14) can be further written as(15)ρl−ρggsinθ1−τlSlA1−Hg−τiSi1A1−Hg+1AHg=0,(16)τg=0,τl=flρlvlf22,τi=fiρgvlf22.(4) Momentum Equation for the Liquid Slug. Since the compression process of air pocket is directly related to the motion of liquid slug, the slug is hereby specified as the controlled volume for which a momentum equation is established. Based on the fact that the resultant force that acts on the controlled volume is equal to the difference of liquid momentum per unit time by subtracting those flowing in the controlled volume from those flowing out, the following equation can be established:(17)Pg−PaA−ΔPfh1sinθ1+h2sinθ2−ρlAh2−h1g=ddtρlAh1sinθ1+h2sinθ2u−ρlQlvlf,where ΔPf is the friction force per unit length and the expression is specified as(18)ΔPf=18πDρλu2,And u is the speed at which the liquid slug advances upward in the upward slope, and it can be calculated using(19)u=1sinθ2dh2dt.During the growth of liquid slug toward the upstream in the downward sloping reach, the liquid rises at a rate equal to the shortening of air pocket:(20)1sinθ1dh1dt+dLgdt=0.


### 2.2.2. Closure Models

(1) Geometric Relationship. Refer to Figure 1; for the air pocket, the following hold: (21)Ag=D24π−δ+12sin2δ,(22)Sg=Dπ−δ,(23)Si=Dsinδ.Based on the above, the void fraction can be calculated as follows:(24)Hg=1−1πδ−12sin2δ.(2) Frictional Coefficient. Since the air pocket is nearly stagnant throughout the filling process, the air frictional coefficient fg is therefore neglected.Spedding and Hand [16] used (25) to calculate the the liquid frictional coefficient between the liquid and pipe wall:(25)fl=24Relf−1,Relf<2100,0.0262HlResl−0.139,Relf>2100.Among which,(26)Relf=Dvlfρlμl,Resl=Dvslρlμl,Hl=1−Hg.Xiao et al. [17] commented on 16 computational methods for interphase frictional coefficient fi calculation after testing them with the multiphase flow database of Tulsa University. The Andritsos–Hanratty correlation is after all recommended when D≤0.127m. When D>0.127m, the corrected Duns and Ros correlation by Baker is preferred. In the situation where the flow regime is to be determined, the recommended interphase frictional coefficient by Shoham and Taitel [18] is adopted, i.e., fi=0.0142, which is the preferable average for interphase wave of small amplitude and is applicable to a diversity of pipe diameters. According to the field pressure drop data of AGA large-diameter pipeline, the computational method for fi recommended by Xiao is well supported.


## 2.2.1. Governing Equations

After the air pocket is formed, the liquid slug at the valley point will start growing towards both upstream and downstream sloping reaches, resulting in the compression of air pocket. In this stage, the entrainment of air pocket is not considered. By specifying the V-shaped segment of gas–liquid two-phase flow as the object of study, the mathematical model for liquid slug growth and air pocket compression is established as shown in Figure2.Figure 2
Schematic diagram of liquid slug growth and air pocket compression.(1) Air Continuity Equation. Since throughout the water filling process, only liquid phase is pumped in from the inlet without any addition of air to the flow, thus, there is no mass exchange between the air pocket and the external environment. By specifying the air pocket segment as the controlled volume, the air continuity equation can be expressed as follows:(5)dmgdt=0.According to the assumptions, the air pocket meets the ideal gas equation. Thus, the mass of air pocket in the controlled volume can be expressed as(6)mg=MgARTPgLgHg.By substituting the above into (5), it can be deduced that(7)dPgLgHgdt=0.(2) Liquid Continuity Equation. By specifying the liquid in the V-shaped segment (the blue section in Figure 2) as the controlled volume, the liquid continuity equation is then given by(8)dmldt=ρlQl.The mass of liquid within the controlled volume is given by (9) based on the assumptions of D/h1≪1 and D/h2≪1:(9)ml=ρlAL1−LgHg+h2sinθ2.By substituting (9) into (8), it can be deduced that(10)1sinθ2dh2dt−dLgHgdt=vsl.(3) Momentum Equation for the Stratified Flow. The most probable flow pattern that could occur along the downward slope is the stratified flow. The quasisteady state force balance equation is used as the momentum equation:(11)ρgAggsinθ1−τgSg+τiSi−AgdPgdx=0,(12)ρlAlgsinθ1−τlSl−τiSi−AldPgdx=0.Combining the two equations above by eliminating the pressure gradient term, the combined momentum equation is then deduced as(13)ρl−ρggsinθ1−τlSlA1−Hg+τgSgAHg−τiSi1A1−Hg+1AHg=0.Among which,(14)τg=fgρgvg22,τl=flρlvlf22,τi=fiρgvg−vlf22.During the filling, the air pocket is compressed under the effect of the incoming flow upstream and liquid slug downstream and after all entrapped at a stagnant position, i.e.,vg=0. In this way, (13) and (14) can be further written as(15)ρl−ρggsinθ1−τlSlA1−Hg−τiSi1A1−Hg+1AHg=0,(16)τg=0,τl=flρlvlf22,τi=fiρgvlf22.(4) Momentum Equation for the Liquid Slug. Since the compression process of air pocket is directly related to the motion of liquid slug, the slug is hereby specified as the controlled volume for which a momentum equation is established. Based on the fact that the resultant force that acts on the controlled volume is equal to the difference of liquid momentum per unit time by subtracting those flowing in the controlled volume from those flowing out, the following equation can be established:(17)Pg−PaA−ΔPfh1sinθ1+h2sinθ2−ρlAh2−h1g=ddtρlAh1sinθ1+h2sinθ2u−ρlQlvlf,where ΔPf is the friction force per unit length and the expression is specified as(18)ΔPf=18πDρλu2,And u is the speed at which the liquid slug advances upward in the upward slope, and it can be calculated using(19)u=1sinθ2dh2dt.During the growth of liquid slug toward the upstream in the downward sloping reach, the liquid rises at a rate equal to the shortening of air pocket:(20)1sinθ1dh1dt+dLgdt=0.


## 2.2.2. Closure Models

(1) Geometric Relationship. Refer to Figure 1; for the air pocket, the following hold: (21)Ag=D24π−δ+12sin2δ,(22)Sg=Dπ−δ,(23)Si=Dsinδ.Based on the above, the void fraction can be calculated as follows:(24)Hg=1−1πδ−12sin2δ.(2) Frictional Coefficient. Since the air pocket is nearly stagnant throughout the filling process, the air frictional coefficient fg is therefore neglected.Spedding and Hand [16] used (25) to calculate the the liquid frictional coefficient between the liquid and pipe wall:(25)fl=24Relf−1,Relf<2100,0.0262HlResl−0.139,Relf>2100.Among which,(26)Relf=Dvlfρlμl,Resl=Dvslρlμl,Hl=1−Hg.Xiao et al. [17] commented on 16 computational methods for interphase frictional coefficient fi calculation after testing them with the multiphase flow database of Tulsa University. The Andritsos–Hanratty correlation is after all recommended when D≤0.127m. When D>0.127m, the corrected Duns and Ros correlation by Baker is preferred. In the situation where the flow regime is to be determined, the recommended interphase frictional coefficient by Shoham and Taitel [18] is adopted, i.e., fi=0.0142, which is the preferable average for interphase wave of small amplitude and is applicable to a diversity of pipe diameters. According to the field pressure drop data of AGA large-diameter pipeline, the computational method for fi recommended by Xiao is well supported.


## 2.3. Air Pocket Entrainment

As the liquid slug grows, the liquid level heights both in the downward and in the upward sloping reaches continue to increase accordingly. This paper assumes that the highest liquid level of upward sloping reach, i.e., the waterfront arrives at the upward summit, marks the end of the slug growth. At this moment, if the liquid level height of the downward sloping reach is lower than that of the upward one, the difference of water levels in upstream and downstream reaches induces backpressure. Furthermore, assuming that only the excess of the backpressure above the air pocket pressure is responsible for the air pocket entrainment, the criterion for the onset of entrainment is thus given by(27)Pb=ρlgL2sinθ2−h1′>Pg.Once entrainment occurs, due to the intaking of large numbers of entrained small bubbles from the pocket tail, the fluid between pocket tail and downstream liquid slug front starts to occur in a form of aerated liquid slug as shown in Figure3.Figure 3
Schematic diagram of air entrainment at the pocket tail.In the figure, “T” represents the tail of air pocket. Bubbles are entrained from the pocket tail at the rate of Φen. Meanwhile, part of the entrained bubbles will flow back and recoalesce with the pocket tail. However, considering the fact that the elongated air pockets in long-distance pipeline are always with great length (normally several kilometers), the amount of entrained small bubbles is usually much more than that of the flow back. Therefore, the effect of bubbles flow back on the entrainment is not considered in this study. All of these bubbles entrained at the rate of Φen will enter the liquid slug downstream to form a bubbly flow.Brauner and Ullmann [19] assume that the entrainment rate Φen is dependent on the balance between the rate of turbulent kinetic energy production and the rate of bubble surface energy production and subsequently deduced that(28)Φenvsl=1400CJd˜critWe−Wec,where(29)We=ρlDvlf−vsl2σ,(30)Wec=100C′d˜crit.In fact,CJ and C′ are tunable constants, where CJ can be tuned based on turbulence intensity and critical bubble size. It is suggested by Brauner and Ullmann that CJ=1,C′=2/3. And this value is found to be in agreement with the experimental data of Nydal and Andreussi [20] in a nearly horizontal 5 cm pipe.The critical size of deformable bubbles is suggested by Brodkey [21]:(31)d˜crit=0.4σρl−ρggcosθ1′D21/2=0.224cosθ1′EoD,EoD≥0.2,0.25,EoD<0.2,Among which,(32)EoD=ΔρgD28σ,θ1′=θ1,θ1≤45∘,90−θ1,θ1>45∘.Regardless of the possible coalescence and buildup of entrained small bubbles in the slug under the effect of buoyancy, it is considered that the entrained small bubbles will simply mix with the liquid downstream to form a bubbly flow. The bubbly flow in the slug region can be described by the drift-flux model. The bubble velocity in liquid slug at pocket tailvblsT is given by(33)vblsT=C0Tvm−vdT,where C0T is the velocity distribution parameter in the near bubble tail region. The rigorous computation of C0T requires the bubble velocity distributions in these regions, which are quite difficult to obtain. However, wall peaking in void fraction distributions generally results in C0T<1, whereas C0T>1 characterizes center peaking. In the case of horizontal or inclined pipeline especially when the liquid velocity is relatively low, C0T<1 is generally considered. This is because the liquid carrying capacity is relatively weak in this situation and the entrained bubbles will normally gather around the upper wall under the buoyancy effect. According to the experimental results by Hernandez and Fabre [22], the recommended value of C0T as 0.5 is adopted in this paper.vm is the velocity of mixture flux which is calculated as follows:(34)vm=vsl+Φen.vdT is the drift velocity. By using the Harmathy model corrected by Wallis [23] for the effect of pipe inclination, the drift velocity is given as follows:(35)vdT=1.53σgρl−ρgρl20.25sinθ11−HglsT1.5.According to the mass conservation, the gas void fraction in liquid slug at the pocket tailHglsT is(36)HglsT=ΦenvblsT.Due to the air entrainment, the pressure, length, and mass of the air pocket decreases constantly. The pressure during the entrainment is given by(37)Pgk+1=Pgk−ΔPg,where k is the current calculation time step and k+1 is the next one. ΔPg is obtained by integrating (11).The change rate of mass of entrained air can be calculated by(38)dmgendt=−ρgTAHglsTvblsT,where ρgT can be deduced from the pressure Pgk+1. And then the remaining air mass and air length after the air pocket is entrained can thus be expressed as(39)mgk+1=mgk−Δmgen,(40)Lgk+1=1−ΔmgenmgkLgk.


## 2.4. Computational Procedure

By specifying the time when the liquid phase is pumped into the pipeline from the inlet as the simulation start time, the initial conditions at the moment of air pocket taking shape are given by(41)Hg0=1−1πδ0−12sin2δ0,(42)Pg0=Pa,(43)Lg0=L1,(44)ρg0=ρa,(45)mg0=ρg0AHg0Lg0,(46)h10=h20=D,(47)u0=Ql2A,where δ0 is iteratively solved by the open-channel flow model in Section 2.1. Then, the air pocket compression model in Section 2.2 is initiated to conduct the calculation. The five ordinary differential equations (ODE set) of (7), (10), (17), (18), and (19) are solved using the forward difference method to obtain Pg,Lg,h1,h2,u and thenHg by (15). The waterfront arriving at the summit of the upward sloping reach marks the end of the slug growth. At this moment, whether the criterion for the onset of entrainment, i.e., (27), is satisfied is judged. If it does not meet the criterion, the simulation stops. Otherwise, the entrainment model in Section 2.3 is activated to continue the calculation as follows. At each calculation time step, (28)–(32) are firstly replaced by the previously solved variables such as Pg,Lg,Hg and other required information as input to obtain Φen and then vm by (34). Secondly, based on the calculated Φen and vm, the remaining three implicit algebraic equations, i.e., (33), (35), and (36), are iteratively solved to obtain HglsT,vblsT,vdT. Finally, the pressure during the entrainment Pgk+1, the remaining air mass mgk+1, and the remaining air length Lgk+1 are deduced by (37), (39), and (40), respectively. After the current time step is finished, the above entrainment process is repeated again to calculate the next time step until the end of the simulation. The flowchart for the entire computational procedure is shown in Figure 4. The time steps were taken as 10 s and 300 s for compression process and entrainment process, respectively.Figure 4
Flow chart of the entire computational procedure.


## 3. Model Verification

Two domestic hilly pipelines which have already been filled with water are chosen as test pipelines to verify the model proposed in this paper. For the pipeline routing, refer to Figure5; for mileage and elevation data, refer to Table 1. The pipeline diameters are 559 mm and 508 mm, respectively. A constant inflow of water at the rate of 600 m3/h and 900 m3/h, respectively, is considered for the two test pipelines.Figure 5
Test pipeline routing: (a) pipeline 1 and (b) pipeline 2.
(a)(b)Table 1
Mileage and elevation data of test pipelines.
PipelineLocationLength (km)Elevation (m)Pipeline 1Initial station0751Peak point116.611345.4Valley point (valve chamber 1)208.91432.7Peak point (valve chamber 2)295.8798Terminal station535.6705Pipeline 2Initial station01921.2Peak point31.82589.03Valve chamber 391.412049.58Valley point145.451560.54Valve chamber 4290.722276.7Peak point402.842827.81Terminal station438.452198.8During water filling the pipelines, four valve chambers, i.e., valve chambers 1 and 2 in pipeline 1 and valve chambers 3 and 4 in pipeline 2, were deployed along the line to monitor the real-time pressure. The monitored pressure values were transmitted to central control room via SCADA. To make a comparison between the field data and the simulation results, the data recorded by SCADA system was exported. The comparison of simulation results and field data is shown in Figure6. It can be seen that the simulated value and measured value are basically consistent. Although some of the data have relatively larger deviations (the maximum deviation does not exceed ±20%), for engineering practice, the overall result is within an acceptable range. That is to say, the mathematical model herein can be used to simulate the evolution of air pocket with good accuracy.Figure 6
Comparison of simulation results and field data: (a) valve chamber 1 in pipeline 1, (b) valve chamber 2 in pipeline 1, (c) valve chamber 3 in pipeline 2, and (d) valve chamber 4 in pipeline 2.
(a)(b)(c)(d)


## 4. Simulation Results

Taking the test pipeline 1 in Section3 as an example, Figure 7 reveals the simulated evolution with time of backpressure Pb, air pocket pressure Pg, and liquid level heights in both downward and upward sloping reaches h1 and h2. It can be seen that the air pocket is formed at about 55.9 h. Since the upstream pressure is not strong enough to push the liquid slug unidirectionally towards the downstream, the slug grows towards both the upstream and downstream, resulting in the gradual increases of liquid level heights both in downward and in upward sloping reach. The air pocket becomes compressed and the pressure increases consequently. At about 96.45 h, the waterfront arrived at the summit of upward sloping reach, i.e., the peak point downstream in Figure 5, which marks the end of the slug growth process. At this moment, the upward and downward liquid level heights are 365.3 m and 261.15 m, respectively, which consequently induces a backpressure of about 1.02 MPa. Comparison with the current gas pressure of 0.14 MPa indicates that the criterion for the onset of air entrainment (i.e., (27)) is exactly satisfied. Therefore, entrainment starts to occur at the tail of air pocket. A large number of small bubbles are entrained from the pocket tail rapidly into the liquid slug downstream and the air pocket pressure decreases constantly. At the same time, the rising liquid level in downward sloping reach leads to a gradual decrease in backpressure. At the time of 641.11 h, the backpressure drops to equal the gas pressure; the entrainment stops along with the end of the whole air removal simulation.Figure 7
Evolution with time ofh1,h2,Pb, and Pg: ① stage 1: air pocket formation; ② stage 2: air pocket compression; and ③ stage 3: air pocket entrainment.In pipeline water filling, field operators are also concerned about the length and location of air pocket in addition to pipeline pressure. This is because the length and location of air pocket plays a critical role in deciding air removal operations. Figure8 gives the evolution with time of air pocket length and mass. At the first stage before the formation of air pocket, the length and mass are all kept constant at 0. After the air pocket is formed and compressed, the length of air pocket decreases gradually but the mass remains constant since no mass exchange exists between the air pocket and the external environment or the liquid film. Then, at the third stage, air entrainment starts to occur, and a large number of small bubbles are entrained from the air pocket to the liquid slug downstream. The length and mass of air pocket decrease consequently. By the end of the entrainment, the length of air pocket is reduced to 56.63 km from the initial 92.3 km. If no further actions are taken, the air pocket with an approximate length of 57 km will not be removed and still remain in the pipeline.Figure 8
Evolution with time of air pocket length and mass.The development of air pocket described above is also reflected in the waterfront movement diagram in Figure9. From the initial station to peak point upstream, the water flows downstream in the form of full pipe flow with a constant velocity of 0.72 m/s. After the waterfront arrives at the peak point upstream, it flows downward to the valley point in the form of slack line flow, where the flow velocity reaches as high as 2.376 m/s. After arriving at the valley point, the water firstly starts accumulating in the elbow and gradually a liquid slug is formed, resulting in the slug velocity drop to a very low value. From the time of about 55.9 h when the air pocket is formed, as the liquid level in the upward sloping reach rises constantly, the slug velocity begins to increase rapidly and remains basically constant at about 0.6 m/s afterwards. After the waterfront arriving at the peak point downstream, the compression of air pocket stops and air entrainment starts. The velocity of the full pipe flow is then accelerated to the initial value, i.e., 0.72 m/s. The same conclusion can also be drawn by observing the tangent lines of the curve of waterfront location variation.Figure 9
Movement of the waterfront.


## 5. Discussion

The terrain category and inlet flow rate of pipeline are thought of having the most important effect on the removal of air pocket. Herein, a dedicated analysis is performed on the two factors as follows.


### 5.1. Effect of Terrain Category

Based on the topographic characteristics of upward and downward sections, terrains are mainly divided into three groups, i.e., terrains 1∼3 shown in Figure10(a). For the corresponding topographic parameters, refer to Table 2. A series of simulations are conducted on the three types of terrains with pipeline inlet flow rate assumed to be 900 m3/h constant. For the simulation results of pressure, length, and mass changes of air pocket, refer to Figures 10(b)–10(d).Figure 10
Three types of V-shaped topographies (a) and simulation results of air pocket: (b) pressure change, (c) length change, and (d) mass change.
(a)(b)(c)(d)Table 2
Topographic parameters for the three types of terrains.
No.Downward sloping reachUpward sloping reachLength (km)Elevation difference (m)Inclination (°)Length (km)Elevation difference (m)Inclination (°)Terrain 110.5612877.0011.936923.33Terrain 25.271241.3510.022901.66Terrain 31.1556.492.8228.89250.830.50From Figures10(b)–10(d), it can be seen that terrain category will have a great influence on the removal of air pocket, which in detail can be summarized as follows:①
Terrain 1: the air pocket is compressed and partly entrained, i.e., the air pocket cannot be totally removed, and a certain amount still remains in the downward sloping reach. From Table2, it can be seen that the upward sloping reach is almost having the same length as the downward one, but with a much smaller elevation difference. Due to this, the water level in the downward sloping reach still keeps rising even after the level in the upward has reached the high elevation point. Consequently, the backpressure causing the entrainment of air pocket decreases. When the backpressure falls below the air pocket pressure, air entrainment ceases. That is why a segment of air pocket with an approximate length of 4.65 km (the initial length of air pocket is 10.56 km) still remains and cannot be removed. For the situation which is met with in Section 4 of this paper, the terrain just falls into Terrain 1②
Terrain 2: the air pocket is compressed and completely removed. In this situation, the condition required for air entrainment is always present. The air pocket with an initial length of 5.27 km is totally entrained into the liquid slug downstream and flows out of the pipeline together with the liquid. It is observed in the simulation that the downward sloping reach will eventually be filled with water and undergo a full pipe flow③
Terrain 3: the air pocket is only compressed without any removal, i.e., after the simulation, only the length of air pocket is shortened, while mass remains unchanged. In Terrain 3, the upward sloping reach is with a very small inclination which is merely 0.5° and can be seen as approximately horizontal. When the liquid slug grows in this sloping reach, it is almost impossible to gain enough backpressure to initiate the entrainment of air pocket. In this situation, the air pocket can only be compressed without any entrainment. After the waterfront arrives at the summit of the upward sloping reach, the simulation stops. The air pocket is compressed to a final pressure of 1.89 MPa.


### 5.2. Effect of Inlet Flow Rate

For different flow rates under Terrain 3, the variations of pressure and length of air pocket during the compression process are shown in Figure11 (“Q” indicates the inlet flow rate). It is obvious that, with the increase of flow rate, the peak pressure (the maximum pressure) increases and the remaining length reduces. Since excessive pressure poses a potential threat to the pipeline, Terrain 3 is hereby defined as the “Dangerous Terrain.” For pipelines located in such terrains as Terrain 3, the inlet flow rate shall be strictly controlled during the filling process.Figure 11
Simulation results of Terrain 3 for different flow rates: (a) variation of air pocket pressure and (b) variation of air pocket length.
(a)(b)In respect of the peak pressure and remaining length of air pocket, Figures12(a) and 12(b), respectively, give the comparison thereof among Terrain 3, Terrain 2, and Terrain 1. It can be concluded that if the terrain falls out of the “Dangerous Terrain” category, i.e., Terrain 1 and Terrain 2, the peak pressure of air pocket becomes relatively lower and less subject to the inlet flow rate. Basically, it remains almost unchanged with the change of inlet flow rate.Figure 12
Simulation results at different flow rates for three types of terrains: (a) peak pressure of air pocket and (b) remaining length of air pocket.
(a)(b)In addition, from the comparison between Terrain 1 and Terrain 3, it is observed that, with the increase of inlet flow rate, the remaining air pocket lengths of both terrains reduce. However, the mass of air pocket of Terrain 3 is constant while that of Terrain 1 decreases gradually (refer to Figure13 where the solid line represents the initial mass of air pocket and the dotted line represents the mass remained at the end of simulation). For Terrain 1, more air is entrained with the increase of flow rate, which implies that a relatively larger flow rate is more beneficial to the removal of air pocket.Figure 13
Initial and remaining air mass at different flow rates for three types of terrains.From Figures12(b) and 13, it is known that under different flow rate conditions the air pocket in Terrain 2 is always completely removed. The remaining length and mass of air pocket both approach 0. In respect of the timing when air pocket is formed, the timing when liquid slug ends growing, the timing when air entrainment ceases, the duration of liquid slug growth, the duration of air entrainment, and the duration of air removal (slug growth + air entrainment), Figure 14 gives the description of the variation thereof under different flow rate conditions. It can be concluded that the larger the flow rate, the earlier the formation of air pocket, the faster the growth of liquid slug, and the earlier the entrainment of air pocket. During the entrainment process, the turbulent flow at the air pocket tail becomes more violent with the increase of flow rate, thus resulting in a much more rapid entrainment rate and shorter entrainment duration. Similar to Terrain 1, an increased flow rate is also beneficial to the removal of air pocket in Terrain 2.Figure 14
Simulation results at different flow rates for Terrain 2: (a) timing variation of each stage and (b) duration variation of each stage.
(a)(b)


## 5.1. Effect of Terrain Category

Based on the topographic characteristics of upward and downward sections, terrains are mainly divided into three groups, i.e., terrains 1∼3 shown in Figure10(a). For the corresponding topographic parameters, refer to Table 2. A series of simulations are conducted on the three types of terrains with pipeline inlet flow rate assumed to be 900 m3/h constant. For the simulation results of pressure, length, and mass changes of air pocket, refer to Figures 10(b)–10(d).Figure 10
Three types of V-shaped topographies (a) and simulation results of air pocket: (b) pressure change, (c) length change, and (d) mass change.
(a)(b)(c)(d)Table 2
Topographic parameters for the three types of terrains.
No.Downward sloping reachUpward sloping reachLength (km)Elevation difference (m)Inclination (°)Length (km)Elevation difference (m)Inclination (°)Terrain 110.5612877.0011.936923.33Terrain 25.271241.3510.022901.66Terrain 31.1556.492.8228.89250.830.50From Figures10(b)–10(d), it can be seen that terrain category will have a great influence on the removal of air pocket, which in detail can be summarized as follows:①
Terrain 1: the air pocket is compressed and partly entrained, i.e., the air pocket cannot be totally removed, and a certain amount still remains in the downward sloping reach. From Table2, it can be seen that the upward sloping reach is almost having the same length as the downward one, but with a much smaller elevation difference. Due to this, the water level in the downward sloping reach still keeps rising even after the level in the upward has reached the high elevation point. Consequently, the backpressure causing the entrainment of air pocket decreases. When the backpressure falls below the air pocket pressure, air entrainment ceases. That is why a segment of air pocket with an approximate length of 4.65 km (the initial length of air pocket is 10.56 km) still remains and cannot be removed. For the situation which is met with in Section 4 of this paper, the terrain just falls into Terrain 1②
Terrain 2: the air pocket is compressed and completely removed. In this situation, the condition required for air entrainment is always present. The air pocket with an initial length of 5.27 km is totally entrained into the liquid slug downstream and flows out of the pipeline together with the liquid. It is observed in the simulation that the downward sloping reach will eventually be filled with water and undergo a full pipe flow③
Terrain 3: the air pocket is only compressed without any removal, i.e., after the simulation, only the length of air pocket is shortened, while mass remains unchanged. In Terrain 3, the upward sloping reach is with a very small inclination which is merely 0.5° and can be seen as approximately horizontal. When the liquid slug grows in this sloping reach, it is almost impossible to gain enough backpressure to initiate the entrainment of air pocket. In this situation, the air pocket can only be compressed without any entrainment. After the waterfront arrives at the summit of the upward sloping reach, the simulation stops. The air pocket is compressed to a final pressure of 1.89 MPa.


## 5.2. Effect of Inlet Flow Rate

For different flow rates under Terrain 3, the variations of pressure and length of air pocket during the compression process are shown in Figure11 (“Q” indicates the inlet flow rate). It is obvious that, with the increase of flow rate, the peak pressure (the maximum pressure) increases and the remaining length reduces. Since excessive pressure poses a potential threat to the pipeline, Terrain 3 is hereby defined as the “Dangerous Terrain.” For pipelines located in such terrains as Terrain 3, the inlet flow rate shall be strictly controlled during the filling process.Figure 11
Simulation results of Terrain 3 for different flow rates: (a) variation of air pocket pressure and (b) variation of air pocket length.
(a)(b)In respect of the peak pressure and remaining length of air pocket, Figures12(a) and 12(b), respectively, give the comparison thereof among Terrain 3, Terrain 2, and Terrain 1. It can be concluded that if the terrain falls out of the “Dangerous Terrain” category, i.e., Terrain 1 and Terrain 2, the peak pressure of air pocket becomes relatively lower and less subject to the inlet flow rate. Basically, it remains almost unchanged with the change of inlet flow rate.Figure 12
Simulation results at different flow rates for three types of terrains: (a) peak pressure of air pocket and (b) remaining length of air pocket.
(a)(b)In addition, from the comparison between Terrain 1 and Terrain 3, it is observed that, with the increase of inlet flow rate, the remaining air pocket lengths of both terrains reduce. However, the mass of air pocket of Terrain 3 is constant while that of Terrain 1 decreases gradually (refer to Figure13 where the solid line represents the initial mass of air pocket and the dotted line represents the mass remained at the end of simulation). For Terrain 1, more air is entrained with the increase of flow rate, which implies that a relatively larger flow rate is more beneficial to the removal of air pocket.Figure 13
Initial and remaining air mass at different flow rates for three types of terrains.From Figures12(b) and 13, it is known that under different flow rate conditions the air pocket in Terrain 2 is always completely removed. The remaining length and mass of air pocket both approach 0. In respect of the timing when air pocket is formed, the timing when liquid slug ends growing, the timing when air entrainment ceases, the duration of liquid slug growth, the duration of air entrainment, and the duration of air removal (slug growth + air entrainment), Figure 14 gives the description of the variation thereof under different flow rate conditions. It can be concluded that the larger the flow rate, the earlier the formation of air pocket, the faster the growth of liquid slug, and the earlier the entrainment of air pocket. During the entrainment process, the turbulent flow at the air pocket tail becomes more violent with the increase of flow rate, thus resulting in a much more rapid entrainment rate and shorter entrainment duration. Similar to Terrain 1, an increased flow rate is also beneficial to the removal of air pocket in Terrain 2.Figure 14
Simulation results at different flow rates for Terrain 2: (a) timing variation of each stage and (b) duration variation of each stage.
(a)(b)


## 6. Conclusions

By giving a comprehensive consideration of air pocket formation, compression, and entrainment, a mathematical method of air pocket evolution during water filling a long-slope pipeline is proposed in this paper. According to the simulation results, it is suggested that after the formation of air pocket, liquid slug starts to grow toward both the upstream and downstream direction, resulting in the air pocket being compressed, pressure being increased, and length being reduced. Under the effect of the backpressure, the air pocket is entrained into small bubbles with its pressure, length, and mass being reduced gradually.The most important factors that affect the removal of air pocket are the terrain category and the inlet flow rate. As revealed by the study herein, the removal of air pocket mainly leads to three kinds of scenarios with each corresponding to one terrain category: the compressed and partly removed scenario, the compressed and completely removed scenario, and the compressed only without any removal scenario. In this paper, the terrain which leads to the last outcome is called the “Dangerous Terrain.” And for the “Dangerous Terrain,” it is of great importance that the inlet flow rate should be strictly confined within a certain level. While for the other two categories of terrains, an increased flow rate is in any respect beneficial to the removal of air pocket.Considering the above, it is suggested that for engineering practice an analysis on the terrain conditions shall be firstly conducted. After the terrain category has been identified, the air removal condition can be determined according to the findings herein. After all these, an appropriate range of inlet flow rate can be specified so as to achieve the maximum removal of air pocket.



---

*Source: 3157850-2020-07-27.xml*