Отправьте статью сегодня! Журнал выйдет 12 июля, печатный экземпляр отправим 16 июля
Опубликовать статью

Молодой учёный

Математическая модель горения пропан-бутановой смеси при недостатке окислителя в диффузионной горелке

14. Общие вопросы технических наук
1682
Поделиться
Библиографическое описание
Польшиков, В. Ю. Математическая модель горения пропан-бутановой смеси при недостатке окислителя в диффузионной горелке / В. Ю. Польшиков. — Текст : непосредственный // Технические науки: теория и практика : материалы I Междунар. науч. конф. (г. Чита, апрель 2012 г.). — Чита : Издательство Молодой ученый, 2012. — С. 130-133. — URL: https://moluch.ru/conf/tech/archive/7/2118/.

Для идентификации процессов протекающих при сжигании топлива в горелочных устройствах необходимо спрогнозировать ход и параметры гидродинамических и тепло-массообменных процессов. Моделирование течения с горением является наиболее сложной задачей движения газа. Данная проблема может эффективно решаться в системе вычислительной гидродинамики FlowVision. FlowVision решает задачи внешней и внутренней аэро-гидродинамики, в частности, задачи диффузионного горения и горения предварительно подготовленных топливных смесей. Последнее обстоятельство является основным при выборе этой системы для моделирования процесса сжигания топливных компонент в диффузионной горелке с «реальными» размерами и расходами газов.

Виртуальный эксперимент в системе вычислительной аэро-гидродинамики, в отличие от натурных испытаний, значительно снижает трудоемкость подготовительных мероприятий, упрощает экспериментальную часть работы, более наглядно и детально раскрывает суть происходящих процессов.

Модель горения в системе FlowVision включает расчет конечно-объемным методом параметров процесса по уравнениям:

- Навье-Стокса (закон сохранения импульса);

- неразрывности (закон сохранения массы сплошной среды);

- закона сохранения энергии, записанного через энтальпию;

- переноса реагентов в приближении брутто-реакции (перенос топлива, окислителя и продуктов сгорания);

- k-&#; модели турбулентности.

Эта модель представляет собой развитие модели течения слабосжимаемой жидкости со всеми ее ограничениями, но при этом учитывается горение газовой смеси. Предполагается, что газовая смесь может быть предварительно перемешана и не перемешана. Ограничения модели предполагают течение с любыми изменениями плотности среды (однако числа Маха все равно малы – т.е. изменения плотности обусловлены только температурными эффектами или влиянием примесей в газе).

При решении задач в FlowVision дифференциальные уравнения аппроксимируются на расчётной сетке в предположении, что каждая ячейка представляет собой конечный объём, в котором скорости изменения физических величин сбалансированы потоками этих величин через грани ячейки. Для расчета используется прямоугольная сетка с локальным измельчением расчетных ячеек у криволинейных границ, кромок и граничных областей малого размера.

Решение конкретных задач во FlowVision предполагает создание трехмерной модели реального объекта во внешней системе автоматизированного проектирования (AutoCAD, ParaSolid, SolidWorks, T-Flex и др.) с последующим импортом геометрии в препроцессор FlowVision. В препроцессоре реализуется интерактивное задание граничных условий на поверхностях модели, диагностика и трансформация геометрии. Далее следуют этапы:

- выбор математической модели (в нашем случае – Модель горения);

- задание граничных условий;

- задание исходной расчетной сетки и критериев её адаптации по решению и по граничным условиям;

- задание параметров расчетного метода;

- проведение расчета (без участия пользователя);

- просмотр результатов расчета по выбранным параметрам в графической форме;

- оценка точности расчетов методом сходимости по сетке.

Блок расчета уравнений и постпроцессор (часть программы, в которой анализируются результаты расчета) объединены и работают одновременно. Это позволяет проводить моделирование и одновременно анализировать результаты, менять граничные условия и параметры математической модели.

Для моделирования процесса в диффузионной горелке с кольцевым соплом для подачи воздуха построена ее 3D-модель в системе T-Flex CAD 3D (рис. 1).

