Abstract

When the tunnel has damages such as deformation and cracking, the lining structure should be set to reinforce the tunnel. In this study, a three-dimensional numerical method is performed to study the aerodynamic effects caused by the train passing through a tunnel with segmented linings at 350 km/h. The influence of the numbers, thickness, location, and types of the segmented lining structure on MPW is analyzed. The results show that the segmented linings located at the middle of the tunnel have a better mitigation effect on MPW than that of the normal linings, which increases with the segment number and the thickness of the linings. When the damage occurs near the entrance of the tunnel, a segmented lining structure with a constricted section slightly away from the tunnel entrance can be used to reduce its adverse effect on MPW.

1. Introduction

The impact of aerodynamic effects caused by high-speed trains entering tunnels increases with the speed-up of trains. Compression waves produced generated when the high-speed train enters a tunnel propagate along the tunnel near to the sonic speed and create a pulse wave at tunnel exits, usually called micro-pressure wave (MPW) [1]. The MPW causes problems such as noise pollution and even a sonic boom, which can lead to structural damage to the buildings surrounding the tunnel exit and negatively affect residents’ normal lives [24]. Wide research has shown that assessing MPW and structural vibrations are associated with safety and discomfort caused in trains [5, 6]. Zou et al. developed an efficient computational model that can accurately predict over-track vibration levels for buildings supported on a transfer structure [7].

The train speed and blockage ratio are significant factors that affect the amplitude of MPW [813]. When the tunnel is under the influence of geology and construction, there are tunnel damages such as structural deformation and cracking [14]. Tunnel damage threatens traffic safety in tunnels and shortens the tunnels’ maintenance cycle and service life. Adding linings technology is a popular way to prevent the occurrence or further expansion of cracks. However, adding linings inside the tunnel changes the clearance area and cross-sectional distribution of the tunnel, which affects the aerodynamic characteristics of the tunnel. The lining is to pour a certain thickness of reinforced concrete in the tunnel to prevent the tunnel from further cracking, damage, and deformation. Jin et al. studied the feature morphology of rubberized concrete by corrosion and expansion of reinforcement in rubberized concrete and the relationship between the deformation of staggered assembled shield tunnel and the disengaging of invert filling [15, 16].

Previously many researches were concentrated on the responses of the tunnel hood to MPW of the tunnel. Tunnel hood is an auxiliary structure specifically designed to reduce the impact of tunnel aerodynamics. Xiang et al. compared different tunnel entrance hoods and optimized hood designs to mitigate the pressure gradient [17]. Murray's theoretical research on the tunnel entrance hoods show that the installation of the tunnel hoods will increase the rise time of the initial compression wave, thereby decreasing the pressure gradient [18]. Studies nowadays have contributed extensively to designing optimized tunnel hoods to reduce the amplitude of MPW at the tunnel exit [1923]. Although related studies have been performed on the influences of tunnel hoods, limited investigations have been conducted on the effects of MPW behaviors of the linings. Based on the real vehicle experiments, dynamic model experiments, and numerical simulations method to analyze the aerodynamic effects of the train passing the tunnel with linings, Liu et al. pointed out that the influence of MPW is the most unfavorable when the linings are added at the tunnel entrance and the pressure amplitude increases with the thickness of the linings [24, 25]. Combined with the test results of the real vehicle, Gao found MPW and the comfort of the tunnel when the train passed through the lined tunnel and proposed that the impact section of the linings should not exceed 95% [26]. Shi et al. investigated the influence of linings length, thickness, and position on the pressure gradient at the tunnel exit through 3D numerical simulation and found that the impact of linings thickness was more significant than that of length [27]. Wang et al. analyzed the transient pressure characteristics of trains passing through tunnels with noncircular linings based on dynamic model experiments. The results show that reasonable distribution of circular linings can reduce the pressure amplitude of initial compression waves [28].

