Analysis of Aerodynamic Noise Characteristics of High-Speed Train Pantograph with Di ff erent Installation Bases

The high-speed-train pantograph is a complex structure that consists of different rod-shaped and rectangular surfaces. Flow phenomena around the pantograph are complicated and can cause a large proportion of aerodynamic noise, which is one of the main aerodynamic noise sources of a high-speed train. Therefore, better understanding of aerodynamic noise characteristics is needed. In this study, the large eddy simulation (LES) coupled with the acoustic finite element method (FEM) is applied to analyze aerodynamic noise characteristics of a high-speed train with a pantograph installed on different configurations of the roof base, i.e. flush and sunken surfaces. Numerical results are presented in terms of acoustic pressure spectra and distributions of aerodynamic noise in near-field and far-field regions under upand down-pantograph as well as flushed and sunken pantograph base conditions. The results show that the pantograph with the sunken base configuration provides better aerodynamic noise performances when compared to that with the flush base configuration. The noise induced by the down-pantograph is higher than that by the up-pantograph under the same condition under the pantograph shape and opening direction selected in this paper. The results also indicate that, in general, the directivity of the noise induced by the down-pantograph with sunken base configuration is slighter than that with the flush configuration. However, for the up-pantograph, the directivity is close to each other in Y-Z or X-Z plane whether it is under flush or sunken roof base condition. However, the sunken installation is still conducive to the noise environment on both sides of the track.


Introduction
Over the last twenty years, high-speed trains have played a major role in community and urban development. It is well known that, as the speed of a high-speed train increases, aerodynamic problems that can be neglected at low speeds become more serious. For such cases, the problem of aerodynamic noise cannot be avoided and must be addressed. Therefore, better understanding of the aerodynamic noise characteristics is needed. To alleviate the aerodynamic noise problems, the shape of complex parts is optimized and appropriate sound barriers are commonly used. Presently, a large number of high-speed trains can reach an operating speed of about 300-350 km/h. With these speeds, the

Computational Methods
The commercial CFD program STAR-CCM+, which is based on the finite-volume method, is used to simulate the flow field around the pantograph and the installation base. The LES method is used to conduct the dynamic flow characteristics, including the fluctuating pressure on the pantograph and the base. Then, the fluctuating pressure in the time domain is transformed into the

Computational Models
The pantograph used in the present study consists of insulators, chassis, upper arm, lower arm rod and double slide, as depicted in Figure 1. The upper arm and lower arm rod can move vertically, and are controlled by the base frame mechanism. On the roof of the train body, the pantograph is installed, either flush or sunken. Namely, the flush surface configuration of the roof is used as the baseline for a comparison with the sunken surface configuration. Scenarios of the up and down pantograph with the different configurations of the pantograph installation on the surface roof are shown in Figure 2.

Computational Methods
The commercial CFD program STAR-CCM+, which is based on the finite-volume method, is used to simulate the flow field around the pantograph and the installation base. The LES method is used to conduct the dynamic flow characteristics, including the fluctuating pressure on the pantograph and the base. Then, the fluctuating pressure in the time domain is transformed into the

Computational Methods
The commercial CFD program STAR-CCM+, which is based on the finite-volume method, is used to simulate the flow field around the pantograph and the installation base. The LES method is used to conduct the dynamic flow characteristics, including the fluctuating pressure on the pantograph and the base. Then, the fluctuating pressure in the time domain is transformed into the frequency domain by the Fast Fourier Transform (FFT) method and the acoustic FEM with commercial program LMS Virtual lab is used to evaluate the acoustic noise of radiation from the pantograph.

Large Eddy Simulation
According to LES, the large-scale eddies in a fully turbulent flow are computed directly and the influence of the small eddies on the large-scale eddies is modeled. The continuity equation and the Favre averaged compressible Navier-Stokes equations, as expressed mathematically in Equations (1) and (2), respectively, are solved implicitly.
where σ i,j and τ i,j are defined as Equations (3) and (4), respectively.
where u i is the resolved filtered velocity vector, σ i,j is the stress tensor caused by molecular viscosity, and τ i,j is the sub-grid scale (SGS) stress tensor, representing the diffusive effect of the sub-grid scale eddies on the resolved ones. The Smagorinsky model [20] is used to model the SGS stresses because of its simplicity.

Acoustic Finite Element Method
The acoustic FEM is used to compute the sound pressure level (SPL). This method can take the reflection effect of the car body structure surface on the noise [21] into consideration. The acoustic FEM procedure starts from calculating the noise propagation from the Hemholtz equation as defined in Equation (5): p(x, y, z) + k 2 p(x, y, z) = f(x, y, z) where p(x, y, z) is the acoustic pressure, k = 2π f /c is the wave-number, and f is the frequency. The corresponding wavelength is computed from Equation (6): where c is the free-stream speed of sound. Then, the FFT method is adopted for spectral analysis. The square of the amplitude of pressure wave is expressed in terms of the summation of sine and cosine functions, as shown in Equation (7): where P n is a data set and represents the fluctuating pressure of the Nth step; n = 0, . . . , N − 1.
Meanwhile, ω k = k/Nδt, where δt is the time step. Finally, the SPL is computed by Equation (8).

Computational Domain and Boundary Conditions
The computation domain is depicted in Figure 3. The flow field is tested to confirm that it is large enough to eliminate boundary effects. The computational domain size is 96.5H × 45H × 20H, where H is the height of the train and H = 3.55 m. The origin of the Cartesian reference system locates at the distance x = 29 m from the train geometry. The inflow velocity is 350 km/h. The Reynolds number, Re, is 2.25 × 10 7 based on inflow velocity and the train's height H. A three-car train model with the pantograph is used to study the effects of the pantograph with different configurations of the installation base on the aerodynamic noise characteristics. The effects of other complicated components, such as bogies and windshields, are neglected to save computational cost.

