CC BY
30
УДК 536.27
Барышева О.Б. - кандидат технических наук, доцент
E-mail: [email protected]
Хабибуллин Ю.Х. - кандидат технических наук, доцент
E-mail: [email protected]
Вагизов А.З. - студент
E-mail: 3 almaz@mail. ru
Казанский государственный архитектурно-строительный университет
Адрес организации: 420043, Россия, г. Казань, ул. Зелёная, д. 1
Математическая модель расчета химически неравновесного состава для утилизации твердых бытовых отходов
Аннотация
Постановка задачи. Для большинства агрегатов, в которых процессы протекают при высокой температуре, свойственны различного рода неравновесные эффекты. Важнейшим из них является неравновесное изменение состава рабочего тела продуктов сгорания, нагрева или разложения, вызванное конечными скоростями протекающих химических реакций, приводящих к существенному влиянию на основные энергетические параметры процесса. Авторами было принято решение создать математическую модель расчета химически неравновесного состава для утилизации твердых бытовых отходов.
Результаты. Авторами рассмотрен процесс утилизации твердых бытовых и промышленных отходов. На основе имеющихся методов расчета разработан альтернативный метод, позволяющий находить состав продуктов сгорания в условиях неравновесности в процессе сгорания строительного мусора. Полученные разработанным методом расчета данные состава и свойств ТБО позволяют прогнозировать выход суперэкотоксикантов в продуктах сгорания.
Выводы. На основе уравнений формальной химической кинетики мы создали альтернативную методику расчета состава продуктов сгорания с учетом скоростей химических реакций. Данная методика заключаются в определении состава рабочего тела по заданному механизму и константам скоростей реакций. Принимается, что все превращения в газовой фазе элементарны, то есть протекают в одну стадию. Любые химические взаимодействия могут быть представлены совокупностью этих элементарных стадий. Рассматривались наиболее вероятны химические реакции трех типов: моно-, би-и тримолекулярные. Данная методика позволяет найти состав продуктов сгорания твердых бытовых отходов еще до процесса сжигания топлива.
Ключевые слова: неравновесные процессы, химическая кинетика, равновесие, состав, твердые бытовые отходы.
Диоксины на частицах дыма, образующегося при сжигании ТБО, обнаружены почти четверть века назад. Однако, несмотря на многочисленные лабораторные и натурные исследования в настоящее время отсутствует полное количественное описание процессов образования и распада полихлорированных диосинов, фуранов, бефинилов. Концентрации диоксинов, образующихся при сжигании ТБО обычно малы, а публикуемые результаты исследований имеют широкий разброс в измеренных концентрациях [1, 2]. Поскольку образование диоксинов является лишь незначительной побочной реакцией в общем процессе распада и синтеза органических соединений, интерпретация экспериментальных результатов является гораздо более сложной процедурой по сравнению с исследованиями основного реакционного процесса.
Важным результатом анализа состава диоксинов на частицах дыма и направляемых в дымовую трубу газов было установление распределения атомов хлора в дибензо-р-диоксинах. В большинстве случаев предпочтительнее образуются более сильно хлорированные соединения, и часто максимум распределения соответствует шести замещений атомами хлора атомов водорода [3, 4].
В некоторых случаях предпочтительными оказываются октохлорзамещенные соединения. Изомер 2,3,7,8-тетрахлордибензо-^-диоксин образуется не в самом большом количестве. Пока нет общепринятой теоретической модели механизма образования диоксинов, воспроизводящей наблюдаемое распределение. Это связано в частности со сложным набором физико-химических процессов, протекающих в установках для сжигания ТБО [5, 6].
Обобщенно процесс преобразования ТБО в конечные продукты можно представить следующим образом. Отходы загружаются различными способами на колосниковую решетку, где они продуваются горячим воздухом. Это ведет к выпариванию воды и пиролитическому разложению органики при низкой температуре в обедненной кислородом области вблизи решетки, что способствует дополнительному образованию воды и низкомолекулярных органических соединений, которые выносятся в газовую фазу [7]. Происходит значительная олигомеризация продуктов, некоторые из которых могут остаться в конденсированной фазе и продолжать пиролизоваться. При контакте органических соединений с кислородом при достаточно высоких температурах начинается процесс горения, который продолжается до конца расходования кислорода. Некоторые из твердых и жидких веществ, находящихся на решетке или в виде мелких частиц, также могут гореть, а те, которые не горят, в результате дальнейших пиролитических реакций преобразуются во все более ненасыщенные структуры. Процесс полимеризации в конечном итоге приводит к образованию углерода. Чтобы минимизировать такие пиролитические процессы при сжигании, естественным является желание осуществлять весь процесс горения при избытке кислорода, а это может способствовать большей эмиссии диоксинов. Непосредственно в камере сгорания, если образуются диоксины, то при прохождении через высокотемпературные зоны происходит их разложение [8]. Поэтому концентрация диоксинов в камере сгорания незначительна. Результаты взятия проб указывают на образование все более хлорированных веществ в направлении от выхода из камеры сгорания к выпускной трубе. Роль самой камеры сгорания состоит в поставке молекул - предвестников, образования диоксинов. Отсюда важным является определение оптимальных условий для минимизации процесса образования этих предвестников, что может ограничить и процесс образования диоксинов.
Анализ литературных источников показывает [9, 10], что образование диоксинов происходит в диапазоне температур 500-800 К. В присутствии частиц дыма, хлора и при избытке кислорода многие ароматические соединения могут образовывать соединения типа диоксинов с соответствующим распределением по семейству хлорзамещенных диоксинов. Наиболее вероятными предвестниками могут быть хлорированные фенолы и высокохлорированные ароматические эфиры. Эксперименты показывают, что образование диоксинов легче всего происходит при наличии пентохлорфенола (С6С15ОН).
С помощью объектно-ориентированного языка программирования С++ нами создана программа расчета состава продуктов сгорания топлива для локального термодинамического равновесия. Однако, процессы, протекающие в мусоросжигательных печах, носят неравновесный характер, поэтому прогнозировать выход диоксинов и суперэкотоксикантов практически невозможно. Можно лишь обнаружить предвестников диоксинов. Поэтому наша цель запрограммировать метод расчета продуктов сгорания с учетом неравновесности.
Рассмотрим метод расчета состава продуктов сгорания ТБО в условиях неравновесности:
Е?= IV 'и Ь 7 = 1.....т, (1)
где В7 (7=1,...,п) - компонент (реагирующая система состоит из п - компонентов, между которыми протекают обратимых химических реакций, то есть протекающих как в прямом, так и в обратном направлении); V ¿у, V¿у - стехиометрические коэффициенты}-й реакции.
Основным законом формальной химической кинетики является закон действующих масс, согласно которому скорость изменения содержания 7-го компонента по прямому направлению}-й реакции есть соотношение:
Щ = М - пЦ) */ш=Ч (2)
Аналогично для обратного направления:
щ = У- п'¡) /с/ге=1 < (з)
где /у",/7" - константы скорости}-й реакции в прямом и обратном направлении.
В формальной химической кинетики константы скоростей как прямых, так и обратных химических реакций обычно принято представлять в аррениусовой форме, так для прямого направления:
к;=А} Т^ехр (-$) , (4)
где А;, п+, Е? - коэффициенты характеризующие прямое направление }-й реакции: предэкспоненциальный множитель, температурный коэффициент и энергию активации. В зависимости (4) размерность определяется порядком реакции, - безразмерная величина, Е^ - имеет размерность обычно Дж/моль. В этом случае = 8,314 Дж
-------- ^-----^-------------^.......... ~ ^— ^ — ^ „ {мольК)
универсальная газовая постоянная. Коэффициенты А}, п^, Е ^ - определяются эмпирически.
Константы скорости прямого и обратного направления реакций связаны через константы равновесия последних:
Кр} = Щ/ к/. (5)
Константа равновесия, являясь функцией температуры, определяет вместе с тем отношение равновесных значений концентраций реагирующих компонентов, т.е. концентраций, при которых = . В этом случае:
кР1 = П * С/п * С. (6)
Константу равновесия целесообразно представить в виде комбинации констант диссоциации, выраженных через парциальные давления:
о К} = П /с (Ф^/П /с (ФУы Я0Т£к(пк] - Ук]), (7)
где К к = рк/ПвР^к, а*к - стехиометрические коэффициенты в реакции диссоциации вещества на атомы.
В формальной химической кинетике предусматривается принцип независимости протекания элементарных химических реакций в реагирующей смеси. То есть, общая скорость изменения любого компонента в смеси равна алгебраической сумме скоростей его изменения в результате протекания всех элементарных реакций, идущих при постоянном объеме:
^ = щ = !?= 1 ^ + £7= 1 щ, (8)
где 7=1,..., п,}=1,., т, п - число индивидуальных веществ; т - количество химических реакций.
С использованием формул (2) и (3) формулу (8) можно записать в следующем виде:
% = г(У'- - уЦ)К П С? + цр=,(У, - У-'/)к- Гк Ср. (9)
Химическая реакция характеризуется порядком, который равен сумме показателей степеней при всех концентрациях, стоящих в правой части кинетического уравнения.
Метод решения уравнений химической кинетики
Рассматривается реагирующая газовая система, в которой между компонентами .4¡(7=1,., п) протекают обратимые элементарные химические реакции (5=1,., т).
Каждая обратимая химическая реакция рассматривается как совокупность двух односторонних. Введем следующие обозначения:
Уа = К - К, па = К,] = 5> 5 = 1'-'^ (ю)
УЧ = У18 - пч = ] = 5+ 5 = 1,-< т
В соответствии с принятыми обозначениями прямые направления химических реакций имеют индексы }=1,., т обратные направления-индексы }=(т+1),., 2т.
Применяя соотношение С1 = п¿/ V и =п¿/Л = С ¿/С = р^р можно уравнение (9) переписать в следующем виде:
= ^уЙШрСО^, (11)
где 7, р=1,.,п, }=1,.,2т, С - мольная концентрация газовой смеси; Ср - мольная парциальная концентрация р -го компонента газовой смеси; т}- - признак участия в у-ой реакции каталитической частицы М(т}=1), реакция протекает с участием этой частицы, в противном случае т}=0); п7 - количество молей 7-го компонента в объеме V, к - константа скорости обратимых химических реакций.
Используя соотношения п7=Ши С7=рг7/(Я0Т), где N - общее число молей в объеме V, получим:
п (12)
] р (12)
=V пРг; - = 1.....-
Число уравнений (12) равно числу компонент газовой смеси. Суммируя по 7 все эти уравнения, получим:
V п с1м + 2 N а п = гам , N а(у 1 гр = 2 2 п ( Р \т]п пр] (13)
1 + 2' ~ = + = 2'пргР , (13)
где ^=1,., п.
Введя обозначение, т^ = ту + V пр ] — 1 и, используя равенство N/V=p/(R0T) получим:
5=V ПА (¿Г Пр ^—„ V 2Р ^ (¿г П г;«. (14)
где 7, р, 9=1,., п, }=1,., 2т.
Соотношение (14) является нелинейным дифференциальным уравнением, которое интегрируется с помощью численных методов, и в этом случае для достаточно малых концентраций г7 ошибки округления могут привести к получению физически нереальных отрицательных концентраций г7. Для предотвращения этого вводится следующая замена:
7г = -1П Г.
Тогда уравнение изменения состава примет следующий вид: 5 = е т' ехр(— ^ пр]Гр) + + 2 ч ^ ^ ехр(- 2 р пр]Гр). (15)
где 7, р, 9=1,., п, }=1,., 2т.
Одно из уравнений (15) можно заменить уравнением нормировки:
2^-^ = 1. (16)
Заменяется уравнение, которое описывает изменение компонента, содержащегося в газовой смеси в максимальной концентрации.
Изменение температуры рабочего тела влияет на его термодинамические и теплофизические свойства и на скорости химических реакций. Поэтому требуется в математическую модель включить соотношение, которое связывает состав рабочего тела с температурой. Это делается на основе формулы для энтальпии смеси:
¿ = 2 ькП^И ¿П. (17)
где 17, т - значения мольной энтальпии и молекулярная масса 7-го компонента.
Значения мольной энтальпии, энтропии и теплоемкость индивидуальных веществ смеси аппроксимируется полиномами базы ИВТАНТЕРМО, отличающихся тем, что рабочий диапазон 298,15.20000 К разбит на шесть температурных диапазонов, каждый из которых описан полиномами с уникальными коэффициентами (а _ 2, а _ а 0, а2, а 3):
1(Т) = [Н(Т) - Я(0)]г + ДЯ/0 г • 10 3, [Н(Т) - Н(0)] = 104[-2а_ 2х - 1 + а_х + аъх + ахх 2 + 2а2х3 + 3 а3 х 4 ], (Г) = -а _2 х" 2 + аь 1п х + а 0 + аь + а 1х2 + 2 а 1 х + 3 а 2 х2 + 4а3х3, ср1 = - 2а_2х~2 + аь + 2а1х + 6а2х2 + 12а3х3, х = Г/10000.
В ходе выполнения расчетов приходится очень часто обращаться к этим соотношениям, поэтому для уменьшения объема вычислений их целесообразно аппроксимировать кусочно-линейными зависимостями типа:
/ = Ь оп + - Гоп), (18)
где / оп - «опорное» значение мольной энтальпии 7-го компонента, Топ - опорное
оп
значение температуры; - опорная мольная теплоемкость, равная среднему значению
на интервале линеаризации. С достаточной для решения прикладных задач точностью этот интервал АТ можно принять 100 К.
При использовании (18) соотношение (17) принимает следующий вид:
i _ & (' °п + с°рП(Т - Топ)) П /(19)
После простых преобразований получим:
т - т°п - & (¿^ - 1°п) / & с°рпп _ 0. (20)
Система уравнений, использованная для определения состава продуктов сгорания твердых бытовых отходов, включает в себя уравнения (15), (16), (20):
^ = -е7; & п]к] ' ехр(- & пр]Ур) + + & & ' ехр(- & пр]Гр),
& ¿е" ( = 1, (21)
1_1 _& / (щ / ~ ^ о)_о
1 1 °п & г°пг. ~ ° Ь1Ср1'1
где 1, р, 4=1,..., п,7=1,., 2га.
Неизвестными в системе (21) являются п-мольных концентраций индивидуальных веществ, составляющих газовую смесь и значение температуры Т. Для индивидуальных веществ составляющих газовую смесь и для всей смеси в целом справедливо уравнение состояния идеального газа:
рУ=Ъ>Т. (22)
Система (21) представляет собой (п-1) обыкновенных дифференциальных уравнений и 1 алгебраическое уравнение.
Специфика решения системы (21) заключается в «жесткости» уравнений изменения состава (15). Получение решений для таких уравнений возможно только при использовании неявных разностных методов. Эффективным неявным методом конечно-разностной аппроксимации дифференциальных уравнений является метод У.Г. Пирумова с последующим использованием метода Ньютона.
Для применения метода Ньютона дифференциальные уравнения предварительно представляются в виде однородных конечно-разностных уравнений вида:
^ = 1 - хъ - Ы№» + (1 - а) к да !»] = о, (23)
где I, &=1,..., 4, 4 - число неизвестных; И - шаг интегрирования; а - параметр аппроксимации (обычно а=0,4); Х£3, ХЦ* 1 - значения неизвестной величины в начале и конце шага интегрирования.
Полученная система представляет собой нелинейные алгебраические уравнения, которые преобразовываются для реализации схемы Ньютона к виду:
8 [-Э ищ _ ил (24)
где [5 Р 1 1 - матрица частных производных (матрица Якоби).
\-дХк5\
Система вида (24) является системой линейных уравнений относительно поправок А , которая решается методом Гаусса. Для нахождения неизвестных применяется итерационная схема:
Х^п+1_хп+1, га + д х пп(25)
где га - номер итерации; АХ™+1т - поправки на т-ой итерации; п+1 - номер шага интегрирования, на котором реализуется схема; 1=1,., д, ц- число неизвестных.
На шаге интегрирования итерации выполняются до получения таких значений ( ), которые отличаются от своих предыдущих значений на заранее заданную
малую величину. Затем осуществляется переход на следующий шаг интегрирования. При этом в качестве схемы интегрирования принята схема с «замороженным» на ряде шагов якобианом, что существенно сокращает время вычислительного процесса.
Наиболее трудоемкой частью алгоритма является вычисление матрицы Якоби. Применение для нее аналитического вычисления позволяет сохранять «замороженный» якобиан на большем числе шагов интегрирования.
Особые требования предъявляются к уравнениям изменения состава: они должны быть вычислены в обобщенном виде. Это позволит создать программу расчета независимую от набора реагирующих веществ и совокупности химических реакций.
Для системы (21), приведенная к конечноразностной форме матрица Якоби имеет
следующие элементы:
dF
= S? + М1 - а)е" Lj Vij(S? - nkj)Wj - Lj vpjnkjО,.,
i = 1,..., (L - 1), (L + 1).....n,
p,k = 1,..., n, j = 1,... ,2m,
где S f - символ Кронекера,
(26)
à-Yi
-¡^ = M1 - a)en^vi}
dF
= exp(-Ki), (27)
Vy W,
1 dAj n.j - m. Ej
Aj дТ T R0T2
г г ГЛ { 1 dAj nj-mj Ej N
- Lv L i VviWi I--+ —-- +--I
^PPJ J \Aj дт T R0T2J
W = k (Îï) ' eXP(- L "pJYPY
(28)
В алгоритме расчета предусмотрено вычисление констант скоростей обратных химических реакций, причем константы скоростей прямых реакций относятся к исходным данным, а константы равновесия находятся как комбинации констант диссоциации.
Константы диссоциации определяются по аппроксимирующим коэффициентам, которые позволяют представлять эти константы в виде усеченной зависимости вида Аррениусовой:
KP=A?exp(-EP/(R 0Т)). (29)
Для сокращения времени расчета весь рабочий температурный интервал (300.4200 К) разбит 20 «опорными» точками на равные температурные интервалы. В этих опорных точках для каждого индивидуального вещества, принимающего участие в расчетах, предварительно рассчитываются и записываются в соответствующий массив «опорные» значения констант диссоциации:
ln K0Dn = (Li aimS° - 5° )/ R0 - (Li a imI? - О/ R0T, (30)
где щт - число атомов «i» в молекуле диссоциирущего вещества «m».
По опорным значениям констант диссоциации рассчитываются аппроксимирующие коэффициенты для определения констант диссоциации на правых и левых границах температурных интервалов:
Г1=Гоп+100 К, Т1=Топ-100 К, (À+D, A-d, E+d, Ed). Этот метод позволяет спрогнозировать выход диоксинов и суперэкотоксикантов еще до самого процесса утилизации ТБО.
Список библиографических ссылок
1. Подавалов Ю. А. Экология нефтегазового производства. М.: Инфра-Инженерия, 2010. 416 с.
2. Савельева В. А., Попова М. Н., Андронова М. М. Экологическая безопасность учреждений УИС в сфере обращения с отходами. М. : ВИПЭ ФСИН России, 2010. 250 с.
3. Бычкова В. А. Коробко В. И. Твердые бытовые отходы. Экономика. Экология. Предпринимательство. М. : Юнити-Дана, 2014. 132 с.
4. Муравьев А. Г., Родин А. А. Другов Ю. С. Экспресс-анализ экологических проб. М. : Бином. Лаборатория знаний, 2010. 424 с.
5. Почекаева Е. И. Экология и безопасность жизнедеятельности. М. : Феникс, 2010. 560 с.
6. Porter M. What is strategy? Harvard : Business Review, 2009. 108 p.
7. Mill R. S., Morrison A. M. Nature and Ecology. Boston : Prentice-Hall, 2009. 230 p.
8. Odum Yu. Ecology. Cambridge : Second Edition, 2008. 256 p.
9. Odum Yu. Ecology. Cambridge : Second Edition, 2010. 338 p.
10. Chapman J. L. Ecology. Cambridge : Second Edition, 2009. 336 p.
Barysheva O.B. - candidate of technical sciences, associate professor
E-mail: [email protected]
Khabibullin Iu.Kh. - candidate of technical sciences, associate professor
E-mail: [email protected]
Vagizov A.Z. - student
E-mail: 3 almaz@mail. ru
Kazan State University of Architecture and Engineering
The organization address: 420043, Russia, Kazan, Zelenaya str., 1
A mathematical model for calculating the chemical nonequilibrium components
for solid waste management
Abstract
Problem statement. For the majority of units in which processes occur at high temperatures characteristic of various types of non-equilibrium effects. The most important of these is the change in non-equilibrium of the working fluid composition of the products of combustion, heating or decomposition caused by finite speed the chemical reactions that lead to a significant impact on the basic energy parameters of the process. The authors decided to create a mathematical model for calculating the chemical nonequilibrium components for solid waste management.
Results. The authors reviewed the process of solid waste and industrial waste. On the basis of calculation methods developed an alternative method to find the composition of the combustion products in the non-equilibrium conditions during the combustion of construction debris. The resulting method of calculation of the data structure and properties of solid waste can be predictive superecotoxicants access to the products of combustion.
Conclusions. On the basis of the formal chemical kinetics equations, we have created an alternative methodology for calculating the composition of the combustion products, taking into account the rates of chemical reactions. This technique consists in the definition of the working fluid composition for a given mechanism and rate constants of reactions. It is assumed that all of the conversions in the gas phase - elementary, i.e. occur in one step. Any chemical interactions may be represented by a set of elementary stages. They were considered the most likely chemical reactions of three types: a mono-, bi- and trimolecular. This technique allows you to find the composition of the products of combustion of municipal solid waste before combustion.
Keywords: nonequilibrium processes, chemical kinetics, balance, composition, municipal solid waste.
References
1. Podavalov Ju. A. Oil and Gas Production Ecology. M. : Infra-Engineering, 2010. 416 p.
2. Saveleva V. A., Popova M. N., Andronova M. M. Environmental security penal institutions in the field of waste management. M. : VIPE FSIN of Russia, 2010. 250 p.
3. Bychkova V. A., Korobko V. I. Municipal solid waste. Economy. Ecology. Entrepreneurship. M. : Unity-Dana, 2014. 132 p.
4. Muravev A. G., Rodin A. A., Drugov Ju. S. Rapid analysis of environmental samples. M. : Binom. Laboratory znanij, 2010. 424 p.
5. Pochekaeva E. I. Ecology and life safety. M. : Fenix, 2010. 560 p.
6. Porter M. What is strategy? Harvard : Business Review, 2009. 108 p.
7. Mill R. S., Morrison A. M. Nature and Ecology. Boston : Prentice-Hall, 2009. 230 p.
8. Odum Yu. Ecology. Cambridge : Second Edition, 2008. 256 p.
9. Odum Yu. Ecology. Cambridge : Second Edition, 2010. 338 p.
10. Chapman J. L. Ecology. Cambridge : Second Edition, 2009. 336 p.
