CC BY
5УДК 504.054
ОПРЕДЕЛЕНИЕ СКОРОСТИ ГОРЕНИЯ УГЛЯ ПО ГЕТЕРОГЕННОЙ
КИНЕТИЧЕСКОЙ МОДЕЛИ
Николай Викторович Бучилин
кандидат технических наук, доцент isk115599@rambler.ru Алексей Васильевич Аксеновский кандидат сельскохозяйственных наук, доцент
noky2002@mail.ru Иван Павлович Криволапов кандидат технических наук, доцент
ivan0068@bk.ru Сергей Юрьевич Щербаков кандидат технических наук, доцент Scherbakov78@yandex.ru Мичуринский государственный аграрный университет
г. Мичуринск, Россия
Аннотация. Проблема возникновения пожаров и контроля процесса горения органических материалов и угля на сегодняшний день остаётся актуальной. Работа посвящена изучению методик расчёта степени превращения (полноты сгорания) угля по гетерогенной кинетической модели реакции. Рассматривается влияние геометрической формы сгораемых частиц на полноту сгорания. Получаемые расчётные значения степени превращения позволяют оценить вероятное время, необходимое для сгорания угля.
Ключевые слова: сжигание угля, неконтролируемое горение, кинетика гетерогенных реакций, константа скорости реакции
Пожарную опасность представляют вещества и материалы, которые в силу своих свойств способны приводить к возникновению или развитию пожара. Математическое моделирование процессов горения позволяет производить оценку пожарной опасности веществ и материалов [1-3]. В проведённой ранее работе приводилось математическое описание процесса горения твёрдых частиц угля идеальной сферической формы [4]. В данной статье рассматривается математическое описание процесса горения твёрдых частиц различной формы.
При горении частиц площадь поверхности убывает со временем протекания реакции [5-7]. Кинетическое уравнение, описывающее этот процесс, имеет следующий вид:
dM
= kFC (1)
dт
где: М - масса частицы в данный момент времени т; к - константа скорости химической реакции; F - площадь поверхности частицы; С -концентрация кислорода в реакционной зоне.
Большинство частиц угля имеют преимущественно окатанную форму, близкую к форме шара [8-9]. Для таких частиц площадь поверхности и масса в данный момент времени сгорания вычисляются по формулам:
4
F = 4тсг2 => М = -пг3р (2)
где: г - радиус частицы; р - плотность частицы.
Следовательно:
2 2 ( 3 \3 2 dM ( 3 \3 2
F = 4rc(-—) хМз => - — = 4ктс(-—) хМзхС => \4пр/ dт \4пр/
Г^М Г ( 1 Л
-I — = к' I dт => 3(М3 - Мз) = к'т (3)
^Мо Мз ^0 V /
где к' - приведённая константа скорости горения продукта при условии постоянной концентрации кислорода в зоне горения:
fc'=44i)3xC (4)
Для частиц в форме куба используются уравнения:
F = 4г3 => М = г3р (5) где: г - длина ребра кубической частицы; р - плотность частицы. Подставляя формулы (3) в формулу (1), получаем:
1 2 /М\3 /Жз
г = (_) => г = 6у (б)
где: г - радиус частицы; р - плотность частицы.
2
ЙМ /М\3 2 / 1 1\
- — = 6к(—) С = ^'М3 => 3 ( М3 - мз) = (7) Получается уравнение, идентичное уравнению для идеального шара.
Для частиц вытянутой цилиндрической формы, у которых длина (l) много больше радиуса (г), можно пренебречь реакцией на торцах. Поэтому в процессе протекания реакции их длина остаётся постоянной (l = const). Получим следующие уравнения:
F = 2nrl => M = пг01р (8)
Подставляя формулы (8) в формулу (1), получаем:
1 i i ( М \2 ( М \2 /пК2 1
Г = Ы => р = 2п1Ы = 2(р) M0 (9)
1
dM i i rMdM
i fMdM Г
, = —| M0C = fc'M2 => -| —r = fc'ldr (10)
i i r2
2(M02 -M2) = fc'x (11)
Как было описано ранее в статье [4], для более простого выражения скорости реакции используется понятие степени превращения (полноты сгораня) а, которая выражается в долях от единицы либо процентах:
М0 - М
^М-т- (12)
Формулы для расчёта степени превращения горящего порошкообразного материала, состоящего из частиц различной формы, и имеющих одинаковый размер, приведены ниже:
Сферические частицы. В формулах г - радиус частицы.
г3 1
а = 1--з => г = г0(1 - а)3 , (13)
Г0
1 ^хСхт 1 - (1 - а)3 =--(14)
ТоХр
Кубические частицы. В формулах г - длина стороны частицы.
г3 1
а = 1--з => г = г0(1 - а)3 , (15)
го
1 2^хСхт
1 - (1 - а)3 =--(16)
тохр
Вытянутые цилиндрические (стержневидные) частицы. В формулах г радиус цилиндрической частицы.
г2 1
а = 1--2 => г = го(1 - а)2 , (17)
го
1 ^хСхт
1 - (1 - а)1=-^хт" (18)
Приведённые уравнения в общем виде позволяют производить математическое описание гетерогенных химических реакций с различной геометрией поверхности раздела. В частности, этими уравнениями описывается
полнота сгорания угля при известных значениях констант скорости реакции горения. Константа зависит от присутствия катализатора, либо ингибитора в реакционной зоне, однако не зависит от геометрии сгораемых частиц. В расчётах константу скорости сгорания угля без использования каталитических добавок принимали равной 0,039. Таким образом, получаются зависимости полноты сгорания частиц угля определённой геометрии от времени, представленные на рисунке 1. В расчётах плотность угля принималась равной 1290 кг/м3, а средний размер частиц (г) - равным 5 см.
3000 3600
т, сек
Рисунок 1 - Зависимость полноты сгорания (а) от продолжительности горения (т)
угля с различной геометрией частиц. * - сферические частицы; - кубические
частицы; - цилиндрические (стержневидные) частицы
Как видно из представленных формул, для частиц идеальной кубической и сферической формы зависимость полноты сгорания от времени пропорциональна кубическому корню, а для частиц стержневидной цилиндрической формы - квадратному корню. Расчётные данные показывают, что за 1 час (3600 секунд) будет происходить полное сгорание частиц угля с любой геометрией. Наибольшей скоростью сгорания обладают частицы идеальной кубической формы - они полностью сгорают через 30 минут (1800
секунд) после начала возгорания. Объясняется это двумя факторами: в расчётах величина (г) соответствует полной длине стороны куба, а для сферических частиц - только половине диаметра, т.е .фактически кубические частицы в длину в 2 раза меньше кубических. Кроме того, в выбранных математических приближениях кубическая форма частиц диктует проявление вдвое большей скорости горения относительно частиц сферической формы (см. формулы 14 и 16).
Наименьшей скоростью горения обладают частицы стержневидной формы: полнота их сгорания достигает 100 % только через час горения. При этом так как для данных частиц зависимость полноты сгорания от времени пропорциональна не кубическому, а квадратному корню, то их сгорание происхолит наиболее плавно.
Представленная методика расчёта может быть использована для прогнозирования последствий горения твёрдых порошкообразных и гранулированных веществ. Для более точных расчётов необходимо учитывать влияние более сложной геометрии реальных частиц, а также влияние скорости движения конвективных потоков воздуха при горении.
Список литературы:
1. Коломиец А.А., Манаенков К.А., Найденов А.А. Оценка показателей надежности автотранспортных средств // Наука и Образование. 2021. Т. 4. № 1. C. 47.
2. Щербаков С.Ю., Криволапов И.П., Стрельников Д.И., Коробельников
A.П. Характеристика методов проведения анализа риска // Наука и Образование. 2019. Т. 2. № 4. С. 253.
3. Клименко Н.Н., Нистратов А.В., Киселева К.И., Делицын Л.М., Сигаев
B.Н. Применение вторичного углеродного волокна для армирования композиционного материала на основе щелочеактивированного доменного шлака // Стекло и керамика. 2020. № 11. С. 28-31.
4. Бучилин Н.В., Аксеновский А.В., Щербаков С.Ю. Кинетика ингибирования процессов горения угля // Наука и Образование. 2022. Т. 5. № 3.
5. Картечина Н.В., Макова Н.Е., Шацкий В.А., Дорохова А.М. Информационная модель учета сельскохозяйственной техники // Наука и Образование. 2021. Т. 4. № 1. C. 40.
6. Ипполитов Е.Г., Артемов А.В., Батраков В.В. Физическая химия. М.: Издательский центр «Академия». 2005. 448 с.
7. Бучилин Н.В., Никитина В.Ю., Луговой А.А., Варрик Н.М., Бабашов В.Г. Получение высокопористых керамических материалов на основе алюмо-магнезиальной шпинели // Стекло и керамика. 2020. № 10. С. 7-14.
8. Строкова Я.А., Клименко Н.Н. Комплексная щелочно-щелочноземельная активация гранулированного доменного шлака // Успехи в химии и химической технологии. 2019. Т. 33. № 4. С. 130-132.
9. Щербаков С.Ю., Аксеновский А.В., Криволапов И.П., Куденко В.Б. Оценка уровня обеспеченности и повышение пожарной безопасности на складах хранения нефтепродуктов предприятий АПК // Сб. научн. Трудов, посвящённый 85-летию Мичуринского государственного аграрного университета. Мичуринск. 2016. Т. 4. С. 110-114.
UDC 504.054
DETERMINATION OF COAL BURNING RATE USING A HETEROGENEOUS KINETIC MODEL
Nikolai V. Buchilin
Candidate of Technical Sciences, Associate Professor
isk115599@rambler.ru Alexey V. Axenowskiy Candidate of Agricultural Sciences, Associate Professor
noky2002@mail.ru
Ivan P. Krivolapov
Candidate of Technical Sciences, Associate Professor
ivan0068@bk.ru Sergey Yu. Sherbakov Candidate of Technical Sciences, Associate Professor
Scherbakov78@yandex.ru Michurinsk State Agrarian University Michurinsk, Russia
Abstract. The problem of the occurrence of fires and control of the combustion processes of organic materials and coal remains relevant today. The work is devoted to the study of methods for calculating the degree of conversion (combustion completeness) of coal using a heterogeneous kinetic model of the reaction. The effect of the geometric shape of combustible particles on the completeness of combustion is considered. The resulting calculated values of the degree of conversion allow us to estimate the probable time, required for the coal combustion.
Keywords: coal burning, uncontrolled combustion, kinetics of heterogeneous reactions, reaction rate constant.
Статья поступила в редакцию 12.11.2022; одобрена после рецензирования 02.12.2022; принята к публикации 20.12.2022.
The article was submitted 05.11.2022; approved after reviewing 02.12.2022; accepted for publication 20.12.2022.
