Flood Simulation Using a Well-Balanced Shallow Flow Model

This work extends and improves a one-dimensional shallow flow model to two-dimensional (2D) for real-world flood simulations. The model solves a prebalanced formulation of the fully 2D shallow water equations, including friction source terms using a finite volume Godunov-type numerical scheme. A reconstruction method ensuring nonnegative depth is used along with a Harten, Lax, and van Leer approximate Riemann solver with the contact wave restored for calculation of interface fluxes. A local bed modification method is proposed to maintain the well-balanced property of the algorithm for simulations involving wetting and drying. Second-order accurate scheme is achieved by using the slope limited linear reconstruction together with a Runge-Kutta time integration method. The model is applicable to calculate different types of flood wave ranging from slow-varying inundations to extreme and violent floods, propagating over complex domains including natural terrains and dense urban areas. After validating against an analytical case of flow sloshing in a domain with a parabolic bed profile, the model is applied to simulate an inundation event in a 36?km2 floodplain in Thamesmead near London. The numerical predictions are compared with analytical solutions and alternative numerical results.     

Divergence Form for Bed Slope Source Term in Shallow Water Equations

A novel technique is presented for the treatment of the bed slope source terms within the numerical solution of the shallow water equations. The proposed method consists of writing the bed slope source term as the divergence of a proper matrix, related to the static force due to bottom slope. Such a technique is founded on analytical reasoning and is physically based, so that it can be easily applied to a wide range of numerical models, as it is completely independent of any adopted discretization technique, and requires a minimum computational effort. Herein, we show an application to a Godunov-type model, second order accurate both in space and time, based on the finite-volume method. The presented technique leads to a strong improvement in the source terms numerical treatment, especially for steady states, in which flux gradients are exactly balanced by source terms. A surprising degree of simplicity is maintained, with respect to different existing methods. The numerical model has been applied to a set of classical test cases and to a selected laboratory experiment, in order to verify its stability, accuracy, and applicability to practical real-world cases.     

Well-Balanced Bottom Discontinuities Treatment for High-Order Shallow Water Equations WENO Scheme

A finite volume well-balanced weighted essentially nonoscillatory (WENO) scheme, fourth-order accurate in space and time, for the numerical integration of shallow water equations with the bottom slope source term, is presented. The main novelty introduced in this work is a new method for managing bed discontinuities. This method is based on a suitable reconstruction of the conservative variables at the cell interfaces, coupled with a correction of the numerical flux based on the local conservation of total energy. Further changes regard the treatment of the source term, based on a high-order extension of the divergence form for bed slope source term method, and the application of an analytical inversion of the specific energy-depth relationship. Two ad hoc test cases, consisting of a steady flow over a step and a surge crossing a step, show the effectiveness of the method of treating bottom discontinuities. Several standard one-dimensional test cases are also used to verify the high-order accuracy, the C-property, and the good resolution properties of the resulting scheme, in the cases of both continuous and discontinuous bottoms. Finally, a comparison between the fourth-order scheme proposed here and a well-established second-order scheme emphasizes the improvement achieved using the higher-order approach.     

Well-Balanced Scheme between Flux and Source Terms for Computation of Shallow-Water Equations over Irregular Bathymetry

Although many numerical techniques such as approximate Riemann solvers can be used to analyze subcritical and supercritical flows modeled by hyperbolic-type shallow-water equations, there are some difficulties in practical applications due to the numerical unbalance between source and flux terms. In this study, a revised surface gradient method is proposed that balances source and flux terms. The new numerical model employs the MUSCL–Hancock scheme and the HLLC approximate Riemann solver. Several verifications are conducted, including analyses of transcritical steady-state flows, unsteady dam break flows on a wet and dry bed, and flows over an irregular bathymetry. The model consistently returns accurate and reasonable results comparable to those obtained through analytical methods and laboratory experiments. The revised surface gradient method may be a simple but robust numerical scheme appropriate for solving hyperbolic-type shallow-water equations over an irregular bathymetry.     

