1. Introductory
Mine fires are a significant concern among scholars worldwide, with 60% of coal mine fires being exogenous fires. Among these, tape fires are the most prevalent, often caused by prolonged friction between the tape and bottom rollers during the conveying process. This friction can lead to temperature increases, igniting any coal residue present. When a fire breaks out in a mine, it typically alters the ventilation resistance of the affected area, impacting the ventilation system and creating turbulence in the underground airflow [
1] (pp. 273–276). Once ignited, toxic gases like carbon monoxide (CO) emitted by the fire can quickly escalate and spread to unaffected areas, posing a threat to the safety of individuals and the overall evacuation process. Scholars both domestically and internationally have extensively researched the evolving characteristics of fire zone resistance. These studies have focused on the contributing factors to fire zone resistance resulting from fires as well as the impact of variables such as fire source intensity, ventilation conditions, elevation of the fire zone tunnel, and other related factors on the magnitude of resistance. Zhou Chuang [
2] (pp. 24–26) and colleagues conducted scaled-down model experiments and numerical simulations to study mine ventilation resistance in the presence of a uniform heat source. By integrating theoretical principles of thermodynamics and the equation of energy conservation, they developed a model for calculating the thermal resistance of a tunnel with a uniform heat source; Li Zongxiang [
3] (pp. 909–913) and other researchers conducted a study on local resistance and resistance coefficient of flames using mathematical derivation methods. They explained the relationship between the size of the cross-sectional area of the fire source and local resistance, provided formulas for local resistance and resistance coefficient, and validated their mathematical models through experiments. In the study of flame shape and cross-section, Ward-Smith [
4] has investigated the one-dimensional flow of smoke generated by fires by means of a one-dimensional treatment of alley fires and has constructed and optimized the differential equations for the viscous losses in the combustion zone of their flow processes by means of theoretical derivations. Fu Peifang et al. [
5,
6,
7] examined the thermal dynamics of fires on roadways, demonstrating the physical significance of thermal dynamics through numerical simulations and scaled-down model experiments. They also determined the relationship between the pressure change gradient due to thermal resistance and the distance from the fire source.
Numerous scholars have discovered that thermal resistance arises from the heating of air by a fire source, causing changes in its density and other factors that impact its resistance [
8,
9]. In the study, a strong correlation was observed between the thermal resistance magnitude and the degree of temperature rise. George and Davood [
10] used numerical simulation to study the heat transfer of the wind flow in a roadway fire, which mainly used Multiflux and CFD numerical simulation software to illustrate the laws of heat and mass transfer of the wind flow. OKA Y et al. [
11,
12,
13] conducted scaled-down model experiments to examine the longitudinal position and temperature variation of the smoke layer in a fire tunnel. They predicted temperature and wind speed variations and proposed a mathematical model for analysis; Kalech et al. [
14] (pp. 1034–1040) conducted a study on the high temperature smoke from tunnel fires in ventilated environments. They utilized FDS numerical simulation software to illustrate the temperature distribution and CO distribution characteristics. In studies conducted by various scholars on local fire zone resistance, a key consideration is the presence of a uniform heat source. However, in real mine fires, the temperature downstream of the fire source varies with distance and is not evenly distributed. Additionally, after a fire, the temperature of the tunnel in its fire zone is strongly related to the ventilation conditions [
15,
16]. The study aims to analyze the temperature change within a tunnel in a fire zone and incorporate it into the calculation of thermal resistance. This approach can enhance the accuracy of calculating the fire zone resistance. Additionally, the research investigates the components of fire zone resistance and develops a calculation model for fire zone resistance. This model can offer theoretical support for determining ventilation resistance and restructuring the wind network during fire incidents.
2. Theoretical Study on the Change of Tunnel Wind Resistance during the Catastrophic Period
When a fire occurs in a tunnel, the total resistance in the tunnel can be broken down into the ventilation resistance on the upwind side of the fire source, the local resistance of the fire source, and the thermal resistance on the downwind side of the fire source, as illustrated in
Figure 1.
Division of the tunnel into 4 intervals, with the numbers at the top of the model representing the interval segments. When the wind speed in the horizontal roadway exceeds the critical wind speed, the high-temperature smoke generated by the fire source has minimal impact on the upwind side, with only local resistance present. In this scenario, the ventilation resistance for Sections 1–2 on the upwind side of the fire source can be calculated using Equation (1), which represents the formula for ventilation resistance.
The ventilation resistance from Section 1 to Section 2 is represented by
h1-2, measured in Pa. The full pressure of Sections 1 and 2 is denoted by
p1 and
p2, respectively, in Pa. The density of the fluid in Sections 1 and 2 is indicated by
ρ1 and
ρ2, measured in Kg/m
3. The fluid flow rate at Sections 1 and 2 is represented by
v1 and
v2, m/s. In
Figure 1, Sections 2–3 discuss the local resistance, which is primarily attributed to the wind flow passing by the fire source. The fire plume creates a wind barrier that hinders the passing wind flow, resulting in an additional loss of local resistance to the thermal bypass flow. Equation (2) shows the expression for the local resistance
h′
e,r.
where
λ is the flame local resistance coefficient,
ρ′ is the density of the wind flow through the exit of the fire zone in kg/m
3, and
v′ is the exit wind speed in m/s.
The expression for the local resistance coefficient is related to the size of the flame section in accordance with the empirical formula for roadway ventilation. By considering the local resistance, we can obtain the following expressions for local resistance and the local resistance coefficient.
The relationship between the local resistance of the flame and the cross-sectional area of the flame can be observed in Equations (3) and (4).
When the horizontal dimensions of the fire source remain constant or change minimally, there is a correlation between the flame height and the rate of heat release [
17,
18], which can be characterized by Equation (5).
where
L is the flame height, m;
D is the equivalent diameter of the fire source, m; and
Q* is the dimensionless heat release rate, which is calculated as Equation (6).
The flame height can be calculated using Equations (5) and (6), and then multiplied by the horizontal size of the fire source to determine the fire stage area. This helps in assessing the resistance to bypassing flow due to the fire source. In these equations, Q represents the heat release rate in Kw, cp is the specific heat capacity of ambient air at constant pressure in Kj/(kg·K), T0 is the ambient temperature in K, g is the acceleration of gravity in m/s2, and H is the hydrodynamic diameter in m.
3. Research on the Construction of Thermal Resistance Model of Fire Zone Roadway
Figure 2 depicts an experiment on one-dimensional inviscid pipeline thermal resistance. The high-quality interface pipeline connects to an infinite container, assuming no heat or work exchange with the environment and no viscosity within the pipeline. This setup enables the derivation of conservation equations, as shown in the experimental schematic diagram below.
In
Figure 2,
P0 is the pressure of the ambient ventilation, pa;
T0 is the temperature of the ambient airflow, K;
ρ0 is the density of the ambient airflow, Kg/m
3;
P1 is the pressure of the airflow into the inlet of the heated tunnel, pa;
T1 is the temperature of the airflow at the inlet of the heated tunnel, K;
ρ1 is the density of the airflow at the inlet of the heated tunnel, Kg/m
3;
P2 is the pressure of the airflow at the outlet of the heated tunnel, pa;
T2 is the temperature of the airflow at the exit of the heated tunnel, K; and
ρ2 is the density of the airflow at the exit of the heated tunnel, Kg/m
3; Section 1 and Section 2 representing the entrance and exit sections of a heated tunnel.
The continuity equation for its tube flow experiment can be expressed by Equation (7):
The momentum equation can be expressed by Equation (8):
In the experimental setup,
P1 and
P2 represent the pressures in Sections 1 and 2 respectively, Pa; while
ρ1 and
ρ2 denote the densities of the fluid in Sections 1 and 2, measured in kilograms per cubic meter, Kg/m
3. Similarly,
v1 and
v2 represent the flow rates in Sections 1 and 2, measured in m/s. By combining Equations (7) and (8), we derive equation (9).
When fluid enters a heated pipe, the physical properties of the fluid cause accelerated thermal expansion, resulting in v2 − v1 > 0, or p1 − p2 > 0. This indicates that the pressure at Section 1 is greater than the pressure at the exit of Section 2. This phenomenon and its calculations can be described as a pressure drop in the direction of fluid flow due to heat in the pipeline, which can be referred to as thermal resistance.
Thermal resistance causes pressure changes at both ends of the fluid, leading to the construction of a thermal resistance coefficient based on the frictional resistance factor, as described in Equation (10).
The energy equation is Equation (11):
where
cp is the specific heat capacity of the fluid at constant pressure, in J/(kg × K);
T1 and
T2 represent the temperatures of Sections 1 and 2 of the wind flow in K;
e is the heat absorption per 1 kg of wind flow passing from Section 1 to Section 2, calculated as
e =
Q/
qm, in J/kg;
Q is the heat transferred to the wind flow from Section 1 to Section 2 per unit time, J; and
qm is the mass flow rate of the fluid, kg/s.
Scholar Yang Shengqiang [
17] (pp. 630) proposed using the equivalent unheated viscous flow control equation with the equivalence condition to derive a factorless criterion for obtaining the formula for thermal resistance, see Equation (12).
where the thermal resistance coefficient
Ct is expressed in Equation (13) after the introduction of the factorless criterion:
The factorless criterion, denoted as He, is defined as He = e/cpT0, where e represents the heat applied and cpT0 is the initial enthalpy of stagnation for a unit mass of fluid flowing from Cross Section 1 to Cross Section 2. Additionally, He1 is defined as g(z1 − z2)/cpT0, where g is the gravitational acceleration, z1 and z2 are the heights of Cross Section 1 and Cross Section 2 respectively, and cpT0 is the initial enthalpy of stagnation.
The thermal resistance coefficient
Ct in Equation (13) can be determined based on its relationship with temperature. The formula for temperature distribution in
Figure 2 serves as a standard, where fluid enters the heated area due to a uniform heat source. In actual mine fires, factors such as wind flow velocity, heat release rate of the fire source, distance from the source of the fire, and other variables impact the temperature of the area. Changes in temperature lead to variations in thermal resistance, with local temperature depending on distance from the fire source. Applying the relevant rules to the calculation of the thermal resistance coefficient can enable the newly constructed model to more accurately reflect the size of the thermal resistance.
5. Experimental and Numerical Simulation of a Reduced Scale Model
5.1. Experimental Platform Construction
During a visit to a mining enterprise in Shanxi Province, a scaled-down model of a roadway section measuring 4 m × 5 m was established. The model was created based on the principles of power similarity and motion similarity in hydrodynamics, ensuring that it was isometrically scaled with the actual model [
22,
23,
24] (pp. 179). The scaling relationship can be described by the formula:
where
λ is the scaling relationship between the scaled model and the full-size model,
l is the feature length,
m is the scaled model corner scale, and
p is the full-size model corner scale.
A scaled model with a geometrically similar scale relationship of approximately 6.5 is employed for fire experiments. The model’s layout is divided into 10 sections to optimize sensor placement. To ensure the precision of experimental results, an adiabatic coating is applied to the model’s roof to prevent heat loss and external heat exchange. Detailed parameters for each section of the model can be found in
Table 1.
Temperature sensors and gas concentration sensors are positioned at the top of the model, with wind speed sensors placed at the exits, entrances, and connections within the model. Additional sensors are located in the middle of the model’s roof. The temperature sensor utilizes a Pt-100 thermocouple, while the gas concentration sensor is part of a multi-gas detection device capable of detecting O2, CO, and CO2 concentrations.
At an ambient temperature of 25 °C, several pt100-type thermocouples and gas concentration sensors were strategically placed on the top of the scale model to monitor temperature and CO2 concentration variations in real-time during a fire event. The fire source was positioned 300 cm away from the model entrance, with the first sensor directly above it. Subsequent sensors were placed at 50 cm intervals in the downwind direction, and a final set of sensors were positioned at the end section of the model.
All sensors transmit real-time detection data to a computer via an RS485 communication line for recording purposes. The model layout is depicted in
Figure 5.
In
Figure 5, the red color represents the location of the fire source, the blue color represents the detection and fan equipment, the arrow shows the direction of the wind flow, the black dots are the sensor deployment locations, and the blue line segments are the communication lines between the monitoring equipment.
The ignition source material, ethanol, with the chemical formula C
2H
6O, was chosen, and its material properties are shown in
Table 2.
The ambient temperature was determined to be 25 °C and the fire experiment was conducted using its model, which was first analyzed for its O
2 and CO
2 variations, which are plotted in
Figure 6.
The graph illustrates that the O
2 concentration initially drops to its lowest point around 40 s, whereas the CO
2 concentration peaks for the first time. Subsequently, both gas concentrations fluctuate within a narrower range but remain relatively stable. This pattern of fluctuation is observed for both gases. Therefore, it can be proved that in about 40 s, its combustion state reaches the fire stability stage. The numerical simulation software Fluent2023 R2 was used to simulate the scaled fire model, and its boundary conditions and the scaled model for the same settings were used to obtain the temperature distribution of each measurement point at 40 s, as shown in
Figure 7.
In the stable phase of the fire, the temperature decay of the two experiments aligns with exponential decay as the distance from the fire source increases, as illustrated in
Figure 7. The temperature in Fluent appears relatively higher compared to the experiments, primarily due to the ideal conditions in Fluent. In the scaled model experiments, despite the use of thermal insulation coatings on the wall, heat exchange with the environment still occurs, leading to heat loss. Consequently, the numerical simulation shows slightly higher temperatures compared to the scale model experiment. The consistency in temperature trends suggests that the simulation of fire on the roadway using numerical simulation software Fluent2023 is highly reliable.
By analyzing the numerical simulated values and the temperature values of the scaled model experiments, the methodology for error calculation can be described using the equation found in the literature [
25]. It is mathematically modelled as Equation (29):
where
δ is relative error, Δ is absolute error, and
L is real value. For this section, the value is taken for the numerical solution. The error is calculated to be about 2.9%. Thus, its simulation accuracy is high.
5.2. Numerical Simulation Models
In fire simulation, the main numerical simulation software includes Fluent2023 and FDS. The choice of turbulence model and combustion model is crucial in fire simulation. Fluent offers a variety of turbulence model options, while FDS primarily focuses on the LES turbulence model. Research by Maele [
26] has shown that the LES turbulence model yields a higher critical wind speed compared to other turbulence models. Ventilation conditions in mines have a significant impact on roadway fires. Mcgrattan [
27] emphasizes that when considering different CFD models for specific fire scenarios, one should focus on model applicability. For mine fires, the numerical simulation software Fluent2023 can be used to conduct full-scale experiments on tunnel fires. SpaceClaim was employed to create a roadway measuring 200 m in length, 4 m in width, and 3.2 m in height. The meshing process was carried out using Fluent, with a rectangular cut-out measuring 1 m in length and 1m in width designated as the fire source positioned 100 m from the entrance. The model diagram can be seen in
Figure 8.
The simulation setup includes a velocity-inlet for the inlet type, a pressure-outlet for the outlet type, a mass-flow-inlet for the fire source, and a wall for the wall boundary. The turbulence model chosen is the RNG k-epsilon model with consideration for buoyancy, and the P1 radiation model is activated. A monitoring point is positioned directly above the fire source and measurements are taken every 5 m downwind to observe changes in temperature, density, pressure, and wind speed in the specified region.
In order to simplify calculations related to fire source temperature and smoke flow, the combustion reaction equations proposed by Huggett were utilized for the fire source setup. These equations can be described as Equations (30) and (31) [
28] (pp. 27–40):
Vi represents the combustible component ratio and S denotes the theoretical oxygen consumption per unit mass of fuel required for complete combustion. Huggett discovered that the rate of heat release during combustion is solely linked to oxygen consumption, with each unit mass of oxygen consumed releasing approximately 1.30 MJ of energy.
Based on the theory mentioned above, when the mass inflow port in Fluent is adjusted to a specific mass flow rate of methane fuel, it results in the consumption of 0.3 kg/m2 of methane per second and the release of 1.5 MW/m2 of heat. The combustion model in Fluent is configured as a PDF non-premixed combustion model. The fire source location is specified as 0.6 kg/m2 of methane fuel per second with a methane mass flow rate of 1. The mixing fraction of methane mass flow into the inlet is set to 1. The oxidant at the ventilation inlet is air, consisting of 0.23 oxygen and 0.77 nitrogen. When the oxidant comes into contact with the fuel, it undergoes rapid combustion, leading to a fast chemical reaction. This model does not account for the specific reaction mechanism, instead focusing on the quick generation of the final combustion product.
5.3. Critical Wind Speed Determination
Critical wind speed is defined as the velocity of wind flow at which the high temperature smoke produced by a fire does not exhibit a backdraft phenomenon. This critical wind speed is typically dependent on factors such as the size of the roadway cross-section, the shape and size of the fire source, among others. In order to further investigate the thermal resistance composition and its changing patterns, numerical simulation experiments are necessary to determine the specific wind speeds involved.
The critical wind speed of a single fire source for this tunnel condition is calculated by the critical wind speed prediction model of WU and Bakar [
29] (pp. 364–380), which is as follows:
Baker’s prediction model for dimensionless wind speed, calculated by dimensionless treatment of the heat release rate from the fire source, is Equation (32):
By calculating the dimensionless wind speed, the critical wind speed can be calculated using the critical wind speed prediction model, which is as in Equation (33):
where
v is the critical wind speed, m/s.
Based on the given parameters, it can be inferred that a mass flow rate of approximately 0.6 kg/m
2 per second of methane fuel results in a heat release rate of 3 MW/m
2. By utilizing Equations (32) and (33), the critical wind speed for this scenario is theoretically around 1.68 m/s. Utilizing Fluent to simulate ventilation speeds of 1.6 m/s, 1.7 m/s, 1.8 m/s, and 1.9 m/s, the critical wind speed was determined to be 1.7 m/s, consistent with the calculated results.
Figure 9 illustrates the temperature distribution at both critical and non-critical wind speeds.
5.4. Analysis of Numerical Simulation Results
Before commencing the numerical simulation, the full domain is initialized along with the longitudinal wind speed of the tunnel being set at 2 m/s at time 0 s, which exceeds the critical wind speed. Subsequently, the wind speeds at various measurement points within the tunnel are examined, specifically at distances of 10 m, 20 m, 30 m, and 40 m from the fire source. The results are then visualized in a variation diagram depicted in
Figure 10.
The figure illustrates that during the initial phase of the simulation up to around 90 s, there was fluctuation in wind speed at each measurement point. Beyond 90 s, the wind speed began to stabilize, with a noticeable trend of higher wind speeds closer to the fire source. The wind flow in the heated area increases the internal energy of the fluid and also performs external expansion work due to its physical properties. As the heated wind flow moves through the environment, it transfers heat, overcomes frictional resistance, and experiences a reduction in flow rate and temperature as internal energy decreases. The wind speed is higher at the lower side of the roadway compared to the entrance due to differences in temperature along the roadway causing thermal expansion of the wind flow. The thermal expansion of the fluid is most clearly observed through the decrease in density on the downwind side of the fire source. This phenomenon is illustrated by the density cloud depicted in
Figure 11.
Therefore, the flow velocity on the downwind side of its fire source is always higher than the flow velocity on the upwind side.
Researchers studied the temperature, wind speed, and fire density in the tunnel during the fire stabilization stage to investigate thermal resistance downstream of the fire.
Figure 12 depicts the variation in total pressure within the tunnel.
In
Figure 12, the deflection of the flame front due to wind speed is evident. Pressure stratification is more pronounced on the upwind side of the fire source, as seen in the same section at the bottom of the roadway where pressure is higher. As vertical height increases, there is a varying degree of pressure reduction. The pressure stratification on the upwind side of the fire source is clearly observable, with higher pressure at the bottom of the roadway and decreasing pressure with height. The flame root creates a wind barrier effect, obstructing the normal wind flow at the bottom of the roadway. This results in the highest pressure at the bottom of the roadway. The continuous area of flame acts like flexible dampers, reducing the movement of wind flow and thereby decreasing ventilation resistance at the bottom. Additionally, as the wind passes through the flame-heated area, the ventilation resistance increases due to thermal resistance, which is temperature-dependent. The high-temperature gas from the fire is influenced by buoyancy force, causing it to accumulate on the downwind side roof of the roadway. Consequently, the roof pressure in the same cross-section is higher than in other areas. The high-temperature smoke flow from the fire diminishes as distance increases, leading to a gradual reduction in thermal resistance and pressure decay, as illustrated in
Figure 13.
The study illustrated in
Figure 13 demonstrates that the thermal resistance caused by the fire decreases as distance increases in a linear fashion. The temperatures, wind speeds, and densities at each monitoring point during the 150 s simulation were collected to validate Equation (28).
The main part of the calculation is the calculation of thermal resistance, which can be expressed by the Formula (34).
The interval section is divided into monitoring points based on the modelling data of
L = 100,
in + 1 =
in + 5. To assess the accuracy of this equation, the thermal resistance decrease gradients are calculated for each interval. The results are illustrated in
Figure 14.
As illustrated in
Figure 14, the numerical simulation results and the thermal resistance trend graph of the constructed mathematical model show a decreasing trend with distance. Upon verification, it was found that the maximum error in the calculated results is within 0.03 pa, representing an error rate of approximately 2.3%. The fit between the results and the numerical simulation is as high as 97.7%, indicating a high level of confidence in the obtained results.