^{1}

^{2}

This present study develops a 2-D numerical scheme to simulate the velocity and depth on the actual terrain by using shallow water equations. The computational approach uses the HLL scheme as a basic building block, treats the bottom slope by lateralizing the momentum flux, then refines the scheme using the Strang splitting to deal with the frictional source term. Besides, a decoupled algorithm is also adopted to compute the aggradation and degradation of bed-level elevation by using the Manning-Strickler formula and Exner’s relationship. The main purpose is to set up the window interface of 2-D numerical model and increase the realization of engineers on these characteristics of hydraulic treatment and maintenance.

Many hazards caused by typhoons in Taiwan have resulted in large damages or pecuniary loss in recent years. The function of local drainage in urban or rural regions usually plays an important role in the task of controlling floods. The detention facility is one of the usual methods to control floods in stream planning. The present research presents a two-dimensional numerical model, DBFVM-2D, which is developed by using the shallow water equations and finite volume method to assess the function of detention facility. The shallow water equations are derived from the Navier-Stokes Equations by depth-averaged integration in vertical direction.

In consideration of frequent flood disasters in Taiwan, hydraulic engineering or drainage system planning need more efficient or accurate hydraulic analyses. Besides, a numerical model is used for inputting disaster conditions for computing the possible flood range before the flood disaster prevention and response operation plan so as to predict possible financial losses or personnel injuries and precede various solutions and strategies with images and other disaster prevention and response data analyses.

Meanwhile, the protected disaster scale should be preset before planning the flood disaster prevention and response operation, and the numerical model inputs related disaster conditions for computing the possible flood range and estimating possible flood range and depth. Based on above disaster scale settings, possible financial losses or personnel injuries therefore could be predicted, and various solutions and strategies are preceded with images and other disaster prevention and response data analyses. It is regarded as the objective of flood disaster potential development.

In the planning of soil and water conservation for slope land development, the runoff increased by watershed sediment getting in the river course or development behaviors could be controlled the sediment or discharge through various soil and water conservation constructions, such as groundsill works, detention basins, and settling basins. Engineers used to precede plan and design according to empirical formula; a numerical model for the situation simulation of allocated equipment could help understand the hydraulic characteristics in order to review the plan and design being played.

Dam-break caused by floods or earthquakes or sudden opening of gates could result in surges on the river flow which are delivered towards downstream. The use of a high-resolution numerical method in 2D shallow water equations to solve surges has been studied in four categories of numerical methods. For instance, MacCormack numerical method and the explicit scheme of Finite Differential Method (FDM) were employed to simulate 2D wave propagations by Garcia & Kahawital and Fennema & Chaudhry; Akanbi & Katopodes applied Finite Element Method (FEM) to simulate flood waves; and, Katopodes & Stelkoff simulated 2D dam-break flows with Method of Characteristics (MOC) [

This study aims to develop a 2D numerical model, in which shallow water equations are regarded as the governing equation and Finite Volume Method is applied to develop user-friendly window interface. Comparing to other numerical methods, Finite Volume Method presents the following advantages:

It could apply irregular grids in Finite Element that it is applicable to various irregular landforms.

Finite Volume Method would discretize conservative differential equations.

The calculation speed of Finite Volume Method equals to Finite Differential Method, but is faster than Finite Element Method that it could reduce large computing time.

Conservative equations of flow continuity and momentum are proceeded integration discretization, when the conservative characteristics are remained, i.e. keeping the conservative properties of some physical quantity, so that it is easier to analyze the flow of discontinuity like seismic waves.

Commonly used inundation models contain 1D steady flow or unsteady flow inundation model, quasi 2D inundation model (or cell model), and 2D inundation model (including Block-Street Network Inundation Model, SOBEK Model, and FLO-2D Model). Each inundation model presents the advantages and disadvantages and generally requires expensive license fee.

Finite Volume Method therefore is utilized for establishing a 2D numerical model and developing the user- friendly window interface by regarding shallow water equations as the governing equation in the present paper. The model covers the functions of dam-break flow, flood inundation simulation, and bed sediment load transport calculation, which could be applied to the hydraulic calculation of groundsill works and detention basins. Moreover, the uncoupled mode is used for calculating the bed load, and sediment continuity equations could simulate the bed scour and deposition. They are beneficial to the application of project planning.

For homogeneous and incompressible flow, the Navier-Stokes equations in a Cartesian coordinate system could be denoted as [

Continuity equation:

Momentum equation:

x direction

y direction

z direction

where x, y , z are the coordinate axes of the Cartesian coordinate system, t is the time, u ,v, w are the components of flow velocity at x, y, z directions,

Preceding the flow depth integration of above equations, the depth-averaged equation could be acquired to derive the following shallow water equations [

Continuity equation:

Momentum equation of water depth-averaged at x direction

Momentum equation of water depth-averaged at y direction

And energy slope

where n is the Manning coefficient, C_{0} = 1 (for SI unit), and C_{0} = 1.49 (for imperial unit).

Explicit uniform grid discretization governing equations in Finite Volume Method are applied to the numerical method, and the Godunov-type scheme is used for solving Riemann problem [

Furthermore, the uncoupled mode is utilized for calculating the bed sediment of hyper concentrated flow. In the beginning (t = 0), the bed-level elevation is assumed known and the relative flow parameters could be acquired by solving shallow water equations, when the sediment transport is not taken into account. Such flow parameters are further used for calculating the bed sediment load of each grid, whose in/out sediment is used for estimating the scour and deposition of each grid. The sediment continuity equation is further utilized for calculating the bed-level elevation:

where p is the bed porosity, z_{b} the bed-level elevation, and q_{s} the volumetric sediment load of per unit width. A new bed profile (

Three bed sediment load equations of Schoklitsch, Meyer-Peter et al., and Einstein sediment load are currently used in this study [

where

Assuming a 200 m × 200 m plane for the simulation, where an asymmetric opening with the width 75 m is in the middle,

The simulation area is a 1.6 m (L) × 0.6 m (W) closed channel, where a 0.12 m × 0.12 m square column is placed in the middle at x = 0.9 m; an initial depth 0.3 m appears before x = 0.4 m and the downstream depth 0.01 m to simulate the 2 cm × 2 cm grid with the Manning n 0.01 and the simulation time 2 sec. The results are shown in

The retrogressive erosion simulation assumes the channel bed profile suddenly appears rapid changes. The simulation data are acquired from the literature [^{2}∙s^{−1},

The functions of the developed model contain dam-break flow, flood inundation simulation, and bed load sediment transport calculation, which could be applied to the hydraulic calculation of groundsill works and detention facilities. Besides, the uncoupled mode is utilized for calculating the bed load, and sediment continuity equations are applied to simulate the bed scour and deposition. Moreover, all programming codes in the numerical model are self-developed that they present high expandability. Meanwhile, Finite Volume Method is utilized for computing the model that it could effectively calculate the surges of dam-break flow. Connecting the uncoupled mode with sediment continuity equations and calculating the bed sediment load with the bed scour and deposition could assist engineers in modifying the construction design.

The paper was supported in part by the Ministry of Science and Technology, Taiwan (NSC 98-2622-B-270-001- CC3) and Opening fund of State Key Laboratory of Geohazard Prevention and Geoenvironment Protection at Chengdu University of Technology (SKLGP2013K008).

Szu-Hsien Peng,Chuan Tang, (2015) Development and Application of Two-Dimensional Numerical Model on Shallow Water Flows Using Finite Volume Method. Journal of Applied Mathematics and Physics,03,989-996. doi: 10.4236/jamp.2015.38121