Трехмерное вычислительное моделирование горения, детонации и переходных процессов на многопроцессорных супер-ЭВМНИР

Three-dimensional numerical simulation of combustion, detonation and deflagration to detonation transition processes on supercomputers

Источник финансирования НИР

грант РФФИ

Этапы НИР

# Сроки Название
1 13 апреля 2017 г.-31 декабря 2017 г. Трехмерное вычислительное моделирование горения, детонации и переходных процессов на многопроцессорных супер-ЭВМ (этап 1, 2017 год)
Результаты этапа: Было рассмотрено моделирование процессов турбулентного горения для двух типов сред: гомогенных горючих смесей газов, в которых горючее и окислитель перемешаны на уровне одной фазы, а также гетерогенных горючих систем, в которых горючее и окислитель предварительно находятся в разных фазах. Методом исследования являлось численное моделирование на основе механики сплошных сред с учетом физических и химических превращений, а также с учетом подсеточных моделей турбулентности. В качестве вычислительных схем были использованы две различные численные схемы сквозного счета: 1) Явная схема 2-го порядка точности по пространству и времени на основе MUSCL-интерполяции переменных на грани при расчете конвективных потоков. Для выбора направления интерполяции и для интерполяции давления применяется метод AUSMP. Схема реализована на регулярной сетке из одинаковых элементов (прямоугольных параллелепипедов), состыковывающихся в произвольной, в общем случае, топологии. Текст программы на языке программирования C. Параллельность исполнения - метод OpenMP. 2) Явная схема 2-го порядка точности по пространству и 2-го по времени на основе схемы Курганова – Тадмора. Использовалание регулярной сетки, состоящей из прямоугольных параллелепипедов. Программа написана на языке CUDA, для ЭВМ гибридной архитектуры. 3) Разработана методика тестирования программ и сравнения с экспериментом. Полученные в 2017 году результаты основывались на научном заделе 2016 года, созданном в рамках работы по гранту 14-03-00171. Результаты были опубликованы в 2017 году со ссылкой на поддержку РФФИ. Участие в научных мероприятиях по тематике Проекта • • 2017 Developing principles for new generation of engines operating cycle (Устный) • Автор: Smirnov N.N. • 1-st Korea-Russia Science and Technology Day, Москва, Россия, 5 июня 2017 • • 2017 Three-dimensional numerical simulation of combustion, detonation and deflagration to detonation transition processes (Устный) • Автор: Smirnov N.N. • GINZBURG CENTENNIAL CONFERENCE on PHYSICS , Москва Физический институт РАН им Лебедева П.Н., Россия, 29 мая - 3 июня 2017 • • 2017 Неравновесные процессы в метастабильных средах (Пленарный) • Автор: Смирнов Н.Н. • Международный симпозиум «Неравновесные процессы в сплошных средах», Пермь, Россия, 15-18 мая 2017
2 1 января 2018 г.-31 декабря 2018 г. Трехмерное вычислительное моделирование горения, детонации и переходных процессов на многопроцессорных супер-ЭВМ (этап 2, 2018 год)
Результаты этапа: 1. Результаты тестирования программ вычислений путем сравнения результатов расчетов с экспериментальными данными по распространению волн в сложных геометрических конфигурациях. При этом проведено сравнение теоретических результатов, полученных с использованием двух разработанных вычислительных кодов с экспериментальными данными и между собой. Проведено сравнение результатов расчетов, при которых использованы две различные численные схемы сквозного счета: 1) Явная схема 2-го порядка точности по пространству и времени на основе MUSCL-интерполяции переменных на грань при расчете конвективных потоков. Для выбора направления интерполяции и для интерполяции давления применяется метод AUSMP. Схема будет реализована на регулярной сетке из одинаковых элементов (прямоугольных параллелепипедов), состыковывающихся в произвольной, в общем случае, топологии. Текст программы на языке программирования C. Параллельность исполнения - метод OpenMP. 2) Явная схема 2-го порядка точности по пространству и 2-го по времени на основе схемы Курганова – Тадмора. Использовалание регулярной сетки, состоящей из прямоугольных параллелепипедов. Программа написана на языке CUDA, для ЭВМ гибридной архитектуры. 2. Произведено уточнение моделей химических взаимодействий и турбулентности на основании сравнения теоретических и экспериментальных данных. 3. Получены новые результаты при исследовании закономерностей переходных процессов при фокусировке ударных волн внутри конуса и клина. Проведено трехмерное вычислительное моделирование горения, детонации и переходных процессов на многопроцессорных супер-ЭВМ. В настоящее время математическое моделирование таких процессов требует высокой степени пространственного разрешения, сравнимого с толщиной зоны реакции. Это делает расчеты реальных процессов развития детонации или ее подавления в областях сложной геометрии трудно выполнимыми в течение разумного времени. Поэтому применение для этой цели параллельных вычислений на суперЭВМ представляется весьма перспективным. С другой стороны моделирование таких сложных взаимозависимых процессов, как взаимодействие химических, физических и газодинамических процессов в турбулентных потоках не поддается распараллеливанию "напрямую", а производительность вычислительного комплекса не прямо пропорциональна числу подключаемых процессоров. Методом исследования являлось численное моделирование на основе механики сплошных сред с учетом физических и химических превращений, а также с учетом подсеточных моделей турбулентности. В качестве вычислительных схем была использована явная схема сквозного счета 2-го порядка точности по пространству и времени на основе MUSCL-интерполяции переменных на грани при расчете конвективных потоков. Для выбора направления интерполяции и для интерполяции давления применяется метод AUSMP. Схема реализована на регулярной сетке из одинаковых элементов (прямоугольных параллелепипедов), состыковывающихся в произвольной, в общем случае, топологии. Текст программы на языке программирования C. Параллельность исполнения - метод OpenMP. Проведено широкомасштабное тестирование разработанных вычислительных комплексов путем сравнения результатов расчетов модельных задач с точными решениями задач о распаде произвольного разрыва в реагирующем газе, а также путем сравнения результатов расчетов с экспериментальными данными по распространению волн в сложных геометрических конфигурациях. Проведено уточнение моделей химических взаимодействий и турбулентности, исходя из сравнения результатов с экспериментальным базисом неравновесных переходных процессов в смесях водорода с воздухом при отражении ударных волн от внутренних поверхностей клина и конуса.
3 1 января 2019 г.-31 декабря 2019 г. Трехмерное вычислительное моделирование горения, детонации и переходных процессов на многопроцессорных супер-ЭВМ (этап 3, 2019 год)
Результаты этапа: 1. Сделан обзор нестационарных переходных режимов распространения самоподдерживающихся экзотермических волн в метастабильных средах различных типов. А именно, рассматривается переход от одного самоподдерживающегося режима к другому. К таким процессам относится переход волн дефлаграции в волны детонации в гомогенных и гетерогенных горючих смесях. Под режимами дефлаграции и детонации в широком смысле будем подразумевать распространение самоподдерживающихся экзотермических волн в среде, когда различие механизмов распространения приводит к существенным различиям структуры волн, дозвуковой и сверхзвуковой скоростям их распространения. 2. Проведено трехмерное вычислительное моделирование горения, детонации и переходных процессов на многопроцессорных супер-ЭВМ. В настоящее время математическое моделирование таких процессов требует высокой степени пространственного разрешения, сравнимого с толщиной зоны реакции. Это делает расчеты реальных процессов развития детонации или ее подавления в областях сложной геометрии трудно выполнимыми в течение разумного времени. Поэтому применение для этой цели параллельных вычислений на суперЭВМ представляется весьма перспективным. Показано, что моделирование таких сложных взаимозависимых процессов, как взаимодействие химических, физических и газодинамических процессов в турбулентных потоках не поддается распараллеливанию "напрямую", а производительность вычислительного комплекса не прямо пропорциональна числу подключаемых процессоров. 3. Проведено широкомасштабное тестирование разработанных вычислительных комплексов путем сравнения результатов расчетов модельных задач с точными решениями задач о распаде произвольного разрыва в реагирующем газе, а также путем сравнения результатов расчетов с экспериментальными данными по распространению волн в сложных геометрических конфигурациях. Проведено уточнение моделей химических взаимодействий и турбулентности, исходя из сравнения результатов с экспериментальным базисом неравновесных переходных процессов в смесях водорода с воздухом при отражении ударных волн от внутренних поверхностей клина и конуса. 4. Разработанная математическая модель перехода горения в детонацию в гомогенных и гетерогенных смесях, позволила изучить все характерные стадии процесса: зажигание, ускоряющееся турбулентное горение, образование горячих точек перед фронтом пламени, возникновение детонационных и ретонационных волн, пересжатые режимы и выход на самоподдерживающуюся детонацию. В отличие от существующих в мире вычислительных кодов достигнуто не только качественное, но и количественное совпадение результатов моделирования переходных процессов с экспериментальными данными. 5. Показано, что возникновение детонации за счет фокусировки сильной ударной волны происходит в режиме пересжатой детонации. Затем детонационная волна замедляется до самоподдерживающегося режима. 6. Переходный режим имеет следующие характерные этапы. Сначала за ударной волной появляется горячая точка, которая впоследствии быстро увеличивается во всех направлениях, тем не менее отставая от ударной волны. Таким образом формируется волна горения. Во-вторых, зона горения увеличивается и фронт волны становится неустойчивым. Морщины образуются на переднем фронте, что особенно заметно в плоскости, ортогональной клину. В-третьих, скорость волны горения резко возрастает, особенно в сморщенной зоне, морщины становятся глубже, образуются карманы несгоревшей смеси, окруженные горячими продуктами реакции. Далее в одном из карманов вблизи оси возникает “взрыв во взрыве”, что ведет к формированию детонационной волны. Детонационная волна начинает отражаться от боковых стенок цилиндра и далее распространяться в виде детонационных и ретонационных волн во всех направлениях.

Прикрепленные к НИР результаты

Для прикрепления результата сначала выберете тип результата (статьи, книги, ...). После чего введите несколько символов в поле поиска прикрепляемого результата, затем выберете один из предложенных и нажмите кнопку "Добавить".