# Shock Wave Speed and Transient Response of PE Pipe with Steel-Mesh Reinforcement

Authors: Wuyi Wan; Xinwei Mao
Journal: Shock and Vibration
 (2016)

---

## Abstract

A steel mesh can improve the tensile strength and stability of a polyethylene (PE) pipe in a water supply pipeline system. However, it can also cause more severe water hammer hazard due to increasing wave speed. In order to analyze the influence of the steel mesh on the shock wave speed and transient response processes, an improved wave speed formula is proposed by incorporating the equivalent elastic modulus. A field measurement validates the wave speed formula. Moreover, the transient wave propagation and extreme pressures are simulated and compared by the method of characteristics (MOC) for reinforced PE pipes with various steel-mesh densities. Results show that a steel mesh can significantly increase the shock wave speed in a PE pipe and thus can cause severe peak pressure and hydraulic surges in a water supply pipeline system. The proposed wave speed formula can more reasonably evaluate the wave speed and improve the transient simulation of steel-mesh-reinforced PE pipes.


---

## Body

## 1. Introduction

Water hammer events can cause noise, vibration, distortion, and fracture in a water supply pipeline system due to sudden increase or decrease of pressure [1, 2]. Numerical simulation and prediction are very important to protect the pipeline from undesired damage due to the water hammer effect [3]. It is the fundamental to control the peak transient pressure by optimal operation [4–6], as well as select and design protection devices in water supply pipelines [7]. Furthermore, it is also applicable to some other pipeline transportation areas, such as oil pipelines [8, 9] and natural gas pipelines [10, 11]. Shock wave speed is a primary factor in various transient simulation models [12–16]. It represents the basic properties of the fluid and pipeline system in transient flow simulations [17]. Usually, wave speed can greatly change the frequency and amplitude of the water hammer waves, as well as the extreme transient pressure distributions along the pipeline [18]. Therefore, the reliability of the result depends heavily on the wave speed in the numerical simulation of the water hammer. Generally, the wave speed is subjected to many factors, such as the density and elastic modulus of the fluid, the material [19, 20] and shape of the pipe [21–23], and the means of fixation of the pipe [24]. Moreover, the temperature, pressure, and gas content can also affect the wave speed in a pipe system [25, 26]. In fact, it is difficult to estimate accurately the shock wave speed. Since the presentation of several basic formulas by Wylie and Streeter [18] for the shock wave speeds, they are widely used in common pipes and conduits. However, it is inapplicable for some specific fluids and composite material pipes, and so some modified approaches have been developed to estimate the shock wave speed for various areas. Sun et al. [24] presented the water hammer speed of fiber-reinforced plastic composite pipes based on three different fixed means. Hachem and Schleiss [27] studied the influence of local wall stiffness decrease on the shock wave speed by experiments, and the result showed the transient pressure change with the wave speed. Han et al. [28] and Zhou et al. [29] studied the shock wave speed of a slurry flow carrying solid particles. Hadj-Taïeb and Lili [30] and Ando et al. [31] studied the shock wave propagation through a deformable tube. Lee and Pejovic [26] studied the influence of air on the similarity of hydraulic transients and vibrations. Soares et al. [22] and Apollonio et al. [23] investigated separately the hydraulic transient properties of a pipeline made of plastic material. Mitosek [32] measured the shock wave speed of PVC pipes. Especially for polyethylene pipe, Covas et al. [19, 33, 34] developed a viscoelastic model to take account of the influences of viscoelasticity on transient processes. Evangelista et al. [20] simulated complex plastic pipes system by viscoelastic model and showed that the viscoelasticity needs to be considered in the transient response of a plastic pipe. These previous achievements are significant, but more accurate simulations are required to prevent the water hammer. However, considering the influence of complicated steel-mesh structures, it is presently still difficult to reasonably estimate the wave speed of a large PE pipe with steel-mesh reinforcement. In order to improve the transient simulation of a PE pipe with steel-mesh reinforcement, we propose an improved wave speed formula for reinforced PE pipes by incorporating the equivalent modulus. Also, we conduct a field test to validate the approach. Based on the proposed approach, the transient responses are simulated and compared for PE pipes with various steel-mesh reinforcement ratios. Then the influence of the steel mesh on the transient response is evaluated and discussed. The result shows that the steel mesh can increase the shock wave speed and cause more serve hydraulic surges in the pipeline. The proposed formula is significant to revise the wave speed and improve the transient simulation of steel-mesh-reinforced PE pipe.


## 2. Basic Shock Wave Speed Formula