Computational Domain and Boundary Conditions
The computation domain is depicted in Figure 3. The flow field is tested to confirm that it is large enough to eliminate boundary effects. The computational domain size is 96.5 45 20 where H is the height of the train and H = 3.55 m. The origin of the Cartesian reference system locates at the distance x = 29 m from the train geometry. The inflow velocity is 350 km/h. The Reynolds number, Re, is 2.25  10 7 based on inflow velocity and the train's height H. A three-car train model with the pantograph is used to study the effects of the pantograph with different configurations of the installation base on the aerodynamic noise characteristics. The effects of other complicated components, such as bogies and windshields, are neglected to save computational cost. The ground is set to be a moving wall, and its velocity is the same as the incoming velocity of the air. Figure 4 shows the observer locations used to evaluate the acoustic pressure spectrum and directivity in the near and far field of the pantograph, namely, two positions at y = 5 and 25 m for the spectrum characteristics, and planes y = 0, 10 and 20 m for the distributions of aerodynamic noise, respectively. Furthermore, five positions in the x-direction i.e. x = 0, 5, 10, 15 and 20 m are monitored to investigate the overall sound pressure level (OASPL) caused by the pantograph with the different base configurations. It should be noted that these five positions are 3.5 m above the ground. The physical time step of the calculation is 10 −4 s and 20 sub-steps are used within each time step. The total physical calculation time is 0.4 s. The Courant-Friedrichs-Lewy number (CFL = u ∆t/∆ x) is less than 1 so that the calculation is converged within each time step. Boundary conditions are set as follows: (1) The surface of the train is a fixed wall with non-slip and penetration conditions. (2) Both lateral sides and the top of the computational domain are given as the far field condition with the characteristic line method, and the air flows in the positive x-direction. (3) The ground is set to be a moving wall, and its velocity is the same as the incoming velocity of the air. Figure 4 shows the observer locations used to evaluate the acoustic pressure spectrum and directivity in the near and far field of the pantograph, namely, two positions at y = 5 and 25 m for the spectrum characteristics, and planes y = 0, 10 and 20 m for the distributions of aerodynamic noise, respectively. Furthermore, five positions in the x-direction i.e. x = 0, 5, 10, 15 and 20 m are monitored to investigate the overall sound pressure level (OASPL) caused by the pantograph with the different base configurations. It should be noted that these five positions are 3.5 m above the ground.

Computational Domain and Boundary Conditions
The computation domain is depicted in Figure 3. The flow field is tested to confirm that it is large enough to eliminate boundary effects. The computational domain size is 96.5 45 20 where H is the height of the train and H = 3.55 m. The origin of the Cartesian reference system locates at the distance x = 29 m from the train geometry. The inflow velocity is 350 km/h. The Reynolds number, Re, is 2.25  10 7 based on inflow velocity and the train's height H. A three-car train model with the pantograph is used to study the effects of the pantograph with different configurations of the installation base on the aerodynamic noise characteristics. The effects of other complicated components, such as bogies and windshields, are neglected to save computational cost. x) t/ u = (CFL   is less than 1 so that the calculation is converged within each time step. Boundary conditions are set as follows: (1) The surface of the train is a fixed wall with non-slip and penetration conditions. (2) Both lateral sides and the top of the computational domain are given as the far field condition with the characteristic line method, and the air flows in the positive x-direction. (3) The ground is set to be a moving wall, and its velocity is the same as the incoming velocity of the air. Figure 4 shows the observer locations used to evaluate the acoustic pressure spectrum and directivity in the near and far field of the pantograph, namely, two positions at y = 5 and 25 m for the spectrum characteristics, and planes y = 0, 10 and 20 m for the distributions of aerodynamic noise, respectively. Furthermore, five positions in the x-direction i.e. x = 0, 5, 10, 15 and 20 m are monitored to investigate the overall sound pressure level (OASPL) caused by the pantograph with the different base configurations. It should be noted that these five positions are 3.5 m above the ground.

Mesh Strategy
The STAR-CCM+ mesher is used to generate the computational grids. On the grid arrangement of flow field computation domain, the first layer thickness is 0.05 mm, the geometric grid growth rate is 1.1, and the number of the grid layers close to the train surface is 10 to resolve flow in the boundary layer. The grid size on the pantograph surface and train body surface is 10 mm and 60 mm, respectively, and mesh refinement rear the train body and pantograph is employed to make the computational grid fine enough to capture flow field physics. An isotropic mesh is used in the whole flow field except near wall boundaries. Three sets of grids, i.e., coarse, medium and fine grids, are used to test the grid sensitivity of the simulations. The total number of grid cells in the three grid configurations are 15, 30 and 60 million cells, respectively. Four mesh refinement zones including small, middle, big and wake region are set up, as shown in Figure 5. The mesh size of each refinement zone is listed in Table 1.

Mesh Strategy
The STAR-CCM+ mesher is used to generate the computational grids. On the grid arrangement of flow field computation domain, the first layer thickness is 0.05 mm, the geometric grid growth rate is 1.1, and the number of the grid layers close to the train surface is 10 to resolve flow in the boundary layer. The grid size on the pantograph surface and train body surface is 10 mm and 60 mm, respectively, and mesh refinement rear the train body and pantograph is employed to make the computational grid fine enough to capture flow field physics. An isotropic mesh is used in the whole flow field except near wall boundaries. Three sets of grids, i.e., coarse, medium and fine grids, are used to test the grid sensitivity of the simulations. The total number of grid cells in the three grid configurations are 15, 30 and 60 million cells, respectively. Four mesh refinement zones including small, middle, big and wake region are set up, as shown in Figure 5. The mesh size of each refinement zone is listed in Table 1.

Grids
Coarse Medium Fine Grid size settings (mm) (small, middle, big, wake) 18, 18, 36, 9 8, 12, 24, 6 6, 9, 18, 6 Total numbers 15 million 30 million 60 million Figure 6 shows the distribution of time-averaged slipstream velocity along a sampling line at a distance of 2 m from the center of the train, 1.3 m above the ground predicted by the coarse, medium and fine meshes. The non-dimensional results of slipstream velocity non U is defined as: where Vx, Vy and Vz are the velocity components along the x-direction, y-direction and z-direction, respectively, as shown in Figure 3.   Figure 6 shows the distribution of time-averaged slipstream velocity along a sampling line at a distance of 2 m from the center of the train, 1.3 m above the ground predicted by the coarse, medium and fine meshes. The non-dimensional results of slipstream velocity U non is defined as: where V x , V y and V z are the velocity components along the x-direction, y-direction and z-direction, respectively, as shown in Figure 3.