The above research has shown that adding a lining structure at the tunnel’s entrance will adversely affect MPW, but the damage may occur anywhere in the tunnel. Therefore, the behaviors of the aerodynamic effect of the segmented lining structure are assessed in this study. In Section 2, the governing equations and mesh motion equations of the train moving in the tunnel are presented. In Section 3, the numerical calculation model of the train tunnel was established and validated with the dynamic model experiments. In Section 4, the influence of the linings’ position, thickness, number of segments, and types on MPW characteristics of the tunnel exit was discussed. The results can provide a fundamental comprehension and reference to the design of the lining structure.

2. Numerical Calculation Method

2.1. Governing Equation

The flow field, produced by trains traveling through the tunnel, is 3D, viscous compressible, and turbulent [9, 2931], which must follow the laws of continuity equations, momentum equations, and energy equations [32]. They can be written aswhere is the density of the air; and are the velocity components in the three directions of the flow field; and are the three components of the coordinate; denotes the pressure; is the shear stress in the flow field; is the thermal conductivity coefficient; is the temperature of the flow field; and is their dynamic viscosity of the air.

The term is the time-averaged Reynolds stresses which represent the turbulent momentum fluxes. Based on the Boussinesq hypothesis, the closed equations under the Boussinesq hypothesis:where is the turbulent viscosity and is the turbulent kinetic energy.

The dynamic mesh method is usually accustomed to simulate the time-varying flow on domain boundaries due to train motion in tunnels. The motion can be a prescribed motion or an unprescribed motion which is determined based on the solution at the current time.

Generally, the integral form of the conservation equation for a general scalar, ϕ, on an arbitrary control volume, V, whose boundaries are moving, can be written as [31, 33]

For solving (3), the differential term is discretized as follows:where and represent the value at the present and next time step. The time step control volume, , is computed from

These equations can then be solved using dynamic meshing methods and updated with new volumes at each time step.

2.2. Computational Model

In this study, the CR400AF vehicle depicted in Figure 1(a) is selected to research the aerodynamic impact of the tunnels with segmented linings. Since the initial compression wave and MPW are mainly associated with the shape of the train nose, four-car marshaling carriages are modeled to improve computational efficiency. The shape of the train is simplified according to the recommendations of the CEN European Standard [34], as the accessory components of the high-speed train (such as pantograph, bogie, doors, and windows) have negligible influence on MPW. According to Chinese high-speed railway tunnel design standards, a double-track tunnel with a cross-sectional area of 100 m2 is adopted for a high-speed train with a design speed of 350 km/h, and the centre-to-centre distance of the two tracks is 5 m. The cross section size of the tunnel is shown in Figure 1(b), and the length of the tunnel is 1000 m. Different numbers of segmented linings with various thicknesses are set within a certain length of 100 m at different locations of the tunnel to investigate the influences of the segmented linings on MPW, as shown in Figure 1(c).

2.3. Computational Domain and Boundary Conditions

The computational domain of the aerodynamic impacts generated by the trains traveling from the tunnel with segmented linings is composed of the atmosphere outside the tunnel portal and the atmosphere inside the tunnel, shown in Figure 2. Considering the acceleration process of the train and the full development of the wake [35], the length of the flow field outside the tunnel is taken as 500 m on both sides of the tunnel.

The train motion in the tunnel is simulated by the sliding mesh model and mesh dynamic layering method, so the entire computational domain is divided into four domains [30, 36]: a motion domain (region-2), two deformation domains (region-1 and region-3), and a fixed domain (region-4). As the train moves at speed as the surrounding motion domain, the lengths of deformation domains are changed to keep the entire computational domain unchanged during computation. Interface boundary conditions are used between the motion domain, the deformation domain, and the fixed domain to realize data exchange. The boundaries of the flow field outside the tunnel portal are defined as pressure inlet and pressure outlet, respectively. No-slip wall is defined as the train surface, ground, and tunnel wall. All walls are defined as Adiabatic wall with an initial temperature of 300 K.

2.4. Grid Generation Scheme