Wave speed can primarily affect the wave propagation and transient pressure, and it is an important parameter for water hammer simulations, as shown in the simplified closed expression equations for water hammer simulation [1]:(1)∂v∂t+v∂v∂x+g∂h∂x+fvv2D=0,(2)∂h∂t+v∂h∂x-vsin⁡β+a02g∂v∂x=0.For PE pipe, considering the additional viscoelastic term [20, 34], the improved continuity equation can be written as(3)∂h∂t+v∂h∂x-vsin⁡β+a02g∂v∂x+2a02gdεrdt=0.Generally, the shock wave speed can be determined by Young’s modulus of the fluid and the wall material for some regular water supply pipes. Based on basic water hammer equation and mass conservation equations, [1] provided a basic shock wave speed formula:(4)a02=AρdpdA.For general elastic pipes, the shock wave speed formula can be expressed as [1](5)a0=K/ρ1+α2R/δK/E.This formula is widely applicable to circular pipes. However, for some complex pipe materials, this formula may need to be improved. As seen in this formula, the modulus of elasticity is an important parameter that determines the wave speed. A reinforced PE pipe has a composite structure, where the steel mesh and polyethylene play different roles. Accordingly, the formula cannot directly determine the wave speed in reinforced PE. Based on the conventional formula, therefore, we will analyze the wave speed for reinforced PE pipe with steel-mesh reinforcement.


## 3. Wave Speed Analysis in Reinforced PE Pipes

### 3.1. Composition of Steel-Mesh-Reinforced PE Pipe

A general PE pipe is made of polyethylene. In order to improve the strength, steel meshes are embedded in the pipe wall. Compared to the common PE pipe, reinforced PE pipes have more strength and are inflexible. Figures1 and 2 illustrate the basic structure of a PE pipe with steel-mesh reinforcement. The pipe wall consists of PE and steel meshes. The modulus of the mesh is much larger than that of PE. Thus, the shock wave speed changes greatly with the steel-mesh reinforcement. It is necessary to revise the shock wave speed due to the steel meshes, in order to obtain a more reasonable transient simulation result.Figure 1
Composition of a PE pipe with steel-mesh reinforcement.Figure 2
Schematic of a steel-mesh-reinforced PE pipe.


### 3.2. Force Analysis of the Pipe Wall under Small Deformation

In order to determine the shock wave speed of a reinforced PE pipe, an equivalent modulus is proposed for the composite reinforced PE pipe wall. Based on the conventional wave speed formula, the equivalent modulus is used as the modified modulus of pipe wall. In fact, stress and strain occur simultaneously when a pipe is subjected to a water hammer wave. Thus, the transient wave propagation represents the variations of fluid and pipe in pressure and deformation.Before deriving the equivalent modulus, we will make two assumptions.(1) The deformation is sufficiently small and can be considered as elasticity in instantaneous time, so that the stress and strain can satisfy Hooke’s law for both steel meshes and PE. (2) For a finite deformation, there is no relative displacement between the steel thread and the PE material. In other words, steel and PE have the same deformation to bear the extra pressure due to the water hammer pressure.We select a unit pipe of specific length, around which a steel thread is wound. As shown in Figure3, the axial length of the pipe is calculated as(6)B=CPcot⁡θ,(7)CP=2πR+0.5δ.Figure 3
Schematic of unfolded pipe wall structure.Figure4 shows the unfolded pipe wall. When an internal hydraulic pressure p acts on the pipe wall, the tensile force is equal to the total horizontal component:(8)P=DBp.Figure 4
Force analysis of the unfolded pipe wall.In the specific length, a spiral stirrup is considered as an equivalent hoop. For a differential deformation, the total circular force in steel mesh can be expressed as(9)Ts=2nAsεEssin⁡θ.Simultaneously, the circular force in the PE material can be expressed as(10)Tp=2Bδ-nAssin⁡θEpε.


### 3.3. Equivalent Modulus of PE Pipe with Steel-Mesh Reinforcement

Based on the above assumptions, the force of pipe wall can be expressed as(11)T=2nAsEsεsin⁡θ+Bδ-nAssin⁡θEpε.If we define an equivalent modulusEe, then the total force can be written as(12)Fe=2BδEeε.Equations (10) and (11) represent the equal force in magnitude to (7) for an equilibrium state analysis. Accordingly (13)2BδεEe=2nAsEsεsin⁡θ+Bδ-nAssin⁡θEpε.Then the equivalent modulus of the pipe wall can be determined as(14)Ee=nAsEssin⁡θ+Bδ-nAs/sin⁡θEpBδ.We define the pipe wall reinforcement ratio as(15)ps=nAsBδ.The reinforcement ratio represents the ratio between the steel areas and the wall section areas. Substituting it into (13), the equivalent elastic modulus can be expressed as(16)Ee=pssin⁡θEs+1-pscsc⁡θ.


### 3.4. Improved Wave Speed Formula for a PE Pipe with Steel-Mesh Reinforcement

