Application of the Alternating Triangular Iterative Method to Solving Shallow Water Hydrodynamics Problems with Graphics Accelerators

Процесс прогнозирования экологических катастроф как техногенного, так и природного характера в настоящее время базируется на достижениях в области математического моделирования. Получение прогноза за короткое время является весьма затруднительным без применения параллельных вычислений и суперкомпьютерных технологий. Большой объем обрабатываемой информации и сложность вычислений приводят к необходимости использования вычислительных кластеров, в состав которых добавляются видеоадаптеры для увеличения производительности вычислительной системы и скорости обработки информации. Статья посвящена разработке программного модуля на языке CUDA C, предназначенного для моделирования гидродинамических процессов мелководных водоёмов, включая решение возникающих при дискретизации систем сеточных уравнений высокой размерности. Для повышения эффективности расчетов часть вычислительной нагрузки передавалась на графические ускорители, для чего был разработан алгоритм параллельных расчетов и его программная реализация в виде модуля. В результате вычислительных экспериментов определена оптимальная двумерная конфигурация потоков в вычислительном блоке, исполняемом на одном потоковом мультипроцессоре. Разработанные алгоритм и программный модуль позволяют более эффективно задействовать вычислительные ресурсы графических ускорителей, используемых для решения вычислительно трудоемких задач гидрофизики. Forecasting environmental disasters, both artificial and natural, is currently based on advances in simulation. Fast forecasting is very difficult without the use of parallel computing and supercomputer technologies. The large volume of information to be processed and the complexity of calculations require to use of computing clusters with graphics cards to increase the computing performance and the data processing rate. This paper deals with the development of a software module in CUDA C to simulate the hydrodynamic processes in shallow water bodies, including the solution of systems of high-dimensional grid equations during discretization. To improve computing performance, graphics accelerators take a part of the load. To enable this, we developed an algorithm for parallel calculations and implemented it as a software module. As a result of computational experiments, the optimal 2D configuration of streams in a computational unit run on a single streaming multiprocessor was determined. The proposed algorithm and software module enables more efficient use of GPU computational resources when solving computationally intensive hydro physics problems.