CC BY
57
Н. Ф. Тимербаев
МАТЕМАТИЧЕСКОЕ ОПИСАНИЕ ПРОЦЕССА ТЕРМИЧЕСКОЙ ПЕРЕРАБОТКИ
ВЛАЖНЫХ ДРЕВЕСНЫХ ОТХОДОВ МЕТОДОМ ПРЯМОТОЧНОЙ ГАЗИФИКАЦИИ
Ключевые слова: газификация, перенос тепла, перенос влаги, пиролиз, термическое разложение, сушка, синтез-газ, биомасса,
математическое описание, граничные условия, система уравнений.
Разработано математическое описание процесса энерготехнологической переработки влажных древесных отходов с применением предварительной сушки отходов отработанными топочными газами. Процесс газификации рассмотрен как совокупность процессов пиролиза, горения и восстановления.
Keywords: gasification, heat carrying over, moisture carrying over, pyrolysis, thermal decomposition, drying, synthesis-gas, a biomass, the mathematical description, boundary conditions, system of the equations.
The mathematical description of process of thermal processing of a damp wood waste a method of directflow gasification.
Одним из эффективных способов термической переработки влажных древесных отходов является их газификация. Процесс прямоточной газификации представляет собой сложный комплекс физико-химических
превращений твердого топлива в горючий газ, включающий в себя следующие основные процессы: сушка, пиролиз, горение и восстановление [1,2].
Для получения математического описания процесса переработки высоковлажных древесных отходов методом прямоточной газификации необходимо рассмотреть и связать между собой все вышеперечисленные процессы.
В общем виде, при описании процесса тепломассопереноса при сушке отходов деревообработки, состоящих преимущественно из технологической щепы и стружки, внутреннюю задачу можно свести к решению уравнения тепломассопереноса для одномерной симметричной пластины:
IT am (Ї)+ am5& Г P7
др
дг
■ = aA
(д2т, "дб2
A
+
rs
о,
au
дг
(1)
(2)
Поле общего давления внутри материала зависит от температурного режима сушки и свойств капиллярно-пористого коллоидного тела, которым является древесины. В частности, для
пиломатериалов из древесины может быть использовано уравнение, полученное Г.С. Шубиным [3]
Фм
дг
С„ц пористость
д 2Рм 0x2
5U
+ s^P„^
дг
Рм ( дТ„
Тм
тЧ -дт I (3)
где пористость древесного определить из выражения
1 W
С„ = 1-
Рб
материала можно
(4)
чРдв 100-Р,
При рассмотрении частного случая, конвективной сушки древесных отходов топочными газами в режиме противотока в сушильном бункере, изменение влагосодержания и температуры
топочного газа по высоте слоя можно определить из уравнений материального и теплового балансов записанных в следующем виде [4]:
j ■ f
В„
dX_=
dh p„(1 -e)w„ L
dT
q
f
(5)
(6)
^ Сгр0г £ • wг
а изменение влажности и температуры внутри древесной частицы, в условиях сушильного бункера, при отсутствии общего градиента давления и фазовых превращений внутри древесной частицы, можно определить с помощью системы дифференциальных уравнений тепломассопереноса:
ди м
1 д
дh xГ дx
5T,, 1 д
ди м
м дx xr X -5T"
+ Sa,
дTм
дx ) I w
1
1
(7)
(8)
аи хГ ах і"'' ах )с„р„\«м
\ / м і и м
Начальные условия для выражений (5)-(8) запишутся в виде:
Хг(0) =Хгк (9)
Тг(0) =Тгк (10)
и„(х,0)= Ым„ (11)
Т (х,0) = Т„ (12)
Влагосодержание топочного газа Хгк на выходе из сушильного бункера можно определить из выражения материального баланса процесса сушки:
X ГК =
ви(Цмн - имк) + LrXn Lr
(1З)
Температура Тгк топочного газа на выходе из сушильного бункера определяется из
технологических соображений при заданной степени насыщения топочного газа водяными парами. Для этого необходимо решить уравнение Р • X • м
ф= Р„ (М.+М.Я ) (14)
относительно температуры при заданной относительной влажности ф, где давление насыщения можно определить с помощью уравнения:
P = exp18,92
5298,3
(T+2735
1)-exp15;
n. 96486,57 94----------^). (15)
(T+27351
Граничные условия для выражений (7), (8) запишутся в виде:
- на поверхности частицы
"аи„
8х
і = РРи(Р„ - Рср ) = РиЭ„
+ 5аТм ах
(16)
Ч - а (ТГ - ТМ ) - ІГ - м
Тм
(17)
■ в центре частицы при условии симметрии
Т
х
и..
х
- 0.
(18)
Равновесное влагосодержание древесины ^ в зависимости от температуры и относительной влажности топочного газа можно определить из уравнения:
и |°131 ехрЦ848ф-ТГ)/ТГ°’451 при и <ипг
!ипг пРИ и >ипг
где предел гигроскопичности для древесины №1г в зависимости от температуры определяется соотношением:
и„
0,314 -1,39 -10-3Тг.
(20)
Теплоту парообразования в граничном условии (17) в зависимости от температуры древесины можно определить выражением:
617443521 -16323,18 - Тм
-103
_ 1- 0,0011- Тм " (21)
При решении задачи термического разложения древесного материала в зоне пиролиза установки газификации изменение температуры материала описывается дифференциальным уравнением теплопереноса [3,4]
с аТсл - аТсл Г. аТсл 1+ Ч . .,,ч
Рсл' " ' ас ах '[ эф - ах ] Чсл (22)
для решения которого необходимо определить сток тепла на прогрев материала, испарение влаги, и приток тепла от химических реакцией на стадии пиролиза
Чсл - Ч„ + Ч„с„ + Чхр. (23)
Для решения уравнения (22)
сформулированы начальные (24) и граничные условия (25)
Тсл(0,х) - Тна
аТ„
-.,ф
= «-(Тг - Тсл).
(24)
(25)
Сток тепла на испарение влаги определяется выражением
р, аи,
Чесї - г- у
ас
(26)
х=0
Изменение влажности определяется совместным решением уравнений (7), (8).
Для определения притока тепла на стадии пиролиза, древесина рассматривается как смесь компонентов: гемицеллюлозы, целлюлозы и
лигнина. А процесс пиролиза древесины как совокупность процессов термического разложения, ее основных компонентов. Механизм термического разложения компонентов древесины, представлен на рис. 1 и предполагает разложение компонентов
древесины на два конечных продукта: древесный уголь и летучие газы.
(И )Летучие
(угц) Промежуточное вещество
ц |(1-Г„.) Летучіші I (Т„.) Уголь I
■| Целлюлоза [ и* V | +
[ (у„) Уголь |
* [ (1-ул) Летучие |
Лигнин
I (у,) Уголь
Рис. 2 - Механизм термического разложения компонентов древесины
Разложение гемицеллюлозы происходит в два этапа: на первом - гемицеллюлоза разлагается на газы и промежуточный остаток, на втором -промежуточный остаток разлагается на газы и уголь. Целлюлоза и лигнин разлагается на газы и уголь в одну стадию. Массовые доли для гемицеллюлозы, целлюлозы и лигнина зависят от породы древесины.
С учетом принятого механизма изменение массы в единице объема для каждого из указанных компонентов запишется уравнениями химической кинетики для гемицеллюлозы: Зт~.
ат„в
для целлюлозы: для лигнина:
аг--к^т
: Угц ■ к гц ■ тгц - к „в ■ т„в
ат
ас
ат л
ас
---кц -тц
(27)
(28)
(29)
(30)
(31)
В результате теплоту термического разложения древесины можно записать в виде:
Чхр - Чгц - (-к гц - тгц + У гц - кгц - тгц - к„в - т„в ) +
+ Чц - (-кц - тц ) + Чл - (-кл - тл ).
Уравнения баланса массы для угля и газа запишутся в виде:
--У„в - к„в - т„в +Уц - (-кц - тц ) + Ул - (-кл - тл ), (32)
ат у
аг
Ф-ті) - а(w„г т
+(1 - Угц) - кгц -тц+(1-у„в)-к1в-ті. +
5с 87 ............... "■ "■ (33)
+(1-Уц)\-т+(1-Ул)-клт-
Предполагая, что при пиролизе температура частицы равна температуре образующихся газов, а в локальном объеме между компонентами газовой и твердой фазы установлено термодинамическое равновесие, уравнение сохранения энергии запишется в виде
. . 8Т 8Т ( 8Т ^ 8Т /■-, л \
(с,цт,ц + с,т, + слтл + с,т, + с„,т„,) -8с = -87 ту]_ тпгс»'«» -87 + ч<р (34)
Коэффициент теплопроводности зависит от доли прореагировавшей древесины, и определяется выражением
. - (1 - л) - .м + Л - .у + Е - .„г +
135 - с - Т3 - ^
(35)
Процесс горения продуктов пиролиза описывается уравнениями, позволяющими
определить расход необходимого воздуха для
х
х-0
х-0
х=0
горения, температуру продуктов сгорания и их компонентный состав.
При рассмотрении процесса
восстановления, вследствие значительного влияния конвективного переноса тепла и массы, теплопроводностью и диффузией по газу пренебрегаем. С учетом принятого допущения уравнения сохранения вещества для каждого компонента газового потока и для угля запишутся в виде:
wг г • — = -к, • С,,
8у 11 (36)
8т у I
wу----- = -ту 'V к,.
у 8у у т=11 (37)
Уравнение сохранения энергии для газового потока и угля соответственно примут вид:
рг-г •сг-г ^ •8у=-ог-г •(Ту -Тr-r)-f+3я •к •(С° -С)),
8у =’ (38)
Ру •су • wу •8у=ог-г •(Ту -Тr_r)•f-яя •к •(Со -С)). (39)
8у ,=1 (39)
Для решения системы уравнений (36-39) приняты граничные условия
Ч=° = и=°; с,|у=° = С, о; т ,1,=° =т,° (40)
Рис. 2 - Алгоритм расчета комплексной
энерготехнологической переработки древесных отходов
На рис. 2 представлен алгоритм расчета комплексной энерготехнологической переработки древесных отходов методом прямоточной газификации. Расчет ведется в следующей последовательности. После ввода постоянных и варьируемых параметров проверяется
необходимость расчета процесса сушки. Если отходы высоковлажные ин>инг, т.е.
влагосодержание древесных отходов больше требуемой начальной влажности отходов перед газификацией, то рассчитывается процесс сушки либо за счет тепла отработанных топочных газов, либо от рекуперативного тепла продуктов
газификации. Если скорость сушки древесных отходов меньше предельной заданной скорости, то рассчитывается дополнительное топочное устройство. Далее проводится расчет процесса газификации. Расчет процесса газификации проводится в соответствии с блоком, представленным на рис 4. Блок расчета процесса газификации включает в себя: расчет зоны термического разложения древесного материала в отсутствие кислорода (пиролиза), расчет процесса горения продуктов пиролиза и расчет восстановительной зоны процесса газификации. Оптимизационная задача расчета процесса
газификации сводится к определению высот зон пиролиза, горения и восстановления, увязанных между собой таким образом, чтобы количество и качество угля образованного в зоне пиролиза удовлетворяло количеству газов получаемых в зоне окисления, с тем, чтобы на выходе из зоны
восстановления не оставалось диоксида углерода, паров воды, либо непрореагировавшего углерода.
Исследования по данной работе выполнены в рамках реализации федеральной целевой программы «Исследования и разработки по приоритетным направлениям научно-технологического комплекса России на 2007-2013 годы по теме: «Создание технологии и опытной установки комплексной переработки древесных отходов лесной
промышленности с получением теплоизоляционного материала», при финансовой поддержке Министерства образования и науки Российской федерации.
Литература
1. Тимербаев, Н.Ф. Моделирование процесса сушки влажных древесных отходов отработанными газами котельных установок / Н.Ф. Тимербаев, А.Н. Грачев, Р.Г. Сафин // Известия высших учебных заведений «Химия и химическая технология». - 2006 г. - Том 49. - Вып. 11. -С. 103-105.
2. Тимербаев, Н.Ф. Моделирование процесса сушки древесных частиц при кондуктивном подводе тепла / Н.Ф. Тимербаев, Р.Г. Сафин, А.Р. Хисамеева // Вестник Казан. технол. ун-та. - 2011. - Т. 14, № 4. -С.84-90.
3. Тимербаев, Н.Ф. Моделирование процесса прямоточной газификации древесных отходов / Н.Ф. Тимербаев, Р.Г. Сафин, А.Р. Садртдинов, А.Р. Хисамеева // Вестник Казан. технол. ун-та. - 2011. - Т. 14,№ 7. -С.75-80.
4. Тимербаев, Н.Ф. Исследование восстановительной зоны
процесса газификации древесных отходов /
Н.Ф. Тимербаев, Р.Г. Сафин, З.Г Саттарова,
Д. А. Ахметова // Вестник Казан. технол. ун-та. -2011. - Т. 14,№ 8. -С.90-96.
5. Тимербаев, Н.Ф. Моделирование процесса пиролиза древесины в установке для производства древесного угля / Н.Ф. Тимербаев, Р.Г. Сафин, И.И Хуснуллин // Вестник Казан. технол. ун-та. - 2011. - Т. 14, № 9. - С.51-57.
6. Тимербаев, Н.Ф. Кондуктивный теплообмен дисперсного материала в установке для производства древесного угля / Н.Ф. Тимербаев, Р.Г. Сафин, И.И Хуснуллин // Вестник Казан. технол. ун-та. - 2011. - Т. 14, № 18. - С.69-76.
© Н. Ф. Тимербаев - канд. техн. наук, доц. каф. переработки древесных материалов КНИТУ, tnail@rambler.ru.
