## 1 Introduction
[PERSON], [PERSON] [19] employed the particle method to calculate the hydrodynamic forces acting on model legs positioned on stairs. [PERSON], [PERSON] [20] put forth an adaptive transmission method to more accurately simulate underground flooding during two-layer connections, while circumventing alterations in grid size due to local model details. A variety of techniques have been utilised by scholars to simulate the process of flood inundation into subway stations, as previously investigated in the literature [21; 22; 23; 24; 25]. In their study, [PERSON], [PERSON] [24] employed a combination of the VOF and \(k\)-\(\varepsilon\) models to simulate turbulent flows, resulting in precise simulations of the free surface, velocity, and pressure on complex step spillways. The simulation outcomes were found to be in close alignment with the measured data, thereby substantiating the efficacy of turbulent numerical simulation for overflow scenarios in intricate step spillways. [PERSON], [PERSON] [21] integrated the renormalization group (RNG) \(k\)-\(\varepsilon\) model with VOF, demonstrating superior simulation capabilities in comparison to the standard k-\(\varepsilon\) model in capturing the flow dynamics of stepped spillways. [PERSON] et al. [23] employed the VOF model in conjunction with the realizable \(k\)-\(\varepsilon\), shear stress transport \(k\)-\(\omega\), v\({}^{2}\)-\(f\), and large eddy simulation (LES) turbulence models, and concluded that the realizable k-\(\varepsilon\) model is the most effective method for simulating the flow of stepped spillways. [PERSON], [PERSON] [22] conducted an analysis of the flow velocity profile of a stepped spillway using a range of turbulence models, including the standard \(k\)-\(\varepsilon\), realizable \(k\)-\(\varepsilon\), RNG \(k\)-\(\varepsilon\), standard \(k\)-\(\omega\), and shear stress transport \(k\)-\(\omega\) models. The results of this analysis demonstrated that all of the aforementioned models yielded satisfactory simulation results. Several studies [21; 22; 23; 25] have highlighted the efficacy of the RNG \(k\)-\(\varepsilon\) model in accurately simulating turbulence in flood scenarios within subway stations. The studies referenced above have demonstrated that numerical simulation has facilitated a deeper comprehension of the flood intrusion process in subterranean environments. The research primarily concentrates on the flooding trend in each subway, the various means of egress for individuals within each subway station, and the limitations of the designated evacuation level. Nevertheless, the current body of research still exhibits certain limitations and gives insufficient consideration to the following areas. (1) The path of floodwater ingress into subway stations; (2) the total volume of water accumulating within stations; and (3) the variation in floodwater flow rates into tunnels. This article presents a case study of Xinxiu Park Station in Nanning City, China. The VOF gas-liquid two-phase flow method was employed in conjunction with the \(k\)-\(\varepsilon\) model to simulate the spatial distribution of floods at varying time intervals within the subway station. Subsequently, the total accumulated water volume within the station was quantified under simulated operational conditions. Furthermore, we analysed the temporal changes in floodwater mass flow rates entering the tunnel and proposed appropriate drainage measures.