Editorial Type:
Article
Article Type:
Research Article

On the Spectral Convergence of the Supercompact Finite-Difference Schemes for the f-Plane Shallow-Water Equations

S. Ghader Institute of Geophysics, University of Tehran, Tehran, Iran

Search for other papers by S. Ghader in
Current site
Google Scholar
PubMed Close
,
A. R. Mohebalhojeh Institute of Geophysics, University of Tehran, Tehran, Iran

Search for other papers by A. R. Mohebalhojeh in
Current site
Google Scholar
PubMed Close
, and
V. Esfahanian Mechanical Engineering Department, University of Tehran, Tehran, Iran

Search for other papers by V. Esfahanian in
Current site
Google Scholar
PubMed Close
Print Publication:
01 Jul 2009
DOI:
https://doi.org/10.1175/2009MWR2824.1
Page(s):
2393–2406
Restricted access

Abstract

For the f-plane shallow-water equations, the convergence properties of the supercompact finite-difference method (SCFDM) are examined during the evolution of complex, nonlinear flows spawned by an unstable jet. The second-, fourth-, sixth-, and eighth-order SCFDMs are compared with a standard pseudospectral (PS) method. To control the buildup of small-scale activity and thus the potential for numerical instability, the vorticity field is damped explicitly by the application of a triharmonic hyperdiffusion operator acting on the vorticity field. The global distribution of mass between isolevels of potential vorticity, called mass error, and the representation of the balance and imbalance are used to assess numerical accuracy. In each of the quantitative measures, a clear convergence of the SCFDM to the PS method is observed. There is no saturation in accuracy up to the eighth order examined. Taking the PS solution as the reference, for the fundamental quantity of potential vorticity the rate of convergence to PS turns out to be algebraic and near-quadratic.

Corresponding author address: Dr. Sarmad Ghader, Institute of Geophysics, University of Tehran, North Karegar Ave., Tehran, Iran. Email: sghader@ut.ac.ir

Abstract

For the f-plane shallow-water equations, the convergence properties of the supercompact finite-difference method (SCFDM) are examined during the evolution of complex, nonlinear flows spawned by an unstable jet. The second-, fourth-, sixth-, and eighth-order SCFDMs are compared with a standard pseudospectral (PS) method. To control the buildup of small-scale activity and thus the potential for numerical instability, the vorticity field is damped explicitly by the application of a triharmonic hyperdiffusion operator acting on the vorticity field. The global distribution of mass between isolevels of potential vorticity, called mass error, and the representation of the balance and imbalance are used to assess numerical accuracy. In each of the quantitative measures, a clear convergence of the SCFDM to the PS method is observed. There is no saturation in accuracy up to the eighth order examined. Taking the PS solution as the reference, for the fundamental quantity of potential vorticity the rate of convergence to PS turns out to be algebraic and near-quadratic.

Corresponding author address: Dr. Sarmad Ghader, Institute of Geophysics, University of Tehran, North Karegar Ave., Tehran, Iran. Email: sghader@ut.ac.ir