Рисунок 1 – Разрез 3D-модели диффузионной горелки



Геометрия из T-Flex CAD экспортировалась через обменный формат Stereo Lithography (stl) и подгружалась в систему FlowVision. В препроцессоре выбиралась осесимметричная задача, устанавливались начальные параметры соответствующие нормальным условиям, из стандартной библиотеки загружались свойства компонентов, подаваемых на горение (воздух и пропан-бутан), и свойства продуктов сгорания для заданной смеси, создавался фильтр одноразового зажигания в параллелепипеде и задавались граничные условия (названия границ указаны стрелками на рис. 2).

В качестве начальных значений и физических параметров модели задавались начальная температура T=293 К; давление 101325 Па; пульсация 0,03; масштаб турбулентности 0,01 м, стехиометрический коэффициент km=15,67 (для топливной смеси пропан-бутан + воздух); температура воспламенения 750 К; диапазон коэффициента избытка окислителя в котором горение поддерживается самопроизвольно &#&#61537;;</FONT>=0,02 &ndash; 3, выбиралась пульсационная модель химической кинетики процесса.</P> <P><A HREF="images/m467f5b5d.gif" TARGET="_blank"><IMG src="https://articles-static-cdn.moluch.orgimages/m467f5b5d.gif" NAME="Объект1" ALIGN=BOTTOM WIDTH=416 HEIGHT=323></A></P> <P><B>Рисунок 2 &ndash; Граничные условия при моделировании процесса в диффузионной горелке</B></P> <P><BR /><BR /> </P> <P>Использовались следующие параметры для границ и источника зажигания в примере расчета форсированных режимов горения (с высоким выходом углеродного депозита при <FONT FACE="Symbol">&#&#61537;;</FONT>=0,4; <SPAN LANG="en-US"><I>k</I></SPAN><SUB><SPAN LANG="en-US"><I>m</I></SPAN></SUB>=6,9):</P> <P>- &laquo;Стенка&raquo; &ndash; тип границы &ndash; &laquo;Стенка&raquo;, логарифмический закон изменения скорости;</P> <P>- &laquo;Вход газа&raquo; &ndash; тип границы &ndash; &laquo;Вход/Выход&raquo;, температура 20 <SUP>о</SUP>С, концентрация 1, нормальная массовая скорость 0,1139<FONT FACE="Symbol">&#&#61655;;</FONT>10<SUP>-4</SUP> кг/(м<SUP>2</SUP><FONT FACE="Symbol">&#&#61655;;</FONT>с);</P> <P>- &laquo;Вход воздуха&raquo; &ndash; тип границы &ndash; &laquo;Вход/Выход&raquo;; температура 20 <SUP>о</SUP>С, концентрация 0, нормальная массовая скорость 0,10834 кг/(м<SUP>2</SUP><FONT FACE="Symbol">&#&#61655;;</FONT>с);</P> <P>- &laquo;Выход&raquo; &ndash; тип границы &ndash; &laquo;Свободный выход&raquo;;</P> <P>- &laquo;Параллелепипед зажигания&raquo; &ndash; для расчета &laquo;холодного&raquo; течения &ndash; &laquo;Неактивный&raquo;, для расчета горения &ndash; &laquo;Одноразовый&raquo;.</P> <P>Для расчетной подобласти генерировалась начальная прямоугольная сетка с разбиением на 50 участков по координатным осям.</P> <P>Процесс моделирования включал два этапа. На первом рассчитывалось &laquo;холодное&raquo; течение газовых компонент без горения. При этом обеспечивалось некоторое распределение коэффициента избытка окислителя <FONT FACE="Symbol">&#&#61537;;</FONT> внутри горелки, при котором горючее достигает границы &laquo;Выход&raquo;.</P> <P>После того как было получено распределение коэффициента избытка окислителя, проводилась инициация процесса горения в &laquo;параллелепипеде зажигания&raquo; с использованием подготовленного фильтра зажигания (переключение режима &laquo;Неактивный&raquo; в режим &laquo;Одноразовый&raquo;).</P> <P>Расчет проводился до стационарного состояния, которое контролировалось по погрешности основных расчетных переменных (точность расчета 10<SUP>-4</SUP>).</P> <P>Для диффузионной горелки с заданным расходом горючего и окислителя получены распределения коэффициента избытка окислителя (рис. 3), скоростей (рис. 4) и температуры (рис. 5).</P> <P><A HREF="images/m1595c04b.jpg" TARGET="_blank"><IMG src="https://articles-static-cdn.moluch.orgimages/m1595c04b.jpg" NAME="Рисунок 3" ALIGN=BOTTOM WIDTH=340 HEIGHT=268 BORDER=0></A></P> <P><B>Рисунок 3 &ndash; Распределение величины коэффициента избытка окислителя </B><B><FONT FACE="Symbol">&#&#61537;;</FONT></B></P> <P><B>в сечении диффузионной горелки</B></P> <P><BR /><BR /> </P> <P><A HREF="images/m7751eddc.jpg" TARGET="_blank"><IMG src="https://articles-static-cdn.moluch.orgimages/m7751eddc.jpg" NAME="Рисунок 2" ALIGN=BOTTOM WIDTH=326 HEIGHT=263 BORDER=0></A></P> <P><B>Рисунок 4 &ndash; Векторное поле модуля скорости в сечении диффузионной горелки</B></P> <P><BR /><BR /> </P> <P><A HREF="images/m36e60025.jpg" TARGET="_blank"><IMG src="https://articles-static-cdn.moluch.orgimages/m36e60025.jpg" NAME="Рисунок 1" ALIGN=BOTTOM WIDTH=350 HEIGHT=264 BORDER=0></A></P> <P><B>Рисунок 5 &ndash; Температурное поле в сечении диффузионной горелки</B></P> <P><BR /><BR /> </P> <P>Анализ результатов моделирования показывает, что в горелке при раздельной подаче горючего и окислителя с заданными параметрами наблюдается неоднородное смесеобразование с коэффициентом избытка окислителя близким к стехиометрическому вблизи границы расширяющегося &laquo;конуса&raquo;, образованного истекающим из центрального сопла горючим. Область богатой топливной смеси, в которой может образовываться углерод в конденсированной фазе, сосредоточена в конусе под центральной частью подложки-саженакопителя.</P> <P>На векторном поле скоростей видно, что по периферии наблюдается циркуляция потока продуктов сгорания (не происходит проскок компонентов топливной смеси без реагирования). Т.е в установившемся режиме можно эффективно регулировать соотношение окислитель/горючее, варьируя их расход на входе в горелку.</P> <P>На подложку-саженакопитель при заданных расходах поступают продукты сгорания с разной температурой.</P> <P>Следует отметить, что при более мягких режимах горения с меньшим расходом компонентов пламя имеет фактически цилиндрическую форму от среза сопла и расширяется на конус малой высоты только вблизи подложки-саженакопителя. При этом диаметр конуса в области сечения подложки-саженакопителя зависит от расходов газовых компонент.</P> <P><BR /><BR /> </P> <P>Литература:</P> <P>1. Польшиков, В.Ю. Разработка технологии и оборудования для синтеза углеродных наноструктурных материалов в диффузионном пламени / В.Ю. Польшиков // Проблемы техногенной безопасности и устойчивого развития: сборник научных статей молодых ученых, аспирантов и студентов. Тамбов: изд-во ГОУ ВПО ТГТУ. &ndash; Вып. II. 2011 &ndash; С .202-206.<P> 2. Борунова, А.Б. Получение углеродных наночастиц при горении метана / А.Б. Борунова, Ю.В. Григорьев, К.Я. Трошин // Горение и взрыв. &ndash; М.: ТОРУС ПРЕСС, 2008. &ndash; С. 10-13.

Можно быстро и просто опубликовать свою научную статью в журнале «Молодой Ученый». Сразу предоставляем препринт и справку о публикации.
Опубликовать статью

Молодой учёный