Математическая модель определения начального участка гетерогенного факела и ее адаптация

УДК 662.95

DOI: 10.14529/power160302

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ОПРЕДЕЛЕНИЯ НАЧАЛЬНОГО УЧАСТКА ГЕТЕРОГЕННОГО ФАКЕЛА И ЕЕ АДАПТАЦИЯ

Е.В. Торопов, К.В. Осинцев

Южно-Уральский государственный университет, г. Челябинск

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

Ключевые слова: факел, начальный участок, воспламенение, математическая модель.

Введение

В теории горения энергетических топлив общепризнанными являются математические модели с критическими условиями самовоспламенения в нестационарной и стационарной постановке, связанные с работами советских академиков Я.Б. Зельдовича и Н.Н. Семенова [1, 2]. Эти модели построены на анализе основных соотношений тепло-переноса в объеме предварительно перемешанных компонентов горения при различных условиях адаптации к реальным схемам. Причем нестационарная теория рассматривает распространение пламени в закрытом объеме, что соответствует в основном нештатным ситуациям в энергетике, а стационарную теорию можно применить к устойчивому процессу горения в движущейся смеси компонентов горения. Обе теории содержат рекомендации по управлению процессом горения при воспламенении от внешнего источника. Изложенный ниже материал дополняет ряд данных по математическому описанию различных элементов процесса горения топлива, начиная с участка, где смесь нагревается до условий воспламенения.

1. Постановка задачи.

Процессы в начальном участке факела

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

смеси в пределах «горючего состава» и наличии локальных актов горения с тепловым эффектом, достаточным для инициирования и распространения горения на весь объем смеси. С точки зрения технологии горения это означает, что смесь находится в интервале воспламенения, между «нижним» и «верхним» пределами, скорость потока сопоставима со скоростью распространения пламени, энтальпия горения достаточна для поддержания процесса.

При классификации способов воспламенения горючих смесей с учетом классических методов можно условно выделить три схемы:

- нагрев объема смеси до температуры самовоспламенения 4в;

- зажигание от контакта с пламенем, имеющим температуру 4л > tсв, но при этом температура смеси может быть ниже 4В;

- зажигание от нагретой твердой поверхности при 4ов > 4в.

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

При определении зависимости температуры факела от его длины Тф(Ьф) посредством предложенной ранее интегральной функции тепловыделения Ф^ф) начало горящего гетерогенного факела размещается в плоскость воспламенения Lф = 0, но эта плоскость отстоит от выходного сопла го-релочного устройства на расстоянии 1н, м. Таким образом, для согласования геометрических размеров факела с функцией тепловыделения необходимо определить протяженность начального участка и сформулировать систему физико-математических соотношений для определения факторов, влияющих на 1н, что и является задачей настоящей работы.

Процесс термического разложения вещества твердого топлива имеет сложную природу, зависящую от исходного состояния топлива и условий нагрева; этот процесс можно разделить на три стадии - бертинирование (0-300 °С), полукоксование (300-600 °С) и коксование (600-1100 °С). На рис. 1 изображены термо-гравиметрические кривые, построенные по данным [3, 4], где показано изменение выхода продуктов термической деструкции природного твердого топлива, относящегося к категории гумолитов, в процентах на сухую беззольную массу.

В первой стадии выделяется небольшое количество газообразных продуктов разложения ~ 3,0 % от исходной массы, с низкой теплотой сгорания. При определении начального участка факела эту стадию можно не учитывать, тем более, что продукты газификации в первой стадии учитываются в составе второй стадии. Температурные функции выделения трех компонентов летучих

веществ - газовой фазы КгТаз, пирогенетической воды Кггаз и смолистых веществ на горючую

(сухую беззольную) массу можно описать следующими аппроксимациями:

Кггаз = 3,0 + 6,17-10"2М ;

кпв = 3,0 +1,5-10"2М ;

Км = 3,7-10"2М.

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

Применение концепции факельного континуума [3] позволяет рассматривать задачу в двух измерениях - с учетом продольной I и поперечной г координат. Причем теплофизические свойства факельной среды считаются постоянными в направлении г, ступенчато изменяющиеся к свойствам топочного газа на границе факела. Задача для начального участка включает не только определение его длины 1н, но и определение возможности применения управляющих воздействий в направлении изменения основных его характеристик.