Mesh Strategy
The STAR-CCM+ mesher is used to generate the computational grids. On the grid arrangement of flow field computation domain, the first layer thickness is 0.05 mm, the geometric grid growth rate is 1.1, and the number of the grid layers close to the train surface is 10 to resolve flow in the boundary layer. The grid size on the pantograph surface and train body surface is 10 mm and 60 mm, respectively, and mesh refinement rear the train body and pantograph is employed to make the computational grid fine enough to capture flow field physics. An isotropic mesh is used in the whole flow field except near wall boundaries. Three sets of grids, i.e., coarse, medium and fine grids, are used to test the grid sensitivity of the simulations. The total number of grid cells in the three grid configurations are 15, 30 and 60 million cells, respectively. Four mesh refinement zones including small, middle, big and wake region are set up, as shown in Figure 5. The mesh size of each refinement zone is listed in Table 1.

Grids
Coarse Medium Fine Grid size settings (mm) (small, middle, big, wake) 18, 18, 36, 9 8, 12, 24, 6 6, 9, 18, 6 Total numbers 15 million 30 million 60 million Figure 6 shows the distribution of time-averaged slipstream velocity along a sampling line at a distance of 2 m from the center of the train, 1.3 m above the ground predicted by the coarse, medium and fine meshes. The non-dimensional results of slipstream velocity non U is defined as: where Vx, Vy and Vz are the velocity components along the x-direction, y-direction and z-direction, respectively, as shown in Figure 3.  Figure 6 shows that the time-averaged non-dimensional slipstream velocity shows larger difference between the coarse and the fine meshes near the train in Region 2 and in the wake propagation Region 3, indicating that the coarse mesh is inadequate analytical accuracy, while the medium mesh exhibits good consistence with the fine mesh. Thus, the medium grid configuration is adopted to calculate characteristics of flow field, as it gives a satisfactory balance between accuracy and computational costs. Figure 7 shows the medium grids around the train surfaces. Because pantograph noise is the main research content, the area near the pantograph is refined. The total number of the computational cells used is about 32 million.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 8 of 21 Figure 6. The distribution non-dimensional time-averaged slipstream velocity. Figure 6 shows that the time-averaged non-dimensional slipstream velocity shows larger difference between the coarse and the fine meshes near the train in Region 2 and in the wake propagation Region 3, indicating that the coarse mesh is inadequate analytical accuracy, while the medium mesh exhibits good consistence with the fine mesh. Thus, the medium grid configuration is adopted to calculate characteristics of flow field, as it gives a satisfactory balance between accuracy and computational costs. Figure 7 shows the medium grids around the train surfaces. Because pantograph noise is the main research content, the area near the pantograph is refined. The total number of the computational cells used is about 32 million.

