Numerical Study on Acoustic Resonance Excitation in Closed Side Branch Pipeline Conveying Natural Gas

Flow-induced acoustic resonance in the closed side branch of a natural gas pipeline can cause intensive vibration which threatens the safe operation of the pipeline. Accurately modeling this excitation process is necessary for a workable understanding of the genetic mechanism to resolve this problem. A realizable k-εDelayed Detached Eddy Simulation (DDES) model was conducted in this study to numerically simulate the acoustic resonance problem. *e model is shown to accurately capture the acoustic resonance phenomenon and self-excited vibration characteristics with low calculation cost. *e pressure pulsation component of the acoustic resonance frequency is gradually amplified and transformed into a narrowband dominant frequency in the process of acoustic resonance excitation, forming a so-called “frequency lock-in phenomenon.”*e gas is pressed into and out of the branch in sinusoidal mode during excitation. *e first-order frequency, single vortex moves at the branch inlet following the same pattern. A quarter wavelength steady standing wave forms in the branch. *e mechanism and characteristics presented in this paper may provide guidelines for developing new excitation suppression methods.


Introduction
Natural gas is a clean, efficient, and high-quality energy resource which is a preferred fuel by many countries across the globe [1,2]. e mileage of natural gas pipeline construction projects has rapidly increased alongside the demand for fuel in recent decades. For example, the length of natural gas pipeline inservice has reached 81,000 km in China; 3,233 km of this length was newly built in 2019. e closed side branch of a natural gas pipeline is a commonly used unit of gas transmission. Flowinduced acoustic resonance occurs as gas flows through the mouth of a closed branch at a certain velocity, causing intensive and continuous gas pressure pulsations that may result severe pipeline vibration and even fatigue damage [3][4][5][6]. ese pulsations are also noisy and intense enough to endanger the physical and mental health of workers. A serious acoustic resonance event occurred in the WKC2 gas transmission station of China in 2012 which fractured a tee connector (Figure 1). ere is substantial practical importance in accurately, comprehensively understanding the mechanism of acoustic resonance and the most effective approaches to suppressing it [7][8][9][10][11][12]. e essential mechanism of acoustic resonance in the closed side branch is the coupling effect of the flow field and sound field. e shear layer detaches from the upstream pipe wall, exciting a sound vortex that travels downstream and is coupled with the gas in the side branch; this causes acoustic vortex convection along the branch. e sound vortex is gradually excited as it detaches. e convection of the vortex (vortices) generates and absorbs sound energy to maintain continuous pressure pulsation [13].
e key to effectively understanding the shear layer detachment process lies in accurately modeling the sound source excitation and hydrodynamic modal in the branch mouth. Deriving a solvable differential equation of sound source excitation analytically is almost impossible, as the turbulence flow in the main pipe is excessively complex. e vortex acoustic theory developed by Powell [14] and Howe [15] has been widely used in the acoustic field, including acoustic resonance problems, but it is still difficult to solve vortex acoustic equations. Michalke [16] used linear stability theory to model a shear layer flow. Karlsson andÅbom [17] developed expressions for the side branch acoustic impedance under various velocities, and investigated the interaction between acoustic impedance and sound source excitation. Considering the nonlinearity of this problem, Graf and Ziada [18] proposed a semiempirical model for predicting the amplitude and phase of pressure fluctuations in tandem side branches and coaxial branches. eir model includes sound source excitation as the nonlinear pressure difference between the upper side and lower side of the detached shear layer in the branch mouth and can predict the acoustic amplitude and frequency of the pressure in the two branches.
Visual flow experiments and computational fluid dynamics (CFD) are the most common approaches to hydrodynamic modal analysis. In visual experiments, the morphology of the acoustic vortex and flow patterns of the fluid during the acoustic resonance can be clearly observed via high-speed camera. Ziada [19], for example, found that the mechanism of acoustic resonance is different from that of nonresonant oscillation of impact shear flow; visual analysis revealed that pressure fluctuation in closed side branches was caused by the coupling of the unstable detachment of shear layer from the upstream wall and standing acoustic waves in the side branch. During the unstable detachment of the shear layer, an acoustic vortex is excited as it is trapped in the mouth of the side branch. Li et al. [20][21][22] studied the flow patterns of acoustic resonance in tandem side branches and coaxial side branches using a high timeresolved particle image velocimetry (PIV) technique. ey captured the velocity distribution in the flow field and constructed a phase delay diagram. Salt et al. [23,24] also used PIV to observe the acoustic energy transfer pattern during an acoustic resonance cycle as per the sound source generation and absorption with shear layer oscillation. Xiao et al. [25] similarly obtained a multimodal image of shear layer detachment under the acoustic resonance of a square single closed side branch.
Flow patterns in closed side branches can be analyzed effectively via CFD methods. For example, Tamura et al. [26] modeled the acoustic resonance in three side branches based on a three-dimensional (3D) model using the finite difference lattice Boltzmann method. Kriesels et al. [5] simulated the gas flow and calculated the sound power generated by a vortex using the vortex-blob method based on a 2D potential flow model. Radavich et al. [27] simulated the 2D flow of a single closed side branch by solving Navier-Stokes equations of unsteady and compressible turbulence. Dequand et al. [28] simulated a 2D compressible coaxial closed side branch flow via Euler equation. Numerical simulations accurately predicted the range of frequency lock-in velocity and flow patterns.
Previous research has mainly centered on the formation conditions, pressure fluctuation amplitudes, and frequencies of acoustic resonance. Few scholars have modeled the acoustic resonance excitation process. Radavich et al. [27] briefly discussed the coupling of sound waves and flow in the excitation process and simulated the variations in pressure pulsation during excitation, but there has been no extensive research on the mechanism and excitation process of acoustic resonance. Such knowledge has great significance in regards to establishing an accurate, descriptive model. is paper presents a numerical simulation of the excitation process of acoustic resonance in a closed side branch of natural gas pipeline using the commercial CFD code FLUENT. Due to safety concerns, researchers often use air instead of natural gas as a fluid medium in physical tests. Our findings show that the ingredients of the fluid medium have a pronounced impact on acoustic resonance characteristics. Figure 1 shows a functional block diagram of the feedback mechanism of the system under analysis here. e red dotted frame marks the self-excitation system given by Ding [29], and the black solid frame is the feedback mechanism of acoustic resonance in the closed side branch. e mean flow kinetic energy in the main pipeline provides constant energy for the system. e closed side branch structure is a regulator. e acoustic particles in the pipe are the main vibration body. e pressure fluctuation (acoustic standing wave) in the closed side branch is fed back to the regulator to control the continuity of the fluctuation. e flow-induced acoustic resonance has self-excited vibration characteristics, so a self-excited vibration model was used in this study. Figure 2(a) shows a typical closed side branch of a natural gas pipeline in which acoustic resonance occurs. Considering its actual dimensions, which are usually larger than 0.5 m in diameter and 2 Shock and Vibration 1.5 m in length, the elbow at the end of the side branch can be neglected. e simplified geometry of closed side branch is shown in Figure 2(b). e model consists of two parts: the main pipe and the closed side branch. e length and inner diameter of the main pipe are 6000.0 mm and 575.0 mm, respectively. e branch pipe is located at the center of the main pipe and has a length and inner diameter of 1712.5 mm and 323.8 mm, respectively. e flow medium is methane (the main component of natural gas), which is input steadily from the inlet side and flows continuously to the outlet side. e coordinate system is shown in Figure 2(b), where the xaxis falls along the center line of the main pipe and the y-axis along the center line of the branch pipe. e pressuremonitoring point p1 (0, −2000, 0) is located at the bottom of the center of the side branch.

