Editorial Type:
Article
Article Type:
Research Article

Implicit–Explicit Multistep Methods for Fast-Wave–Slow-Wave Problems

Dale R. Durran University of Washington, Seattle, Washington

Search for other papers by Dale R. Durran in
Current site
Google Scholar
PubMed Close
and
Peter N. Blossey University of Washington, Seattle, Washington

Search for other papers by Peter N. Blossey in
Current site
Google Scholar
PubMed Close
Print Publication:
01 Apr 2012
DOI:
https://doi.org/10.1175/MWR-D-11-00088.1
Page(s):
1307–1325
Corrigendum to this article
Restricted access

Abstract

Implicit–explicit (IMEX) linear multistep methods are examined with respect to their suitability for the integration of fast-wave–slow-wave problems in which the fast wave has relatively low amplitude and need not be accurately simulated. The widely used combination of trapezoidal implicit and leapfrog explicit differencing is compared to schemes based on Adams methods or on backward differencing. Two new families of methods are proposed that have good stability properties in fast-wave–slow-wave problems: one family is based on Adams methods and the other on backward schemes. Here the focus is primarily on four specific schemes drawn from these two families: a pair of Adams methods and a pair of backward methods that are either (i) optimized for third-order accuracy in the explicit component of the full IMEX scheme, or (ii) employ particularly good schemes for the implicit component. These new schemes are superior, in many respects, to the linear multistep IMEX schemes currently in use.

The behavior of these schemes is compared theoretically in the context of the simple oscillation equation and also for the linearized equations governing stratified compressible flow. Several schemes are also tested in fully nonlinear simulations of gravity waves generated by a localized source in a shear flow.

Corresponding author address: Dale R. Durran, Dept. of Atmospheric Sciences, University of Washington, Box 351640, Seattle, WA 98195-1640. E-mail: drdee@uw.edu

Abstract

Implicit–explicit (IMEX) linear multistep methods are examined with respect to their suitability for the integration of fast-wave–slow-wave problems in which the fast wave has relatively low amplitude and need not be accurately simulated. The widely used combination of trapezoidal implicit and leapfrog explicit differencing is compared to schemes based on Adams methods or on backward differencing. Two new families of methods are proposed that have good stability properties in fast-wave–slow-wave problems: one family is based on Adams methods and the other on backward schemes. Here the focus is primarily on four specific schemes drawn from these two families: a pair of Adams methods and a pair of backward methods that are either (i) optimized for third-order accuracy in the explicit component of the full IMEX scheme, or (ii) employ particularly good schemes for the implicit component. These new schemes are superior, in many respects, to the linear multistep IMEX schemes currently in use.

The behavior of these schemes is compared theoretically in the context of the simple oscillation equation and also for the linearized equations governing stratified compressible flow. Several schemes are also tested in fully nonlinear simulations of gravity waves generated by a localized source in a shear flow.

