CC BY
34
DOI 10.23859/1994-0637-2018-1-83-6 УДК 674.047.3
© Синицын Н.Н., Телин Н.В., Грызлов В.С., Андреев А.С., Виноградова М.С., Гаркавченко Э.В., Гневашева Т.В., Кузнецова В.П., Павлова А. И., 2018
Синицын Николай Николаевич
Доктор технических наук, профессор, Череповецкий государственный университет (Череповец, Россия) E-mail: [email protected]
Телин Николай Владимирович
Доктор технических наук, профессор, Череповецкий государственный университет (Череповец, Россия) E-mail: [email protected]
Грызлов Владимир Сергеевич
Доктор технических наук, профессор, Череповецкий государственный университет (Череповец, Россия) E-mail:[email protected]
Андреев Александр Сергеевич
Кандидат технических наук, преподаватель, Череповецкое высшее военное инженерное училище радиоэлектроники (Череповец, Россия) E-mail:[email protected]
Sinitsyn Nikolay Nikolaevich
Doctor of Technical Sciences, professor, Cherepovets State University (Cherepovets, Russia) E-mail: [email protected]
Telin Nikolay Vladimirovich Doctor of Technical Sciences, professor, Cherepovets State University (Cherepovets, Russia) E-mail: [email protected]
Gryzlov Vladimir Sergeevich
Doctor of Technical Sciences, professor, Cherepovets State University (Cherepovets, Russia) E-mail:[email protected]
Andreev Alexander Sergeevich
PhD in Technical Sciences, teacher, of the Cherepovets Higher Military School of Engineering and Radioelectronics (Cherepovets, Russian) E-mail:[email protected]
Виноградова Мария Сергеевна
Инженер, ООО «Севесталь-проект»
(Череповец, Россия)
E-mail: [email protected]
Гаркавченко Элина Вадимовна
Инженер, ПАО «Северсталь» (Череповец, Россия) E-mail: [email protected]
Гневашева Татьяна Вадимовна
Инженер, ПАО «Северсталь»
(Череповец, Россия)
E-mail: [email protected]
Кузнецова Виктория Павловна
Инженер, ПАО «Северсталь» (Череповец, Россия) E-mail: [email protected]
Павлова Алина Игоревна
Инженер, ООО «АВПоиск»
(Череповец, Россия)
E-mail: [email protected]
Vinogradova Maria Sergeyevna
Engineer, OOO "Severstal-project" (Cherepovets, Russia) E-mail: [email protected]
Garkavchenko Elina Vadimovna
Engineer, PJSC "Severstal" (Cherepovets, Russia) E-mail: [email protected]
Gnevacheva Tatyana Vadimovna Engineer, PJSC "Severstal" (Cherepovets, Russia) E-mail: [email protected]
Kuznetsova Victoria Pavlovna
Engineer, PJSC "Severstal" (Cherepovets, Russia) E-mail: [email protected]
Pavlova Alina Igorevna
Engineer, OOO "AVSearch" (Cherepovets, Russia) E-mail: [email protected]
МОДЕЛИРОВАНИЕ ДИНАМИКИ ВЫХОДА ВЛАГИ И ЛЕТУЧИХ ВЕЩЕСТВ В ПРОЦЕССЕ НАГРЕВА ДРЕВЕСНОЙ БИОМАССЫ (КОРЫ)1
STUDY ON OPPORTUNITIES AND EVALUATION OF TIME OF WOOD CRUDE DRYING WITH WASTE HIGH-TEMPERATURE SMOKE GASES OF BOILER UNITS
1 Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта №17-48-350810.
Аннотация. В работе предложена модель процесса горения на стадиях сушки и выхода летучих веществ из куска высоковлажной коры древесины, состоящая из системы нестационарных дифференциальных уравнений. Получены зависимости числа Фурье от критерия Коссовича и критерия температурных отношений.
Ключевые слова: модель, древесная кора, сушка, влага, летучие вещества, фазовый переход влаги
Abstract. In this article a model of the combustion process at the stages of drying and exit of volatile substances from the bush of the high-moisture crust of wood, consisting of a system of non-stationary differential equations is suggested. Correlation of the Fourier number, the Kosovich test and the temperature relationship criterion are obtained.
Keywords: model, tree bark, drying, moisture, volatile substances, phase transition of moisture.
Введение
Технологии получения энергии из древесных отходов непрерывно развиваются и совершенствуются. Основной особенностью твердых природных топлив является их способность при нагревании выделять продукты термического разложения органической массы (летучие). Точный расчет состава продуктов термолиза твердых природных топлив в зависимости от температуры и времени представляет собой сложную задачу, причем кинетические характеристики, необходимые для этого расчета, еще не обобщены. В ряде инженерных задач не требуется знать состав продуктов термолиза, но нужны ответы на более простые вопросы: какое количество летучих может выделиться из топлива, за какое время и в каком процентном отношении произойдет разложение топлива при его термолизе в условиях данного промышленного агрегата. Определяющим фактором процесса горения является динамика выхода влаги и летучих веществ из древесной биомассы при различных температурах. Для организации процесса горения побочных продуктов заготовки, обработки и переработки древесины необходимо знать кинетические характеристики этого процесса. Решение этой проблемы позволит обеспечить безотходность переработки древесины на деревоперерабатывающих предприятиях.
Основная часть
В упрощенном виде процесс горения древесной коры может быть разделен на четыре зависящих от температуры стадии: сушка; выход летучих компонентов; возгорание газообразных летучих веществ и горение твердого углерода. После попадания топлива в топку его температура сразу начинает увеличиваться. Испарение влаги происходит при достаточно низкой температуре (50...100 °С) [2]-[4]. На испарение воды, содержащейся в исходном топливе, уходит часть энергии, выделяемой в процессе горения топлива, что уменьшает температуру в топке и замедляет процесс его сжигания. Процесс выделения летучих веществ (прежде всего углеводородов) начинается при достижении температуры топлива 100.105 °С [2]-[4]. В результате испарения влаги и выделения летучих веществ топливные частички становится пористыми. Скорость выхода летучих веществ увеличивается с повышением температуры. Сначала происходит разложение гемицеллюлозы и затем, при более высокой температуре, разложение целлюлозы. Полное выделение большинства летучих веществ происходит при температуре до 400 °С [2]-[4]. С дальнейшим повышением температуры скорость выхода летучих компонентов резко возрастает. Однако при температуре 400.500 °С может наблюдаться низкая скорость выхода летучих веществ, определяемая скоростью разложения лигнина, увеличивающейся с повышением температуры. Возгорание газообразных летучих веществ происходит около поверхности
твердого слоя топлива при температуре около 600 °С. Газы сгорают в виде пламени, которое нагревает твердые частицы. При сжигании древесины приблизительно 80 % энергии выделяется в виде газов, а оставшаяся часть в виде древесного угля.
Горение твердого углерода происходит в горящем слое при температуре свыше 800 °С. Основное тепловыделение при сжигании древесной коры ввиду высокого выхода горючих веществ у этих топлив происходит в топочном объеме, а не в слое топлива. Необходимым условием использования древесной коры в качестве экологически приемлемого вида топлива является создание условий ее эффективного и полного сгорания. Это достигается пространственной и временной организацией процесса горения на всех его стадиях. Целью данной работы является моделирование динамики выхода влаги и летучих веществ в процессе нагрева древесной биомассы (коры).
На основе анализа экспериментальных данных, оттестированных методик расчета отдельных процессов составлена модель процесса горения на стадиях сушки и выхода летучих веществ из куска высоковлажной коры древесины, состоящая из системы нестационарных дифференциальных уравнений: уравнения теплопроводности, уравнения испарения влаги и выхода летучих компонентов, учитывающих перемещение фронтов испарения влаги и выхода летучих веществ.
Прогрев частицы коры на стадиях сушки и выхода летучих компонентов описывается сквозным уравнением теплопроводности с переменными граничными условиями третьего рода, учитывающими теплообмен конвекцией и тепловым излучением. Для периода сушки и выхода летучих веществ система уравнений имеет вид:
Сэф(^(Г) Ы И (1)
интегрируемое в области: 0 < х < 5, 0 <г<гк начальное условие: при: т = 0 Т| = Т0;
граничное условие: при х = 0
-А(Т)^т = «(Т - т )
ах
при х =
-А.(Т) ^ = 0 . ах
где р - плотность материала; а - коэффициент теплоотдачи; 5 - половина толщины пластины; Тср - температура среды; Т - начальная температура материала; 1 - коэффициент теплопроводности.
При этом выделение теплоты фазового перехода в уравнении (1) учитывают с помощью эффективной теплоемкости сэф, задаваемой выражением:
Ч (Т) ,Т > Тл;
С (Тс) )-у + с(Тл)-(1 -у) + ^ , Тл < Т < Тс с2 (Т), Т <Тс
С =
коэффициент теплопроводности и плотность определяют по формулам:
А1, Т < Т с
X-у + А 2-(1 -у), Т с < Т < Т л
А2, Т > Тл
[р1, Т < Тс; Р = ]р1 -У + Р 2 '(1 -V), Т с < Т < Т
[Р2, Т > Тл ,
где Тл = Тф +19,5 , Тс = Тф - 37 - фиктивные температуры начала и окончания фазового перехода воды; с(Т) - теплоемкость материала; с1 и с2 - теплоемкость сухого и влажного слоев материала; р! и р2 - плотность сухого и влажного слоев материала;
и Х2 - коэффициенты теплопроводности сухого и влажного слоев материала; g -доля влаги в элементарном объеме материала; - половина толщины слоя материала; Ь - удельная теплота фазового перехода влаги; ¥ - доля влажного материала; Тф - температура фазового перехода воды. Величина ¥ определяется по формуле:
V =
1, Т < Т,
Т л - Т Т -Т '
Т с < Т < Т
0, Т > Тл.
а = ак +ал ,
где ак - коэффициент теплоотдачи конвекцией; ал - коэффициент теплоотдачи излучением.
Процесс термолиза твердого топлива представляет собой разрыв физико-химических связей внутри его составляющих с последующим весьма быстрым образованием конечных продуктов разложения. Поэтому можно выделить отдельные группы связей, которые разрушаются при определенном условии с определенной скоростью, примерно одинаковой для данной группы связей. В этом случае скорость выделения продуктов термолиза, образующихся в результате разрушения такой группы связей, можно записать следующим образом [5]:
• ^ = к, • (1 - У)* ,
с ¡и <1 с0(
где Уг - доля продуктов термолиза 1-й группы в общем количестве продуктов термолиза; т - время разрушения связей г-й группы; с0г - доля продуктов термолиза г-й
группы в общем количестве продуктов термолиза, получающаяся при полном разрушении связей данной группы (0 < у < с0г); пр - порядок реакции (ниже принимаемый равным 1, что для реакций разложения в первом приближении допустимо);
к. = к„. • е
Е
где к - константа скорости реакции, характеризующей разрушение связи данной группы; Ео - знергия активации, МДж/кмоль; Я - универсальная газовая постоянная;
Т - температура материала в узле; к. - предъэкспоненциальный множитель, 1/с;
т - время.
При проведении расчетов принято, что максимальный выход летучих веществ в расчете на сухую массу коры древесины составляет V¡Cax = 83,6%; E = 59,9 МДж /кмоль ; k¡ = 38,3 с"1 [5]. В формуле (1) принято, что сы = 1.
В этом случае скорость выделения продуктов термолиза, образующихся в результате разрушения такой группы связей, описывается уравнением (2):
дУ , — = ехр а г
" Е КТА
(1" V), е< х < 5. (2)
Изменение массы образца коры древесины в течении времени прогрева за счет выхода влаги и летучих веществ описывается уравнением
атч аш ^ ау (3)
а г а г а г
Решение системы уравнений (1)-(3) получено численным методом. Уравнение (1) решено методом конечных разностей по явной схеме аппроксимации производных. Уравнения (2), (3) решены методом Рунге-Кутты четвертого порядка.
На рис. 1 представлена зависимость изменения температуры поверхности образца от времени прогрева. Начиная с температуры 180 оС, температура поверхности резко увеличивается и из поверхностных слоев увеличивается выход летучих веществ.
На рис. 2 представлен взаимосвязь чисел Бо от Ко [1] при отсутствии излучающего объема Ш/ = 85 %. Из рисунка видно, что зависимость чисел Бо и Ко имеет
вид параболы, где Ео = "22 г - число Фурье, Ко = - ^ ё Р2
§2 ' С1Р1 (Тф "То )
Рис. 1. Зависимость температуры поверхности куска коры древесины от времени при Тср= 352 °С
На рис. 3 представлено сравнение расчетных и экспериментальных данных по выходу влаги летучих веществ. Начиная с момента времени т = 270 с, масса образца изменяется незначительно, а начиная с момента времени - 300 с масса образца начинает снова уменьшаться за счет выхода летучих веществ. Оценка погрешности отклонения расчетных, по предложенной методике, и экспериментальных данных допустима.
Fo
Рис. 2. Зависимость числа Бо от Ко при отсутствии излучающего объема при W¡' = 85 %.
ш, 10-5 кг
Рис. 3. Зависимость массы прогреваемого куска коры от времени 1 - расчетные данные, 2 - экспериментальные данные
Таким образом, предложенная математическая модель прогрева влажной коры древесины достаточно точно описывает процесс горения на стадиях сушки и выхода летучих веществ. Отклонение расчетных данных от экспериментальных не превышает 12 %. Математическую модель можно использовать при разработке технических устройств, обеспечивающих сушку коры древесины и выход летучих веществ, например, при сжигании высоковлажной коры.
Выводы
Разработана математическая модель прогрева высоковлажной коры древесины с учетом фазового перехода влаги и выхода летучих веществ. Получены графические зависимости безразмерного времени прогрева высоковлажной коры древесины от характеризующих отношение количества теплоты, выделяющегося при испарении влаги, к количеству теплоты идущего на прогрев материала до температуры греющей среды и критерия Кт.
Литература
1. Лыков А.В. Теория теплопроводности. М.: Высшая школа, 1967. 600 с.
2. Любов В.К., Любов С.В. Повышение эффективности энергетического использования биотоплива. Архангельск, 2010. 496 с.
3. Любов В.К. Энергетическое использование биотоплива. Архангельск: АГТУ, 2007. 156 с.
4. Марьяндышев П.А., Чернов А.А., Любов В.К. Анализ термогравиметрических и кинетических данных различных видов древесного биотоплива северо-западного региона Российской Федерации // Лесной журнал. 2016. №1. С. 167-172.
5. Померанцев В.В. и др. Основы практической теории горения. Л.: Энергоатомиздат, 1986. 312. с.
References
1. Lykov A.V. Teoriia teploprovodnosti [Theory of heat conductivity]. Moscow: Vysshaia shkola, 1967. 600 p.
2. Lyubov V.K., Lyubov S.V. Povysheniie effektivnosti energeticheskogo ispol'zovaniia biotop-liva [Increase the effectiveness of energy use of biofuels]. Arkhangelsk, 2010. 496 p.
3. Lyubov V.K. Energeticheskoie ispol'zovaniye biotopliva [Energy use of biofuel]. Arkhangelsk: AGTU, 2007. 156 p.
4. Mar'iandyshev P.A., Chernov A.A., Liubov V.K. Analiz termogravimetricheskikh i kineti-cheskikh dannykh razlichnykh vidov drevesnogo biotopliva severo-zapadnogo regiona Rossiyskoi Federatsii [Analysis of thermogravimetric and kinetic data of different types of wood biofuel in the northwestern region of the Russian Federation]. Lesnoi zhurnal [Forest Journal], 2016, no. 1, pp. 167-172.
5. Pomerantsev V.V. Osnovy prakticheskoi teorii goreniia [Bases of the practical theory of burning]. Leningrad: Energoatomizdat, 1986. 312 p.
Для цитирования: Синицын Н.Н., Телин Н.В., Грызлов В.С., Андреев А.С., Виноградова М.С., Гаркавченко Э.В., Гневашева Т.В., Кузнецова В.П., Павлова А.И. Моделирование динамики выхода влаги и летучих веществ в процессе нагрева древесной биомассы (коры) // Вестник Череповецкого государственного университета. 2018. №2(83). С. 47-53. DOI 10.23859/ 1994-0637-2018-1-83-6
For citation: Sinitsyn N.N., Telin N.Vl., Gryzlov V.S., Andreev A.S., Vinogradova M.S., Gar-kavchenko E.V., Gnevacheva T.V., Kuznetsova V.P., Pavlova A.I. Study on opportunities and evaluation of time of wood crude drying with waste high-temperature smoke gases of boiler units. Bulletin of the Cherepovets State University, 2018, no. 2 (83), pp. 47-53. DOI 10.23859/1994-06372018-1-83-6
