Технология решения задач физики плазмы на суперЭВМ
(По материалам докторской диссертации)
Семинар: Информационно-вычислительные технологии
Начало заседания: 16:00
Дата выступления: 18 Октябрь 2016
Организация: ИВМиМГ СО РАН (Новосибирск)
Авторы: к.ф.-м.н. Снытников Алексей Владимирович
В работе представлена суперкомпьютерная технология, созданная для решения актуальных задач физики плазмы, таких как моделирование плазменно-пучкового взаимодействия, моделирование взаимодействия лазерного импульса с плазмой, солнечные вспышки II и III типа и др., в рамках кинетического подхода. Основными элементами технологии являются эйлерово-лагранжева декомпозиция расчетной области, реализация наиболее трудоемких расчетных процедур на ускорителях вычислений, параметризованная форма реализации метода частиц в ячейках, и методика межархитектурного переноса программ.
Эйлерово-лагранжева декомпозиция расчетной области представляет собой сочетание эйлеровой декомпозиции, когда пространство моделирования разделяется на части между процессорами, и лагранжевой декомпозиции, когда модельные частицы распределяются между процессорами вне зависимости от положения в расчетной области. Для комбинированной эйлерово-лагранжевой декомпозиции вначале происходит разделение пространства моделирования на подобласти вместе с частицами, а затем частицы каждой подобласти разделяются между процессорами. Таким образом достигнута эффективность распараллеливания (в слабом смысле) 92 % для 500 графических ускорителей Nvidia Tesla на суперкомпьютере «Ломоносов», НИВЦ МГУ.
Алгоритм, реализующий метод частиц в ячейках, и включающий в себя расчет электромагнитного поля по методу FDTD и расчет движения частиц с использованием схемы Бориса, был перенесен на графические ускорители. При этом было достигнуто ускорение для стадии расчета движения модельных частиц в 40 раз для графического ускорителя Nvidia Kepler по сравнению с 4-ядерным процессором Intel Xeon.
Создана параметризованная реализация метода частиц в ячейках для суперЭВМ гибридной архитектуры. Тем самым достигнута возможность замены различных вычислительных процедур метода частиц, а также адаптации кода к особенностям архитектуры суперЭВМ, возможность изменения размерности задачи, а также учета ее физических особенностей.
Методика межархитектурного переноса программ обеспечивает возможность переноса алгоритма, реализованного для графических ускорителей с применением технологии CUDA на ускорители Intel Xeon Phi. Данная методика позволяет проводить расчеты по моделированию динамики плазмы без изменения кода на суперЭВМ всех архитектур, представленных в первой десятке наиболее мощных компьютеров мира.