Vestnik КRAUNC. Fiz.-Mat. nauki. 2024. vol. 49. no. 4. P. 65-84. ISSN 2079-6641
INFORMATION AND COMPUTING TECHNOLOGIES
https://doi.org/10.26117/2079-6641-2024-49-4-65-84
Research Article
Full text in Russian
MSC 68-04, 45K05, 76W05
Geodynamo Simulations Software Package Suite Based on Spectral Hereditary Models
G. M. Vodinchar^{\ast}, E. A. Kazakov, L. K. Feshchenko
Institute of Cosmophysical Research and Radio Wave Propagation FEB RAS, 684034, Paratynka, Kamchatka, Mirnaya str., 7, Russia
Abstract. The study of the geodynamics problem is often carried out on the basis of spectral models, when the fields are fully or partially expanded in eigenfields (eigenmodes) of suitable spectral problems. The most meaningful from the physical point of view are spectral problems of free oscillations or free decay of fields. The compilation of spectral models first of all requires calculating the parameters of the basic modes, and then the model coefficients. Most often, these are the Galerkin coefficients. Then the problem of actually solving the model equations numerically arises. The paper describes a software package developed by the authors that allows solving such problems. It includes modules for calculating the mode parameters, a module for calculating the Galerkin coefficients, two modules for numerically solving the system, and a noise generation module. The package allows calculating the model with hereditary quenching of the α-effect by the field energy. Two types of the quenching functional kernel are provided, requiring different difference schemes. These schemes are implemented in two numerical solution modules. Random noise simulates the effect of spontaneous synchronization of small-scale field components, which is absent on average. The calculation of the parameters of the basic modes and Galerkin coefficients is performed using combined symbolic-numerical computations, so the corresponding modules are implemented in the Maple package. The need for symbolic computations is associated with the great complexity of the expressions of the modes themselves and the integrands when calculating the Galerkin coefficients. Therefore, the task arises first of all to form the necessary expressions. This is done using symbolic calculations. The remaining modules are implemented in C++. The developed package can be useful for specialists studying the geodynamo problem based on spectral models and memory effects in this problem.
Key words: geodynamo, spectral models, memory models, simulation software packages, symbolicnumerical computations.
Received: 05.11.2024; Revised: 19.11.2024; Accepted: 26.11.2024; First online: 28.11.2024
For citation. Vodinchar G. M., Kazakov E. A., Feshchenko L. K. Geodynamo simulations software package suite based on spectral hereditary models. Vestnik KRAUNC. Fiz.-mat. nauki. 2024, 49: 4, 65-84. EDN: QYHFFU. https://doi.org/10.26117/2079-6641-2024-49-4-65-84.
Funding. This research was funded by Russian Science Foundation grant number 22-11-00064 https://rscf.ru/project/22-11-00064/.
Competing interests. The authors declare that there are no conflicts of interest regarding authorship and publication.
Contribution and Responsibility. All authors contributed to this article. Authors are solely responsible for providing the final version of the article in print. The final version of the manuscript was approved by all authors.
^{\ast}Correspondence: E-mail: gvodinchar@ikir.ru
The content is published under the terms of the Creative Commons Attribution 4.0 International License
© Vodinchar G. M., Kazakov E. A., Feshchenko L.K., 2024
© Institute of Cosmophysical Research and Radio Wave Propagation, 2024 (original layout, design, compilation)
References
- Glatzmaier G. A., Roberts P. H. A three-dimensional self-consistent computer simulation of a geomagnetic field reversal, Nature, 1995, vol. 377, pp. 203–209. DOI: 10.1038/377203a0
- Glatzmaier G. A., Roberts P. H. A three-dimensional convective dynamo solution with rotating and finitely conducting inner core and mantle, Physics of the Earth and Planetary Interiors, 1995, vol. 91, no. 1–3, pp. 63–75. DOI: 10.1016/0031-9201(95)03049-3
- Merril R., McElhinny M., McFadden P. The Magnetic Field of the Earth: Paleomagnetism, the Core, and the Deep Mantle. New York, Acad. Press, 1996, 532 p.
- Gledzer E. B., et. al. Sistemy gidrodinamicheskogo tipa i ih primenenie [Hydrodynamic type systems and their application]. Moscow, Nauka, 1981, 368 p. (In Russian).
- Moinin A. B. Teoreticheskie osnovy geofizicheskoj gidrodinamiki [Theoretical foundations of geophysical fluid dynamics]. Leningrad, Gidrometeoizdat, 1988, 424 p. (In Russian).
- Vodinchar G. M., Feshchenko L. K. Model of Geodynamo Driven by Six-jet Convection in the Earth’s Core, Magnetohydrodynamics, 2016, vol. 52, no. 1, pp. 287–299.
- Vodinchar G., Feshchenko L. Computational Technology for the Basis and Coefficients of Geodynamo Spectral Models in the Maple System, Mathematics, 2023, vol. 11, no. 13, 3000.
- Krause F., R¨adler K.-H. Mean-filed magnetohydrodynamics and dynamo theory. Berlin, Academic-Verlag, 1980, 284 p.
- Zeldovich Ya. B., Ruzmaikin A. A., Sokoloff D. D. Magnetic Fields in Astrophysics. New York, Gordon and Breach, 1983, 365 p.
- Brandenburg A. Memory effects in turbulent transport, Astr. Jour., 2009, 706: 1, 712.
- Tihonov A. N., Samarskii A. A. Equations of Mathematical Physics. New York, Dover Publications, 2013, 800 p.
- L’Ecuyer P. Random Number Generation. In: Handbook of Computational Statistics: Concepts and Methods. Berlin–Heidelberg, Springer, 2012, pp. 35-72-112
- Vodinchar G., Feshchenko L. Fractal Properties of the Magnetic Polarity Scale in the Stochastic Hereditary αω-Dynamo Model, Fractal Fract, 2022, vol. 6, no. 6, 328.
- Amosov A. A., et. al. Vychislitel’nye metody dlya inzhenerov [Computational Methods for Engineers]. Moscow, Vysshaya shkola, 1994, 544 p. (In Russian).
- Vodinchar G. M., Kazakov E. A. Elimination of the integral term in the equations of one hereditary system related to the problem of hydromagnetic dynamo, Vest. KRAUNC. Fiz.-Mat. Nauki, 2023, 42:1, 180–190. DOI: 10.26117/2079-6641-2023-42-1-180-190 (In Russian).
Information about the authors
Vodinchar Gleb Mikhailovich – PhD (Phys & Math), Associate Professor, Leading Researcher, Institute of Cosmophysical Research and adio Wave Propagation FEB RAS, Paratunka, Kamchatka, Russia, ORCID 0000-0002-5516-1931.
Kazakov Evgeny Anatolyevich – Junior Researcher, Institute of Cosmophysical Research and adioWave Propagation FEB RAS, Paratunka, Kamchatka, Russia, ORCID 0000-0001-7235-4148.
Feshchenko Liybov Konstantinovna – PhD (Phys & Math), Researcher, Institute of Cosmophysical Research and adioWave Propagation FEB RAS, Paratunka, Kamchatka, Russia, ORCID 0000-0001-5970-7316.