• Colella, P., , and P. R. Woodward, 1984: The Piecewise Parabolic Method (PPM) for gas-dynamical simulations. J. Comput. Phys., 54, 174201.

    • Search Google Scholar
    • Export Citation
  • Friedrich, O., 1998: Weighted essentially non-oscillatory schemes for the interpolation of mean values on unstructured grids. J. Comput. Phys., 144, 194212.

    • Search Google Scholar
    • Export Citation
  • Golub, G. H., , and C. Reinsch, 1970: Singular value decomposition and least squares solutions. Numer. Math., 14, 403420.

  • Gross, E. S., , L. Bonaventura, , and G. Rosatti, 2002: Consistency with continuity in conservative advection schemes for free-surface models. Int. J. Numer. Methods Fluids, 38, 307327.

    • Search Google Scholar
    • Export Citation
  • Harris, L. M., , P. H. Lauritzen, , and R. Mittal, 2011: A flux-form version of the conservative semi-Lagrangian multi-tracer transport scheme (CSLAM) on the cubed sphere grid. J. Comput. Phys., 230, 12151237.

    • Search Google Scholar
    • Export Citation
  • Heikes, R., , and D. A. Randall, 1995: Numerical integration of the shallow-water equations on a twisted icosahedral grid. Part I: Basic design and results of tests. Mon. Wea. Rev., 123, 18621880.

    • Search Google Scholar
    • Export Citation
  • Iske, A., , and T. Sonar, 1996: On the structure of function spaces in optimal recovery of point functionals for ENO-schemes by radial basis functions. Numer. Math., 74, 177202.

    • Search Google Scholar
    • Export Citation
  • Lee, J.-L., , and A. E. MacDonald, 2009: A finite-volume icosahedral shallow-water model on a local coordinate. Mon. Wea. Rev., 137, 14221437.

    • Search Google Scholar
    • Export Citation
  • Leonard, B. P., , M. K. MacVean, , and A. P. Lock, 1993: Positivity-preserving schemes for multidimensional advection. NASA Tech. Memo. 106055/ICOMP-93-05, Institute for Computational Mechanics in Propulsion, Lewis Research Center, Cleveland, OH, 62 pp.

  • Lipscomb, W. H., , and T. D. Ringler, 2005: An incremental remapping transport scheme on a spherical geodesic grid. Mon. Wea. Rev., 133, 23352350.

    • Search Google Scholar
    • Export Citation
  • Majewski, D., and Coauthors, 2002: The operational global icosahedral-hexagonal gridpoint model GME: Description and high-resolution tests. Mon. Wea. Rev., 130, 319338.

    • Search Google Scholar
    • Export Citation
  • Masuda, Y., , and H. Ohnishi, 1986: An integration scheme of the primitive equations model with an icosahedral–hexagonal grid system and its application to the shallow water equations. Short- and Medium-Range Numerical Weather Prediction, T. Matsuno, Ed., Japan Meteorological Society, 317–326.

  • Mittal, R., , H. C. Upadhyaya, , and O. P. Sharma, 2007: On near-diffusion-free advection over spherical geodesic grids. Mon. Wea. Rev., 135, 42144225.

    • Search Google Scholar
    • Export Citation
  • Miura, H., 2007a: A fourth-order-centered finite-volume scheme for regular hexagonal grids. Mon. Wea. Rev., 135, 40304037.

  • Miura, H., 2007b: An upwind-biased conservative advection scheme for spherical hexagonal–pentagonal grids. Mon. Wea. Rev., 135, 40384044.

    • Search Google Scholar
    • Export Citation
  • Miura, H., , and M. Kimoto, 2005: A comparison of grid quality of optimized spherical hexagonal–pentagonal geodesic grids. Mon. Wea. Rev., 133, 28172833.

    • Search Google Scholar
    • Export Citation
  • Miura, H., , M. Satoh, , T. Nasuno, , A. T. Noda, , and K. Oouchi, 2007: A Madden-Julian oscillation event realistically simulated by a global cloud-resolving model. Science, 318, 17631765.

    • Search Google Scholar
    • Export Citation
  • Nair, R. D., , and P. H. Lauritzen, 2010: A class of deformational flow test cases for linear transport problems on the sphere. J. Comput. Phys., 229, 88688887.

    • Search Google Scholar
    • Export Citation
  • Niwa, Y., , H. Tomita, , M. Satoh, , and R. Imasu, 2011: A three-dimensional icosahedral grid advection scheme preserving monotonicity and consistency with continuity for atmospheric tracer transport. J. Meteor. Soc. Japan, 89, 255268.

    • Search Google Scholar
    • Export Citation
  • Rasch, P. J., 1994: Conservative shape-preserving two-dimensional transport on a spherical reduced grid. Mon. Wea. Rev., 122, 13371350.

    • Search Google Scholar
    • Export Citation
  • Ringler, T. D., , and D. A. Randall, 2002: A potential enstrophy and energy conserving numerical scheme for solution of the shallow-water equations on a geodesic grid. Mon. Wea. Rev., 130, 13971410.

    • Search Google Scholar
    • Export Citation
  • Ringler, T. D., , R. P. Heikes, , and D. A. Randall, 2000: Modeling the atmospheric general circulation using a spherical geodesic grid: A new class of dynamical cores. Mon. Wea. Rev., 128, 24712490.

    • Search Google Scholar
    • Export Citation
  • Satoh, M., , T. Matsuno, , H. Tomita, , H. Miura, , T. Nasuno, , and S. Iga, 2008: Nonhydrostatic icosahedral atmospheric model (NICAM) for global cloud resolving simulations. J. Comput. Phys., 227, 34843514.

    • Search Google Scholar
    • Export Citation
  • Skamarock, W. C., , and M. Menchaca, 2010: Conservative transport schemes for spherical geodesic grids: High-order reconstructions for forward-in-time schemes. Mon. Wea. Rev., 138, 44974508.

    • Search Google Scholar
    • Export Citation
  • Skamarock, W. C., , and A. Gassmann, 2011: Conservative transport schemes for spherical geodesic grids: High-order flux operators for ODE-based time integration. Mon. Wea. Rev., 139, 29622975.

    • Search Google Scholar
    • Export Citation
  • Skamarock, W. C., , J. B. Klemp, , M. G. Duda, , L. Fowler, , S.-H. Park, , and T. D. Ringler, 2012: A multiscale nonhydrostatic atmospheric model using centroidal Voronoi tesselations and C-grid staggering. Mon. Wea. Rev., 140, 30903105.

    • Search Google Scholar
    • Export Citation
  • Stuhne, G. R., , and W. R. Peltier, 1996: Vortex erosion and amalgamation in a new model of large scale flow on the sphere. J. Comput. Phys., 128, 5881.

    • Search Google Scholar
    • Export Citation
  • Thuburn, J., 1995: Dissipation and cascades to small scales in numerical models using a shape-preserving advection scheme. Mon. Wea. Rev., 123, 18881903.

    • Search Google Scholar
    • Export Citation
  • Thuburn, J., 1996: Multidimensional flux-limited advection schemes. J. Comput. Phys., 123, 7483.

  • Thuburn, J., 1997: A PV-based shallow-water model on a hexagonal–icosahedral grid. Mon. Wea. Rev., 125, 23282347.

  • Tomita, H., , M. Tsugawa, , M. Satoh, , and K. Goto, 2001: Shallow water model on a modified icosahedral geodesic grid by using spring dynamics. J. Comput. Phys., 174, 579613.

    • Search Google Scholar
    • Export Citation
  • van Leer, B., 1977: Towards the ultimate conservative difference scheme. IV. A new approach to numerical convection. J. Comput. Phys., 23, 276299.

    • Search Google Scholar
    • Export Citation
  • Weller, H., , H. G. Weller, , and A. Fournier, 2009: Voronoi, Delaunay, and block-structured mesh refinement for solution of the shallow-water equations on the sphere. Mon. Wea. Rev., 137, 42084224.

    • Search Google Scholar
    • Export Citation
  • Williamson, D. L., , J. B. Drake, , J. J. Hack, , R. Jakob, , and P. N. Swarztrauber, 1992: A standard test set for numerical approximations to the shallow water equations in spherical geometry. J. Comput. Phys., 102, 211224.

    • Search Google Scholar
    • Export Citation
  • Yeh, K.-S., 2007: The streamline subgrid integration method: I. Quasi-monotonic second-order transport schemes. J. Comput. Phys., 225, 16321652.

    • Search Google Scholar
    • Export Citation