A hybrid grid scheme is adopted to divide the entire computational domain into blocks accounting for computational accuracy and speed [37]. The unstructured grid is used in the motion domain surrounding the train as the complicated shape of the train, and the geometrically complex parts of the train surface are locally refined. The deformation domains before and after the motion domain can be divided into high-quality structured grids, and the adjacent regions are connected by setting common nodes on the interface with the motion domain. The structured grid is also used in the fixed domain, and the regions near the segmented linings are refined to capture the influence of the local flow field on the MPW.

In order to evaluate the mesh density on the numerical results, four different grid numbers were taken into consideration, which are 3.98 × 106, 5.19 × 106, 6.60 × 106, and 8.12 × 106, respectively. The details of grids and the amplitude values of MPW calculated by 4 different numbers of grids are shown in Table 1. The results show that the amplitude values of MPW are very little when the number of the grids exceeds 5.19 × 106, which means that the computational accuracy and efficiency can be met when the mentioned mesh is used. In such cases, the grid of the train head and segmented linings are shown in Figure 3.

2.5. Measuring Points

The evaluation standard of MPW at the tunnel exit of a high-speed railway is that the amplitude values of MPW at 20 m and 50 m away from the tunnel exit are less than 50 Pa and 20 Pa, respectively [22].

Studies show that the magnitude of MPW emitting from the tunnel exit is almost proportional to pressure gradients of the wavefront of the initial compression wave reaching the tunnel exit and inversely proportional to the distance to the tunnel exit plane [38], which can be expressed as follows:where is the cross-sectional area of the tunnel exit, denotes the solid angle around the tunnel exit portal, is the speed of sound, is the distance from the measuring point to the tunnel entrance, and demotes the pressure gradient at the tunnel exit.

Both the pressure gradient near the tunnel exits and MPW outside the tunnel are monitored simultaneously in the calculation. Two groups of monitoring points 1.2 m above the ground are arranged symmetrically on both sides of the tunnel centreline, as shown in Figure 4. The monitoring point located at 50 m inside the tunnel is marked as T1 and T2, and M1,2 and M3,4 denote the monitoring points located at 20 m and 50 m outside the tunnel, respectively.

3. Numerical Method Validation

To verify the accuracy of the presented method in the simulation, the numerical results are compared with that of the moving model test conducted by Central South University [39]. The model ratio in the test is 1 : 20, and the original lengths of the train and the double tracks tunnel are 51.7 m and 1000 m, respectively. A ring tunnel portal with a slope value of 1 : 1.25 is set at tunnel entrance. Both the transient pressure at 65.5 m from the entrance in the tunnel and MPW at 20 m from the exit outside the tunnel are measured in the test when the train passes the tunnel at 350 km/h. The numerical simulation is performed in the same condition as the moving model test.

ANSYS Fluent is used to explain the governing equation of flow motion. The ideal gas model is adopted in the calculation, and the RNG k-ε turbulence model is used due to the high turbidity. The coupling velocity and pressure are dealt with the SIMPLE algorithm, the gradient term is discretized by the Green-Gaussian Cell-Based method. The pressure term is discretized using the standard format, the other spatial terms are discretized by the second-order upwind format and the second-order implicit format is used in the time discretization.

The initial compression wave and MPW acquired by the numerical simulation and moving model tests are shown in Figure 5. It can be seen that the results calculated by the numerical method agree well with those of the moving model test. The amplitude values of MPW and the pressure gradient of the initial compression wave are also listed in Table 2. Because the numerical calculation model has assumptions and errors in the calculation process, the results of the numerical method are not exactly the same as those of the moving model test. The differences between the results from the numerical method and that of the moving model test are 2.86% and 2.60%, respectively. The results show that the numerical methods adopted can meet the requirements of calculation accuracy and efficiency. Thus, the mentioned method in this paper can be used to study the characteristics of MPW caused by the train passing through the tunnel.

4. Results and Discussion

4.1. Propagation Characteristics of Initial Compression Wave

