Комплекс программ моделирования геодинамо на основе спектральных эредитарных моделей
Authors: Г.М. Водинчар , Е.А. Казаков , Л.К. Фещенко 
Published: November 2024

Исследование задачи геодинамо часто проводят на основе спектральных моделей, когда поля задачи полностью или частично раскладываются по собственным полям (модам) подходящих спектральных задач. Наиболее осмысленными с физической точки зрения являются спектральные задачи о свободных колебаниях или свободном затухании полей. Составление спектральных моделей прежде всего требует расчета параметров базисных мод, а затем и коэффициентов модели. Чаще всего это коэффициенты Галеркина. Затем возникает задача собственно численного решения уравнений модели. В работе описывается разработанный авторами комплекс программ, позволяющий решать такие задачи. Он включает в себя модули расчета параметров мод, модуль расчета коэффициентов Галеркина, два модуля численного решения системы, модуль генерации шума. Комплекс позволяет вести расчет модели с эредитарным подавлением α-эффекта энергией поля. Предусмотрено два типа ядер функционала подавления, требующих различных разностных схем. Эти схемы реализованы в двух модулях численного решения. Случайный шум имитирует влияние спонтанной синхронизации мелкомасштабных компонент полей, которое в среднем отсутствует. Расчет параметров базисных мод и коэффициентов Галеркина выполняется с помощью комбинированных символьно-численных вычисления, поэтому соответствующие модули реализованы в пакете Maple. Необходимость символьных вычислений связана с большой сложностью выражений самих мод и подынтегральных выражений при расчете коэффициентов Галеркина. Поэтому возникает задача прежде всего формировать необходимые выражения. Это и делается с помощью символьных вычислений. Остальные модули реализованы на C++. Разработанный комплекс может быть полезен специалистам, изучающим задачу геодинамо на основе спектральных моделей и эффекты памяти в этой задаче.

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.