Quasi-static Flow Model for Predicting the Extreme Values of Air Pocket Pressure in Draining and Filling Operations in Single Water Installations

: Inertial models have been used by researchers to simulate the draining and filling processes in water pipelines, based on the evolution of the main hydraulic and thermodynamic variables. These models use complex differential equations, which are solved using advanced numerical codes. In this study, a quasi-static flow model is developed to study these operations in hydraulic installations. The quasi-static flow model represents a simplified formulation compared with inertial flow models, in which its numerical resolution is easier because only algebraic equations must be addressed. Experimental measurements of air pocket pressure patterns were conducted in a 4.36 m long single pipeline with an internal diameter of 42 mm. Comparisons between measured and computed air pocket pressure oscillations indicate how the quasi-static flow model can predict extreme values of air pocket pressure for experimental runs, demonstrating the possibility of selecting stiffness and pipe classes in actual pipelines using this model. Two case studies were analysed to determine the behaviour of the quasi-static flow model in large water pipelines.


Introduction
The interaction of water columns with entrapped air pockets in water pipelines is difficult to understand because hydraulic and thermodynamic formulations need to be analysed [1,2]. Furthermore, the understanding of water and the air phase is crucial to determine the evolution of involved variables [3].
The water phase can be analysed using flow models depending on fluid and pipe characteristics, where two groups of flow models have been used: (i) inertial models, which consider the fluid inertia, and its resolution is conducted using either the method of characteristics for solving the elastic model (water hammer) [4,5] or numerical methods when the analysis is addressed with a rigid column model (also named the mass oscillation equation) [6,7]; (ii) non-inertial models (or quasi-static flow models), where the fluid inertia is neglected, and its numerical resolution can be computed using a system composed of algebraic equations [1].
• The polytropic model describes the thermodynamic behaviour of the air phase [5,10].
• The air-water interface is simulated using a piston flow model [2,16]. Sections 2.1 and 2.2 present the mathematical formulations of the draining and filling processes using the quasi-static flow model, respectively.

