## Abstract

This paper proposes a linear model that relates the pressure head variations at the downstream end of an inverted siphon to the flow rate variations at two ends. It divides the pressure head variations in the inverted siphon into low-frequency part and high-frequency part. The two parts are caused by the deformation of the siphon wall and the reflection of acoustic wave, respectively. In order to build a simplified relation between wall deformation and low-frequency pressure head variations, the Preissmann slot method (PSM) is adopted in this paper. The linear model can also be used in other forms of structures, such as pipes and tunnels, where a pressurized flow condition is present. In comparison with simulation results using the finite volume method, the linear model shows an *L*_{2} norm of 0.177 for a large-scale inverted siphon and 0.044 for a PVC pipe. To this end, the linear model is adopted to model a large-scale inverted siphon in a virtual water delivery system. Simulation results show that the inverted siphon can reduce water fluctuations. An equation to quantify this effect is proposed based on the linear model.

## NOTATION

*A*_{f}wetted area of full cross-section in m

^{2};*a*speed of acoustic wave in m/s;

*B*_{sl}width of Pressimann slot in m;

*h*relative water elevation in m;

*h*_{f}frictional head loss of quasi-steady flow in m;

*h*_{i}inertia head in m;

*Imag*resonance frequency in rad/s;

*q*relative discharge in m

^{3}/s;*s*Laplace variable in s

^{−1};*y*water depth in m;

*τ*delay time of propagation downstream in s.

## INTRODUCTION

Under the present changing climate, the conflict between water demand and water supply has become more and more critical. In order to balance the water resources in western and eastern areas of China, some long distance water distribution systems are being built. However, a large amount of water is wasted due to the low management level in these systems. Canal automatic control is an efficient tool to maintain stable flows, improve water-use efficiency and fulfill users' demands in a timely manner. It is achieved by automatic operations of control structures based on periodical observations of canal systems (Malaterre *et al.* 1998; Clemmens *et al.* 2005; Clemmens & Strand 2010). Commonly, a water level monitor is fixed at a certain point in an open-channel; this point is referred to as the target point. The water level at the target point is sent to a feedback controller periodically and a control algorithm is adopted to decide the opening of check gates so that the water surface at the target point can keep a designed level. For example, when the water surface at the target point is lower than the designed level, the check gate upstream is raised and more water is let into the open-channel so that the water surface goes up, while the feedback controller determines the movement of check gates according to the water level at the target point. A sketch to illustrate this system is shown in Figure 1.

After the check gate is raised, there is a delay time before the water level at the target point starts to rise due to the wave movement in the open-channel. Therefore, the feedback controller needs to predict the water level variations at the target point corresponding to the movement of the check gate (Clemmens *et al.* 2015; Wahlin & Zimbelman 2018). Water in open-channels can be simulated using the Saint-Venant (S-V) equations (Cunge *et al.* 1980; Rogers *et al.* 1993). However, when it comes to the design of feedback controller, the S-V equations are not easily implementable due to the large number of parameters involved. Based on the understanding of wave movements in open-channels, researchers have proposed a linear model as a simplification of the S-V equations (Schuurmans *et al.* 1999). Schuurmans's model contains a delay section that represents the delay time of wave movement in the open-channel and an integer section that reflects the effect of water volume in the open-channel on the water level at the target point; therefore, this model is referred to as the ID (integer-delay) model. The parameters of the ID model can be derived using the size of open-channels or field data, the latter of which makes it more robust (Foo *et al.* 2012, 2014). Furthermore, the ID model can be used in optimization of the feedback controller due to its simplicity (Overloop *et al.* 2005). Litrico & Fromion (2004) improved the accuracy of the ID model by adding two new parameters, which became the integer-delay-zero (IDZ) model.

