Fourth order transport model on Yin-Yang grid by multi-moment constrained finite volume scheme

Abstract A fourth order transport model is proposed for global computation with the application of multi-moment constrained finite volume (MCV) scheme and Yin-Yang overset grid. Using multi-moment concept, local degrees of freedom (DOFs) are point-wisely defined within each mesh element to build a cubic spatial reconstruction. The updating formulations for local DOFs are derived by adopting multi moments as constraint conditions, including volume-integrated average (VIA), point value (PV) and first order derivative value (DV). Using Yin-Yang grid eliminates the polar singularities and results in a quasi-uniform mesh over the whole globe. Each component of Yin-Yang grid is a part of the LAT-LON grid, an orthogonal structured grid, where the MCV formulations designed for Cartesian grid can be applied straightforwardly to develop the high order numerical schemes. Proposed MCV model is checked by widely used benchmark tests. The numerical results show that the present model has fourth order accuracy and is competitive to most existing ones.

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

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

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

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

[5]  Nigel Wood,et al.  SLICE‐S: A Semi‐Lagrangian Inherently Conserving and Efficient scheme for transport problems on the Sphere , 2004 .

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

[7]  Paul A. Ullrich,et al.  A conservative semi-Lagrangian multi-tracer transport scheme (CSLAM) on the cubed-sphere grid , 2010, J. Comput. Phys..

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

[9]  Y. Baba,et al.  Dynamical Core of an Atmospheric General Circulation Model on a Yin-Yang Grid , 2010 .

[10]  C. Jablonowski,et al.  Moving Vortices on the Sphere: A Test Case for Horizontal Advection Problems , 2008 .

[11]  F. Xiao,et al.  A Multimoment Finite-Volume Shallow-Water Model On The Yin-Yang Overset Spherical Grid , 2008 .

[12]  Feng Xiao,et al.  A global shallow water model using high order multi-moment constrained finite volume method and icosahedral grid , 2010, J. Comput. Phys..

[13]  Feng Xiao,et al.  High order multi-moment constrained finite volume method. Part I: Basic formulation , 2009, J. Comput. Phys..

[14]  F. Mesinger,et al.  A global shallow‐water model using an expanded spherical cube: Gnomonic versus conformal coordinates , 1996 .

[15]  Feng Xiao,et al.  Conservative constraint for a quasi‐uniform overset grid on the sphere , 2006 .