Improved double Fourier series on a sphere and its application to a semi-implicit semi-Lagrangian shallow-water model
<p>One way to reduce the computational cost of a spectral model using spherical harmonics (SH) is to use double Fourier series (DFS) instead of SH. The transform method using SH usually requires <span class="inline-formula">O(<i>N</i><sup>3</sup>)</sp...
| Main Author: | |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
Copernicus Publications
2022-03-01
|
| Series: | Geoscientific Model Development |
| Online Access: | https://gmd.copernicus.org/articles/15/2561/2022/gmd-15-2561-2022.pdf |
| _version_ | 1811273038711226368 |
|---|---|
| author | H. Yoshimura |
| author_facet | H. Yoshimura |
| author_sort | H. Yoshimura |
| collection | DOAJ |
| description | <p>One way to reduce the computational cost of a spectral
model using spherical harmonics (SH) is to use double Fourier series (DFS)
instead of SH. The transform method using SH usually requires
<span class="inline-formula">O(<i>N</i><sup>3</sup>)</span> operations, where <span class="inline-formula"><i>N</i></span> is the truncation
wavenumber, and the computational cost significantly increases at high
resolution. On the other hand, the method using DFS requires only
<span class="inline-formula">O(<i>N</i><sup>2</sup>log <i>N</i>)</span> operations. This paper proposes a new
DFS method that improves the numerical stability of the model compared with
the conventional DFS methods by adopting the following two improvements: a
new expansion method that employs the least-squares method (or the Galerkin
method) to calculate the expansion coefficients in order to minimize the
error caused by wavenumber truncation, and new basis functions that satisfy
the continuity of both scalar and vector variables at the poles. Partial
differential equations such as the Poisson equation and the Helmholtz
equation are solved by using the Galerkin method. In the semi-implicit
semi-Lagrangian shallow-water model using the new DFS method, the Williamson
test cases and the Galewsky test case give stable results without the
appearance of high-wavenumber noise near the poles, even without horizontal
diffusion and without a zonal Fourier filter. In the Eulerian advection
model using the new DFS method, the Williamson test case 1, which simulates
a cosine bell advection, also gives stable results without horizontal
diffusion but with a zonal Fourier filter. The shallow-water model using the
new DFS method is faster than that using SH, especially at high resolutions
and gives almost the same results, except that very small oscillations near
the truncation wavenumber in the kinetic energy spectrum appear only in the
shallow-water model using SH.</p> |
| first_indexed | 2024-04-12T22:52:32Z |
| format | Article |
| id | doaj.art-0a00a089766449f194c9d498d63f955a |
| institution | Directory Open Access Journal |
| issn | 1991-959X 1991-9603 |
| language | English |
| last_indexed | 2024-04-12T22:52:32Z |
| publishDate | 2022-03-01 |
| publisher | Copernicus Publications |
| record_format | Article |
| series | Geoscientific Model Development |
| spelling | doaj.art-0a00a089766449f194c9d498d63f955a2022-12-22T03:13:19ZengCopernicus PublicationsGeoscientific Model Development1991-959X1991-96032022-03-01152561259710.5194/gmd-15-2561-2022Improved double Fourier series on a sphere and its application to a semi-implicit semi-Lagrangian shallow-water modelH. Yoshimura<p>One way to reduce the computational cost of a spectral model using spherical harmonics (SH) is to use double Fourier series (DFS) instead of SH. The transform method using SH usually requires <span class="inline-formula">O(<i>N</i><sup>3</sup>)</span> operations, where <span class="inline-formula"><i>N</i></span> is the truncation wavenumber, and the computational cost significantly increases at high resolution. On the other hand, the method using DFS requires only <span class="inline-formula">O(<i>N</i><sup>2</sup>log <i>N</i>)</span> operations. This paper proposes a new DFS method that improves the numerical stability of the model compared with the conventional DFS methods by adopting the following two improvements: a new expansion method that employs the least-squares method (or the Galerkin method) to calculate the expansion coefficients in order to minimize the error caused by wavenumber truncation, and new basis functions that satisfy the continuity of both scalar and vector variables at the poles. Partial differential equations such as the Poisson equation and the Helmholtz equation are solved by using the Galerkin method. In the semi-implicit semi-Lagrangian shallow-water model using the new DFS method, the Williamson test cases and the Galewsky test case give stable results without the appearance of high-wavenumber noise near the poles, even without horizontal diffusion and without a zonal Fourier filter. In the Eulerian advection model using the new DFS method, the Williamson test case 1, which simulates a cosine bell advection, also gives stable results without horizontal diffusion but with a zonal Fourier filter. The shallow-water model using the new DFS method is faster than that using SH, especially at high resolutions and gives almost the same results, except that very small oscillations near the truncation wavenumber in the kinetic energy spectrum appear only in the shallow-water model using SH.</p>https://gmd.copernicus.org/articles/15/2561/2022/gmd-15-2561-2022.pdf |
| spellingShingle | H. Yoshimura Improved double Fourier series on a sphere and its application to a semi-implicit semi-Lagrangian shallow-water model Geoscientific Model Development |
| title | Improved double Fourier series on a sphere and its application to a semi-implicit semi-Lagrangian shallow-water model |
| title_full | Improved double Fourier series on a sphere and its application to a semi-implicit semi-Lagrangian shallow-water model |
| title_fullStr | Improved double Fourier series on a sphere and its application to a semi-implicit semi-Lagrangian shallow-water model |
| title_full_unstemmed | Improved double Fourier series on a sphere and its application to a semi-implicit semi-Lagrangian shallow-water model |
| title_short | Improved double Fourier series on a sphere and its application to a semi-implicit semi-Lagrangian shallow-water model |
| title_sort | improved double fourier series on a sphere and its application to a semi implicit semi lagrangian shallow water model |
| url | https://gmd.copernicus.org/articles/15/2561/2022/gmd-15-2561-2022.pdf |
| work_keys_str_mv | AT hyoshimura improveddoublefourierseriesonasphereanditsapplicationtoasemiimplicitsemilagrangianshallowwatermodel |