Study on Numerical Simulation of Fire Danger Area Division in Mine Roadway

High-temperature poisonous smoke produced by coal mine roadway fire seriously affects miners’ lives and safety. Studying the development law of high-temperature smoke in the process of mine roadway fire and then exploring the danger of roadway are of great significance to personnel safety and post-disaster rescue. In order to study this problem, the CFD numerical simulation method is used to establish a fire calculation model based on ANSYS Fluent software in the development stage of mine fire. .e high-temperature flue gas flow in the roadway during the development stage of mine fire is simulated, and the variation law of temperature field and gas concentration field with time and space position under different levels of roadway in the development stage of fire is revealed. .e variation rules of environmental parameters, such as temperature, CO, and CO2, are obtained by numerical calculation. Based on these, the danger zones of smoke spread in fire development stage are divided by the critical values of high-temperature smoke and toxicity evaluation index, and the mathematical fitting analysis of the evolution of the dangerous area with time is carried out. .e research results have certain theoretical guiding significance for reducing underground environmental pollution and ensuring the personal safety of workers and rescuers.


Introduction
As one of the five possible major mine disasters, mine fires are a serious threat to the safety of underground workers and industrial property [1][2][3]. Once a fire breaks out underground, the mine tunnel would form a high-temperature flue gas environment, which contains toxic asphyxiating gases (such as CO, CO 2 , NO x , HCL, Cl 2 , H 2 S, SO 2 , and NH 3 ) and harmful smoke particles [4][5][6]. As the fire intensity increased during the fire development process, the mine ventilation network was destroyed, and the air flow temperature and pressure were changed [7,8]. Because the mine tunnel is a confined space, the smoke and heat generated during the fire are difficult to discharge in time. erefore, the high-temperature smoke diffuses rapidly in the tunnel, polluting the underground working environment and seriously threatening the lives of workers [9,10]. However, due to the inability to scientifically and rationally define dangerous areas, the evacuation of underground personnel and the implementation of emergency rescue operations is very difficult during underground mine fires [11,12].
Many researchers have studied the spread of high-temperature smoke during mine fires to determine the associated danger, and they have obtained fruitful results. Since 1980s, scholars all over the world have studied the flow of smoke during tunnel fires and the impact of fire on mine ventilation systems through experimental tests and numerical simulations [13][14][15]. Xu et al. [16], Egan [17], Li [18], Li and Ingason [19], and Shafee and Yozgatligil [20] and other researchers have established a full-scale tunnel and different small-scale tunnels to carry out extensive fire testing and studied the evolution of environmental parameters such as temperature, CO, and CO 2 in the tunnel during the catastrophe, and a series of quantitative studies were carried out on the relationship between the dimensionless smoke retreat length and the roadway wind speed and fire source. It is widely used in the numerical simulation of tunnel fire because the numerical simulation method has the characteristics of short test period, high efficiency, resource saving, and strong flexibility. For example, to study the distribution law of high-temperature gas and the influence of throttling in tunnel fire, Zhang et al. [15], Li et al. [21], Zhu et al. [22], Vaitkevicius and Carvel [23] considered different initial conditions and geometric factors of tunnel (height, slope, and radius of curvature). Fan et al. [1] and He et al. [24] employed the FDS numerical simulation to study the development characteristics of high-temperature smoke flow in tunnels under stable fire source power of 1 MW-20 MW and different ventilations. Chen et al. [25] regarded the fire source as a high-temperature flue gas release source with a fixed temperature of 2000 K and a fixed amount of gas release, and the distribution law of the temperature and toxic gases in the tunnel were explored to the danger of the tunnel. At present, numerous research works have been carried out on the development law of tunnel fire smoke flow, throttling effect, and other fire evolution characteristics. However, there are the relatively few studies on the distribution law of mine roadway fire smoke flow and the distribution characteristics of tunnel dangerous areas in the condition of unsteady fire source.
According to the investigation and analysis of the accident, the main causes of casualties caused by underground tunnel fire are toxic gas and high-temperature [26][27][28]. In this paper, the ANSYS Fluent15.0 was employed to simulate calculation of the changes of roadway temperature and smoke flow characteristics during the fire evolution in the small-scale straight tunnels, which used the t2 fire model to treat the fire source as a source term that changes with time, and set the toxic asphyxiant gas release source term. In addition, based on the foundation of previous studies, standards for dividing tunnel dangerous areas were established to monitor the changes of temperature and gas concentration in the tunnel and determined the distribution characteristics law of the fire dangerous area in the tunnel.

