## Abstract

The entrapped air pockets in a pressurized water supply pipe may produce a significant transient pressure, resulting in pipe deformation or even pipe rupture. This article focuses on the maximum pressure of entrapped air pockets in a rapid filling pipe and its influencing factors. In this paper, the rigid water column theory and the ideal gas equation are applied to the air vessel-pressured pipe–air pocket system, and the theoretical formula representing the relationship between the maximum air pressure (MAP) and air comprehensive coefficient is presented. In addition, a numerical model of the system is established by the method of characteristics. The rapid filling process is simulated according to the experimental setup and the results are in good agreement with the experimental data. The influence of the air comprehensive coefficient on the MAP is also studied. Both the theoretical formula and numerical simulation show that the smaller air comprehensive coefficient will result in larger MAP (i.e., the MAP decreases with the increase of the polytropic exponent and initial air pressure, and increases with the increase of the initial air length).

## NOMENCLATURE

### Symbols

*Z*_{u}Water level of upper steel air vessel

*Z*_{t}Top elevation of air chamber

*P*Pressure of the air pocket

*P*_{u}Air pressure in steel air vessel

*γ*Unit weight of water

*z*Absolute water level variation in air chamber

*L*Length of the diversion pipeline

*m*Polytropic exponent

*Q*Flow rate in the diversion pipe

*α*Total head loss coefficient in the diversion system

*A*Cross-sectional areas of the diversion pipe

*F*Cross-sectional areas of the air chamber

*g*Gravitational acceleration

*l*_{0}Initial length of the air pocket

*P*_{0}Initial absolute pressure of the air pocket

*σ*Air comprehensive coefficient

*f*Darcy–Weisbach friction factor

*θ*Angle of the pipeline inclined with the horizontal

*x*Distance along the pipeline

*ξ*Head loss coefficient of the spherical valve

*z*_{c}Stable water level variation in air chamber

*z*_{m}Maximum water level variation in air chamber

### Acronym

## INTRODUCTION

Air bubbles are easily entrapped in pipes forming air pockets in undulating water supply systems. If the pipeline system cannot discharge the air before the opening of the valve, the phenomenon of pressurized water flow-impacted entrapped air pockets in the pipeline system will inevitably appear. The pressure generated in this process is large and may cause damage to the pipeline system.

The prediction of pressure transient process in pipes with air is complicated, as it is related to air pocket size and movement, flow velocity and depth, pipe size and slope, venting and boundary conditions. The influence of entrapped air on abnormal transient pressures is often ambiguous because the compressibility of the air pocket permits the flow to accelerate but also partly cushions the system; the balance of these tendencies have been associated with the initial void fraction of the air pocket (Zhou *et al.* 2011). Compared with a single-phase water flow, free air increases the elasticity of the water and reduces the wave velocity, thus reducing the pressure pulsation during the transient process (Wylie *et al.* 1993). Experimental results show that an air pocket in pipeline system will produce a great pressure fluctuation in the hydraulic transient process (Vasconcelos & Wright 2011). For rapid filling scenarios, hydraulic bores, air–water interactions, and boundary reflections may lead to air pocket entrapment (Hamam & McCorquodale 1982; Li & McCorquodale 1999), which in turn can influence the transient pressure in the system because of compression, expansion, migration, and expulsion of the air pocket (Martin 1976; Aimable & Zech 2003; Zhou *et al.* 2011; Ferreri *et al.* 2014; Besharat *et al.* 2016). Many researchers have studied the rapid filling process in pipes with open ends, closed ends, orifice plates, or small air valves. Liou & Hunt (1996) analyzed the velocity variation of water in a rapid filling process, undulating, empty pipe with an open end. Zhou *et al.* (2002) presented experimental measurements from a horizontal pipeline in which the pipe length was constant and the inlet pressure was more than twice the atmospheric pressure; the work shows that when less air was present, the maximum pressure of the air pocket increased as the cushioning effect of the air pocket decreased. Chaudhry & Reddy (2011) analyzed the initial filling of an empty tunnel with a closed gate at the end and without an air vent before the gate. Zhou *et al.* (2011) carried out experiments to study pipes containing one or more air pockets and set up different ports for the tail of the pipe. Vasconcelos & Wright (2011) presented a model based on lumped inertia analysis which was developed by Lewis & Wright (2012). These models were successful in representing the experimental conditions presented by the authors, but failed in attempting to describe the severity of actual geysering episodes. Hou *et al.* (2013) investigated the filling velocity in rapid filling pipes with an open end. Martins *et al.* (2017) used a three-dimensional computational fluid dynamics model to simulate the rapid filling process and describe the maximum pressures of air pockets for a given range of initial conditions and compared the results with the experimental data. Muller *et al.* (2017) investigated air-related geysers that provide insight into the mechanisms of air release and the displacement of water in vertical shafts.