При определении начального участка гетерогенного факела существенную роль играет динамика выхода летучих веществ из топливных частиц при нагреве. Суммарный эффект выхода летучих веществ определяется сложением уравнений (1)-(3), а перевод в измерение в процентах от рабочей смеси - умножением на (100 - Ар - ^)/100, в результате получается уравнение

Vр = 3,480 + 6,595 -10"2 М

или в массовых долях

Vр = 0,0348 + 0,6595 • 10"3 М .

100

>>

о о cd

2 г 2 Он

о

80-

60-

40-

§ 20

Й

PQ

1

200 400 600

Температура, оС

800

1000

Рис. 1. Выход продуктов термического разложения в зависимости от температуры: 1 - газ; 2 - пирогенетическая вода; 3 - смола; 4 - твердый остаток

0

Уравнения (4)-(5) связывают количество выделяемых летучих веществ с уровнем нагрева частиц топлива в интервале температур 300-600 °С, причем Дt = t - 300. При получении уравнений (4)-(5) приняты средние значения балластных составляющих бурого угля = 7,0 %, Ц^ = 35,0 %.

При выборе алгоритма определения начального участка факела от выхода дисперсного потока из горелочного устройства до начала воспламенения 1н необходимо выделить несколько схем по принципу учета влияния определяющих факторов, на основе идентификации которых разрабатываются математические модели расчета 1н. Предлагается следующая схема классификации:

- расчет на основе определения температуры интенсивного выхода летучих веществ tштт:, которая косвенно соответствует температуре воспламенения;

- расчет на основе определения концентрации летучих веществ в смеси с воздухом ц^, которая формируется динамикой нагрева топливных частиц и соответствует нижнему пределу воспламенения газо-воздушной смеси цсм =

- применение стандартной методики расчета взрывобезопасных концентраций топливной пыли в системе подготовки и транспорта топливной пыли [3];

- использование экспериментальных данных по времени воспламенения летучих веществ в неподвижном окислителе [4];

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

Все перечисленные алгоритмы относятся к

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

2. Математическая модель определения 1н на основе расчета

Для учета динамики нагрева частиц топлива размером 5, м, можно применить критериальное уравнение для термически тонких тел

1п © = -3Bi • Fo,

где Ы = а5А. - число Био, Бо = ах/5 - число Фурье, 0 = 0л/00 - безразмерная температура топливной частицы. На рис. 2 представлена номограмма, построенная по уравнению (6), которая позволяет быстро определить безразмерное время процесса нагрева Бо для частиц размером 5, причем область больших значений числа Бо > 10 относится к частицам меньшего размера 5.

Зависимость (6) преобразуется с учетом коэффициентов массивности т = 1,003-1,067 для самых мелких и самых крупных частиц топливной пыли и формы £ = 3,0:

т = [5рст/

где 00 = ^ - 4ач - начальный температурный перепад; 0л = ^ - 4 - температурный перепад соответствующий выделению летучих веществ в количестве, достаточном для воспламенения при нижнем

1п6>

Число Го

Рис. 2. Номограмма по определению безразмерного времени процесса нагрева Ро для частиц размером 5 по числу В^ 1 - Б1 = 0,18; 2 - Б1 = 0,091; 3 - Bi = 0,045; 4 - Bi = 0,0364; 5 - Bi = 0,0182; 6 - Б1 = 0,01164; 7 - Б1 = 0,0091; 8 - Bi = 0,0082; 9 - Bi = 0,0045

пределе цн; также 4, 4ач, 4 - температура топочных газов, начальная температура топливных частиц и температура воспламенения летучих веществ соответственно, °С. Температуре 4 присваивается смысл такого нагрева топливной пыли, который приводит к выделению летучих веществ в количестве, достаточном для воспламенения.

Применение рис. 2 ограничено условиями внешнего нагрева частиц пыли с интенсивностью результирующей теплоотдачи а, Вт/м2 К. После воспламенения результирующий тепловой поток определяется с учетом собственного излучения частицы.

Выделение летучих веществ и их смешение с воздухом в количестве Уа = авК0, м3/кг топлива, приводит к образованию газо-воздушной смеси, состав которой можно определить в соответствии с (5) в массовых долях по формуле

^см = Vр/

Этой концентрации газо-воздушной смеси соответствует динамика выделения летучих веществ, определяемая с учетом (5) по формуле

^см = т^-3-;

(0,0348 + 0,6595 • 10-3 М + рв Vа ) при получении обратной функции

М =цсм

0,6595•Ю-3 (1 — цсм) "

При достижении газо-воздушной смесью концентрации, соответствующей нижнему пределу воспламенения цсм = цвс, температура соответствует температурному перепаду воспламенения

0,5

ДtВc = tвс - 300, а 0л = 0вс = tт - Дtвс и с заменой концентрации и температуры на таковые при воспламенении можно в (7) подставить вместо 0л = 0вс, что дает возможность определить время воспламенения при т ~ 1,0 и £ = 3,0:

ТвС =[0,4 •Ю6 5/6а] 1п

где 5 - размер топливной частицы, м; а - коэффициент суммарной теплоотдачи на поверхность частицы топливной пыли, Вт/м2К, 00 = tт - tнaч.

С учетом (6), (11), а также аэромеханического коэффициента увеличения поперечного сечения факела ^ = (1 + 0,211н /й0)-2, можно получить

1н = Щ kF V

и далее

% +9,5238й0/12 +Ц3379Й3н = 56^0й02/а, (13) где принято: 4 = 1220 °С, 4ач = 20 °С, 4с = 650 °С.

Уравнение (13) можно решить методом Кар-дано или численно-графическим методом с помощью построения соответствующего графика. На рис. 3 изображен график, который позволяет определить /н в соответствии с уравнением (13), значения левой части которого нанесены на ось ординат в функции искомого /н, а значения правой части уравнения (13), которые не зависят от /н, соответствуют абсциссе.

График (см. рис. 3) построен для значения й0 = 0,20 м, если на практике определяется длина начального участка для горелочного устройства с диаметром выходного сопла йпр, отличным от й0, к определенному по графику значению /нгр вводится поправка, равная 4,375йф, то есть, /нф = 4,375йф/нгр. Если значение правой части уравнения (13) не совпадает с нанесенными на линии числами, то можно

~Г

0,3 Комплекс

График для определения начального участка факела 4 по значению комплекса кка = 2,2576 п>0 /а

0,4-s 0,3-

cd X

s

ц

4 0,2-од-

0

Рис. 3.

применить интерполяцию, или рассчитать индивидуальное значение числа правой части для конкретных условий. При определении 1н для тестовых значений w0 = 10 м/с, d0 = 0,20 м, а = 150 Вт/м2К получаем по графику 1н = 0,212 м, что удовлетворительно согласуется с экспериментальными данными.

3. Математическая модель определения 1н

на основе расчета

Температура начала выхода летучих веществ зависит от степени метаморфизма исходного угля и изменяется от 170 °С для бурых углей до 380 °С для антрацита [4]; завершается выход летучих веществ при температуре 800-1000 °С. Приняв, как и ранее, интервал 300-600 °С основным для определения состава газовой фазы и пределов воспламенения, можно отметить, согласно рис. 1, что процесс нагрева частиц угольной пыли в этом интервале температур аналогичен начальному периоду нагрева угольной шихты при получении кокса.

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

Пределы воспламенения смесей различных газов с окислителем по методике ЭНИН [5] можно определить, поделив сумму объемных концентраций индивидуальных газов на сумму отношений этих концентраций к пределам воспламенения индивидуальных газов. При разбавлении газовой смеси инертными примесями С02) расчет

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

Изменение пределов воспламенения смеси с инертными компонентами определяется по экспериментальным графикам ЭНИН [5]. Как показывают расчеты пределов воспламенения по методике ЭНИН, при переходе от промышленного коксового газа к факельному газу, который считается аналогичным первичному коксовому газу, нижний предел воспламенения дн в объемных процентах возрастает в 1,352 раза - от 5,70 до 7,71 %, а верхний предел воспламенения снижается в 0,816 раз -от 29,4 до 24,0 %. Перевод объемных процентов в массовые доли производится согласно простой зависимости = 7,71рг /(100рсм) = 0,07. При

