In this paper, the concentration distribution of CO emitted by motor vehicles in Beijing Yuntong Tunnel is simulated and analyzed based on Fluent numerical simulation software. Firstly, the physical model of tunnel geometry is established. Secondly, the suitable control equations and boundary conditions are selected to describe the diffusion distribution of pollutants in the tunnel. At the same time, the reliability of the model is verified according to the measured data, and the distribution law of CO concentration in the tunnel is observed through the simulation results. The tunnel is unorganized at the entrance of the tunnel. In the simulation results, the CO concentration at the entrance of the tunnel is far more than the CO unorganized emission limit in Beijing. The CO emission of the tunnel should be purified to prevent pollution of the urban environment.
As a major engineering building to alleviate urban traffic pressure, urban tunnels are already very common in cities with high vehicle ownership. Beijing is the mega city. The urban tunnels are developing longer and larger. The Beijing Yuntong Tunnel is a one-way three lanes and length of 6310 m to become a representative of the city’s long tunnels. Different from outdoor, the space in the tunnel is closed, and the pollutants are easy to accumulate. If the accumulated pollutants in the urban tunnel are discharged into the atmosphere without treatment, it will affect the surrounding environment of the tunnel entrance. The tunnel ventilation and air purification system can ensure the safety of driving and the health of the surrounding residents [
Simulation experiment, numerical analysis method and actual measurement method are the main methods to study the emission and diffusion of motor vehicle pollutants in tunnels [
The control of motor vehicle exhaust pollution in big cities is very important, which is particularly important for improving the effect of environmental protection, ensuring the health of urban residents and promoting social development [
Beijing Yuntong Tunnel is located in the sub-center of Beijing and officially opened to traffic in January 2021. The tunnel passes under the whole line, the design speed is 80 km/h, the lane is two-way six lanes, and the outer lane is a large-capacity bus lane. The length of the closed section of the tunnel is 6310 m, the cross-sectional area of the tunnel is 82.3 m2, the height of the tunnel is 6.8 m, the width of the tunnel is 12.25 m, the maximum slope of the road section is +3.5%, the longest section is 993 m, and the slope is +0.17%.
In July 2021, the research team carried out field measurement on the Yuntong tunnel. A total of three measuring points were set up. The layout positions were 1500 m, 3350 m and 4500 m from the entrance of the tunnel in the direction out of Beijing (right line). The measuring points were located at the ventilation of the tunnel distribution room, close to the wall surface, and about 1.5 - 3 m high; the measured contents include: traffic flow, traffic wind speed, CO concentration in the tunnel. During the measured period, the tunnel fan runs in different time periods, and the data is obtained when the fan is not running.
The establishment of physical model is the basis of clear research object. The physical model can determine the spatial scope and geometric shape of the research, and the mathematical model can be established after the preliminary definition of the research object. The two models are mutually constrained and the appropriate numerical method is used to solve the research object more accurately. This section describes the process of establishing the physical model of Yuntong tunnel. Through reasonable simplification, analysis and summary of model establishment ideas, combined with model establishment ideas and engineering drawings, the physical model of the tunnel is established, and the physical model is meshed by software.
In order to ensure the accuracy of the simulation results, through reasonable simplification and analysis, the setting method of model boundary conditions is determined. The specific simplification of the model and the setting of boundary conditions are as follows:
① The Beijing Yuntong Tunnel was officially opened to traffic in January 2021. The measured period belongs to the trial operation stage of the tunnel, with fewer motor vehicles. The measured data is set as the boundary condition of the simulation model;
② The distribution of CO concentration is used to represent the distribution of pollutant concentration in the tunnel;
③ The distance between the vehicle exhaust pipe and the ground is compared with the height of the tunnel. It can be assumed that the vehicle exhaust pipe is close to the ground of the road. When the model is established, the vehicle in the tunnel is simplified into a flat plate to release the pollution source in the closed space;
④ Without considering chemical reaction.
In this paper, the 1:1 geometric model of Beijing Yuntong Tunnel is established by using Rhino3D drawing software. When establishing the model, the decoration of the inner wall of the tunnel is neglected and the smaller slope is neglected. The tunnel is simplified as a hexagon with a total length of 6310 m, a section of 12.25 m at the bottom, 10.12 at the bottom, and a height of 6.8 m. The simulation calculation focuses on the diffusion of pollutants.
In this study, Fluent Mashing meshing software was used. Because the tunnel was long and the cross-sectional area was small, the pollution source was simplified to an outer-swept plate for pollutant diffusion calculation. In order to simplify the calculation time and ensure the accuracy of the calculation, the hexahedral unstructured grid was used to mesh the model to reduce the number of model nodes. At the same time, according to the third item in the modeling idea, the vehicle pollution source discharge port is merged with the lower surface of the model, and the lower surface area is too large without mesh encryption.
In this paper, computational fluid dynamics technology is used to study the distribution of pollutant concentration in the tunnel, and the calculation results are verified by experimental fluid mechanics. The research of computational fluid dynamics usually follows the following steps: first, determine the geometric description of the flow region, including the geometric shape of the flow field, the flow conditions and the requirements of numerical simulation; secondly, the control equations and boundary conditions are selected. Then determine the meshing strategy and numerical method; then the corresponding CFD software is used to solve the fluid dynamics equation and obtain the numerical solution. Finally, the numerical solution is post-processed to display and analyze the flow field.
In this paper, Fluent 2020 is used to calculate the wind speed field and pollutant
concentration distribution in the tunnel. Determining the research problem is the first step in fluid mechanics calculation, including the determination of flow field geometry, flow conditions and numerical simulation boundary conditions. Secondly, the suitable control equations and boundary conditions are selected to describe the pollutant diffusion distribution in the tunnel model. The N-S equation can describe the movement of Newtonian fluid. The N-S equation includes mass conservation equation, momentum conservation equation and energy conservation equation.
Considering that the fluid flow inside the tunnel is mainly the turbulent flow of ideal fluid, and this paper mainly studies the distribution of pollutant concentration. The most widely used turbulence model is selected to describe the turbulent flow of pollutants in the tunnel. The standard K-ε two-equation is used as the control equation of the gas flow state in the tunnel. The advantage is that the standard K-ε model has strong adaptability to the velocity field. The pressure-velocity coupling SIMPLEC algorithm is used to accelerate the convergence.
a) Law of conservation of mass
The integral type is:
∂ ∂ t ∭ Ω ρ d Ω + ∯ s ρ V ⋅ n d S = 0
where Ω is the control of the integral equation, dΩ = dxdydz, S = ∂Ω is the surface of the control body, ρ is the density, and V is the velocity vector of the fluid motion.
The differential equation is:
∂ ρ ∂ t + ∇ ⋅ ρ V = 0
b) Law of conservation of momentum
The integral type is:
∂ ∂ t ∭ Ω ρ V d Ω + ∯ s ρ V V ⋅ n d S = ∭ Ω ρ F d Ω + ∯ s ρ ( − P I ¯ + ζ ¯ ¯ ) ⋅ n d S
The differential equation is:
∂ ( ρ V ) ∂ t + ∇ ⋅ ρ V = ρ F + ∇ ⋅ ( − P I ¯ + ζ ¯ ¯ )
where I ¯ is the unit tensor, p is the pressure, and ζ ¯ ¯ is the viscous stress tensor.
c) Law of conservation of energy
When there is no heat source in the flow field, the integral energy conservation equation is:
∂ ∂ t ∭ Ω ρ E d Ω + ∯ s ρ E V ⋅ n d S = ∭ Ω ρ F V d Ω + ∯ s [ ( − P I ¯ + ζ ¯ ¯ ) ⋅ V − q ] ⋅ n d S
The differential equation is:
∂ ( ρ E ) ∂ t + ∇ ⋅ ρ E V = ρ F ⋅ V − ∇ ⋅ q + ∇ ⋅ ( − P I ¯ + ζ ¯ ¯ ) ⋅ V
where E is the total energy and q is the heat flux.
d) K-ε model:
K-ε is the most widely used turbulence model. The standard K-ε model used in this study has wide adaptability to conventional velocity field and heat transfer model.
The K equation is:
∂ ( ρ k ) ∂ t + ∂ ( ρ k μ i ) ∂ x j = ∂ ∂ x j [ ( μ + μ i σ k ) ∂ k ∂ x j ] + G k + G b − ρ ε − Y M + S K
The ε equation is:
∂ ( ρ ε ) ∂ t + ∂ ( ρ ε μ i ) ∂ x = ∂ ∂ x j [ ( μ + μ i σ k ) ∂ ε ∂ x j ] + C 1 e ε k ( G k + G 3 e G b ) − C 2 e ρ ε 2 k + S K
where Gk represents the turbulent kinetic energy generated by the laminar velocity gradient; Gb is the turbulent kinetic energy generated by buoyancy; C1e, C2e, C3e are constants.
The boundary condition of the pollution source is set to the velocity inlet boundary condition, so the pollutant emission rate represented by the mass flow rate needs to be converted into velocity for calculation.
PIARC 2019 [
Q = q e x ( v , i ) ⋅ f h ⋅ f t ⋅ f e (1)
In the formula, Q is the single vehicle emissions of CO and NOx, g/(h); qex(v,i) is the single vehicle emission factor of CO of PC in Class A area, which is related to vehicle speed and road slope. The vehicle speed of the control model is 50 km/h, g/(h); fh is the altitude correction coefficient, 1 in Beijing; ft is the annual correction coefficient, and the CO correction coefficient in 2020 is 0.91. fe is the influencing factor of other technical standards, and 1 is taken for the calculation of yuntong tunnel.
The results obtained from the PIARC2019 version are the mass flow results. Through the mass fraction of each pollutant component in the automobile exhaust, the exhaust emission rate of a single car under the blocking condition is calculated.
The number of vehicles in the tunnel is calculated by Equation (2), and taking into Equation (1) to calculate the total emission rate of all motor vehicles in the tunnel.
n = N × L × 0.001 v (2)
In the formula, N is traffic flow, veh/h; L is tunnel length, m; v is the speed, km/h.
The total vehicle exhaust emission rate formula is as follows:
Q = M exhaustemission × n (3)
Because the Velocity-inlet boundary condition is assigned to velocity, not mass flow. Therefore, it is necessary to convert it into volume flow rate, and then obtain the corresponding exhaust outlet speed, from formula (4)-(6).
ρ = M V = M O 2 , CO 2 , CO , H 2 O , N 2 ρ ′ O 2 , CO 2 , CO , H 2 O , N 2 (4)
G = Q ρ (5)
v ′ = G ρ × A = Q ρ × ρ × n (6)
where, ρ is the density of pollutants, and the exhaust density of single vehicle is 1.23 kg/m3; V is the total volume of exhaust gas, m3; M is the mass of each component, O2 = 0.0224, CO2 = 0.1689, CO = 0.0049, H2O = 0.0754, N2 = 0.7284; the density of each component is O2 = 1.429 kg/m3, CO2 = 1.977 kg/m3, CO = 1.25 kg/m3, H2O = 0.6 kg/m3, N2 = 1.25 kg/m3. G is the volume flow, m3/s; v ′ , is the exhaust outlet velocity, m/s; A is the area of pollutant emission source.
In this study, the model volume is large, the ideal grid size is small and the calculation result is accurate. Affected by the actual conditions, the larger the grid density and the more grids, the longer the calculation time. Therefore, this simulation study simplifies the grid into a hexahedral unstructured grid, with a total of 1.2 million grids.
The tunnel entrance and exit are set as the velocity inlet boundary condition. According to the measured data, the average wind speed is 3.6 m/s when there are few vehicles in the tunnel. Therefore, the tunnel entrance boundary condition is set to V = 3.6 m/s, and the airflow direction is the driving direction. The tunnel outlet boundary is set as the pressure outlet to suppress the backflow, and the tunnel wall is set as the adiabatic and fixed wall boundary conditions. The measured tunnel traffic flow is 1400 veh/h. The lower surface of the model is used as the emission source, and the emission direction is positive in the Z axis. According to Chapter 3.2, the tunnel pollutant emission rate can be calculated to be 1.037 × 10−6 m/s.
The above simulation results are all natural conditions. The mechanical ventilation system inside the tunnel is all closed. The CO distribution map at 2 m height in the tunnel is shown in
From the longitudinal profile of
1) Yuntong tunnel is located in the sub center of Beijing city, special geographical location in the tunnel construction without ventilation shaft emissions
| Measuring point | wind velocity (m/s) | relative error | co concentration (mg/m3) | relative error | ||
|---|---|---|---|---|---|---|
| measured value | value of simulation | measured value | value of simulation | |||
| Measuring point 1 | 1.88 | 1.92 | 2.23% | 5.36 | 5.84 | 8.09% |
| Measuring point 2 | 1.94 | 2.18 | 10.97% | 7.45 | 8.43 | 11.59% |
| Measuring point 3 | 2.20 | 1.90 | 13.51% | 8.38 | 9.98 | 16.09% |
of pollutants, only through the tunnel exit for fugitive emissions, according to the simulation calculation of tunnel exit CO concentration value reached 12 mg/m3, far more than the Beijing fugitive emissions limit, need to open the purification device in the tunnel exit.
2) The ‘Highway Tunnel Ventilation Design Rules’ 2014 edition [
3) By comparing the measured data and simulation data to verify the reliability of the model established in this paper, the model can be used to predict the CO concentration of the tunnel or to simulate the diffusion of NOx, another major pollutant of motor vehicles (
The authors declare no conflicts of interest regarding the publication of this paper.
Zheng, S.L., Li, Q. and Xue, H.W. (2023) Simulation Analysis of CO Concentration Distribution in Beijing Yuntong Tunnel. Open Journal of Civil Engineering, 13, 210-220. https://doi.org/10.4236/ojce.2023.132016