1923 Gleno Dam Break: Case Study and Numerical Modeling

On the morning of December 1, 1923, the Gleno Dam (located in the Central Italian Alps) suddenly collapsed a few days after the first complete reservoir filling. Nearly 4.5×106??m3 of water was released. The consequent inundation caused significant destruction along the downstream valley and a death toll of at least 356 lives. This failure is the only historical case of dam break caused by structural deficiencies that has occurred in Italy. As a result, it has deeply influenced the evolution of Italian regulations regarding dam design and hydraulic risk evaluation. However, in spite of its relevance, this event has never been characterized from a hydraulic standpoint. This paper reports the main information obtained from the analysis of a vast amount of historical documents regarding the Gleno Dam break to set up a case study useful for validating dam-break models in mountain settings. Moreover, it presents the main results of one-dimensional (1D) modeling of the dam break wave propagation accomplished with a first-order finite volume numerical scheme recently proposed in the literature for field applications. The overall effectiveness and reliability of the model are evaluated for this case characterized by very irregular topography. Finally, the practical relevance of several choices that the numerical reconstruction of this kind of event demands is tested.     

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.     

Two-Dimensional Fast-Response Flood Modeling: Desktop Parallel Computing and Domain Tracking

Emergency flood management is enhanced by using models that can estimate the timing and location of flooding. Typically, flood routing and inundation prediction is accomplished by using one-dimensional (1D) models. These have been the models of choice because they are computationally simple and quick. However, these models do not adequately represent the complex physical processes present for shallow flows located in the floodplain or in urban areas. Two-dimensional (2D) models developed on the basis of the full hydrodynamic equations can be used to represent the complex flow phenomena that exist in the floodplain and are, therefore, recommended by the National Research Council for increased use in flood analysis studies. The major limitation of these models is the increased computational cost. Two-dimensional flood models are prime candidates for parallel computing, but traditional methods/equipment (e.g., message passing paradigm) are more complex in terms of code refactoring and hardware setup. In addition, these hardware systems may not be available or accessible to modelers conducting flood analyses. This paper presents a 2D flood model that implements multithreading for use on now-prevalent multicore computers. This desktop parallel computing architecture has been shown to decrease computation time by 14 times on a 16-processor computer and, when coupled with a wet cell tracking algorithm, has been shown to decrease computation by as much as 310 times. These accomplishments make high-fidelity flood modeling more feasible for flood inundation studies using readily available desktop computers.     

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.     

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.     

Case Study: Flood Mitigation of the Muda River, Malaysia

The 2003 flood of the Muda River reached 1,340?m3/s at Ladang Victoria and adversely impacted 45,000 people in Malaysia. A flood control remediation plan proposed a levee height based on a 50-year discharge of 1,815?m3/s obtained from hydrologic models. This design discharge falls outside the 95% confidence intervals of the flood frequency analysis based on field measurements. Instream sand and gravel mining operations also caused excessive riverbed degradation, which largely off sets apparent benefits for flood control. Pumping stations have been systematically required at irrigation canal intakes. Several bridge piers have also been severely undermined and emergency abutment protection works were needed in several places. Instream sand and gravel mining activities should be replaced with offstream mining in the future.     

Modeling Evapotranspiration of Two Land Covers Using Integrated Hydrologic Model