этом значении дн по формуле (10) получаем Двс = 570 °С и 0вс = 1220 - 570 = 650 °С и соответствие с предыдущим расчетом, отличие заключается в определении температуры воспламенения через расчет нижнего предела воспламенения по данным ЭНИН.

4. Адаптация стандартной методики

определения взрывобезопасных

концентраций топливной пыли

В теплоэнергетике принято определять взры-воопасность пыли природных твердых топлив при наличии источника зажигания, которая зависит от состава топлив, химической активности компонентов летучих веществ, содержания кокса золы. Для оценки взрывоопасности природных твердых топ-лив используется критерий взрываемости Кт, определяемый расчетным путем как отношение объема летучих веществ на сухую массу топлива Vе к нижнему пределу воспламеняемости с учетом нелетучего остатка (кокса и золы) днлб /Кт = К7днлб. По величине Кт энергетические твердые топлива разделены на четыре группы от Кт < 1,0 до Кт > 3,5.

Методика определения нижнего предела воспламенения основана на расчете теплоты сгорания летучих веществ на сухую беззольную массу, которая также корректируется с учетом кислорода топлива и называется «кислородной массой». Теплота сгорания летучих веществ определяется по разности между теплотой сгорания исходного топлива и теплотой сгорания углерода кокса.

Нижний предел воспламенения определяется в двух вариантах: первый расчет достаточно обоснован, так как днл определяется делением экспериментальной константы 1,26105 на теплоту сгорания летучих веществ. Второй расчет производится с учетом твердого нелетучего остатка (кокса и золы), что вызывает сомнения, так как оставшийся после выхода летучих веществ твердый остаток почти не влияет на воспламенение смеси летучих веществ с воздухом.

Эти данные не предоставляют возможности использовать определенные в [6] Кт для расчета начального участка факела, так как в основе определения критерия взрываемости лежит не расчет нижнего предела воспламенения смеси летучих веществ с воздухом, а расчет днлб, что не соответствует физико-химической природе воспламенения.

Но применить результаты работы [6] можно в отношении цнл, но при этом необходимо произвести пересчет с сухой беззольной массы на рабочую массу. Такой расчет для тестовых условий (Наза-ровский бурый уголь марки Б2) дает 7,23 мас. %, что хорошо согласуется с данными, рассчитанным по методике ЭНИН. Подстановка значения днл = 7,23 рг /(100 рсм) в формулу (10) дает возможность определить по описанной выше методике и по графику рис. 3: 1н. = 0,21 м.

5. Определение 1н по экспериментальным данным для времени нагрева топливной пыли до воспламенения летучих веществ

В работе [7] опубликованы результаты экспериментального исследования характеристик широкой гаммы энергетических твердых топлив в отношении воспламенения и выгорания частиц топливной пыли. В методическом и физико-техническом плане опыты тщательно проработаны и исполнены с применением специально изготовленной аппаратуры. В части определения длины начального участка факела приводится обобщенная зависимость времени воспламенения летучих веществ тв.л в секундах как функция ряда факторов: температуры топочной среды Тг, К, размера частиц топлива 5, м, и опытного коэффициента пропорциональности &в.л, индивидуального для каждой марки угля [7]:

тв.л = 5,3 -1014 £в.лТг-480,8.

При расчете по формуле (14) для тестовых значений Назаровского угля марки Б2 при &в.л = 1,0 получаем твл = 0,04 с. При переходе к формуле (12) и учете аэромеханического коэффициента расширения факельной струи в виде ^ = (1 + 0,21/н ^0)-2 получаем 1н = 0,25 м. Отклонение от аналогичных данных в п. 2 составляет +15...+18 %, что легко объясняется тем фактом, что в экспериментальных данных [7] приводятся данные, полученные при относительной скорости потока окислителя около 0,1 м/с, в то время как в реальном факеле, особенно в начальном участке, относительная скорость значительно выше. Поэтому согласование вышеприведенной методики с практическими данными [7] следует признать удовлетворительным.