Draining Process
The draining process in water pipelines produces sub-atmospheric pressure conditions because the initial air pocket size is expanded during the hydraulic event. Figure 1 illustrates the draining process in a single hydraulic installation, where the main hydraulic and thermodynamic variables are shown. The draining process begins when a drain valve (installed on the downstream end) is opened; subsequently, the water column exits by a drain valve, and an entrapped air pocket increases its initial size. The draining process cannot be concluded based on the opening percentage of the drain valve. This shows that part of the water column remains inside the considered pipeline, and the hydraulic installation should be under a sub-atmospheric condition. The draining process can be described as follows: • The steady-state equation. The mass oscillation equation describes the water column movement during a draining operation, and its formulation for a single pipeline is as follows: where v = water velocity, t = time, * = air pocket pressure, * = atmospheric pressure, = water density, = water column size, = pipe elevation at point i, = friction factor, D = internal pipe diameter, A = cross-sectional area of pipe and Rv = resistance coefficient of a drain valve. The terms used in the mass oscillation equation are as follows: ⁄ ( / ) is the variation of water velocity over time; ( * − * )/ is the difference in pressure energies between points 1 and 2; ( − ) is the difference in potential energies between points 1 and 2; | |/ (2 ) represents the friction losses from points 1 to 2; | | represents the minor losses in the valve located at point 2.
Neglecting the inertial term (dv/dt = 0) in Equation (1), the steady-state equation (or Bernoulli's equation) can be expressed mathematically as follows [1]: Equation (2) is used to compute the evolution of water velocity in pipe systems and is categorised as a non-inertial flow. This equation is applicable for slow transient flows.

•
The piston flow formulation. The proposed model considers a perpendicular air-water interface to the main direction of the flow; consequently, it can be applied for water pipelines with a significant slope in their axis. Although real pipelines present a horizontal air-water interface, the considered hypothesis has been implemented by many researchers to obtain a suitable approximation of this phenomenon [1,4,8]. The air-water interface can be modelled as follows [5,16]: Applying the finite backward difference approximation, Equation (3) becomes where = air volume, k = polytropic coefficient and = time at 0 s (the hydraulic system is at rest). When the polytrophic coefficient is 1.0, the air pocket exhibits an isothermal thermodynamic behaviour, and the formulation is named as Boyle´s law; if k = 1.4, the air pocket presents an adiabatic evolution [18]. In real water pipelines, intermediate values of the polytropic coefficient can occur between 1.0 and 1.4.
Considering a cross-sectional area of a pipe and an air pocket size ( = − , ), an explicit formulation can be obtained with the involved variables [5]: where = total length of a pipeline.

Filling Process
The filling operation is characterised by the compression of an entrapped air pocket, which generates pressure surges. Figure 2 shows all involved variables to describe this process. The operation starts when a regulating valve is opened; subsequently, the pipeline installation begins to be filled by an energy source (pump or tank). As the dead end of a pipeline is closed, an entrapped air pocket is compressed, producing peaks of air pocket pressure. The filling operation can be described using the following formulations: • The steady-state equation. The Bernoulli's equation is applied to the water column from points 1 to 2 (see Figure 2), which should be considered for emptying operations: where * = absolute pressure supplied by a pump or tank, and = resistance coefficient of a regulating valve.
• The piston flow formulation. The air-water interface can be described as follows: The polytropic equation of an air pocket. Equation (6) can be used to represent the behaviour of an entrapped air pocket during a filling operation.

Numerical Resolution
The resolution of draining processes in single pipelines using the quasi-static model involves the determination of three unknown variables: , , and , * for each time step (∆ ). The resolution of the system implies an extended period simulation of single water pipelines. Figure 3 shows the flowchart to solve the set of three algebraic formulations presented in Sections 2.1 and 2.2 for draining and filling processes, respectively. The first stage is the reading of water parameters, topographical characteristics of a single pipeline, location of the initial position of the air-water interface and valve coefficients. The pipe slopes were obtained using the formulation, sin = ∆ ( − ) ⁄ . Next, when the hydraulic system is at rest ( = 0), the water velocity is null, and the air pocket pressure is at the atmospheric condition. For the remaining time steps, the three unknown variables are computed by solving Equations (2), (4) and (6) for each time step for a draining process; meanwhile, a filling process can be computed using Equations (6), (7) and (8). At the end of transient phenomena, no changes occur between and in the estimation of water velocity, air pocket pressure or length of the draining column as the selection of the steady-state model omits the inertial term ( ⁄ = 0).  (2), (4), and (6)) Filling operation (Equations (6), (7), and (8)) , , , , * +1 = ∆ + End of the extended period simulation End Yes No Figure 4). A manual ball valve (MBV1) was positioned at the upstream end to fill the pipe installation. A pressure transducer was located at the highest point of the hydraulic installation to record air pocket pressure oscillations. In addition, another manual ball valve (MBV2) was placed at the downstream end, where the draining operation was performed. During all experiments, air pocket sizes of 0.205 m (0.67 ft), 0.340 m (1.11 ft) and 0.450 m (1.48 ft) were configured in the pipeline, as well as pipe slopes of 0.457 rad and 0.515 rad. Different valve travels in the MBV2 were performed, which were defined through marks on a parallel white surface with regard to the position of the MBV2. Table 1 shows the characteristics of all performed runs, which have been previously discussed by the authors [9].

Model Validation
This section presents a comparison between computed and measured air pocket pressure patterns for 12 runs, as shown in Figure 5. The analysed pipeline installation was composed of an inclined and a vertical pipe branch measuring 4.16 m (13.65 ft) and 0.2 m (0.66 ft), respectively. For all runs, the water column remained inside the inclined pipe branch. The gravity term ( − ) , ⁄ (see Equation (2)) can be calculated as follows: During the 12 runs, a constant friction factor (f) of 0.018 and a permanent polytropic coefficient (k) of 1.4 (adiabatic condition) were obtained based on experimental results, as mentioned in a previous publication by the authors [9], as the hydraulic event occurred rapidly. The quasi-static model is suitable for representing Run 6 and Runs 10 to 12, as presented in Figure 5, where the trends of the experimental results exhibit no oscillations for the highest values of resistance coefficient in MBV2; however, it is not appropriate for simulating air pocket pressure oscillations from Runs 1 to 5 and Runs 7 to 9, where low values of resistance coefficient of MBV2 exist. The quasi-static model cannot simulate these absolute pressure oscillations as Equation (2) neglects the inertial term (dv/dt = 0).
Typically, engineers size a drain valve orifice between 1/3 and 1/4 of an internal pipe diameter in actual pipelines [12]. Therefore, the quasi-static flow model can be used to study the behaviour of the main hydraulic and thermodynamic variables during a draining operation because low values of flow factor in drain valves are obtained in real hydraulic installations.
For all runs, the quasi-static model can predict the drop of sub-atmospheric pressure patterns, which is critical because the minimum sub-atmospheric pressure is required for selecting a stiffness class pipe. For instance, for Run No. 1, the minimum value recorded by the pressure transducer of an absolute pressure head was 8.03 m, and using the quasi-static flow model, a value of 8.42 m was obtained; by contrast, for Run No. 12, a sub-atmospheric pressure head of 8.30 m was obtained, which was close to that of the quasi-static model (8.36 m).
A comparison between the inertial model developed by the authors [6,9] and the quasi-static model is shown in Figure 6 for Runs 1 and 11, which are two representative experimental tests. The remaining tests exhibit similar air pocket pressure patterns. The comparison indicates that the inertial model can accurately predict air pocket pressure patterns for all runs as it considers the inertial term ( ⁄ ≠ 0). In Run No. 1, the inertial model predicts not only the drop of the sub-atmospheric pressure head (8.03 m), but also the absolute pressure pattern. Run No. 11 was predicted using both the inertial and quasi-static flow models. The mean square error (MSE) of the minimum air pocket pressure head for the 12 runs was analysed. Using the quasi-static model, the MSE was 2.494%, implying a good estimation of extreme conditions. As expected, the inertial model agrees well with the calculations of air pocket pressure head oscillations (0.402 %) compared with the quasi-static model. Although the inertial model could predict the behaviour of sub-atmospheric pressure patterns more accurately, both models provided good results in terms of extreme values, which enables a stiffness pipe class to be selected using both models.

