МОДЕЛИРОВАНИЕ ХИМИКО-ТЕХНОЛОГИЧЕСКИХ ПРОЦЕССОВ КАК ОБЪЕКТОВ УПРАВЛЕНИЯ УЧЕБНО-ИССЛЕДОВАТЕЛЬСКИЙ ПРАКТИКУМ
.pdfКонтрольные вопросы
1.Частные уравнения интенсивности источников (стоков) тепла
впотоках: теплообмен.
2.ММ проточного аппарата идеального смешения с учетом источников тепла.
3.ММ проточного аппарата идеального вытеснения с учетом источников тепла.
4.ММ теплообменника типа смешение – смешение.
5.ММ теплообменника типа вытеснение – смешение.
6.ММ теплообменника типа вытеснение – вытеснение.
21
3. МОДЕЛИРОВАНИЕ КАТАЛИТИЧЕСКОГО РЕАКТОРА КОЖУХОТРУБЧАТОГО ТИПА
Цель работы: изучение особенностей моделирования каталитических химических процессов, сопровождающихся выделением или поглощением тепла, на примере реактора кожухотрубчатого типа с катализатором, помещенным в трубное пространство.
Теоретические сведения
Многие процессы химической технологии осуществляются в аппаратах кожухотрубчатого типа. В таких аппаратах чаще всего проводят химические реакции, характеризирующиеся необходимостью подвода тепла извне, либо охлаждения (экзо- и эндотермические реакции). При этом процессы могут идти в присутствии катализатора либо в газофазном потоке. Общая схема такого аппарата представлена на рис. 3.1.
T 0, CA0
Tx0
TC |
PC |
Рис. 3.1. Каталитический реактор кожухотрубчатого типа
В качестве моделируемого процесса в работе рассмотрен процесс риформинга, который в схематичном виде можно представить следующим образом:
22
А → продукты реакции + Q.
Реакция сопровождается выделением тепла, следовательно, для поддержания требуемого температурного режима необходим постоянный отвод тепла. Однако переохлаждение реакционного пространства недопустимо, так как в этом случае реакция затухает.
Для составления математического описания реактора приняты следующие допущения:
1)трубное пространство аппарата следует модели реактора идеального вытеснения;
2)в кожухотрубном пространстве температура постоянна. Скорость химической реакции представляет собой функцию
состава реакционной массы, температуры, давления и других факторов.
Пусть выражение для скорости реакции имеет вид
− |
q |
|
E |
|
|
|
w = KCA = K0e T CA , q = |
, |
(3.1) |
||||
|
||||||
|
|
|
R |
|
где CA – концентрация исходного компонента, кмольм3 ; K – константа скорости реакции, 1мин;
T – температура смеси, К;
R – универсальная газовая постоянная, Джкмоль К; E – энергия активации, Джкмоль; ;
K0i – предэкспоненциальный множитель, 1мин.
Дифференциальное уравнение, описывающее изменение концентрации компонента А, имеет следующий вид:
F |
∂ CA |
= −υ ∂ |
|
CA |
+ FqA , |
(3.2) |
|
|
|
||||
|
∂ t |
∂ l |
|
где F – площадь теплообмена, м2 ;
t– астрономическое время;
υ– объемный расход реакционной массы, м3мин;
qA – интенсивность расходования исходного вещества, кмольм3 мин.
23
Разделив уравнение (3.2) на F и с учетом того, что |
qA = −w, |
||||||||
получим |
|
|
|
|
|
|
|
|
|
|
∂ CA |
= −υ ∂ |
|
CA |
|
− |
q |
|
|
|
|
− |
K0e T CA . |
(3.3) |
|||||
|
∂ t |
|
|||||||
|
F |
|
∂ l |
|
|
|
|
Как упоминалось выше, рассматриваемая реакция является экзотермической (идет с выделением тепла), поэтому в математическом описании необходимо учесть тепловой эффект реакции и теплообмен реагирующего в трубном пространстве потока со средой, находящейся в межтрубном пространстве.
Уравнение теплового баланса имеет вид
|
|
|
|
Fcp |
∂ T |
|
= −υ |
|
cp ∂ |
|
|
T |
+ |
F (−∆ H ) w+ |
|
|
KT F |
(T−х T ), |
(3.4) |
|||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
∂ t |
|
∂ |
|
l |
|
|
|
|
|
|
L |
|
||||||||||
где |
cp |
– удельная теплоемкость исходного реагента, Дж м3 K ; |
||||||||||||||||||||||||||
|
∆ H – изменение энтальпии реакции, |
Дж кмоль; |
|
|||||||||||||||||||||||||
|
KT |
|
– коэффициент теплопередачи, Вт м2 К; |
|
||||||||||||||||||||||||
|
Tх |
– температура хладагента, подаваемого в межтрубное про- |
||||||||||||||||||||||||||
странство, К; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
F |
– суммарная площадь поперечного сечения трубного про- |
||||||||||||||||||||||||||
странства, м2 ; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
L – длина труб, м. |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
Рассмотрим установившийся режим работы реактора, для кото- |
|||||||||||||||||||||||||||
рого |
|
∂ T |
= 0 , |
∂ CA |
= 0 . Тогда уравнения (3.3) и (3.4) приобретают |
|||||||||||||||||||||||
|
|
|
|
|||||||||||||||||||||||||
вид |
|
∂ t |
|
|
∂ t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dCA |
|
|
|
|
F |
− |
q |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
= − |
|
|
|
K0e T CA , |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
dl |
|
|
|
|
v |
|
|
|
|
|
|
|
|
(3.5) |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
dT |
= |
|
F |
(−∆ H ) w+ |
KT F |
(Tх− T ). |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
υ cp |
|
|
|
Lυ cp |
|
|||||||||||||||
|
|
|
|
|
|
dl |
|
|
|
|
|
24
Введя обозначения |
|
F |
K0 = K0′ , |
1 |
(−∆ H )= a1 , |
KT F |
= a2 , сис- |
||||||||
|
υ |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
cp |
Lυ cp |
|||
тему (3.5) можно переписать в следующем виде: |
|
|
|||||||||||||
dCA |
|
|
− |
q |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
= −K0′e T CA , |
|
|
|
|
||||||
dl |
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
(3.6) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
dT |
|
|
|
F |
|
− |
q |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
= |
|
|
a1K0e T CA + a2 (Tх − T ). |
|
|
|||||||
|
|
|
υ |
|
|
||||||||||
dl |
|
|
|
|
|
|
|
|
|
|
|
|
Применительно к кожухотрубчатому каталитическому реактору изменение температуры и концентрации исходного реагента нагляднее изучать не по длине аппарата, а по времени пребывания элементов потока в аппарате (реакционной зоне):
|
= |
LF |
, |
(3.7) |
|
τ |
|||||
υ |
|||||
|
|
|
|
где τ – среднее время пребывания вещества в аппарате.
Выразим текущее время пребывания как функцию конструктивных характеристик аппарата и расхода потока:
τ = |
lF |
, |
(3.8) |
|
υ |
||||
|
|
|
где τ – время пребывания вещества в аппарате.
Из (3.8) можно получить, что d τ = |
|
|
F |
dl |
и |
||||||||||||||||
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
υ |
|
|
|
|
|
|
|
|
dl = |
|
υ |
d τ |
. |
|
|
(3.9) |
|||||||||
|
|
|
|
|
|
F |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
После подстановки (3.9) в (3.6) получим |
|
||||||||||||||||||||
dCA |
|
K0′ |
υ |
− |
q |
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
= − |
|
|
|
e |
|
|
T CA , |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||
|
τ |
|
|
F |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(3.10) |
|
d |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
dT |
|
|
− |
q |
|
|
|
|
|
|
υ |
(Tх − T ). |
|||||||||
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
= a1K0e T CA |
+ |
|
|
|
a2 |
|||||||||||||
|
|
|
|
|
|||||||||||||||||
d τ |
|
|
|
|
|
|
|
|
|
|
|
|
|
F |
|
|
25
Для упрощения системы (3.10) все постоянные коэффициенты можно сгруппировать и переобозначить:
K0′υ |
|
= K0 , |
||
F |
||||
(3.11) |
||||
|
υ |
|||
a2 |
|
= a2 . |
||
F |
|
|||
|
|
|
С учетом (3.11) система (3.10) приобретает окончательный вид
|
dCA |
|
− |
|
q |
|
|
||
|
= −K0e T CA , |
|
|||||||
|
|
||||||||
|
d |
τ |
|
|
|
|
|
|
(3.12) |
|
|
|
|
|
|
|
|||
dT |
|
− |
q |
|
(Tх − T ). |
||||
|
|
||||||||
|
|
|
= a1K0e T CA + a2 |
||||||
|
|
||||||||
d τ |
|
|
|
|
|
|
|
|
Система (3.12) решается любым из известных методов численного или аналитического решения систем дифференциальных урав-
нений при наличии начальных условий (CA (0) = CA0 ; T (0) = T 0 )
с использованием данных о конструктивных характеристиках аппарата, физических свойствах смеси, тепловом эффекте реакции.
Задание и порядок выполнения работы
1.Ознакомиться с основными теоретическими сведениями.
2.В соответствии с вариантом задания (табл. 3.1) составить систему дифференциальных уравнений вида (3.12).
3.Используя один из методов решения систем дифференциальных уравнений программы Pr&Co&Rk.exe, пакетов Matlab или
MathCad, получить кривые состояния реактора по температуре
иконцентрации исходного вещества.
4.Для заданной температуры и концентрации реагентов на входе подобрать температуру хладагента, при которой температура по длине реактора не превышает заданного значения Тогр.
5.Сделать выводы по работе.
26
Содержание отчета
1.Теоретические сведения.
2.Система математического описания процесса.
3.Исходные данные и результат работы (найденное значе-
ние Тх).
|
4. |
Таблицы зависимостей CA = f (τ ), |
T = |
f (τ ). |
|
|
|
|
||||
|
5. |
Графики зависимостей CA = f (τ ), |
T = |
f (τ ). |
|
|
|
|
||||
|
6. |
Выводы по работе. |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
Таблица 3.1 |
|
|
|
|
Исходные данные для моделирования |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
№ |
q, K |
K0, 1/мин |
a1, |
a2, |
|
|
CA0, |
Т |
0 |
, K |
Тогр, K |
|
K м3/кмоль |
1/мин |
|
кмоль/м3 |
|
||||||||
1 |
9000 |
1,00E+11 |
400 |
1 |
|
|
0,1 |
400 |
530 |
|||
2 |
10 000 |
5,00E+10 |
350 |
1,5 |
|
|
0,25 |
425 |
540 |
|||
3 |
9500 |
7,50E+10 |
375 |
1,25 |
|
|
0,2 |
410 |
520 |
|||
4 |
9750 |
1,25E+11 |
425 |
1,1 |
|
|
0,3 |
400 |
545 |
|||
5 |
9000 |
5,00E+10 |
350 |
0,75 |
|
|
0,25 |
425 |
520 |
|||
6 |
10 000 |
7,50E+10 |
375 |
1,25 |
|
|
0,2 |
410 |
530 |
|||
7 |
9500 |
1,25E+11 |
425 |
1,1 |
|
|
0,3 |
400 |
540 |
|||
8 |
9750 |
1,25E+11 |
350 |
0,75 |
|
|
0,25 |
425 |
520 |
|||
9 |
10 000 |
5,00E+10 |
375 |
1,25 |
|
|
0,2 |
410 |
545 |
|||
10 |
9500 |
7,50E+10 |
425 |
1,1 |
|
|
0,3 |
425 |
520 |
|||
11 |
9750 |
5,00E+10 |
350 |
0,75 |
|
|
0,2 |
410 |
530 |
|||
12 |
9000 |
7,50E+10 |
375 |
1 |
|
|
0,3 |
425 |
540 |
|||
13 |
10 000 |
1,25E+11 |
350 |
1,5 |
|
|
0,25 |
410 |
530 |
|||
14 |
9500 |
1,25E+11 |
375 |
1,25 |
|
|
0,2 |
425 |
540 |
|||
15 |
9750 |
5,00E+10 |
425 |
1,1 |
|
|
0,3 |
410 |
520 |
|||
16 |
10 000 |
7,50E+10 |
400 |
0,75 |
|
|
0,25 |
425 |
545 |
|||
17 |
9500 |
1,25E+11 |
350 |
1,25 |
|
|
0,2 |
410 |
520 |
|||
18 |
9750 |
5,00E+10 |
375 |
1,1 |
|
|
0,3 |
400 |
545 |
|||
19 |
9000 |
5,00E+10 |
425 |
0,75 |
|
|
0,25 |
425 |
520 |
27
Окончание табл. 3.1
|
|
|
a1, |
a2, |
CA0, |
0 |
|
Тогр, K |
№ |
q, K |
K0, 1/мин |
K м3/кмоль |
1/мин |
кмоль/м3 |
Т |
, K |
|
20 |
10 000 |
7,50E+10 |
350 |
1 |
0,2 |
410 |
545 |
|
21 |
9500 |
1,25E+11 |
375 |
1,5 |
0,3 |
400 |
520 |
|
22 |
9750 |
1,25E+11 |
425 |
1,25 |
0,25 |
425 |
530 |
|
23 |
9750 |
5,00E+10 |
350 |
1,1 |
0,2 |
410 |
540 |
|
24 |
9000 |
7,50E+10 |
375 |
0,75 |
0,3 |
400 |
520 |
|
25 |
10 000 |
1,25E+11 |
400 |
1,25 |
0,25 |
425 |
545 |
|
26 |
9500 |
5,00E+10 |
350 |
1,1 |
0,2 |
410 |
520 |
|
27 |
9750 |
7,50E+10 |
375 |
0,75 |
0,2 |
400 |
530 |
|
28 |
9750 |
1,25E+11 |
425 |
1,25 |
0,3 |
425 |
520 |
|
29 |
9750 |
1,25E+11 |
350 |
1,1 |
0,25 |
410 |
545 |
|
30 |
9000 |
5,00E+10 |
375 |
0,75 |
0,2 |
410 |
520 |
|
31 |
10 000 |
7,50E+10 |
425 |
1,25 |
0,3 |
400 |
530 |
|
32 |
9500 |
1,25E+11 |
350 |
1,1 |
0,25 |
425 |
540 |
|
33 |
9750 |
7,50E+10 |
375 |
0,75 |
0,2 |
410 |
520 |
Контрольные вопросы
1.ММ проточного аппарата идеального вытеснения с учетом источников вещества и тепла.
2.Частные уравнения интенсивности источников (стоков) вещества и тепла в потоках: химические реакции и теплообмен.
28
4. МОДЕЛИРОВАНИЕ СИСТЕМЫ УПРАВЛЕНИЯ ХИМИЧЕСКИМ РЕАКТОРОМ
Цели работы: получение практических навыков математического моделирования процессов, происходящих в химическом реакторе; получение с помощью вычислительного эксперимента на модели настроек регулятора; моделирование системы управления реактором. Задание на выполнение работы может служить темой курсовой работы студентов по дисциплине «Моделирование систем».
Исходной базой для проектирования любой системы управления являются математические модели типовых технологических процессов. Построение математической модели всегда начинается с составления формализованного описания процессов, происходящих в объекте моделирования. Основными процессами, протекающими в химическом реакторе, являются химические превращения, сопровождающиеся выделением или поглощением тепла (экзотермические и эндотермические реакции). Как правило, наиболее часто встречаются реакции следующих типов:
•линейные, когда происходит последовательное превращение одного вещества в другое: А→ В→ С→ D (например, реакции полимеризации);
•обратимые, когда происходит одновременное протекание как прямого, так и обратного превращения: A B (большинство химических реакций, например окисление сернистого ангидрида в производстве серной кислоты);
•реакции синтеза, когда из двух или более реагентов образуется один целевой продукт реакции: А→ В← С (например, получение нитроглицерина);
•реакции разложения, когда из одного вещества образуются два или более продукта: А← В→ С (например, процесс кальцинации).
На практике в чистом виде такие реакции, как правило, не встречаются. В промышленных химических процессах обычно протекают реакции смешанных типов, например:
29
В C
А |
|
D |
A B |
|
С |
D |
|||
|
|
|||
|
|
|
Скорость химической реакции представляет собой функцию состава реакционной массы, температуры, давления и других факторов. В частности, для реакций на приведенных схемах выражение для скоростей можно представить в виде
|
|
|
|
wi |
= Ki f (C j ), |
|
|
|
(4.1) |
||
где |
i |
– номер реакции; |
|
|
|
|
|
|
|
||
|
j |
– номер реагирующего компонента; |
|
|
|
|
|||||
|
Ki – константа скорости i-й реакции; |
|
|
|
|
||||||
|
f (C j ) – функция от концентрации j -го компонента. |
|
|||||||||
|
В дифференциальной форме уравнение скоростей превращения |
||||||||||
компонентов реагирующей смеси имеет вид |
|
|
|
|
|||||||
|
|
|
dC j |
= ∑ ν ij wi , |
i = 1,l; |
j = 1, n , |
(4.2) |
||||
|
|
|
|
l |
|
|
|
|
|
|
|
|
|
|
dt |
i =1 |
|
|
|
|
|
|
|
где |
ν ij |
– стехиометрический коэффициент |
j -го компонента в i-й |
||||||||
реакции; |
|
|
|
|
|
|
|
|
|||
|
l |
– количество реакций; |
|
|
|
|
|
|
|
||
|
n – количество реагирующих компонентов. |
|
|||||||||
|
Например, для сложной химической реакции типа |
|
1 2
A P S
3
C
выражения для скоростей простых реакций, предполагая для них первый порядок, можно представить в следующем виде:
w1 = K1CA , |
|
|
w2 |
= K2CP , |
(4.3) |
|
= K3CA . |
|
w3 |
|
30