A parallel well-balanced finite volume method for shallow water equations with topography on the cubed-sphere

A finite volume scheme for the global shallow water model on the cubed-sphere mesh is proposed and studied in this paper. The new cell-centered scheme is based on Osher's Riemann solver together with a high-order spatial reconstruction. On each patch interface of the cubed-sphere only one layer of ghost cells is needed in the scheme and the numerical flux is calculated symmetrically across the interface to ensure the numerical conservation of total mass. The discretization of the topographic term in the equation is properly modified in a well-balanced manner to suppress spurious oscillations when the bottom topography is non-smooth. Numerical results for several test cases including a steady-state nonlinear geostrophic flow and a zonal flow over an isolated mountain are provided to show the flexibility of the scheme. Some parallel implementation details as well as some performance results on a parallel supercomputer with more than one thousand processor cores are also provided.

[1]  Shian-Jiann Lin,et al.  Finite-volume transport on various cubed-sphere grids , 2007, J. Comput. Phys..

[2]  Akio Arakawa,et al.  Integration of the Nondivergent Barotropic Vorticity Equation with AN Icosahedral-Hexagonal Grid for the SPHERE1 , 1968 .

[3]  A. Kageyama,et al.  ``Yin-Yang grid'': An overset grid in spherical geometry , 2004, physics/0403123.

[4]  R. Sadourny.  Conservative Finite-Difference Approximations of the Primitive Equations on Quasi-Uniform Spherical Grids , 1972 .

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

[6]  J. Hack,et al.  Spectral transform solutions to the shallow water test set , 1995 .

[7]  Stephen J. Thomas,et al.  Semi-Implicit Spectral Element Atmospheric Model , 2002, J. Sci. Comput..

[8]  Jan S. Hesthaven,et al.  Nodal high-order discontinuous Galerkin methods for the spherical shallow water equations , 2002 .

[9]  S. Osher,et al.  Upwind difference schemes for hyperbolic systems of conservation laws , 1982 .

[10]  James A. Rossmanith,et al.  A wave propagation method for hyperbolic systems on the sphere , 2006, J. Comput. Phys..

[11]  J. Greenberg,et al.  A well-balanced scheme for the numerical processing of source terms in hyperbolic equations , 1996 .

[12]  N SwarztrauberPaul,et al.  A standard test set for numerical approximations to the shallow water equations in spherical geometry , 1992 .

[13]  John M. Dennis,et al.  Towards an efficient and scalable discontinuous Galerkin atmospheric model , 2005, 19th IEEE International Parallel and Distributed Processing Symposium.

[14]  Chi-Wang Shu,et al.  Strong Stability-Preserving High-Order Time Discretization Methods , 2001, SIAM Rev..

[15]  H. Tufo,et al.  Development of a Scalable Global Discontinuous Galerkin Atmospheric Model , 2005 .

[16]  P. Paolucci,et al.  The “Cubed Sphere” , 1996 .

[17]  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..

[18]  Stephen J. Thomas,et al.  A Discontinuous Galerkin Global Shallow Water Model , 2005, Monthly Weather Review.

[19]  P. Swarztrauber,et al.  A standard test set for numerical approximations to the shallow water equations in spherical geometry , 1992 .

[20]  David L. Williamson,et al.  Integration of the barotropic vorticity equation on a spherical geodesic grid , 1968 .

[21]  William Gropp,et al.  Newton-Krylov-Schwarz Methods in CFD , 1994 .

[22]  Mark A. Taylor,et al.  High-Resolution Mesh Convergence Properties and Parallel Efficiency of a Spectral Element Atmospheric Dynamical Core , 2005, Int. J. High Perform. Comput. Appl..

[23]  Feng Xiao,et al.  Shallow water model on cubed-sphere by multi-moment finite volume method , 2008, J. Comput. Phys..

[24]  Stanley Osher,et al.  Upwind schemes and boundary conditions with applications to Euler equations in general geometries , 1983 .

[25]  William Putman,et al.  The finite-volume dynamical core on the cubed-sphere , 2006, SC.

[26]  Chao Yang,et al.  A Fully Implicit Domain Decomposition Algorithm for Shallow Water Equations on the Cubed-Sphere , 2010, SIAM J. Sci. Comput..