Inﬂuences of High-Speed Train Speed on Tunnel Aerodynamic Pressures

: To comprehensively investigate the characteristics of aerodynamic pressures on a tunnel caused by the whole tunnel passage of a high-speed train at different speeds, we conduct a series of three-dimensional numerical simulations. Based on the ﬁeld test results obtained by other researchers, the input parameters of our numerical simulation are determined. The process of a high-speed train travelling through a railway tunnel is divided into three stages according to the spatial relationship between the train and tunnel. Stage I: before train nose enters the entrance; Stage II: while the train body runs inside the tunnel; Stage III: after the train tail leaves the exit. The inﬂuences of high-speed train speed on the tunnel aerodynamic pressures of these three stages are systematically investigated. The results show that the maximum peak pressure value decreases with increasing distance from the entrance and increases with increasing train speed in Stage I. There is an approximately linear relationship between the three types of maximum peak pressure (positive peak, negative peak, and peak-to-peak pressures) and the power of the train speed in Stage II. These three types of maximum peak pressure values of the points near tunnel portals increase with increasing train speed in Stage III. Moreover, these three types of maximum peak pressure in the tunnel’s middle section at different train speeds are more complex than those near the tunnel portals, and there is one or more turning points due to the superimposed effects of different pressure waves.


Introduction
After a high-speed train travels through a railway tunnel, the aerodynamic pressure waves are generated and exert loads on the tunnel lining. The interactions between the running train and the air inside the tunnel produce pressure waves that propagate along the tunnel at the speed of sound. Moreover, the pressure field surrounding the running train also exerts pressure on tunnel lining at the same speed of train [1]. The pressure waves are partially reflected at tunnel portals and other discontinuities, the interactions of which produce a complex pressure field on tunnel [2][3][4][5][6]. The adverse effects associated with the aerodynamic pressures increase with increasing train speed, because of the approximately linear relationship between the amplitude of pressure change and the square of train speed [7]. For example, the aerodynamic fatigue damage of tunnel lining under long-term cycles of aerodynamic pressures [8]. Therefore, the research on the influence of train speed on tunnel aerodynamic pressures has a practical significance for tunnel structural safety during operation.
The aerodynamic pressures associated with a high-speed railway (HSR) tunnel have been investigated by field test, moving model rig, theoretical analysis, and numerical simulation. The related studies mainly focus on the train speed, train marshalling form, train nose length, tunnel length, tunnel hood configuration, and blockage ratio. Using the field test, the influences of train speed and marshalling form on tunnel aerodynamic pressures were investigated [9,10]. Adopting the moving model rig, the effects of tunnel hood configuration and train speed on tunnel aerodynamic pressures were studied [11,12]. Employing the theoretical analysis, the influences of hood configuration and Mach number on initial aerodynamic pressure were examined [13,14]. Applying the numerical simulation, the influences of train nose shape and length on tunnel aerodynamic pressures were explored [15,16]. Moreover, the problem of two trains running in opposite directions inside a tunnel was examined by numerical simulation [17] and field test [11]. Although considerable researchers have studied the characteristics of tunnel aerodynamic pressures, they mostly focused on the behaviours during the high-speed train running inside a railway tunnel. The behaviours before the high-speed train enter the tunnel and after it exited the tunnel have not been further explored.
In recent years, the propagation process of initial pressure wave produced by an HSR train through numerical simulation was studied [18]. The results showed that tunnel aerodynamic pressures were generated before the train nose entering the entrance. a series of field tests on the Beijing-Shanghai HSR line were conducted [19]. The results showed that tunnel aerodynamic pressures did not stop immediately after train tail leaving the exit but presented a periodic change. However, these phenomena and their influencing factors have not been substantially studied.
Therefore, the behaviours of the whole tunnel passage of the high-speed train are assessed, i.e., considering not only the train running inside the tunnel, but also the other two series before the train nose entering the entrance and after the train tail leaving the exit. In this research, numerous three-dimensional numerical calculations are performed to explore the characteristics of tunnel aerodynamic pressures. A high-speed train travelling through a 1000 m long railway tunnel with a speed from 275 km/h to 400 km/h is simulated. First, the input parameters in the simulation are verified by field test results. Then, the time history of aerodynamic pressures is classified into three stages according to the spatial relationship between train and tunnel. Finally, the influences of high-speed train speed on tunnel aerodynamic pressures of these three stages are systematically analyzed.

