Скачиваний:
692
Добавлен:
07.01.2014
Размер:
207.36 Кб
Скачать

Российский химико-технологический университет им. Д.И.Менделеева

кафедра общей химической технологии

ОТЧЁТ

по лабораторной работе:

" Окисление SO2 в SO3 "

выполнили : Благодёр А.

Косарев А.

Берестов В.

Барыкин С.

Субботин Ю.

группа : O-44

МОСКВА -2005

1. Цель работы

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

Анализ и выбор оптимальных условий функционирования химико-технологиче-

ской системы.

Исследование включает следующие этапы:

Определение равновесной степени превращения.

Определение скоpости pеакции W(T) пpи заданном X.

Определение оптимальной темпеpатуpы Tоpt (X)

Определение вpемени контакта в слое t (X) (pежимы ИС и ИВ)

2. Общие сведения

Скорость реакции

SO2 + 0,5O2 = SO3 (2.1)

описывается уравнением Борескова-Иванова :

, (2.2)

где - парциальные давления SO2, SO3 и O2;

k, KP - константы скорости реакции и равновесия.

Используя степень превращения SO2 - x, уравнение (2.2) преобразуется:

, (2.3)

где a, b - начальные концентрации SO2 и O2, об.доли.

Р - общее давление;

Константа равновесия KP зависит от температуры T [K]:

lg KP = 4905/T - 4,6455 (2.4)

Зависимость констаны скорости реакции k от T в аррениусовских координатах "lnK - 1/T" приведен на рис.2.1.

В реакционной среде фазовое состояние и химический состав активного компонента ванадиевого катализатора изменяется с температурой. Поэтому зависимость k(T) будет различна в разных температурных областях. На аррениусовском графике “ln k - 1/T” это представлено ломаной линией (рис.2.1), каждый из участков которой характеризуется своим значением энергии активации (E1, E2, E3)

Рис. 2.1. Зависимость константы скорости k реакции окисления SO2 на ванадиевом катализаторе от температуры Т.

Пояснения в тексте

Алгоритм расчета. При значениях температур от Тн шагом Т по формуле (2.3) рассчитываются значения W(х). Если температура превышает Тр ( при этом W(х) < 0), то расчет прекращается.

Pавновесные степени пpевpащения xp(T)

Равновесие в реакции (2.1) определяется уравнением

(2.5)

Здесь ()равн - парциальные давления компонентов в равновесии

Перейдя к равновесной степени превращения SO2 xр, получим

(2.6)

Алгоритм расчета. Для каждого значения температуры Т из (2.4) рассчитывют КP и решают уравнение (2.6) относительно xр методом половинного деления. Первоначальный интервал поиска xр 0  1 или 0 2b/a, если 2b/a < 1. Точность расчета xр - 10-5

Оптимальные темпеpатуpы Tопт(х)

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

dW(х,T)/dT = 0 (2.7)

Представим W(х,T) и входящие в нее КP(Т) и k(Т) в таком виде, как следует из (2.3)-(2.5):

W(x,T) = К f(x) [1 - (x)/]

; КP = К exp(Q/RT); k = kо exp(-E/RT)

Из условия (2.7) определения оптимальной температуры получим

(2.8)

Алгоритм расчета. При значениях степеней превращения от xн с шагом x до xк по формуле (2.8) расcчитываются значения Топт

Вpемя контакта в слое (х)(pежимы идеального смешения и вытеснения)

Математически процесс описывается уравнениями:

в слое идеального вытеснения: dx/d = W(x,T), при  = 0: x = xн (2.9)

в слое идеального смешения: (x - xн)/ = W(x,T) (2.10)

В адиабатическом режиме :

T = Tн + Tад(x - xн), (2.11)

В изотермическом режиме T = Const, что получим, задавая T = Тн и Tад = 0.

Здесь Tад = Qрa/Cр - величина адиабатического разогрева;

Qр - тепловой эффект реакции;

Ср - теплоемкость реакционной смеси.

Из (2.9) и (2.10) следуют расчетные уравнения (2.12) и (2.13) соответственно:

(2.12)

 = (x - xн)/W(x,T) (2.13)

Алгоритм расчета. Значение  для слоя идеального смешения вычисляется непосредственно по формуле (2.13). Интеграл из (2.12) вычисляется методом Симпсона с постоянным шагом х. Если f(х) = 1/W(х), то приращение интеграла i наi-м шаге

i = x/6· [f(xi) + 4f(xi + x/2) + f(xi +x)]

Принятый шаг интегрирования x = (xк - xн)/200 обеспечивает достаточную точность интегрирования (не хуже 0,1%отн.).

3.Исходные данные

Состав газа: 10% SO2,14% O2, остальное-азот.

Расход газа: 20 л/час

Степень конверсии SO2 : 0.987

ΔΤad=279 К

Катализатор: Classic

Давление : 1 атм.

4.Расчётная часть

Рассчёт равновеcных степеней превращения в интервале от 680 до 850 К с шагом 17 К.

T,K

X*

Kp

680

0,993

369,6

697

0,9895

246,49

714

0,9845

167,59

731

0,978

116,01

748

0,9691

81,66

765

0,9571

58,38

782

0,9423

42,35

799

0,9234

31,15

816

0,9001

23,2

833

0,8721

17,49

850

0,8394

13,34

График 4.1 X* = f(T)

График 4.2 Kp = f(T)

Расчет скорости реакции при различных степенях конверсии(Х=0,2÷0,95)

T,K

W при Х =

0,5

0,6

0,7

0,8

0,9

0,95

680

0,427

0,340

0,254

0,169

0,083

0,039

697

0,818

0,652

0,487

0,322

0,156

0,068

714

1,065

0,848

0,632

0,416

0,194

0,068

731

1,366

1,087

0,807

0,524

0,223

0,028

748

1,729

1,372

1,011

0,638

0,214

-0,118

765

2,158

1,702

1,234

0,734

0,098

-

782

2,211

1,723

1,210

0,628

-0,0238

-

799

2,175

1,658

1,088

0,384

-

-

816

2,111

1,540

0,869

-0,056

-

-

833

1,998

1,332

0,483

-

-

-

850

1,805

0,977

-0,176

-

-

-

Построение ЛОТ

Конверсия

Торт

0,02

1355,5

0,12

1102,7

0,22

1028,6

0,32

981,7

0,42

945,3

0,52

913,6

0,62

880,4

0,72

852,2

0,82

816,1

0,92

764,1

График 4.3 ЛОТ

Модель идеального смешения

№ слоя

Xн

Xk

τad

1

0

0,6783

0,243

2

0,6783

0,8652

0,234

3

0,8652

0,9434

0,416

4

0,9434

0,9748

0,675

5

0,9748

0,9870

0,865

Σ τad=2,453

Расчёт объёма контактной массы(ν):

V=20 л/час/ 3600 с/час=0,005556 л/с

νi=V· τadi

№ слоя

τad,c

νi, л

1

0,243

0,00135

2

0,234

0,0013

3

0,416

0,002311

4

0,675

0,00375

5

0,865

0,004806

Σ τad=2,453

Σ νi=0,013629

Вывод:

1)Понижение температуры повышает равновесную степень превращения и конс-танту равновесия.

2)Оптимальная температура снижается с ростом конверсии.

3) Суммарное время контакта растёт с увеличением степени превращения.

Суммарное время контакта 20 л/час смеси ,состоящей из 10% диоксида серы , 18,9% кислорода и 71,1% азота при конверсии 0,987 составило 2,453 с, числом слоёв катализатора Classic- 5, общим объёмом - 0,013629 литра

8