Equation (15) provides the equivalent Young’s modulus in terms of the reinforcement ratio of the pipe wall. That is to say, we can deal with the reinforced PE pipe wall as a composite material with an equivalent elastic modulus. Accordingly, the wave speed can be expressed as(17)a0=K/ρ1+α2R/δK/pssin⁡θEs+1-pscsc⁡θEp.Provided we know the number and winding angle of the spiral steel wire in pipe cross section, the improved formula can calculate the shock wave speed for a PE pipe with steel-mesh reinforcement.


### 3.5. Validation of the Wave Speed Formula

In fact, it is difficult to conduct an experiment to obtain the shock wave speed for a large PE pipe with steel-mesh reinforcement in the field. Fortunately, a practical project supports the field to validate the proposed formula. As a result, the experiment is subject to the practical scale and layout of the project. Figure5 shows the schematic of pipeline and stations in the field test, where s1 is a 20.16 km PE pipe, s2 is a 0.4 km steel pipe, and s3 is a 0.1 km steel pipe. Moreover, the PE pipe is 0.232 m in internal radius and 0.018 m in thickness.Figure 5
Schematic of wave speed measurement in the field.In the experiment, the control valve will give rise to a negative pressure wave by sudden opening and discharging water from the main pipe at initial timet0, and then the negative wave propagates along the pipe to the stations. Stations 1 and 2 detect the wave at times t1 and t2, respectively. According to the time and wave speed, the following equations can be written:(18)s1ap+s2as=t1-t0,(19)s1ap+s2as+s3as=t2-t0.According tot0, t1, and t2, the experimental wave speed can be obtained as follows:(20)ap=s1s3t1-t0s2+s3-t2-t0s2.Figure6 shows the pressure response at station 1 and station 2, respectively. The beginning time of the impulse t0 is 29.61 s, and the response times t1 and t2 are, respectively, 83.109 and 83.214 s at station 1 and station 2. According to (19), the experimental wave speed can be obtained.Figure 6
Pressure surge in the stations.As shown in Table1, the experimental wave speed is 379.8 m/s. According to the steel-mesh density and winding angle of the spiral steel wire, the proposed theoretical formula yields a collective wave speed 385.0 m/s. The error is about 1.37%. The comparison shows that the computational result agrees well with the measurement. The proposed formula can provide a reasonable result for the PE pipe with steel-mesh reinforcement.Table 1
Comparison of test and calculation.
R (mm) δ (mm) p s (%) θ E p (GPa) E s (GPa) E e (GPa) a n (m/s) a m (m/s) ε a 232 18 1.48% 42.26° 1.43 207 3.43 385.0 379.8 1.37%


### 3.6. Influence of Reinforcement Ratio on Wave Speeds of a PE Pipe

As is well known, the elastic modulus can greatly affect the shock wave speed. Since the reinforcement can primarily affect the elastic modulus, it can also change the shock wave speed. Next, the influence of the steel-mesh density on the wave speed is studied, including the reinforcement ratios and winding angles. Based on the proposed formula, Figure7 shows the influence of the reinforcement ratio on shock wave speed for three different PE pipe walls. For a given pipe diameter, wall thickness, and winding angle of the steel wire, the wave speed increases with the reinforcement ratio. If we define reinforcement ratio as 0, it will indicate actually the original PE material pipe without steel meshes, and the wave speed is only about 250 m/s. Obviously, the wave speed increases with the reinforcement ratio. Moreover, for a given diameter, wall thickness, and reinforcement ratio, the wave speed also increases with the winding angle of steel wire, as shown in Figure 8. In fact, for a specified reinforcement ratio in the section, the density of the steel mesh increases with the winding angle along the axial direction. It shows that the wave speed changes because of the steel mesh enhancing the elastic modulus of the PE pipe. In practice, the range of the reinforcement ratio is 1%-2%, and the winding angle of a rhombic metal mesh is about 42.26°; accordingly, the wave speed in the PE pipe may increase to 300–400 m/s.Figure 7
Influence of reinforcement ratio on the wave speed.Figure 8
Influence of spiral winding angle on the wave speed.


## 3.1. Composition of Steel-Mesh-Reinforced PE Pipe

A general PE pipe is made of polyethylene. In order to improve the strength, steel meshes are embedded in the pipe wall. Compared to the common PE pipe, reinforced PE pipes have more strength and are inflexible. Figures1 and 2 illustrate the basic structure of a PE pipe with steel-mesh reinforcement. The pipe wall consists of PE and steel meshes. The modulus of the mesh is much larger than that of PE. Thus, the shock wave speed changes greatly with the steel-mesh reinforcement. It is necessary to revise the shock wave speed due to the steel meshes, in order to obtain a more reasonable transient simulation result.Figure 1
Composition of a PE pipe with steel-mesh reinforcement.Figure 2
Schematic of a steel-mesh-reinforced PE pipe.


