scholarly journals Towards fast large-scale flood simulations using 2D Shallow water modelling with depth-dependant porosity

Author(s):  
Vita Ayoub ◽  
Carole Delenne ◽  
Patrick Matgen ◽  
Pascal Finaud-Guyot ◽  
Renaud Hostache

<p><span>In hydrodynamic modelling, the mesh resolution has a strong impact on run time and result accuracy. Coarser meshes allow faster simulations but often at the cost of accuracy. Conversely, finer meshes offer a better description of complex geometries but require much longer computational time, which makes their use at a large scale challenging. In this context, we aim to assess the potential of a two-dimensional shallow water model with depth-dependant porosity (SW2D-DDP) for flood simulations at a large scale. This modelling approach relies on nesting a sub-grid mesh containing high-resolution topographic and bathymetric data within each computational cell via a so-called depth-dependant storage porosity. It enables therefore faster simulations on rather coarse grids while preserving small-scale topography information. The July 2007 flood event in the Severn River basin (UK) is used as a test case, for which hydrometric measurements and spatial data are available for evaluation. A sensitivity analysis is carried out to investigate the porosity influence on the model performance in comparison with other classical parameters such as boundary conditions.</span></p>

Water ◽  
2021 ◽  
Vol 13 (6) ◽  
pp. 750
Author(s):  
Antonio Pasculli ◽  
Jacopo Cinosi ◽  
Laura Turconi ◽  
Nicola Sciarra

The current climate change could lead to an intensification of extreme weather events, such as sudden floods and fast flowing debris flows. Accordingly, the availability of an early-warning device system, based on hydrological data and on both accurate and very fast running mathematical-numerical models, would be not only desirable, but also necessary in areas of particular hazard. To this purpose, the 2D Riemann–Godunov shallow-water approach, solved in parallel on a Graphical-Processing-Unit (GPU) (able to drastically reduce calculation time) and implemented with the RiverFlow2D code (version 2017), was selected as a possible tool to be applied within the Alpine contexts. Moreover, it was also necessary to identify a prototype of an actual rainfall monitoring network and an actual debris-flow event, beside the acquisition of an accurate numerical description of the topography. The Marderello’s basin (Alps, Turin, Italy), described by a 5 × 5 m Digital Terrain Model (DTM), equipped with five rain-gauges and one hydrometer and the muddy debris flow event that was monitored on 22 July 2016, were identified as a typical test case, well representative of mountain contexts and the phenomena under study. Several parametric analyses, also including selected infiltration modelling, were carried out in order to individuate the best numerical values fitting the measured data. Different rheological options, such as Coulomb-Turbulent-Yield and others, were tested. Moreover, some useful general suggestions, regarding the improvement of the adopted mathematical modelling, were acquired. The rapidity of the computational time due to the application of the GPU and the comparison between experimental data and numerical results, regarding both the arrival time and the height of the debris wave, clearly show that the selected approaches and methodology can be considered suitable and accurate tools to be included in an early-warning system, based at least on simple acoustic and/or light alarms that can allow rapid evacuation, for fast flowing debris flows.


2009 ◽  
Vol 137 (10) ◽  
pp. 3339-3350 ◽  
Author(s):  
Ramachandran D. Nair

Abstract A second-order diffusion scheme is developed for the discontinuous Galerkin (DG) global shallow-water model. The shallow-water equations are discretized on the cubed sphere tiled with quadrilateral elements relying on a nonorthogonal curvilinear coordinate system. In the viscous shallow-water model the diffusion terms (viscous fluxes) are approximated with two different approaches: 1) the element-wise localized discretization without considering the interelement contributions and 2) the discretization based on the local discontinuous Galerkin (LDG) method. In the LDG formulation the advection–diffusion equation is solved as a first-order system. All of the curvature terms resulting from the cubed-sphere geometry are incorporated into the first-order system. The effectiveness of each diffusion scheme is studied using the standard shallow-water test cases. The approach of element-wise localized discretization of the diffusion term is easy to implement but found to be less effective, and with relatively high diffusion coefficients, it can adversely affect the solution. The shallow-water tests show that the LDG scheme converges monotonically and that the rate of convergence is dependent on the coefficient of diffusion. Also the LDG scheme successfully eliminates small-scale noise, and the simulated results are smooth and comparable to the reference solution.


Author(s):  
Feng Jie Zheng ◽  
Fu Zheng Qu ◽  
Xue Guan Song

