Математическая модель определяется лимитирующей стадией процесса. В случае, когда лимитирующей стадией является химическая реакция, математическая модель будет включать в себя дифференциальные уравнения первого порядка, описывающие изменение концентрации каждого вещества во времени:
где [x] – концентрация интересующего нас вещества,
t – время,
rx – скорость изменения концентрации данного вещества во времени, имеющая, следующий вид:
, где - предэкспоненциальный множитель, энергия активации процесса,R – универсальная газовая постоянная,
С – концентрации реагентов
m, n – порядок реакции по веществам a и b соответственно. В данной работе все порядки приняты равными единице.
Если же лимитирующей стадией процесса является внутренняя или внешняя диффузия, математическая модель будет состоять из дифференциальных уравнений второго порядка. В эти уравнения входят также такие параметры, как скорость движения потока, размер пор катализатора, коэффициент диффузии. Такие зависимости позволяют определить концентрацию вещества в зависимости от длины реактора и расстояния от его оси. Вид этих уравнений представлен ниже.
Для нахождения концентрации вещества в потоке:
Для нахождения концентрации вещества на поверхности и внутри катализатора:
Скорость реакции разложения оксида азота описывается при помощи кинетической модели по механизму Лэнгмюра-Хиншельвуда:
Перечень символов:
- концентрация в потоке; - концентрация в твердой фазе (на поверхности катализатора и внутри его пор; - коэффициент эффективной диффузии;r – радиальная координата;
- коэффициент молекулярной диффузии;k – константа скорости
ra – скорость реакции;
Z – аксиальная координата.
U – скорость газового потока;
1.5. Лимитирующие стадии гетерогенного каталитического процесса.
Как известно, у гетерогенной каталитической реакции может быть несколько лимитирующих стадий. Лимитирующей стадией может являться:
Внешняя диффузия реагентов или продуктов реакции. В этом случае скорость процесса будет определяться скоростью, с которой частицы реагентов из ядра газового потока будут попадать на поверхность катализатора или тем, насколько быстро образовавшиеся частицы будут уходить с поверхности катализатора в ядро газового потока.
Внутренняя диффузия реагентов или продуктов реакции. В этом случае скорость процесса определяется тем, насколько быстро молекулы реагента проникают в поры катализатора или же тем, как быстро молекулы продукта реакции освобождают пространство пор.
Адсорбция или же десорбция на поверхности катализатора. В такой ситуации скорость процесса будет определяться тем, насколько быстро происходит насыщение активных центров катализатора молекулами веществ или насколько быстро десорбируются полученные вещества с поверхности катализатора.
Непосредственно химическая реакция – в этом случае скорость процесса определяется скоростью взаимодействия веществ активных центрах катализатора.
В литературе существует ряд работ, показывающих, что большинство реальных химических реакций являются сложными и идут через посредство активных промежуточных продуктов. Реакция между двумя стабильными молекулами требует большой энергии активации. Поэтому скорость такой реакции мала. Реально наблюдаемые нами реакции идут обычно "окольными" путями, позволяющими обойти этот высокий энергетический барьер.
В случае гомогенных реакций первой стадией является обычно образование из исходных веществ тех или иных активных промежуточных продуктов. Природа этих продуктов в настоящее время еще недостаточно хорошо изучена, но, надо полагать, она может быть различной. Во многих случаях роль активных промежуточных продуктов, несомненно, играют свободные радикалы или даже свободные атомы. В других случаях это могут быть и довольно сложные и относительно стабильные молекулы, обладающие, однако, по тем или иным причинам, повышенной реакционной способностью, например, органические перекиси.
Активные промежуточные продукты, реагируя с исходными веществами, переводят их в конечные продукты реакции. Эти процессы требуют сравнительно малой энергии активации (особенно когда активными продуктами являются свободные радикалы или атомы) и идут с большей скоростью. Но первичное образование активных продуктов из стабильных исходных молекул требует большой энергии активации, и поэтому с большой скоростью проходить не может.
Для того чтобы реакция посредством активных продуктов могла протекать достаточно быстро, необходимо, чтобы активные продукты регенерировались при реакции, т.е. чтобы при взаимодействии активных продуктов с исходными веществами получались не только стабильные конечные продукты реакции, но и новые молекулы активных промежуточных продуктов. Такие реакции, в которых имеет место регенерация активных промежуточных продуктов, называются цепными(3, стр 17) Как было установлено Семеновым, Хиншельвудом и др., большинство реальных гомогенных сложных реакций являются цепными.
Обозначим концентрацию активного продукта как x. Изменение этой величины со временем для цепных реакций подчиняется следующему кинетическому уравнению:
где no – скорость зарождения цепей;
f – константа скорости процесса разветвления цепей (в общем случае не совпадает с общепринятым понятием константы скорости элементарной реакции
g - константа скорости обрыва цепей.
Под зарождением цепей подразумевают начальный процесс образования продукта из исходных веществ, под разветвлением цепей – процесс, в котором одна молекула активного продукта, реагируя с исходными веществами, вызывает образование двух или нескольких молекул активного продукта; под обрывом цепей – процесс, при котором активный продукт безвозвратно уничтожается.
Актуальность проблем очистки газовых выбросов от оксидов азота подтверждается значительным количеством тематических публикаций. Из литературных данных видно, что интерес исследователей к этой проблеме не угасает. Наиболее перспективными каталитическими методами удаления оксидов азота являются процессы селективного каталитического восстановления NOx с использованием в качестве восстановителей аммиака и углеводородов.
Проблема очистки газовых выбросов от NOх газов дизельных двигателей, в которых сгорание производится на “бедных” смесях, содержащих значительный избыток кислорода до сих пор не решена, т. к. в таких смесях оксиды азота обычно подвергаются окислению, а не восстановлению.
Наиболее эффективными катализаторами являются катализаторы блочного типа. Однако процессы, происходящие на таких катализаторах, по-прежнему изучены слабо. Математическое моделирование дает возможность более детально рассмотреть механизм данного процесса, оценить влияние различных стадий, а также параметров химического процесса на итоговую конверсию оксида азота, сделать предположения об условиях, в которых целесообразно будет проведение процесса.
2. ПОСТАНОВКА ЗАДАЧИ И ПОСТРОЕНИЕ МОДЕЛИ.
2.1. Экспериментальные данные.
Рис 3. Зависимость конверсии оксида азота от площади катализатора.
Учитывая наше допущение о том, что процесс восстановления оксида азота является радикальным, логично предположить, что: