Parallel simulation of kinetic processes by Monte Carlo method (dedicated to the memory of the Chief Theoretiсian of Cosmonautics academician M.V. Keldysh in the 60th anniversary of the launch of the first artificial sputnik)
Authors: Марченко М.А. , Сушкевич Т.А. 
Published: August 2021

В 2017 году мировая общественность отмечает 60-летний юбилей запуска 4 октября 1957 года в СССР первого искусственного спутника Земли, положившего начало космической эры. Баллистические расчеты проводились на первой серийной ЭВМ Стрела в Институте Келдыша. При решении сложнейших задач создания ракетно-ядерного щита были заложены основы новых направлений в математике - вычислительной математики и математического моделирования. И в СССР, и в США параллельно разрабатывались детерминированные и статистические численные методы. Методы Монте-Карло (ММК) как инструмент для прямого статистического моделирования были разработаны в США в рамках Манхэттонского проекта создания ядерного оружия. Джон фон Нейман первым предложил использовать саму ЭВМ для генерации случайных чисел. В 1949 году Джон фон Нейман и Станислав Улам предложили первый алгоритм получения псевдослучайных величин, который впоследствии был назван ММК и послужил основой для развития методики генерации псевдослучайных чисел с использованием ЭВМ. Разработка ММК и эффективность его применения всегда начинается с разработки генератора случайных или псевдослучайных чисел, который зависит от класса решаемых задач и конкретной структуры и архитектуры ЭВМ. Методы Монте-Карло стали массово применять на всех архитектурах вычислительных систем с параллельными и распределенными вычислениями. Сейчас в эпоху супервычислений преобладают ММК как следствие простоты их
реализации. Но эта простота обманчива. В статье представлен разработанный отечественный комплексный методический подход, в котором на примере трех сложных больших задач, описывающих
пространственно-неоднородные кинетические процессы диффузии, коагуляции и переноса заряженных частиц, системно рассматриваются теория методов и алгоритмов ММК и практика их реализации в формате не просто программ, а также параллельных генераторов псевдослучайных чисел, библиотек программ, средств обработки
данных, управляющих программ и т.д., т.е. все этапы создания цифрового продукта. На примере вероятностных моделей для численного моделирования кинетических процессов диффузии, коагуляции и переноса заряженных частиц, когда ансамбли траекторий или частиц содержат по 10 в 7-13 степени элементов, продемонстрированы возможности и эффективность новых параллельных алгоритмов и распределенных вычислений ММК для решения больших и сложных задач не только для расчета отдельных функционалов или оценок, но и для всего фазового объема задачи. Это важнейшее достижение, которое повышает конкурентность ММК с детерминированными конечно-разностными и сеточными методами при параллельном моделировании.

In 2017, the global community celebrates the 60th anniversary of the launch on 4 October 1957 in USSR of the first Earth sputnik, which ushered the space age. The ballistic calculations were conducted on the first production computer Strela in the Keldysh Institute. When solving the most complicated problems of creating a nuclear missile shield were laid the foundations for new directions in mathematics - computational mathematics and mathematical modeling. The USSR and the USA simultaneously developed the deterministic and statistical numerical methods. Methods Monte Carlo (MMC) as a tool for direct statistical simulation has been developed in the United States under the Manhattan project to create nuclear weapons. John von Neumann first proposed to use the computers to generate random numbers. In 1949, John von Neumann and Stanislaw Ulam proposed the first algorithm generating pseudo-random values, which was subsequently named MMC and was the basis for the development of the technique of generating pseudo-random numbers using a computer. The development of MMC and the effectiveness of its application always begins with the development of random or pseudo-random numbers, which depends on the class of tasks and the specific structure and architecture of computers. Methods Monte Carlo began to apply on all architectures of computing systems with parallel and distributed computing. Now in the era of superficially dominated by MMC as a result of the simplicity of their implementation. But this simplicity is deceptive. Developed a complex methodological approach on the example of three complicated big problems describing spatially inhomogeneous kinetic processes of diffusion, coagulation, and migration of charged particles, systematically discusses the theory of methods and algorithms for MMC and practice implementing them in the format of not just programs, as well as the parallel random number generators, libraries, data processing, control programs, etc., i.e., all the stages of creating a "digital product". For example, probabilistic models for numerical simulation of the kinetic processes of diffusion, coagulation, and migration of charged particles when the ensembles of trajectories or particles contain 10 in 7-13 degrees of elements, and demonstrated the potential and effectiveness of the MMC to resolve big and complex tasks not only for the calculation of particular functionals or estimates, but for the whole phase volume of the task. This is the most important achievement that increases MMC's competitiveness with the deterministic finite-difference and grid methods for parallel simulation.