Numerical Simulation Method
Fluent and ICEM, as embedded modules of ANSYS third-party commercial software, are widely used in fluid calculation, heat transfer, and chemical reaction simulation. In this paper, the geometric model is established by ICEM, the fire source model is set by user-defined function (UDF), and the model is imported into FLUENTsoftware to calculate the development law of high-temperature toxic smoke in roadway fire process.

Physical Model and Fire Source.
e physical model is based on the geometric prototype of the small flat roadway fire test bench of Xi'an University of Science and Technology [18]. A geometric model with a ratio of 1 : 1 was established. e experimental roadway model is 11 m long, 0.6 m wide, and 0.4 m high.
e roadway section has a continuous, constant arch. e fire source was set to a circular imprinted surface with a radius of 0.06 m, located bottom of roadway 1 m from the air inlet. Ten measuring points were laid along the roadway model to monitor the simulation data results. ese points were arranged in two layers along the central axis of the roadway and numbered as levels L 1 and L 2 , respectively, which were 0.1 m and 0.2 m below the top of the tunnel arch. e concrete model is shown in Figure 1. e commonly used formulas for calculating heat release rate of unsteady fire sources are t 2 , complete development model, MRFC heat release rate calculation formula, and FFB index calculation model. t 2 fire is the most famous unsteady fire source, which is widely used in North America to describe the relationship between heat release rate and time during the development of fire [29][30][31]. In this paper, the t 2 model was used as the fire source term to calculate the temperature change during the fire process, and the fire source power is embedded in Fluent 15.0 [32] through UDF code of fire source power [33] for unsteady state numerical simulation calculation: where Q is the heat release rate in kilowatts, b is the development coefficient, taken as 0.0469 kW/s −2 , and t is the fire development time in seconds. UDF code of fire source power [33]: e release process of toxic asphyxiant gas was realized by a fixed source term. Studies have shown that, after coal reaches its ignition point (673 K), the CO 2 concentration of the outlet reaches a maximum of about 350,000 ppm, and the CO concentration reaches a maximum of about 60,000 ppm [34]. In this paper, the stable release of flue gas was achieved by setting CO and CO 2 gas inlets at the fire source.