Reservoir-pipe-valve (RPV) systems are widely used in many industrial process. The pressure in an RPV system plays an important role in the safe operation of the system, especially during the sudden operation such as rapid valve opening/closing. To investigate the pressure especially the pressure fluctuation in an RPV system, a multidimensional and multiscale model combining the method of characteristics (MOC) and computational fluid dynamics (CFD) method is proposed. In the model, the reservoir is modeled by a zero-dimensional virtual point, the pipe is modeled by a one-dimensional MOC, and the valve is modeled by a three-dimensional CFD model. An interface model is used to connect the multidimensional and multiscale model. Based on the model, a transient simulation of the turbulent flow in an RPV system is conducted, in which not only the pressure fluctuation in the pipe but also the detailed pressure distribution in the valve are obtained. The results show that the proposed model is in good agreement with the full CFD model in both large-scale and small-scale spaces. Moreover, the proposed model is more computationally efficient than the CFD model, which provides a feasibility in the analysis of complex RPV system within an affordable computational time.


Author(s):  
Hilary Weller

The shallow water equations are solved using a mesh of polygons on the sphere, which adapts infrequently to the predicted future solution. Infrequent mesh adaptation reduces the cost of adaptation and load-balancing and will thus allow for more accurate mapping on adaptation. We simulate the growth of a barotropically unstable jet adapting the mesh every 12 h. Using an adaptation criterion based largely on the gradient of the vorticity leads to a mesh with around 20 per cent of the cells of a uniform mesh that gives equivalent results. This is a similar proportion to previous studies of the same test case with mesh adaptation every 1–20 min. The prediction of the mesh density involves solving the shallow water equations on a coarse mesh in advance of the locally refined mesh in order to estimate where features requiring higher resolution will grow, decay or move to. The adaptation criterion consists of two parts: that resolved on the coarse mesh, and that which is not resolved and so is passively advected on the coarse mesh. This combination leads to a balance between resolving features controlled by the large-scale dynamics and maintaining fine-scale features.


2018 ◽  
Author(s):  
LMD

We show how the two-layer moist-convective rotating shallow water model (mcRSW), which proved to be a simple and robust tool for studying effects of moist convection on large-scale atmospheric motions, can be improved by including, in addition to the water vapour, precipitable water, and the effects of vaporisation, entrainment, and precipitation. Thus improved mcRSW becomes cloud-resolving. It is applied, as an illustration, to model the development of instabilities of tropical cyclone-like vortices.


2019 ◽  
Vol 141 (5) ◽  
Author(s):  
Feng Jie Zheng ◽  
Chao Yong Zong ◽  
William Dempster ◽  
Fu Zheng Qu ◽  
Xue Guan Song

Reservoir-pipe-valve (RPV) systems are widely used in many industrial processes. The pressure in an RPV system plays an important role in the safe operation of the system, especially during the sudden operations such as rapid valve opening or closing. To investigate the pressure response, with particular interest in the pressure fluctuations in an RPV system, a multidimensional and multiscale model combining the method of characteristics (MOC) and computational fluid dynamics (CFD) method is proposed. In the model, the reservoir is modeled as a zero-dimensional virtual point, the pipe is modeled as a one-dimensional system using the MOC, and the valve is modeled using a three-dimensional CFD model. An interface model is used to connect the multidimensional and multiscale model. Based on the model, a transient simulation of the turbulent flow in an RPV system is conducted in which not only the pressure fluctuation in the pipe but also the detailed pressure distribution in the valve is obtained. The results show that the proposed model is in good agreement when compared with a high fidelity CFD model used to represent both large-scale and small-scale spaces. As expected, the proposed model is significantly more computationally efficient than the CFD model. This demonstrates the feasibility of analyzing complex RPV systems within an affordable computational time.


2013 ◽  
Vol 21 (04) ◽  
pp. 1350017
Author(s):  
RAMIN KAVIANI ◽  
VAHID ESFAHANIAN ◽  
MOHAMMAD EBRAHIMI

