An Upwind-Biased Transport Scheme Using a Quadratic Reconstruction on Spherical Icosahedral Grids

AbstractSeveral transport schemes developed for spherical icosahedral grids are based on the piecewise linear approximation. The simplest one among them uses an algorithm where the tracer distribution in the upwind side of a cell face is reconstructed using a linear surface. Recently, it was demonstrated that using second- or fourth-order reconstructions instead of the linear one produces better results. The computational cost of the second-order reconstruction method was not much larger than the linear one, while that of the fourth-order one was significantly larger. In this work, the authors propose another second-order reconstruction scheme on the spherical icosahedral grids, motivated by some ideas from the piecewise parabolic method. The second-order profile of a tracer is reconstructed under two constraints: (i) the area integral of the profile is equal to the cell-averaged value times the cell area and (ii) the profile is the least squares fit to the cell-vertex values. The new scheme [the second u...

[1]  G. R. Stuhne,et al.  Vortex Erosion and Amalgamation in a New Model of Large Scale Flow on the Sphere , 1996 .

[2]  Masaki Satoh,et al.  Nonhydrostatic icosahedral atmospheric model (NICAM) for global cloud resolving simulations , 2008, J. Comput. Phys..

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

[4]  Todd D. Ringler,et al.  Modeling the Atmospheric General Circulation Using a Spherical Geodesic Grid: A New Class of Dynamical Cores , 2000 .

[5]  H. Miura.  An Upwind-Biased Conservative Advection Scheme for Spherical Hexagonal–Pentagonal Grids , 2007 .

[6]  Peter H. Lauritzen,et al.  A class of deformational flow test cases for linear transport problems on the sphere , 2010, J. Comput. Phys..

[7]  B. P. Leonard,et al.  Positivity-preserving numerical schemes for multidimensional advection , 1993 .

[8]  J. Thuburn.  Dissipation and Cascades to Small Scales in Numerical Models Using a Shape-Preserving Advection Scheme , 1995 .

[9]  B. V. Leer,et al.  Towards the Ultimate Conservative Difference Scheme , 1997 .

[10]  Luca Bonaventura,et al.  Consistency with continuity in conservative advection schemes for free‐surface models , 2002 .

[11]  Aimé Fournier,et al.  Voronoi, Delaunay, and Block-Structured Mesh Refinement for Solution of the Shallow-Water Equations on the Sphere , 2009 .

[12]  Michael Buchhold,et al.  The Operational Global Icosahedral-Hexagonal Gridpoint Model GME: Description and High-Resolution Tests , 2002 .

[13]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[14]  Kao-San Yeh,et al.  The streamline subgrid integration method: I. Quasi-monotonic second-order transport schemes , 2007, J. Comput. Phys..

[15]  O. Friedrich,et al.  Weighted Essentially Non-Oscillatory Schemes for the Interpolation of Mean Values on Unstructured Grids , 1998 .

[16]  Peter H. Lauritzen,et al.  A flux-form version of the conservative semi-Lagrangian multi-tracer transport scheme (CSLAM) on the cubed sphere grid , 2011, J. Comput. Phys..

[17]  Yoshinobu Masuda,et al.  An Integration Scheme of the Primitive Equation Model with an Icosahedral-Hexagonal Grid System and its Application to the Shallow Water Equations , 1986 .

[18]  William C. Skamarock,et al.  Conservative Transport Schemes for Spherical Geodesic Grids: High-Order Flux Operators for ODE-Based Time Integration , 2011 .

[19]  A. Iske,et al.  On the structure of function spaces in optimal recovery of point functionals for ENO-schemes by radial basis functions , 1996 .

[20]  Alexander E. MacDonald,et al.  A Finite-Volume Icosahedral Shallow-Water Model on a Local Coordinate , 2009 .

[21]  J. Thuburn.  A PV-Based Shallow-Water Model on a Hexagonal-Icosahedral Grid , 1997 .

[22]  Hiroaki Miura,et al.  A Comparison of Grid Quality of Optimized Spherical Hexagonal–Pentagonal Geodesic Grids , 2005 .

[23]  Hiroaki Miura,et al.  A Fourth-Order-Centered Finite-Volume Scheme for Regular Hexagonal Grids , 2007 .

[24]  Hiroaki Miura,et al.  A Madden-Julian Oscillation Event Realistically Simulated by a Global Cloud-Resolving Model , 2007, Science.

[25]  R. Heikes,et al.  Numerical Integration of the Shallow-Water Equations on a Twisted Icosahedral Grid , 1995 .

[26]  Todd D. Ringler,et al.  A Multiscale Nonhydrostatic Atmospheric Model Using Centroidal Voronoi Tesselations and C-Grid Staggering , 2012 .

[27]  J. Thuburn.  Multidimensional Flux-Limited Advection Schemes , 1996 .

[28]  Gene H. Golub,et al.  Singular value decomposition and least squares solutions , 1970, Milestones in Matrix Computation.

[29]  William H. Lipscomb,et al.  An Incremental Remapping Transport Scheme on a Spherical Geodesic Grid , 2005 .

[30]  William Skamarock,et al.  Conservative Transport Schemes for Spherical Geodesic Grids: High-Order Reconstructions for Forward-in-Time Schemes , 2010 .

[31]  Ryoichi Imasu,et al.  A Three-Dimensional Icosahedral Grid Advection Scheme Preserving Monotonicity and Consistency with C , 2011 .

[32]  P. Rasch.  Conservative Shape-Preserving Two-Dimensional Transport on a Spherical Reduced Grid , 1994 .

[33]  Hirofumi Tomita,et al.  Shallow water model on a modified icosahedral geodesic grid by using spring dynamics , 2001 .

[34]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. IV. A new approach to numerical convection , 1977 .

[35]  Rashmi Mittal,et al.  On Near-Diffusion-Free Advection over Spherical Geodesic Grids , 2007 .

[36]  D. Randall,et al.  A Potential Enstrophy and Energy Conserving Numerical Scheme for Solution of the Shallow-Water Equations on a Geodesic Grid , 2002 .