No abstract

No abstract

In the present work, a porosity-based numerical scheme for the Shallow Water Equations is presented. With the aim of accounting for the presence of storage areas, such as gardens, yards and dead zones, and for preferential flow pathways, both an isotropic storage porosity parameter and anisotropic friction are adopted. Particularly, the anisotropic effects due to the building alignments are evaluated defining conveyance porosities along principal directions and using them to express the friction losses in tensor form. The storage and conveyance porosities are evaluated from the geometry of the urban layout at a district scale and then assigned to computational cells rather than to cell sides, thus avoiding oversensitivity to the mesh design. The proposed formulation guarantees the C-property also in presence of wet-dry fronts. Model testing is performed analyzing schematic and idealized urban layouts, and against experimental data as well. The results obtained by the proposed anisotropic scheme are similar to a high-resolution model with resolved buildings, also in the presence of low-friction regimes, meanwhile with a remarkable reduction of the computational times.

Abstract. With the aim of improving resilience to flooding and increasing preparedness to face levee-breach-induced inundations, this paper presents a methodology for creating a wide database of numerically simulated flooding scenarios due to embankment failures, applicable to any lowland area protected by river levees. The analysis of the detailed spatial and temporal flood data obtained from these hypothetical scenarios is expected to contribute both to the development of civil protection planning and to immediate actions during a possible future flood event (comparable to one of the available simulations in the database) for which real-time modelling may not be feasible. The most relevant criteria concerning the choice of mathematical model, grid resolution, hydrological conditions, breach parameters and locations are discussed in detail. The proposed methodology, named RESILIENCE, is applied to a 1100 km2 pilot area in northern Italy. The creation of a wide database for the study area is made possible thanks to the adoption of a GPU-accelerated shallow-water numerical model which guarantees remarkable computational efficiency (ratios of physical to computational time up to 80) even for high-resolution meshes (2.5–5 m) and very large domains (>1000 km2).

The capability of a GPU-parallelized numerical scheme to perform accurate and fast simulations of surface runoff in watersheds, exploiting high-resolution digital elevation models (DEMs), was investigated. The numerical computations were carried out by using an explicit finite volume numerical scheme and adopting a recent type of grid called Block-Uniform Quadtree (BUQ), capable of exploiting the computational power of GPUs with negligible overhead. Moreover, stability and zero mass error were ensured, even in the presence of very shallow water depth, by introducing a proper reconstruction of conserved variables at cell interfaces, a specific formulation of the slope source term and an explicit discretization of the friction source term. The 2D shallow water model was tested against two different literature tests and a real event that recently occurred in Italy for which field data is available. The influence of the spatial resolution adopted in different portions of the domain was also investigated for the last test. The achieved low ratio of simulation to physical times, in some cases less than 1:20, opens new perspectives for flood management strategies. Based on the result of such models, emergency plans can be designed in order to achieve a significant reduction in the economic losses generated by flood events.