In this work, the study of smoke fire transportation inside of a subway passenger platform is presented. The study includes a set of numerical simulations to observe the behavior of the smoke inside the platform. Two smoke transport simulations using the FDS program are also included. Subsequently, the development of a 3D - 1:100 scale model is described and it was used to perform an experimental observation of the phenomenon. The model was built by using a 3D printer which allowed to include more architectural details of the real scenario. The inclusion of these details allowed to observe qualitative similarity between the results of the simulation and the experimental work. Although there are clear differences between what could happen in a real scenario and what was observed in the scale model, it was identified that the model is an important complement to the simulations. In addition to the simulations, the use of this type of 3D models allows the observation of the phenomenon by different specialists such as firefighters, policeman, medical personnel, etc., in the same place and its intention is to provide a more interactive tool to the observation group, increasing the time devoted to the development of contingency actions and reducing the costs associated with the logistics of a real simulacrum. The model allows to better identify the strengths, opportunities, weaknesses and threats of the contingency procedures developed by the safety and hygiene groups and to make their corresponding adjustments if necessary.
Keywords: CFD, FDS, ANSYS-FLUENT, smoke, critical velocity, 3D print
During a fire emergency in a subway platform, several specialized teams have to interact each other to ensure the safety of the passengers. At least 25 tasks and resources have been identified during an empirical study of emergency response . They have to deal with a situation in which, the architectural distribution creates important differences in the behavior of the fire and smoke distribution  loss of visibility [3, 4], as well as the ventilation characteristics . The reaction procedure of the staff is intended to adequately manage the time that elapses until the last person is evacuated [6, 7]. Passenger reactions to an emergency in subways have shown to be independent of their individual characteristics such as gender, educational level and for example; carrying-on luggage . These kind of disasters involve a great human cost because in subterranean emergency scenarios, the possibilities for the passengers to leave or scape out of a critical location is quite compromised. The limited space for mobility, low visibility, presence of fresh air, etc., are some of the restrictions that a passenger can face in these situations. Therefore we propose the use of computer simulations together with observations in 3D scale models to increase the scope in the study of these phenomenon.
When designing subterranean installations, connections to atmospheric pressure are one of the most important issues that have to be considered. The location and size of the ventilation ports, external and internal weather conditions among others create set of variable air conditions. In the case of subway tunnels, convective flow currents enhance the generation of vortex in random locations even when the air velocity can be maintained constant along the tunnel by natural conditions or by means of a mechanical ventilation system. The emission of smoke close to a ventilation duct connected to the atmosphere does not always guarantee a safety condition because the atmospheric pressure creates an equilibrium between the internal and external pressures affected by the temperature and the resulting direction of the air flow in consequence is not easy to predict.
This study is composed by several steps with the purpose of observing in a separate way the events occurring in a fire scenario. In these events, several physical situations are presented simultaneously, by dividing the sequence of them, interesting observations can be made to understand the complexity of the phenomenon when the equilibrium condition is reached. The temperature changes generated inside the platform will have a straight interaction with the ducts connected to atmospheric pressure, and finally addressing the transportation of the smoke. Therefore, the study was divided in three stages. The first stage was developed to validate the critical velocity calculation. The second one to simulate the smoke transportation phenomena and the third one describes the creation of the scaled model to observe and compare the displacement of the smoke with respect to the simulation.
The first stage considers the analysis of the flow velocity inside the platform when the critical velocity is injected from one side of the tunnel. There is not fire source in this case but only the air flow velocity of the calculation. The NFPA 130-7 provides important information to carry out the calculation of the critical velocity which is the velocity required to provide adequate ventilation conditions and to produce adequate smoke transportation in case of fire emergencies . This critical velocity avoids also the backlayering effect in which the smoke in certain cases tends to move in a different direction of the one desired. To observe the flow pattern, a simplified theoretical 2D model was implemented in FORTRAN. This flow pattern was compared with the results obtained in ANSYS-FLUENT and in FDS programs by extracting only a sliced plane of each one. The observation of these results are very important during the selection and location process of commercial fans for subterranean installations.
The second stage considers a 180 seconds fire simulation in FDS . In this case, a totally developed fire of 5MW is located close to the center of the platform as in . The evacuation of the smoke trough the natural ventilation ports is another interesting issue to observe in the simulation  and to notice the presence or not of the backlayering effect with the injection of the critical velocity [13-15]. It is not easy to predict the displacement of the smoke because at least 4 ventilation ducts connected to atmospheric pressure are close to the fire source. In one simulation, the air conditions are allowed to move freely. In a second simulation, the previously calculated critical velocity is injected from the left side of the platform and both results are compared.
In the third stage, a 3D printed model was constructed to create a suitable observation experiment setup by using a Mendelmax printer and Polylactic Acid (PLA) as the deposition material with 0.2mm thickness layers. The stage has been developed with the purpose of re-creating the scenario to observe the smoke spread. Once the smoke output is constant, an exchange of temperatures within the platform is generated, these differences between the hot and cold gases, after some time, reach an equilibrium behavior creating a flow distribution (under free motion conditions). This equilibrium is directly related with the distance between the ventilation ducts because they are communicated through the platform and are directly affected by the atmospheric pressure. Once the ventilation system has been turned on, it will be possible to observe if the behavior is similar to the one observed in the simulations (stage 2), although in a real scenario many other factors will be presented and may differ from this experiment. Table 1, shows a list of all the analyzed cases in each stage of this study.
|Case No.||Program||Simulation type||Objective|
|1||Fortran||Air flow injection at critical velocity||Validation|
|2||ANSYS-FLUENT||Air flow injection at critical velocity||Validation|
|3||FDS||Air flow injection at critical velocity||Validation|
|4||FDS||Smoke transported at free motion condition||Observe and measure the free smoke movement|
|5||FDS||Smoke transported with the air flow at critical velocity||Observe and measure the influenced movement of the smoke|
|6||Experimental setup||Smoke transported in free motion||Observe experimental behavior with respect to the FDS simulations|
|7||Experimental setup||Smoke transported with the fans on||Observe experimental behavior with respect to the FDS simulations|
In the presence of fire inside underground facilities, different parameters should be considered to handle the emergency situation and provide a suitable breading atmosphere for the users, it must be kept sustainable conditions for life. However, every subway systems have particular characteristics, they require specific solutions which relate internal and external conditions of the tunnel as an enclosed place with limited escape routes.
The analyzed platform considered has 15m depth from the sidewalk level and 150m in length, it has four ventilation ports (ducts) at the ceiling connected to atmospheric pressure and three passenger access towards the street. Figure 1, shows the dimensions and geometry considered for the 3D model. They were obtained by direct measurements of the platform creating a simplified geometry.
During the backlayering effect, the ventilation conditions are not suitable to transport the smoke to a desired direction or exit. The critical velocity has become in one of the prime criteria for the designing of ventilation systems, especially in tunnel section . In this study, the Kennedy model (equations 1 up to 3) were used to calculate the critical velocity in the subterranean scenario, where the value of the critical velocity is determined as a solution of the three equations of the Kennedy model [16-19].
For the model, K1= 0.605 is a dimensionless constant value given by equation 3, the Froude number considered for the calculation takes a value of 4.5 according to [20-22] and K2=1 is a dimensionless slope factor; the average burning temperature was TF=369.34 °K (96.2°C) which resulting values corresponds to the ones presented in . The 5MW fire source was considered according with [24, 25], and the remaining values considered for the critical velocity calculation are shown in Table 2 .
|Gravity acceleration (g)||9.81|
|Tunnel height (H)||5.2 m|
|Air density ( )||1.176|
|Specific heat of air (Cp)||1.007|
|Transverse section area (A)||41.6|
|Air temperature ( )||300.15 °K (97 °C)|
|Critical velocity (Vc)||1.464|
The present analysis is based on the pressure-velocity coupling method, known as vorticity-stream function. The governing equations for the analysis consider an incompressible and Newtonian flow (2D), the mass conservation and momentum balance are given in equations 4 and 5. The analysis considers a procedure step by step with respect to time, variable velocities (), density () and dynamic viscosity (μ):
The vorticity - stream function ( ) equations were discretized using the finite differences method where each node is used to describe the computer dominion (equations 6 and 7) :
Here, U and V corresponds to the velocities in the x and y directions (equation 8):
After discretization of equation 6, the following is obtained:
In which, the coefficients are:
And for equation 7:
The coefficients are:
The model was constructed by considering a rectangular shape with 50 divisions in width and 275 divisions in length representing the top view of the passenger platform, the discretization was defined by considering a regular mesh with an equivalent and of 0.2m. This value was obtained by considering a physical real dimension of 10.5m for the width and 58m for the length of the model and used for the analysis of case 1. On the other hand, the resulting Reynolds number was in this case 4.898 x 10-5.
Figure 2 a) shows the distribution of boundary conditions in the 2D model (case 1). The resulting mesh arrangement is shown in figure 2 b). The figure represents the top view of an idealized duct; the control area, it is represented by a horizontal layer crossing along the center of the passenger platform. The location of each boundary condition is applied on each side by considering the Thom formulation and an idealized velocity profile .
This analysis (case 2) was carried out only to observe the flow stream along the platform at the critical velocity and compare with case 1. It considers the Finite Volume formulation allowing the discretization and solution of the differential equations. The resulting 3D model is shown in figure 3, it is composed by a large rectangular shape representing the passenger platform, and four rectangular blocks were used to represent the natural ventilation ports. Another three slender blocks (tilted) were used to represent the stairs. The location of the areas connected to atmospheric pressure were identified as condition type 1. These includes the natural ventilation ports and passenger accesses. For the type 2 condition, 1.46m/s of air velocity was applied at the end of the tunnel (left side).
Different mesh methods were tested until an excellent mesh quality range (qualification) was obtained according to . The resulting model contains 284,816 nodes and 247,974 brick type elements as shown in figure 4.
Table 3 shows the main parameters considered for the ANSYS-FLUENT simulation.
|Solution type||Pressure based|
|Viscous model||k-ε standard|
|Air density||1.164 kg/m3|
|Dynamic viscosity||1.86 e-05 kg/m.s|
|Inlet velocity||1.464 m/s|
Fire Dynamic Simulator (FDS) was used to simulate the behavior of the calculated critical velocity in the platform [30-32]. The CAD modelling for this case was created considering the same dimensions of the previous case and the resulting model is shown in figure 5. In this model, three case studies were considered, one was used to verify the development of the critical velocity named case 3. Case 4, includes the fire source in free motion and case 5 includes the fire source and the application of the critical velocity (location 2).
For the FDS meshing, the total number of rectangular cells was 2,474,693, a simplified view of the resulting brick mesh is shown in figure 6 a) isometric view and a close view in figure 6 b).
|Figure 6. Resulting mesh in FDS. (a) Isometric view of the passenger station. (b) Close view of mesh layers at the center of the passenger platform|
The gathering of results was made by considering a horizontally oriented rectangular plane crossing along the model, located at 3.4m (height) from the floor level, in the middle of the platform, this plane was used to observe the flow contour plot similar to the one presented in  but in this case horizontally oriented.
Figure 7 a), shows a dimensionless scale of the resulting air velocity distribution. This results were scaled to the real dimensions and tracked through a crossing line at the middle of the plane as shown in figure 7 b).
Figure 8, shows the resulting velocity contour obtained in ANSYS-FLUENT. The red frame represents the portion of the model analyzed in FORTRAN. An increment of the velocity is present at both ending sides of the model due to the reduction of the cross section area (-70m and 80m) between the tunnel and the passenger platform.
The results observed in figure 8 shows that the maximum velocity is achieved at approximate 10m away from the location of the air injection. After that, the air reduces its velocity specially when achieving the passenger entrance area which is also in contact with the exit stairs. When the air crosses trough the natural ventilation ducts this velocity is also reduced. An interesting effect is then observed following this direction path. The air velocity suddenly increases when entering the tunnel at approximate 5m. This can be produced by the abrupt reduction of the cross section area between the passenger platform and the tunnel. This phenomenon indicates that special considerations have to be taken when stablishing the switching on and off of the mechanical ventilation fans because the architectural details of the construction directly affect the behavior of the flow.
Figure 9, shows the velocity contour obtained in FDS for case 3. The results present a similar behavior as in ANSYS-FLUENT.
Case 4 corresponds to the simulation of the smoke evolution in free motion condition. The smoke moves only according to the influence of the temperature changes and atmospheric interactions. Case 5 presents the evolution of the smoke affected with an air velocity of 1.46 m/s. The evolution of the simulations is presented in figure 10. From a) to c) corresponds to case 4 and figures d) to f) corresponds to case 5.
|Evolution of smoke for Case 4||Evolution of smoke for Case 5|
|(a) 30 sec.||(d) 30 sec.|
|(b) 90 sec.||(e) 90 sec.|
|(c) 180 sec.||(f) 180 sec.|
|Figure 10. Evolution of smoke transportation for cases 4 and 5.|
According to the calculation described in , the resulting flame height is in this case 6.07m which is higher than the top ceiling level as observed in figure 11, and a safety distance for the passengers to be away from the fire source goes beyond 10m [35, 36].
The difference of heat distribution for cases 4 and 5 after 3 minutes of fire evolution is presented in figure 12, a) and b) for free motion, and c) and d) with the fans switched on. When the critical velocity is injected from the left side of the platform.
|Temperature distribution for case 4||Temperature distribution for case 5|
|(b) 180 sec.||(d) 180 sec.|
The distribution of the heat is maintained above 35°C for more than 80m long in the free motion case with a maximum of 63°C. When the critical velocity is acting on the scenario, a notorious local increment in the temperature reaches more than 220°C. After approximately 10 meters this temperature drops to 100°C and 40meters away it falls down below 60°C. As in many fire situations, an increment in the air velocity can produce enrichment of the fire reaction.
The use of the experimental setup is important because it allows to include more physical senses during the study of the phenomenon. Although the creation of the scaled model is not a straight way work, instead it has presented a series of challenges, because of that, a series of recommendations are described below.
Figure 13 shows the distribution of the experimental setup, the scaled model is mounted on a support table. Apart from the 3D printing of the model, extra support equipment will be required to locate the camera systems depending on the shots or pictures that are required to obtain. Another important factor that must be carefully considered is the treatment of the filming (or photographs) in scaled scenarios. For this reason, an adequate lighting system will allow to obtain better pictures during the experiment. In this case, different attempts had to be made to obtain a suitable observation sequences.
For the creation of the scaled model, a series of modifications were considered in the original 3D CAD model in order to create the printable file. The modifications included the elimination of the ceiling and selected walls also. These layers will be next substituted with transparent glass to allow the observation of the smoke. The final CAD file is shown in figure 14.
It is possible to observe also that the model was divided in many sections to be finally assembled together. A series of mechanical supports were also included to provide rigid support and to hold in position the transparent layers. Previous prints were made to obtain the model shown in figure 15.
In this case, this model allow us to validate the assembling process between all the sections and make corrections if needed. For convenience of this picture, the bottom stairs were located on top of the model upside down (number 11 in the picture). Once the entire model was validated, a 1:100 scaled model was finally printed. It is recommended to create a numbered map and list of all the pieces to avoid confusion during the assembling. In this form, the model can be easily transported to a preferred location for the experiments. Figure 16 shows the resulting assembled model, now the bottom stairs are properly located at the bottom face of the model.
Finally, the transparent layers of the ceiling were installed, figure 17 shows the final scenario constructed for the experiment. Prior to the experiment it was verified the room temperature (21°C), that there were no smoke leaks around the model and no external air currents also. An LCD digital gauge meter measure anemometer (GM8908) together with a regulated power supply connected to the fan motor were used to control the velocity of the air injected to the scaled model. After several trails, a velocity of 0.245 m/s was achieved. This value was calculated by using the critical velocity equation but replacing the corresponding scaled values (area, height, fire power, etc.).
Besides of the top views captured during the experiments, a close lateral view in the location of the smoke source is also considered and presented in figure 18. In this case, the ceiling reflect is observed on the top side of the picture. This is one of the issues that can appear during the close-up shots and have to be identified during the interpretation of results.
Once the smoke has been injected into the setup, a free displacement of smoke was allowed. After 60 seconds, the internal equilibrium of the system was clearly achieved generating a flow pattern like the one shown in figure 19. In the right side of the figure, the tendency of the smoke to move out of the platform in the same ventilation duct as the one observed in the FDS simulation can be observed.
In figure 20 a) a top view of the smoke movement was captured when the fans were switched on. This is compared with the result obtained in the FDS simulation in figure 20 b). Similar displacement of the smoke was observed in both scenarios showing that the smoke leaves out in the same ventilation duct despite of the differences appreciated between the scaled model and the simulation.
Figure 21 shows a lateral close view in the smoke emission zone in two different cases, a) when the fans are turned off and b) when the fans are turned on.
Several differences with respect a real scenario can be observed also. For example, in a real fire source, the temperature of hot gases can maintain the smoke plume close to the ceiling for a longer distance than the one observed in picture 21 a). The smoke displacement on the bottom stairs is another important difference. According to the fire source location, there exist the possibility of the passengers to use this bottom stairs to evacuate towards the opposite side of the platform. In the experiment the presence of smoke was observed there creating a dangerous situation in which passengers taking the bottom stairs could be facing the presence of smoke on top and bottom sides of this particular route as observed in the simulations.
In figure 10 a) which corresponds to 30 seconds in the simulation, the smoke moves enough to reach the closest ventilation duct and leave the platform, which corresponds similarly in the experimental setup. There, the smoke starts to leave the platform in the same duct after 40 seconds from the beginning of the emission. It was observed that, even when the smoke reaches the same longitudinal distance in a similar time, it took at least 10 more seconds to the smoke of the experiment to ascend through the ventilation ducts, due to the wall roughness and size relation between simulation and experiment and this an important drawback to consider. Likewise, the reduced dimensions of the setup have an important influence in the temperature of the smoke. In the simulation, the smoke travels close to the ceiling of the platform and in the experimental setup, the smoke travels close to the floor which can be observed by comparing figure 11 with figure 21 a). After 180 seconds of simulation, most of the smoke presents stagnation in the platform. In the experimental work, this was also observed even before the 180 seconds. The stagnation continued with the same behaviour for more than 5 minutes. The amount of smoke leaving the platform in both analyses of free motion is not enough to identify a smoke-free route towards the exit stairs. When the ventilation fans were switched on in the simulation, no flashback was observed but in the experimental setup, this behaviour is difficult to observe because the amount of air enclosed in the model is very sensitive to the smallest changes of air movements. The amount of smoke being emitted and the velocity of air injected to the model are two main sensitive variables not easy to control and they can create totally different flow patterns inside the platform.
The development of the 3D printed setup combined with computer simulations, allows to have better interaction with the phenomena but also with the development of emergency procedures. The implementation of a contingency plan is complex and expensive because it requires the time of coordinated interaction with different specialized personnel. For example, one of the initial steps during an emergency contingency considers the communication procedure between the security personnel and to the passengers advising the need of evacuation. It requires the interaction of internal and external groups affecting not only the subterranean service but also (in the case of Mexico City) the activation of the external emergency transport service, which can face issues related to the urban traffic or sudden obstructions of the streets . Another important step in emergency contingency is related to the evacuation process. It would be easy to generate external obstructions in the streets near the stations during the evacuation process. External groups like fire fighters, police patrols, ambulances, etc., require in just that moment a free access to the stations. Mistakes in this procedures can conduct to several injuries or demands to the transport service. The setup helps also to create a better understanding of the subterranean distribution of the building. It reduces the chance to have situations in which, the emergency personal get lost inside the subterranean installation specially in long transfer platforms or the situations in which the arrival of the crew is not the one expected for example in an opposite or different doorway of the same platform. The amount of people on each platform, agglomeration of people at peak hours, size and complexity of the station, internal and external accessibility to the platform, location of extinguishers, etc., are elements that must be discussed and studied before the appearance of the emergency situations. Better designed 3D scaled models are a suitable complement because they create an interactive tool in which several crew members can observe and discuss about common topics and to identify weaknesses, strengths, opportunities, and threats of previously designed contingency plans. A wide variety of analysis scenarios can be observed in a very short time by changing the location of the smoke source (hose location) and probably not only for fire contingency. Several future work implementations can be applied in the setup from sophisticated control systems to special illumination for pattern recognition, image analysis algorithms, thermography or other technologies for intelligent support of firefighters that can be implemented also in real platforms.
According to the simulation, 1.46m/s is required to produce adequate smoke transportation of the 5MW fire source; no backlash was observed. This was an idealized case in which the velocity of air can be altered for example by the presence of the piston effect produced by trains traveling along the tunnels. These results consider just a limited case of study, for example in real cases, an increment of the velocity caused by external factors could produce an increment in the combustion also due oxygen enrichment towards the fire source. The results presented here correspond qualitatively to the ones presented in the literature with differences related to the power of fire considered in each investigation. The time it takes the smoke to occupy the entire platform can be used to calculate the available time the passengers have to evacuate the platform, in these situations the critical velocity plays an important role for the transportation of the smoke.
Authors acknowledge the Instituto Politécnico Nacional (IPN) and the Consejo Nacional de Ciencia y Tecnología (CONACYT) México for the support provided for this academic research.
 Qiming Li, Yongliang Deng, Cong Liu, Qingtian Zeng, Ying Lu, Modeling and analysis of subway fire emergency response: An empirical study, Safety Science 84 (2016) 171–180, http://dx.doi.org/10.1016/j.ssci.2015.12.003 0925-7535
 Ran Gao, Angui Li, Ying Zhang, Na Luo, How domes improve fire safety in subway stations, Safety Science 80 (2015) 94–104, http://dx.doi.org/10.1016/j.ssci.2015.07.015 0925-7535
 Kai Kang, A smoke model and its application for smoke management in an underground mass transit station, Fire Safety Journal 42 (2007) 218–231, doi:10.1016/j.firesaf.2006.10.003
 Manabu Tsukahara, Yusuke Koshiba, Hideo Ohtani, Effectiveness of downward evacuation in a large-scale subway fire using Fire Dynamics Simulator, Tunnelling and Underground Space Technology 26 (2011) 573–581, doi:10.1016/j.tust.2011.02.002
 Na Luo, Angui Li, Ran Gao, Zhenguo Tian, Zhipei Hu, Smoke confinement utilizing the USME ventilation mode for subway station fire, Safety Science 70 (2014) 202–210, http://dx.doi.org/10.1016/j.ssci.2014.06.002 0925-7535
 Wang Zhilei, Hua Min, Xu Dayong, Pan Xuhai, Simulation research on human evacuation in subway with a single point fire scenario, Procedia Engineering 84 ( 2014 ) 595 – 602, doi: 10.1016/j.proeng.2014.10.472
 Jae Seong Roh, Hong Sun Ryou, Won Hee Park, Yong Jun Jang, CFD simulation and assessment of life safety in a subway train fire, Tunneling and Underground Space Technology 24 (2009) 447–453, doi:10.1016/j.tust.2008.12.003
 Jing-Hong Wang, Wen-yu Yan, You-ran Zhi, Jun-cheng Jiang, Investigation of the panic psychology and behaviors of evacuation crowds in subway emergencies, Procedia Engineering 135 ( 2016 ) 128 – 137, doi: 10.1016/j.proeng.2016.01.091
 NFPA 130, Standard for Fixed Guideway Transit & Passenger Rail Systems, (2007 edition)
 Jae Seong Roh, Hong Sun Ryou, Won Hee Park, Yong Jun Jang, CFD simulation and assessment of life safety in a subway train fire, Tunnelling and Underground Space Technology 24 (2009) 447–453, doi:10.1016/j.tust.2008.12.003
 Dong-Ho Rie, Myung-Whan Hwang, Seong-Jung Kim, Sung-Wook Yoon, Jae-Woong Ko, Ha-Yong Kim, A study of optimal vent mode for the smoke control of subway station fire, Tunnelling and Underground Space Technology 21 (2006) 300–301, doi:10.1016/j.tust.2005.12.157
 Shaogang Zhang, Yongzheng Yao, Kai Zhu, Kaiyuan Li, Ruifang Zhang, Song Lu, Xudong Cheng, Prediction of smoke back-layering length under different longitudinal ventilations in the subway tunnel with metro train, Tunnelling and underground space technology 53 (2016) 13-21, http://dx.doi.org/10.1016/j.tust.2015.12.013 0886-7798.
 Shaogang Zhang, Xudong Cheng, Yongzheng Yao, Kai Zhu, Kaiyuan Li, Song Lu, Ruifang Zhang, Heping Zhang, An experimental investigation on blockage effect of metro train on the smoke back-layering in subway tunnel fires, Applied Thermal Engineering 99 (2016) 214–223, http://dx.doi.org/10.1016/j.applthermaleng.2015.12.085 1359-4311
 Miao-cheng Weng, Xin-ling Lu, Fang Liu, Xiang-peng Shi, Long-xing Yu, Prediction of backlayering length and critical velocity in metro tunnel fires, Tunnelling and Underground Space Technology 47 (2015) 64–72, http://dx.doi.org/10.1016/j.tust.2014.12.010 0886-7798
 Y.Z. Li, B. Lei, H. Ingason, Study of Critical velocity and backlayering length in longitudinally ventilated tunnel fires, Fire Safety Journal 45 (2010) 361–370
 Ahmed Kashef, Zhongyuan Yuan, Bo Lei, Ceiling temperature distribution and smoke diffusion in tunnel fires with natural ventilation, Fire Safety Journal 62 (2013) 249–255, http://dx.doi.org/10.1016/j.firesaf.2013.09.019
 W. K. Chow, On smoke control for tunnels by longitudinal ventilation, Fire and life safety, Tunnelling and Underground Space Technology, vol. 13 (1998) 3, 271–275, http://dx.doi.org/10.1016/S0886-7798(98) 00061-3
 Y. Wu, M.Z.A. Bakar, Control of smoke in tunnel fires using longitudinal ventilation systems - a study of the critical velocity, Fire Safety Journal 35 (2000) 363-390, http://dx.doi.org/10.1016/S0379-7112(00)00031-X
 Falin Chen, Smoke propagation in road tunnels, Applied Mechanics Reviews, 53 (2000) 8, 207-218, doi:10.1115/1.3097350
 W.D. Kennedy, J.A. Gonzalez, J.G. Sanchez, Derivation and application of the SES critical velocity equations, Transactions-American Society of Heating Refrigerating and Air Conditioning Engineers (ASHRAE), vol. 102, pt. 2, (1996), USA ISSN 0001-2505
 B. Giachetti, D. Couton, F. Plourde, Smoke spreading analyses in a subway fire scale model, Tunnelling and Underground Space Technology 70 (2017) 233-239, http://dx.doi.org/10.1016/j.tust.2017.08.008
 F. Tang, L. J. Li, M. S. Dong, Q. Wang, F. Z. Mei, L. H. Hu, Characterization of the bouyant flow stratification behaviours by Richardson (Froude) number in a tunnel fire with complex combination of longitudinal ventilation and ceiling extraction, Applied Thermal Engineering 110 (2017) 1021-1028, http://dx.doi.org/10.1016/j.applthermaleng.2016.08.224
 J. H. Klote, J. A. Milke, P. G. Turnbull, A. Kashef, M. J. Ferreira, Handbook of Smoke Control Engineering. Atlanta. ASHRAE, ICC, SFPE & NFPA, (2012) USA, ISBN 978-1-936504-24-4
 E. K. Stefopoulos, D. G. Damigos, Design of emergency ventilation system for an underground storage facility, Tunnelling and Underground Space Technology, 22 (3), 293–302, (2007), http://dx.doi.org/10.1016/j.tust.2006.07.002
 Morgan J. Hurley, SFPE Handbook of Fire Protection Engineering, 5th edition, Springer (2016) USA, ISBN-13: 978-1493925643
 NFPA 502, Standard for Road Tunnels, Bridges, and Other Limited Access Highways, (2011 edition).
 Suhas V. Patankar, Numerical Heat Transfer and Fluid Flow. Series in computational methods in mechanics and thermal sciences, McGraw-Hill, (1980) USA, ISBN: 0-07-048740-5
 A. Thom, The flow past circular cylinders at low speeds, Proc. Roy. Soc. London A 141, 651-669, (1933).
 ANSYS FLUENT Meshing User's Guide, ANSYS Inc., Release 15.0 (2013) USA, http://www.ansys.com
 K. McGrattan, B. Klein, S. Hostikka, J. Floyd, Fire Dynamics Simulator User’s Guide (Version 5), NIST Special Publication 1019-5, (2007) USA
 W. Lei, A. Li, C. Tai, The effect of different fire source locations on the environment of a subway station, Procedia engineering 205 (2017) 3721-3726, 10.1016/j.proeng.2017.10.305
 W. Lei, A. Li, C. Tai, The effect of heat release rate on the environment of a subway station, Procedia engineering 205 (2017) 3717-3720, 10.1016/j.proeng.2017.10.301
 You HZ, Faeth GM, An investigation of fire impingement on a horizontal ceiling, NBS-GCR-81-304, National Bureau of standards, Washington D. C., USA (1981)
 Wang Binbin, Comparative Research on FLUENT and FDS’s Numerical Simulation of Smoke Spread in Subway Platform Fire, Procedia Engineering 26 (2011) 1065 – 1075, doi:10.1016/j.proeng.2011.11.2275
 J. H. Klote, J. A. Turnbull, P. G. Kashef, M. J. Ferreira, Handbook of Smoke Control Engineering, Atlanta ASHRAE, (2012)
 Deng Jun, Ma Li, Wang Zhen-ping, Xing Zhen, Wang Wei-feng, Simulation Study on Critical Velocity of Longitudinal Ventilation Tunnel Fire, Procedia Engineering 52 (2013) 67 – 71, doi: 10.1016/j.proeng.2013.02.107
 L. A. Flores-Herrera, J. M. Sandoval-Pineda, U. S. Silva-Rivera, P. A. Tamayo-Meza, R. Rivera-Blas, CFD Simulation of obstructed ventilation ports in a subway tunnel section, International journal of simulation modeling, 16, (2017) 3, 386-398. ISSN 1726-4529, https://doi.org/10.2507/IJSIMM16(3)2.380