CFD Turbulence Model Configuration.
Reynolds Average Navier-Stokes (RANS), Detached Eddy Simulation (DES), and Large Eddy Simulation (LES) [30,31] are the most commonly used CFD models. Among them, RANS is most often applied due to its minimal computational cost. However, it is often unable to accurately predict pressure fluctuations in transient flow fields. We found that the pressure fluctuation amplitude of acoustic resonance in a closed side branch, as calculated by RANS, gradually decreases over time and ultimately attenuates to 0. We considered a hybrid model of the subregional computation of RANS and LES, in this study, so DES was a workable solution. DES solves near-wall flow field by RANS and other flow fields by LES, thus avoiding the need for an overly fine grid near the LES wall and striking a balance between numerical calculation cost and accuracy [32]. In near-wall flow fields, a realizable k-ε model can accurately simulate boundary-free shear flows [33]. A realizable k-ε DDES model [34][35][36] was used in this study to investigate acoustic resonance in the closed side branch.

Conservation
Equations. e mass, momentum, and energy conservation equations as well as ideal gas law that the acoustic resonance process complies with can be written as follows: where ρ is the gas density, t is the time, p is the pressure, τ is the viscous stress tensor, E is the sum of internal energy and kinetic energy, T is the temperature, h is the gas enthalpy, J is the diffusion energy, k eff is the effective thermal conductivity, τ eff is the effective viscous dissipation coefficient, R is the universal gas constant, and M is the molecular weight.