Governing Equations
The flow field, produced by an HSR train travelling through the tunnel, is threedimensional, viscous, compressible, and turbulent [20][21][22][23]. Computational fluid dynamics is commonly used to solve the problems involving fluid flows. The key governing equations in fluid dynamics are continuity equations, momentum equations, and energy equations.
Continuity equations: ∂ρ ∂t where ρ and t are air density and time, respectively. u i is air speed of x i -axis. The subscript i denotes the ith axis of the Cartesian coordinate. Momentum equations: where p is air static pressure. g i is gravitational acceleration of x i -axis. τ ij is the stress tensor, and can be expressed as follows [24]: where µ is molecular viscosity. δ ij is the kronecker symbol. The definitions of the subscripts j and k are similar to i. Appl where T is temperature. k and c p are thermo-conductivity and specific heat capacity of fluid, respectively. S T is a viscous dissipative term [25]. The airflows around the train inside the tunnel are considered compressible, so the state equation of ideal gas is satisfied. That is: where R is gas constant.

Numerical Simulation Method
The software FLUENT 17.2 is adopted in the simulation. The Reynolds-averaged Navier-Stokes (RANS) k-ε turbulence model is adopted to simulate the flow field inside the tunnel [26,27]. The coupled relationship between pressure and velocity is dealt with by the Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) algorithm. The upwind scheme and implicit scheme are set as second-order, which produce better results than the linear model [28,29]. The residual values of continuity, momentum, and energy equations are set as 10 −5 , 10 −5 , and 10 −7 , respectively. The time step is set as 0.004 s. The maximum iteration in each time step is set as 50. The calculations are conducted at the High-Performance Computing Center of Beijing Jiaotong University.

Tunnel and Train Models
A typical double-track HSR tunnel in China is adopted in our study. The configuration of the tunnel is shown in Figure 1. The center-to-center distance of two tracks is 5.0 m. The internal cross-sectional area and the length of the tunnel are 100 m 2 and 1000 m, respectively. A CRH380A vehicle is modelled ( Figure 2). As the accessory components of the HSR train (such as roof pantographs and handlebars) have negligible influences on the simulation results; they are not considered to improve computational efficiency. The constant cross-sectional area of the train body is 11.22 m 2 . The blockage ratio is 0.11. The train speed (v) in our simulation varies from 275 km/h to 400 km/h.

Computational Domain
The computational domain of our simulation includes a tunnel domain and two outer domains, as shown in Figure 3. The two outer domains are identical cuboids. The dimensions of these two outer domains are large enough to ensure that airflows inside the tunnel can be fully developed [30]. The starting location of the HSR train is placed 50 m in front of the tunnel entrance [31,32].

Computational Domain
The computational domain of our simulation includes a tunnel domain and two outer domains, as shown in Figure 3. The two outer domains are identical cuboids. The dimensions of these two outer domains are large enough to ensure that airflows inside the tunnel can be fully developed [30]. The starting location of the HSR train is placed 50 m in front of the tunnel entrance [31,32].

Computational Domain
The computational domain of our simulation includes a tunnel domain and two outer domains, as shown in Figure 3. The two outer domains are identical cuboids. The dimensions of these two outer domains are large enough to ensure that airflows inside the tunnel can be fully developed [30]. The starting location of the HSR train is placed 50 m in front of the tunnel entrance [31,32].

Boundary Conditions
The boundary conditions of our study are shown in Figure 4. The computational domain is divided into stationary part 1 and sliding part 2 to realize the relative movement between train and tunnel. The common surfaces of stationary part 1 and sliding part 2 are set as the interface. The top, two sides, front, and back surfaces of stationary part 1 are set as the pressure far-field. The bottom surface of stationary part 1 is set as the no-slip solid-wall [33]. The surfaces of stationary part 1, nearby the tunnel entrance and exit, are also set as the no-slip solid-wall. The top and two sides surfaces of sliding part 2 are set as the interface. The front and back surfaces are set as the pressure far-field, and the bottom surface is set as the no-slip solid-wall. The surfaces of tunnel wall and train body are set as the no-slip solid-wall. The pressure far-field allows pressure waves to penetrate without reflection. The interface allows the exchange of information between stationary part 1 and sliding part 2. The no-slip solid-wall allows pressure waves to reflect back from the wall. wall [33]. The surfaces of stationary part 1, nearby the tunnel entrance and exit, are also set as the no-slip solid-wall. The top and two sides surfaces of sliding part 2 are set as the interface. The front and back surfaces are set as the pressure far-field, and the bottom surface is set as the no-slip solid-wall. The surfaces of tunnel wall and train body are set as the no-slip solid-wall. The pressure far-field allows pressure waves to penetrate without reflection. The interface allows the exchange of information between stationary part 1 and sliding part 2. The no-slip solid-wall allows pressure waves to reflect back from the wall.

