NAJul 4, 2016
A large time-step and well-balanced Lagrange-Projection type scheme for the shallow-water equationsChristophe Chalons, Pierre Kestener, Samuel Kokh et al.
This work focuses on the numerical approximation of the Shallow Water Equations (SWE) using a Lagrange-Projection type approach. We propose to extend to this context recent implicit-explicit schemes developed in the framework of compressibleflows, with or without stiff source terms. These methods enable the use of time steps that are no longer constrained by the sound velocity thanks to an implicit treatment of the acoustic waves, and maintain accuracy in the subsonic regime thanks to an explicit treatment of the material waves. In the present setting, a particular attention will be also given to the discretization of the non-conservative terms in SWE and more specifically to the well-known well-balanced property. We prove that the proposed numerical strategy enjoys important non linear stability properties and we illustrate its behaviour past several relevant test cases.
NASep 22, 2017
Experimenting with the p4est library for AMR simulations of two-phase flowsFlorence Drui, Alexandru Fikl, Pierre Kestener et al.
Many physical problems involve spatial and temporal inhomogeneities that require a very fine discretization in order to be accurately simulated. Using an adaptive mesh, a high level of resolution is used in the appropriate areas while keeping a coarse mesh elsewhere. This idea allows to save time and computations, but represents a challenge for distributed-memory environments. The MARS project (for Multiphase Adaptative Refinement Solver) intends to assess the parallel library p4est for adaptive mesh, in a case of a finite volume scheme applied to two-phase flows. Besides testing the library's performances, particularly for load balancing, its user-friendliness in use and implementation are also exhibited here. First promising 3D simulations are even presented.
NAFeb 4, 2019
An all-regime and well-balanced Lagrange-projection type scheme for the shallow water equations on unstructured meshesChristophe Chalons, Samuel Kokh, Maxime Stauffert
In this work, we focus on the numerical approximation of the shallow water equations in two space dimensions. Our aim is to propose a well-balanced, all-regime and positive scheme. By well-balanced, it is meant that the scheme is able to preserve the so-called lake at rest smooth equilibrium solutions. By all-regime, we mean that the scheme is able to deal with all flow regimes, including the low-Froude regime which is known to be challenging when using usual Godunov-type finite volume schemes. At last, the scheme should be positive which means that the water height stays positive for all time. Our approach is based on a Lagrange-projection decomposition which allows to naturally decouple the acoustic and transport terms. Numerical experiments on unstructured meshes illustrate the good behaviour of the scheme.