In a pipeline system with entrapped air pockets, the hydraulic transients between different conditions is complicated as the air pressure is a nonlinear restoring force, and the influencing factors on the maximum air pressure (MAP) and their relationship are unclear. Particularly, the analytical expression that addresses the effect of influencing factors on the MAP is not sufficient, which is more general than specific numerical simulation. In this paper, the rigid water column theory and the ideal gas equation are applied to the pipeline system, and the theoretical formula representing the relationship between the maximum pressure of air pocket and air comprehensive coefficient is deduced, thereby contributing to a more general understanding of the dynamic characteristics of entrapped air in a rapid filling pipe. In addition, the numerical model of the system is established based on the method of characteristics (Karney & Ghidaoui 1997; Miao *et al.* 2017; Liu *et al.* 2017), the rapid filling process is simulated based on the experimental setup and the results are compared with the experimental data. The influence of air comprehensive coefficient on the MAP is also studied. Numerical simulation of different operation conditions are carried out in order to clarify the influence of the polytropic exponent, the initial air pressure and the initial air length on the MAP. Finally, the conclusions can provide a reference for the design and operation of water supply systems.

## THEORETICAL ANALYSIS OF THE MAP AND ITS INFLUENCING FACTORS

### Mathematical model

The ideal gas equation shows that as the length of the air pocket reaches its minimum value, the air pocket pressure is at maximum. Therefore, the investigation on the air pocket pressure can be transformed into the study of the variation of the length of the air pocket, then the air pocket pressure can be obtained by the ideal gas equation. Martins (2012) carried out physical experiments to study the dynamic characteristics of an entrapped air pocket. In this paper, the pipeline system is established based on a physical experiment setup (Martins *et al.* 2017), and the layout of the pipeline system is shown in Figure 1. Pipe 3, which contains an air pocket, can be treated as an air chamber. In the process of the theoretical analysis on MAP, the following assumptions are made: (a) ignoring the water inertia in the air chamber, (b) the water in the pipeline is incompressible, and (c) the friction in the pipeline is quasi steady.

The basic equations of the pipeline system with an entrapped air pocket are presented as follows (Chaudhry 2014).

*Z*is the water level of the upper steel air vessel.

_{u}*Z*is the top elevation of the air chamber.

_{t}*P*and

*P*are the pressure of the air pocket and the air pressure in the steel air vessel, respectively.

_{u}*γ*is the unit weight of water.

*l*is the instantaneous equivalent air pocket length obtained by

*l*

*=*

*l*-

_{0}*z*, in which

*z*is the absolute water level variation in the air chamber, measured positive upwards.

*L*is the length of the pipeline.

*Q*is the discharge in the pipeline.

*α*is the total head loss coefficient in the diversion system, which is a constant during the transient process.

*F*and

*A*are the cross-sectional areas of the air chamber and diversion pipes, respectively.

*P*and

_{0}*l*are the initial absolute pressure and the length of the air pocket before the valve opens.

_{0}*m*is the polytropic exponent, taken as 1.0 in the isothermal process and 1.4 in the adiabatic process.

*g*is the acceleration due to gravity.

*t*is time. Subscript 0 refers to the initial steady value.

*σ*is defined as the air comprehensive coefficient,

*σ*

*=*

*1*

*+*

*mP*, and

_{0}/γl_{0}*μ*