Case studies
This section presents two case studies regarding draining and filling operations in water distribution systems, which are also used to validate the quasi-static flow model.

Draining process
The simulation of a large draining process in a single pipe was performed using the following parameters: LT = 600 m, D = 0.35 m, k = 1.2, f = 0.018, ∆z/Lw,t = sin θ = 0.025 m/m, Rv = 0.06 m/s/m 6 and t = 10 s. The initial conditions were described by the following parameters: x0 = 200 m, v(0) = 0 m/s and * = * = 101325 Pa. Figure 7 presents the results of water column length, water velocity and air pocket pressure. The behaviours of the draining process using the quasi-static flow model are summarised as follows: • The quasi-static model can predict the evolution of the water column length during the draining process, as shown in Figure 7a.    Figure 9 shows a comparison between the quasi-static and inertial models [7] of the main hydraulic and thermodynamic variables. The following can be deduced regarding filling based on the results shown in Figure 9: • The quasi-static model is not suitable for predicting oscillation patterns of water column position, water velocity, and air pocket pressure because it neglects the internal term (dv/dt = 0); however, extreme values can be predicted by the quasi-static model, which are used to select the pipe class. • The quasi-static flow and inertial models can predict the final position of the air-water interface; as shown in Figure 9a, the water column approaches 382 m (1253.3 ft) from 120 s to the end of the hydraulic event.

Conclusions
In this study, a new one-dimensional mathematical model was developed to compute the main hydraulic and thermodynamic variables during draining and filling processes with an entrapped air pocket in water installations, which was based on the quasi-static flow model, the polytropic formulation of an air pocket and the piston flow model to represent the air-water interface.
The quasi-static model represents an alternative formulation, which can be used by engineers to model filling and draining operations in water single installations, as it can adequately predict extreme values of air pocket pressure patterns. The quasi-static model can be programmed using typical applications (e.g., Microsoft Excel) because complex equations such as the elastic water model or mass oscillation equation are not necessitated, as they are replaced by Bernoulli's equation. Consequently, the numerical resolution of the proposed model is much easier to use compared with inertial models. The mathematical model was validated in an experimental facility housing a 4.36-m-long (14.3ft) methacrylate pipeline with an internal pipe diameter of 42 mm (1.65 inches). Results indicated that the proposed model was suitable for predicting the drop in air pocket pressure patterns in a draining process, which was critical as this extreme value was used to select pipeline characteristics. The minimum drop in sub-atmospheric pressure of the proposed model was compared with the attained values of the inertial model (rigid water column model), in which MSE values of 2.494% and 0.402% were obtained, respectively. The validation of a filling process was not conducted in the experimental facility as these formulations could be applied similarly as a draining process.
Finally, the behaviour of the inertial and proposed models was compared using two cases studies of large diameter (0.35 m or 1.15 ft) for filling and draining operations. The quasi-static model yielded similar results of extreme values of air pocket pressure oscillations compared with the inertial model. The time step did not significantly affect the extreme values of air pocket pressure in the proposed model. Therefore, in a design stage, engineers can apply the quasi-static model to select a stiffness pipe and a class pipe for draining and filling processes, respectively, in water single installations.
The next step is to apply the quasi-static model for complex water pipeline systems.