Improved multislope MUSCL reconstruction on unstructured grids for shallow water equations

In shallow water flow and transport modeling, the monotone upstream-centered scheme for conservation laws (MUSCL) is widely used to extend the original Godunov scheme to second-order accuracy. The most important step in MUSCL-type schemes is the MUSCL reconstruction, which calculate extrapolates the values of independent variables from the cell center to the edge. The monotonicity of the scheme is preserved with the help of slope limiters that prevent the occurrence of new extrema during the reconstruction. On structured grids, the calculation of the slope is straightforward and usually based on a two-point stencil that uses the cell centers of the neighbor cell and the so-called far-neighbor cell of the edge under consideration. On unstructured grids, the correct choice for the upwind slope becomes non-trivial. In this work, two novel TVD schemes are developed based on different techniques for calculating the upwind slope and downwind slope. An additional treatment that stabilizes the scheme is discussed. The proposed techniques are compared to two existing MUSCL reconstruction techniques and a detailed discussion of the results is given. It is shown that the proposed MUSCL reconstruction schemes obtain more accurate results with less numerical diffusion and higher efficiency

An efficient unstructured MUSCL scheme for solving the 2D shallow water equations

The aim of this paper is to present a novel monotone upstream scheme for conservation law (MUSCL) on unstructured grids. The novel edge-based MUSCL scheme is devised to construct the required values at the midpoint of cell edges in a more straightforward and effective way compared to other conventional approaches, by making better use of the geometrical property of the triangular grids. The scheme is incorporated into a two-dimensional (2D) cell-centered Godunov-type finite volume model as proposed in Hou et al. (2013a,c) to solve the shallow water equations (SWEs). The MUSCL scheme renders the model to preserve the well-balanced property and achieve high accuracy and efficiency for shallow flow simulations over uneven terrains. Furthermore, the scheme is directly applicable to all triangular grids. Application to several numerical experiments verifies the efficiency and robustness of the current new MUSCL scheme

A depth-averaged non-cohesive sediment transport model with improved discretization of flux and source terms

This paper presents novel flux and source term treatments within a Godunov-type finite volume framework for predicting the depth-averaged shallow water fl ow and sediment transport with enhanced the accuracy and stability. The suspended load ratio is introduced to differentiate between the advection of the suspended load and the advection of water. A modified Harten, Lax and van Leer Riemann solver with the contact wave restored (HLLC) is derived for the fl ux calculation based on the new wave pattern involving the suspended load ratio. The source term calculation is enhanced by means of a novel splitting-point implicit discretization. The slope effect is introduced by modifying the critical shear stress, with two treatments being discussed. The numerical scheme is tested in five examples that comprise both fixed and movable beds. The model predictions show good agreement with measurement, except for cases where local three-dimensional effects dominate.China Scholarship Counci

Second-order discontinuous Galerkin flood model: comparison with industry-standard finite volume models

Finite volume (FV) numerical solvers of the two-dimensional shallow water equations are core to industry-standard flood models. The second-order Discontinuous Galerkin (DG) alternative is well-known to perform better than first- and second-order FV to capture sharp flow fronts and converge faster at coarser resolutions, but DG2 models typically rely on local slope limiting to selectively damp numerical oscillations in the vicinity of shock waves. Yet flood inundation events are smooth and gradually-varying, and shock waves play only a minor role in flood inundation modelling. Therefore, this paper investigates two DG2 variants - with and without local slope limiting - to identify the simplest and most efficient DG2 configuration suitable for flood inundation modelling. The predictive capabilities of the DG2 variants are analysed for a synthetic test case involving advancing and receding waves representative of flood-like flow. The DG2 variants are then benchmarked against industry-standard FV models over six UK Environment Agency scenarios. Results indicate that the DG2 variant without local slope limiting closely reproduces solutions of the commercial models at twice as coarse a spatial resolution, and removing the slope limiter can halve model runtime. Results also indicate that DG2 can capture more accurate hydrographs incorporating small-scale transients over long-range simulations, even when hydrographs are measured far away from the flooding source. Accompanying details of software and data accessibility are provided

Two-dimensional discontinuous Galerkin shallow water model for practical flood modelling applications