6. Определение 1н по экспериментальным данным для Питтсбургского угля [8]

В работе [8] приводятся данные экспериментального изучения показателей процесса горения угольной пыли из Питтсбургского битуминозного угля - по терминологии РФ из угля марки Ж. Опыты проводились на специальной установке, оборудованной устройством для отбора проб твердой и газообразной фаз. Содержание летучих веществ в исходном угле составляло 36 % на сухую массу, выводы акцентировались на частицы размером 200 мкм, несмотря на то, что основное число частиц располагалось в области 0-60 мкм.

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

движения от горелки по известному распределению скорости потока.

Содержание летучих веществ в массовой доле от начальных 100 % определялось по содержанию влаги, золы и расчетному углероду; на временной абсциссе также фиксировалось расположение фронта воспламенения твл = 0,051 с. Эти данные позволяют адаптировать результаты работы [8] к параметрам решаемой задачи.

В п. 5 была решена аналогичная задача для условий работы [7]; в эти данные необходимо внести изменения, связанные в марками углей: для марки Б2 было принято къл = 1,0, для марки Ж необходимо ввести къл = 0,85. Поэтому фактическое скорректированное время до воспламенения летучих веществ угля марки Ж составляет твл = = 0,051 0,85 = 0,04335 с и длина начального участка /н = 0,25 • 0,04335/0,04 = 0,27 м. Расхождение следует признать удовлетворительным с учетом всех погрешностей, связанных с экспериментами в работах [7, 8], а также того факта, что битуминозные угли по составу ближе к сапропелитам, а не к гумолитам, для которых построен рис. 1.

Выводы

Сформулирована математическая модель определения начального участка гетерогенного факела, позволяющая связать интегральную зависимость выгорания топлива с пространственными координатами факельного континуума Ф(и) = F(Lф).

Разработанная математическая модель адаптирована к результатам теоретических и экспериментальных работ в этом направлении. Сравнение численных решений и экспериментальных данных с предложенной моделью показало достаточно удовлетворительное согласование, что подтверждает достоверность предложенной математической модели и надежность расчетов по ней для определения основных характеристик факела в зоне интенсивного горения котельных агрегатов.

Литература

1. Математическая теория горения и взрыва / Я.Б. Зельдович, Г.И. Баренблатт, В.Б. Либрович, Г.М. Махвиладзе. -М.: Наука, 1980. - 478 с.

2. Семенов, Н.Н. О некоторых проблемах химической кинетики и реакционной способности / Н.Н. Семенов. - М.: Изд-во АН СССР, 1954. - 350 с.

3. Торопов, Е.В. Концепция факельного континуума для зоны интенсивного горения котельного агрегата / Е.В. Торопов, К.В. Осинцев // Вестник ЮУрГУ. Серия «Энергетика». - 2015. -Т. 15, № 3. - С. 5-10. DOI: 10.14529/power150301

4. Белосельский, Б.С. Технология топлива и энергетических масел / Б.С. Белосельский. - М.: Изд-во МЭИ, 2003. - 340 с.

5. Спейшер, В.А. Повышение эффективности использования газа и мазута в энергетических

установках /В.А. Спейшер, А.Д. Горбаненко. - М.: Энергия, 1974. - 208 с.

6. РД 153-34.1-03.352-99. Правила взрывобе-зопасности топливоподачи и установок для приготовления и сжигания пылевидного топлива. -М.: Энергия, 2000. - 37 с.

7. Бабий, В.И. Горение угольной пыли и расчет пылеугольного факела / В.И. Бабий, Ю.Ф. Ку-ваев. - М.: Энергоатомиздат, 1986. - 208 с.

8. Jones, J.C. Combustion Science: Principles and Practice. - Australia: Millennium Books, 1993. -306p. DOI: 10.1016/0010-2180(95)00222-7

Торопов Евгений Васильевич, д-р техн. наук, профессор, заслуженный деятель науки и техники РФ, профессор кафедры «Промышленная теплоэнергетика», Южно-Уральский государственный университет, г. Челябинск; evtor@mai1 ги.

