Geodynamo simulations software package suite based on spectral hereditary models

G. M. Vodinchar   ; E. A. Kazakov   ; L. K. Feshenko

Vestnik KRAUNC. Fiziko-matematičeskie nauki, Tome 49 (2024) no. 4, pp. 65-84

Résumé

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 $\alpha$-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.

Détail

Keywords: geodynamo, spectral models, memory models, simulation software packages, symbolicnumerical computations.

Bibliographie

[1] Glatzmaier G. A., Roberts P. H., “A three-dimensional self-consistent computer simulation of a geomagnetic field reversal”, Nature, 377 (1995), 203–209 DOI: 10.1038/377203a0 | DOI

[2] 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, 91:1–3 (1995), 63–75 DOI: 10.1016/0031-9201(95)03049-3 | DOI

[3] Merril R., McElhinny M., McFadden P., The Magnetic Field of the Earth: Paleomagnetism, the Core, and the Deep Mantle, Acad. Press, New York, 1996, 532 pp.

[4] Gledzer E. B., Dolzhanskii F. V., Obukhov A. M., Sistemy gidrodinamicheskogo tipa i ikh primenenie, Nauka, M., 1981, 368 pp. | MR

[5] Monin A. B., Teoreticheskie osnovy geofizicheskoi gidrodinamiki, Gidrometeoizdat, L., 1988, 424 pp.

[6] Vodinchar G. M., Feshchenko L. K., “Model of Geodynamo Driven by Six-jet Convection in the Earth's Core”, Magnetohydrodynamics, 52:1 (2016), 287–299 | DOI | MR

[7] Vodinchar G., Feshchenko L., “Computational Technology for the Basis and Coefficients of Geodynamo Spectral Models in the Maple System”, Mathematics, 11:13 (2023), 3000 DOI: 10.3390/math11133000 | DOI

[8] Krause F., Radler K.-H., Mean-filed magnetohydrodynamics and dynamo theory, Academic-Verlag, Berlin, 1980, 284 pp. | MR

[9] Zeldovich Ya. B., Ruzmaikin A. A., Sokoloff D. D., Magnetic Fields in Astrophysics, Gordon and Breach, New York, 1983, 365 pp.

[10] Brandenburg A., “Memory effects in turbulent transport”, Astrophysical Journal, 706:1 (2009), 712 DOI: 10.1088/0004-637X/706/1/712 | DOI

[11] Tikhonov A. N., Samarskii A. A., Uravneniya matematicheskoi fiziki, Nauka, M., 1972, 736 pp. | MR

[12] L’Ecuyer P., “Random Number Generation”, Handbook of Computational Statistics: Concepts and Methods, Springer, Berlin–Heidelberg, 2012, 35-72 | MR

[13] Vodinchar G., Feshchenko L., “Fractal Properties of the Magnetic Polarity Scale in the Stochastic Hereditary $\alpha\omega$-Dynamo Model”, Fractal and Fractional, 6:6 (2022), 328 DOI: 10.3390/fractalfract6060328 | DOI

[14] Amosov A. A., Dubinskii Yu. A., Kopchenova N. V., Vychislitelnye metody dlya inzhenerov, Vysshaya shkola, M., 1994, 544 pp.

[15] Vodinchar G. M., Kazakov E. A., “Isklyuchenie integralnogo chlena v uravneniyakh odnoi ereditarnoi sistemy, svyazannoi s zadachei gidromagnitnogo dinamo”, Vestnik KRAUNTs. Fiz.-mat. nauki, 42:1 (2023), 180–190 DOI: 10.26117/2079-6641-2023-42-1-180-190 | DOI | MR | Zbl