In previous studies, canal automatic control is limited to open-channels because only linear models for open-channels have been proposed (Guan *et al.* 2011; Shahdany *et al.* 2016). Pressure distribution in an inverted siphon or pipe is influenced by the boundary conditions, the flow rate distribution, and the material of the siphon (or pipe) wall (Sadeghi *et al.* 2012; Okhravi *et al.* 2017). The wave movements in an inverted siphon are very different from that in free-surface flows. The speed of acoustic wave can reach 1,000 m/s which is two orders of magnitude higher than the wave speed in open-channels, indicating that the delay time between check gate movement and water level variations at the target point is much smaller (Brito *et al.* 2017; Li *et al.* 2018). Moreover, unlike in the open-channel where wave dissipates quickly due to wall friction and gravity, the wave in an inverted siphon may reflect for many times before it is completely damped. In addition, the influence of the wall deformation on pressure head variations must be taken into consideration (Zhu *et al.* 2018).

In order to apply canal automatic control to complex water distribution systems consisting of a large-scale inverted siphon, a new linear model that is especially tailored to describe the wave movements in an inverted siphon is required for the design of feedback controller. Due to the high speed of acoustic wave, the pressure head variations are very frequent and hard to analyze in the time domain. Frequency analysis has been shown as an efficient tool for evaluating the relation between pressure head variations and flow rate variations. It has been widely used in leakage detection in pressurized pipelines (Duan *et al.* 2011; Wang & Ghidaoui 2018; Lin *et al.* 2019). In this paper, frequency analysis is used to evaluate pressure head variations at the downstream end of an inverted siphon using the flow rate at two ends. The Preissmann slot method (PSM) is used in this paper to build a simplified relation between wall deformation and pressure head variations for the first time. The linear model proposed in this paper can be used for other forms of structures, like pipes and tunnels, where a pressurized flow condition is present. It can be used to construct a coupled linear model of open-channels and inverted siphons, which can then optimize controller design and improve the management level of complex water distribution systems.

This paper consists of five sections. The first section introduces the ID model and the IDZ model, as well as the PSM. The following section investigates characteristics of pressurized flows in inverted siphon using frequency analysis. The linear model is proposed and validated in the next section, then a section describes the linear model being used to model an inverted siphon in a virtual water distribution system. Conclusions are drawn in the final section.

### Preissmann slot method

The PSM has been widely used in modelling transitions between free-surface flow and pressurized flow (Dazzi *et al.* 2016; An *et al.* 2018; Maranzoni & Mignosa 2018). It permits simulation of two flow regimes using one set of governing equations, i.e., the S-V equations, by assuming a narrow slot on the top of pressurized flow and forming a virtual free-surface flow, shown in Figure 2.

*B*is computed as: where

_{sl}*A*is the wetted area of the full cross-section;

_{f}*a*is the speed of acoustic wave which can be computed according to the Korteweg-Joukowsky formula: where

*ρ*is the density of the liquid;

*K*is the bulk modulus of elasticity of the liquid;

*E*is the modulus of elasticity of wall material;

*b*is the thickness of wall. With the PSM, the pressurized flow can be simulated using the modified S-V equations: where

*H*is pressure head;

*v*is flow velocity;

*A*is wetted area;

*g*is acceleration of gravity;

*S*=

_{f}*n*

^{2}

*v*

^{2}/

*R*

^{4/3}is hydraulic slope;

*n*is Manning coefficient;

*R*is hydraulic radius;

*S*is bed slope.

_{b}*A*and pressure head

*H*follows:

As shown by Equation (7), the width of the virtual slot *B _{sl}* can be regarded as the ratio of wall deformation to pressure head variations, it is a parameter which reflects the wall material, wall size, and wave speed. Equation (7) is a simplified relation because it ignores the time-dependent nature of stress-strain behavior of the wall material (Bertaglia

*et al.*2018). In this paper, it is used to derive the linear model of inverted siphon for its simplicity.

### ID model and IDZ model of open-channel

*Q*(

_{0}*x*) and

