Spectral interpolation in semi-implicit semi-Lagrangian methods for shallow water equations on the sphere
Fuente:
arXiv
Guardado en:
| Autores principales: | , , |
|---|---|
| Formato: | Preprint |
| Publicado: |
2026
|
| Materias: | |
| Acceso en línea: | |
| Etiquetas: |
Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
|
| _version_ | 1866911640404885504 |
|---|---|
| author | Chiwere, Michael Fortunato, Daniel Wright, Grady B. |
| author_facet | Chiwere, Michael Fortunato, Daniel Wright, Grady B. |
| contents | Semi-implicit semi-Lagrangian (SISL) methods are commonly used for the shallow water equations (SWE) because they allow for larger time steps than those permitted by the Courant-Friedrichs-Lewy (CFL) stability condition in Eulerian schemes. In these methods, the semi-Lagrangian treatment of advection is typically performed using lower-order interpolation, such as tensor-product Lagrange interpolation with cubic or quintic polynomials. However, operational SISL schemes routinely employ spectrally accurate spatial discretizations, such as spherical harmonics or the double Fourier sphere (DFS) method, for computing horizontal derivatives of the prognostic variables. This creates a mismatch in numerical accuracy, making the use of low-order interpolation less clearly justified. In this work, we present the first numerical investigation of spectrally accurate interpolation in SISL schemes for the SWE. Our approach builds upon the recently developed DFS-based SWE model, incorporating a spectral interpolation scheme that is accelerated using the nonuniform fast Fourier transform (NUFFT) to maintain the same overall computational complexity as the original model. Using several standard SWE test cases, we evaluate the accuracy, conservation, and numerical diffusion of the new model, particularly over long integration times. Compared to an equivalent SISL model with low-order interpolation, the new model achieves higher accuracy, improved mass and energy conservation, and reduced numerical diffusion, demonstrating the potential benefits of incorporating spectrally accurate interpolation into SISL schemes. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2605_01313 |
| institution | arXiv |
| publishDate | 2026 |
| record_format | arxiv |
| spellingShingle | Spectral interpolation in semi-implicit semi-Lagrangian methods for shallow water equations on the sphere Chiwere, Michael Fortunato, Daniel Wright, Grady B. Numerical Analysis 65T40, 65D05, 65M25, 42A15 Semi-implicit semi-Lagrangian (SISL) methods are commonly used for the shallow water equations (SWE) because they allow for larger time steps than those permitted by the Courant-Friedrichs-Lewy (CFL) stability condition in Eulerian schemes. In these methods, the semi-Lagrangian treatment of advection is typically performed using lower-order interpolation, such as tensor-product Lagrange interpolation with cubic or quintic polynomials. However, operational SISL schemes routinely employ spectrally accurate spatial discretizations, such as spherical harmonics or the double Fourier sphere (DFS) method, for computing horizontal derivatives of the prognostic variables. This creates a mismatch in numerical accuracy, making the use of low-order interpolation less clearly justified. In this work, we present the first numerical investigation of spectrally accurate interpolation in SISL schemes for the SWE. Our approach builds upon the recently developed DFS-based SWE model, incorporating a spectral interpolation scheme that is accelerated using the nonuniform fast Fourier transform (NUFFT) to maintain the same overall computational complexity as the original model. Using several standard SWE test cases, we evaluate the accuracy, conservation, and numerical diffusion of the new model, particularly over long integration times. Compared to an equivalent SISL model with low-order interpolation, the new model achieves higher accuracy, improved mass and energy conservation, and reduced numerical diffusion, demonstrating the potential benefits of incorporating spectrally accurate interpolation into SISL schemes. |
| title | Spectral interpolation in semi-implicit semi-Lagrangian methods for shallow water equations on the sphere |
| topic | Numerical Analysis 65T40, 65D05, 65M25, 42A15 |
| url | https://arxiv.org/abs/2605.01313 |