The aerodynamics caused by the train traveling through the tunnel without linings at 350 km/h is simulated by the mentioned method. The comparison of the pressure and pressure gradient of two symmetrical measuring points in the same plane of the tunnel is shown in Figure 6.

It can be found that the results at T1 are basically the same as that at T2, which indicates that the initial compression wave presents a one-dimensional plane wave characteristic propagating in the tunnel [40]. Due to the relationship with the pressure gradient of the initial compression wave at the tunnel exit, the MPW outside the tunnel exit is also symmetrical, as shown in Figure 7. So, the results at the monitoring points on only one side (T1, M1, and M3) are considered in the subsequent discussion.

4.2. Influence of Linings Position and Quantity

The compress wave and MPW characteristics caused by a train passing through tunnels with normal linings and various numbers of segmented linings at different positions with a speed of 350 km/h are studied by numerical simulation.

For the normal linings located at the tunnel entrance, middle, and exit, the maximum values of pressure gradient at T1 and MPW at M1 and M3 are shown in Table 3. It can be found that both the pressure gradient and MPW increase compared to that without linings, which has an adverse effect on the aerodynamics of the tunnel. When the linings are located at tunnel middle and exit, the pressure gradient and MPW decrease compared to that without linings and a favorable effect is induced. The results also show that the mitigation effect on MPW is better when the linings are located near the tunnel exit.

The MPW at M1 and M3 when the trains travel through the tunnel with different numbers of segmented linings at different positions are shown in Figures 8 and 9, respectively.

It can be seen that the influence of the lining’s position on MPW is basically the same for both normal linings and the different number of segmented linings. That is, the linings close to the tunnel middle and exit have a favorable effect on the mitigation of MPW, while an adverse effect will be induced when the linings locate near the entrance of the tunnel.

The amplitude values of MPW at 20 m and 50 m away from the tunnel exit are listed in Tables 4 and 5, respectively. It can be found that when the linings are located at the tunnel entrance and middle, the MPW generated by the tunnel with segmented linings decrease compared to that with normal linings. For the segmented linings locate at the middle and exit of the tunnel, the mitigation effects on the micropressure wave increase with the segment numbers due to the increscent dissipation effect caused by the linings. The amplitude values of MPW at M1 and M3 are 82.172 Pa and 35.096 Pa when a 30-segment lining is set at the tunnel exit, which is 4.789% and 5.277% lower than that of the tunnel without linings.

However, when the segmented linings are located at the tunnel entrance, the amplitude value of MPW increases gradually with the number of segments. The reason is that the shrinking section of the lining gets closer to the tunnel entrance as the increase of the segmented number and the compression effect of the air is dominant compared to the dissipation effect, which has an adverse impact on MPW. Therefore, the adverse effect of the linings located at the tunnel entrance can be reduced by a small number of segments. When a 10-segment lining is set at the entrance of the tunnel, the amplitude values of MPW at M1 and M3 are 91.225 Pa and 39.181 Pa, which are 5.701% and 5.749% higher than that of the tunnel without linings but decrease by 6.884% and 6.699% compared to that of the tunnel with normal linings.

4.3. Influence of Linings Thickness

The influence of the thickness of the linings on MPW is also studied. The aerodynamics caused by a train passing through the tunnel with normal and 10-segment linings of various thicknesses located at the middle of the tunnel is simulated by numerical methods.

The variations of MPW at 20 m and 50 m away from the tunnel exit are shown in Figures 10 and 11, respectively. The results show that MPW decrease slightly with the increment of linings thickness in a certain range for both normal linings and segmented linings due to the increscent dissipation effect.

The amplitude values of MPW at M1 and M3 are also given in Tables 6 and 7. It can be found that the amplitude values of MPW at two monitoring points drop by 0.629% and 0.707% when the thickness of the segmented linings changes from 0.4 m to 0.6 m.