Sliding Mesh Method
The sliding mesh method has been widely adopted to simulate the running train in the tunnel [34]. Compared with the dynamic mesh method, the sliding mesh method updates the mesh without mesh-deformation and mesh-regeneration during the simulation. As a result, the calculation time required for the sliding mesh method is less. Thus, the sliding mesh method is applied in our simulation. The information exchange of airflows inside the computational domain can be accomplished by the interfaces. The schematic diagram of information exchange on the interface of the sliding mesh method is shown in Figure 5. When the sliding grid part 1 moves toward the right, the relative motion between grid part 1 and part 2 forms a common interface. The information on airflows is transmitted from zone 1-I, 1-II, and 1-III to zone 2-I, 2-II, and 2-III, respectively. The number of intersectional grids on the common interface changes when the interface 1 of part 1 moves relative to the interface 2 of part 2. The information of airflows passing through the common interface is calculated by the intersectional grids on the common interface. Thus, the information on airflows can be exchanged between adjacent parts.

Sliding Mesh Method
The sliding mesh method has been widely adopted to simulate the running train in the tunnel [34]. Compared with the dynamic mesh method, the sliding mesh method updates the mesh without mesh-deformation and mesh-regeneration during the simulation. As a result, the calculation time required for the sliding mesh method is less. Thus, the sliding mesh method is applied in our simulation. The information exchange of airflows inside the computational domain can be accomplished by the interfaces. The schematic diagram of information exchange on the interface of the sliding mesh method is shown in Figure 5. When the sliding grid part 1 moves toward the right, the relative motion between grid part 1 and part 2 forms a common interface. The information on airflows is transmitted from zone 1-I, 1-II, and 1-III to zone 2-I, 2-II, and 2-III, respectively. The number of intersectional grids on the common interface changes when the interface 1 of part 1 moves relative to the interface 2 of part 2. The information of airflows passing through the common interface is calculated by the intersectional grids on the common interface. Thus, the information on airflows can be exchanged between adjacent parts.

Grid Generation
The surface meshes of train nose and tunnel portal are shown in Figure 6. The hexahedral mesh is adopted in grid generation. The quality of the near-wall mesh can be evaluated by the dimensionless distance y + . The y + values of train surface and tunnel wall in our simulation are controlled between 30 and 180 for the train speed varies from 275 km/h to 400 km/h, which meet the requirements of the standard [35]. The elements nearby the train streamline nose are refined manually as these regions have significant influences on the simulation results.

Grid Generation
The surface meshes of train nose and tunnel portal are shown in Figure 6. The hexahedral mesh is adopted in grid generation. The quality of the near-wall mesh can be evaluated by the dimensionless distance y + . The y + values of train surface and tunnel wall in our simulation are controlled between 30 and 180 for the train speed varies from 275 km/h to 400 km/h, which meet the requirements of the standard [35]. The elements nearby the train streamline nose are refined manually as these regions have significant influences on the simulation results.

Grid Generation
The surface meshes of train nose and tunnel portal are shown in Figure 6. The hexahedral mesh is adopted in grid generation. The quality of the near-wall mesh can be evaluated by the dimensionless distance y + . The y + values of train surface and tunnel wall in our simulation are controlled between 30 and 180 for the train speed varies from 275 km/h to 400 km/h, which meet the requirements of the standard [35]. The elements nearby the train streamline nose are refined manually as these regions have significant influences on the simulation results.  To study the influences of mesh density on the calculation results, two kinds of mesh densities, coarse mesh (16 million) and fine mesh (20 million), are used. In our simulation, To study the influences of mesh density on the calculation results, two kinds of mesh densities, coarse mesh (16 million) and fine mesh (20 million), are used. In our simulation, the train speed and the tunnel length are 350 km/h and 1000 m, respectively. The comparations of aerodynamic pressure time histories calculated by coarse mesh and fine mesh are shown in Figure 7. The differences between the results obtained by the two mesh densities are negligible. The difference between the maximum positive peaks (maximum negative peaks) of these two mesh densities is only 1.0% (0.2%). Thus, the coarse mesh density is fair enough, and is used in the subsequent numerical simulations.  Figure 7. The differences between the results obtained by the two mesh densities are negligible. The difference between the maximum positive peaks (maximum negative peaks) of these two mesh densities is only 1.0% (0.2%). Thus, the coarse mesh density is fair enough, and is used in the subsequent numerical simulations.

Layout of Monitoring Points
The layouts of monitoring points in the longitudinal section and cross section of the tunnel are shown in Figure 8. A total of 19 monitoring sections are arranged at different

Layout of Monitoring Points
The layouts of monitoring points in the longitudinal section and cross section of the tunnel are shown in Figure 8. A total of 19 monitoring sections are arranged at different locations along the tunnel. There are nine monitoring points on each monitoring section.

Layout of Monitoring Points
The layouts of monitoring points in the longitudinal section and cross section of the tunnel are shown in Figure 8. A total of 19 monitoring sections are arranged at different locations along the tunnel. There are nine monitoring points on each monitoring section.  We use the labels in the form of "1-II" to indicate the locations of monitoring points on the tunnel wall. The Arabic numerals in the labels show the cross-section ID, and the Roman numerals show the ID of the monitoring points in each monitoring section. For example, monitoring point 1-II indicates the second monitoring point on the first monitoring cross section.