The affordable grid resolutions in conventional large-eddy simulations (LESs) of high Reynolds jet flows are unable to capture the sound generated by fluid motions near and beyond the grid cut-off scale. As a result, the frequency spectrum of the extrapolated sound field is artificially truncated at high frequencies. In this paper, a new method is proposed to account for the high frequency noise sources beyond the resolution of a compressible flow simulation. The large-scale turbulent structures as dominant radiators of sound are captured in LES, satisfying filtered Navier–Stokes equations, while for small-scale turbulence, a Kolmogorov's turbulence spectrum is imposed. The latter is performed via a wavelet-based extrapolation to add randomly generated small-scale noise sources to the LES near-field data. Further, the vorticity and instability waves are filtered out via a passive wavelet-based masking and the whole spectrum of filtered data are captured on a Ffowcs-Williams/Hawkings (FW-H) surface surrounding the near-field region and are projected to acoustic far-field. The algorithm can be implemented as a separate postprocessing stage and it is observed that the computational time is considerably reduced utilizing a hybrid of many-core and multi-core framework, i.e. MPI-CUDA programming. The comparison of the results obtained from this procedure and those from experiments for high subsonic and transonic jets, shows that the far-field noise spectrum agree well up to 2 times of the grid cut-off frequency.


Water ◽  
2020 ◽  
Vol 12 (2) ◽  
pp. 439 ◽  
Author(s):  
Sara Venturi ◽  
Silvia Di Francesco ◽  
Martin Geier ◽  
Piergiorgio Manciola

This work compares three forcing schemes for a recently introduced cascaded lattice Boltzmann shallow water model: a basic scheme, a second-order scheme, and a centred scheme. Although the force is applied in the streaming step of the lattice Boltzmann model, the acceleration is also considered in the transformation to central moments. The model performance is tested for one and two dimensional benchmarks.


2020 ◽  
Author(s):  
Giulia Mazzotti ◽  
Richard Essery ◽  
Johanna Malle ◽  
Clare Webster ◽  
Tobias Jonas

<p>Forest canopies strongly affect snowpack energetics during wintertime. In discontinuous forest stands, spatio-temporal variations in radiative and turbulent fluxes create complex snow distribution and melt patterns, with further impacts on the hydrological regimes and on the land surface properties of seasonally snow-covered forested environments.</p><p>As increasingly detailed canopy structure datasets are becoming available, canopy-induced energy exchange processes can be explicitly represented in high-resolution snow models. We applied the modelling framework FSM2 to obtain spatially distributed simulations of the forest snowpack in subalpine and boreal forest stands at high spatial (2m) and temporal (10min) resolution. Modelled sub-canopy radiative and turbulent fluxes were compared to detailed meteorological data of incoming irradiances, air and snow surface temperatures. These were acquired with novel observational systems, including 1) a motorized cable car setup recording spatially and temporally resolved data along a transect and 2) a handheld setup designed to capture temporal snapshots of 2D spatial distributions across forest discontinuities.</p><p>The combination of high-resolution modelling and multi-dimensional datasets allowed us to assess model performance at the level of individual energy balance components, under various meteorological conditions and across canopy density gradients. We showed which canopy representation strategies within FSM2 best succeeded in reproducing snowpack energy transfer dynamics in discontinuous forests, and derived implications for implementing forest snow processes in coarser-resolution models.</p>


2020 ◽  
Vol 77 (9) ◽  
pp. 3119-3137
Author(s):  
Marcin J. Kurowski ◽  
Wojciech W. Grabowski ◽  
Kay Suselj ◽  
João Teixeira

Abstract Idealized large-eddy simulation (LES) is a basic tool for studying three-dimensional turbulence in the planetary boundary layer. LES is capable of providing benchmark solutions for parameterization development efforts. However, real small-scale atmospheric flows develop in heterogeneous and transient environments with locally varying vertical motions inherent to open multiscale interactive dynamical systems. These variations are often too subtle to detect them by state-of-the-art remote and in situ measurements, and are typically excluded from idealized simulations. The present study addresses the impact of weak [i.e., O(10−6) s−1] short-lived low-level large-scale convergence/divergence perturbations on continental shallow convection. The results show a strong response of shallow nonprecipitating convection to the applied weak large-scale dynamical forcing. Evolutions of CAPE, mean liquid water path, and cloud-top heights are significantly affected by the imposed convergence/divergence. In contrast, evolving cloud-base properties, such as the area coverage and mass flux, are only weakly affected. To contrast those impacts with microphysical sensitivity, the baseline simulations are perturbed assuming different observationally based cloud droplet number concentrations and thus different rainfall. For the tested range of microphysical perturbations, the imposed convergence/divergence provides significantly larger impact than changes in the cloud microphysics. Simulation results presented here provide a stringent test for convection parameterizations, especially important for large-scale models progressing toward resolving some nonhydrostatic effects.


Sign in / Sign up

Export Citation Format

Share Document