*Y*(

_{0}*x*), where

*Q*is flow rate and

*Y*is water depth. Assume

*Q*(

*x*,

*t*) =

*Q*(

_{0}*x*) +

*q*(

*x*,

*t*), and

*Y*(

*x*,

*t*) =

*Y*(

_{0}*x*) +

*y*(

*x*,

*t*): where

*v*

_{0},

*c*

_{0},

*B*

_{0},

*S*

_{f0},

*Y*

_{0},

*A*

_{0},

*T*

_{0}, and

*P*

_{0}(wetted perimeter) are functions of

*x*;

*c*

_{0}is the speed of gravity wave; and 0 denotes the value of each term at the initial state (Schuurmans

*et al.*1995). In the ID model, an open-channel is divided into uniform flow area and backwater area. The flow profile is parallel to the bed in uniform flow area and almost horizontal in the backwater area, as shown in Figure 3.

*q*

_{1}=

*Q*(0,

*t*) −

*Q*

_{0}(0),

*y*

_{1}=

*Y*(0,

*t*) −

*Y*

_{0}(0),

*q*

_{2}=

*Q*(

*L*,

*t*) −

*Q*

_{0}(

*L*),

*y*

_{2}=

*Y*(

*L*,

*t*) −

*Y*

_{0}(

*L*) and

*L*is the length of open-channel. Schuurmans proved that the wave traveling upstream dampens exponentially with distance whereas the wave moving downstream always reaches the downstream end in open-channels. The variations of water depth in the backwater area can be computed using the water volume variations in the backwater area: where

*A*is the surface area of backwater;

_{s}*τ*is the travelling time of wave from upstream end to downstream end (Schuurmans

*et al.*1995, 1999). Litrico & Fromion (2004) proposed to add two new parameters to improve the accuracy of the ID model: where

*y*(

*t*) is an intermediate variable;

*b*and

_{u}*b*are parameters that reflect the instantaneous variations of

_{d}*y*

_{2}when wave peak is present at the downstream end. Figure 4 is a sketch of

*h*

_{2}computed using the ID model, the IDZ model, and the realistic variations. Before

*t*=

*t*

_{0}, a steady flow condition is ensured, then a step change

*q*

_{1}at the upstream end arises. During the time the wave is travelling to the downstream end,

*y*

_{2}equals zero. As soon as it reaches the backwater area, the water depth at the downstream end rises. In spite of the wave peak, the rising speed of

*y*

_{2}equals

*q*

_{1}/

*A*, as shown in Figure 4.

_{s}The ID model and the IDZ model are not appropriate for inverted siphons due to the different flow characteristics. In this paper, the authors propose a linear model that describes the variations of pressure head in inverted siphon as a combination of a low-frequency part and high-frequency part.

## FREQUENCY DOMAIN ANALYSIS

*H*is pressure head and

*H*(

*x*,

*t*) =

*H*

_{0}(

*x*) +

*h*(

*x*,

*t*),

*h*

_{1}=

*H*(0,

*t*) −

*H*

_{0}(0) and

*h*

_{2}=

*H*(

*L*,

*t*) −

*H*

_{0}(

*L*);

*L*is length of inverted siphon. The derivation of Equations (12) and (13) is listed in the Appendix. The state-transition matrix can be obtained by applying a Laplace transform to Equations (12) and (13): where

*q*

_{1}(

*s*) and

*q*

_{2}(

*s*) are the Laplace transforms of

*q*

_{1}(

*t*) and

*q*

_{2}(

*t*),

*h*

_{1}(

*s*) and

*h*

_{2}(

*s*) are Laplace transforms of

*h*

_{1}(

*t*) and

*h*

_{2}(

*t*), respectively.

**(**

*Φ**L*,

*s*) is state-transition matrix of

*q*

_{1}(

*s*),

*h*

_{1}(

*s*) and

*q*

_{2}(

*s*),

