, , , ,

"Решение кинетического уравнения Больцмана на графических ускорителях". Сергей Олегович Полешкин, магистратура НГУ, 23.7.2014

nsu_sopoleshkin_201407.pdf

Состав коллектива

Научное содержание работы

'Признанным методом расчета течений разреженного газа является метод прямого статистического моделирования (ПСМ), предложенный в работах Г. Берда и заметно усовершенствованный в работах М. С. Иванова и С. В. Рогазинского, обосновавших, так называемую, схему «мажорантной частоты». На основе метода ПСМ были разработаны программные комплексы, эффективно применяемые при расчетах высокоскоростных течений. Однако при расчетах медленных течений и нестационарных течений этот метод является неэффективным из-за статистического шума. В связи с этим, представляется важным подход к моделированию течений разреженного газа, основанный на детерминированных методах решения уравнения Больцмана. Начало разработки этих методов было положено в конце 70-х годов в работах В. В. Аристова и Ф. Г. Черемисина. Но только в последние 5-10 лет детерминированный подход стал бурно развиваться. Это связано с удешевлением и широким распространением мощных вычислительных систем, позволяющих решать нетривиальные задачи детерминированными методами, которые требуют очень больших вычислительных затрат. Узким местом при численном решении кинетического уравнения Больцмана является вычисление интеграла столкновений, оно занимает более 90 процентов времени выполнения программы. Диссертационная работа посвящена алгоритмической и программной оптимизации вычисления интеграла столкновений в одной из схем детерминированного метода решения уравнения Больцмана. Ввиду вышеизложенного понятна актуальность темы работы. Предложен и программно реализован адаптивный метод вычисления интеграла столкновений на основе сокращения количества узлов расчетной сетки в скоростном пространстве, отбираемых в качестве узлов квадратуры. Им были проведены тестовые расчеты для различных классов задач, которые показали примерно трехразовое ускорение. Проведенная оптимизация позволила в сокращенные сроки провести серию расчетов профиля ударной волны при больших числах Маха. Эти расчеты подтвердили существование в профиле особенности, локального максимума температуры вниз по потоку. Возможность существования этой особенности долгое время обсуждалась в литературе и полученный результат можно считать фундаментальным.'

Публикации