## 3.2. Force Analysis of the Pipe Wall under Small Deformation

In order to determine the shock wave speed of a reinforced PE pipe, an equivalent modulus is proposed for the composite reinforced PE pipe wall. Based on the conventional wave speed formula, the equivalent modulus is used as the modified modulus of pipe wall. In fact, stress and strain occur simultaneously when a pipe is subjected to a water hammer wave. Thus, the transient wave propagation represents the variations of fluid and pipe in pressure and deformation.Before deriving the equivalent modulus, we will make two assumptions.(1) The deformation is sufficiently small and can be considered as elasticity in instantaneous time, so that the stress and strain can satisfy Hooke’s law for both steel meshes and PE. (2) For a finite deformation, there is no relative displacement between the steel thread and the PE material. In other words, steel and PE have the same deformation to bear the extra pressure due to the water hammer pressure.We select a unit pipe of specific length, around which a steel thread is wound. As shown in Figure3, the axial length of the pipe is calculated as(6)B=CPcot⁡θ,(7)CP=2πR+0.5δ.Figure 3
Schematic of unfolded pipe wall structure.Figure4 shows the unfolded pipe wall. When an internal hydraulic pressure p acts on the pipe wall, the tensile force is equal to the total horizontal component:(8)P=DBp.Figure 4
Force analysis of the unfolded pipe wall.In the specific length, a spiral stirrup is considered as an equivalent hoop. For a differential deformation, the total circular force in steel mesh can be expressed as(9)Ts=2nAsεEssin⁡θ.Simultaneously, the circular force in the PE material can be expressed as(10)Tp=2Bδ-nAssin⁡θEpε.


## 3.3. Equivalent Modulus of PE Pipe with Steel-Mesh Reinforcement

Based on the above assumptions, the force of pipe wall can be expressed as(11)T=2nAsEsεsin⁡θ+Bδ-nAssin⁡θEpε.If we define an equivalent modulusEe, then the total force can be written as(12)Fe=2BδEeε.Equations (10) and (11) represent the equal force in magnitude to (7) for an equilibrium state analysis. Accordingly (13)2BδεEe=2nAsEsεsin⁡θ+Bδ-nAssin⁡θEpε.Then the equivalent modulus of the pipe wall can be determined as(14)Ee=nAsEssin⁡θ+Bδ-nAs/sin⁡θEpBδ.We define the pipe wall reinforcement ratio as(15)ps=nAsBδ.The reinforcement ratio represents the ratio between the steel areas and the wall section areas. Substituting it into (13), the equivalent elastic modulus can be expressed as(16)Ee=pssin⁡θEs+1-pscsc⁡θ.


## 3.4. Improved Wave Speed Formula for a PE Pipe with Steel-Mesh Reinforcement

Equation (15) provides the equivalent Young’s modulus in terms of the reinforcement ratio of the pipe wall. That is to say, we can deal with the reinforced PE pipe wall as a composite material with an equivalent elastic modulus. Accordingly, the wave speed can be expressed as(17)a0=K/ρ1+α2R/δK/pssin⁡θEs+1-pscsc⁡θEp.Provided we know the number and winding angle of the spiral steel wire in pipe cross section, the improved formula can calculate the shock wave speed for a PE pipe with steel-mesh reinforcement.


## 3.5. Validation of the Wave Speed Formula

In fact, it is difficult to conduct an experiment to obtain the shock wave speed for a large PE pipe with steel-mesh reinforcement in the field. Fortunately, a practical project supports the field to validate the proposed formula. As a result, the experiment is subject to the practical scale and layout of the project. Figure5 shows the schematic of pipeline and stations in the field test, where s1 is a 20.16 km PE pipe, s2 is a 0.4 km steel pipe, and s3 is a 0.1 km steel pipe. Moreover, the PE pipe is 0.232 m in internal radius and 0.018 m in thickness.Figure 5
Schematic of wave speed measurement in the field.In the experiment, the control valve will give rise to a negative pressure wave by sudden opening and discharging water from the main pipe at initial timet0, and then the negative wave propagates along the pipe to the stations. Stations 1 and 2 detect the wave at times t1 and t2, respectively. According to the time and wave speed, the following equations can be written:(18)s1ap+s2as=t1-t0,(19)s1ap+s2as+s3as=t2-t0.According tot0, t1, and t2, the experimental wave speed can be obtained as follows:(20)ap=s1s3t1-t0s2+s3-t2-t0s2.Figure6 shows the pressure response at station 1 and station 2, respectively. The beginning time of the impulse t0 is 29.61 s, and the response times t1 and t2 are, respectively, 83.109 and 83.214 s at station 1 and station 2. According to (19), the experimental wave speed can be obtained.Figure 6
Pressure surge in the stations.As shown in Table1, the experimental wave speed is 379.8 m/s. According to the steel-mesh density and winding angle of the spiral steel wire, the proposed theoretical formula yields a collective wave speed 385.0 m/s. The error is about 1.37%. The comparison shows that the computational result agrees well with the measurement. The proposed formula can provide a reasonable result for the PE pipe with steel-mesh reinforcement.Table 1
Comparison of test and calculation.
R (mm) δ (mm) p s (%) θ E p (GPa) E s (GPa) E e (GPa) a n (m/s) a m (m/s) ε a 232 18 1.48% 42.26° 1.43 207 3.43 385.0 379.8 1.37%