Modeling evapotranspiration (ET) distribution in shallow water table environments is of great importance for understanding and reproducing other hydrologic fluxes such as runoff and recharge. Unfortunately, ET distribution can be the most difficult hydrologic process to analyze. The partitioning of ET into upper zone ET, lower zone ET, and groundwater ET is complex because it depends on land cover and subsurface characteristics. One comprehensive distributed parameter model, integrated hydrologic model (IHM), builds on an improved understanding and characterization of ET partitioning between surface storages, vadose zone storage, and saturated groundwater storage. It provides a smooth transition to satisfy ET demand between the vadose zone and the deeper saturated groundwater. In this paper, the IHM was used to analyze ET contribution from different regions of the vadose zone and saturated zone. Rigorous testing was done on two distinct land covers, grass land and forest land, at a study site in West-Central Florida. Sensitivity analysis on the key parameters was investigated and influence of parameters on ET behavior was also discussed. Statistics with the root mean square error and mean bias error for forest total ET were about 1.46 and 0.04 mm/day, respectively, and 1.61 and 1.07 mm/day for grass total ET. Modeling results further proved that ET distributions from the upper and lower soil and water table, while incorporating field-scale variability of soil and land cover properties, can be predicted reasonably well using IHM model.     

Lagrangian Modeling of Weakly Nonlinear Nonhydrostatic Shallow Water Waves in Open Channels

A Lagrangian, nonhydrostatic, Boussinesq model for weakly nonlinear and weakly dispersive flow is presented. The model is an extension of the hydrostatic model—dynamic river model. The model uses a second-order, staggered grid, predictor-corrector scheme with a fractional step method for the computation of the nonhydrostatic pressure. Numerical results for solitary waves and undular bores are compared with Korteweg-de Vries analytical solutions and published numerical, laboratory, and theoretical results. The model reproduced well known features of solitary waves, such as wave speed, wave height, balance between nonlinear steepening and wave dispersion, nonlinear interactions, and phase shifting when waves interact. It is shown that the Lagrangian moving grid is dynamically adaptive in that it ensures a compression of the grid size under the wave to provide higher resolution in this region. Also the model successfully reproduced a train of undular waves (short waves) from a long wave such that the predicted amplitude of the leading wave in the train agreed well with published numerical and experimental results. For prismatic channels, the method has no numerical diffusion and it is demonstrated that a simple second-order scheme suffices to provide an efficient and economical solution for predicting nonhydrostatic shallow water flows.     

Case Study: Design of Flood Control Systems on the Vara River by Numerical and Physical Modeling

In the present paper, we investigate the effectiveness of a flood defense project based on storage reservoirs, presently under study for the Magra River and Vara River (Italy). We have focused the analysis on two detention reservoirs and studied their response to different hydrological scenarios mostly in terms of flood mitigation efficiency, leaving aside sediment transport issues. The analysis has been carried out with the aid of a physical model and one-dimensional numerical simulations. Experimental and numerical simulations have been performed spanning a wide range of hydrological conditions. Some of the results can be generalized for different applications where similar flood control systems are employed.     

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.     

Flood Inundation Modeling with an Adaptive Quadtree Grid Shallow Water Equation Solver

Flood risk studies require hydraulic modeling in order to estimate flow depths and other hydraulic variables in the floodplain for a wide range of input conditions. Currently there is a need to improve the computational efficiency of fully two-dimensional numerical models for large-scale flood simulation. This paper describes an adaptive quadtree grid-based shallow water equation solver and demonstrates its capability for flood inundation modeling. Due to the grid dynamically adapting to dominant flow features such as steep water surface gradients and wet-dry fronts, the approach is both efficient and accurate. The quadtree model is applied to a realistic scenario of flood inundation over an urban area of 36?km2, resulting from the flood defenses breaching at Thamesmead on the River Thames, United Kingdom. The results of the simulation are in close agreement with alternative predictions obtained using the commercially available software TUFLOW.     

Hydraulic Jet Control for River Junction Design of Yuen Long Bypass Floodway, Hong Kong