Thus, the thickness of the lining structure can be designed according to the requirements of structural strength in engineering practice and a thick lining can be used as the cost allows to alleviate the MPW problem outside the tunnel exit.

4.4. Influence of Linings Type

The previous studies show that an effect is induced by the segmented linings located at the tunnel entrance and exit. The aerodynamics of the segmented linings with different types located at the tunnel middle and exit are simulated to study the influences of the linings type on MPW.

As shown in Figure 12, Type 1 denotes the segmented linings start with a constant section compared to the tunnel, while Type 2 denotes the segmented linings start with a constricted section, and 10 segments are set for both types with a thickness of 0.5 m.

When a train travels through the tunnels with different types of segmented linings, the MPW at 20 m and 50 m away from the tunnel portal are given in Figures 13 and 14, respectively.

It has shown quite a difference between the two types of segmented linings located at the entrance of the tunnel. For the segmented linings that start with a constant section (Type 1), the rising of the initial compression wave is divided into two stages, resulting in decreases in pressure gradient and MPW. While the segmented linings start with a constricted section (Type 2), the compression effect is enhanced due to the reduction of the clearance area of the tunnel entrance, which has an adverse on the micropressure wave. The amplitude values of MPW at the M1 and M3 are 95.787 Pa and 41.06 Pa, respectively, which are 4.999% and 4.796% higher than that with Type 1.

However, the results show little difference between the two types of segmented linings located at the exit of the tunnel. For the segmented linings end with constricted section (Type 1), the mitigation effect is better due to a smaller clearance area at the tunnel exit, as shown in equation (6). The amplitude values of MPW at the M1 and M3 are 83.204 Pa and 35.569 Pa, which are 0.375% and 0.829% smaller than that with Type 2.

The results show that the type of segmented linings located at tunnel entrance has a greater impact on MPW. Therefore, a segmented lining structure can be used if damages occur close to the tunnel entrance, and the constricted section of the lining should be slightly away from the tunnel entrance.

5. Conclusions

In this paper, a 3D numerical method is used to research the aerodynamic impacts produced by trains traveling through a tunnel with linings at 350 km/h. The influence of the numbers, thickness, location, and types of the segmented lining structure on MPW are analyzed. The following can be concluded from the results:(1)Both normal linings and segmented linings located at the entrance of the tunnel have an adverse effect on MPW, while a favorable effect is induced when the linings locate at the middle and exit of the tunnel and the mitigation effect is better for the linings at the tunnel exit.(2)When the segmented linings are located at the tunnel entrance, the amplitude of MPW increases with the increase of the number of segments. When the segmented linings are located at the middle and exit of the tunnel, the amplitude of MPW decreases with the increase of the number of segments. In addition, when the linings are located in the tunnel middle, the segmented linings have a better mitigation effect on MPW than the normal linings, which increases with the segment number of the linings. When a 30-segment lining structure with a thickness of 0.5 m is set near the exit of the tunnel, the amplitude values of MPW at 20 m and 50 m outside the tunnel decrease by 4.789% and 5.277%, respectively;(3)The mitigation effects of both normal linings and segmented linings increase slightly with the thickness of the linings.(4)When the damage occurs near the entrance of the tunnel, a segmented lining structure can be used to reduce its adverse effect on MPW and the constricted section of the lining structure should be slightly away from the tunnel entrance.(5)The problem of tunnel damage can be repaired by adding segmented linings and reducing MPW. Other aerodynamic effects caused by train tunnel (such as train wind, transient pressure, passenger comfort, and so on) will be studied in future work.

Data Availability

All data generated or analyzed during this study are included in this article and are available upon request by contact with the corresponding author.

Conflicts of Interest

The author declares that there are no conflicts of interest.

Acknowledgments

This work was supported by the State Key Laboratory for Track Technology of High-Speed Railway (Grant no. 2021YJ056) and Natural Science Foundation of Shandong Province, China (Grant nos. ZR2020MA057 and ZR2021QE108), and Science and Technology R&D Program of China Railway (P2021G053).