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

Numerical simulation of dispersed mixtures combustion incorporating condensed phase of multiple fractions

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

грант РФФИ

Этапы НИР

# Сроки Название
1 1 февраля 2018 г.-31 декабря 2018 г. Вычислительное моделирование горения диспергированных топлив с многофракционной конденсированной фазой (этап 1, 2018 год)
Результаты этапа: За первый год выполнения проекта: 1. Были разработаны замкнутые вычислительные модели и отлажены компьютерные коды для расчета турбулентного горения в потоке газообразного окислителя капель жидкого углеводородного горючего с содержащим ися в них дисперсными частицами ме таллов с использованием редуцированных моделей макрокинетики горения и зажигания, учитывающей физические и химические превращения и турбулентное перемешивание. 2. Была разработана модель зажигания струй в камере. 3. Была проведена отладка и тестирование компьютерных кодов. В проекте разработана математическая модель горения капли жидкости с металлическими включениями. Разработаны математические модели горения одиночной капли с включением мелких частиц металла в случае смачивающей жидкости с учетом повышенной теплопроводности металла и существенно более высокой температуры его воспламенения нежели температура испарения жидкости. Разработанная модель предполагает, что при попадании такой комбинированной капли в камеру сгорания вначале происходит испарение жидкого горючего, воспламенение его паров и горение в газофазном режиме с образованием диффузионного микропламени. При этом температура поверхности остается на уровне температуры фазового перехода жидкого горючего, что не позволяет воспламениться металлическим частицам. На этом этапе роль частиц сказывается в изменении средней теплопроводности и теплоемкости комбинированной капли. По мере выгорания жидкого горючего концентрация твердых частиц вблизи поверхности возрастает. При достижении максимальной поверхностной плотности частиц их прогрев приводит к интенсификации поверхностного испарения жидкости, которая отступает вглубь. Агломерации частиц формируют на поверхности пористое тело, температура которого начинает повышаться, что приводит к воспламенению и формированию гетерогенной реакции. Выделяющееся тепло передается внутрь капли, что способствует дальнейшему испарению жидкости, выходу паров через пористое тело и поддержанию, тем самым, газофазной реакции. Окислитель при этом расходуется на поддержание как газофазной, так и гетерогенной реакций. Также в окислении металла может участвовать и углекислый газ, который является продуктом газофазной реакции. Оригинальность результатов состоит в разработке новых моделей диспергирования, нагрева и горения гетерогенных горючих, содержащих смесь жидких углеводородов и конденсированных мелкодисперсных порошков металлов. При этом учтены особенности воспламенения и горения таких систем, а именно: воспламенение и горение углеводородов в газовой фазе после испарения, преобладание поверхностных гетерогенных реакций при горении металлов, газообразные продукты горения углеводородов, конденсированные продукты горения металлов. Принципиально новым моментом в подходе к изучению зажигания и, в частности, его пределов, является учет всех видов неравновесности: механической (взаимодействие с газом и фрагментация частиц), фазового перехода, химической кинетики, состава капель, включающих жидкую и конденсировнную фазу.
2 1 января 2019 г.-31 декабря 2019 г. Вычислительное моделирование горения диспергированных топлив с многофракционной конденсированной фазой (этап 2, 2019 год)
Результаты этапа: В 2019 году было сделано: 1. Разработаны базовые модели для введения дополнений в код по учету многофазности и полидисперсности смеси в камере сгорания и особенностей прогрева и горения многофазных частиц, содержащих жидкие и конденсированные горючие компоненты: 1. Разработана модель с дополнительным ансамблем частиц на фоне поля течения газовой фазы. 2. Разработаны модели с учетом в газодинамических уравнениях источниковых членов в результате влияния ансамбля частиц при размере ячеек сетки больше размеров частиц. 3. Разработана модель работы форсунки при заданной функции распределения частиц по размерам и заданном расходе в условиях наличия обратных токов. 4. Создана физическая и математическая модели горения агломератов капель, и изучено влияние многофазности топлива (жидкий углеводород с угольными включениями) на условия зажигания и распространение горения в неравномерной полидисперсной среде. Учтено влияние скученности капель на скорость испарения каждой из них. Проведено компьютерное моделирование распространения пламени при горении полидисперсной распыленного капельно-жидкого топлива с твердыми угольными включениями. 5. Получено, что при наличии в капле твердого конденсированного топлива испарение жидкости происходит быстрее из-за дополнительного выделения энергии на поверхности за счет гетерогенных реакций. Когда от капли остается лишь твердый остаток, скорость массообмена резко уменьшается, но затем снова возрастает за счет повышения температуры остатка и соответствующего ускорения гетерогенных реакций. 6. Показано, что при распространении пламени в полидисперсной смеси толщина зоны горения может составлять десятки сантиметров. При распространении горения из центра закрытой камеры вначале толщина зоны горения остается постоянной. В дальнейшем по достижении фронтом пламени стенок камеры эта толщина увеличивается, форма же становится неустойчивой вследствие обратных потоков в сторону центра камеры. При этом общая интенсивность горения увеличивается, но максимум интенсивности снижается. 7. Проведенные аналитические исследования, а также имеющиеся экспериментальные данные, показали, что существуют два основных режима горения капель углеводорода: высокотемпературный и низкотемпературный. Причем полученные аналитические решения показывают, что существуют две поверхности пламени, причем первая поверхность пламени может быть зафиксирована измерительными приборами, а вторая поверхность пламени невидима. Для описания поведения первой и второй поверхностей пламени были получены аппроксимационные аналитические формулы.
3 1 февраля 2020 г.-31 декабря 2020 г. Вычислительное моделирование горения диспергированных топлив с многофракционной конденсированной фазой (этап 3, 2020 год)
Результаты этапа: При выполнении проекта были разработаны физическая, математическая и численная модели горения капли с включением частицы металла в случае смачивающей жидкости. С помощью полученных моделей были созданы компьютерные коды для расчета турбулентного горения в потоке газообразного окислителя капель жидкого углеводородного горючего с содержащимися в них дисперсными частицами металлов. Кинетика процесса горения и зажигания моделировалась с помощью редуцированных механизмов, учитывающих физические и химические превращения и турбулентное перемешивание. Вычислительный алгоритм, использующийся в программе, имеет опцию параллельной реализации. В 2020 году были получены следующие результаты: 1. Разработан вычислительный алгоритм параллельной реализации задачи горения многофракционных частиц в турбулентном потоке. 2. Изучена немонотонная зависимость ускорения вычислений при подключении дополнительных параллельных процессов. Расчеты производились на компьютере АПК-1 при использовании от 1 до 48 процессоров для различного числа нитей исполнения при методе OpenMP. Был выявлен немонотонный характер зависимости времени выполнения расчетов от числа нитей при числе процессов больше 40, то есть близкой к предельно возможному числу 48 для одной платы памяти. Причиной немонотонности может быть как характер пересылок, так и вмешательство конкурирующих системных процессов, замедляющих расчет при дефиците свободных исполнителей. Также было получено, что ускорение расчета возрастает линейно до 33 процессоров, а затем линейность нарушается. Максимум достигается на 43 и 47 процессорах. 3. Определено оптимальное количество параллельных процессов, при котором будет достигнуто максимальное ускорение работы алгоритма. Было показано, что эффективность на первых 6 процессорах быстро падает от 1 до 0.8. В целом после 33 процессоров эффективность в среднем уменьшается, но все же остается около 0.65 для максимального (для одной платы памяти) числа процессоров 48.

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

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