Well-balanced methods for the shallow water equations in spherical coordinates

Abstract The goal of this work is to obtain a family of explicit high order well-balanced methods for the shallow water equations in spherical coordinates. Application of shallow water models to large scale problems requires the use of spherical coordinates: this is the case, for instance, of the simulation of the propagation of a Tsunami wave through the ocean. Although the PDE system is similar to the shallow water equations in Cartesian coordinates, new source terms appear. As a consequence, the derivation of high order numerical methods that preserve water at rest solutions is not as straightforward as in that case. Finite Volume methods are considered based on a first order path-conservative scheme and high order reconstruction operators. Numerical methods based on these ingredients have been successfully applied previously to the nonlinear SWEs in Cartesian coordinates. Some numerical tests to check the well-balancing and high order properties of the scheme, as well as its ability to simulate planetary waves or tsunami waves over realistic bathymetries are presented.

[1]  Michael Dumbser,et al.  Finite volume schemes of very high order of accuracy for stiff hyperbolic balance laws , 2008, J. Comput. Phys..

[2]  Marc de la Asunción,et al.  Efficient GPU implementation of a two waves TVD-WAF method for the two-dimensional one layer shallow water system on structured meshes , 2013 .

[3]  Manuel Jesús Castro Díaz,et al.  High Order Extensions of Roe Schemes for Two-Dimensional Nonconservative Hyperbolic Systems , 2009, J. Sci. Comput..

[4]  E. Toro,et al.  Restoration of the contact surface in the HLL-Riemann solver , 1994 .

[5]  Enrique D. Fernández Nieto,et al.  A family of stable numerical solvers for the shallow water equations with source terms , 2003 .

[6]  Manuel Jesús Castro Díaz,et al.  High order finite volume schemes based on reconstruction of states for solving hyperbolic systems with nonconservative products. Applications to shallow-water systems , 2006, Math. Comput..

[7]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[8]  Yulong Xing,et al.  High-order well-balanced finite volume WENO schemes for shallow water equation with moving water , 2007, J. Comput. Phys..

[9]  Manuel Jesús Castro Díaz,et al.  Well-Balanced High Order Extensions of Godunov's Method for Semilinear Balance Laws , 2008, SIAM J. Numer. Anal..

[10]  F. Bouchut.  Nonlinear Stability of Finite Volume Methods for Hyperbolic Conservation Laws: and Well-Balanced Schemes for Sources , 2005 .

[11]  Carlos Parés,et al.  On the well-balance property of Roe?s method for nonconservative hyperbolic systems , 2004 .

[12]  Manuel Jesús Castro Díaz,et al.  A Class of Computationally Fast First Order Finite Volume Solvers: PVM Methods , 2012, SIAM J. Sci. Comput..

[13]  Randall J. LeVeque,et al.  Balancing Source Terms and Flux Gradients in High-Resolution Godunov Methods , 1998 .

[14]  Yulong Xing,et al.  High order well-balanced finite volume WENO schemes and discontinuous Galerkin methods for a class of hyperbolic systems with source terms , 2006, J. Comput. Phys..

[15]  Alfredo Bermúdez,et al.  Upwind methods for hyperbolic conservation laws with source terms , 1994 .

[16]  Matteo Semplice,et al.  On the Accuracy of WENO and CWENO Reconstructions of Third Order on Nonuniform Meshes , 2015, Journal of Scientific Computing.

[17]  Mária Lukácová-Medvid'ová,et al.  Well-balanced finite volume evolution Galerkin methods for the shallow water equations , 2015, J. Comput. Phys..

[18]  Y. Okada.  Surface deformation due to shear and tensile faults in a half-space , 1985 .

[19]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[20]  Manuel Jesús Castro Díaz,et al.  A second order PVM flux limiter method. Application to magnetohydrodynamics and shallow stratified flows , 2014, J. Comput. Phys..

[21]  C. Parés.  Numerical methods for nonconservative hyperbolic systems: a theoretical framework. , 2006 .

[22]  Eleuterio F. Toro,et al.  Centred TVD schemes for hyperbolic conservation laws , 2000 .

[23]  José Miguel Mantas,et al.  Two-Dimensional Compact Third-Order Polynomial Reconstructions. Solving Nonconservative Hyperbolic Systems Using GPUs , 2011, J. Sci. Comput..

[24]  Emmanuel Audusse,et al.  A Fast and Stable Well-Balanced Scheme with Hydrostatic Reconstruction for Shallow Water Flows , 2004, SIAM J. Sci. Comput..

[25]  José M. Mantas,et al.  An Introduction to GPU Computing for Numerical Simulation , 2016 .

[26]  Manuel Jesús Castro Díaz,et al.  On a well-balanced high-order finite volume scheme for shallow water equations with topography and dry areas , 2007, J. Comput. Phys..

[27]  Y. Xing.  Numerical Methods for the Nonlinear Shallow Water Equations , 2017 .

[28]  P. Lax,et al.  On Upstream Differencing and Godunov-Type Schemes for Hyperbolic Conservation Laws , 1983 .

[29]  Eleuterio F. Toro,et al.  A weighted average flux method for hyperbolic conservation laws , 1989, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[30]  Antonio Marquina,et al.  Local Piecewise Hyperbolic Reconstruction of Numerical Fluxes for Nonlinear Scalar Conservation Laws , 1994, SIAM J. Sci. Comput..

[31]  Jostein R. Natvig,et al.  Well-balanced finite volume schemes of arbitrary order of accuracy for shallow water flows , 2006, J. Comput. Phys..

[32]  Kun Xu,et al.  A gas-kinetic scheme for shallow-water equations with source terms , 2004 .

[33]  Carlos Parés,et al.  Well-Balanced Schemes and Path-Conservative Numerical Methods , 2017 .

[34]  Tan Weiyan,et al.  Shallow Water Hydrodynamics: Mathematical Theory and Numerical Solution for a Two-dimensional System of Shallow Water Equations , 2012 .

[35]  Michael Dumbser,et al.  Building Blocks for Arbitrary High Order Discontinuous Galerkin Schemes , 2006, J. Sci. Comput..

[36]  Chi-Wang Shu,et al.  Total variation diminishing Runge-Kutta schemes , 1998, Math. Comput..

[37]  S. Osher,et al.  Regular ArticleUniformly High Order Accurate Essentially Non-oscillatory Schemes, III , 1997 .