The Yuen Long Bypass Floodway (YLBF) was designed to collect flows from the Sham Chung River (SCR) and the San Hui Nullah (SHN) and to serve as a diversion channel of the Yuen Long Main Nullah (YLMN). Under a 200-year return period design condition, the floodway was designed (1) to divert a flow of approximately 38?m3/s from the supercritical YLMN flow and (2) to convey a total combined flow of 278?m3/s to downstream within acceptable flood levels. The success of the design depends critically on complicated junction flow interactions that cannot be resolved by 1D unsteady flow models. These features include the supercritical-subcritical flow transition at the San Hui-Floodway (SHN-YLBF) junction and the diversion of part of the supercritical flow from the Main Nullah (YLMN). A laboratory Froude scale physical model was constructed to study water stages and flow characteristics in the floodway and to investigate optimal design arrangements at channel junctions and transitions. This paper summarizes the main features of the unique river junction network, in particular the use of the hydraulic jet principle at the SHN-YLBF junction to lower flood levels. In addition, a numerical flow model is employed to study flow details at the river junctions. The model is based on the general 2D shallow water equations in strong conservation form. The equations are discretized using the total variation diminishing finite-volume method which captures the discontinuity in hydraulic jumps. The numerical model predictions are well supported by the laboratory data, and the theoretical and experimental results offer useful insights for the design of urban flood control schemes under tight space constraints.     

Dynamic Modeling of Storm Surge and Inland Flooding in a Texas Coastal Floodplain

The purpose of this study is to analyze the combined effects of storm surge and inland rainfall on the floodplain of a coastal bayou in the Houston area by using dynamic hydraulic modeling. Most existing floodplains in the Houston area are defined using only rainfall as an input into steady-state hydraulic models and do not consider the impact of hurricane-induced storm surge on the floodplain. HEC-RAS, a one-dimensional flow model, was run for both steady- and unsteady-states to analyze the additional effect storm surge has on the coastal floodplain. Storm surge and rainfall data from Hurricane Ike were utilized to run an unsteady hydraulic model on Horsepen Bayou near Galveston Bay. The dynamic model generated a good match between the modeled hydrograph and measured data in the watershed. Additionally, a timing sensitivity analysis was completed by shifting the timing of the storm surge both earlier and later in time. The dynamic model revealed that the timing of both rainfall and storm surge play a significant role in the magnitude of inland flooding.     

3D Unsteady RANS Modeling of Complex Hydraulic Engineering Flows. I: Numerical Model

A general-purpose numerical method is developed for solving the full three-dimensional (3D), incompressible, unsteady Reynolds-averaged Navier-Stokes (URANS) equations in natural river reaches containing complex hydraulic structures at full-scale Reynolds numbers. The method adopts body-fitted, chimera overset grids in conjunction with a grid-embedding strategy to accurately and efficiently discretize arbitrarily complex, multiconnected flow domains. The URANS and turbulence closure equations are discretized using a second-order accurate finite-volume approach. The discrete equations are integrated in time via a dual-time-stepping, artificial compressibility method in conjunction with an efficient coupled, block-implicit, approximate factorization iterative solver. The computer code is parallelized to take full advantage of multiprocessor computer systems so that unsteady solutions on grids with 106 nodes can be obtained within reasonable computational time. The power of the method is demonstrated by applying it to simulate turbulent flow at R ? 107 in a stretch of the Chattahoochee River containing a portion of the actual bridge foundation located near Cornelia, Georgia. It is shown that the method can capture the onset of coherent vortex shedding in the vicinity of the foundation while accounting for the large-scale topographical features of the surrounding river reach.     

National Assessment of Pesticide Runoff Loads from Grass Surfaces

Pesticide runoff loads from grass surfaces were estimated through simulation experiments for 37 chemicals registered for use on U.S. lawns and golf courses. Simulation runs were made for each chemical and surface (lawns, greens, fairways) using 100-year weather records generated for nine U.S. cities. Results were summarized as mean annual and 1-in-10?year annual maximum daily pesticide loads. These loads varied greatly with pesticide, grass surface, and city, ranging from less than one to over 400??g/ha for mean annual loads and from less than one to over 500??g/ha for 1-in-10?year maximum daily loads. Mean annual loads averaged over the 37 chemicals and three grass surfaces were found to be closely related to growing season precipitation. Variations among the nine cities were well-captured by three general climate categories: humid, represented by Atlanta and Houston; mesic, as with Albany, Columbus, Madison, and Olympia; and dry, represented by Bismarck, Fresno, and Roswell. Mean annual pesticide runoff was 19, 6, and 2??g/ha in the humid, mesic, and dry regions, respectively.     

