Discontinuous Galerkin scheme for the spherical shallow water equations with applications to tsunami modeling and prediction

JOURNAL OF COMPUTATIONAL PHYSICS

卷:362

Article

Bonev, Boris¹  Hesthaven, Jan S.¹  Giraldo, Francis X.²  Kopera, Michal A.³ 

[1] Ecole Polytech Fed Lausanne, CH-1015 Lausanne, Switzerland

[2] Naval Postgrad Sch, Monterey, CA 93943 USA

[3] Univ Calif Santa Cruz, Santa Cruz, CA 95064 USA

关键词: Shallow water equations;    Well-balanced scheme;    Wetting and drying;    Discontinuous Galerkin method;    Curved mesh;    Non-conforming mesh;   

DOI  :  10.1016/j.jcp.2018.02.008

来源: Elsevier

【 摘 要 】

We present a novel high-order discontinuous Galerkin discretization for the spherical shallow water equations, able to handle wetting/drying and non-conforming, curved meshes in a well-balanced manner. This requires a well-balanced discretization, that cannot rely on exact quadrature, due to the curved mesh. Using the strong form of the discontinuous Galerkin discretization, we achieve a splitting of the well-balanced condition into individual problems for the flux and volume terms, which has significant advantages: It allows for the construction of non-conforming, well-balanced flux discretizations, i.e. we can perform non-conforming mesh refinement while preserving the well-balanced property of the scheme. More importantly, this approach enables the development of a new method for handling wet/dry transitions. In contrast to other wetting/drying methods, it is well-balanced and able to handle wetting/drying robustly at any polynomial order, without the introduction of physical model assumptions such as viscosity, artificial porosity or cancellation of gravity. We perform a series of one-dimensional tests and analyze the properties of our scheme. In order to validate our method for the simulation of large-scale tsunami events on the rotating sphere, we perform numerical simulations of the 2011 Tohoku tsunami and compare our results to real-world buoy data. The method is able to predict arrival times and wave amplitudes accurately even over long distances. This indicates that our method accurately captures all physical phenomena relevant to the long-term evolution of tsunami waves. (C) 2018 Elsevier Inc. All rights reserved.

【 授权许可】

Free