Осинцев Константин Владимирович, канд. техн. наук, доцент, заведующий кафедрой «Промышленная теплоэнергетика», Южно-Уральский государственный университет, г. Челябинск; osintsev2008@ yandex.ru.

Поступила в редакцию 5 июля 2016 г.

DOI: 10.14529/power160302

MATHEMATICAL MODEL FOR DETERMINING INITIAL SECTION OF HETEROGENEOUS FLAME AND ITS ADAPTATION

E.V. Toropov, evtor@mail.ru,

K.V. Osintsev, osintsev2008@yandex.ru

South Ural State University, Chelyabinsk, Russian Federation

This paper describes a new mathematical model of aerodynamics and heat transfer at the initial section of heterogeneous flame continuum, which largely determine the fuel utilization efficiency of the boiler unit. It provides an analytical description of volatile organic emissions depending on the temperature. This temperature relates to heating coal-dust particles; at that, the ignition of air and volatile-matter mixture is considered both at the known temperature and at the known lower threshold of mixture explosion. These two approaches lead to the same results at determination of an initial flame section. Model reliability is ensured by the known experimental data obtained at determination of concentration limits of air and volatile-matter mixture ignition, period to mixture ignition and other data.

Keywords: flame, initial section, ignition, mathematical model.

References

1. Zel'dovich Ya.B., Barenblatt G.I., Librovich V.B., Makhviladze G.M. Matematicheskaya teoriya goreniya i vzryva [Mathematical Theory of Combustion and Explosion]. Moscow, Nauka Publ., 1980. 478 p.

2. Semenov N.N. O nekotorykh problemakh khimicheskoy kinetiki i reaktsionnoy sposobnosti [About Some Problems of Chemical Kinetics and Reactivity]. Moscow, AN SSSR Publ., 1954. 350 p.

3. Toropov E.V., Osintsev K.V. [Concept of Flame Continuum for Zone of Intense Burning of Boiler Unit]. Bulletin of South Ural State University. Ser. Power Engineering, 2015, vol. 15, no. 3, pp. 5-10. (in Russ.) DOI: 10.14529/power150301

4. Belosel'skiy B.S. Tekhnologiya topliva i energeticheskikh masel [Fuel and Fuel Oil Technology]. Moscow, MEI Publ., 2003. 340 p.

5. Speysher V.A., Gorbanenko A.D. Povyshenie effektivnosti ispol'zovaniya gaza i mazuta v energeticheskikh ustanovkakh [Improving Efficiency of Using Gas and Fuel Oil in Power Plants]. Moscow, Energiya Publ., 1974. 208 p.

6. RD 153-34.1-03.352-99. Pravila vzryvobezopasnosti toplivopodachi i ustanovok dlya prigotovleniya i szhiganiya pylevidnogo topliva [Regulations on Explosion Safety of Fuel Supply and Pulverized Coal Conditioning and Bruning Systems]. Moscow, Energiya Publ., 2000. 37 p.

7. Babiy V.I., Kuvaev Yu.F. Gorenie ugol'noy pyli i raschet pyleugol'nogo fakela [Burning Coal Dust and Coal-dust Flame Calculation]. Moscow, Energoatomizdat Publ., 1986. 208 p.

8. Jones, J.C. Combustion Science: Principles and Practice. Australia, Millennium Books, 1993. 306 p. DOI: 10.1016/0010-2180(95)00222-7

Received 5 July 2016

ОБРАЗЕЦ ЦИТИРОВАНИЯ

Торопов, Е.В. Математическая модель определения начального участка гетерогенного факела и ее адаптация / Е.В. Торопов, К.В. Осинцев // Вестник ЮУрГУ. Серия «Энергетика». - 2016. - Т. 16, № 3. -С. 15-22. DOI: 10.14529/power160302

FOR CITATION

Toropov E.V., Osintsev K.V. Mathematical Model for Determining Initial Section of Heterogeneous Flame and Its Adaptation. Bulletin of the South Ural State University. Ser. Power Engineering, 2016, vol. 16, no. 3, pp. 15-22. (in Russ.) DOI: 10.14529/power160302