3D Unsteady RANS Modeling of Complex Hydraulic Engineering Flows. II: Model Validation and Flow Physics

A chimera overset grid flow solver is developed for solving the unsteady Reynolds-averaged Navier-Stokes (RANS) equations in arbitrarily complex, multiconnected domains. The details of the numerical method were presented in Part I of this paper. In this work, the method is validated and applied to investigate the physics of flow past a real-life bridge foundation mounted on a fixed flat bed. It is shown that the numerical model can reproduce large-scale unsteady vortices that contain a significant portion of the total turbulence kinetic energy. These coherent motions cannot be captured in previous steady three-dimensional (3D) models. To validate the importance of the unsteady motions, experiments are conducted in the Georgia Institute of Technology scour flume facility. The measured mean velocity and turbulence kinetic energy profiles are compared with the numerical simulation results and are shown to be in good agreement with the numerical simulations. A series of numerical tests is carried out to examine the sensitivity of the solutions to grid refinement and investigate the effect of inflow and far-field boundary conditions. As further validation of the numerical results, the sensitivity of the turbulence kinetic energy profiles on either side of the complex pier bent to a slight asymmetry of the approach flow observed in the experiments is reproduced by the numerical model. In addition, the computed flat-bed flow characteristics are analyzed in comparison with the scour patterns observed in the laboratory to identify key flow features responsible for the initiation of scour. Regions of maximum shear velocity are shown to correspond to maximum scour depths in the shear zone to either side of the upstream pier, but numerical values of vertical velocity are found to be very important in explaining scour and deposition patterns immediately upstream and downstream of the pier bent.     

Coupled Mechanical and Hydraulic Modeling of Geosynthetic-Reinforced Column-Supported Embankments

Geosynthetic-reinforced column-supported (GRCS) embankments have increasingly been used in the recent years for accelerated construction. Numerical analyses have been conducted to improve understanding and knowledge of this complicated embankment system. However, most studies so far have been focused on its short-term or long-term behavior by assuming an undrained or drained condition, which does not consider water flow in saturated soft soil (i.e., consolidation). As a result, very limited attention has been paid to a settlement-time relationship especially postconstruction settlement, which is critical to performance of pavements on embankments or connection between approach embankments and bridge abutments. To investigate the time-dependent behavior, coupled two-dimensional mechanical and hydraulic numerical modeling was conducted in this study to analyze a well-instrumented geotextile-reinforced deep mixed column-supported embankment in Hertsby, Finland. In the mechanical modeling, soils and DM columns were modeled as elastic-plastic materials and a geotextile layer was modeled using cable elements. In the hydraulic modeling, water flow was modeled to simulate generation and dissipation of excess pore water pressures during and after the construction of the embankment. The numerical results with or without modeling water flow were compared with the field data. In addition, parametric studies were conducted to further examine the effects of geosynthetic stiffness, column modulus, and average staged construction rate on the postconstruction settlement and the tension in the geosynthetic reinforcement.     

Hydraulic Modeling of an Automatic Upstream Water-Level Control Gate

This paper proposes an efficient mathematical model of an automatic upstream water-level control gate, called a Begemann or flap gate. This automatic gate controls the upstream level close to a reference level for free gate flow, using a counterweight to compensate for the hydraulic pressure on the gate. The proposed gate model is designed to be included in a hydraulic simulation model. A discharge law for the gate is first derived using simple physical assumptions. Then a method to compute the static equilibrium is obtained by modeling the opening force exerted by the water on the gate. This mathematical model is validated on experimental data from a small-scale gate and on other data from the literature in order to show the ability of the model to simulate various gates.