*=*

*Z*

_{u}*+*

*P*

_{u}/γ-Z_{t}*+*

*l*

_{0}*−*

*P*. As the pressure on the left side of the spherical valve is greater than the right side as the valve is closed, the value of

_{0}/γ*μ*is greater than zero.

### Formula derivation

*φ*

*=*

*2gAF/L*,

*φ*is the constant greater than zero.

*C*is an integral constant.

_{1}*z*, where

_{c}*z*is the absolute value between the new steady water level and the initial chamber water level, and the actual value of

_{c}*z*is positive. Taking

_{c}*Q*= 0,

*z*=

*z*into Equation (8) gives:

_{c}*Q*

*=*0,

*z*

*=*

*z*), where

_{m}*z*is the absolute value between the maximum water level and the initial chamber water level, and the actual value of

_{m}*z*is positive. Taking

_{m}*Q*= 0,

*z*

*=*

*z*into Equation (10) gives: when the system reaches a new steady state, the discharge in the pipe becomes zero and the water level variation in the air chamber becomes

_{m}*z*. Taking

_{c}*Q*= 0,

*z*

*=*

*z*into Equation (5) can draw

_{c}*z*=

_{c}*μ*/

*σ*, then Equation (11) can be simplified to:

*ɛ*on both sides of Equation (13) gives: when the system reaches a new steady state, the water level variation in the air chamber will be less than the maximum water level variation (i.e.,

*z*

_{c}*−*

*z*< 0), therefore 1 − e

_{m}*≠ 0. The following formula is available:*

^{−φαzm}As illustrated in Equation (16), the air comprehensive coefficient, a parameter reflecting the initial state of the air pocket in the system, has an important effect on the maximum pressure of the air pocket during the transient process. In the case of keeping other parameters constant, the increase of the air comprehensive coefficient will lead to the decrease of MAP. Thus, the smaller air comprehensive coefficient will result in a larger MAP. It can be seen from the ideal gas equation and the expression of the air comprehensive coefficient that when the initial water level in the air chamber and other parameters is constant, the increase of the polytropic exponent and the initial air pressure will lead to the increase of the air comprehensive coefficient, resulting in a decrease in the MAP. The air comprehensive coefficient decreases with the increase of the initial air pocket length, leading to the increase of the MAP. Thus, keeping other parameters constant, larger initial air pocket length will result in a higher MAP.

## NUMERICAL INVESTIGATION

### Experimental model

An experiment was conducted by Martins *et al.* (2012) to study the maximum pressure of an air pocket in a pipeline system. The model layout diagram and the main parameters of the system are shown in Figure 1 and Table 1. The physical experimental rig mainly consists of the following parts: steel air vessel, water pipe, spherical valve, water pipe, closed end of vertical pipe with an air pocket. Pipes are made of polyvinyl chloride (PVC), with an inner diameter of 0.0536 m and a wall thickness of 4.7 mm. The spherical valve links the pressure steel air vessel and the water pipes with an inner diameter of 0.0536 m and opened totally in 0.23 s.

Pipe | Length/m | Diameter/m | Steel air vessel | |
---|---|---|---|---|

1 | 0.24 | 0.0536 | Diameter/m | 0.80 |

2 | 2.00 | 0.0536 | Volume/m^{3} | 1.00 |

3 | 1.38 | 0.0536 | Water Volume/m^{3} | 0.40 |

Pipe | Length/m | Diameter/m | Steel air vessel | |
---|---|---|---|---|

1 | 0.24 | 0.0536 | Diameter/m | 0.80 |

2 | 2.00 | 0.0536 | Volume/m^{3} | 1.00 |

3 | 1.38 | 0.0536 | Water Volume/m^{3} | 0.40 |

