Abstract

A time-stepping scheme is proposed. It is based on the leapfrog method and a fourth-order time filter. The scheme requires only one evaluation per time step and uses an implicit filter, but the effort needed to implement it in an explicit manner is trivial. Comparative tests demonstrate that the proposed scheme produces numerical approximations that are more stable and highly accurate compared to the standard Robert-Asselin (RA) and the Robert-Asselin-Williams (RAW) filtered leapfrog scheme, even though both methods use filter coefficients that are tuned such that the 2Δt modes are damped at the same rate. Formal stability analysis demonstrates that the proposed method generates amplitude errors of O[(Δt)4], implying third-order accuracy. This contrasts with the O[(Δt)2] errors produced by the standard RA and RAW filtered leapfrog. A second scheme that produces amplitude errors of O[(Δt)6] is also presented. The proposed scheme is found to do well at controlling numerical instabilities arising in the diffusion equation and in nonlinear computations using Lorenz's system and the global shallow-water spectral model. In addition to noticeably improving the resolution of the physical modes, the proposed method is simple to implement and has a wider region of stability compared to the existing time-filtered leapfrog schemes. This makes the proposed method a potential alternative for use in atmospheric, oceanic, and climate modeling.

Original languageEnglish (US)
Pages (from-to)2545-2560
Number of pages16
JournalMonthly Weather Review
Volume142
Issue number7
DOIs
StatePublished - 2014

Fingerprint

numerical method
filter
stability analysis
climate modeling
shallow water
method

Keywords

  • Climate models
  • Mesoscale models
  • Model comparison
  • Numerical weather prediction/forecasting

ASJC Scopus subject areas

  • Atmospheric Science

Cite this

@article{a5fd9a3b540d4d6a85aebf3e47551174,
title = "A numerical method based on leapfrog and a fourth-order implicit time filter",
abstract = "A time-stepping scheme is proposed. It is based on the leapfrog method and a fourth-order time filter. The scheme requires only one evaluation per time step and uses an implicit filter, but the effort needed to implement it in an explicit manner is trivial. Comparative tests demonstrate that the proposed scheme produces numerical approximations that are more stable and highly accurate compared to the standard Robert-Asselin (RA) and the Robert-Asselin-Williams (RAW) filtered leapfrog scheme, even though both methods use filter coefficients that are tuned such that the 2Δt modes are damped at the same rate. Formal stability analysis demonstrates that the proposed method generates amplitude errors of O[(Δt)4], implying third-order accuracy. This contrasts with the O[(Δt)2] errors produced by the standard RA and RAW filtered leapfrog. A second scheme that produces amplitude errors of O[(Δt)6] is also presented. The proposed scheme is found to do well at controlling numerical instabilities arising in the diffusion equation and in nonlinear computations using Lorenz's system and the global shallow-water spectral model. In addition to noticeably improving the resolution of the physical modes, the proposed method is simple to implement and has a wider region of stability compared to the existing time-filtered leapfrog schemes. This makes the proposed method a potential alternative for use in atmospheric, oceanic, and climate modeling.",
keywords = "Climate models, Mesoscale models, Model comparison, Numerical weather prediction/forecasting",
author = "Mohamed Moustaoui and Alex Mahalov and Eric Kostelich",
year = "2014",
doi = "10.1175/MWR-D-13-00073.1",
language = "English (US)",
volume = "142",
pages = "2545--2560",
journal = "Monthly Weather Review",
issn = "0027-0644",
publisher = "American Meteorological Society",
number = "7",

}

TY - JOUR

T1 - A numerical method based on leapfrog and a fourth-order implicit time filter

AU - Moustaoui, Mohamed

AU - Mahalov, Alex

AU - Kostelich, Eric

PY - 2014

Y1 - 2014

N2 - A time-stepping scheme is proposed. It is based on the leapfrog method and a fourth-order time filter. The scheme requires only one evaluation per time step and uses an implicit filter, but the effort needed to implement it in an explicit manner is trivial. Comparative tests demonstrate that the proposed scheme produces numerical approximations that are more stable and highly accurate compared to the standard Robert-Asselin (RA) and the Robert-Asselin-Williams (RAW) filtered leapfrog scheme, even though both methods use filter coefficients that are tuned such that the 2Δt modes are damped at the same rate. Formal stability analysis demonstrates that the proposed method generates amplitude errors of O[(Δt)4], implying third-order accuracy. This contrasts with the O[(Δt)2] errors produced by the standard RA and RAW filtered leapfrog. A second scheme that produces amplitude errors of O[(Δt)6] is also presented. The proposed scheme is found to do well at controlling numerical instabilities arising in the diffusion equation and in nonlinear computations using Lorenz's system and the global shallow-water spectral model. In addition to noticeably improving the resolution of the physical modes, the proposed method is simple to implement and has a wider region of stability compared to the existing time-filtered leapfrog schemes. This makes the proposed method a potential alternative for use in atmospheric, oceanic, and climate modeling.

AB - A time-stepping scheme is proposed. It is based on the leapfrog method and a fourth-order time filter. The scheme requires only one evaluation per time step and uses an implicit filter, but the effort needed to implement it in an explicit manner is trivial. Comparative tests demonstrate that the proposed scheme produces numerical approximations that are more stable and highly accurate compared to the standard Robert-Asselin (RA) and the Robert-Asselin-Williams (RAW) filtered leapfrog scheme, even though both methods use filter coefficients that are tuned such that the 2Δt modes are damped at the same rate. Formal stability analysis demonstrates that the proposed method generates amplitude errors of O[(Δt)4], implying third-order accuracy. This contrasts with the O[(Δt)2] errors produced by the standard RA and RAW filtered leapfrog. A second scheme that produces amplitude errors of O[(Δt)6] is also presented. The proposed scheme is found to do well at controlling numerical instabilities arising in the diffusion equation and in nonlinear computations using Lorenz's system and the global shallow-water spectral model. In addition to noticeably improving the resolution of the physical modes, the proposed method is simple to implement and has a wider region of stability compared to the existing time-filtered leapfrog schemes. This makes the proposed method a potential alternative for use in atmospheric, oceanic, and climate modeling.

KW - Climate models

KW - Mesoscale models

KW - Model comparison

KW - Numerical weather prediction/forecasting

UR - http://www.scopus.com/inward/record.url?scp=84905653143&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84905653143&partnerID=8YFLogxK

U2 - 10.1175/MWR-D-13-00073.1

DO - 10.1175/MWR-D-13-00073.1

M3 - Article

AN - SCOPUS:84905653143

VL - 142

SP - 2545

EP - 2560

JO - Monthly Weather Review

JF - Monthly Weather Review

SN - 0027-0644

IS - 7

ER -