КАРТОЧКА ПРОЕКТА ФУНДАМЕНТАЛЬНЫХ И ПОИСКОВЫХ НАУЧНЫХ ИССЛЕДОВАНИЙ,
ПОДДЕРЖАННОГО РОССИЙСКИМ НАУЧНЫМ ФОНДОМ

Информация подготовлена на основании данных из Информационно-аналитической системы РНФ, содержательная часть представлена в авторской редакции. Все права принадлежат авторам, использование или перепечатка материалов допустима только с предварительного согласия авторов.

 

ОБЩИЕ СВЕДЕНИЯ


Номер 22-21-00902

НазваниеАсимптотическое и численное моделирование динамики детонационных волн

РуководительКасимов Аслан Рамазанович, кандидат наук (признаваемый в РФ PhD)

Организация финансирования, регион Автономная некоммерческая образовательная организация высшего образования «Сколковский институт науки и технологий», г Москва

Период выполнения при поддержке РНФ 2022 г. - 2023 г. 

Конкурс№64 - Конкурс 2021 года «Проведение фундаментальных научных исследований и поисковых научных исследований малыми отдельными научными группами».

Область знания, основной код классификатора 01 - Математика, информатика и науки о системах, 01-310 - Газовая динамика

Ключевые словаДетонация, устойчивость, слабонелинейные волны, асимптотический анализ, редуцированные модели, нелинейная динамика, разложение по динамическим модам, машинное обучение, нейронные сети.

Код ГРНТИ30.51.23


СтатусЗакрыт досрочно


 

ИНФОРМАЦИЯ ИЗ ЗАЯВКИ


Аннотация
Проект посвящен теоретическому и численному анализу многомерных уравнений Навье-Стокса для сжимаемого и реагирующего газа с целью описания динамики детонационных волн. Предполагается получение упрощенной асимптотической модели для слабонелинейных волн, в которой будут учитываться как эффекты сжимаемости, так и вязкости, диффузии и теплопроводности. Будет проведен анализ возможных автомодельных решений асимптотической модели, их линейной устойчивости и нелинейной динамики. С помощью численных расчетов решений асимптотической модели будет проведено исследование перехода от быстрого дефлаграционного режима к детонационному режиму в слабонелинейном приближении. Предложенный анализ имеет фундаментальное значение для разработки теории сложной нелинейной динамики волновых процессов в химически реагирующих потоках и для понимания механизмов перехода дефлаграции в детонацию. Также будет разработан подход, основанный на методах машинного обучения с целью построения нейронной сети для решения задач динамики детонации. Успешная реализации подобного подхода позволит многократно ускорить расчеты детонации.

Ожидаемые результаты
Примечание: здесь и далее даются ссылки на уравнения (xx) и публикации [yy], приведённые в приложенных к проекту дополнительных файлах. Математическая модель (31-33) будет исследована как в стационарном, так и в нестационарном случаях. Будет проведен анализ существования и численный расчет решений в виде бегущих волн на основе геометрической теории сингулярных возмущений [16, 18, 19]. Первые численные расчеты бегущих волн уже проведены в [15], но требуется их более детальный анализ, в особенности в том, что касается строгого математического обоснования существования решений, а также их эффективного численного расчета. Ввиду сильно сингулярного характера уравнений численный счет связан с определенными трудностями. Например, в случае расчета перехода от дефлаграции к детонации, процесс сопровождается явлениями с существенно разными временными масштабами. Фронт дефлаграции распространяется со скоростью, которая, в основном, определяется диффузионными процессами, в то время как ускорение этого фронта вызывает волны сжатия, которые распространяются со значительно большими скоростями. Численный метод должен быть в состоянии рассчитать как и тот, так и другой процесс. Потребуется разработка соответствующей численной схемы и её эффективная имплементация. Будет сделан строгий вывод асимптотической модели путем анализа уравнений Навье-Стокса для сжимаемых реагирующих течений с целью заменить или обосновать качественную модель (31-33). При этом аналитические и численные исследования как модели (31-33), так и полной системы уравнений Навье-Стокса будут способствовать выявлению основных факторов, которые играют важную роль в динамике волн быстрого горения и детонации. Будет проведен анализ существования и численный расчет решений в виде бегущих волн на основе геометрической теории сингулярных возмущений. Будет разработан новый алгоритм машинного обучения и построена нейронная сеть для анализа как прямой, так и обратной задач детонации. Такой подход ранее не применялся к решению задач детонации; его успешная реализация может оказать прорывное влияние на расчеты динамики многомерной детонации. Асимптотический анализ, проведенный совместно с численными расчетами, позволит прояснить физические механизмы, ответственные за сложную нелинейную динамику детонации и, в особенности, за процесс перехода быстрого горения в детонацию. Это будет сделано на основе упрощенной модели, полученной строго из полных уравнений движения. Применение методов машинного обучения и разработка нейронной сети для расчетов детонационных волн может дать в перспективе инструменты для намного более быстрых расчетов детонационных волн, чем можно делать в настоящее время. Как хорошо известно, волны детонации являются исключительно требовательными с точки зрения вычислительных ресурсов. Современные методы машинного обучения успешно применяются для решения многих сложных задач механики и физики и их адаптация для задач данного проекта представляется перспективной в области расчетов детонационных волн и вообще в вычислительной газодинамике.


 

ОТЧЁТНЫЕ МАТЕРИАЛЫ