Validation
To determine the values of input parameters in our simulation, the numerical results are compared with the field test results. A series of field tests were conducted on the Beijing-Shanghai high-speed railway line [36]. The length of CRH2C train and doubletrack tunnel were 200 m and 1005 m, respectively. The cross-sectional area of tunnel was 100 m 2 . The pressure sensors were installed on tunnel wall 1.5 m above the floor and 140 m from the entrance. The train speed of field tests was 300 km/h. The numerical simulation is performed in the same conditions as those of the field test. The aerodynamic pressures obtained by numerical simulation and field test are shown in Figure 9. The aerodynamic pressures obtained by the numerical simulation agree well with the field test results. The differences in the maximum positive peak, maximum negative peak, and maximum peak-peak values between these two methods are 9.8%, 2.1%, and 6.9%, respectively. As many factors, such as, temperature, humidity, and discontinuities, are not considered in the numerical simulation, there is a small difference between the time histories of aerodynamic pressures obtained by numerical simulation and field test. Generally, the numerical simulation can reasonably predict the field monitoring data. Therefore, the values of input parameters in the numerical simulation are suitable, and they can be used in the subsequent analysis. obtained by numerical simulation and field test are shown in Figure 9. The aerodynamic pressures obtained by the numerical simulation agree well with the field test results. The differences in the maximum positive peak, maximum negative peak, and maximum peakpeak values between these two methods are 9.8%, 2.1%, and 6.9%, respectively. As many factors, such as, temperature, humidity, and discontinuities, are not considered in the numerical simulation, there is a small difference between the time histories of aerodynamic pressures obtained by numerical simulation and field test. Generally, the numerical simulation can reasonably predict the field monitoring data. Therefore, the values of input parameters in the numerical simulation are suitable, and they can be used in the subsequent analysis.

Typical Stages of Time History of Aerodynamic Pressures
To facilitate the illustration, the process of a single high-speed train travelling through a railway tunnel is divided into three stages ( Figure 10). Stage I: before the train nose enters the entrance; Stage II: while the train body runs inside the tunnel; Stage III: after the train tail leaves the exit. Noted, the length of tunnel is longer than that of train [8].

Typical Stages of Time History of Aerodynamic Pressures
To facilitate the illustration, the process of a single high-speed train travelling through a railway tunnel is divided into three stages ( Figure 10). Stage I: before the train nose enters the entrance; Stage II: while the train body runs inside the tunnel; Stage III: after the train tail leaves the exit. Noted, the length of tunnel is longer than that of train [8]. In the case of a single high-speed train travels through a 1000 m long railway tunnel with a speed of 350 km/h (97.22 m/s), the development of aerodynamic pressures with time at some specific monitoring points are shown in Figure 11. It takes 0.5 s from starting for the HSR train nose to arrive at the tunnel entrance. It takes 12.9 s from starting for the train tail to leave the exit. After the train tail leaves the exit, the tunnel aerodynamic pressures last for 67.   In the case of a single high-speed train travels through a 1000 m long railway tunnel with a speed of 350 km/h (97.22 m/s), the development of aerodynamic pressures with time at some specific monitoring points are shown in Figure 11. It takes 0.5 s from starting for the HSR train nose to arrive at the tunnel entrance. It takes 12.9 s from starting for the train tail to leave the exit. After the train tail leaves the exit, the tunnel aerodynamic pressures last for 67. The peak pressures and the duration of Stage I are negligible in comparison to those of Stage II and Stage III. Therefore, the previous studies did not pay attention to the aerodynamic pressures of Stage I. However, considering the entire tunnel life cycle, the total duration of aerodynamic pressures and the total number of peak pressures of Stage I cannot be ignored. For example, a total of 195 pairs of high-speed trains travel on the Guangzhou-Shenzhen HSR line every day in China. The design life of an HSR tunnel is 100 years [8]. The total duration of aerodynamic pressures and the total amount of peak pressures of Stage I are 1384 h and 1.42 × 10 7 , respectively.
The maximum peak pressure of Stage I is much lower than the fracture strength of tunnel lining, and the total number of peak pressures is much higher than 10 6 . Thus, the fatigue damage of the tunnel lining produced by the aerodynamic pressures of Stage I can be analyzed by the method of high-cycle fatigue [37]. To comprehensively investigate the long-term influences of aerodynamic effects on tunnel lining, the aerodynamic pressures of Stage I should be considered, especially near tunnel portals.