Save
Cover Monthly Weather Review
Monthly Weather Review

  • Asselin, R., 1972: Frequency filter for time integrations. Mon. Wea. Rev., 100 , 487–490.

  • Chung, T. J., 2002: Computational Fluid Dynamics. Cambridge University Press, 1012 pp.

  • Dritschel, D. G., L. M. Polvani, and A. R. Mohebalhojeh, 1999: Contour-advective semi-Lagrangian algorithm for the shallow-water equations. Mon. Wea. Rev., 127 , 1151–1165.

    • Search Google Scholar
    • Export Citation

  • Durran, D. R., 1999: Numerical Methods for Wave Equations in Geophysical Fluid Dynamics. Springer, 465 pp.

  • Esfahanian, V., S. Ghader, and K. Ashrafi, 2004: Accuracy analysis of super compact scheme in nonuniform grid with application to parabolized stability equations. Int. J. Numer. Methods Fluids, 46 , 485–505.

    • Search Google Scholar
    • Export Citation

  • Esfahanian, V., S. Ghader, and A. R. Mohebalhojeh, 2005: On the use of the super compact scheme for spatial differencing in numerical models of the atmosphere. Quart. J. Roy. Meteor. Soc., 131 , 2109–2129.

    • Search Google Scholar
    • Export Citation

  • Ford, R., M. E. McIntyre, and W. A. Norton, 2000: Balance and the slow quasimanifold: Some explicit results. J. Atmos. Sci., 57 , 1236–1254.

    • Search Google Scholar
    • Export Citation

  • Fu, D., and Y. Ma, 1995: High resolution schemes. Computational Fluid Dynamic Review, M. Hafez and K. Oshima, Eds., John Wiley and Sons, 234–250.

    • Search Google Scholar
    • Export Citation

  • Galewsky, J., R. K. Scott, and L. M. Polvani, 2004: An initial-value problem for testing numerical models of the global shallow-water equations. Tellus, 56A , 429–440.

    • Search Google Scholar
    • Export Citation

  • Ghader, S., 2000: Numerical simulation of the stability of the laminar boundary layer flow over a flat plate by PSE equations using super compact method (in Persian). M.S. thesis, Dept. of Mechanical Engineering, University of Tehran, 91 pp.

  • Ghader, S., 2005: Super compact spatial differencing for the numerical simulation of fluid dynamics problems with application to the shallow water equations (in Persian). Ph.D. thesis, University of Tehran, 197 pp.

  • Holdaway, D., J. Thuburn, and N. Wood, 2007: On the relation between order of accuracy, convergence rate and spectral slope for linear numerical methods applied to multiscale problems. Int. J. Numer. Methods Fluids, 56 , 1297–1303.

    • Search Google Scholar
    • Export Citation

  • Hoskins, B. J., M. E. McIntyre, and A. W. Robertson, 1985: On the use and significance of isentropic potential-vorticity maps. Quart. J. Roy. Meteor. Soc., 111 , 877–946.

    • Search Google Scholar
    • Export Citation

  • Iserles, A., 1996: A First Course in the Numerical Analysis of Differential Equations. Cambridge University Press, 378 pp.

  • Lele, S. K., 1992: Compact finite difference schemes with spectral-like resolution. J. Comput. Phys., 103 , 16–42.

  • Ma, Y., and D. Fu, 1996: Super compact finite difference method (SCFDM) with arbitrary high accuracy. Comput. Fluid Dyn. J., 5 , 259–276.

    • Search Google Scholar
    • Export Citation

  • Merilees, P. E., and S. A. Orzag, 1979: The pseudo-spectral method. Numerical Methods Used in Atmospheric Models, GARP Publication Series, Vol. II, No. 17, World Meteorological Organization, 276–299.

    • Search Google Scholar
    • Export Citation

  • Mohebalhojeh, A. R., and D. G. Dritschel, 2000: On the representation of gravity waves in numerical models of the shallow water equations. Quart. J. Roy. Meteor. Soc., 126 , 669–688.

    • Search Google Scholar
    • Export Citation

  • Mohebalhojeh, A. R., and D. G. Dritschel, 2001: Hierarchies of balance conditions for the f-plane shallow water equations. J. Atmos. Sci., 58 , 2411–2426.

    • Search Google Scholar
    • Export Citation

  • Mohebalhojeh, A. R., and D. G. Dritschel, 2004: Contour-advective semi-Lagrangian algorithms for many-layer primitive equation models. Quart. J. Roy. Meteor. Soc., 130 , 347–364.

    • Search Google Scholar
    • Export Citation

  • Mohebalhojeh, A. R., and D. G. Dritschel, 2007: Assessing the numerical accuracy of complex spherical shallow water flows. Mon. Wea. Rev., 135 , 3876–3894.

    • Search Google Scholar
    • Export Citation

  • Naulin, V., and A. H. Nielsen, 2003: Accuracy of spectral and finite difference schemes in 2D advection problems. SIAM J. Sci. Comput., 25 , 104–126.

    • Search Google Scholar
    • Export Citation

  • Polvani, L. M., R. K. Scott, and S. J. Thomas, 2004: Numerically converged solutions of the global primitive equations for testing the dynamical core of atmospheric GCMs. Mon. Wea. Rev., 132 , 2539–2552.

    • Search Google Scholar
    • Export Citation

  • Ritchie, H., 1988: Application of the semi-Lagrangian method to a spectral model of the shallow water equations. Mon. Wea. Rev., 116 , 1587–1598.

    • Search Google Scholar
    • Export Citation

  • Robert, A., 1966: The integration of a low order spectral form of the primitive meteorological equations. J. Meteor. Soc. Japan, 44 , 237–245.

    • Search Google Scholar
    • Export Citation

  • Salmon, R., 1998: Lectures on Geophysical Fluid Dynamics. Oxford University Press, 378 pp.

  • Smolarkiewicz, P. K., and L. G. Margolin, 1994: Variational solver for elliptic problems in atmospheric flows. Appl. Math. Comp. Sci., 4 , 527–551.

    • Search Google Scholar
    • Export Citation

  • Smolarkiewicz, P. K., and L. G. Margolin, 2000: Variational methods for elliptic problems in fluid models. Proc. Workshop on Developments in Numerical Methods for Very High Resolution Global Models, Reading, United Kingdom, ECMWF, 49–68.

    • Search Google Scholar
    • Export Citation
All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 1653 1494 120
PDF Downloads 118 30 1