All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 24 24 8
PDF Downloads 21 21 5

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

View More View Less
  • 1 Atmosphere and Ocean Research Institute, The University of Tokyo, Chiba, Japan
  • | 2 Mesoscale and Microscale Meteorology Division, National Center for Atmospheric Research, Boulder, Colorado
© Get Permissions
Restricted access

Abstract

Several 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 upwind-biased quadratic approximation (UQA-2)] is more accurate than the preceding second-order reconstruction scheme [the first upwind-biased quadratic approximation (UQA-1)] in most of the tests in this work. Solutions of UQA-2 are sharper than those of UQA-1, although with slightly larger phase errors. The computational cost of UQA-2 is comparable to UQA-1.

Corresponding author address: Hiroaki Miura, Department of Earth and Planetary Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan. E-mail: h_miura@eps.s.u-tokyo.ac.jp

Abstract

Several 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 upwind-biased quadratic approximation (UQA-2)] is more accurate than the preceding second-order reconstruction scheme [the first upwind-biased quadratic approximation (UQA-1)] in most of the tests in this work. Solutions of UQA-2 are sharper than those of UQA-1, although with slightly larger phase errors. The computational cost of UQA-2 is comparable to UQA-1.

Corresponding author address: Hiroaki Miura, Department of Earth and Planetary Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan. E-mail: h_miura@eps.s.u-tokyo.ac.jp
Save