An all Froude high order IMEX scheme for the shallow water equations on unstructured Voronoi meshes

被引:7
|
作者
Boscheri, Walter [1 ]
Tavelli, Maurizio [2 ]
Castro, Cristobal E. [3 ]
机构
[1] Univ Ferrara, Dept Math & Comp Sci, Ferrara, Italy
[2] Free Univ Bozen, Fac Comp Sci, Bolzano, Italy
[3] Univ Tarapaca, Fac Ingeniena, Dept Ingeniena Mecan, Arica, Chile
关键词
IMEX schemes; Finite volume and discontinuous Galerkin; methods; High order in space and time; Asymptotic preserving; Shallow water equations; All Froude flows; NAVIER-STOKES EQUATIONS; DISCONTINUOUS GALERKIN METHOD; HYPERBOLIC CONSERVATION-LAWS; RUNGE-KUTTA SCHEMES; FINITE-DIFFERENCE METHODS; GODUNOV-TYPE SCHEMES; CENTRAL WENO SCHEMES; FREE-SURFACE FLOWS; MACH NUMBER LIMIT; VELOCITY RECONSTRUCTION;
D O I
10.1016/j.apnum.2022.11.022
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
We propose a novel numerical method for the solution of the shallow water equations in different regimes of the Froude number making use of general polygonal meshes. The fluxes of the governing equations are split such that advection and acoustic-gravity sub-systems are derived, hence separating slow and fast phenomena. This splitting allows the nonlinear convective fluxes to be discretized explicitly in time, while retaining an implicit time marching for the acoustic-gravity terms. Consequently, the novel schemes are particularly well suited in the low Froude limit of the model, since no numerical viscosity is added in the implicit solver. Besides, stability follows from a milder CFL condition which is based only on the advection speed and not on the celerity. High order time accuracy is achieved using the family of semi-implicit IMEX Runge-Kutta schemes, while high order in space is granted relying on two discretizations: (i) a cell-centered finite volume (FV) scheme for the nonlinear convective contribution on the polygonal cells; (ii) a staggered discontinuous Galerkin (DG) scheme for the solution of the linear system associated to the implicit discretization of the pressure sub-system. Therefore, three different meshes are used, namely a polygonal Voronoi mesh, a triangular subgrid and a staggered quadrilateral subgrid. The novel schemes are proved to be Asymptotic Preserving (AP), hence a consistent discretization of the limit model is retrieved for vanishing Froude numbers, which is given by the so-called "lake at rest" equations. Furthermore, the novel methods are well-balanced by construction, and this property is also demonstrated. Accuracy and robustness are then validated against a set of benchmark test cases with Froude numbers ranging in the interval Fr approximate to [10-6; 5], hence showing that multiple time scales can be handled by the novel methods.(c) 2022 IMACS. Published by Elsevier B.V. All rights reserved.
引用
收藏
页码:311 / 335
页数:25
相关论文
共 50 条
  • [1] A flux-vector splitting scheme for the shallow water equations extended to high-order on unstructured meshes
    Toro, Eleuterio F.
    Castro, Cristobal E.
    Vanzo, Davide
    Siviglia, Annunziato
    INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, 2022, 94 (10) : 1679 - 1705
  • [2] The FVC Scheme on Unstructured Meshes for the Two-Dimensional Shallow Water Equations
    Ziggaf, Moussa
    Boubekeur, Mohamed
    Kissami, Imad
    Benkhaldoun, Fayssal
    El Mahi, Imad
    FINITE VOLUMES FOR COMPLEX APPLICATIONS IX-METHODS, THEORETICAL ASPECTS, EXAMPLES, FVCA 9, 2020, 323 : 455 - 465
  • [3] Asymptotic preserving scheme for the shallow water equations with source terms on unstructured meshes
    Duran, A.
    Marche, F.
    Turpault, R.
    Berthon, C.
    JOURNAL OF COMPUTATIONAL PHYSICS, 2015, 287 : 184 - 206
  • [4] Asymptotic Consistency of the RS-IMEX Scheme for the Low-Froude Shallow Water Equations: Analysis and Numerics
    Zakerzadeh, Hamed
    THEORY, NUMERICS AND APPLICATIONS OF HYPERBOLIC PROBLEMS II, 2018, 237 : 665 - 675
  • [5] An All Speed Second Order IMEX Relaxation Scheme for the Euler Equations
    Thomann, Andrea
    Zenk, Markus
    Puppo, Gabriella
    Klingenberg, Christian
    COMMUNICATIONS IN COMPUTATIONAL PHYSICS, 2020, 28 (02) : 591 - 620
  • [6] ADER scheme on unstructured meshes for shallow water: simulation of tsunami waves
    Castro, Cristobal E.
    Toro, Eleuterio F.
    Kaeser, Martin
    GEOPHYSICAL JOURNAL INTERNATIONAL, 2012, 189 (03) : 1505 - 1520
  • [7] NUMERICAL DISCRETIZATIONS FOR SHALLOW WATER EQUATIONS WITH SOURCE TERMS ON UNSTRUCTURED MESHES
    Duran, Arnaud
    Marche, Fabien
    Berthon, Christophe
    Turpault, Rodolphe
    HYPERBOLIC PROBLEMS: THEORY, NUMERICS, APPLICATIONS, 2014, 8 : 541 - 549
  • [8] Well-balanced high-order centered schemes on unstructured meshes for shallow water equations with fixed and mobile bed
    Canestrelli, Alberto
    Dumbser, Michael
    Siviglia, Annunziato
    Toro, Eleuterio F.
    ADVANCES IN WATER RESOURCES, 2010, 33 (03) : 291 - 303
  • [9] A high order semi-implicit discontinuous Galerkin method for the two dimensional shallow water equations on staggered unstructured meshes
    Tavelli, Maurizio
    Dumbser, Michael
    APPLIED MATHEMATICS AND COMPUTATION, 2014, 234 : 623 - 644
  • [10] High order finite volume schemes with IMEX time stepping for the Boltzmann model on unstructured meshes
    Boscheri, Walter
    Dimarco, Giacomo
    COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2021, 387