Implicit kinetic schemes for the Saint-Venant system
EL HASSANIEH, Chourouk
Laboratoire Jacques-Louis Lions [LJLL (UMR_7598)]
Numerical Analysis, Geophysics and Ecology [ANGE]
Laboratoire Jacques-Louis Lions [LJLL (UMR_7598)]
Numerical Analysis, Geophysics and Ecology [ANGE]
SAINTE-MARIE, Jacques
Laboratoire Jacques-Louis Lions [LJLL (UMR_7598)]
Numerical Analysis, Geophysics and Ecology [ANGE]
Laboratoire Jacques-Louis Lions [LJLL (UMR_7598)]
Numerical Analysis, Geophysics and Ecology [ANGE]
EL HASSANIEH, Chourouk
Laboratoire Jacques-Louis Lions [LJLL (UMR_7598)]
Numerical Analysis, Geophysics and Ecology [ANGE]
Laboratoire Jacques-Louis Lions [LJLL (UMR_7598)]
Numerical Analysis, Geophysics and Ecology [ANGE]
SAINTE-MARIE, Jacques
Laboratoire Jacques-Louis Lions [LJLL (UMR_7598)]
Numerical Analysis, Geophysics and Ecology [ANGE]
< Reduce
Laboratoire Jacques-Louis Lions [LJLL (UMR_7598)]
Numerical Analysis, Geophysics and Ecology [ANGE]
Language
en
Document de travail - Pré-publication
This item was published in
2023-03-28
English Abstract
Explicit (in time) kinetic schemes applied to the nonlinear shallow water equations have been extensively studied in the past. The novelty of this paper is to investigate an implicit version of such methods in order to ...Read more >
Explicit (in time) kinetic schemes applied to the nonlinear shallow water equations have been extensively studied in the past. The novelty of this paper is to investigate an implicit version of such methods in order to improve their stability properties. In the case of a flat bathymetry we obtain a fully implicit kinetic solver satisfying a discrete entropy inequality and keeping the water height non negative without any restriction on the time step. Remarkably, a simplified version of this nonlinear implicit scheme allows to express the update explicitly which we implement in practice. The case of varying bottoms is then dealt with through an iterative solver combined with the hydrostatic reconstruction technique. We show that this scheme preserves the water height non-negativity under a CFL condition and satisfies a discrete entropy inequality without error term, which is an improvement over its explicit version. An extension of the implicit and iterative methods to the two dimensional case is also discussed. Finally we perform some numerical validations underlining the advantages and the computational cost of our strategy.Read less <
English Keywords
Shallow water equations
Kinetic solver
Fully discrete entropy inequality
Well-balanced schemes
Hydrostatic reconstruction
Origin
Hal imported