3.1. Back-Layering Length and Critical Velocity
Back-layering length is the important parameter in roadway fires. The back-layering flow is the most fatal contamination to workers who are blocked in upstream of the fire source.
Figure 3 shows the relationship between back-layering length and longitudinal ventilation velocity under different heat release rates.
Figure 3 indicates that the back-layering length decreases constantly as the ventilation velocity increases and presents a linear trend. Moreover, under the same ventilation velocity, the back-layering length increases as the heat release rate increases. However, when the heat release rate is more than 600kw, the back-layering length increases quite slowly. This indicates that the influence of heat release rate on the back-layering length is gradually weakened as the heat release rate increases.
Critical ventilation velocity is also a key parameter for ensuring the safety and proper emergency evacuation of workers in roadway fires. The critical ventilation velocity in single-tube tunnel fires has been researched thoroughly. Wu and Bakar [
12] established a mathematical model for dimensionless critical ventilation velocity and dimensionless heat release rate in a single-tube tunnel, as shown in Equation (2):
where
represents the dimensionless critical ventilation velocity,
represents the critical ventilation velocity (m/s),
represents the gravity acceleration (m/s
2),
represents the hydraulic tunnel height (m),
represents the dimensionless heat release rate,
,
represents the heat release rate (kw),
represents the ambient air density (kg/m
3),
represents the specific heat capacity (kJ/(kg⋅K)), and
represents the environment temperature (K).
Li et al. [
13] also acquired the correlation between dimensionless critical ventilation velocity and dimensionless heat release rate in a single-tube tunnel, which can be expressed in Equation (3):
where
and
represents the tunnel height (m).
In this paper, the critical ventilation velocity was determined by the X direction velocity vector beneath the ceiling in the main roadway. The critical ventilation velocities are 1.5 m/s, 1.8 m/s, 2 m/s, and 2.1 m/s when the heat release rates (HRR) are 300 kw, 600 kw, 900 kw, and 1200 kw, respectively. The critical ventilation velocity predicted by CFD simulation under different HRRs is shown in
Figure 4a, and the critical ventilation velocities calculated by equations proposed by Wu [
12] and Li [
13] are also presented in
Figure 4a for comparison. It can be seen that the critical ventilation velocity increases with the increase in heat release rate and the increased value decreases slowly. Also, we can find that the critical ventilation velocity is higher than Wu’s and Li’s model, which means the critical ventilation velocity in the T-shaped bifurcated roadway is higher than in the single-tube roadway when the fire is located upstream of the T-junction. This is because, when the airflow passes through the T-junction, a portion of ventilation mass flow is pushed into the branched roadway, so the actual ventilation mass flow in the main roadway of the T-shaped bifurcated roadway is lower than in the single-tube roadway.
Figure 4b presents the correlation between dimensional critical ventilation velocity and dimensional heat release rate. According to
Figure 4b,
is directly proportional to
; hence, the prediction model of dimensional critical ventilation velocity in the T-shaped bifurcated roadway can be expressed as:
It can be found that Equation (4) is similar to Equations (2) and (3) and is more close to Li’s model.
To further validate the accuracy of the numerical model adopted in this paper, the critical ventilation velocities predicted by CFD were compared with data from Li’s small-scale tests [
13]. Li’s small-scale experimental tests were conducted in a 12 m long model tunnel. The fire source, which was a 100 mm diameter porous bed burner, was located in the center of the tunnel model. We built the same simulated tunnel model as Li’s, and the numerical model was the same as mentioned above (see
Section 2.2). The fire source was set as a volumetric source with a height of 0.25 m and a diameter of 0.1 m.
Table 2 presents the critical ventilation velocity predicted by CFD under different cases in Li’s experiment, and the comparison between simulation values of critical ventilation velocity and tested results by Li is shown in
Figure 5. According to
Table 2 and
Figure 5, there is a reasonable agreement between the simulation and the experiment tests.
3.2. Longitudinal Temperature Profile in the Main and Branched Roadway
Inhaling high-temperature smoke is the main factor to cause casualties in fire accidents, so the temperature distribution at breathing zone height (Z = 1.6 m) was analyzed in this paper. The tolerance time in a fire environment is 12 min for a person [
28], so the escape temperature that a person can escape from the fire environment successfully is set to be 60 °C.
Figure 6 presents the longitudinal temperature profile at breathing zone height of the main roadway when heat release rates are 300 kw, 600 kw, and 900 kw, respectively.
Figure 6a illustrates that, when the longitudinal ventilation velocity is 1 m/s, the maximum temperature is located at the fire source and the value is 101.636 °C, and the influence range of the fire source on upstream temperature is 21.4 m, which can be explained by smoke backflow and thermal convection. When the longitudinal ventilation velocities are 1.5 m/s, 2 m/s, 2.5 m/s, and 3 m/s, respectively, the influence range of fire source on upstream temperature is 0 m and the smoke temperature at the fire source increases dramatically from 25 °C to 41.16 °C, 37.13 °C, 34.73 °C, and 33.07 °C, respectively; then, the temperature rises slightly and becomes steady along the longitudinal direction (X direction). The smoke temperatures in the main roadway I are all lower than escape temperature when the ventilation velocities are 1 m/s~3 m/s. As shown in
Figure 6b, when the longitudinal ventilation velocities are 1 m/s and 1.5 m/s, the values of maximum temperature are 179 °C and 140 °C, respectively, and the influence ranges of fire source on upstream temperature are 24.2 m and 14 m, respectively. When the longitudinal ventilation velocities are 2 m/s, 2.5 m/s, and 3 m/s, respectively, the smoke temperature at the fire source increase dramatically from 25 °C to 37.366 °C, 34.869 °C, and 33.2 °C, respectively; then, the temperature rises slightly and becomes steady along the longitudinal direction (X direction). The smoke temperatures in the main roadway I are all lower than escape temperature when the ventilation velocities are 2 m/s, 2.5 m/s, and 3 m/s. According to
Figure 6c, when the longitudinal ventilation velocities are 1 m/s and 1.5 m/s, the values of maximum temperature are 296.96 °C and 138.8 °C, respectively, and the influence ranges of fire source on upstream temperature are 27 m and 18 m, respectively. When the longitudinal ventilation velocities are 2 m/s, 2.5 m/s, and 3 m/s, respectively, the smoke temperature at the fire source increase dramatically from 25 °C to 61.90 °C, 54.35 °C, and 49.68 °C, respectively; then, the temperature rises slightly and becomes steady along the longitudinal direction (X direction). The smoke temperatures in the main roadway I are all lower than escape temperature when the ventilation velocities are 2.5 m/s and 3 m/s.
Then, the effect of longitudinal ventilation velocity and heat release rate on average temperature at breathing zone height inside the bifurcated roadway is discussed thoroughly.
Figure 7 presents the average temperature at Z = 1.6 m inside the main roadway I and branched roadway under different ventilation velocities and different heat release rates.
It can be observed from
Figure 7 that the average temperatures in the main roadway I and branched roadway are lower than 60 °C when the HRR is 300 kw and the ventilation velocities are 1 m/s, 1.5 m/s, 2 m/s, 2.5 m/s, and 3 m/s and when the HRR is 600 kw and the ventilation velocities are 2 m/s, 2.5 m/s, and 3 m/s and when the HRR is 900 kw and the ventilation velocities are 2.5 m/s and 3 m/s. The average temperatures in the branched roadway under different cases are all slightly lower than that in the main roadway I. When the heat release rate is 300 kw, the maximum temperature difference between the main roadway I and the branched roadway is 3.1 °C and the minimum temperature difference is 0.7 °C. When the heat release rate is 600 kw, the maximum temperature difference is 2.12 °C and the minimum temperature difference is 0.585 °C. When the heat release rate is 900 kw, the maximum temperature difference is 6.47 °C and the minimum temperature difference is 1.43 °C. When the heat release rate is 1200 kw, the maximum temperature difference is 8 °C and the minimum temperature difference is 0.3 °C.
Figure 7 also shows that, under different heat release rates, the average temperatures in the main roadway I and branched roadway decrease as the ventilation velocity increases. This occurs because the thermal convection between airflow and smoke is intensified with the increase in ventilation velocity. And the temperature decay fits the power function, as shown in Equation (5):
where
represents the average temperature at Z = 1.6 m (°C),
represents longitudinal ventilation velocity (m/s), and
and
represent dimensionless coefficient, respectively.
The fitting curves are shown with solid lines and dashed lines in
Figure 7. The values of
and
are displayed in
Table 3. It can be seen from
Figure 7 and
Table 3 the values of R-squares are all above 0.98; it can be concluded that the variation in average temperature in the main roadway I and branched roadway according to ventilation velocity can be accurately describe by Equation (5).
The variation in average temperature with the increasing heat release rate in the bifurcated roadway is shown in
Figure 8. It can be observed that the average temperatures in the main roadway I and branched roadway increase constantly as the heat release rate increases, and the predicted data of average temperature can be correlated to HRR with the following equation:
where
represents the average temperature at Z = 1.6 m (°C),
represents the heat release rate (kw), and
and
represent the dimensionless coefficient, respectively.
The fitting lines are shown with solid lines and dashed lines in
Figure 8. The values of
and
are displayed in
Table 4. It can be seen from
Figure 8 and
Table 4 the values of R-squares are all above 0.99, which indicates that the variation in average temperature in the main roadway I and branched roadway according to heat release rate can be accurately predicted by Equation (6).
3.3. The Profile of CO Concentration in the Main and Branched Roadway
Figure 9 presents the smoke propagation in the T-shaped roadway when the HRR is 600 kw and longitudinal ventilation velocity is 1.5 m/s. According to
Figure 9, it can be seen that a high concentration of CO mainly gathers near the top plate of the fire source and spreads upstream of the fire source at a certain distance. Under the influence of longitudinal ventilation, the CO generated by the fire gradually spreads downstream of the main roadway and inside the branched roadway and then is rapidly diluted. The CO concentration in the main roadway I and branched roadway are much lower than near the fire source.
Based on previous fire accident analyses, we know that high-temperature smoke contains CO, and casualties were partly caused by inhaling CO, so the distribution of CO concentration at breathing zone height is thoroughly analyzed and discussed in this section.
Figure 10 shows the concentration distribution of CO at breathing zone height in the main roadway I and branched roadway under different longitudinal ventilation velocities when the HRR is 600 kw.
We can see clearly from
Figure 10 that the concentration of CO at the height of the breathing zone in the main roadway I and branched roadway gradually decreases and tends to be of uniform distribution as the longitudinal ventilation velocity increases. When the ventilation velocities are 1 m/s and 1.5 m/s, the concentration of CO at the height of the breathing zone in the main roadway I is higher than in the branched roadway. And, when the ventilation velocities are 2 m/s, 2.5 m/s, and 3 m/s, respectively, the concentration of CO in the branched roadway is slightly higher than in the main roadway I.
The CO volume concentration in the safe evacuation passage is less than 500 ppm in the fire environment [
28], so 500 ppm is chosen to be the critical concentration to escape safely for a person. The average concentration of CO at the breathing zone height inside the main roadway I and branched roadway is shown in
Figure 11.
According to
Figure 11, it can be seen that the average concentrations of CO in the main roadway I and branched roadway are all lower than 500 PPM. When the HRR is 900 kw, the maximum concentration difference between the main roadway I and branched roadway is 84.444PPM and the minimum concentration difference is 5.676 PPM. When the HRR is 600 kw, the maximum concentration difference is 19.72 PPM and the minimum concentration difference is 4.776 PPM. When the HRR is 300 kw, the maximum concentration difference is 54.094 PPM and the minimum concentration difference is 2.686 PPM. In addition, the concentration differences of CO between the main roadway I and the branched roadway decrease as the ventilation velocity increases. And we also can find that the average concentrations of CO inside the main roadway I and branched roadway decrease constantly with the increase in longitudinal ventilation velocity, which fits the power function. The fitting curves are shown in
Figure 11 and the fitting functions are shown in
Table 5. The values of R-squares are all above 0.95 according to
Figure 10 and
Table 5; the equations describing the variation in average concentration of CO are shown in
Table 5.
The effect of HRR on average concentration of CO inside the main roadway I and the branched roadway are discussed then.
Figure 12 shows the variation in the average concentration of CO according to the heat release rate. It can be seen that, as the heat release rate increases, the average concentration of CO inside the main roadway I and branched roadway increase linearly. And it can be described by the linear function; the fitting lines are shown in
Figure 12 and the equations are presented in
Table 6.
Table 6 illustrates that the value of R-squares is 0.89406 when the ventilation velocity is 1 m/s and the type of roadway is branched roadway. However, the values of R-squares are all above 0.94 in other cases. It can be explained that, when the ventilation velocity is 1 m/s, the volume of fresh airflow pushed into the branched roadway is much lower than in the main roadway I, so the mixing of smoke and fresh airflow is uneven, which leads to the disordered distribution of smoke in the branched roadway. Therefore, it can be concluded that the average concentration of CO inside the main roadway I and branched roadway can be predicted better by this kind of equation.