Aerodynamic Pressures in Stage I
When the train nose arrives at the entrance at a speed of 350 km/h, the aerodynamic pressure time histories at some specific monitoring points are shown in Figure 12  The peak pressures and the duration of Stage I are negligible in comparison to those of Stage II and Stage III. Therefore, the previous studies did not pay attention to the aerodynamic pressures of Stage I. However, considering the entire tunnel life cycle, the total duration of aerodynamic pressures and the total number of peak pressures of Stage I cannot be ignored. For example, a total of 195 pairs of high-speed trains travel on the Guangzhou-Shenzhen HSR line every day in China. The design life of an HSR tunnel is 100 years [8]. The total duration of aerodynamic pressures and the total amount of peak pressures of Stage I are 1384 h and 1.42 × 10 7 , respectively.
The maximum peak pressure of Stage I is much lower than the fracture strength of tunnel lining, and the total number of peak pressures is much higher than 10 6 . Thus, the fatigue damage of the tunnel lining produced by the aerodynamic pressures of Stage I can be analyzed by the method of high-cycle fatigue [37]. To comprehensively investigate the long-term influences of aerodynamic effects on tunnel lining, the aerodynamic pressures of Stage I should be considered, especially near tunnel portals.

Aerodynamic Pressures in Stage I
When the train nose arrives at the entrance at a speed of 350 km/h, the aerodynamic pressure time histories at some specific monitoring points are shown in Figure 12. It takes 0.18 s (60 m / 340 m/s = 0.18 s) for the initial pressure wave produced by running of train to reach the monitoring point 1-V. When the train nose arrives at the entrance, the propagation distance of the initial pressure wave is 174.86 m (340 m/s × (50 m/97.22 m/s) = 174.86 m). The airflows from the tunnel entrance to the initial wave surface are disturbed by the initial pressure wave. As a result, the aerodynamic pressures are recorded from monitoring section 1 to monitoring section 6. Meanwhile, the airflows from the initial wave surface to the tunnel exit are undisturbed. The distributions of the maximum peak pressures on the tunnel vault along the tunnel at different train speeds are shown in Figure 13. When the train nose arrives at the entrance, the propagation distances of the initial pressure waves at different train speeds are shown in Table 1. The peak pressure values of the monitoring points located at 200 m from tunnel entrance are not available when the train speed is faster than 300 km/h. The characteristics of maximum peak pressures on the tunnel vault along the tunnel's longitudinal axis at different train speeds are similar. The maximum peak pressure decreases with increasing distance from the entrance. The maximum peak pressure increases with increasing train speed.  The distributions of the maximum peak pressures on the tunnel vault along the tunnel at different train speeds are shown in Figure 13. When the train nose arrives at the entrance, the propagation distances of the initial pressure waves at different train speeds are shown in Table 1. The peak pressure values of the monitoring points located at 200 m from tunnel entrance are not available when the train speed is faster than 300 km/h. The characteristics of maximum peak pressures on the tunnel vault along the tunnel's longitudinal axis at different train speeds are similar. The maximum peak pressure decreases with increasing distance from the entrance. The maximum peak pressure increases with increasing train speed. The distributions of the maximum peak pressures on the tunnel vault along the tunnel at different train speeds are shown in Figure 13. When the train nose arrives at the entrance, the propagation distances of the initial pressure waves at different train speeds are shown in Table 1. The peak pressure values of the monitoring points located at 200 m from tunnel entrance are not available when the train speed is faster than 300 km/h. The characteristics of maximum peak pressures on the tunnel vault along the tunnel's longitudinal axis at different train speeds are similar. The maximum peak pressure decreases with increasing distance from the entrance. The maximum peak pressure increases with increasing train speed.   The maximum peak pressures at different train speeds of monitoring point 1-V and the related fitting curve are shown in Figure 14. The 'x' and 'y' show the train speed and the maximum peak pressure, respectively. The maximum peak pressure on tunnel vault increases linearly with increasing train speed, and the value of the adjusted coefficient of determination R 2 is 0.99.

Aerodynamic Pressures in Stage II
Most previous studies in the field of the aerodynamic pressure waves inside tunnel focused on the behaviours in Stage II. The reflection of pressure waves from tunnel portals was widely investigated [38]. However, the reflection and transmission of pressure waves associated with the running train inside the tunnel have not been substantially explored, which have notable influences on the tunnel lining pressure. To cater to this issue, we systematically analyze the development of the pressure waves associated with the train body in this research.
The schematic diagram of propagation and reflection of a pressure wave between tunnel wall and train body is shown in Figure 15. When a compression wave Pa arrives at train tail, the effective sectional area of the airspace inside the tunnel is suddenly contracted due to the train. Then, most of the pressure wave Pa continues to transmit as a compression wave, and the magnitude of transmitted wave Pta is larger than Pa. The pressure wave Pa is partially reflected back as a compression wave, and the magnitude of the reflected wave Pra is smaller than Pa. When the compression wave Pta arrives at train nose, the effective sectional area of the airspace is suddenly enlarged, and the magnitude of transmitted wave Ptta is decreased. The reflected wave Prta turn into a rarefaction wave. Noted, the reflected wave is inverted (turning upside down) if the incident pressure wave disturbs a more massive medium, and the reflection is upright if the incident wave disturbs a less massive medium. The transmission of a pressure wave from the train nose to tail is similar to that from tail to nose. In fact, a pressure wave is partially reflected and partially transmitted at the boundary between media (i.e., tunnel portals, trains inside the tunnel). Therefore, numerous pressure waves are generated and propagate inside the tunnel, of which the superposition results in a complex pressure field on tunnel lining. Because of the friction effects of tunnel wall and train body, the intensity of pressure waves gradually decreases and ultimately disappears.

