Abstract. To reduce the computational time for modelling large tunnel networks, we used FDS6.1 to investigate the solver’s capabilities for multi-scale modelling by combining 1D and full CFD models in the same domain. This study is to investigate the capabilities of FDS6.1 for simulating a tunnel ventilation system (cold and fire) using the multi-scale method. Firstly, we replicated the physical condition for cold flow using the multi-scale method and we validated the predictions against field measurements in the Dartford Tunnel, UK. This study demonstrated both the multi-scale and full models have good correlation with each other, and with the field measurements. Depending on the ratio of 1D and full CFD plus the scenarios modelled, the reduction in the computational time for the multi-scale method varies from 20% to 220% compared to the full CFD method. Secondly, we introduced a fire up to 75 MW in the multi-scale and full CFD models. Both models exhibit oscillatory mass flows along the tunnel similar to the pulsation effect observed in a tunnel fire. In the multi-scale method, the models significantly under or over-predicted the mass flows.
With exponential increase in population growth and development of ever larger densely populated cities, tunnels are often the only viable option to expand the transport infrastructure connecting towns and cities. Construction of modern road and rail tunnels increasingly becomes longer, larger and intricately connected. For example, Laerdal Tunnel is currently the world's longest tunnel at 24.51 km.
With an increasing number of long tunnels used for transport infrastructure, these tunnels need to be adequately ventilated to ensure tunnel users will not be negatively impacted by fumes from vehicles or by smoke in a fire. Although there are various options to provide ventilation in a tunnel, these methods can be broadly categorises into longitudinal, transverse or semi-transverse ventilation. In design and engineering, the ventilation system in tunnels are currently designed using both 1D and full CFD modelling methods.
1D is used at the early stages of a project where time and resources are constrained, whereas full CFD modelling is used at the later stages to validate the designs informed by the 1D modelling method. Although full CFD modelling provides detailed results, this method is significantly constrained by the computational time which can take from days to weeks. In addition, to ensure the computational time is kept to a realistic time frame for the design and engineering process, only a section instead of the larger tunnel network is modelled. This shows the current limitations for the practical application of the full CFD method.
With the multi-scale method, this is a hybrid approach combining CFD for complex regions and 1D for far-field regions of a tunnel section where an area-average representation of the variables is acceptable. This approach allows designers to simulate a significantly longer tunnel section that is impractical using only either the full CFD, or the 1D method due to the constraints of computational time.
In this research, we have used FDS6.1 [3] and the built-in HVAC feature to implement the multi-scale modelling method. The tunnel model in this research has been based on the Dartford West Tunnel and the onsite cold flow ventilation data (velocity) collected by Collela et al [2]. The onsite cold flow data is used to determine whether the multi-scale modelling method has a good correlation with both the field data and the predictions from the full CFD method. This is important to demonstrate the predictions obtained from the multi-scale modelling method correlates with real world observation. The computational times from modelling the cold flow ventilation of the Dartford West Tunnel have also been used to assess the reduction in computational time between the full CFD method and the multi-scale modelling method.
Following on that, we have used the same tunnel models but this time included a fire up to 75 MW in both the full CFD and multi-scale models to consider the feasibility of modelling a tunnel fire using the multi-scale modelling method.
For the full CFD models, a grid sensitivity study and a non-dimensionalised validation study have also been considered. We have benchmarked the full CFD models against the ArupFire Tunnel Fire Experiments documented in the FDS Validation Guide [4]. The results of the validation study has not been presented here.
Note that we have previously published the cold flow validation study in a separate peer reviewed journal [1]. The results from the multi-scale fire modelling are first published here.
The Dartford Tunnel is 1.5 km long with an effective area of 41 m2. There are 14 jet fans pairs in total where 7 is located at the South portal (100 m) and 7 is located at the North portal (50 m). These jet fans pairs are spaced 50 m apart in series, and the individual fans is spaced at 1.2 m in parallel and are mounted at 5.5 m above the road. Each fan has a diameter of 0.5 m2 and generates a volume of rate of 8.9 m3/s. See Figure 1.
We have used FDS6.1 (current at the time the work was carried out) and the built in HVAC feature to implement the multi-scale modelling method. Briefly, multi-scale modelling can be divided into two methods, namely the direct and indirect method [1].
A key consideration for a multi-scale model is the length of the CFD domain (near field), i.e. regions where CFD is needed to resolve the flow's characteristic. Once this is decided, the remaining areas are treated as far field regions where these can be represented in 1D.
For a cold flow only scenario, the CFD domain is the region where the jet fans are located, with the far field domain modelled in 1D. When a fire is introduced, the fire will also be located in a CFD domain given its complex reaction. By modelling the jet fans and the fires in CFD regions together with the regions in between represented with 1D, this essentially couples together the jet fans and the fires, allowing the effect of this coupling to be simulated.
See Figure 2 for the implementation of the multi-scale model. Note that LP is 50 m, LJF,DW is 130 m and LJF,UP is 35 m. We have derived these distances by calibrating the CFD models.
To compare the numerical model for cold flow multi-scale modelling against field data, we have carried out modelling for five scenarios including one with only background velocity (no jet fans active), and models with one, four, seven and fourteen jet fans pairs activated [1].
Figure 3 shows the comparison of the on-site and both 1D and full CFD modelling average velocities at the centre of the tunnel. Although the velocities at 20 m to 60 m have weak correlations to the field data, the velocity profiles further downstream of the jet fans correlate with the field data. This is within expectation as the numerical modelling of the jet fans are highly sensitive, particularly when considering the velocities profile near to the jet fans (highly turbulent and sensitive to the configuration of the jet fans). As the results show, as the flow from the jet fans becomes fully developed, the velocity profiles are within the field measurements.