Transport
Equations. e flows discussed here are characterized by high Reynolds number turbulence, as their Reynolds numbers range from 2.03 E+07 to 1.42 E+08. Additional transport equations are needed to construct the turbulence closure models. e DDES model adopts the following dual realizable k-ε equations [34]: where Y M � ρk 3/2 /l des , l de s � l rke − f d max(0, l rke − C de s Δ max ) C 1ε � max 0.43, (η/(η + 5)), C 2ε � 1.9, C 3ε � tanh|v/u| , η � Sk/ε, G k denotes the turbulent kinetic energy due to the mean velocity, Pr t � 0.85 in the realizable k-ε model, C des � 0.61, and Δ max is the grid spacing, which in the case of a rectilinear hexahedral cell is the maximum edge length. In the region far away from the wall, r d is far less than 1 and f d is equal to 1, then the governing equation changes to LES; in the region near the wall, f d is equal to 0, and then the governing equation is in the RANS form of the realizable k-e model. In short, Y M thus controls the transformation of RANS to LES.

Grid Partition.
A hexahedral structured grid was selected to improve the accuracy and convergence of the computation. All block grids in the pipeline were made O-shaped to improve the quality of the grid, as shown in Figure 3. e grid was refined in the side branch mouth and parts adjacent to the main pipe as the acoustic vortices are mainly located in these areas when acoustic resonance occurs. e thickness of the first layer near the wall was set to 0.04 mm, and the internal area grid increased inward at a growth rate of 1.2.

Boundary Conditions and Solver.
e inlet boundary condition was considered a velocity inlet. e inlet velocities of all cases were configured with a mean flow and without fluctuation (Table 1)

and represented by Strouhal number S t .
e Strouhal number is a dimensionless number that characterizes the unsteady similarity of fluid flow and is calculated as S t � fd/v, where d is the mouth diameter of the side branch, v is the mean flow velocity of the main pipe, and f is the first-order acoustic resonance frequency calculated as follows [37]: n � 1, 2, 3, . . . , is the order of the frequency, c is the flow field sound velocity, L is the branch length, and Le is correction factor, 0.4-0.425 times the diameter of the branch. e outlet boundary condition was considered a pressure outlet with a pressure of 0 Pa and operating pressure of 6 MPa. e methane gas (CH 4 , the main component of natural gas) in the pipe was assumed to be compressible ideal gas. e wall temperature was set to 293.15 K. e walls were designed as rigid and nonpermeable, so heat does not transfer through them. It should be noted that although there was no heat exchange between the fluid and the environment through the wall in the simulation, energy conversion did exist in the flow (e.g., the conversion of pressure potential energy to heat energy, molecular friction leading to heat generation), so equation (3) is a necessary component of the CFD model. e PISO algorithm coupled with neighbour and skewness corrections was adopted as the conservation equation solver. e conservation equations were spatially discretized by the finite volume method with the discretization schemes shown in Table 2. e parameters of the mean flow velocity, pressure, and flow time were normalized as the dimensionless parameters Strouhal number St, normalized pressure p * , and normalized time t * , respectively. ese parameters are calculated as follows: where f is the first-order frequency of acoustic resonance, d is the diameter of the side branch, v is the mean flow velocity, p is the pressure amplitude, ρ is the gas density, t is the flow time, and l is the distance from the inlet to the branch.