## 3.6. Influence of Reinforcement Ratio on Wave Speeds of a PE Pipe

As is well known, the elastic modulus can greatly affect the shock wave speed. Since the reinforcement can primarily affect the elastic modulus, it can also change the shock wave speed. Next, the influence of the steel-mesh density on the wave speed is studied, including the reinforcement ratios and winding angles. Based on the proposed formula, Figure7 shows the influence of the reinforcement ratio on shock wave speed for three different PE pipe walls. For a given pipe diameter, wall thickness, and winding angle of the steel wire, the wave speed increases with the reinforcement ratio. If we define reinforcement ratio as 0, it will indicate actually the original PE material pipe without steel meshes, and the wave speed is only about 250 m/s. Obviously, the wave speed increases with the reinforcement ratio. Moreover, for a given diameter, wall thickness, and reinforcement ratio, the wave speed also increases with the winding angle of steel wire, as shown in Figure 8. In fact, for a specified reinforcement ratio in the section, the density of the steel mesh increases with the winding angle along the axial direction. It shows that the wave speed changes because of the steel mesh enhancing the elastic modulus of the PE pipe. In practice, the range of the reinforcement ratio is 1%-2%, and the winding angle of a rhombic metal mesh is about 42.26°; accordingly, the wave speed in the PE pipe may increase to 300–400 m/s.Figure 7
Influence of reinforcement ratio on the wave speed.Figure 8
Influence of spiral winding angle on the wave speed.


## 4. Effects of Steel Mesh on Transient Responses of the PE Pipe

### 4.1. Selection of Transient Simulation Model

Considering the viscoelasticity of PE materials, the viscoelastic model is needed to simulate the hydraulic transient processes. Covas et al. [34] developed a conventional viscoelastic model, and the model is also well verified in complex plastic pipes system [20]. Given the wave speed, (1) and (3) form a closed system for the water hammer simulations. Considering the viscoelastic model [34] and coupling the method of characteristics [1], the modified MOC equation systems can be expressed as(21)ga0dhdt+dvdt+fvv2D+2a0∂εr∂t=0dx/dt=+a0,-ga0dhdt+dvdt+fvv2D+2a0∂εr∂t=0dx/dt=-a0.Integral along the characteristic line, the system of equations can be converted as(22)hi,t+Δt+a0gvi,t+Δt+2a02Δtg∂εr∂t-CF=0,hi,t+Δt-a0gvi,t+Δt+2a02Δtg∂εr∂t-CR=0.In the equations,CF=hi-1,t+vi-1,ta0/g-sgn⁡(v)fΔxvi-1,t2/(2gD) and CR=hi+1,t-vi+1,ta0/g+sgn⁡(v)fΔxvi+1,t2/(2gD).About the term∂εr/∂t, Covas et al. [34] have proposed an applied model. According to the model, the term can be expressed as(23)∂εri,t+Δt∂t=CShi,t+Δt+CT,where CS=(αρa02D/δ)∑k=1nJk(1-ek),ek=e-Δt/τk,CT=2(a02Δt/g)∑k=1n(Jkek/τk)(αD/2δ)(hi,t-hi,0)-(Jk1-ek/Δt)(αD/2δ)ρghi,t-ekε~rk(i,t),  Fi,t=(αD/2δ)ρghi,t-hi,0,  and  ε~rk(i,t)=JkFi,t-JkekFi,t-Δt-(Jkτk1-ek/Δt)(Fi,t-Fi,t-Δt)+ekε~rk(i,t-Δt).For a present step, all values at initial timet are known. A system of closed equations is obtained. Then the solution of the simultaneous equations can be written as(24)hi,t+Δt=0.5CF+CR+2CT1+CS,vi,t+Δt=hi,t+Δt+CS+CT-CRga0.To simulate the transient process with viscoelastic properties, a six-element Kelvin-Voigt model and basic creep functions of polyethylene [34] are adopted. In order to analyze the influence of the elastic steel mesh, an extra spring is added. Figure 9 shows the improved six-element Kelvin-Voigt model. In the model, the instantaneous elastic is modified to the compound of the PE instantaneous elastic and the steel instantaneous elastic according to their contribution factors. In the figure, Ees=pssin⁡θEs and Eep=(1-pscsc⁡θ)Ep≈Ep.Figure 9
Six-element Kelvin-Voigt viscoelastic model with steel mesh.For viscoelastic solid, the basic Kelvin-Voigt creep function can be written as [35](25)Jt=J0+∑k=1NJk1-e-t/τk.Covas et al. [34] have fixed several creep functions for typical polyethylene materials. In order to describe the behavior of the polyethylene, especially, the basic creep coefficient is approximated by the six-element Kelvin-Voigt model [34], which has the corresponding instantaneous creep compliance J0=0.7×10-9Pa-1. For PE pipe with steel mesh, the improved model is modified to J0=1/(Eep+Ees), as shown in Figure 9. Considering the influence of the steel mesh, Figure 10 shows the modified creep functions for different reinforcement ratios.Figure 10
Modified creep function for PE pipe with steel mesh.


