Azimuthal averaging–reconstruction filtering techniques for finite-difference general circulation models in spherical geometry
-
Published:2021-02-10
Issue:2
Volume:14
Page:859-873
-
ISSN:1991-9603
-
Container-title:Geoscientific Model Development
-
language:en
-
Short-container-title:Geosci. Model Dev.
Author:
Dang Tong,Zhang Binzheng,Lei Jiuhou,Wang Wenbin,Burns Alan,Liu Han-li,Pham Kevin,Sorathia Kareem A.
Abstract
Abstract. When solving hydrodynamic equations in spherical or cylindrical geometry using explicit finite-difference schemes, a major difficulty is that the time step is greatly restricted by the clustering of azimuthal cells near the pole due to the Courant–Friedrichs–Lewy condition. This paper adapts the azimuthal averaging–reconstruction (ring average) technique to finite-difference schemes in order to mitigate the time step constraint in spherical and cylindrical coordinates. The finite-difference ring average technique averages physical quantities based on an effective grid and then reconstructs the solution back to the original grid in a piecewise, monotonic way. The algorithm is implemented in a community upper-atmospheric model, the Thermosphere–Ionosphere Electrodynamics General Circulation Model (TIEGCM), with a horizontal resolution up to 0.625∘×0.625∘ in geographic longitude–latitude coordinates, which enables the capability of resolving critical mesoscale structures within the TIEGCM. Numerical experiments have shown that the ring average technique introduces minimal artifacts in the polar region of general circulation model (GCM) solutions, which is a significant improvement compared to commonly used low-pass filtering techniques such as the fast Fourier transform method. Since the finite-difference adaption of the ring average technique is a post-solver type of algorithm, which requires no changes to the original computational grid and numerical algorithms, it has also been implemented in much more complicated models with extended physical–chemical modules such as the Coupled Magnetosphere–Ionosphere–Thermosphere (CMIT) model and the Whole Atmosphere Community Climate Model with thermosphere and ionosphere eXtension (WACCM-X). The implementation of ring average techniques in both models enables CMIT and WACCM-X to perform global simulations with a much higher resolution than that used in the community versions. The new technique is not only a significant improvement in space weather modeling capability, but it can also be adapted to more general finite-difference solvers for hyperbolic equations in spherical and polar geometries.
Highlights.
The ring average technique is adapted to solve the issue of clustered grid cells in polar and spherical coordinates with a finite-difference method. The ring average technique is applied to develop a 0.625∘×0.625∘ high-resolution TIEGCM and more complicated geoscientific models with polar and spherical coordinates as well as finite-difference numerical schemes. The high-resolution TIEGCM shows good capability in resolving mesoscale structures in the ionosphere–thermosphere (I–T) system.
Funder
Fundamental Research Funds for the Central Universities China Postdoctoral Science Foundation National Postdoctoral Program for Innovative Talents National Natural Science Foundation of China
Publisher
Copernicus GmbH
Reference48 articles.
1. Basu, S., Basu, S., Sojka, J. J., Schunk, R. W., and MacKenzie, E.: Macroscale modeling and mesoscale observations of plasma density structures in the polar cap, Geophys. Res. Lett., 22, 881–884, 1995. a 2. Bouaoudia, S. and Marcus, P.: Fast and accurate spectral treatment of
coordinate singularities, J. Comput. Phys., 96, 217–223, 1991. a 3. Codrescu, M. V., Fuller-Rowell, T. J., and Foster, J. C.: On the importance of E-field variability for Joule heating in the high-latitude thermosphere,
Geophys. Res. Lett., 22, 2393–2396, 1995. a 4. Colella, P. and Woodward, P. R.: The Piecewise Parabolic Method (PPM) for
gas-dynamical simulations, J. Comput. Phys., 54, 174–201, 1984. a 5. Courant, R., Friedrichs, K., and Lewy, H.: Über die partiellen
Differenzengleichungen der mathematischen Physik, Math. Annal., 100, 32–74, 1928. a
Cited by
20 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献
|
|