*h*

_{2}(

*s*); for the derivation of

**(**

*Φ**L*,

*s*), readers are referred to Litrico & Fromion (2004). The transfer function between

*h*

_{2}(

*s*) and

*q*

_{1}(

*s*) is noted as

*p*

_{21}. The transfer function between

*h*

_{2}(

*s*) and

*q*

_{2}(

*s*) is noted as

*p*

_{22}. where

*λ*

_{1}(

*s*) and

*λ*

_{2}(

*s*) are eigenvalues of

**(**

*Φ**L*,

*s*).

In this paper, two scenarios are investigated. Scenario A is an inverted siphon in the Dianzhong water distribution system. It consists of five parallel closed conduits and each one has a circular cross-section, shown as Figure 5. The parameters of scenario A are listed in Table 1. The speed of acoustic wave is 1,014 m/s according to the Korteweg-Joukowsky formula. The slot width is computed using Equation (1) and is 5.9885 × 10^{−4} m.

Diameter (m) | 4 |

Number of conduits | 5 |

Length (m) | 516 |

Acoustic wave speed (m/s) | 1,014 |

Designed flow rate (m^{3}/s) | 120 |

Bed slope | −1/460 |

Manning coefficient | 0.014 |

Slot width (m) | 5.9885 × 10^{−4} |

Diameter (m) | 4 |

Number of conduits | 5 |

Length (m) | 516 |

Acoustic wave speed (m/s) | 1,014 |

Designed flow rate (m^{3}/s) | 120 |

Bed slope | −1/460 |

Manning coefficient | 0.014 |

Slot width (m) | 5.9885 × 10^{−4} |

Scenario B is an experimental horizontal polyvinyl chloride (PVC) pipe, shown as Figure 6. In order to avoid the vibration of PVC pipe, it is fixed to a concrete foundation. Its parameters are listed in Table 2. The speed of acoustic wave which is validated by experiments is 348 m/s, and the slot width computed using Equation (1) is 3.5218 × 10^{−6} m.

Diameter (m) | 0.2354 |

Length (m) | 271.3 |

Acoustic wave speed (m/s) | 348 |

Designed flow rate (m^{3}/s) | 0.007 |

Bed slope | 0 |

Manning coefficient | 0.010 |

Slot width (m) | 3.5218 × 10^{−6} |

Diameter (m) | 0.2354 |

Length (m) | 271.3 |

Acoustic wave speed (m/s) | 348 |

Designed flow rate (m^{3}/s) | 0.007 |

Bed slope | 0 |

Manning coefficient | 0.010 |

Slot width (m) | 3.5218 × 10^{−6} |

### Bode plots of transfer function

Bode plot is a common method to study transfer functions in frequency domain. Its abscissa is frequency (Freq rad/s) of input and its ordinate is magnification (Gain dB) and phase difference (Phase degree) of output. Figures 7–10 are bode plots of *p*_{21} and *p*_{22} which contain 500 frequency points that range from 0.001 rad/s to 100 rad/s.

At low-frequency domain where *Freq* < 1 rad/s, the magnification line of *p*_{21} is straight with a slope of −20 dB while the phase line is horizontal. Thses lines are the same as that of an integrator section. The phase line is below zero in low-frequency domain, which means there is a delay in *p*_{21}. Thus *p*_{21} can be estimated by an integrator section and a delay section. The magnification line of *p*_{22} is similar to that of *p*_{21} at low-frequency domain. However, the phase line is above zero indicating that there is no delay time in *p*_{22}. Thus, *p*_{22} can be estimated by an integrator section only. For *Freq**>* 1 rad/s, the magnification lines of *p*_{21} and *p*_{22} change in an oscillatory manner indicating that waves are superimposed in high-frequency domain. Several frequencies of the magnification peaks are studied in the following text.

### Poles of transfer function

*Imag*, the authors assume

*σ*and