Finite Volume Model for Shallow Water Equations with Improved Treatment of Source Terms

A simple yet precise relation between the flux gradient and the bed slope source term is presented, which produces a net force within the cell with an inclined water surface, but ensures still water condition when there is no flow across the boundaries. The proposed method consists of calculating the pressure term based on the water depths at the cell vertices, which may be computed by a higher order scheme and the bed slope source term by a centered discretization technique. The methodology is demonstrated with a Godunov-type upwind finite volume formulation. The inviscid fluxes are calculated using Roe’s approximate Riemann solver and a second-order spatial accuracy is obtained by implementing multidimensional gradient reconstruction and slope limiting techniques. The accuracy and applicability of the numerical model is verified with a couple of test problems and a real flow example of tidal water movement in a stretch of River Hooghly in India.     

Modeling Landslide Dambreak Flood Magnitudes: Case Study

Landslide dams typically comprise unconsolidated and poorly sorted material and are vulnerable to rapid failure and breaching, resulting in significant and sudden flood risk downstream. Hence they constitute a serious natural hazard, and rapid assessment of the likely peak flow rate is required to enable preparation of adequate mitigation strategies. To determine the relative utility and accuracy of dambreak flood forecasts, field estimates of peak outflow rates from the failure of the Poerua landslide dam in October 1999 were compared with estimates from physical laboratory modeling, empirical methods, and computer modeling. There was reasonable agreement among the field estimates, laboratory modeling, and computer modeling. Some empirical estimates were less reliable. Reasonably reliable estimates of peak outflow can be obtained from computer model routines sufficiently rapidly to be of use in an emergency management situation. The laboratory modeling demonstrated the effect of dam batter slopes and valley bed slope on peak outflow; this information could be used to refine empirical or numerical estimates of peak outflow.     

Finite-Volume Model for Shallow-Water Flooding of Arbitrary Topography

A model based on the finite-volume method is developed for unsteady, two-dimensional, shallow-water flow over arbitrary topography with moving lateral boundaries caused by flooding or recession. The model uses Roe’s approximate Riemann solver to compute fluxes, while the monotone upstream scheme for conservation laws and predictor-corrector time stepping are used to provide a second-order accurate solution that is free from spurious oscillations. A robust, novel procedure is presented to efficiently and accurately simulate the movement of a wet/dry boundary without diffusing it. In addition, a new technique is introduced to prevent numerical truncation errors due to the pressure and bed slope terms from artificially accelerating quiescent water over an arbitrary bed. Model predictions compare favorably with analytical solutions, experimental data, and other numerical solutions for one- and two-dimensional problems.     

Discontinuous Galerkin Finite-Element Method for Transcritical Two-Dimensional Shallow Water Flows

A total variation diminishing Runge Kutta discontinuous Galerkin finite-element method for two-dimensional depth-averaged shallow water equations has been developed. The scheme is well suited to handle complicated geometries and requires a simple treatment of boundary conditions and source terms to obtain high-order accuracy. The explicit time integration, together with the use of orthogonal shape functions, makes the method for the investigated flows computationally as efficient as comparable finite-volume schemes. For smooth parts of the solution, the scheme is second order for linear elements and third order for quadratic shape functions both in time and space. Shocks are usually captured within only two elements. Several steady transcritical and transient flows are investigated to confirm the accuracy and convergence of the scheme. The results show excellent agreement with analytical solutions. For investigating a flume experiment of supercritical open-channel flow, the method allows very good decoupling of the numerical and mathematical model, resulting in a nearly grid-independent solution. The simulation of an actual dam break shows the applicability of the scheme to nontrivial bathymetry and wave propagation on a dry bed.     

Finite Volume Model for Two-Dimensional Shallow Water Flows on Unstructured Grids