Aerodynamic Pressures in Stage II
Most previous studies in the field of the aerodynamic pressure waves inside tunnel focused on the behaviours in Stage II. The reflection of pressure waves from tunnel portals was widely investigated [38]. However, the reflection and transmission of pressure waves associated with the running train inside the tunnel have not been substantially explored, which have notable influences on the tunnel lining pressure. To cater to this issue, we systematically analyze the development of the pressure waves associated with the train body in this research.
The schematic diagram of propagation and reflection of a pressure wave between tunnel wall and train body is shown in Figure 15. When a compression wave P a arrives at train tail, the effective sectional area of the airspace inside the tunnel is suddenly contracted due to the train. Then, most of the pressure wave P a continues to transmit as a compression wave, and the magnitude of transmitted wave P ta is larger than P a . The pressure wave P a is partially reflected back as a compression wave, and the magnitude of the reflected wave P ra is smaller than P a . When the compression wave P ta arrives at train nose, the effective sectional area of the airspace is suddenly enlarged, and the magnitude of transmitted wave P tta is decreased. The reflected wave P rta turn into a rarefaction wave. Noted, the reflected wave is inverted (turning upside down) if the incident pressure wave disturbs a more massive medium, and the reflection is upright if the incident wave disturbs a less massive medium. The transmission of a pressure wave from the train nose to tail is similar to that from tail to nose. In fact, a pressure wave is partially reflected and partially transmitted at the boundary between media (i.e., tunnel portals, trains inside the tunnel). Therefore, numerous pressure waves are generated and propagate inside the tunnel, of which the superposition results in a complex pressure field on tunnel lining. Because of the friction effects of tunnel wall and train body, the intensity of pressure waves gradually decreases and ultimately disappears.
tail is similar to that from tail to nose. In fact, a pressure wave is partially reflected and partially transmitted at the boundary between media (i.e., tunnel portals, trains inside the tunnel). Therefore, numerous pressure waves are generated and propagate inside the tunnel, of which the superposition results in a complex pressure field on tunnel lining. Because of the friction effects of tunnel wall and train body, the intensity of pressure waves gradually decreases and ultimately disappears.  (Note: P a : pressure wave, P ra : reflected wave of P a , P ta : transmitted wave of P a , P rta : reflected wave of P ta , P tta : transmitted wave of P ta ).
The relationship between the wave diagram and the pressure variations of monitoring point 10-V in Stage II is shown in Figure 16. The train speed is 350 km/h. In Figure 16a, the solid lines (with letters C and c) and dotted lines (with letters R and r) denote compression and rarefaction waves, respectively. The letters N and T denote the running paths of train nose and tail inside the tunnel, respectively. The paths C1 and R1 denote the initial pressure waves produced by the entering of train nose and tail, respectively. Except for C1 and R1, the uppercase letters (C and R) denote the reflected waves from tunnel portals. The lowercase letters (c and r) denote the reflected waves from the train nose and tail. The Arabic numerals denote the sequence of the reflected waves from a boundary (portal or train body) in the form of compression or rarefaction waves. As the intensity of the reflected wave from the train body is smaller than that of the incident wave, the subsequent reflections of this reflected wave are ignored. As a result, multiple pressure waves move between the tunnel portals during a train travelling in the tunnel. The recorded aerodynamic pressures of monitoring point 10-V, 500 m from both tunnel portals, are shown in Figure 16b. The vertical lines represent the moments at which the train nose, train tail, and pressure waves arriving at the monitoring section.
When the compression waves and the train tail arrive at the monitoring section, positive leaps of the monitoring pressure are produced. When the rarefaction waves and the train nose arrive at the monitoring section, negative drops are generated. The aerodynamic pressure time history of a point on tunnel lining is a superposition of numerous pressure waves [39]. The aerodynamic pressure field of the tunnel lining is complex, which varies with location and time in Stage II. When the compression waves and the train tail arrive at the monitoring section, positive leaps of the monitoring pressure are produced. When the rarefaction waves and the train nose arrive at the monitoring section, negative drops are generated. The aerodynamic pressure time history of a point on tunnel lining is a superposition of numerous pressure waves [39]. The aerodynamic pressure field of the tunnel lining is complex, which varies with location and time in Stage II.
The distributions of three types of maximum peak pressures (positive peak, negative peak, and peak-to-peak pressures) along the tunnel vault in Stage II are shown in Figure  17. The three types of maximum peak pressure values of the points near tunnel middle section are generally larger than those of the points near tunnel portals. The relationships between the maximum peak pressure value of tunnel vault and the train speed of three The distributions of three types of maximum peak pressures (positive peak, negative peak, and peak-to-peak pressures) along the tunnel vault in Stage II are shown in Figure 17. The three types of maximum peak pressure values of the points near tunnel middle section are generally larger than those of the points near tunnel portals. The relationships between the maximum peak pressure value of tunnel vault and the train speed of three typical tunnel cross section in Stage II are shown in Figure 18. The maximum peak pressure value increases with increasing train speed. There is an approximately linear relationship between each of the three types of maximum peak pressure value and the train speed raised to the power from 1.12 to 2.33 except for the boundary condition 0 kPa at 0 km/h. When the compression waves and the train tail arrive at the monitoring section, positive leaps of the monitoring pressure are produced. When the rarefaction waves and the train nose arrive at the monitoring section, negative drops are generated. The aerodynamic pressure time history of a point on tunnel lining is a superposition of numerous pressure waves [39]. The aerodynamic pressure field of the tunnel lining is complex, which varies with location and time in Stage II.
The distributions of three types of maximum peak pressures (positive peak, negative peak, and peak-to-peak pressures) along the tunnel vault in Stage II are shown in Figure  17. The three types of maximum peak pressure values of the points near tunnel middle section are generally larger than those of the points near tunnel portals. The relationships between the maximum peak pressure value of tunnel vault and the train speed of three typical tunnel cross section in Stage II are shown in Figure 18. The maximum peak pressure value increases with increasing train speed. There is an approximately linear relationship between each of the three types of maximum peak pressure value and the train speed raised to the power from 1.12 to 2.33 except for the boundary condition 0 kPa at 0 km/h.