Validation of Computational Methods
Many structures of the pantograph could be treated as cylinders. The acoustic finite element method is first validated by a cylindrical flow calculation, as shown in Figure 8. The cylinder has a diameter of 10 mm, the length of the cylinder is D, and the radius of the flow field is 15D. The trimmer meshing method is employed, and the first layer thickness is 0.0035 mm with a growth ratio of 1.1 in the wall-normal direction. The value of normalized wall-normal distance y+ is defined as: where y is the distance from the wall,  u is the friction velocity (  is the air density and is the air kinetic viscosity. The values of y+ are below 1 in downstream of the front facing stagnation line on the cylinder surface, which is adequate for the LES model. The cylindrical surface grid size is 0.25 mm. The total number of cells is about 6 million. Figure 8 shows the boundary conditions, namely, the pressure far-field boundary with a uniform inlet velocity of 72 m/s and a gauge pressure of 0 Pa, and the parallel plane (on the two ends of a cylinder) is set as periodic.

Validation of Computational Methods
Many structures of the pantograph could be treated as cylinders. The acoustic finite element method is first validated by a cylindrical flow calculation, as shown in Figure 8. The cylinder has a diameter of 10 mm, the length of the cylinder is πD, and the radius of the flow field is 15D. The trimmer meshing method is employed, and the first layer thickness is 0.0035 mm with a growth ratio of 1.1 in the wall-normal direction. The value of normalized wall-normal distance y+ is defined as: where y is the distance from the wall, u * is the friction velocity (u * = τ w ρ 1/2 ), τ w is wall shear stress, ρ is the air density and ν is the air kinetic viscosity. The values of y+ are below 1 in downstream of the front facing stagnation line on the cylinder surface, which is adequate for the LES model. The cylindrical surface grid size is 0.25 mm. The total number of cells is about 6 million. Figure 8 shows the boundary conditions, namely, the pressure far-field boundary with a uniform inlet velocity of 72 m/s and a gauge pressure of 0 Pa, and the parallel plane (on the two ends of a cylinder) is set as periodic.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 8 of 21 Figure 6. The distribution non-dimensional time-averaged slipstream velocity. Figure 6 shows that the time-averaged non-dimensional slipstream velocity shows larger difference between the coarse and the fine meshes near the train in Region 2 and in the wake propagation Region 3, indicating that the coarse mesh is inadequate analytical accuracy, while the medium mesh exhibits good consistence with the fine mesh. Thus, the medium grid configuration is adopted to calculate characteristics of flow field, as it gives a satisfactory balance between accuracy and computational costs. Figure 7 shows the medium grids around the train surfaces. Because pantograph noise is the main research content, the area near the pantograph is refined. The total number of the computational cells used is about 32 million.

Validation of Computational Methods
Many structures of the pantograph could be treated as cylinders. The acoustic finite element method is first validated by a cylindrical flow calculation, as shown in Figure 8. The cylinder has a diameter of 10 mm, the length of the cylinder is D, and the radius of the flow field is 15D. The trimmer meshing method is employed, and the first layer thickness is 0.0035 mm with a growth ratio of 1.1 in the wall-normal direction. The value of normalized wall-normal distance y+ is defined as: where y is the distance from the wall,  u is the friction velocity (  is the air density and is the air kinetic viscosity. The values of y+ are below 1 in downstream of the front facing stagnation line on the cylinder surface, which is adequate for the LES model. The cylindrical surface grid size is 0.25 mm. The total number of cells is about 6 million. Figure 8 shows the boundary conditions, namely, the pressure far-field boundary with a uniform inlet velocity of 72 m/s and a gauge pressure of 0 Pa, and the parallel plane (on the two ends of a cylinder) is set as periodic.  The LES method is used to calculate the flow field. The pressure-velocity coupling method is addressed by the Semi-Implicit Method for Pressure Linked Equations (SIMPLE) algorithm. The discrete algebraic equation is solved by the Gauss-Seidel iteration technique. Figure 9 shows the generation, development, and detachment of the leeward vortices downstream of the cylinder. The vortices alternate axially in the downstream direction of the flow, forming a von Karman vortex street. It can be seen that the leeward side vortices are mainly generated when the airflow flowing from the cylindrical wall surface leaves the leeward side wall surface. The generation, development and shedding of the vortices lead to noise generation. The acoustic FEM is used to solve the far-field noise at an inflow Mach number is 0.21. In this study, Kato equation [22] is used to correct the SPL obtained from numerical simulations because the spanwise length of the LES computational domain is limited by the available computational resources. As a result, the simulation spanwise length, L sim , of the geometry is shorter than the experimental spanwise length, L exp , reported in Jacob et al. [23]. To compare the predicted SPL with experimental data, the value obtained from numerical simulations must be corrected. The SPL is monitored 185D above directly the cylinder. Figure 10 shows the validations of the aerodynamic results and SPL results of 185D. It can be seen that the pressure coefficient (Cp) and the corresponding Strouhal number (St) of the vortex shedding predicted by the present study agrees with the experimental results and other numerical simulation solutions obtained in [23][24][25][26]. Therefore, it is reasonable to use the FEM for the far-field noise characteristic prediction in the following analysis.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 9 of 21 The LES method is used to calculate the flow field. The pressure-velocity coupling method is addressed by the Semi-Implicit Method for Pressure Linked Equations (SIMPLE) algorithm. The discrete algebraic equation is solved by the Gauss-Seidel iteration technique. Figure 9 shows the generation, development, and detachment of the leeward vortices downstream of the cylinder. The vortices alternate axially in the downstream direction of the flow, forming a von Karman vortex street. It can be seen that the leeward side vortices are mainly generated when the airflow flowing from the cylindrical wall surface leaves the leeward side wall surface. The generation, development and shedding of the vortices lead to noise generation. The acoustic FEM is used to solve the far-field noise at an inflow Mach number is 0.21. In this study, Kato equation [22] is used to correct the SPL obtained from numerical simulations because the spanwise length of the LES computational domain is limited by the available computational resources. As a result, the simulation spanwise length, Lsim, of the geometry is shorter than the experimental spanwise length, Lexp , reported in Jacob et al. [23]. To compare the predicted SPL with experimental data, the value obtained from numerical simulations must be corrected. The SPL is monitored 185D above directly the cylinder. Figure 10 shows the validations of the aerodynamic results and SPL results of 185D. It can be seen that the pressure coefficient (Cp) and the corresponding Strouhal number (St) of the vortex shedding predicted by the present study agrees with the experimental results and other numerical simulation solutions obtained in [23][24][25][26]. Therefore, it is reasonable to use the FEM for the far-field noise characteristic prediction in the following analysis.  This paper seeks to show that the numerical method is able to model the evolution of the structures in the gap between the solid bodies. A comparison between downstream flow characteristics of the cylinder obtained by measurements and simulation using the same mesh level changes in the cylinder wake as used for the pantograph simulation is performed. As shown in Figure  11, the cylinder has a diameter of 0.05715 m, the separation distance, L, between the cylinders is 3.7D. The length is 3D. The trimmer meshing method is employed, the first layer thickness is 0.001 mm The LES method is used to calculate the flow field. The pressure-velocity coupling method is addressed by the Semi-Implicit Method for Pressure Linked Equations (SIMPLE) algorithm. The discrete algebraic equation is solved by the Gauss-Seidel iteration technique. Figure 9 shows the generation, development, and detachment of the leeward vortices downstream of the cylinder. The vortices alternate axially in the downstream direction of the flow, forming a von Karman vortex street. It can be seen that the leeward side vortices are mainly generated when the airflow flowing from the cylindrical wall surface leaves the leeward side wall surface. The generation, development and shedding of the vortices lead to noise generation. The acoustic FEM is used to solve the far-field noise at an inflow Mach number is 0.21. In this study, Kato equation [22] is used to correct the SPL obtained from numerical simulations because the spanwise length of the LES computational domain is limited by the available computational resources. As a result, the simulation spanwise length, Lsim, of the geometry is shorter than the experimental spanwise length, Lexp , reported in Jacob et al. [23]. To compare the predicted SPL with experimental data, the value obtained from numerical simulations must be corrected. The SPL is monitored 185D above directly the cylinder. Figure 10 shows the validations of the aerodynamic results and SPL results of 185D. It can be seen that the pressure coefficient (Cp) and the corresponding Strouhal number (St) of the vortex shedding predicted by the present study agrees with the experimental results and other numerical simulation solutions obtained in [23][24][25][26]. Therefore, it is reasonable to use the FEM for the far-field noise characteristic prediction in the following analysis.  This paper seeks to show that the numerical method is able to model the evolution of the structures in the gap between the solid bodies. A comparison between downstream flow characteristics of the cylinder obtained by measurements and simulation using the same mesh level changes in the cylinder wake as used for the pantograph simulation is performed. As shown in Figure  11, the cylinder has a diameter of 0.05715 m, the separation distance, L, between the cylinders is 3.7D. The length is 3D. The trimmer meshing method is employed, the first layer thickness is 0.001 mm This paper seeks to show that the numerical method is able to model the evolution of the structures in the gap between the solid bodies. A comparison between downstream flow characteristics of the cylinder obtained by measurements and simulation using the same mesh level changes in the cylinder wake as used for the pantograph simulation is performed. As shown in Figure 11, the cylinder has a diameter of 0.05715 m, the separation distance, L, between the cylinders is 3.7D. The length is 3D.
The trimmer meshing method is employed, the first layer thickness is 0.001 mm with a growth ratio of 1.2 in the wall-normal direction. The cylindrical surface grid size is 1 mm. The total number of cells is about 8.5 million. The pressure outlet is set as a gauge pressure of 0 Pa, the inlet condition is set as a uniform velocity of 44 m/s, and the other parallel planes are given as symmetry.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 10 of 21 with a growth ratio of 1.2 in the wall-normal direction. The cylindrical surface grid size is 1 mm. The total number of cells is about 8.5 million. The pressure outlet is set as a gauge pressure of 0 Pa, the inlet condition is set as a uniform velocity of 44 m/s, and the other parallel planes are given as symmetry. Figure 11. Geometry and computational domain for flow field of two cylinders in tandem. Figure 12a shows the generation, development and detachment of the leeward vortices downstream of the two cylinders. The predicted root-mean-square (rms) of the pressure coefficient (Cp) on the cylinder downstream cylinder surfaces is compared with the previous work [27], in Figure 12b. The angle "Theta" is measured from the upstream stagnation point and is positive in the clockwise direction. The predicted results show that this method is suitable to obtain the flow characteristics in the gap and downstream between the solid bodies.   Figure 12a shows the generation, development and detachment of the leeward vortices downstream of the two cylinders. The predicted root-mean-square (rms) of the pressure coefficient (Cp) on the cylinder downstream cylinder surfaces is compared with the previous work [27], in Figure 12b. The angle "Theta" is measured from the upstream stagnation point and is positive in the clockwise direction. The predicted results show that this method is suitable to obtain the flow characteristics in the gap and downstream between the solid bodies.

Flow Characteristics of Pantograph Area
Appl. Sci. 2019, 9, x FOR PEER REVIEW 10 of 21 with a growth ratio of 1.2 in the wall-normal direction. The cylindrical surface grid size is 1 mm. The total number of cells is about 8.5 million. The pressure outlet is set as a gauge pressure of 0 Pa, the inlet condition is set as a uniform velocity of 44 m/s, and the other parallel planes are given as symmetry. Figure 11. Geometry and computational domain for flow field of two cylinders in tandem. Figure 12a shows the generation, development and detachment of the leeward vortices downstream of the two cylinders. The predicted root-mean-square (rms) of the pressure coefficient (Cp) on the cylinder downstream cylinder surfaces is compared with the previous work [27], in Figure 12b. The angle "Theta" is measured from the upstream stagnation point and is positive in the clockwise direction. The predicted results show that this method is suitable to obtain the flow characteristics in the gap and downstream between the solid bodies.

Flow Characteristics of Pantograph Area
In this section, the flow characteristics around the pantograph and the installation bases are discussed. As shown in Figures 13 and 14, the vorticity and the instantaneous iso-surface vorticity as shown by the normalized Q-criterion, and the velocity streamlines around the pantograph are presented, respectively. In Figure 13, the flow separates over the pantograph surface. A series of vortices detaches from the leeward side of each rod and they interact with each other downstream. The hairpin vortices of different size are observed. In addition to the common characteristics of the flow field mentioned above, Figure 13a,b also shows that the flow field characteristics of the pantograph are different due to the opening directions of the pantograph. For the down-pantograph with flush installation base, the incoming flow first impacts on the corner of the upper and lower arm rods, and then strikes the rear of the carbon skateboard and the base in the downstream area. The turbulent upstream flow field hits downstream complex structure like the carbon skateboard and base again, which makes the flow field in the wake region more complex. Because the contact area between the pantograph and the incoming flow is concentrated, the eddy size under the larger contact area is also larger. For the up-pantograph, the distance between different structures on the pantograph is larger, and the inflow flows directly to the downstream area after impact with the carbon skateboard. Compared with the down-pantograph, the rear eddy scale is smaller. For Figure 13c, when the pantograph is placed in the cavity, the flow field above the cavity shielding part is almost the same as when it is installed on the flush surface. However, for the case that the pantograph is installed in the sunken base, the flow field characteristics have something in common with the general flow field in a cavity, as previously studied in [28][29][30]. However, the presence of the pantograph affects the pressure fluctuation on the cavity surface, and then it further affects the noise induced by the pantograph. For a pantograph that is installed in the sunken cavity, the base and insulator are basically located in the cavity, avoiding direct collision with the incoming flow. Generally, the flow past a cavity exhibits a strong relationship with the cavity length (L) to height (D) ratio as studied by Lawson and Barakos [19]. For the present study, the ratio of L/D of the cavity with pantograph is about 4.7. Therefore, according to the criterion of Lawson and Barakos [19], it is considered as the open cavity. In addition, the pressure field in the cavity is coupled with the shear layer shed from the cavity upstream edge, and the pressure fluctuation intensity in the cavity depends on the characteristics of the shear layer and on the evolution of momentum and vorticity during the impact of the shed shear layer with the cavity rear bulkhead as studied by Ouyang et al. [31]. Because the presence of the pantograph, the typical open cavity flow field structure is destroyed. As shown in Figure 14, one can observe that the air-flow acts on the pantograph first, and then vortices shed by the pantograph impinge on the surface downstream. A part of the air flow returns to the front of the cavity and continues to act on the pantograph. This differs from an empty cavity flow in which the shear layer shed at the cavity upstream edge impacts the cavity downstream edge, unobstructed. The pressure wave (feedback pressure wave) propagating upstream after acting on the surface of the cavity collides with the pressure wave reflected from the front wall when it arrives upstream at the last time. The air circulation inside and outside the cavity exchanges, resulting in pressure fluctuation in the cavity. The pressure fluctuation on the base and insulator surface even as well as the overall pressure fluctuation on the surface of the pantograph is affected by the flow in the cavity.

Aerodynamic Noise Source
In this section, sources of aerodynamic noise are discussed as it is expected that the flow field of the different base configurations of the pantograph may have significant effects on noise spectrum characteristics and on the noise directivity, which could lead to a change in the perceived noise along the train track. Therefore, the aerodynamic pressure amplitude distributions on the pantograph at 200 Hz and 1000 Hz are compared. As shown in Figure 15, it can be observed that along the direction of the flow velocity, the aerodynamic pressure amplitude in the front regions of the pantograph is lower than that in the rear regions of the pantograph. This may be explained by the fact that the front parts always collide with the incoming flow sharply, forming a windward stagnation area. As a result, the fluctuating pressure is relatively small. Flow with higher levels of vorticity, as shown in Figure 13, passes over the rear parts of the pantograph. As a result, the amplitude of aerodynamic pressure fluctuation in the rear region is larger. Meanwhile, the amplitude of the aerodynamic pressure decreases when the frequency increases from 200 Hz to 1000 Hz. Where the flush configuration is used under the down-pantograph condition, the maximum aerodynamic pressure fluctuation at 1000 Hz is smaller than that at 200 Hz by as much as 18 dB. The phenomenon that the maximum aerodynamic pressure fluctuation decreases with the increasing frequency also can be observed by Zhang et al. [15]. At the low frequency, the SPL distribution of the aerodynamic pressure on the surface of the pantograph is non-uniform. The largest amplitude aerodynamic pressure mainly features on the supporting slide, the rotating shaft, the lower arm, and the leeward surface of the installation base. At 1000 Hz, the aerodynamic pressure amplitude distribution is more uniform than that at 200 Hz. It seems that, although the aerodynamic pressure amplitude distribution over the pantograph surfaces are different at these two frequencies, the largest amplitudes are always located on the downstream base frame, the connection between the upper and lower arms and the downstream position of the double skateboard bow. These aerodynamic pressure fluctuations act as dipole noise sources.

Aerodynamic Noise Source
In this section, sources of aerodynamic noise are discussed as it is expected that the flow field of the different base configurations of the pantograph may have significant effects on noise spectrum characteristics and on the noise directivity, which could lead to a change in the perceived noise along the train track. Therefore, the aerodynamic pressure amplitude distributions on the pantograph at 200 Hz and 1000 Hz are compared. As shown in Figure 15, it can be observed that along the direction of the flow velocity, the aerodynamic pressure amplitude in the front regions of the pantograph is lower than that in the rear regions of the pantograph. This may be explained by the fact that the front parts always collide with the incoming flow sharply, forming a windward stagnation area. As a result, the fluctuating pressure is relatively small. Flow with higher levels of vorticity, as shown in Figure 13, passes over the rear parts of the pantograph. As a result, the amplitude of aerodynamic pressure fluctuation in the rear region is larger. Meanwhile, the amplitude of the aerodynamic pressure decreases when the frequency increases from 200 Hz to 1000 Hz. Where the flush configuration is used under the down-pantograph condition, the maximum aerodynamic pressure fluctuation at 1000 Hz is smaller than that at 200 Hz by as much as 18 dB. The phenomenon that the maximum aerodynamic pressure fluctuation decreases with the increasing frequency also can be observed by Zhang et al. [15]. At the low frequency, the SPL distribution of the aerodynamic pressure on the surface of the pantograph is non-uniform. The largest amplitude aerodynamic pressure mainly features on the supporting slide, the rotating shaft, the lower arm, and the leeward surface of the installation base. At 1000 Hz, the aerodynamic pressure amplitude distribution is more uniform than that at 200 Hz. It seems that, although the aerodynamic pressure amplitude distribution over the pantograph surfaces are different at these two frequencies, the largest amplitudes are always located on the downstream base frame, the connection between the upper and lower arms and the downstream position of the double skateboard bow. These aerodynamic pressure fluctuations act as dipole noise sources.

Aerodynamic Noise Radiation and Attenuation Characteristics
To study aerodynamic noise radiation and attenuation characteristics induced by the up-and down-pantograph with the different base configurations, the contours of spatial distribution in nearfield and far-field noise at three planes in the y-direction i.e. y = 0 (near-field), 10 and 20 m (in farfield) are compared, as depicted in Figure 16. In this analysis, the acoustic finite element method is directly used to deal with the whole field noise, including near-field y = 0 plane noise, without considering the effect of quadrupole noise and corresponding convection effect. The main reason is that Mach number of the incoming flow is only 0.278 in this calculation. As it is known that the ratio of quadrupole noise to dipole noise is proportional to the square of Mach number [32], the influence of the quadrupole noise at this speed is relatively small, and it is reasonable to ignore it in the present study. As the pantograph is the noise source, the sound pressure decreases with the y distance. The contours also indicate that the regions with higher sound levels for both frequencies are in the vicinity of the pantograph. These distributions are more uniform at the high frequency. It is shown that, as frequency increases, the sound pressure amplitude decreases.

Aerodynamic Noise Radiation and Attenuation Characteristics
To study aerodynamic noise radiation and attenuation characteristics induced by the up-and down-pantograph with the different base configurations, the contours of spatial distribution in near-field and far-field noise at three planes in the y-direction i.e. y = 0 (near-field), 10 and 20 m (in far-field) are compared, as depicted in Figure 16. In this analysis, the acoustic finite element method is directly used to deal with the whole field noise, including near-field y = 0 plane noise, without considering the effect of quadrupole noise and corresponding convection effect. The main reason is that Mach number of the incoming flow is only 0.278 in this calculation. As it is known that the ratio of quadrupole noise to dipole noise is proportional to the square of Mach number [32], the influence of the quadrupole noise at this speed is relatively small, and it is reasonable to ignore it in the present study. As the pantograph is the noise source, the sound pressure decreases with the y distance. The contours also indicate that the regions with higher sound levels for both frequencies are in the vicinity of the pantograph. These distributions are more uniform at the high frequency. It is shown that, as frequency increases, the sound pressure amplitude decreases. To obtain the noise characteristics in the far field, two positions in the y-direction, i.e. at y = 5 m and y = 25 m from the center line of the train, as shown previously in Figure 4, are evaluated, as plotted in Figure 17. Some phenomena are observed. Firstly, the noise amplitude at y = 25 m is lower than that at y = 5 m at all frequencies. Next, it is shown that the noise caused by the down-pantograph is higher than that by the up-pantograph. Finally, the noise from the pantograph with the flush base configuration is of higher amplitude than that with the sunken base configuration irrespective of whether the pantograph is under the up or down condition. These phenomena could be explained by the fact that, for the sunken configuration, the vortices generated in the rear of the insulator and the base move inside the cavity and the movement to the surroundings is blocked by the cavity surface. This may explain the predicted weakening of the far-field noise. Additionally, one can see that, for frequencies below 200 Hz, the noise caused by the sunken type configuration is almost the same as that from the flush. However, at higher frequencies, the noise caused by the sunken configuration of the pantograph is lower than that by the flush configuration in general. This is because the large vortices produced by the sunken insulator are blocked by the cavity, therefore their To obtain the noise characteristics in the far field, two positions in the y-direction, i.e. at y = 5 m and y = 25 m from the center line of the train, as shown previously in Figure 4, are evaluated, as plotted in Figure 17. Some phenomena are observed. Firstly, the noise amplitude at y = 25 m is lower than that at y = 5 m at all frequencies. Next, it is shown that the noise caused by the down-pantograph is To obtain the noise characteristics in the far field, two positions in the y-direction, i.e. at y = 5 m and y = 25 m from the center line of the train, as shown previously in Figure 4, are evaluated, as plotted in Figure 17. Some phenomena are observed. Firstly, the noise amplitude at y = 25 m is lower than that at y = 5 m at all frequencies. Next, it is shown that the noise caused by the down-pantograph is higher than that by the up-pantograph. Finally, the noise from the pantograph with the flush base configuration is of higher amplitude than that with the sunken base configuration irrespective of whether the pantograph is under the up or down condition. These phenomena could be explained by the fact that, for the sunken configuration, the vortices generated in the rear of the insulator and the base move inside the cavity and the movement to the surroundings is blocked by the cavity surface. This may explain the predicted weakening of the far-field noise. Additionally, one can see that, for frequencies below 200 Hz, the noise caused by the sunken type configuration is almost the same as that from the flush. However, at higher frequencies, the noise caused by the sunken configuration of the pantograph is lower than that by the flush configuration in general. This is because the large vortices produced by the sunken insulator are blocked by the cavity, therefore their streamwise growth and spread are constrained. As a result, large eddies are not easily broken into small ones. To further study the trend of the overall sound pressure levels (OASPL) at different distances along x (parallel to the train) and y (normal to the train) directions, Figure 18a,b shows the OASPL in the far field along the x and y directions, respectively. Figure 18a indicates that the OASPL is inversely proportional to the distance of sound source for all cases. To predict the OASPL in the region from the pantograph to the far field, it is found that the logarithm function, OASPL = OASPL0 -b  ln(x + c), is fitted appropriately using the values of the OASPL under the relation between the OASPL and the distance in the range of 0 ≤ y ≤ 30 m. Besides, under the same configuration base, the OASPL obtained by the down-pantograph is higher than that by the up-pantograph. To evaluate the OASPL caused by the pantograph in the far-field region, it is acceptable if the distance from the sound source is more than 25 m, the sound source is assumed as a point source [13]. Therefore, the sound pressure propagation can be considered as the spherical surface from the point source center and the attenuation value ∆L of the OASPL from the point source in the far-field regions is computed by the equation: ∆L = 20log(r/r0). This can be used To further study the trend of the overall sound pressure levels (OASPL) at different distances along x (parallel to the train) and y (normal to the train) directions, Figure 18a,b shows the OASPL in the far field along the x and y directions, respectively. Figure 18a indicates that the OASPL is inversely proportional to the distance of sound source for all cases. To predict the OASPL in the region from the pantograph to the far field, it is found that the logarithm function, OASPL = OASPL 0 − b × ln(x Besides, under the same configuration base, the OASPL obtained by the down-pantograph is higher than that by the up-pantograph. To evaluate the OASPL caused by the pantograph in the far-field region, it is acceptable if the distance from the sound source is more than 25 m, the sound source is assumed as a point source [13]. Therefore, the sound pressure propagation can be considered as the spherical surface from the point source center and the attenuation value ∆L of the OASPL from the point source in the far-field regions is computed by the equation: ∆L = 20log(r/r0). This can be used to the predict the geometric attenuation of noise. Figure 18b shows the OASPL along the x direction at y = 25 m. The down-pantograph always produces higher noise than the up-pantograph for both configuration bases. Besides, Figure 18a,b indicates that the OASPL in the far field along the x and y directions are always found that for all x and y positions the down-pantograph with the flush base provides the highest OASPL. The lowest OASPL is obtained by the up-pantograph with the sunken base. In addition, the OASPL caused by the up-pantograph with the flush base is higher than that by the down-pantograph with the sunken base.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 17 of 21 to the predict the geometric attenuation of noise. Figure 18b shows the OASPL along the x direction at y = 25 m. The down-pantograph always produces higher noise than the up-pantograph for both configuration bases. Besides, Figure 18a,b indicates that the OASPL in the far field along the x and y directions are always found that for all x and y positions the down-pantograph with the flush base provides the highest OASPL. The lowest OASPL is obtained by the up-pantograph with the sunken base. In addition, the OASPL caused by the up-pantograph with the flush base is higher than that by the down-pantograph with the sunken base.

Aerodynamic Noise Directivity
In this section, the aerodynamic noise directivity around the pantograph with the different configurations is presented. According to the ISO3095-2013 standard [33], the far-field noise receivers should be located 7.5 or 25 m away from the center line of the track and 1.2 or 3.5 m above the ground. The present study investigates the far-field noise characteristics around the pantograph with a distance of 7.5 m in the radial direction and 3.5 m above the ground. The OASPL in the X-Y, Y-Z and X-Z planes is discussed, as shown in Figure 19, and the OASPL is considered at every five degrees. It is found that in all the planes the noise directivity induced by the sunken pantograph is slightly smaller than the flush pantograph in general. Only at some angles, the noise induced by the sunken pantograph is equal to or larger than that by the flush pantograph. Besides, some interesting phenomena are observed in each plane. First, in the X-Y plane, when the pantograph is used with the flush surface, a relatively high OASPL is found in the rear arc of the pantograph (−5 ≤  ≤ 5), and along the sideline over the ranges 85 ≤  ≤ 125 and 235 ≤  ≤ 275. When the pantograph is installed in the sunken surface, the rear arc still has the higher OASPL. However, the high OASPL values are also predicted at other angles. In the down-pantograph, a large part of the structure is located inside the cavity. The corresponding amplitude of the four corners of the cavity is slightly higher than that of other positions. This may be because of the effects of the sunken configuration. Specifically, during the process of noise propagation, the sound wave can reflect within the cavity of the sunken surface, four corners can easily become the convergence area of reflected sound waves. One can see that the directivity of noise is not completely symmetrical the X-Y plane. Similar instances of incomplete symmetry are also reported in the literature [13]. The main reason is that the far-field noise monitors are usually placed at every 5-10 degrees around the pantograph in different planes. This may omit part of the amplitude between the monitors, resulting a certain impact on the results. In addition, in the calculation of the sound field, the frequency resolution has a slight impact on the results. However, further detailed analysis of the up-and down-pantograph in the cavity has shown that there is almost no difference when the frequency resolution is increased from 10 Hz to 5 Hz, namely, the overall noise difference is less than 1 dB. In the Y-Z plane, the noise distributions obtained by the two

Aerodynamic Noise Directivity
In this section, the aerodynamic noise directivity around the pantograph with the different configurations is presented. According to the ISO3095-2013 standard [33], the far-field noise receivers should be located 7.5 or 25 m away from the center line of the track and 1.2 or 3.5 m above the ground. The present study investigates the far-field noise characteristics around the pantograph with a distance of 7.5 m in the radial direction and 3.5 m above the ground. The OASPL in the X-Y, Y-Z and X-Z planes is discussed, as shown in Figure 19, and the OASPL is considered at every five degrees. It is found that in all the planes the noise directivity induced by the sunken pantograph is slightly smaller than the flush pantograph in general. Only at some angles, the noise induced by the sunken pantograph is equal to or larger than that by the flush pantograph. Besides, some interesting phenomena are observed in each plane. First, in the X-Y plane, when the pantograph is used with the flush surface, a relatively high OASPL is found in the rear arc of the pantograph (−5 • ≤ θ ≤ 5 • ), and along the sideline over the ranges 85 • ≤ θ ≤ 125 • and 235 • ≤ θ ≤ 275 • . When the pantograph is installed in the sunken surface, the rear arc still has the higher OASPL. However, the high OASPL values are also predicted at other angles. In the down-pantograph, a large part of the structure is located inside the cavity. The corresponding amplitude of the four corners of the cavity is slightly higher than that of other positions. This may be because of the effects of the sunken configuration. Specifically, during the process of noise propagation, the sound wave can reflect within the cavity of the sunken surface, four corners can easily become the convergence area of reflected sound waves. One can see that the directivity of noise is not completely symmetrical the X-Y plane. Similar instances of incomplete symmetry are also reported in the literature [13]. The main reason is that the far-field noise monitors are usually placed at every 5-10 degrees around the pantograph in different planes. This may omit part of the amplitude between the monitors, resulting a certain impact on the results. In addition, in the calculation of the sound field, the frequency resolution has a slight impact on the results. However, further detailed analysis of the up-and down-pantograph in the cavity has shown that there is almost no difference when the frequency resolution is increased from 10 Hz to 5 Hz, namely, the overall noise difference is less than 1 dB. In the Y-Z plane, the noise distributions obtained by the two installation bases are nearly symmetric and the higher sound level area is found in in the range of 60 • ≤ θ ≤ 120 • . It also shows that, in this range of theta with the pantograph under the down configuration, the OASPL obtained by the flush installation base is higher than that by the sunken base configuration. However, under the up-pantograph condition, the OASPL values obtained by the two configuration bases are close to each other. In the X-Z plane, the down-pantograph with the flush base produces higher noise than the sunken base in all directions. For the up-pantograph, the noise amplitudes of the flush and sunken bases are very close to each other in general. Nonetheless, it is seen that at some angles, the noise induced by the sunken base is higher than the flush base. In both Y-Z and X-Z plane, for the up-pantograph with the flush base or sunken base, it can be seen that the noise amplitudes are very close to each other in general, especially in the Z direction directly above the body. This may be due to the effect of sound wave reflection from the cavity floor and that the up-pantograph itself is a large source of noise in this direction.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 18 of 21 installation bases are nearly symmetric and the higher sound level area is found in in the range of 60 ≤  ≤ 120. It also shows that, in this range of theta with the pantograph under the down configuration, the OASPL obtained by the flush installation base is higher than that by the sunken base configuration. However, under the up-pantograph condition, the OASPL values obtained by the two configuration bases are close to each other. In the X-Z plane, the down-pantograph with the flush base produces higher noise than the sunken base in all directions. For the up-pantograph, the noise amplitudes of the flush and sunken bases are very close to each other in general. Nonetheless, it is seen that at some angles, the noise induced by the sunken base is higher than the flush base. In both Y-Z and X-Z plane, for the up-pantograph with the flush base or sunken base, it can be seen that the noise amplitudes are very close to each other in general, especially in the Z direction directly above the body. This may be due to the effect of sound wave reflection from the cavity floor and that the up-pantograph itself is a large source of noise in this direction.

Conclusions
This work presents the numerical study of the aerodynamic noise characteristics induced by a pantograph over a simplified train geometry with the different configurations of the installation base using the large eddy simulation (LES) with finite element method (FEM). The numerical results are carried out in terms of the flow field around the pantograph and its installation base, the spatial distribution and spectrum characteristics, and the noise directivity in near and far fields under the

Conclusions
This work presents the numerical study of the aerodynamic noise characteristics induced by a pantograph over a simplified train geometry with the different configurations of the installation base using the large eddy simulation (LES) with finite element method (FEM). The numerical results are carried out in terms of the flow field around the pantograph and its installation base, the spatial distribution and spectrum characteristics, and the noise directivity in near and far fields under the up and down conditions at the train speed of 350 km/h. Through the discussion, the following interesting phenomena are drawn: (1) The complex shape of the pantograph is the main reason for its induced aerodynamic noise. A series of vortices detaches from the leeward side of each rod and they interact with each other and with downstream solid structures. The aerodynamic pressure fluctuation is generated on the rods induces aerodynamic noise. The main aerodynamic noise sources are located on the downstream base frame and on the double skateboard bow, at the connection between the upper and lower arms. The aerodynamic noise sources over the pantograph surfaces are different at 200 Hz and at 1 kHz.
(2) For the spatial distributions, along the direction of the flow velocity, the aerodynamic pressure fluctuation amplitude in the front regions of the pantograph is lower than that in the rear regions of the pantograph. These distributions are more uniform at 1 kHz than at 200 Hz. In addition, the amplitude of the noise is lower at the higher frequency.
(3) In the far-field regions along the x and y directions, for all x and y positions under the pantograph shape and opening direction selected in this paper, the down-pantograph with the flush configuration provides the highest OASPL. The up-pantograph with the sunken configuration produces the lowest OASPL. In addition, the OASPL caused by the up-pantograph with the flush base is higher than that by the down-pantograph with the sunken base. Thus, the pantograph with sunken configuration is a better choice in a practical engineering application.
(4) On the directivity, the noise induced by the pantograph with the sunken surface is of lower amplitude than that with the flush surface, especially in the down-pantograph. For up-pantograph, the amplitude of the noise from the sunken surface configuration is close to that from the flush surface configuration in Y-Z plane and in the X-Z plane. In the X-Z plane, the results even show noise amplitude from the sunken surface configuration is predicted to be higher than that from the flush configurations at some angles in some angles. However, this still does not affect the fact that the sunken installation is conducive to a quieter noise environment on both sides of the road.
It should be noted that effects of noise transmission to the exterior and interior of high-speed train need to be further studied, especially in the case of the pantograph with the sunken configuration. The present study considers only the aerodynamic noise generated by the wall fluctuating pressure of the pantograph in the exterior noise transmission to the environment. In fact, there are strong vortices in the regions near the pantograph and spatial distribution of near-field, such as Y = 0 m, may obtain a more accurate result. These are the quadrupole noise sources that contribute to aerodynamic noise. Since the operating speed is not too high (Ma = 0.278), the proportion of noise generated by the quadrupole source, which scales by u 8 [34], is smaller than that from dipole sources, which scales as u 6 [35], where u is the train velocity. Thus, the induced noise by the quadrupole source is not included in this study.