In the initial operating condition, the spherical valve is closed. The pressure on the left side of the valve is greater than the pressure on the right side, so there is a pressure difference between the two sides. When the valve is opened, the water flows from part A to part B through a vertical pipe. As the water level in part B rises, the air pocket is compressed and the air pocket pressure is increased. When the water level reaches the maximum in part B, the air pocket pressure reaches its maximum value as well. At this moment, the air pressure in part B is greater than the air pressure in the left steel air vessel. The water begins to flow in the opposite direction due to the pressure difference. Due to the presence of friction in the system, water level fluctuations are always diminishing, but the water level will reach a new stable state after a while. To understand the dynamical behaviors of the filling process and the water–air interaction, one-dimensional (1D) numerical simulation is conducted in the following sections.

### method of characteristics model

#### Pressurized pipe model

The following equations are used to describe the transient flow in pipes (Wylie *et al.* 1993).

*H*is the piezometric head.

*V*is the flow velocity.

*a*is the wave speed.

*g*is the acceleration due to gravity.

*f*is the Darcy–Weisbach friction factor.

*θ*is the angle of the pipeline inclined with the horizontal.

*x*is the distance along the pipeline measured positive in the downstream direction.

*D*is the diameter of the pipe.

*t*is time.

The momentum and continuity equations governing transient flow in closed conduits are classified as hyperbolic partial differential equations for which no analytical solutions are available for a general pipe system. The method of characteristics is widely used to solve the set of partial differential equations. The method of characteristics is the most commonly used numerical method for solving hydraulic transient process in a pipeline system. The method of characteristics is a kind of approximate calculation method for solving the hyperbolic partial differential equation based on the characteristic theory of the partial differential equation (Wylie *et al.* 1993). If the problem is simple, the analytical solution or approximate analysis solution can be solved by this method. If the problem is complicated, a numerical solution of high accuracy can be obtained. The method of characteristics has many advantages: (1) stability criteria can be established; (2) boundary conditions are easy to program; (3) it is applicable to all kinds of pipeline hydraulic transient analysis; and (4) it has better precision than all the other different methods.

#### Spherical valve model

*H*and

_{P1}*H*are the piezometric heads at each side of the spherical valve, respectively.

_{P2}*ξ*is the head loss coefficient of the spherical valve.

*Q*is the flow rate through the valve.

_{P}*A*is the area of the valve. The head loss coefficient

_{m}*ξ*under different valve opening

*τ*is provided by the manufacturer, and the diameter of the spherical valve used in this paper is 0.0536 m. Substituting compatibility equations (Wylie

*et al.*1993) into Equation (19) results in: in which the coefficients

*C*,

_{P}*B*,

_{P}*C*, and

_{M}*B*are constants, and

_{M}*Q*is the discharge of the previous time step.

_{P0}*Q*is the calculated discharge.

_{P}#### Mathematical model parameters

In this case, the method of characteristics is used to simulate the pipeline transition process with the closed air pocket under the condition of valve opening, and the results are compared with the experimental results. The layout of the numerical simulation system is the same as the experimental system, shown in Figure 1. The lengths of pipe 1, pipe 2, and pipe 3 are 0.24 m, 2.00 m, and 1.38 m, respectively. The diameter of all the pipes is 0.0536 m, and the pipe roughness is 0.012. The wave velocity of water in the pipeline is 450 m/s. The initial pressure of the air pocket is the atmospheric pressure (i.e., 10.33 m). The bottom elevation of the pipeline is 0.0 m. The spherical valve is closed under the initial condition and opens totally in 0.23 s with one-stage linear open law.

### Comparison with experimental results

The selected conditions are that the initial pressure of the air in the steel air vessel is 20.66 m and the initial air pocket length in part B is different under different conditions. The polytropic exponent is taken as 1.4. The maximum air pocket pressure difference and air pocket pressure variation are shown in Figure 2 and Table 2. It can be seen that although the maximum pressure of the air pocket is different in both cases, the differences are smaller relative to the MAP obtained by the physical experiment. In addition, the pressure variations of the air pocket are in agreement with each other in different conditions. Therefore, the numerical results have certain accuracy and can satisfy the accuracy requirements of the following research.

Condition | Maximum experiment pressure/m | Maximum simulation pressure/m | Pressure difference/m | Relative error/% |
---|---|---|---|---|

a | 36.4 | 36.0 | 0.4 | 1.10 |