Control Equation.
e high-temperature smoke spreading process in roadway fire is a complex three-dimensional unsteady fluid flow including heat and mass transfer. e numerical simulation satisfies the following conservation equations.
(1) Continuity equation [35]: In the formula, ρ is the fluid density, t is the time, and u i is the velocity component of the velocity vector u in the direction of i. (2) Momentum conservation equation [36,37]: where μ is the dynamic viscosity (absolute viscosity), S i is the generalized source term in the direction of i, and gradu i is the velocity gradient in the direction of i. (3) Energy conservation equation [28]: where T is the temperature of the microelement, a is the heat transfer coefficient of the fluid, c p is the specific heat capacity of the fluid, and S T is the viscous dissipation term, including the internal heat source of the fluid and the heat energy generated by the viscous interaction. At the fire source location, S T is the fire source power source term including the t 2 model in equation (1). (4) Component mass conservation equation [38]: where c i is the volume fraction of component i in the fluid, D i is the diffusion coefficient of component i in the fluid, and S i is the volume fraction i in unit time. At the fire source, S j is the fixed source term of CO and CO 2 . (5) Turbulence equation: incompressible fire smoke flows in the tunnel as a three-dimensional turbulent flow processes, in which the governing equations of the fluid area and the three-dimensional flux in the unsteady and turbulent state were the time-averaged mass Navier-Stokes combined with the k-e realizable model [39,40]. In the simulation process, the standard k-epsilon turbulence model equation can be used to describe, including the turbulent kinetic energy equation and the turbulent dissipation rate equation, as shown in equations (6) and (7), respectively: where k is the turbulent kinetic energy, ε is the dissipation rate, xi is the three-dimensional direction, μ t is the turbulent viscosity coefficient, σ k is the turbulent Prandtl number of the turbulent kinetic energy equation, 1.0, G k is the average velocity gradient turbulent kinetic energy, G b is the turbulent kinetic energy generated by buoyancy, σ ε is the turbulent  to carry out comprehensive numerical simulation of the development and change pattern of mine fire according to the actual physical phenomena in the fire combustion process. erefore, the combustion model needs to be reasonably simplified. Based on previous research, our simulation used the following assumptions [24,41]: (1) e air flow in the roadway is regarded as an ideal gas, ignoring the impact of the air heat expansion in the lane, and no chemical reaction occurs between the gas components there (2) Assuming that the roadway model is a horizontal smooth roadway, the wall of the roadway is regarded as one with adiabatic constant temperature, no slippage, the air velocity and temperature at the inlet are not affected by the fire, and both the wall and inlet air flow temperature are 293 K (3) e main components of the smoke generated by fire are CO and CO 2 , and other toxic gases are not considered (4) During the development of the fire, the location of the fire source does not move, and the shape of the tunnel remains intact without deformation (5) Assuming that all the heat in the fire process comes from the heat source model, no other heat is generated, and there is no mass loss of products

Boundary Conditions and Parameter
Settings. e numerical simulation is a transient calculation. In the simulation calculation, an open energy model, the RNG kepsilon (2-equ) turbulence model, was chosen, and gravity and buoyancy effects were considered. Species model selects species transportation model and opens the P1 radiation model. e initial temperature of the fire source is set to 673 K. e inlet boundary condition was set as a velocity inlet, and the inlet air flow temperature was set to 298 K. Research shows that the optimum wind speed to promote the development of fire is 0.4 m/s in the same size test-bed [42]. erefore, the air velocity value was selected at the entrance boundary. e turbulence intensity was 4.15%, and hydraulic diameter was 0.73 m. e boundary conditions of the air outlet were set to fully develop the outflow. e wall of roadway is a nonslip wall and surface velocity and diffusion coefficient were set to zero regardless of wall penetration. e SIMPLE pressure-velocity coupling method is used to solve the problem iteratively.

Grid Determination.
In roadway fire simulation, the grid size needed to be reasonably selected according to computer resources and fire characteristic parameters, and the fire source characteristic size D * was the most relevant and important parameter [43]. Scholars have conducted a series of research on grid rules; for example, Hwang and Edwards [44] used a grid size of 0.06 D * -0.12 D * for fire simulation, and the result is scientific and reasonable. Li [45] considered a grid size of 0.075 D * can satisfy the calculation conditions after analyzing the grid size accuracy. D * can be calculated by where Q is the heat release power, KW, ρ 0 is the air density, ρ 0 � 1.2043 kg/m 3 , when the temperature is 293 K, C p is the constant pressure specific heat, C p � 1.004 KJ/(Kg·K), T a is the ambient temperature, 293 K, and g is the acceleration of gravity, 9.81 m/s 2 . When the calculation time is 60 s, according to formula (1), the heat release rate of the fire source reaches 168.84 KW, and D * is 0.47 m. erefore, according to the above scholars' research results, the grid size reaches 0.075 D * � 0.035 m which can be required in the calculation. Combined with computer performance, this paper chose a grid size of 0.01 m to calculate the result. In addition, the grid is encrypted at the fire source location, and the expansion layer is added at the tunnel wall, as shown in Figure 2. e minimum size of the grid after partial densification was 3.97 × 10 −5 m, and the whole grid model contained 654,060 cells.