*γ*to be zero, which is a frictionless and horizontal case. Under this assumption, Equation (18) becomes:

*Imag* of the first five poles and frequencies of the first five magnification peaks in two scenarios are listed in Table 3. Where *k* is the serial number of pole, *j*^{2} = −1. It is shown that the frequencies are close to the imaginary part of poles with a difference of less than 0.2 rad/s which is introduced due to the selected finite frequencies.

K | Imag of scenario A (j) | Frequencies of magnification peaks of scenario A (rad/s) | Imag of scenario B (j) | Frequencies of magnification peaks of scenario B (rad/s) |
---|---|---|---|---|

0 | 0 | 0 | ||

1 | 6.1736 | 6.132 | 4.0298 | 4.048 |

2 | 12.3471 | 12.25 | 8.0595 | 8.088 |

3 | 18.5207 | 18.56 | 12.0893 | 11.97 |

4 | 24.6943 | 24.48 | 16.1190 | 16.16 |

… | … | … | … | … |

K | Imag of scenario A (j) | Frequencies of magnification peaks of scenario A (rad/s) | Imag of scenario B (j) | Frequencies of magnification peaks of scenario B (rad/s) |
---|---|---|---|---|

0 | 0 | 0 | ||

1 | 6.1736 | 6.132 | 4.0298 | 4.048 |

2 | 12.3471 | 12.25 | 8.0595 | 8.088 |

3 | 18.5207 | 18.56 | 12.0893 | 11.97 |

4 | 24.6943 | 24.48 | 16.1190 | 16.16 |

… | … | … | … | … |

Note that *L*/*α* + *L*/*β* is the duration time of a wave moving in an inverted siphon back and forth once. It can be derived that for *q*_{1} and *q*_{2}, whose frequencies are multiples of wave reflection, frequency causes larger amplitude of *h*_{2} as a resonance. To validate this conclusion, a series frequencies of *q*_{1} around 6.132 rad/s are simulated for scenario A, and a series frequencies of *q*_{1} around 4.048 rad/s are simulated for scenario B. The simulated frequencies are shown in Figures 11 and 12. The amplitudes of *q*_{1} for scenario A and scenario B are 0.01 and 0.001 m^{3}/s, respectively. The results are shown in Figures 13 and 14.

For scenario A, *q*_{1} with a frequency of 6.132 rad/s (solid line) causes the largest amplitude in variations of *h*_{2} among all the tested frequencies. This conclusion can also be applied to scenario B. In order to avoid the large amplitude of pressure head variations in inverted siphons and pipes, flow rate change whose frequencies are close to the frequencies computed using Equation (20) should be avoided.

## LINEAR MODEL OF INVERTED SIPHON

*B*. Considering that wall deformation is a gradual process, the deformation at the outlet can be approximated using the average deformation along the siphon and pipe. Neglecting the compressibility of water, the average wall deformation can be quantified using the variations of water volume in the siphon: where

_{sl}*V*is the total volume of water in the siphon or pipe. In addition, the variations of water volume can be computed using the flow rate at inlet and outlet: where

*τ*is the delay time of acoustic wave that travels to downstream end from upstream end, it can be computed approximately using: Applying Equations (23) and (21) to Equation (22), the low-frequency variations of pressure head can be computed using: where

*B**

_{sl}*L*corresponds to the integrator section in

*p*

_{21}and

*p*

_{22};

*L*/

*a*corresponds to the delay section in

*p*

_{21}.

*Freq*> 1 rad/s, the magnification line of

*p*

_{21}and

*p*

_{22}shows obvious oscillations. Previous research has proven that the high-frequency variations of pressure head are caused by reflection of acoustic wave (Skulovich

*et al.*2015; Mitosek & Szymkiewicz 2016). In this paper, the authors use the water hammer equation to descrbe the high-frequency variations. Whereas a positive

*q*

_{1}increases

*h*

_{2}, a positive

