A coupled 1D–2D hydrodynamic model for flood simulation in flood detention basin

Qiang Liu,Yi Qin,Yang Zhang,Ziwen Li

DOI: https://doi.org/10.1007/s11069-014-1373-3

IF: 3.158

2014-09-05

Natural Hazards

Abstract:AbstractA coupled 1D–2D hydrodynamic model linking the channel and flood detention basin for flood simulation with complex topography and irregular boundary was presented in this paper. The 1D Saint-Venant equations were used for governing flows in channel, and the four-point implicit Preissmann scheme was used for discretization. To simulate floods in flood detention basin, the 2D shallow equations were adopted as governing equations, and a well-balanced, unstructured finite-volume method was proposed for the numerical resolution. The 2D model is based on a new formulation of the classical shallow water equations in hyperbolic conservation form. The numerical fluxes are computed by HLLC algorithm, and the MUSCL–Hancock predictor–corrector scheme is used to achieve high-accuracy and high-resolution results. A simple and efficient method was proposed to reflect the coupled hydraulic connection between the channel and the flood detention basin. The novelties of the proposed model include (1) a robust method for wetting/drying treatment and (2) an efficient way to link the 1D and 2D models. The 1D, 2D, and coupled models are all tested through some benchmark cases, and numerical results validated the accuracy of the models. Furthermore, the coupled 1D–2D model was used for a real flood simulation in Jiakouwa flood detention basin, China. The flood-risk information including flood arrival time and maximal water depth was mapped using GIS. Those flood-risk maps can be used as an important decision-making basis of flood control and rescue for the flood control departments at all levels.

geosciences, multidisciplinary,water resources,meteorology & atmospheric sciences

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to conduct flood simulation in flood detention areas with complex terrains and irregular boundary conditions. Specifically, the author proposes a coupled one - dimensional (1D) and two - dimensional (2D) hydrodynamic model to simulate the flow of flood between rivers and flood detention areas. This model aims to predict the behavior of flood in these areas more reasonably and accurately, especially when floods occur simultaneously in rivers and flood detention areas, and the hydraulic connection between them needs to be considered.

### Main features of the model:
1. **1D model**: The Saint - Venant equations are used to describe the water flow in the river, and the four - point implicit Preissmann scheme is used for discretization.
2. **2D model**: The shallow water equations are used as the control equations, and the unstructured finite - volume method is used for numerical solution. The numerical flux is calculated by the HLLC algorithm, and the MUSCL - Hancock prediction - correction scheme is used to achieve high - precision and high - resolution results.
3. **Coupling method**: A simple and effective method is proposed to reflect the hydraulic connection between the river and the flood detention area. This method is achieved by connecting the 1D river model and the normal flow boundary conditions of the 2D domain through source terms.

### Innovation points:
1. **Wet/dry treatment**: A robust method is proposed to deal with the wet/dry interface.
2. **1D - 2D model coupling**: An efficient method is provided to link the 1D and 2D models.

### Verification and application:
- **1D model verification**: The accuracy of the 1D model is verified through a 50 - kilometer - long channel case.
- **2D model verification**: The accuracy of the 2D model is verified through a two - dimensional circular dam - break flood simulation case.
- **1D - 2D coupled model verification**: The effectiveness of the coupled model is verified through a simple test benchmark.
- **Practical application**: The coupled 1D - 2D model is applied to the actual flood simulation in the Jiakouwa flood detention area of Daqing River in China, and risk maps of flood arrival time and maximum water depth are generated. These maps can be used as an important basis for flood control departments at all levels to make flood control and rescue decisions.

### Formula display:
- **1D model control equations**:
\[
\frac{\partial z}{\partial t}+\frac{1}{B}\frac{\partial Q}{\partial x} = q
\]
\[
\frac{\partial Q}{\partial t}+\frac{gA}{C_0}\frac{\partial z}{\partial x}+\frac{2Q}{A}\frac{\partial Q}{\partial x}=\frac{Q^2}{A^2}\frac{\partial A}{\partial x}-gQ\left|Q\right|\frac{1}{Ac^2R}
\]
where \(z\) is the water level, \(B\) is the river width, \(Q\) is the flow rate, \(q\) is the lateral flow rate, \(A\) is the wetted perimeter area, \(g\) is the gravitational acceleration, \(R\) is the hydraulic radius, and \(c\) is the Chezy coefficient.

- **2D model control equations**:
\[
\frac{\partial U}{\partial t}+\frac{\partial E}{\partial x}+\frac{\partial G}{\partial y}=S
\]
where:
\[
U=\begin{pmatrix}h\\hu\\hv\end{pmatrix},\quad E = \begin{pmatrix}hu\\hu^{2}+g(h^{2}-b^{2})/2\\huv\end{pmatrix},\quad G=\begin{pmatrix}hv\\huv\\hv^{2}+g(h^{2}-b^{2})/