A fully 3-D numerical model to predict flood wave propagation and assess efficiency of flood protection measures

The execution of hydraulic models at large spatial scales has yielded a step change in our understanding of flood risk. Yet their necessary simplification through the use of coarsened terrain data results in an artificially smooth digital elevation model with diminished representation of flood defense structures. Current approaches in dealing with this, if anything is done at all, involve either employing incomplete inventories of flood defense information or making largely unsubstantiated assumptions about defense locations and standards based on socioeconomic data. Here, we introduce a novel solution for application at scale. The geomorphometric characteristics of defense structures are sampled, and these are fed into a probabilistic algorithm to identify hydraulically relevant features in the source digital elevation model. The elevation of these features is then preserved during the grid coarsening process. The method was shown to compare favorably to surveyed U.S. levee crest heights. When incorporated into a continental-scale hydrodynamic model based on LISFLOOD-FP and compared to local flood models in Iowa (USA), median correspondence was 69% for high-frequency floods and 80% for low-frequency floods, approaching the error inherent in quantifying extreme flows. However, improvements versus a model with no defenses were muted, and risk-based deviations between the local and continental models were large. When simulating an event on the Po River (Italy), built and tested with higher quality data, the method outperformed both undefended and even engineering-grade models. As such, particularly when employed alongside model components of commensurate quality, the method here generates improved-accuracy simulations of flood inundation. Plain Language SummaryTraditional flood risk assessments are carried out using computer models built with local data, but their spatial coverage is impaired by how expensive and time-consuming they are to produce. Recent advances in data availability, understanding of necessary physical process representation, and computational capacity have enabled hydraulic models of the entire globe to be built in an automated fashion at a fraction of the financial and human cost. However, their accuracy can be significantly impaired by a lack of information on flood defenses. As the model is built, elevation data are coarsened to reduce the number of calculations required to simulate flooding over such wide areas. This results in flood defense structures being smoothed out of the terrain information used in the model. Publicly available defense inventories are of insufficient coverage to ameliorate this issue. In this paper, a method is presented, which automatically detects levee-like features in high-resolution elevation data and accurately represents their heights during this necessary coarsening process. Simulating flood inundation over this "defended" topography results in high correspondence between local models and observations for test cases in the United States and Italy, with improve...

Section: Methodsmentioning

confidence: 99%

A New Automated Method for Improved Flood Defense Representation in Large‐Scale Hydraulic Models

Wing

Bates

Neal

et al. 2019

“…The same solver was previously used to conduct unsteady RANS and DES simulations in our group for complex turbulent flows. These studies (e.g., Cheng et al, 2018;Horna-Munoz & Constantinescu, 2018Wu et al, 2020) also include detailed validation with experimental data.…”

Section: Viscous Flow Solver and Turbulence Modelmentioning

confidence: 99%

Flow and Entrainment Mechanisms Around a Freshwater Mussel Aligned With the Incoming Flow

Constantinescu

Zeng

2020

Self Cite

The paper uses information from eddy-resolving simulations to characterize flow and turbulence around freshwater mussels at the organism scale. The focus is on the simplest case of a partially burrowed, isolated mussel aligned with the incoming flow in an open channel. The wake structure and the capacity of the flow to displace the mussel from the bed substrate and to induce local scour around the mussel are investigated as a function of the ratio between the height of the exposed part of the mussel, h, and the total mussel height, d, and as a function of the filtration velocity ratio (VR) between the incoming channel flow velocity, U 0 , and the mean velocity inside the excurrent siphon, U e. As opposed to the flow past most surface-mounted obstacles where the bed shear stresses are reduced in the wake of the body, the capacity of the flow to induce bed erosion behind isolated mussels aligned with the flow is relatively large because of strong downwelling motions inside the horizontal separated shear layers. These flow features are associated with the formation of counterrotating base vortices of unequal coherence that induce upwash behind the mussel. The total circulation of the base vortices increases with increasing h/d and VR. Though both symmetrical and antisymmetrical shedding are observed in the wake, the antisymmetrical mode dominates, and its strength increases with increasing VR. The nondimensional streamwise force acting on the emerged part of the mussel's shell increases with increasing h/d and VR. Finally, the paper discusses the effects of varying h/d and VR on the dynamics and dilution of the jet of filtered water originating in the excurrent siphon, which is important to understand how mussels affect mixing and water quality (e.g., nutrient availability and phytoplankton concentration) in natural streams.

“…A related question is whether the shallow-water equations are well suited to deal with dam-break flows that experience substantial three-dimensional effects (for instance, because of the valley narrowings and widenings). Horna Munoz and Constantinescu (2018) conducted a comparison between the shallow-water equations and three-dimensional Reynolds-averaged Navier-Stokes equations and found differences of about 10% in the computed hydrographs along the Iowa River. As this value is below the uncertainty range in our computations, we conclude that the errors made in neglecting three-dimensional effects are consistent with our general computational framework.…”

Section: Relevance Of the Computational Frameworkmentioning

confidence: 99%

Hydraulic Reconstruction of the 1818 Giétro Glacial Lake Outburst Flood

Ancey

Bardou

Funk

et al. 2019

In the spring of 1818, ice avalanches from the Giétro Glacier created an ice dam, which in turn formed a glacial lake in the Drance Valley (Canton of Valais, Switzerland). Today, its maximum volume is estimated to have been 25×106 m3. Cantonal authorities commissioned an engineer named Ignaz Venetz to mitigate the risk of the ice dam's failure. He supervised the construction of a tunnel through which a large volume of water was drained as the lake rose (9×106 m3 according to his estimates and 11×106 m3 according to our model). After 2.5 days of slow drainage, the ice dam failed on 16 June 1818 and caused major flooding in the Drance Valley up to 40 km downstream, resulting in about 40 deaths. Venetz's lake monitoring notes, numerous testimonies gathered in the disaster's aftermath, and our field survey have made it possible to collect a wealth of information on this event, which is one of the world's major documented glacial lake outburst floods. Reconstructing major outburst floods remains challenging because not only do they involve enormous volumes of water spreading over long distances but they are also associated with additional physical processes such as massive erosion; intense transport of ice, sediment, and debris; and damage to vegetation and buildings. This paper attempts to reconstruct the 1818 Giétro flood by focusing on its water component. We develop a simple model to estimate the initial hydrograph during the slow drainage and failure phases. The flood's features are deduced by solving the shallow‐water equations numerically. The computational framework involves six free parameters, of which five are constrained by physical considerations. Using iterative manual parameter adjustments, we matched the numerical simulations to the historical data. We found that the peak discharge was close to 14,500 m3/s, the flood's front velocity was about 6 m/s, and flow depth varied considerably along the River Drance's bed (from 30 m just downstream of the ice dam to 2 m on the alluvial fan, 24 km west of the dam). To achieve a good agreement between computations and historical data, we had to select a high value for the Manning friction coefficient n (with n as large as 0.08 s/m1/3). As the Drance Valley is narrow, high flow resistance caused the flood's leading edge to behave like a plug, moving at a fairly constant velocity, with little dependence on what happened behind it. This result may explain why a simple flood routing model is able to reproduce the flood's features, because in an Alpine valley, a lateral spreading of the water volume is limited.