YU Wangyang1,2, MA Jianming3, YIN Yueming4, YU Haijun1,2, WU Binbin1,2, MU Jie1,2
1. China Institute of Water Resources and Hydropower Research, Beijing 100038; 2. Research Center on Flood & Drought Disaster Reduction of the Ministry of Water Resources, Beijing 100038; 3. National Institute of Natural Hazards, Beijing 100085; 4. Sichuan Water Resources Planning and Research Institute, Chengdu 610072
Abstract:In order to simulate the flood propagation and inundation, a numerical model was established adapting the finite volume method scheme based on adaptive grid. The model takes the 2D shallow water equation as the governing equation, uses HLL scheme to approximately solve the Riemann solution to obtain the interface flux, and uses the WAF scheme to make weighted corrections to the calculated flux, so as to ensure the second order accuracy in time and space, satisfy the TVD condition, and reduce the numerical oscillation. The bottom slope source term is reconstructed by water depth difference, and the friction source term is semi-implicit solved, which ensures the harmony and stability of numerical calculation. The adaptive grid is quadtree grid, and the mesh refinement is carried out in the region with rapid water level gradient. The numerical example shows that the adaptive grid technology can automatically identify the region with large parameter gradient and the boundary between dry and wet, adjust the grid size, accurately simulate the dynamic change of water flow propagation, and can be applied to the actual flood simulation.
[1] 王志力, 耿艳芬, 金生. 具有复杂计算域和地形的二维浅水流动数值模拟[J]. 水利学报, 2005(4):439-444. [2] 马利平, 侯精明, 张大伟. 耦合溃口演变的二维洪水演进数值模型研究[J]. 水利学报, 2019, 50(10):1253-1267. [3] 宋利祥, 徐宗学. 城市暴雨内涝水文水动力耦合模型研究进展[J]. 北京师范大学学报(自然科学版), 2019, 55(5):581-587. [4] Kim B, Sanders B F, Schubert J E, et al.Mesh type tradeoffs in 2D hydrodynamic modeling of flooding with a Godunov-based flow solver[J].Advances in Water Resources, 2014, 68(3):42-61. [5] Borthwick A G L, S. Cruz León, J. Józsa.The shallow flow equations solved on adaptive quadtree grids[J].International Journal for Numerical Methods in Fluids, 2015, 37(6):691-719. [6] Liang Q, Borthwick A.Adaptive quadtree simulation of shallow flows with wet-dry fronts over complex topography[J].Computers & Fluids, 2009, 38(2):221-234. [7] Liang Q, Borthwick A, Stelling G.Simulation of dam-and dykebreak hydrodynamics on dynamically adaptive quadtree grids[J]. International Journal for Numerical Methods in Fluids, 2004, 46(2):127-162. [8] Liang Q,Zang J, Borthwick A G L, et al.Shallow flow simulation on dynamically adaptive cut cell quadtree grids[J].International Journal for Numerical Methods in Fluids, 2010, 53(12):1777-1799. [9] 王亮, 毛科峰, 陈希, 等.Gerris数值方案及其在海洋数值模拟中的应用[J]. 南京信息工程大学学报(自然科学版), 2015, 7(1):92-96. [10] Fraccarollo L, Toro E F.Experimental and numerical assessment of the shallow water model for two-dimensional dam-break type problems[J].Journal of Hydraulic Research, 1995, 33(6):843-864. [11] Rogers B, Fujihara M, Borthwick A G L.Adaptive Q-tree Godunovtype scheme for shallow water equations[J].International Journal for Numerical Methods in Fluids, 2001, 35(3):247-280.