Grid Independent Solution Validation.
It is necessary to verify the grid independence to prevent grid quality from affecting the accuracy of the numerical computation. For grid independence verification, the number of adjacent grids is usually set to be approximately twice as large as that of adjacent grids [38,39]. Four sets of grids were prepared to compute with the boundary condition with a 40 m/s inlet velocity. e cell number of grids was about 159151, 354435, 790423, and 1494845, respectively. All the cases were simulated to 2 s, at which point the quasisteady state of acoustic resonance was reached. e pressure fluctuation at the monitoring point p1 of the four sets of grids is shown in Figure 4(a). e dominant frequency of all pressure fluctuations was equal to 50.6 Hz by FFT analysis, which shows that the four sets of grids did not affect the pressure fluctuation frequency in this case.
As shown in Figure 4(a), the pressure fluctuation waveforms of 159151 grids have larger errors than those of other grids while the pressure fluctuation results of 354435, 790423, and 1494845 grids are basically the same, with nearly coincident curves. e amplitude and phase of the dominant frequency pressure fluctuation of four sets of grids were obtained by FFT as shown in Figure 4

Experimental Validation.
It is difficult to do experiments under pressure as high as 6 MPa. We compared our simulation results of 0.3 MPa with the experimental results of Ziada Shine [11]. As shown in Figure 5, the CFD simulation pressure is generally higher than the experimental pressure, as the wall was rigid and impenetrable in the CFD simulation. e pipe absorbed vibration energy in the experiment, however, so the acoustic energy could propagate outward through the pipe wall. e simulation result is basically consistent with the experimental result at the Strouhal number (St) range of acoustic resonance, which indicates that our model can effectively simulate the acoustic resonance in the closed side branch. Figure 6 shows the time history of pressure fluctuation in monitoring point p1 at different mean flow velocities. All of these cases were initialized from the inlet of the main pipe at the corresponding velocity and were calculated until a stable acoustic resonance formed. e acoustic resonance excitation process initially attenuated at various flow velocities, presenting cluttered signals. e difference was whether acoustic resonance finally formed. For St > 1, as shown in Figures 6(a) and 6(b)), acoustic resonance did not occur and the pressure fluctuation amplitude was low. For St < 1, as shown in Figures 6(c)-6(h)), acoustic resonance did occur and the amplitude of pressure fluctuation was higher. It appears to take some time for stable acoustic resonance to form, less time when the mean flow velocity is relatively large.