b | 110.0 | 107.4 | 2.6 | 2.36 |

c | 130.6 | 124.5 | 6.1 | 4.67 |

Condition | Maximum experiment pressure/m | Maximum simulation pressure/m | Pressure difference/m | Relative error/% |
---|---|---|---|---|

a | 36.4 | 36.0 | 0.4 | 1.10 |

b | 110.0 | 107.4 | 2.6 | 2.36 |

c | 130.6 | 124.5 | 6.1 | 4.67 |

### Influencing factors on MAP

Under the guidance of the analytical expression between the air comprehensive coefficient and the maximum pressure of air pocket, this section investigates the effects of initial air pressure, polytropic exponent, and initial air length on the maximum pressure of the air pocket by using a control variate method. It is generally believed that the process of air transition in the pipes is between the adiabatic process and the isothermal process (i.e., m = 1.0 ∼ 1.4). Therefore, the polytropic exponents take *m* as 1.0, 1.2 and 1.4 for comparison, respectively. The initial air pocket pressure is 10.33 m and the initial air length is 0.102 m. The calculated results are shown in Figure 3.

As shown in Figure 3, with the other parameters remaining unchanged, the maximum air pocket pressure is decreasing with the increase of the polytropic exponent. The smaller polytropic exponent resulted in the larger MAP. When the polytropic exponent is taken as 1.4, the difference between the calculated maximum pressure of the air pocket and the experimental results is small. With the decrease of the polytropic exponent, the difference between the experimental results and the numerical results began to increase. When the polytropic exponent is taken as 1.0, the numerical simulation result of MAP is 125.02 m, which is 15.00 m higher than the experimental data. The polytropic exponent is an uncertain factor in the experiment, and the value of polytropic exponent in each experiment may not be the same due to the experimental instrument and the experimental environment. The MAP decreases with the increase of the polytropic exponent, which agrees with the theoretical analysis.

It can be seen from the theoretical analysis (Equation (16)) that the increase of the initial air pressure will change the value of the air comprehensive coefficient, affecting the MAP during the transient process. To verify this conclusion, the rapid filling process with different initial air pressure is simulated, and the effects of initial air pressure on maximum pressure of the air pocket are also analyzed. In this simulation, the polytropic exponent is taken as 1.4, the initial air length is taken as 0.102 m, and the initial water level in the air chamber is taken as the same value under different conditions. The initial air pressures are 10.33, 12.33, 14.33, and 16.33 m, respectively, and the corresponding calculation results are shown in Table 3 and Figure 4.

Initial air pressure/m | Air comprehensive coefficient | Initial water level/m | Highest water level/m | MAP/m |
---|---|---|---|---|

10.33 | 142.78 | 1.278 | 1.361 | 107.40 |

12.33 | 170.24 | 1.357 | 103.40 | |

14.33 | 197.69 | 1.354 | 100.89 | |

16.33 | 225.14 | 1.352 | 99.22 |

Initial air pressure/m | Air comprehensive coefficient | Initial water level/m | Highest water level/m | MAP/m |
---|---|---|---|---|

10.33 | 142.78 | 1.278 | 1.361 | 107.40 |

12.33 | 170.24 | 1.357 | 103.40 | |

14.33 | 197.69 | 1.354 | 100.89 | |

16.33 | 225.14 | 1.352 | 99.22 |

It can be seen from Table 3 and Figure 4 that with the increase of the initial air pressure, the air comprehensive coefficient also increased, and the maximum water level and the MAP are decreasing as the initial water level remains the same. The results are consistent with the theoretical analysis, which shows that the maximum pressure of the air pocket decreases as the initial air pressure increases when the other parameters are kept constant.

As discussed in the previous section, the MAP increases with the increase of initial air pocket length. Experimental results conducted by Martins *et al.* (2017) proved this conclusion. In the experimental study of Martins *et al.* (2017), the initial air pocket length was chosen to be 1.9 times and 9.3 times the diameter of the vertical pipe, respectively, and other parameters were kept constant. *L _{air0}* is the initial air pocket length and

*D*is the diameter of the vertical pipe. The experimental results are shown in Figure 5.