Aerodynamic Pressures in Stage III
The wave diagram and the aerodynamic pressure time history of monitoring point 10-V in Stage III are shown in Figure 19. The train speed is 350 km/h. After train tail leaves the exit, the pressure field of tunnel lining is produced by the reflections of pressure waves from tunnel portals. The solid and dotted lines denote compression and rarefaction waves, respectively. The red and magenta lines denote pressure waves produced by the entering of train nose and tail, respectively. The blue and green lines denote pressure waves produced by the leaving of train nose and tail, respectively. As the energy of the pressure waves reflected from the train nose and tail is relatively small and has been reduced by the friction effect in Stage II, these reflected waves are ignored in Stage III.

Aerodynamic Pressures in Stage III
The wave diagram and the aerodynamic pressure time history of monitoring point 10-V in Stage III are shown in Figure 19. The train speed is 350 km/h. After train tail leaves the exit, the pressure field of tunnel lining is produced by the reflections of pressure waves from tunnel portals. The solid and dotted lines denote compression and rarefaction waves, respectively. The red and magenta lines denote pressure waves produced by the entering of train nose and tail, respectively. The blue and green lines denote pressure waves produced by the leaving of train nose and tail, respectively. As the energy of the pressure waves reflected from the train nose and tail is relatively small and has been reduced by the friction effect in Stage II, these reflected waves are ignored in Stage III.
The aerodynamic pressure time history presents a periodic variation due to the periodic reflections of the pressure waves from tunnel portals. The period of the aerodynamic pressures is mainly determined by tunnel length and wave speed. The pressure wave period can be defined as T ac = 2 L/c (L: tunnel length, c: sonic speed). The energy of pressure waves is gradually reduced by the dissipations at tunnel portals and the frictions on boundary surfaces (e.g., tunnel wall and track). Therefore, the amplitude of the aerodynamic pressures decreases with time. The amplitude of the aerodynamic pressures is increased after two or more pressure waves, with the disturbances along the same line, arriving at the monitoring section. Moreover, the amplitude of the aerodynamic pressures is decreased after two pressure waves, with the disturbances along the different lines, arriving at the monitoring section. The aerodynamic pressure time history presents a periodic variation due to the periodic reflections of the pressure waves from tunnel portals. The period of the aerodynamic pressures is mainly determined by tunnel length and wave speed. The pressure wave period can be defined as Tac = 2 L/c (L: tunnel length, c: sonic speed). The energy of pressure waves is gradually reduced by the dissipations at tunnel portals and the frictions on boundary surfaces (e.g., tunnel wall and track). Therefore, the amplitude of the aerodynamic pressures decreases with time. The amplitude of the aerodynamic pressures is increased after two or more pressure waves, with the disturbances along the same line, arriving at the monitoring section. Moreover, the amplitude of the aerodynamic pressures is decreased after two pressure waves, with the disturbances along the different lines, arriving at the monitoring section.
The distributions of the maximum peak pressures along the tunnel vault are shown in Figure 20. The maximum peak pressure values of the points near tunnel middle section are generally larger than those near tunnel portals.  Figure 20. The maximum peak pressure values of the points near tunnel middle section are generally larger than those near tunnel portals.
The relationships between the maximum peak pressure and the train speed of specific monitoring points in Stage III are shown in Figure 21. The three types of maximum peak pressure values of the points near tunnel portals increase with increasing train speed. The three types of maximum peak pressure of tunnel middle section are more complex than those near tunnel portals. There are one or more turning points of these three types of maximum peak pressure of tunnel middle section. The relationships between the maximum peak pressure and the train speed of specific monitoring points in Stage III are shown in Figure 21. The three types of maximum peak pressure values of the points near tunnel portals increase with increasing train speed. The three types of maximum peak pressure of tunnel middle section are more complex than those near tunnel portals. There are one or more turning points of these three types of maximum peak pressure of tunnel middle section. The relationships between the maximum peak pressure and the train speed of specific monitoring points in Stage III are shown in Figure 21. The three types of maximum peak pressure values of the points near tunnel portals increase with increasing train speed. The three types of maximum peak pressure of tunnel middle section are more complex than those near tunnel portals. There are one or more turning points of these three types of maximum peak pressure of tunnel middle section. As shown in Figure 19, when the train speed is 350 km/h, the aerodynamic pressure of the point near tunnel middle section achieves the maximum positive peak at 14.1 s. After three consecutive rarefaction waves pass through the monitoring section, the aerodynamic pressure achieves the maximum negative peak at 16.9 s. When the train speed is 375 km/h, the aerodynamic pressure of the point near tunnel middle section achieves the maximum positive peak at 14.2 s (Figure 22). After two rarefaction waves and one compression wave passe through the monitoring section, the aerodynamic pressure achieves the maximum negative peak at 17.0 s. As a result of the superimposed effect of pressure waves, the maximum negative peak pressure value of the point near tunnel middle section decreases after train speed increases from 350 km/h to 375 km/h. The existence of other turning points can also be explained by the superimposed effect of pressure waves. As shown in Figure 19, when the train speed is 350 km/h, the aerodynamic pressure of the point near tunnel middle section achieves the maximum positive peak at 14.1 s. After three consecutive rarefaction waves pass through the monitoring section, the aerodynamic pressure achieves the maximum negative peak at 16.9 s. When the train speed is 375 km/h, the aerodynamic pressure of the point near tunnel middle section achieves the maximum positive peak at 14.2 s (Figure 22). After two rarefaction waves and one compression wave passe through the monitoring section, the aerodynamic pressure achieves the maximum negative peak at 17.0 s. As a result of the superimposed effect of pressure waves, the maximum negative peak pressure value of the point near tunnel middle section decreases after train speed increases from 350 km/h to 375 km/h. The existence of other turning points can also be explained by the superimposed effect of pressure waves.