Pressure Fluctuation and FFT Analysis.
For main pipe mean flow velocity of v > 40 m/s, the pressure fluctuation shows beat vibration characteristics indicative of unstable self-excited vibration, though continuous fluctuations were maintained in this case.
According to the FFT analysis (           (Figure 7(g)), unlike other cases where acoustic resonance occurred, the frequencies of pressure fluctuation contained two frequencies: a dominant frequency of 50.3 Hz and a secondary frequency of 185.6 Hz. ese results indicate that the mean flow velocity v � 60 m/s fell in the transitional velocity range from the first-order acoustic resonance to the second. e order of the acoustic resonance frequency is defined by the parameter n in equation (7).

Frequency and Amplitude Characteristics of Excitation
Process.
e dominant frequency and its pressure amplitude are presented in Figure 8 as function of St and velocity.
e dominant frequencies in the cases of 0.3 < St < 1 showed a frequency lock-in phenomenon, which has been found in previous experimental studies. In the cases of St > 1, the dominant frequency and its amplitude were lower and the pressure fluctuation amplitude could be neglected. e pressure fluctuation amplitude reached its maximum when St � 0.4, which is consistent previous experimental results wherein the peak pressure fluctuation amplitude appears around St of 0.4 or 0.45 [37]. Previously published experiments have generally shown that acoustic resonance occurs in the range of 0.2 < St < 0.6 when the fluid medium is air. In the present study, acoustic resonance still occurred when St reaches 1.0 for the fluid medium of methane. is indicates that the fluid medium affects the value of St at which acoustic resonance occurs. e fluid medium of methane acoustic resonance occurs at a higher St (lower mean flow velocity) than air, so the natural gas pipeline is a riskier environment in terms of acoustic resonance in closed side branches.    e pressure fluctuations in other cases were analyzed using the same method to find that the excitation processes follow the same mode as the case of v � 50 m/s. As shown in Figures 5(a) and 5(b)), in the case of the mean flow velocity v � 10 m/s and 15 m/s, the pressure fluctuation still had a component of approximately 50 Hz but no acoustic resonance occurred. Comparison with Figure 9(a) indicates that in the case of mean flow velocity v � 10 m/s and 15 m/s, the power of the 50 Hz component was too small to excite the gas acoustic resonance in the closed side branch.  Figure 10, the mass flow rate into the closed side branch was sinusoidal; the amplitude increased gradually with time until a stable acoustic resonance formed. e pressure at the bottom of the branch appears to be positively correlated with the mass flow rate into the branch. e pressure increased when the mass flow rate was positive, and vice versa. Acoustic resonance excitation is essentially a process in which the gas is pressed into and out of the branch in accordance with the sinusoidal mode.

Vortex Modes.
e 200 time steps shown in Figure 10 were evenly divided into eight parts corresponding to the 0T/8-8T/8 of the period. e vorticity on the z � 0 plane of the dividing time points (0T/8-8T/8) is shown in Figures 11-13 for the three time intervals and the velocity curl: e vorticity around the mouth of the branch increased with the flow time (Figures 11-13). e vorticity shows three distinct vortex modes at the three time intervals. At t1 � 0.1-0.12 s, the detached shear layer beat the air at the branch mouth with low amplitudes; the vortex mode at most moments was similar to that of backward-facing step flow. At t2 � 0.6123-0.6323 s, the detached shear layer formed a large vortex at the branch mouth which then impinged on the downstream wall. Small vortices were transmitted to the bottom of the branch at this time in a mode similar to shallow cavity flow. At t3 � 1.7-1.72 s, the detached shear layer formed a larger vortex that was pushed out of the branch mouth without impingement on the downstream wall. is is a unique vortex mode.
We conclude that the acoustic resonance excitation process in the closed side branch can be divided into three steps: (1) the shear layer is free to detach from the upstream wall, (2) the detached shear layer interacts with the downstream wall at the branch mouth, and (3) the detached shear layers are coupled to the gas column resonance to form a periodic large vortex, which is excited in the branch as the shear layer no longer interacts with the downstream wall of the branch.
At the time (6T/8 ⟶ 7T/8 ⟶ 8T/8 ⟶ 0T/8 ⟶ 1T/8), the shear layer gradually detached from the upstream wall to form a large vortex. is large vortex propagated from the main pipe to the branch, pressing the gas into the branch. We define this as the "former vortex." At the time (2T/8 ⟶ 3T/ 8 ⟶ 4T/8 ⟶ 5T/8 ⟶ 6T/8), the large vortex was gradually pushed to the main pipe by the compressed gas in the branch, accompanied by the gas moving out of the branch. e "later detached shear layer" from the upstream wall was brought upwards by the former vortex.
At 4T/8 (Figure 13(e)), the later detached shear layer was almost perpendicular to the flow direction. At 5T/8 ( Figure 13(f )), the later detached shear layer was broken into an independent vortex (the later detached vortex) by the former detached vortex. At 6T/8 ( Figure 13(g)), the vorticity of the later detached vortex began to exceed the former detached vortex. At 7T/8 ( Figure 13(h)), the trajectory of the later vortex and former vortex intersected; the two vortices integrated to a certain extent. In this unique mode, the vortex mode repeats periodically so that a stable acoustic resonance forms. We found that after this formation, the center of the large vortex is located in the upstream half around the branch mouth and the large eddy does not impinge the downstream wall.