Finite volume (FV) numerical solvers to the two-dimensional shallow water equations are the foundation of the current state-of-the-practice, industry-standard flood models. The second-order Discontinuous Galerkin (DG2) alternative show a promising way to improve current FV-based flood model formulations, but is yet under-studied and rarely utilised to support flood modelling applications. This is contributed by the mathematical complexity constructed within the DG2 formulation that could lead to large computational costs and compromise its stability and robustness when used for practical modelling. Therefore, this PhD research aims to develop a new flood model based on simplified DG2 solver that is improved for flood modelling practices. To achieve this aim, three objectives have been formed and addressed through analyses involving academic and experimental test cases, as well as test cases that are recommended by the UK Environment Agency to validate 2D flood model capabilities, whilst benchmarking the simplified DG2 solver against four FV-based industrial models. Key research findings indicate that the simplified DG2 solver can equally retain conservative properties and provide second-order accurate predictions as the standard DG2 solver whilst offering around 2.6 times runtime speed up. Additionally, the simplified DG2 solver can be reliably efficient to provide predictions close to the outputs of the industrial models, in simulating flood scenarios covering large catchment-scale areas and at a grid resolution greater or equal than 5 m, particularly when the local limiting is disabled. However, the local limiting is still needed by the simplified DG2 solver when modelling detailed velocity fields at sub-metre grid resolutions, particularly in regions of highly active wave-structure interactions as commonly encountered in urban flooding around steep-sloped building structures

SERGHEI (SERGHEI-SWE) v1.0: a performance-portable high-performance parallel-computing shallow-water solver for hydrology and environmental hydraulics

The Simulation EnviRonment for Geomorphology, Hydrodynamics, and Ecohydrology in Integrated form (SERGHEI) is a multi-dimensional, multi-domain, and multi-physics model framework for environmental and landscape simulation, designed with an outlook towards Earth system modelling. At the core of SERGHEI's innovation is its performance-portable high-performance parallel-computing (HPC) implementation, built from scratch on the Kokkos portability layer, allowing SERGHEI to be deployed, in a performance-portable fashion, in graphics processing unit (GPU)-based heterogeneous systems. In this work, we explore combinations of MPI and Kokkos using OpenMP and CUDA backends. In this contribution, we introduce the SERGHEI model framework and present with detail its first operational module for solving shallow-water equations (SERGHEI-SWE) and its HPC implementation. This module is designed to be applicable to hydrological and environmental problems including flooding and runoff generation, with an outlook towards Earth system modelling. Its applicability is demonstrated by testing several well-known benchmarks and large-scale problems, for which SERGHEI-SWE achieves excellent results for the different types of shallow-water problems. Finally, SERGHEI-SWE scalability and performance portability is demonstrated and evaluated on several TOP500 HPC systems, with very good scaling in the range of over 20 000 CPUs and up to 256 state-of-the art GPUs

A robust high-resolution hydrodynamic numerical model for surface water flow and transport processes within a flexible software framework

Paralleltitel: Ein robustes hochauflösendes hydrodynamisch-numerisches Modell für Oberflächenabfluss- und Transportprozesse innerhalb eines flexiblen Software-Framework

La méthode MOOD Multi-dimensional Optimal Order Detection : la première approche a posteriori aux méthodes volumes finis d'ordre très élevé

