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

УДК 66.092 - 977

Д. Ф. Зиатдинова

МАТЕМАТИЧЕСКОЕ ОПИСАНИЕ ТЕРМОХИМИЧЕСКИХ ПРОЦЕССОВ ПЕРЕРАБОТКИ

ДРЕВЕСНЫХ МАТЕРИАЛОВ, СОПРОВОЖДАЮЩИХСЯ ВЫБРОСАМИ ПАРОГАЗОВОЙ ФАЗЫ

Ключевые слова: испарение, конденсация, термическое разложение, перенос тепла, перенос массы, гомогенные процессы,

гетерогенные процессы, граничные условия.

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

Keywords: evaporation, condensation, thermal expansion, heat transfer, mass transfer, homogeneous processes, heterogeneous

processes, boundary conditions.

Developed a mathematical description of the thermochemical processes of processing of wood, fitted with a combined gas cleaning system released gas mixture.

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

ы

VIs

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

Структурная схема (рис. 1) включает в себя [2] участок подготовки древесного сырья 1 (сушка до необходимого влагосодержания, измельчение до требуемой фракции, прогрев до определенной температуры, предварительная химическая обработка), реактор для проведения термохимического процесса 18 (термомодификация, пластификация, варка, термическая переработка), участок разгрузки переработанного древесного материала 16 (выдувной резервуар, разгрузочная площадка), систему газоочистки (конденсацию, абсорбцию, адсорбцию, каталитическую очистку).

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

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

да из камеры образующейся парогазовоздушной смеси 2. В результате осуществления технологической операции, сопровождаемой иногда химической реакцией между компонентами композиции, в камере образуются потоки паров и неконденсирующихся газов. Для обеспечения нормальных санитарных условий в помещении (концентрация выделяющихся паров и газов в воздухе рабочей зоны должна поддерживаться на уровне, не превышающем предельно допустимой концентрации - ПДК) камера должна быть сообщена устройством откачки парогазовой смеси с технологического оборудования, также в помещении должна быть приточно-вытяжная вентиляция 9, 19. Для очистки отходящих газов перед выбросом их в атмосферу на линии устанавливаются устройства газоочистки: конденсатор - 5,7, абсорбер, 11 - адсорбер или их комбинация. Согласно представленной схеме осуществляются существующие процессы переработки древесных материалов на большинстве предприятий ЛПК. Анализ способов осуществления процессов в деревообрабатывающей промышленности по данной схеме позволил выявить несколько возможных вариантов их реализации.

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

Для исключения выделения образующихся паров и газов в рабочую зону объемная производительность системы откачки должна быть больше

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

В соответствии с формализацией процесса, древесину рассматриваем как смесь, состоящую из трех главных составляющих: гемицеллюлозы, целлюлозы и лигнина с известными массовыми долями [3,4],. С учетом принятых допущений изменение массы в единице объема для ьго компонента запишется в общем виде

дт.

-Г- _ ктг (1)

дт

Изменение массы газовой фазы в единице объема представлена как сумма массовых потоков за счет конвекции и химических реакций:

Фт)_ %г-т)

дт

где е - пористость, определяемая из выражения

(

а _

1 --

т

др

т

И1 -адр )

(2)

(3)

др 0

Скорость газового потока определяется по закону Дарси:

^пг

АР

~Г

(4)

где цпг - динамическая вязкость газа; Кч - газовая проницаемость

частицы Кч рассчитывается по формуле:

Кч _ (1 -V). Кр + v -Ку, (5)

где п - степень химического превращения, определяемая по формуле

V _ 1 --

т

др

т

др 0

(6)

Давление газа определяется из уравнения состояния идеальных газов:

Р _ Рш ^ Тпг („)

Мк ' (7)

где М - молекулярный вес; Я - универсальная газовая постоянная.

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

. дТ