### 4.2. Pipe and Material Parameters

As shown in the previous section, the steel mesh can greatly increase the shock wave speed; meanwhile, it can also cause more severe transient pressure. In order to analyze the influence of the steel mesh on the transient response, various reinforcement ratios are considered in the same scale pipeline system. As shown in Figure11, the system is composed of a booster pump, a control valve, the main pipe, and the downstream reservoir. In the example, the pipe is 0.5 m in external diameter, 0.018 m in wall thickness, and 2500 m in length. Four types of reinforcements are simulated to analyze the influence of the steel mesh on the transient response. The reinforcement ratios are 0% (the pure PE materials), 0.37%, 0.74%, and 1.48% separately. Table 2 shows the equivalent elastic modulus and shock wave speed parameters for next transient simulation according to the proposed method.Table 2
Wave speeds of PE pipes with various reinforcement ratios.
PE R (mm) δ (mm) p s (%) θ (°) E p (GPa) E s (GPa) E e (GPa) a 0 (m/s) Φ500 232 18 0.00 — 1.43 207 1.43 253.04 Φ500 232 18 0.37 42.26 1.43 207 1.94 292.82 Φ500 232 18 0.74 42.26 1.43 207 2.44 327.10 Φ500 232 18 1.48 42.26 1.43 207 3.45 384.95Figure 11
Basic model of a water supply PE pipe system.


### 4.3. Influence of Steel-Mesh Reinforcement on Transient Response of PE Pipe

In the example, the hydraulic transients of pump failure were numerically simulated and compared for the water supply system. Figure12 compares the transient response by the classic model and the viscoelastic model. The result shows that viscoelasticity has a great effect on transient process; it needs to take account of viscoelasticity in transient simulation of PE pipe with steel mesh. Figure 13 shows the transient pressure waves of the same scale PE pipe with different reinforcement ratios. As seen in the figure, the negative pressure is the largest when the PE pipe has the reinforcement ratio of 1.48%. Conversely, the negative pressure is the least for the PE pipe without the steel mesh. Obviously, the transient intensities increase with the reinforcement ratio, as well as the wave frequency and amplitude. Accordingly, Figure 14 shows the distribution of extreme transient pressure along the pipeline. Compared with the original extreme pressures, Table 3 shows that the steel-mesh reinforcement has increased the amplitudes of the transient pressure along the entire pipeline.Table 3
Extreme pressure with various reinforcement ratios (viscoelastic model).
p s (%) a (m/s) h 0 (m) h m i n (m) h m a x (m) A m (m) 0.00 253.04 49.95 37.84 56.75 12.11 0.37 292.82 49.95 36.88 56.91 13.07 0.74 327.10 49.95 36.16 56.92 13.79 1.48 384.95 49.95 35.20 56.78 14.75Figure 12
Comparison of pressure waves by different transient models.Figure 13
Pressure waves for various reinforcement ratios.Figure 14
Extreme pressure distributions along the pipeline.Figure15 shows the changes of the extreme pressures and amplitudes with the reinforcement ratio. As seen in this figure, the pressure surge increases with the reinforcement ratio. The result shows that steel-mesh reinforcement can cause more severe transient response due to increase in the shock wave speed of the PE pipes. Therefore, the modified wave speed formula and viscoelastic model can improve the transient simulation of the PE pipe with steel-mesh reinforcement for a water supply system.Figure 15
Transient intensions for various reinforcement ratios.


## 4.1. Selection of Transient Simulation Model