Corresponding author address: Dale R. Durran, Dept. of Atmospheric Sciences, University of Washington, Box 351640, Seattle, WA 98195-1640. E-mail: drdee@uw.edu
Save
Cover Monthly Weather Review
Monthly Weather Review

  • Ascher, U. M., S. J. Ruuth, and B. T. R. Wetton, 1995: Implicit-explicit methods for time-dependent partial differential equations. SIAM J. Numer. Anal., 32, 797–823.

    • Search Google Scholar
    • Export Citation

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

  • Bannon, P. R., 1996: On the anelastic approximation for a compressible atmosphere. J. Atmos. Sci., 53, 3618–3628.

  • Benoit, R., M. Desgagné, P. Pellerin, S. Pellerin, Y. Chartier, and S. Desjardins, 1997: The Canadian MC2: A semi-Lagrangian, semi-implicit wideband atmospheric model suited for finescale process studies and simulation. Mon. Wea. Rev., 125, 2382–2415.

    • Search Google Scholar
    • Export Citation

  • Dahlquist, G. G., 1963: A special stability problem for linear multistep methods. BIT Numer. Math., 3, 27–43.

  • Durran, D. R., 1989: Improving the anelastic approximation. J. Atmos. Sci., 46, 1453–1461.

  • Durran, D. R., 1991: The third-order Adams–Bashforth method: An attractive alternative to leapfrog time differencing. Mon. Wea. Rev., 119, 702–720.

    • Search Google Scholar
    • Export Citation

  • Durran, D. R., 2008: A physically motivated approach for filtering acoustic waves from the equations governing compressible stratified flow. J. Fluid Mech., 601, 365–379.

    • Search Google Scholar
    • Export Citation

  • Durran, D. R., 2010: Numerical Methods for Fluid Dynamics: With Applications in Geophysics. Springer, 516 pp.

  • Evans, K. J., M. A. Taylor, and J. B. Drake, 2010: Accuracy analysis of a spectral element atmospheric model using a fully implicit solution framework. Mon. Wea. Rev., 138, 3333–3341.

    • Search Google Scholar
    • Export Citation

  • Fornberg, B., and T. A. Driscoll, 1999: A fast spectral algorithm for nonlinear wave equations with linear dispersion. J. Comput. Phys., 155, 456–467.

    • Search Google Scholar
    • Export Citation

  • Frank, J., W. Hundsdorfer, and J. G. Verwer, 1997: On the stability of implicit-explicit linear multistep methods. Appl. Numer. Math., 25, 193–205.

    • Search Google Scholar
    • Export Citation

  • Giraldo, F. X., 2005: Semi-implicit time-integrators for a scalable spectral element atmospheric model. Quart. J. Roy. Meteor. Soc., 131, 2431–2454, doi:10.1256/qj.03.218.

    • Search Google Scholar
    • Export Citation

  • Giraldo, F. X., M. Restelli, and M. Läuter, 2010: Semi-implicit formulations of the Navier–Stokes equations: Application to nonhydrostatic atmospheric modeling. SIAM J. Sci. Comput., 32, 3394–3425.

    • Search Google Scholar
    • Export Citation

  • Karniadakis, G. E., M. Israeli, and S. A. Orszag, 1991: High-order splitting methods for the incompressible Navier-Stokes equations. J. Comput. Phys., 97, 414–443.

    • Search Google Scholar
    • Export Citation

  • Kennedy, C. A., and M. H. Carpenter, 2003: Additive Runge-Kutta schemes for convection-diffusion-reaction equations. Appl. Numer. Math., 44, 139–181, doi:10.1016/S0168-9274(02)00138-1.

    • Search Google Scholar
    • Export Citation

  • Klemp, J. B., and R. Wilhelmson, 1978: The simulation of three-dimensional convective storm dynamics. J. Atmos. Sci., 35, 1070–1096.

  • Kwizak, M., and A. J. Robert, 1971: A semi-implicit scheme for grid point atmospheric models of the primitive equation. Mon. Wea. Rev., 99, 32–36.

    • Search Google Scholar
    • Export Citation

  • Lipps, F., and R. Hemler, 1982: A scale analysis of deep moist convection and some related numerical calculations. J. Atmos. Sci., 39, 2192–2210.

    • Search Google Scholar
    • Export Citation

  • Nevanlinna, O., and W. Liniger, 1978: Contractive methods for stiff differential equations Part I. BIT Numer. Math., 18, 457–474.

  • Ogura, Y., and N. A. Phillips, 1962: Scale analysis for deep and shallow convection in the atmosphere. J. Atmos. Sci., 19, 173–179.

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

    • Search Google Scholar
    • Export Citation

  • Tapp, M. C., and P. W. White, 1976: A non-hydrostatic mesoscale model. Quart. J. Roy. Meteor. Soc., 102, 277–296.

  • Tatsumi, Y., 1983: An economical explicit time integration scheme for a primitive model. J. Meteor. Soc. Japan, 61, 269–287.

  • Ullrich, P., and C. Jablonowski, 2012: Operator-split Runge–Kutta–Rosenbrock methods for nonhydrostatic atmospheric models. Mon. Wea. Rev., 140, 1257–1284.

    • Search Google Scholar
    • Export Citation

  • Varah, J. M., 1980: Stability restrictions on second order, three level finite difference schemes for parabolic equations. SIAM J. Numer. Anal., 17, 300–309.

    • Search Google Scholar
    • Export Citation

  • Wicker, L. J., and W. C. Skamarock, 2002: Time-splitting methods for elastic models using forward time schemes. Mon. Wea. Rev., 130, 2088–2097.

    • Search Google Scholar
    • Export Citation

  • Williams, P. D., 2009: A proposed modification to the Robert–Asselin time filter. Mon. Wea. Rev., 137, 2538–2546.

  • Williams, P. D., 2011: The RAW filter: An improvement to the Robert–Asselin filter in semi-implicit integrations. Mon. Wea. Rev., 139, 1996–2007.

    • Search Google Scholar
    • Export Citation
All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 684 259 68
PDF Downloads 554 215 24