A numerical model based upon a second-order upwind finite volume method on unstructured triangular grids is developed for solving shallow water equations. The HLL approximate Riemann solver is used for the computation of inviscid flux functions, which makes it possible to handle discontinuous solutions. A multidimensional slope-limiting technique is employed to achieve second-order spatial accuracy and to prevent spurious oscillations. To alleviate the problems associated with numerical instabilities due to small water depths near a wet/dry boundary, the friction source terms are treated in a fully implicit way. A third-order total variation diminishing Runge–Kutta method is used for the time integration of semidiscrete equations. The developed numerical model has been applied to several test cases as well as to real flows. Numerical tests prove the robustness and accuracy of the model.     

Analysis of Debris Wave Development with One-Dimensional Shallow-Water Equations

The objective of this contribution is to analyze the formation of debris waves in natural channels. Numerical simulations are carried out with a 1D code, based on shallow-water equations and on the weighted averaged flux method. The numerical code represents the incised channel geometry with a power-law relation between local width and flow depth and accounts for all source terms in the momentum equation. The debris mixture is treated as a homogeneous fluid over a fixed bottom, whose rheological behavior alternatively follows Herschel-Bulkley, Bingham, or generalized viscoplastic models. The code is first validated by applying it to dam-break tests on mudflows down a laboratory chute and verifying its efficiency in the simulation of rapid transients. Then, following the analytical method developed by Trowbridge, the stability of a uniform flow for a generalized viscoplastic fluid is examined, showing that debris flows become unstable for Froude numbers well below 1. Applications of the code to real debris flow events in the Cortina d’Ampezzo area (Dolomites) are presented and compared with available measured hydrographs. A statistical analysis of debris waves shows that a good representation of wave statistics can be obtained with a proper calibration of rheological parameters. Finally, it is shown that a minimum duration of debris event and channel length are required for waves showing up, and an explanation, confirmed both by field data and numerical simulations, is provided.     

Role of Artificial Dissipation Scaling and Multigrid Acceleration in Numerical Solutions of the Depth-Averaged Free-Surface Flow Equations

A numerical model is developed for solving the depth-averaged, open-channel flow equations in generalized curvilinear coordinates. The equations are discretized in space in strong conservation form using a space-centered, second-order accurate finite-volume method. A nonlinear blend of first- and third-order accurate artificial dissipation terms is introduced into the discrete equations to accurately model all flow regimes. Scalar- and matrix-valued scaling of the artificial dissipation terms are considered and their effect on the accuracy of the solutions is evaluated. The discrete equations are integrated in time using a four-stage explicit Runge–Kutta method. For the steady-state computations, local time stepping, implicit residual smoothing, and multigrid acceleration are used to enhance the efficiency of the scheme. The numerical model is validated by applying it to calculate steady and unsteady open-channel flows. Extensive grid sensitivity studies are carried out and the potential of multigrid acceleration for steady depth-averaged computations is demonstrated.     

Numerical Modeling of Basin Irrigation with an Upwind Scheme

In recent years, upwind techniques have been successfully applied in hydrology to simulate two-dimensional free surface flows. Basin irrigation is a surface irrigation system characterized by its potential to use water very efficiently. In basin irrigation, the field is leveled to zero slope and flooded from a point source. The quality of land leveling has been shown to influence irrigation performance drastically. Recently, two-dimensional numerical models have been developed as tools to design and manage basin irrigation systems. In this work, a finite volume-based upwind scheme is used to build a simulation model considering differences in bottom level. The discretization is made on triangular or quadrilateral unstructured grids and the source terms of the equations are given a special treatment. The model is applied to the simulation of two field experiments. Simulation results resulted in a clear improvement over previous simulation efforts and in a close agreement with experimental data. The proposed model has proved its ability to simulate overland flow in the presence of undulated bottom elevations, inflow hydrographs, and colliding fronts.     

Discretization of Integral Equations Describing Flow in Nonprismatic Channels with Uneven Beds

