=Paper=
{{Paper
|id=Vol-1482/611
|storemode=property
|title=Решение задач газовой динамики горения на гибридных вычислительных системах
(Solution of problems of combustion gas dynamics on hybrid CPU/GPU computing system)
|pdfUrl=https://ceur-ws.org/Vol-1482/611.pdf
|volume=Vol-1482
}}
==Решение задач газовой динамики горения на гибридных вычислительных системах
(Solution of problems of combustion gas dynamics on hybrid CPU/GPU computing system)==
https://ceur-ws.org/Vol-1482/611.pdf
Суперкомпьютерные дни в России 2015 // Russian Supercomputing Days 2015 // RussianSCDays.org
Решение задач газовой динамики горения на гибридных
вычислительных системах
Л.И. Стамов, Е.В. Михальченко
Федеральное государственное учреждение «Федеральный научный центр Научно-
исследовательский институт системных исследований Российской академии наук»
Задачи горения и детонации являются чрезвычайно интересными и важными для дальней-
шего развития современного мира, особенно интересны переходы между этими режимами. Та-
кого рода задача о переходе горения в детонацию является до сих пор наименее изученной, а ее
решение, например, может способствовать созданию двигателей нового типа. Сложность изу-
чения такого типа процессов представляет отсутствие возможности решить эти задачи анали-
тически. В связи с чем часто для моделирования поведения таких систем применяют вычисли-
тельный эксперимент.
В данной работе рассмотрена возможность применения ряда схем высокого порядка точно-
сти для решения многомерных задач горения и детонации. Рассмотрены численные схемы типа
PPM (piece parabolic method) [1] и ряд центральных схем [2-4].Особенностью методов типа PPM
является использование параболы в качестве интерполяционной функций, что позволяет полу-
чать более точное представление о поведении пространственных градиентов, а также, лучше
описывает разрывы, особенно контактные разрывы. Главной особенностью второго семейства
схем является их простота. Выбранные центрированные схемы не требует каких-либо решате-
лей Римана, которые используются в противопотоковых схемах высокого порядка точности, а
также в них нет необходимости вычислять точный якобиан системы уравнений.
Для моделирования протекания химических реакций была выбрана модель нестационарно-
го процесса горения в химически реагирующей однородной газовой среде [5]. В расчетах ис-
пользовался кинетический механизм горения кислородно-воздушной смеси Мааса и Варнаца
[6]. Решение системы химических уравнений осуществлялось с помощью четырех стадийного
метода типа Розенброка [7].
Все численные схемы были протестированы с помощью ряда тестовых задач, были сравне-
ны с аналитическими решениями, где это представлялось возможным. Результаты показали
хорошее совпадение численных решений с аналитическими в рамках выбранных схем. Однако,
стоит отметить, что, например, схема PPM [1] показала меньшее размазывание ударных и кон-
тактных поверхностей по сравнению с рассмотренными центральными схемами [2-3]. В то же
время, данная схема требует значительно большего числа вычислений, чем центрированные
схемы.
В работе для сокращения времени моделирования были построены и рассмотрены парал-
лельные реализации выбранных численных схем. Для их написания использовались технологии
OpenMP [8] и CUDA [9]. Распараллеливание производилось на основе геометрической деком-
позиции, т.е. расчетная область делилась поровну между имеющимися в системе вычислитель-
ными ресурсами. Было рассмотрено влияние размера вычислительной сетки и размеров блоков
на производительность. Было получено, что с использованием графических процессоров, в за-
висимости от их технических характеристик и от количества узлов расчетной сетки, можно по-
лучить более чем 10-кратное ускорение по сравнению с центральным процессором. Было пока-
зано, что с использованием выбранных технологий параллельного программирования можно
значительно ускорить время проведения моделирования задач физико-химической газовой ди-
намики на многопроцессорных системах и системах с графическими сопроцессорами.
Литература
1. Colella P., Woodward P. R. The Piecewise Parabolic Method (PPM) for Gas-Dynamical Simula-
tions // Journal of Computational Physics. 1984. Vol. 54, pp. 174-201.
611
Суперкомпьютерные дни в России 2015 // Russian Supercomputing Days 2015 // RussianSCDays.org
2. Jiang G.-S., Tadmor E. Nonoscillatory central schemes for multidimensional hyperbolic conserva-
tion laws // SIAM J. SCI. COMPUT. 1988. Vol. 19, pp. 1892-1917.
3. Kurganov A., Tadmor E. New High-resolution central schemes for nonlinear conservation laws
and convection–diffusion equations // Journal of Computational Physics. 2000. Vol. 160, pp. 241-
282.
4. Kurganov A., Levy D. A third-order semidiscrete central scheme for conservation laws and con-
vection-diffusion equations // SIAM J. SCI. COMPUT. 2000. Vol. 22, No. 2, pp. 1461-1488.
5. Cмирнов Н.Н., Никитин В.Ф., Алиари Шурехдели Ш. Переходные режимы распростране-
ния волн в метастабильных системах // Физика горения и взрыва. 2008. 44, 5. c. 25-37.
6. Maas U., Warnatz J. Ignition Processes in Hydrogen-Oxygen Mixtures // Combustion and Flame.
1988. Vol. 74, pp. 53-69.
7. Хайрер Э., Ваннер Г. Решение обыкновенных дифференциальных уравнений. Жесткие и
дифференциально-алгебраические задачи. Пер. с англ. М., 1999.
8. Антонов А.С. Технологии параллельного программирования MPI и OpenMP: Учеб. посо-
бие. Предисл.: В.А.Садовничий. М., 2012.
9. NVIDIA CUDA. Programming Guide. 2015. URL: http://developer.nvidia.com/cuda-downloads (дата
обращения: 27.05.2015).
612
Суперкомпьютерные дни в России 2015 // Russian Supercomputing Days 2015 // RussianSCDays.org
Solution of problems of combustion gas dynamics on hybrid
CPU/GPU computing system
Lyuben Stamov and Elena Mikhalchenko
Keywords: combustion, detonation, gas dynamics, parallel algorithms, GPU
In this paper an application of modern GPU containing computing system for computation of
combustion problems is study. Some test problems of initiation of detonation in hydrogen-air
mixtures were considered. Several parallel algorithms based on GPU-computing were
developed. High-order difference schemes for gas dynamics and chemistry equations were
used. Good performance for a certain range of parameters was obtained.