Standing Wave.
e pressure appeared to form a standing wave in the branch (Figure 14), so we further examined this aspect of the excitation process. e pressure along the branch axis at three time intervals is shown in Figure 15. At t1, there is no obvious standing wave phenomenon but the wave node is clearly located near the inner side of the branch mouth. At t2, the wave node obviously falls near the inner side of the branch mouth and the antinode is at the bottom of the branch. At t3, the antinode is located at the bottom of the branch and there are no obvious nodes; all are located near to the outer side of the branch mouth. As shown in Figure 15 e reason for the irregularity of the standing wave node is observable in Figure 13. During the positive pressure time, the vortex moved toward the branch pipe, which contributed to the pressure drop; during the negative pressure time, the vortex moved toward the main pipe. e motion of the vortex prevented the sound source from being fixed. e distance from the branch mouth to the main pipe was greater than that to the branch, which is why equation (7) needs to contain the correction coefficient Le.
e frequency validation result f 1c � 50.8 Hz, which is smaller than the result predicted by equation (7), is is because the acoustic vortex (sound source) moved toward the main pipe, so that the wavelength of the standing wave increased and the frequency decreased significantly. At the time intervals of t1 and t2, the standing wave length was shorter due to the effect of the 65 Hz pressure fluctuation component.  [40] pointed out that the energy transfer from the flow field to the sound field can overcome the energy taken away by sound damping and radiation, maintaining high-pressure fluctuations. e mechanism of self-excited resonance in the closed side branch is similar to that around inline cylinders. Figure 16 shows the phase trajectories of all the cases in the phase plane, where it appears that whether acoustic resonance is formed determines the final shape of the phase trajectory. When acoustic resonance can be formed, the limit cycle exists and the phase trajectory eventually converges to the limit cycle. When acoustic resonance cannot form, the limit cycle does not exist and the phase trajectory converges to the equilibrium point. In the case of the mean flow velocity v � 10 m/s and 20 m/s, the phase trajectory tended toward the coordinate origin (0, 0) from the outside to the center, indicating that the equilibrium point fell at the approximate coordinate origin (Figures 16(a) and 16(b)). In the other cases, the phase trajectories converged to a limit cycle from the center to the outside. A curve direction toward the center (limit cycle) reflects the consumption (production) of source energy. Again, whether acoustic resonance occurs depends on whether the sound source energy is greater than the energy taken away by the sound damping and radiation, and the area of this curve illustrates the strength of the source.