It can be seen from Figure 5 that the larger initial air pocket length results in the larger MAP. When the initial air pocket length was chosen to be 1.9 times the diameter of the vertical pipe, the MAP appeared at 0.15 s, and the value of MAP was 109.99 m. When the initial air length was chosen to be 9.3 times the diameter of the vertical pipe, the MAP appeared at 0.22 s, and the value of MAP was 130.59 m. The MAP increases with the increasing initial air pocket length, which verified the theoretical analysis.

The influence of the different air pocket volumes on the MAP has been verified by the experimental model, by which the initial air length influenced the MAP indirectly. However, the volume of the air pocket is constant in some real cases, and the effect of the shape parameters on the MAP is not understood. As the MAP increases with the increasing initial air length, it is possible to reduce the MAP by reducing the length of the initial air pocket. The diameter of the air pocket will increase as the length of the air pocket decreases when the initial volume of the air pocket is constant. Therefore, the purpose of this section is to study the effect of different diameters of air pocket on the MAP with the constant initial volume and initial pressure of entrapped air pocket. The polytropic exponent is taken as 1.4, the initial air pressure is taken as 10.33 m, the diameter of the pipes are taken as 0.038, 0.054, and 0.121 m and the initial air length are taken as 1.00, 0.50, and 0.10 m, respectively. The initial parameters and corresponding results are shown in Figure 6.

As shown in Figure 6, when the initial volume and initial pressure of the air pocket are constant, the initial air length increases as the diameter of the air pocket decreases, resulting in the increase of the MAP. When the initial air length is 0.10 m, the MAP is 107.40 m and the maximum water level variation is 0.083 m. When the initial air length is 1.00 m, the MAP is 108.87 m and the maximum water level variation is 0.814 m. As the numerical model parameters in this case are relatively small, the MAP difference is not significant, while the change law of the MAP with the increase of the initial air length is confirmed. The experimental data and numerical simulation results verified the theoretical derivation. In addition, the numerical results illustrate that the MAP decreases with the increase of the diameter of the pipe containing the entrapped air when the initial volume of the air pocket is constant.

## CONCLUSIONS

The entrapped air pockets in a pressurized water supply pipe may produce a significant transient pressure, resulting in pipe deformation or even pipe rupture. In this paper, the mathematical model of the air vessel-pressured pipe–air pocket system based on the rigid water column theory and the ideal gas equation is established, and the theoretical formula reflecting the relationship between the MAP and air comprehensive coefficient is presented. In addition, the numerical model of the system is established based on the method of characteristics. The rapid filling process of a pipeline system is simulated based on the experimental setup, and the results have been compared with the experimental data collected by Martins *et al.* (2017). A good agreement compared with the experimental data proved that the numerical model of the system has good precision. The influence of air comprehensive coefficient on the maximum pressure of the air pocket is studied, and both the theoretical formula and numerical simulation showed that the air comprehensive coefficient has a great influence on the maximum pressure of air pocket. This study found that the maximum pressure of entrapped air pocket and the fluctuation of water level in the air chamber will decrease with the increase of the air comprehensive coefficient during the rapid filling process. Air comprehensive coefficient reflects the initial characteristics of the entrapped air pocket, whose value is related to polytropic exponent, initial air pressure, and initial air length. With the increase of polytropic exponent and initial air pressure, the MAP is decreased; when the initial volume and initial pressure of the air pocket are constant, when the initial air length increased, the diameter of the air pocket decreased, resulting in the increase of the MAP.

Currently, some of the dynamic characteristics of the entrapped air pocket are understood theoretically, but more in-depth research needs to be carried out to translate these laws into effective mitigation methods.

## ACKNOWLEDGEMENTS

This paper was supported by the National Key Research and Development Program of China (Grant No. 2016YFC0401810) and the Fundamental Research Funds for the Central Universities (Grant Nos 2016B10814, 2016B41714 and 2016B04914).

## REFERENCES

*.*

*Pressurization Dynamics of Hydraulic Systems with Entrapped Air*

*PhD thesis*

*.*