Секция 1
33
Голуба - Кахана. Описываются алгоритмы с односторонним и двусторонним предобусловливанием СЛАУ. Исследуются обобщения метода бидиагонализации Голуба - Кахана - Ариоли для решения алгебраических систем с матрицами седлового типа, возникающими в актуальных задачах электромагнетизма, фильтрации и других приложениях. Обсуждаются вопросы масштабируемого распараллеливания средствами гибридного программирования на многопроцессорных вычислительных системах с распределенной и иерархической общей памятью.
The stability properties of symmetric multi-step methods in the predictor-corrector mode
E. D. Karepova, Yu. V. Shan'ko, I. R. Adaev Institute of Computational Modeling SB RAS Email: e.d.karepova@icm.krasn.ru DOI: 10.24412/cl-35065-2022-1-00-19
The orbital motion is described by the system of second-order ordinary differential equations. J. Lambert and I. Watson proposed the symmetric methods [1], that possess a periodicity property when the product of the step-size and the angular frequency lies within a certain interval called the interval of periodicity. The numerical integration of orbit by the symmetric methods with the step-size from the interval of periodicity gives the longitude error which increases linearly, whereas the energy error remains roughly constant, which distinguishes symmetric methods from other one for the better.
The symmetric methods are not uniquely determined even if their order and explicitness are specified. We construct and investigate the high-order symmetric explicit and implicit methods in the "Predict-Evaluate-Cor-rect-Evaluate" (PECE) mode. We propose a technology for testing the stability of predictor-corrector pairs and obtaining their intervals of periodicity or stability. Note that the symmetry is preserved only in the case when the order of the corrector is one greater than the order of the predictor, otherwise the pair should be examined for stability.
This work is supported by the Krasnoyarsk Mathematical Center and financed by the Ministry of Science and Higher Education of the Russian Federation in the framework of the establishment and development of regional Centers for Mathematics Research and Education (Agreement No. 075-02-2022-873).
References
1. J. D. Lambert, I. A. Watson. Symmetric Multistep Methods for Periodic Initial Value Problems // J. Inst. Maths Applics. 1976. V. 18. P. 189-202.
On convergence of shock-capturing schemes inside the shocks influence area
O. A. Kovyrkina12, V. V. Ostapenko12, V. F. Tishkin3
1Lavrentyev Institute of Hydrodynamics SB RAS
2Novosibirsk State University
3Keldysh Institute of Applied Mathematics RAS
Email: olyana@ngs.ru
DOI: 10.24412/cl-35065-2022-1-00-20
We show that in the Nonlinear Flux Correction (NFC) schemes CABARETM [1] and WENO5 [2] (in contrast to the Rusanov scheme [3]), there is no local convergence of the difference solution inside the shock influence area when calculating the shallow water dam break problem. This is due to the fact that the numerical solutions obtained by these schemes have persistent oscillations inside the domain of a constant flow between a shock and a centered rarefaction wave. In this case, taking into account the Lax - Wendroff theorem [4], the numerical solutions obtained by the NFC schemes converge to the exact solution inside the shock influence area in the weak sense only.
The reported study was funded by RFBR and NSFC, project number 21-51-53012.
34 Методы вычислительной алгебры и решения уравнений математической физики
References
1. Kovyrkina O.A., Ostapenko V.V. Monotonicity of the CABARET scheme approximating a hyperbolic system of conservation laws // Comput. Math. Math. Phys. 2018. V. 58 (9). P. 1435-1450.
2. Jiang G.S., Shu C.W. Efficient implementation of weighted ENO schemes // J. Comput. Phys. 1996. V. 126. P. 202-228.
3. Rusanov V.V. Difference schemes of the third order of accuracy for the forward calculation of discontinuous solutions // Dokl. Akad. Nauk SSSR. 1968. V. 180. P.1303-1305.
4. Lax P.D., Wendroff B. Systems of conservation laws // Commun. Pure Appl. Math. 1960. V. 13 (2). P. 217-237.
Использование сеточно-характеристических методов в задачах железнодорожной безопасности транспорта
А. А. Кожемяченко, Е. А. Песня, А. В. Фаворская Московский физико-технический институт Email: anton-kozhemyachenko@yandex.ru DOI: 10.24412/cl-35065-2022-1-00-21
Целью работы является исследование задач, связанных с железнодорожной безопасностью транспорта [1, 2] в условиях тяжеловесного и высокоскоростного движения на различных участках пути и применения в них сеточно-характеристического метода на структурированных сетках. Используя явные се-точно-характеристические методы, рассмотрено движение поезда по балластному и безбалластному мостовому полотну. В результате исследований были получены различные волновые картины и картины динамического распределения компонент тензора напряжений Коши при движении состава по железнодорожному пути, представленного в виде комбинации линейно-упругих сред [3]. Рассматривается конструирование неявного сеточно-характеристического метода первого и второго порядка [4] и возможность его применения для моделирования динамических процессов в акустической и изотропной линейно-упругой средах. В модельных постановках наблюдались как акустические волны с большой скоростью распространения, так и макроволны с малой скоростью распространения.
Исследование выполнено при финансовой поддержке Российского научного фонда (проект № 20-71-10028). Список литературы
1. В. Ю. Поляков, Д. Н. Тхань. Безопасность движения и динамические свойства мостового полотна на ВСМ. Вестник научно-исследовательского института железнодорожного транспорта. 2018. Т. 77, № 6. С. 357-363
2. В. Ю. Поляков, Д. Н. Тхань. Ударное взаимодействие колеса и рельса на мостах высокоскоростных магистралей. Интернет-журнал "Транспортные сооружения". 2019. № 1.
3. Новацкий В. Теория упругости. М.: Мир, 1975.
4. Kozhemyachenko A.A., Pesnya E., Favorskaya A.V. Application of Implicit Grid-Characteristic Methods for Modeling Wave Processes in Linear Elastic Media. In: Czarnowski I., Howlett R.J., Jain L.C. (eds) Intelligent Decision Technologies. Smart Innovation, Systems and Technologies. Springer, Singapore. 2021. V. 238.
Математическое моделирование протяженных электронно-оптических систем
А. Н. Козырев, В. М. Свешников
Институт вычислительной математики и математической геофизики СО РАН
Email: victor@lapasrv.sscc.ru
DOI: 10.24412/cl-35065-2022-1-00-25
Математическое моделирование протяженных электронно-оптических систем с интенсивными пучками приводит к решению самосогласованной нелинейной задачи, включающей в себя расчет электрических и магнитных полей, траекторий заряженных частиц и объемного заряда. Под протяженной