Model Validation and Calculation Time Determination.
e continuity is set in the numerical simulation process, the convergence criteria of k and epsilon are both set to 10 −5 , the convergence criteria of other variables are set to 10 −3 , the calculation time step is 0.1 s, the number of time steps is 600, and the simulation result is verified. Since few reports on tunnel fire experiments on the t 2 model and the t 2 fire source power model was a function of time, the heat release was not affected by other variables. erefore, refer to the experimental test method in [42], which has set the fire source as a stable fire source of 3.14 KW, and other parameters were calculated to prove the correctness of the parameter settings in the model. e corresponding point (define the point as P v -#3) with the measuring point #3 (5.4 m from the air inlet, 0.3 m height) of the experimental platform in the physical model to compare the results is selected. rough monitoring, the numerical simulation results of measuring point P v -#3 at different times are obtained, as shown in Figure 3. It can be seen that the numerical simulation results are in good agreement with the experimental results. In the initial stage, the temperature rises slowly and starts to rise rapidly at 15 s; meanwhile, the growth rate of the simulation results is higher than that of experiment. e final temperature is stable at 301-301.5 K, and the results of simulation are higher than that of experiment as a whole, which may be caused by the heat loss of high-temperature smoke flow in the experimental process. e comparison shows that the model can better reflect the flow characteristics of the smoke flow in roadway and meet the basic requirements of roadway fire risk area division.
rough simulation calculation, under the calculation time of 60 s, the high-temperature flue gas can fill the entire tunnel. Consider the evacuation of people during the mine tunnel fire, and the escape time of miners should be considered when determining the numerical simulation calculation time. In [46], scholars have carried out experiments on the escape speed of people under different roadway types, slopes, and cross-section conditions. In addition, based on the experimental results, the model of the escape speed of miners during a mine fire was established. e results are as follows: where V 1 and V 2 are the miners' escape speed, m/s, S is the roadway slope,°, A is the tunnel section, m 2 , and R t is the tunnel type, which including track and trackless types. According to formula (9), combined with the case of this paper and physical similar conditions [10], the maximum escape speed of miners in a straight trackless tunnel was 4.342 m/s, and the escape distance within 60 s was 260.52 m, so the calculation time of 60 s was selected.

Results and Discussion
3.1. Distribution of Temperature Field. As the simulation proceeded, the simulated gas flow process was monitored by 10 measurement points at levels L 1 and L 2 , and the temperature data at each point were recorded and analyzed. e variation pattern of temperature by roadway length at different time points is shown in Figure 4. e variation pattern of temperature at level L 1 is shown in Figure 4(a). It can be seen from the figure that the temperature at each measuring point has a similar pattern for each time step over the length of the roadway. e highest temperature is 1.4 m downwind of the fire source, the temperature rises gradually over time, and the highest temperature reaches about 319 K. With a gradual increase in the length of the fire source, the temperature shows a downward trend, and at 10 s and 20 s, the temperature is still low. is is because the heat generated by the fire source in the initial stage has not yet been transmitted to the downwind side, resulting in a temperature change that is not obvious, but a stable heat transfer phase is achieved in 20-30 s of development.
e variation pattern of temperature at level L 2 is shown in Figure 4(b). It can be seen from the figure that the temperature variation of the L 2 horizontal measuring points has a similar pattern. e highest temperature point is 0.4 m downwind from the fire source, and the temperature gradually increases with time. e maximum temperature is about 360 K. With the gradual increase of the length of the fire source, the temperature shows a rapid downward trend, and the temperature drop pattern is obvious at the time of 10 s and 20 s. Again, this is because the heat generated by the fire source in the initial stage has not been transmitted to the downwind side, resulting in only a small temperature change, but stable heat transfer starts at 20-30 s.
A temperature cloud map across the vertical section of the axis of the roadway at different time points is shown in Figure 5. It can be seen that the high-temperature flue gas at the fire source flows upward due to buoyancy, and it is affected by the horizontal wind flow of the roadway so that the flue gas flows along toward the downwind side. e fire pressure and high-temperature smoke generated by the fire disturb the steady wind flow in the roadway and create turbulence. As downwind of the fire source, the horizontal flow is at an acute angle with respect to the top of the roadway. After the flue gas generated by the fire source rises to the top of the roadway, most of the flue gas flows downwind from the fire source, but some flue gas flows along the top of the roadway to the upwind side of the fire source to form a countercurrent flue gas layer. It can be seen that the counterflow length gradually increases with time. After the high-temperature smoke from the fire spreads to the downwind side of the roadway, the average temperature in the roadway gradually increases with time but decreases with increasing distance from the fire source. Since the smoke generated by the fire accumulates near the ceiling along the entire roadway, the flue gas stratifies, expands, and  Mathematical Problems in Engineering spreads stably in the roadway, resulting in the temperature of the vertical space of the roadway ranging from high to low, with higher temperatures at the ceiling. It can be seen from the temperature cloud map that the high-temperature wind flow basically spreads along the entire roadway at 30 s, and by 60 s, the outlet temperature reaches 310 K.

Distribution of Smoke Concentration Field.
During the development of the fire, the smoke generated is entrained by the wind and also spreads along the roadway. rough calculation and postprocessing, the distribution patterns of CO and CO2 gas concentrations in the longitudinal axial planes of the roadway at 10 s, 20 s, 30 s, 40 s, 50 s, and 60 s were obtained, as shown in Figure 6.
It can be seen from Figure 6(a) that the highest CO concentration is about 5% at the source of the fire. During the diffusion of the smoke stream, the CO concentration near the fire source remains basically unchanged and the flue gas generated by the fire diffuses to the roadway ceiling under the action of buoyancy. As the flue gas continues to diffuse in the wind flow through the roadway, the CO concentration decreases with increasing distance from the fire source. e CO gas concentration contours in the roadway are layered in the horizontal direction. It can be seen from the distribution characteristics of the CO gas concentration contour in the vertical direction that the smoke flows along the top of the roadway due to buoyancy. e CO concentration at the roadway floor is lower than that of the top of the roadway. At about 30 s, smoke from the fire has spread along the entire roadway.
It can be seen from Figure 6(b) that the concentration and variation of CO 2 concentration in the flue gas have a pattern similar to that of CO. e CO 2 concentration at the fire source is the highest, about 30%. e CO 2 concentration near the fire source did not change much during the development of the fire. It can be seen from the contours of CO 2 concentrations that the smoke near the fire source spreads to the top of the roadway due to the heat created by the fire. As with the CO pattern, downwind flow through the     roadway diffuses the smoke, and the CO 2 concentration gradually decreases as the distance from the fire source increases. In addition, due to buoyancy, the smoke accumulates at the top of the roadway, and the CO 2 gas concentration contour is hierarchically distributed in the vertical direction; at 30 s, the CO 2 is diffused along the entire roadway.

Classification of Fire Danger Zones
In a coal mine fire, the fire source will release a large amount of heat and toxic gases. Due to the wind flow, it spreads widely through the limited space of the roadway. After a certain time, the high-temperature flue gas products are gradually distributed throughout the roadway space, certainly forming an oxygen-deficient and poisonous environment, which poses grave danger to the life and safety of underground workers [47]. erefore, based on the simulation results, the authors used the temperature and toxic gas distribution patterns at the measurement points to classify the fire danger areas of the roadway based on the selected fire risk assessment indicators.

Assessment Indicators and reshold
Value. e main cause of casualties in mine fire accidents is suffocation and poisoning by high-temperature smoke. Previous survey results showed that among the main combustion products of fire-CO, CO 2 , SO 2 , and other high-temperature organic gases-the largest cause of casualties is CO and CO 2 in hightemperature flue gas [48]. erefore, this paper selects the temperature, CO, and CO 2 as indicators. According to the research results of scholars [5,25,[49][50][51][52][53][54], the impact of different indexes on human body and critical values are shown in Table 1. According to the degree of influence of high-temperature flue gas on the human body, thresholds were selected to classify risk in the flue gas-spreading area.

Flue Gas
Temperature. According to the influence of different temperatures on the human body, 303.15 K, 313.15 K, and 338.15 K are used as the temperature thresholds to classify the danger zone according to heat injury in a mine fire. e downwind side of the roadway was divided into four areas: safe area (I), slightly dangerous area (II), moderately dangerous area (III), and highly dangerous area (IV). e same classes were used for all three indicators. e specific classification thresholds are shown in Table 2.

CO Concentration.
According to the influence of different CO gas concentrations on the human body, 0.005%, 0.32%, and 0.64% are used as the CO concentration thresholds for hazard classification. ese thresholds are shown in Table 3.

CO 2 Concentration.
Carbon dioxide itself is nontoxic, but it asphyxiates people by displacing oxygen, and it is a major cause of fire casualties, especially in enclosed spaces. Its long-term allowable concentration is 5,000 ppm. According to the influence of different CO 2 gas concentrations on the human body, 0.5%, 3%, and 6% are used as the CO 2 concentration thresholds to create the hazard classifications. e thresholds are shown in Table 4.

High-Temperature Hazardous Areas.
We used the simulation results for the roadway distribution pattern of the flue gas temperature at levels L 1 and L 2 for each time step to draw two critical temperature classification lines at 303.15 K and 313.15 K. Furthermore, the horizontal temperature field danger zone was also classified, and the result is shown in Figure 7.
It can be seen that, during the first 60 s of fire development, the maximum temperature at level L 1 reaches 360 K, which exceeds the critical value of danger class IV. e highest temperature at the L 2 level reaches 320 K, which is within the danger class III range. During this time, both the level L 1 and L 2 locations having danger class I gradually disappear on the downwind side of the fire source and become class II zones. e critical temperature classification was applied to the results of the roadway horizontal distribution of different high-temperature danger zones at each time step, and the results are shown in Table 5  initial fire development (10 s), although a few similar areas exist downwind of the fire source at this time. Two danger class II areas and additional class III areas are distributed downwind of the fire source. As the fire develops, the danger class I areas within the roadway remain essentially unchanged, but the range of danger class II areas is constantly changing. Figure 8 uses the distribution result for danger level classification of the roadway to illustrate the variation of the high-temperature danger zone range of measurement levels L 1 and L 2 over time. It shows that the area identified danger class I is the largest before 10 s; class I areas are rapidly reduced from 10 s to 20 s, and the size of the danger class II areas increase. After 20 s, the danger class I areas of both levels L 1 and L 2 remained unchanged, but the danger class II area began to decrease linearly with time, with reduction rates of 0.17 m/s and 0.06 m/s for levels L 1 and L 2 , respectively. e danger class III areas for levels L 1 and L 2 increased linearly, with growth rates of 0.17 m/s and 0.05 m/s, respectively, and the linear growth rate of the danger class IV area of level L 2 was 0.01 m/s. In general, the unsafe area determined by the L 1 horizontal temperature is larger than that determined by level L 2 .

Smoke Spread Danger Zone.
All the gaseous components of the flue gas lead to poisonous asphyxiation and death. erefore, when smoke spreads through the roadway, it is necessary to evaluate both the CO and CO 2 levels to classify areas with asphyxiation hazards. From the simulation results for smoke spreading, it can be seen that, during the mine fire development, the gas generated by the fire source will gather near the ceiling and spread downwind Lose consciousness after several breaths 1-2 minutes, resulting in death    [55,56], the gas concentrations of CO and CO 2 at the L 2 level were selected to classify the dangerous areas in the roadway. e smoke toxicity risk assessment index was used to draw the line representing the CO and CO 2 concentration curves at the classification thresholds, as shown in Figure 9. It can be seen from Figure 9 that, during fire development, as the CO and CO 2 gases generated by the fire source flowed downwind from the fire source, their concentrations in the roadway gradually increased with time. Areas closer to the fire source experienced greater increases in gas concentration, which rapidly crossed the thresholds for successively more dangerous classifications. e distribution results were used to classify the smoke danger at each horizontal point of the roadway over time, as shown in Table 6. It can be seen from Table 6 that some areas of the roadway, classified as toxic according to CO and CO 2 hazard, are coincident. erefore, in the final result determination, the CO and CO 2 gas toxicity classification results were combined and further divided based on the results of the high toxicity level.    e results of the reclassification are shown in Table 7. It can be seen from Table 7 that, in the development stage of a mine fire, the toxic environment in the roadway can be divided into four areas. In the early stage of the fire, the safe area in the roadway is mainly on the upwind side of the fire source and the roadway exit. As the fire continues to develop, the flue gas spreads to the downwind side, resulting in a gradual increase in the risk level throughout the tunnel, and the danger class II area is distributed on the downwind side of the fire source. As the fire develops, the safety zone in the roadway remains basically unchanged for more than 30 s, but the range of the more hazardous areas is constantly changing.
e result of the classification was processed to view how the range of the danger zone of the roadway changes over time, as shown in Figure 10. As can be seen from Figure 10, the safe area of the roadway is the largest before 25 s. With the development of the fire, the safe area downwind of the fire source rapidly shrinks over time, and the range of areas with danger classes II and III are initially small and then increase and decrease with time, with a small rate of change.
e range of danger class IV increases rapidly at first, and then, its growth rate slows. e change of the range of danger class I and IV has obvious regularity. e following relationship between the range and time of the dangerous areas is obtained by fitting: , where X D is the is the length of the danger area, m, and t is time, s.

Results of Comprehensive Analysis on Hazardous Areas.
As mine fires develop, the amount of generated heat and smoke increases greatly, and the safe areas in the roadways are simultaneously threatened by high temperatures and toxins. During the fire development process, the unsafe (class II) high-temperature area of the roadway near the tunnel ceiling (level L 1 ) is larger than that at the middle height of the tunnel (level L 2 ). e horizontal high-temperature class II area range linearly decreases over time, with a reduction rate of 0.17 m/s, and the class III area range increases linearly at a rate of 0.17 m/s. However, compared with the high-temperature danger areas, the toxicity of the flue gas creates larger areas with higher danger classification zones. e area safe from a toxic atmosphere shrinks rapidly over time, from 0.15 m/s to 0.38 m/s. e areas with danger classes II and III are small, and the most dangerous class IV area grows rapidly, up to 0.32 m/s.

Conclusions
e numerical simulation results show that, during the mine fire in a straight roadway, the temperature and gas concentration form a gradient from high to low in the downward vertical direction, and part of the flue gas flows along the top of the roadway toward the upwind side of the fire source to form a flue gas counterflow. e numerical simulation results are in good agreement with the experimental results of relevant scholars.
In this study, the temperature and toxicity of flue gas are taken as evaluation indexes, and the risk assessment of flue gas temperature and toxicity is carried out according to the harm degree of different indexes to human body. Based on previous research foundations, the classification standards of dangerous areas were further excavated, and the critical value of classification of dangerous areas was determined, and the smoke spread area was divided into four classifications, ranging from safe (I) to highly dangerous (IV).
During fire development, areas made unsafe by high temperatures are larger near the tunnel ceiling (level L 1 ) than  those at the middle height (level L 2 ). e danger class II area at level L 1 linearly decreased over time, and the area with danger class III increased linearly. Compared with the range of high-temperature danger areas, areas with dangerous levels of toxic flue gas were larger and more dangerous. Combined with high-temperature and toxic gas index, the relationship between the range of dangerous area and time is determined. rough the research of numerical simulation, the idea and method are provided for the division of dangerous area and the judgment of rescue danger of mine tunnel fire. According to the numerical simulation results and the physical similarity simulation criteria, the dangerous area of the roadway during the mine fire under actual site conditions can be calculated. However, because the simulation in this paper is based on the calculation model established by the small-scale tunnel fire experiment platform, there may be some deficiencies in the results of physical similarity calculation, which needs to be further improved in the future research.

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

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.