Application of the finite-volume method in one dimension for open channel flow predictions mandates the direct discretization of integral equations for mass conservation and momentum balance. The integral equations include source terms that account for the forces due to changes in bed elevation and channel width, and an exact expression for these source term integrals is presented for the case of a trapezoidal channel cross section whereby the bed elevation, bottom width, and inverse side slope are defined at cell faces and assumed to vary linearly and uniformly within each cell, consistent with a second-order accurate solution. The expressions may be used in the context of any second-order accurate finite-volume scheme with channel properties defined at cell faces, and it is used here in the context of the Monotone Upwind Scheme for Conservation Laws (MUSCL)-Hancock scheme which has been adopted by many researchers. Using these source term expressions, the MUSCL-Hancock scheme is shown to preserve stationarity, accurately converge to the steady state in a frictionless flow test problem, and perform well in field applications without the need for upwinding procedures previously reported in the literature. For most applications, an approximate, point-wise treatment of the bed slope and nonprismatic source terms can be used instead of the exact expression and, in contrast to reports on other finite-volume-based schemes, will not cause unphysical oscillations in the solution.     

Creep Modeling in Excavation Analysis of a High Rock Slope

Based on the distinct element method, a numerical procedure is presented for simulation of creep behavior of jointed rock slopes due to excavation unloading. The Kelvin model is used to simulate viscous deformation of joints. A numerical scheme is introduced to create incremental contact forces, which are equivalent to producing creep deformation of a rock-joint system. The corresponding displacement of discrete blocks due to creep deformation of contact joints can be calculated by equilibrium iteration. Comparisons of results between the numerical model and theoretical solutions of a benchmark example show that the presented model has excellent accuracy for analysis of creep deformation of rock-joint structures. As an application of the model, residual deformations of the high rock slopes of the Three Gorges shiplock due to excavation unloading and creep behavior are investigated. By simulating the actual excavation process, the deformation history of a shiplock slope is studied. Good agreement has been achieved between numerical prediction and field measurements. It demonstrates the effectiveness of the presented model in analysis of the creep deformation due to excavation unloading of high rock slopes.     

Treatment of Natural Geometry in Finite Volume River Flow Computations

A method is proposed for the treatment of irregular bathymetry in one-dimensional finite volume computations of open-channel flow. The strategy adopted is based on a reformulation of the Saint-Venant equations. In contrast with the usual treatment of topography effects as source terms, the method accounts for slope and nonprismaticity by modifying the momentum flux. This makes it possible to precisely balance the hydrostatic pressure contributions associated with variations in valley geometry. The characteristic method is applied to the revised equations, yielding topographic corrections to the numerical fluxes of an upwind scheme. Further adaptations endow the scheme with an ability to capture transcritical sections and wetting fronts in channels of abrupt topography. To test the approach, the scheme is first applied to idealized benchmark problems. The method is then used to route a severe flood through a complex river system: the Tanshui in Northern Taiwan. Computational results compare favorably with gauge records. Discrepancies in water stage represent no more than a fraction of the magnitude of typical bathymetry variations.     

边坡复杂地质结构三维可视化及数值模型构建

地质模型的精确性是露天矿边坡稳定性数值分析可靠性的决定因素。针对地质体间相互交叉、穿切的复杂性和无规律性,研究了基于非均匀有理B样条NURBS（Non-Uniform Rational B-Spline）的自由曲面生成技术。以矿山普遍采用的二维地质勘探剖面图为基础数据,通过对其进行矢量化处理提取出岩土体边界离散点信息,建立起NURBS三维地质体分界曲面,据此以MIDAS软件为平台制定了边坡工程复杂地质结构三维可视化及数值模型构建方案,并选取江西城门山露天矿边坡工程进行了应用实践。结果表明,运用该技术构建复杂边坡地质体模型,方法简单可靠,可极大提高模型的精度和建模效率,不仅使复杂地质体得以形象再现,同时可为边坡稳定性分析提供可靠的有限元数值计算模型。     