Considering the viscoelasticity of PE materials, the viscoelastic model is needed to simulate the hydraulic transient processes. Covas et al. [34] developed a conventional viscoelastic model, and the model is also well verified in complex plastic pipes system [20]. Given the wave speed, (1) and (3) form a closed system for the water hammer simulations. Considering the viscoelastic model [34] and coupling the method of characteristics [1], the modified MOC equation systems can be expressed as(21)ga0dhdt+dvdt+fvv2D+2a0∂εr∂t=0dx/dt=+a0,-ga0dhdt+dvdt+fvv2D+2a0∂εr∂t=0dx/dt=-a0.Integral along the characteristic line, the system of equations can be converted as(22)hi,t+Δt+a0gvi,t+Δt+2a02Δtg∂εr∂t-CF=0,hi,t+Δt-a0gvi,t+Δt+2a02Δtg∂εr∂t-CR=0.In the equations,CF=hi-1,t+vi-1,ta0/g-sgn⁡(v)fΔxvi-1,t2/(2gD) and CR=hi+1,t-vi+1,ta0/g+sgn⁡(v)fΔxvi+1,t2/(2gD).About the term∂εr/∂t, Covas et al. [34] have proposed an applied model. According to the model, the term can be expressed as(23)∂εri,t+Δt∂t=CShi,t+Δt+CT,where CS=(αρa02D/δ)∑k=1nJk(1-ek),ek=e-Δt/τk,CT=2(a02Δt/g)∑k=1n(Jkek/τk)(αD/2δ)(hi,t-hi,0)-(Jk1-ek/Δt)(αD/2δ)ρghi,t-ekε~rk(i,t),  Fi,t=(αD/2δ)ρghi,t-hi,0,  and  ε~rk(i,t)=JkFi,t-JkekFi,t-Δt-(Jkτk1-ek/Δt)(Fi,t-Fi,t-Δt)+ekε~rk(i,t-Δt).For a present step, all values at initial timet are known. A system of closed equations is obtained. Then the solution of the simultaneous equations can be written as(24)hi,t+Δt=0.5CF+CR+2CT1+CS,vi,t+Δt=hi,t+Δt+CS+CT-CRga0.To simulate the transient process with viscoelastic properties, a six-element Kelvin-Voigt model and basic creep functions of polyethylene [34] are adopted. In order to analyze the influence of the elastic steel mesh, an extra spring is added. Figure 9 shows the improved six-element Kelvin-Voigt model. In the model, the instantaneous elastic is modified to the compound of the PE instantaneous elastic and the steel instantaneous elastic according to their contribution factors. In the figure, Ees=pssin⁡θEs and Eep=(1-pscsc⁡θ)Ep≈Ep.Figure 9
Six-element Kelvin-Voigt viscoelastic model with steel mesh.For viscoelastic solid, the basic Kelvin-Voigt creep function can be written as [35](25)Jt=J0+∑k=1NJk1-e-t/τk.Covas et al. [34] have fixed several creep functions for typical polyethylene materials. In order to describe the behavior of the polyethylene, especially, the basic creep coefficient is approximated by the six-element Kelvin-Voigt model [34], which has the corresponding instantaneous creep compliance J0=0.7×10-9Pa-1. For PE pipe with steel mesh, the improved model is modified to J0=1/(Eep+Ees), as shown in Figure 9. Considering the influence of the steel mesh, Figure 10 shows the modified creep functions for different reinforcement ratios.Figure 10
Modified creep function for PE pipe with steel mesh.


## 4.2. Pipe and Material Parameters

As shown in the previous section, the steel mesh can greatly increase the shock wave speed; meanwhile, it can also cause more severe transient pressure. In order to analyze the influence of the steel mesh on the transient response, various reinforcement ratios are considered in the same scale pipeline system. As shown in Figure11, the system is composed of a booster pump, a control valve, the main pipe, and the downstream reservoir. In the example, the pipe is 0.5 m in external diameter, 0.018 m in wall thickness, and 2500 m in length. Four types of reinforcements are simulated to analyze the influence of the steel mesh on the transient response. The reinforcement ratios are 0% (the pure PE materials), 0.37%, 0.74%, and 1.48% separately. Table 2 shows the equivalent elastic modulus and shock wave speed parameters for next transient simulation according to the proposed method.Table 2
Wave speeds of PE pipes with various reinforcement ratios.
PE R (mm) δ (mm) p s (%) θ (°) E p (GPa) E s (GPa) E e (GPa) a 0 (m/s) Φ500 232 18 0.00 — 1.43 207 1.43 253.04 Φ500 232 18 0.37 42.26 1.43 207 1.94 292.82 Φ500 232 18 0.74 42.26 1.43 207 2.44 327.10 Φ500 232 18 1.48 42.26 1.43 207 3.45 384.95Figure 11
Basic model of a water supply PE pipe system.