Nous introduisons et développons dans cette thèse un nouveau type de méthodes Volumes Finis d'ordre très élevé pour les systèmes hyperboliques de lois de conservations. Appelée MOOD pour Multidimensional Optimal Order Detection, elle permet de réaliser des simulations très précises en dimensions deux et trois sur maillages non-structurés. La conception d'une telle méthode est rendue délicate par l'apparition de singularités dans la solution (chocs, discontinuités de contact) pour lesquelles des phenomènes parasites (oscillations, création de valeurs non physiques...) sont générés par l'approximation d'ordre élevé. L'originalité de cette thèse réside dans le traitement de ces problèmes. A l'opposé des méthodes classiques qui essaient de contrôler ces phénomènes indésirables par une limitation a priori, nous proposons une approche de traitement a posteriori basée sur une décrémentation locale de l'ordre du schéma. Nous montrons en particulier que ce concept permet très simplement d'obtenir des propriétés qui sont habituellement difficiles à prouver dans le cadre multi-dimensionel non-structuré (préservation de la positité par exemple). La robustesse et la qualité de la méthode MOOD ont été prouvées sur de nombreux tests numériques en 2D et 3D. Une amélioration significative des ressources informatiques (CPU et stockage mémoire) nécessaires à l'obtention de résultats équivalents aux méthodes actuelles a été démontrée.We introduce and develop in this thesis a new type of very high-order Finite Volume methods for hyperbolic systems of conservation laws. This method, named MOOD for Multidimensional Optimal Order Detection, provides very accurate simulations for two- and three-dimensional unstructured meshes. The design of such a method is made delicate by the emergence of solution singularities (shocks, contact discontinuities) for which spurious phenomena (oscillations, non-physical values creation, etc.) are generated by the high-order approximation. The originality of this work lies in a new treatment for theses problems. Contrary to classical methods which try to control such undesirable phenomena through an a priori limitation, we propose an a posteriori treatment approach based on a local scheme order decrementing. In particular, we show that this concept easily provides properties that are usually difficult to prove in a multidimensional unstructured framework (positivity-preserving for instance). The robustness and quality of the MOOD method have been numerically proved through numerous test cases in 2D and 3D, and a significant reduction of computational resources (CPU and memory storage) needed to get state-of-the-art results has been shown

Schémas numériques pour la simulation de l'explosion

In nuclear facilities, internal or external explosions can cause confinement breaches and radioactive materials release in the environment. Hence, modeling such phenomena is crucial for safety matters. The purpose of this thesis is to contribute to the creation of efficient numerical schemes to solve these complex models. The work presented here focuses on two major aspects: first, the development of consistent schemes for the Euler equations which model the blast waves, then the buildup of reliable schemes for the front propagation, like the flame front during the deflagration phenomenon. Staggered discretization is used in space for all the schemes. It is based on the internal energy formulation of the Euler system, which insures its positivity and the positivity of the density. A discrete kinetic energy balance is derived from the scheme and a source term is added in the discrete internal energy balance equation to preserve the exact total energy balance. High order, MUSCL-like interpolators are used in the discrete momentum operators. The resulting scheme is consistent (in the sense of Lax) with the weak entropic solutions of the continuous problem. We use the properties of Hamilton-Jacobi equations to build a class of finite volume schemes compatible with a large number of meshes to model the flame front propagation. These schemes satisfy a maximum principle and have important consistency and monotonicity properties. These latter allow to derive a convergence result for the schemes based on Cartesian grids.Dans les installations nucléaires, les explosions, qu’elles soient d’origine interne ou externe, peuvent entrainer la rupture du confinement et le rejet de matières radioactives dans l’environnement. Il est donc fondamental, dans un cadre de sûreté de modéliser ce phénomène. L’objectif de cette thèse est de contribuer à l’élaboration de schémas numériques performants pour résoudre ces modèles complexes. Les travaux présentés s’articule autour de deux axes majeurs : le développement de schémas volumes finis consistants pour les équations d’Euler compressible qui modélise les ondes de choc et celui de schémas performants pour la propagation d’interfaces comme le front de flamme lors d'une déflagration. La discrétisation spatiale est de type mailles décalées pour tous les schémas développés. Les schémas pour les équations d'Euler se basent sur une formulation en énergie interne qui permet de préserver sa positivité ainsi que celle de la masse volumique. Un bilan d'énergie cinétique discret peut être obtenu et permet de retrouver un bilan d'énergie totale par l'ajout d'un terme de correction dans le bilan d'énergie interne. Le schéma ainsi construit est consistant au sens de Lax avec les solutions faibles entropiques des équations continues. On utilise les propriétés des équations de type Hamilton-Jacobi pour construire une classe de schémas volumes finis performants sur une large variété de maillages modélisant la propagation du front de flamme. Ces schémas garantissent un principe du maximum et possèdent des propriétés importantes de monotonie et consistance qui permettent d'obtenir un résultat de convergence