Simulation of a 3D Numerical Viscous Wave Tank

A numerical scheme was developed to solve the unsteady three-dimensional (3D) Navier–Stokes equations and the fully nonlinear free surface boundary conditions for simulating a 3D numerical viscous wave tank. The finite-analytic method was used to discretize the partial differential equations, and the marker-and-cell method was extended to treat the 3D free surfaces. A piston-type wave generator was incorporated in the computational domain to generate the desired incident waves. This wave tank model was applied to simulate the generation and propagation of a solitary wave in the wave tank and the diffraction of periodic waves by a semiinfinite breakwater. The computation was carried out by a PC cluster established by connecting several personal computers. The message passing interface (MPI) parallel language and MPICH software were used to write the computer code for parallel computing. High consistency between the numerical results and the theoretical solutions for the wave and velocity profiles confirms the accuracy of the proposed wave tank model.     

Upwind Conservative Scheme for the Saint Venant Equations

An upwind conservative scheme with a weighted average water-surface-gradient approach is proposed to compute one-dimensional open channel flows. The numerical scheme is based on the control volume method. The intercell flux is computed by the one-sided upwind method. The water surface gradient is evaluated by the weighted average of both upwind and downwind gradients. The scheme is tested with various examples, including dam-break problems in channels with rectangular and triangular cross-sections, hydraulic jump, partial dam-break problem, overtopping flow, a steady flow over bump with hydraulic jump, and a dam-break flood case in a natural river valley. Comparisons between numerical and exact solutions or experimental data demonstrated that the proposed scheme is capable of accurately reproducing various open channel flows, including subcritical, supercritical, and transcritical flows. The scheme is inherently robust, stable, and monotone. The scheme does not require any special treatment, such as artificial viscosity or front tracking technique, to capture steep gradients or discontinuities in the solution.     

Numerical Simulation of Flushing of Trapped Salt Water from a Bar-Blocked Estuary

Results of a numerical simulation investigating the complicated flushing process of an isolated trapped volume of salt water from a bar-blocked estuary are presented. A multiphase model, a part of the commercial code FLUENT 6.2, is applied. The governing equations together with initial and boundary conditions and the numerical scheme are described. The time-dependent salt-wedge position, vertical-density distribution, and proportion of total input kinetic energy converted into potential energy are examined for various incoming flow densimetric Froude number and estuary bed slope. The vertical position and thickness of the interfacial mixed layer between freshwater and salt water as well as the local gradient Richardson number are determined from simulated density profiles and velocity fields. The good agreement between the simulated and measured results indicates that the numerical model can be successfully applied to investigate the complex flushing process involving stratified flow.     

Eulerian-Lagrangian Numerical Scheme for Simulating Advection, Dispersion, and Transient Storage in Streams and a Comparison of Numerical Methods

Past applications of one-dimensional advection, dispersion, and transient storage zone models have almost exclusively relied on a central differencing, Eulerian numerical approximation to the nonconservative form of the fundamental equation. However, there are scenarios where this approach generates unacceptable error. A new numerical scheme for this type of modeling is presented here that is based on tracking Lagrangian control volumes across a fixed (Eulerian) grid. Numerical tests are used to provide a direct comparison of the new scheme versus nonconservative Eulerian numerical methods, in terms of both accuracy and mass conservation. Key characteristics of systems for which the Lagrangian scheme performs better than the Eulerian scheme include: nonuniform flow fields, steep gradient plume fronts, and pulse and steady point source loadings in advection-dominated systems. A new analytical derivation is presented that provides insight into the loss of mass conservation in the nonconservative Eulerian scheme. This derivation shows that loss of mass conservation in the vicinity of spatial flow changes is directly proportional to the lateral inflow rate and the change in stream concentration due to the inflow. While the nonconservative Eulerian scheme has clearly worked well for past published applications, it is important for users to be aware of the scheme’s limitations.