*q*

_{2}decreases

*h*

_{2}; the high-frequency variations can be expressed as:

Although Equations (26) and (27) are similar to the IDZ model in form, the meaning of the parameters involved are completely different. In the IDZ model, the low-frequency variations of water depth are caused by water volume variations in the backwater area, and the high-frequency variations are caused by wave peaks. In the linear model for inverted siphon, the low-frequency variations are caused by wall deformation and the high-frequency variations are caused by reflection of acoustic wave.

### Comparison of linear model with p_{21} and p_{22}

Figures 15–18 are comparisons of *p*_{21} and *p*_{22} with the linear model in frequency domain. Results show that the linear model is accurate in the low-frequency domain but less accurate in the high-frequency domain since it cannot present the oscillations. However, pressure head variations of very high frequencies dampen very quickly and rarely appear. Thus they are less important in the effectiveness of the linear model for practical use.

### Comparison of linear model with simulation and experimental results

*et al.*2008): where the flow variable vectors

**U**, the flux vectors

**F**, and the source term vectors

**S**are: where

*p*is the pressure that acts on the center of the cross-section. In this paper, a quasi-steady friction model is adopted (Ghidaoui Mohamed

*et al.*2002). A first-order Godunov-type scheme with a Roe solver is used for numerical simulation. The boundary condition employs the ghost cells proposed by Toro (2002).

^{3}/s, and a water tank is placed at the upstream end to provide a constant water head of 21.4 m. The ghost cell at the upstream boundary applies a constant pressure head of 21.4 m while the discharge adopts a transmissive condition. The ghost cell at the downstream boundary applies a constant discharge of 120 m

^{3}/s while the pressure head adopts a transmissive condition. The distance step of simulation is 25.8 m, and the time step is 0.0242 s. The linear model of scenario A is written in incremental form as: where Δ

*t*is the time step. Figure 19 shows the pressure at the downstream end (

*P*) in the simulation results (

_{bv}*FVM*) and linear model (

*Res*). The values of

*q*

_{1}and

*q*

_{2}used in Equation (30) are obtained from simulation results of finite volume method.

*et al.*(2011). In the experiment, a water tank is placed at the inlet to provide a constant water head of 21.4 m, the initial flow rate is 0.007 m

^{3}/s, and a ball valve is set at the outlet to control the outflow rate. At the beginning of the experiment, the ball valve is closed immediately. In simulation, the distance step is 13.5 m, the time step is 0.0369 s, and the Courant number is 0.95. The ghost cell at the upstream boundary applies a constant pressure head of 21.4 m while the discharge adopts a transmissive condition. The ghost cell at the downstream boundary applies a constant discharge of 0 m

^{3}/s while the pressure head adopts a transmissive condition. The linear model of scenario B can be written in incremental form as:

Figure 20 shows the pressure at ball valve (*P _{bv}*),

*FVM*represents

*P*obtained using finite volume method,

_{bv}*Exp*represents

*P*obtained by experiments,

_{bv}*Res*represents

*P*obtained by Equation (31). The values of

_{bv}*q*

_{1}and

*q*

_{2}used in Equation (31) are obtained from simulation results of finite volume method.

According to *Exp*, when the valve is closed a wave is generated at the ball valve. The *P _{bv}* changes in a periodical manner due to the reflection of acoustic wave between the inlet and the valve. At the first 5 seconds, some variations which are of higher frequency than wave reflection exist in the experimental data and then quickly disappear. From

*t*= 5 s to

*t*= 20 s, two or three frequencies of variations are present in

*P*while from

_{bv}*t*= 20 s to 40 s only one frequency of variations remains. In this process, the amplitude of pressure gradually decreases due to the stress-strain behavior.

*L*,

*A*, and

_{f}*a*. Among them, only

*a*is an implicit value which needs to be calibrated by experiments. In order to evaluate the influence of wave speed on the accuracy of the linear model, different values of wave speed are adopted in Equations (24) and (25). Herein, the

