- •В.А. Тихонов математическое моделирование химико-технологических процессов
- •Ведение
- •1 Математическое моделирование хтп
- •1.1. Математическое моделирование
- •1.2. Химико-технологический процесс как объект математического моделирования
- •1.3. Основные стадии математического моделирования
- •1.4. Основные группы уравнений, входящих в математическое описание процесса
- •1.5. Блочный принцип построения математической модели
- •Математические модели химических реакторов
- •.1. Реактор идеального смешения
- •Vвх ρвх сРвх (Твх – Тну) dt – Vвых ρ ср (т – Тну) dt –
- •.1.1 Стационарный режим работы рис
- •Vвх ≠ Vвых
- •Vвх ≈ Vвых ≈ V
- •Vвх ≠ Vвых
- •.1.2 Динамический режим работы рис
- •.2. Реактор идеального вытеснения
- •.2.1. Решение математического описания
- •Литература
- •Математическое моделирование химико-технологических процессов
.2. Реактор идеального вытеснения
Исследуется химическая реакция в аппарате идеального вытеснения (РИВ) (рис.):
Рис. Реактор идеального вытеснения
Полное математическое описание (детерминированная математическая модель) процесса будет представлено покомпонентным материальным балансом и тепловым балансом элементарной ячейки с объёмом dV для малого промежутка времени dt, так как параметры потока меняются по длине аппарата и во времени.
Материальный баланс элементарной ячейки реактора РИВ:
mi dt – приход вещества с потоком,
mi dt = Vвх Сi dt = u S Ci dt ;
(mi + dmi) dt – расход вещества с потоком,
(mi + dmi) dt = (Vвх Ci + d(Vвх Ci)) dt ;
dV ri dt – изменение вещества в химической реакции,
dMi – накопление вещества,
dMi = dCi dV ;
mi dt – (mi + dmi) dt + dV ri dt = dCi dV ,
где ri – скорость изменения концентрации вещества i в результате химического превращения, моль/(м3 с),
mi = Vвх Сi = u S Ci – мольный поток вещества i, моль/с,
u – линейная скорость потока, м/с.
Преобразуя, поделим все члены уравнения на (dV dt). Конечный вид уравнения материального баланса реактора идеального вытеснения для вещества i в размерности моль/(м3 с) имеет вид:
(2.1)
Решением этого уравнения будет являться функция Ci = f (Vвх, t).
В случае стационарного режима работы реактора:
,
.
(2.2)
Решение его – функция Ci = f (Vвх).
При Vвх = Vвых = V (мономолекулярная реакция) Твх = Твых = Т:
.
(2.3)
Решение этого уравнения – функция Ci = f (z, t).
В случае стационарного режима работы реактора:
,
или
, (2.4)
где
,
при u,V
– const
по длине реактора.
Решение уравнения (2.4) – функция Ci = f (z) или Ci = f ().
Тепловой баланс элементарной ячейки РИВ.
Vвх ρ СР (Т – Тну) dt – теплосодержание входного потока;
[Vвх ρ СР (Т – Тну) + d(Vвх ρ СР Т)] dt – теплосодержание выходного потока (аналогично (mi + dmi) dt в материальном балансе);
– тепло химического
превращения;
KF (T – TS) L dz dt – теплоперенос через стенку;
d(ρ CP T) dV – накопление тепла (аналогично dCi dV в материальном балансе);
где (ρ CP T), Дж/м3 – аналог концентрации Ci, моль/м3,
(Vвх ρ СР Т), Дж/с – аналог мольного потока mi, моль/с.
j = 1, m – количество стадий реакции;
НР – энтальпия химического превращения (тепловой эффект стадии с обратным знаком), Дж/моль;
KF – коэффициент теплопередачи через стенку, Дж/(м2 с K);
L dz – поверхность теплообмена, м2;
TS – температура хладоагента, K.
Тепловой баланс элементарной ячейки:
Vвх ρ СР (Т – Тну) dt – [Vвх ρ СР (Т – Тну) + d(Vвх ρ СР Т)] dt –
– – KF (T – TS) L dz dt = d(ρ CP T) dV
Преобразуя, поделим все члены уравнения на (dV dt). Конечный вид уравнения теплового баланса реактора идеального вытеснения в размерности Дж/(м3 с):
(2.5)
Уравнение для адиабатического реактора упрощается за счет того, что:
Практически это возможно при KF = 0 (футеровка, теплоизоляция) или T = TS (специально организованный температурный режим у внешней поверхности реактора).
Параметры потока Vвх, ρ, СР, T, Ci, согласно уравнению (2.5), меняются по длине реактора идеального вытеснения и во времени.
В практических расчётах можно пользоваться упрощённой формой уравнения теплового баланса адиабатического реактора:
(2.6)
Решением этого уравнения является функция T = f (Vвх, t).
Стационарный режим работы реактора:
(2.7)
Его решением является функция T = f (Vвх).
При использовании в качестве независимой переменной линейной координаты
уравнение (2.6) преобразуется к виду:
(2.8)
Решением этого уравнения является функция T = f (z, t).
Для стационарного режима работы:
(2.9)
Его решением является функция T = f (z).