Conclusions
To study the influences of high-speed train speed on tunnel aerodynamic pressures, we carry out a series of three-dimensional numerical simulations. The input parameters of our simulations are verified by field monitoring data. The process of a single high-speed train travelling through a railway tunnel is divided into three stages according to the spatial relationship between train and tunnel. The influences of high-speed train speed on tunnel aerodynamic pressures of these three stages are systematically studied. The main results are as follows: (1) The peak pressure values at different train speeds are much lower than the fracture strength of tunnel lining, and the total number of peak pressures during the entire life cycle of tunnel in Stage I is much higher than 10 6 . Thus, the aerodynamic fatigue damage of tunnel lining of Stage I can be calculated by the method of high-cycle fatigue.
(2) The characteristics of maximum peak pressure values acting on the tunnel vault along the tunnel's longitudinal axis at different train speeds are similar in Stage I. The maximum peak pressure value decreases with increasing distance from the entrance and increases with increasing train speed.
(3) Three types of maximum peak pressure values (positive peak, negative peak, and peak-to-peak pressures) increase with increasing train speed in Stage II. There is an approximately linear relationship between each of these three types of maximum peak pressure values and high-speed train speed raised to the power from 1.12 to 2.33.
(4) The three types of maximum peak pressure values of the points near tunnel portals increase with increasing train speed in Stage III. These three types of maximum peak pressure of the tunnel middle section at different train speeds are more complex than those near the tunnel portals, and there is one or more turning point of the curves due to the superimposed effects of different pressure waves.