*L*

_{2}norm is defined as an indicator of accuracy in comparison with FVM: where

*p*is pressure computed using the linear model,

_{Res}*p*is pressure in simulation results of finite volume method, and Δ

_{FVM}*t*is the time step. The

*L*

_{2}norms of different wave speed are shown in Figures 21 and 22.

The abscissa is the wave speed in multiples of the correct wave speed *a*. It can be seen that the minimum *L*_{2} norm is achieved when *a* is used; the minimum *L*_{2} norms are 0.177 and 0.044 for scenario A and B, respectively. The linear model is sensitive to the wave speed, when a false speed is used, the *L*_{2} norm can reach to 6 and 7 for scenario A and B, respectively. Thus, a correct value of *a* is critical in the accuracy of the linear model. The value of *a* may be affected by the upstream structures in practice (Mitosek & Szymkiewicz 2016), so it needs to be calibrated by experiments.

## INFLUENCE ON WATER FLUCTUATION OF INVERTED SIPHON

Finally, the linear model is applied to the model inverted siphon in a conceptualized water distribution system, shown as Figure 23. It consists of two open-channels and one inverted siphon. Open-channels 1 and 2 are prismatic with a trapezoidal cross-section while the inverted siphon has a circular cross-section. The wave speed in the inverted siphon is set at 1,014 m/s, and the slot width is 2.9943 × 10^{−5} m according to Equation (1). A gate is placed at the upstream end of open-channel 1 to control the inflow rate of the entire system. A weir is placed at the downstream end of open-channel 2. Model parameters are listed in Table 4. This model is simulated using the finite-difference method with a time step of 0.2 s and a distance step of 10 m for two open-channels and the inverted siphon.

Length (m) | Slope | Manning coefficient | |
---|---|---|---|

Open-channel 1 | 1,000 | 5/10,000 | 0.015 |

Open-channel 2 | 1,000 | 5/10,000 | 0.015 |

Inverted siphon | 2,000 | 1/2,000 | 0.014 |

Length (m) | Slope | Manning coefficient | |
---|---|---|---|

Open-channel 1 | 1,000 | 5/10,000 | 0.015 |

Open-channel 2 | 1,000 | 5/10,000 | 0.015 |

Inverted siphon | 2,000 | 1/2,000 | 0.014 |

The initial flow rate is 3 m^{3}/s, then at *t* = 5 min, the flow rate at the gate increases to 8 m^{3}/s and remains constant. Figure 24 shows that, at *t* = 9 min, *h*_{1} and *h*_{2} of the inverted siphon start to change. The variations of *h*_{1} show obvious oscillations while the variations of *h*_{2} are smaller and mild. This finding indicates that the inverted siphon can filter the high-frequency variations and decrease the amplitude of variations.

*h*

_{2}consist of low- and high-frequency parts. At the low-frequency domain, the deformation of siphon wall at the two ends are similar, the difference between

*h*

_{1}and

*h*

_{2}is frictional head loss

*h*. Because of the transmissive condition in the inverted siphon, only the wave that moves from the upstream end to the downstream end is present. Thus, the high-frequency variations of

_{f}*h*

_{2}can be computed using the first term at the right side of Equation (25):

*h*

_{1}can be determined similar to the second term at the right side of Equation (25):

*τ*is very short, the following transformation can be made: where

*τ*can be approximated by

*L*/

*a*. Then, Equation (36) becomes:

In this paper, the second term at the right side of Equation (37) is noted as *h _{i}*. The frictional head loss

*h*is computed using the flow rate at the upstream end of the inverted siphon under the assumption of a quasi-steady flow condition. The comparisons between

_{f}*h*

_{1}–

*h*,

_{f}*h*

_{1}–

*h*–

_{f}*h*and

_{i}*h*

_{2}are shown in Figure 25. The