Phase Plane Characteristics. Shaaban and Mohany
However, in the cases of convergence toward the limit cycle (Figure 16 velocity v � 60 m/s or 70 m/s, the oval long axis is on the right side. e limit cycle is a regular oval in the case of the mean flow velocity v � 20 m/s and 30 m/s. e limit cycle is relatively regular but oval in the case of the mean flow velocity v � 40 m/s and an irregular oval in the case of the mean flow velocity v � 50 m/s and 60 m/s. e limit cycle is also relatively oval in the case of the mean flow velocity v � 70 m/s. is phase shape ovality was expected due to harmonics, as is also shown in Figure 7. Figure 16 with Figure 7 together show that the shape of the limit cycle is affected by the second-order dominant frequency component. e limit cycle appears to grow more irregular the more obvious the second-order dominant frequency component is. Furthermore, the influence of the second-order dominant frequency intensifies when closer to the mean flow velocity; the limit cycle also grows more irregular. e phase trajectories were similar to the shape of the final limit cycle in all cases except for mean flow velocity v � 60 m/s (Figure 16(g)). In the case of the mean flow velocity v � 60 m/s (Figure 16(g)), the phase trajectory in the phase plane first shows an oval long axis along the direction of Y, and then ultimately the oval long axis falls along the direction of X. ese two different trends illustrate that the excitation process of the oval long axis is relatively unstable. e amount of regularity should indicate whether the coupling is stable or not. e stability of the above process is related to the approximation of the second-order acoustic resonance velocity. e phase trajectory appears to grow less stable as the mean flow velocity grows closer to the second-order acoustic resonance velocity. As shown in Figure 16(g), the phase trajectory was extremely unstable in the case of the mean flow v � 60 m/s. is indicates that the instability factor increased with flow velocity during the transition from firstorder acoustic resonance to second-order acoustic resonance. However, the process recovered a great deal of its stability in the case of mean flow v � 70 m/s.

Self-Excited Vibration Characteristics of Acoustic
Resonance.
e acoustic resonance in the closed side branch presented self-excited vibration (SEV) in our simulation. e SEV and acoustic resonance characteristics discussed above suggest that acoustic resonance is an SEV system in which the pressure fluctuation of the gas is a vibration variable. In an actual system, the gas conveyed in the pipeline and the pipe can also be regarded as an SEV system. We did not build vibration of the pipeline into our simulation, so this is not discussed here. e mathematical model of the SEV system must be a nonlinear differential equation [29]. As an SEV system, the acoustic resonance in the closed side branch shows pure-            tone characteristics which approximate the simple harmonic function, so it is a quasiharmonic SEV system (weakly nonlinear autonomous system). e natural frequency of a weakly nonlinear SEV system is ω 0 , so the equation of motion is as follows: where f(x, _ x) is a nonlinear function that does not contain the linear term of x, and μ is a small constant. A larger μ value makes the phase trajectory tend toward the limit cycle (or the equilibrium point) more quickly. Our simulation results show that the mean flow velocity in the main pipe v affects the speed of the phase trajectory approaching the limit cycle as well as the shape of the limit cycle. We conclude that the mean flow velocity v in the main pipe should be within the nonlinear damping force system of SEV system equation (12) and positively correlated with μ. e length of the closed side branch and the acoustic velocity of the fluid determine the frequency of acoustic resonance (equation (7)), so the length of the branch and the acoustic velocity determine the ω 0 in equation (12).

Conclusion
We simulated acoustic resonance in a single closed side branch of a natural gas pipeline based on a realizable k-ε DDES turbulence model in this study. We investigated flow characteristics of the acoustic resonance excitation process accordingly via the phase plane method from the SEV perspective. Our conclusions can be summarized as follows: (1) In the excitation process, the dominant acoustic resonance frequency of the closed side branch can be selected and amplified while other components gradually disappear. e acoustic resonance frequency ultimately dominates the pure-tone characteristic.  (2) A large vortex forms around the branch mouth which brings the gas in to (out of ) the branch pipe with a sinusoidal periodic mass. e vortex mode is unique from the background step flow and shallow cavity flow, but the excitation process successively experiences the characteristics of both flows. (3) e excitation time and stability of acoustic resonance are related to the mean flow velocity. e excitation time is negatively correlated with the mean flow velocity; it is shorter when the mean flow velocity is relatively large. e stability of acoustic resonance decreases with the mean flow velocity during the transition from the first-order acoustic resonance to the second order, and then the secondorder acoustic resonance recovers the stability. (4) Acoustic resonance is an SEV system. We recommend modeling the acoustic resonance as a weakly nonlinear SEV equation.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.