Аннотация результатов, полученных в 2022 году
Проведен численный анализ новой математической модели, описывающей процессы дефлаграции и детонации в горючих газовых смесях. Особенностью рассмотренной модели является одновременное присутствие как диссипативных эффектов (диффузии, теплопроводности и вязкости), так и эффектов слабой нелинейности и сжимаемости, что позволяет исследовать в рамках одной модели и явление детонации с учетом диссипативных эффектов, и явление дефлаграции с учетом эффектов сжимаемости. Для численного решения модели используется метод расщепления по физическим процессам, что позволяет рассмотреть сложные явления, происходящие на разных масштабах в пространстве и времени. Результаты расчетов для начальной задачи инициирования указывают на соответствие решений модели реальным процессам в волнах дозвукового и сверхзвукового горения, а также на способность модели описывать переходные режимы, важные с теоретической и практической точек зрения. Рассмотрена задача об инициировании в трубе от горячего источника, расположенного у закрытого конца трубы. Обнаружено, что при достаточно низкой температуре инициатора образуются быстрая акустическая волна-прекурсор и медленная волна горения, которая распространяется по среде, возмущенной прекурсором. Если же температура инициатора достаточно высока, то волна горения может резко ускориться и перейти в детонационный режим. Таким образом, явление перехода медленного горения в детонацию впервые описывается в рамках упрощённой математической модели. Нами также проведен асимптотический анализ полной системы уравнений Навье-Стокса для сжимаемой реагирующей среды с тем, чтобы получить полностью рационально обоснованную упрощенную модель перехода горения в детонацию. Предсказательные свойства упрощенных моделей будут изучены на основе численного кода для решения полной системы уравнений Навье-Стокса для сжимаемых реагирующих течений, который разработан нами в среде OpenFOAM.

 

Публикации


Аннотация результатов, полученных в 2023 году
Проведен численный и теоретический анализ новой упрощенной математической модели, описывающей процессы дефлаграции и детонации в горючих газовых смесях. Особенностью и преимуществом модели является её простота по сравнению с полной системой уравнений Навье-Стокса для сжимаемых реагирующих течений и, в то же время, учет в ней как диссипативных эффектов (диффузии, теплопроводности и вязкости), так и эффектов слабой нелинейности и сжимаемости. Эти свойства позволяют исследовать в рамках данной относительно простой модели как явления детонации с учетом диссипативных эффектов, явления дефлаграции с учетом эффектов сжимаемости, а также - что представляет особый интерес - явление перехода горения в детонацию (ПГД). Для численного решения уравнений модели используется метод расщепления по физическим процессам, что позволяет рассмотреть сложные явления, происходящие на разных масштабах в пространстве и времени. Результаты расчетов для начальной задачи инициирования указывают на соответствие решений модели реальным процессам в волнах дозвукового и сверхзвукового горения, а также на способность модели описывать переходные режимы, важные с теоретической и практической точек зрения. Численные эксперименты указывают на наличие ПГД в достаточно длинных трубах с закрытым левым концом и подавление ПГД в случае открытого конца. В качестве физического механизма ПГД выявлено образование слабых пульсирующих волн давления, образующихся у закрытого горячего конца трубы, которые догоняют зону реакции и фронт пламени, дополнительно прогревая смесь перед пламенем за счет сжатия. Это приводит к ускорению пламени и спонтанному переходу к детонации при достаточно интенсивном прогреве смеси. Также нами исследованы аналитически стационарные решения указанной модельной системы и проведено качественное сравнение с поведением решений полной системы уравнений Навье–Стокса. В рамках проведенной работы разработан алгоритм для численного построения бегущих волн модельной системы. Соответствующий программный комплекс реализован на языке C++ с использованием библиотеки конечных элементов с открытым исходным кодом deal.II. Проведен асимптотический анализ для построения решений системы в приближении слабой диссипации и получены аналитические решения. Классифицированы и построены четыре типа стационарных решений, с описанием их свойств, характерных масштабов и ограничений на параметры системы. Проведено сравнение этих решений с результатами, полученными численно, а также проведено качественное сравнение нестационарных решений модельной системы и системы Навье–Стокса. Для численного решения уравнений Навье-Стокса для реагирующих течений разработана программа в OpenFOAM. Также, впервые для задач детонации применен один из самых перспективных методов машинного обучения для решения уравнений в частных производных, а именно метод нейронного оператора Фурье. На примере модельного уравнения из теории детонации продемонстрированы основные свойства данного метода и его преимущества. Для двумерных волн детонации, описываемых системой уравнений Эйлера, распространяющихся в среде с поперечными периодическими неоднородностями нами впервые продемострирована возможность синхронизации и захвата мод детонационных ячеек. Показано, что относительно слабые возмущения в среде приводят к перестройке ячеистой структуры с ее возможной регуляризацией, а сильные возмущения способствуют образованию стационарных течений разной степени стабильности.

 

Публикации

1. Голдин А.Ю., Касимов А.Р. Detonation cell locking in periodically nonuniform mixtures Proceedings of the Combustion Institute, - (год публикации - )

2. Голдин А.Ю., Магомедов Ш.М., Фариа Л.М., Касимов А.Р. Study of a qualitative model for combustion waves: flames, detonations, and deflagration-to-detonation transition Computers and Fluids, - (год публикации - )

3. Шахзеб Аамир, Аслан Касимов Application of the Fourier Neural Operator to the reactive Burgers equation of detonation theory Mathematics, - (год публикации - )


Возможность практического использования результатов
Идеи регуляризации и стабилизации детонационных волн через создание специальных условий в исходной смеси или в канале имеют непосредственное отношение к технологиям создания детонационных двигателей. Эти идеи были выдвинуты впервые руководителем проекта и активно развиваются, начиная с пионерской работы, представленной в 2019 году на конференции в Китае (Kasimov & Gonchar, ICDERS, 2019). Для практического применения требуется развитие этих идей и исследование детонационных процессов в многомерной постановке. Это планируется сделать в дальнейшем на основе тех инструментов, включая программу в OpenFOAM для трехмерных расчетов, которые были разработаны в рамках данного проекта.