*L*norm between

_{2}*h*

_{1}–

*h*and

_{f}*h*

_{2}is 1.206 while the

*L*

_{2}norm between

*h*

_{1}–

*h*–

_{f}*h*and

_{i}*h*

_{2}is 0.026. This result shows that an inverted siphon reduces the fluctuation of water level in open canals. The effect can be quantified using Equation (37).

## CONCLUSIONS

A linear model relating the pressure head variations at the downstream end to the flow rate variations at two ends of an inverted siphon or pipe is proposed in this paper. The linear model divides the pressure head variations into low-frequency part and high-frequency part, which are caused by wall deformation and wave reflection, respectively. In the high-frequency domain, specific frequencies of flow rate variations at the two ends may cause a larger amplitude in pressure head variations as a resonance. An equation is proposed to compute these frequencies, i.e., Equation (20). For the low-frequency variations, the Preissmann slot method is adopted to derive a simplified relation between wall deformation and pressure head variations. While for the high-frequency variations, the water hammer equation is adopted to evaluate the pressure head variations caused by acoustic wave reflection. The parameters of the linear model are computed using *L*, *a*, and *A _{f}*. Two scenarios are studied in this paper, one is a large-scale inverted siphon and the other one is a PVC pipe. The accuracy of the linear model is verified in frequency domain using Bode plot in two scenarios. Then the pressure head computed using the linear model is compared with the simulation results of FVM in two scenarios. The accuracy of this model is based on an accurate speed of acoustic wave. The minimum

*L*

_{2}norms are achieved when a correct wave speed is used. The minimum

*L*

_{2}norms are 0.177 and 0.044 for the inverted siphon and the PVC pipe, respectively. When a false speed is used, the

*L*

_{2}norms can reach 6 and 7 for the inverted siphon and the PVC pipe, respectively. Thus, a correct speed of acoustic wave is important to guarantee the accuracy of the linear model. The linear model is subsequently applied to describe an inverted siphon in a conceptualized water distribution system. Simulation results show that the inverted siphon can reduce water fluctuations. This effect is quantified by a frictional head loss

*h*under the assumption of a quasi-steady flow and an inertia head multipled by

_{f}*L*/(

*A*). On the basis of the aforementioned analysis, the following conclusions can be drawn:

_{f}gAn equation to compute the frequencies of flow rate variations at the two ends that may cause resonant pressure head variations in an inverted siphon or pipe is proposed. These frequencies are multiples of the wave reflection frequency between the two ends. In order to avoid damage of the siphon or pipe wall caused by pressure variations, flow rate variations around these frequencies should be avoided.

The pressure head variations in an inverted siphon or pipe can be divided into a low-frequency part and high-frequency part. The low-frequency part is caused by wall deformation while the high-frequency part is caused by the reflection of acoustic wave.

A linear model is proposed to calculate the pressure head variations at the downstream end of the inverted siphon using the flow rate variations at two ends. This model is applicable to other structures, such as pipes and tunnels, where pressurized flows are present. The low-frequency pressure head variations are related to the flow rate at two ends using a virtual slot width

*B*, which represents the ratio of wall deformation to pressure head variations. The high-frequency variations are computed using the water hammer equation. The linear model is sensitive to the speed of acoustic wave, but in practice, it needs to be calibrated by experiments to ensure the accuracy of the linear model._{sl}An equation to quantify the reduction effect of an inverted siphon on the water fluctuations is proposed. The reduction effect can be computed using a frictional head loss under the assumption of a quasi-steady flow condition and an inertia head multiplied by

*L*/(*A*)._{f}g

The linear model proposed in this study can be used to build a coupled model of open-channel and pressured flows (i.e., syphon and tunnel), which can then optimize the feedback controller. These findings may be a reference for large water delivery systems.

## ACKNOWLEDGEMENTS

The authors acknowledge the support of the NSFC grant 51979202.