## 4.3. Influence of Steel-Mesh Reinforcement on Transient Response of PE Pipe

In the example, the hydraulic transients of pump failure were numerically simulated and compared for the water supply system. Figure12 compares the transient response by the classic model and the viscoelastic model. The result shows that viscoelasticity has a great effect on transient process; it needs to take account of viscoelasticity in transient simulation of PE pipe with steel mesh. Figure 13 shows the transient pressure waves of the same scale PE pipe with different reinforcement ratios. As seen in the figure, the negative pressure is the largest when the PE pipe has the reinforcement ratio of 1.48%. Conversely, the negative pressure is the least for the PE pipe without the steel mesh. Obviously, the transient intensities increase with the reinforcement ratio, as well as the wave frequency and amplitude. Accordingly, Figure 14 shows the distribution of extreme transient pressure along the pipeline. Compared with the original extreme pressures, Table 3 shows that the steel-mesh reinforcement has increased the amplitudes of the transient pressure along the entire pipeline.Table 3
Extreme pressure with various reinforcement ratios (viscoelastic model).
p s (%) a (m/s) h 0 (m) h m i n (m) h m a x (m) A m (m) 0.00 253.04 49.95 37.84 56.75 12.11 0.37 292.82 49.95 36.88 56.91 13.07 0.74 327.10 49.95 36.16 56.92 13.79 1.48 384.95 49.95 35.20 56.78 14.75Figure 12
Comparison of pressure waves by different transient models.Figure 13
Pressure waves for various reinforcement ratios.Figure 14
Extreme pressure distributions along the pipeline.Figure15 shows the changes of the extreme pressures and amplitudes with the reinforcement ratio. As seen in this figure, the pressure surge increases with the reinforcement ratio. The result shows that steel-mesh reinforcement can cause more severe transient response due to increase in the shock wave speed of the PE pipes. Therefore, the modified wave speed formula and viscoelastic model can improve the transient simulation of the PE pipe with steel-mesh reinforcement for a water supply system.Figure 15
Transient intensions for various reinforcement ratios.


## 5. Discussion

For the composite structure of a PE pipe with steel-mesh reinforcement, an equivalent instantaneous elastic modulus is used to determine the water hammer properties of the reinforced system. Then an improved shock wave speed formula is proposed for the reinforced PE pipe based on the reinforcement ratios and winding angles of spiral steel threads. Accordingly, the formula is validated by a field measurement in a submarine PE pipe with steel-mesh reinforcement. Finally, transient simulation is applied to analyze the influence of the steel mesh on the hydraulic transients of the reinforced PE pipe. Compared to the original uniform PE material, the steel-mesh-reinforced PE material has a higher elastic modulus and a larger shock wave speed. Therefore, it can cause more severe transient pressure. For example, the negative pressure, wave amplitudes, and frequencies greatly increase with the steel-mesh reinforcement ratio. Viscoelastic model especially is needed to take account of the viscoelasticity of the PE pipe. Obviously, a higher wave speed increases the water hammer risks due to the steel-mesh reinforcement in a PE pipe. It is necessary to evaluate the water hammer pressure more seriously because of the effects of the steel mesh in PE pipe system. The viscoelastic model and the proposed shock wave speed can yield a more reasonable transient simulation in the PE pipe with steel-mesh reinforcement.


## 6. Conclusion

A steel mesh can greatly affect the transient intensity in a PE pipe water supply system, because it increases the shock wave speed in the PE pipe. Generally, a steel mesh can cause more serious water hammer peak pressure and hydraulic fluctuation; thus more consideration on the reinforcement mesh is necessary to protect the pipeline from water hammer damage. It is worth reevaluating the shock wave speed according to the density and winding angle of the steel mesh. The equivalent elastic modulus is presented to indicate the material properties of steel-mesh-reinforced PE pipe wall. Accordingly, an improved wave speed formula is proposed to evaluate reasonably the water hammer speed for the reinforced PE pipe. It agrees well with a field measurement in a submarine reinforced PE pipeline system. Given the reinforcement ratio and winding angle of the spiral steel threads, the improved formula can conveniently be applied to the calculation of water hammer speed. Moreover, viscoelasticity has a great effect on transient process; it needs to take account of viscoelasticity in transient simulation of PE pipe with steel mesh. The comparison of various reinforcement ratios shows that the steel mesh can increase the transient pressure, as well as the wave frequency and amplitude. Consequently, it is significant to evaluate reasonably the wave speed by the proposed formula and improve the transient simulation by viscoelastic model for water hammer prediction and prevention in a reinforced PE pipe system.



---

*Source: 8705031-2016-05-25.xml*