(Сгитги + Сц+ Слтл + Спгтпг _

(8)

дт

дТ Л дП дТ

¥ Г ¥ |-^ ^ + ЯХ

Яхр _ Ягц ■ (-к.

т

гц гц

)

+

+ Яц -(-кц-тц) + Ял -(-кл ■тл)

(9)

где ^ - коэффициент теплопроводности частицы древесного материала; т - масса вещества в единице объема; q - энергия, выделяемая (поглощаемая) химической реакцией; к - константа скорости химической реакции, определяется по закону Аррениуса:

к1 _ к10 ■ ехР

(

Е а

Л

ЯТ

(10)

Дифференциальное уравнение переноса энергии для жидкой фазы имеет вид:

ссШ$Т+с1Лт+щ£Лм_ (11)

г=1 Н Н Н

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

Количество тепла, выделяющееся в ходе реализации экзотермической химической реакции, определится выражением вида

V у

q хр _ жг ж у q хр^-

1 _1 с т

(12)

Дифференциальное уравнение переноса массы можно записать соотношением:

атсм _ЩзСт-Ь^сСт-Ущ^сС^ Опат. (13)

1_1 1_1 1_1

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

— _ -

пг

Изменение теплоемкости смеси определяется из соотношения:

¿Ссм = 2 ^Х; . 1=1

(14)

Концентрация 1-го компонента меняется за счет испарения, за счет химической реакции и за счет загрузки:

¿Х1 = ёХ1 исп + ¿Х1Х.р. + з .

(15)

Для реальных жидкостей изменение концентрации за счет испарения имеет вид:

=

ар

у1 ехр А1 -

-Х^ат-^ау/ (16)

Т У;

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

<т. Ь =-

тм,<&и

см гисп

(17)

¥ ■ <т ¥ ■ <т Скорость химической реакции численно может быть определена изменением весовой концентрации 1-го компонента в единицу времени :

ах

Ю: = -

1 х .р.

(18)

а т

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

ас п

— = к П ат 11 г

(19)

для гетерогенных систем

<С п

^ = К¥ П Сг

< т ¡=1

(20)

где К - коэффициент массопередачи, осложненный химическим превращением.

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

Р3 + ЪвР2 гвРв = V,

¿А, ат

- по пару

ЪпРз ^п А = V

ат

(21)

(22)

В уравнении (21) первый член левой части характеризует приток газовых компонентов; второй член - натекание воздуха через имеющиеся неплотности и технологические отверстия; третий - отвод

газовоздушной смеси системой откачки. В уравнении (22) первый член левой части характеризует приток паровых компонентов смеси; второй член -отвод образующихся паров системой откачки паров. Правая часть уравнений (21) и (22) представляет собой изменение в объеме аппарата массы газовоздушной и паровой смесей соответственно.

Уравнение теплового баланса паровой смеси определяется по соотношению:

рсмссХат=с^ Хъ к - ОкРсмссмтат.

V 1=1

(23)

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

Алгоритм расчета состоит из четырех блоков. В первом блоке вводятся теплофизические характеристики исходных веществ, константы уравнений, варьируемые переменные, величина временного шага и исходные данные для расчета первого шага, задаются условия расчета. Далее представлен цикл, в котором оператором перехода является наличие химической реакции Ю Ф 0 , вид химического разложения древесины.. При термическом разложении древесины (Т>Тр) расчет ведется на базе системы уравнения (1^10). При химическом разложении древесины (Т<Тр), за счет реагентов (например белого щелока), расчет ведется по соотношением (11^20) с учетом скоростей химических реакций: для гомогенных систем по уравнению (19) или гетерогенных систем по уравнению (20).

При отсутствии химических превращений или после завершения химических процессов (Юг = 0) ведется расчет паровой фазы в выдувном

резервуаре с использование упрощенных уравнений теплового и материального баланса (21^23), в которых отсутствуют члены, описывающие химические процессы.

При учете в паровой фазе неконденсирующихся газов Сг Ф 0 расчет ведется с введением III

блока с системой уравнений, включающей выражение (21).

В четвертом блоке осуществляется выбор системы газоочистки. При отсутствии в выбросах

газовой составляющей (Сг = 0) улов токсичных или ценных компонентов из паровой смеси осуществляется конденсацией.

При этом, если температура конденсации паров не превышает 300С, то требуется подготовка хладагента с применением холодильных устройств.

При улавливании ценных Ээфф>0 или токсичных паров Сг > ПДК1 целесообразно применять поверхностные конденсаторы.

Для улавливания парогазовых выбросов

(Сг Ф 0) необходимо устанавливать абсорбера, а для достижения более глубокой очистки при с1>ПДК1 рекомендуются адсорбционные устройства. При выполнении всех заданных условий происходит завершение работы программы.

газоочистки [5,6].

Литература

1. Сафин Р.Г. Технологические процессы и оборудование деревообрабатывающих производств: Учебное пособие. 2-е изд., испр. и доп. - М.: МГУЛ, 2002. - 688 с.

2. Зиатдинова, Д.Ф. Совершенствование технологий переработки древесных материалов, сопровождающихся выделением парогазовой фазы: Монография [Текст] / Д.Ф. Зиатдинова - Казань: КГТУ, 2009. - 140 с.

3. Зиатдинова, Д.Ф. Математическая модель технологических процессов, сопровождающихся локальными выбросами [Текст] / Д.Ф. Зиатдинова, Р.Г.Сафин, В.Н.Башкиров, Н.И.Горбачевский // Известия вузов: Химия и химическая технология. - 2004 г. - Т.47 - №10 - С. 129-132.

4. Зиатдинова, Д.Ф. Разработка технологии производства целлюлозы и побочных продуктов в виде спиртов, дрожжей, фурфурола из низкокачественной древесины [Текст] / Д.Ф. Зиатдинова // Вестник Московского государственного университета леса - Лесной Вестник. - 2010. - №4 (73). - С.84-87.

5. Зиатдинова, Д. Ф. Совершенствование системы газоочистки при сульфатной варке целлюлозы [Текст] / Д.Ф. Зиатдинова, Р.Г. Сафин, М.А. Мазохин, Р.Р. Зиатдинов, Д.А. Ахметова // Вестник Казанского технологического университета. - 2010. - №10 - С. 166-170.

6. Зиатдинова, Д.Ф. Моделирование процесса выгрузки варочного котла при переработки древесных отходов [Текст] / Д.Ф. Зиатдинова, Р.Г. Сафин, Д.Ш. Гайнуллина, М.А. Мазохин, Р.Р. Зиатдинов // Вестник Казанского технологического университета. - 2011. Т.14 - №18. -С. 76-80.

Рис. 2 - Блок-схема системы газоочистки по обобщенной математической модели

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

© Д. Ф. Зиатдинова - докторант каф. переработки древесных материалов КНИТУ, ziatdinova@rambler.ru.