Исследование кинетики гомогенных химических реакций

Исследование кинетических закономерностей протекания химической реакции методом математического моделирования заключается в определении изменения концентраций реагирующих веществ во времени при заданной температуре.

Пусть протекают химические реакции

k1 k2

 
 
A ® B +2C ® D.

На основании закона действующих масс запишем уравнения скоростей химических реакций и составим кинетическую модель:

; ; ; ; ,  

где СА, СВ, СС - концентрации веществ, моль/л;

ki - константа скорости i-й химической реакции первого порядка, с-1; (для реакций второго порядка размерность константы ;для реакций третьего порядка размерность константы ).

Wi - скорость i-й химической реакции, моль/л×с;

t - время реакции, с.

Систему обыкновенных дифференциальных уравнений первого порядка можно решить с использованием численного метода Эйлера, алгоритм которого записывается по уравнению (15).

Блок-схема расчета кинетики гомогенной химической реакции методом Эйлера приведена на рис.1.2.

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

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

Необходимо учитывать, что скорость химической реакции зависит от температуры, поэтому, чтобы использовать кинетическую модель (1.35) для исследования процесса при различных температурах, необходимо ввести зависимость константы скорости химической реакции от температуры по уравнению Аррениуса (1.9).

 

 
 

 

 


Рис. 1.2. Блок - схема расчета кинетики гомогенной химической

реакции методом Эйлера

 


Варианты заданий

Таблица 1.1

Варианты заданий

№ задания Уравнение химической реакции Начальные концентра- ции, моль/л Значение констант скоростей при Т=500 К Значения энергии активации Е, Дж/моль Значение предэкспоненциального множителя
k1 k2 k3 Е1 Е2 Е3 k01 k02 k03
  k1 2А«2В+С k2   CA0=0,5   0.2   0,15   -   9,305∙104   10,1∙104   -   1,376∙108   5,870∙108   -
  k1 2A®B+C k2 k3 D+E®B   CA0=0,3; CE0=0,2; CD0=0,2     0,2     0,4     -     11,514∙104     9,524∙104     -     1,686∙1010     4,214∙108     -
  k1 A«2B+C k2 k3 B+D®2C   CA0=0,3; CD0=0,4     0,2     0,2     0,1     12,444∙104     14,897∙104     14,1∙104     1,306∙1011     3,181∙1012     2,463∙1012
  k1 A+B«C k2 k3 2C+D®E   CA0=0,2; CB0=0,3; CD0=0,1       0,2     0,2     0,15     14,305∙104     12,594∙104     13,511∙104     7,714∙1012     1,853∙1011     9,949∙1011
  k1 2A+B«C k2   CA0=0,08       0,2   0,4   -   11,514∙104   10,386∙104   -   1,723∙1010   3,016∙109   -

Окончание табл.1.1

 

№ задания Уравнение химической реакции Начальные концентра- ции, моль/л Значение констант скоростей при Т=500 К Значение энергии активации,Е, Дж/моль Значение предэкспоненциального множителя
k1 k2 k3 Е1 Е2 Е3 k01 k02 k03
  k1 A+2C«B k2 k3 B®2D   CA0=0,2; CC0=0,1       0,42     0,2     0,25     11,354∙104     9,592∙104     9,431∙104     2,529∙1010     2,622∙108     2,260∙108
  k1 A®E+C k2 2C«D k3     CA0=0,3     0,3     0,25     0,18     14,18∙104     13,344∙104     13,762∙104     8,803∙1012     1,153∙1012     1,849∙1012
  k1 A«B+C k2 k3 2A®D     CA0=0,4     0,5     0,2     0,1     15,852∙104 15,016∙104       14,389∙104 5,693∙1013     4,100∙1012 4,747∙1012
k1 А«В k2 CA0=0,2; CB=0,05   0,37   0,16   -   12,717∙104     13,727∙104   -   4,266∙1011   1,643∙1012   -
k1 В А k2 С CA0=0,07; CВ0=СС0=0   0,4   0,2   -   10,351∙104   9,180∙104   -   2,424∙109   9,322∙108   -
k1 В А k2 k3 С CA0=0,06; CВ0=СС0=0   0,19   0,2   0,12   11,183∙104   12,103∙104   10,264∙104   7,888∙109   6,209∙1010   5,843∙108
k1 k2 А® В « С k3 CA0=0,2; CВ0=0,1; СС0=0   0,5   0,2   0,1   10,0∙104   9,8∙104   9,5∙104   1,572∙109   4,061∙108   1,054∙108

Порядок выполнения работы

1. Составить в соответствии с заданной схемой реакций на основании закона действующих масс кинетическую модель.

2. Разработать алгоритм расчёта составленной кинетической модели с использованием численного метода Эйлера.

3. Разработать программу расчета кинетики с учетом температурной зависимости констант скорости гомогенной химической реакции.

4. Обсудить результаты. Сделать выводы по работе.

5. Составить отчет.

 

ЛИТЕРАТУРА

 

1. Кафаров В.В. Методы кибернетики в химии и химической технологии. - М.: Химия,1988. - 489с.

2. Кравцов А.В., Новиков А.А., Коваль П.И. Методы анализа химико-технологических процессов. - Томск: изд-во. ТПУ,1994.-76с.

3. Кафаров В.В., Глебов М.В. Математическое моделирование основных процессов химических производств. - М.; Высш. шк.,1991.-400с.

4. Мойзес О.Е., Коваль П.И., Баженов Д.А., Кузьменко Е.А. Информатика / Учеб. пособие. В 2-х ч.- Томск: 1999, - 150 с.

5. Турчак Л.И.Основы численных методов. – М: Наука, 1987. -320 с.

 



>2