УДК 66.10 - 503.4.001.57 С. И. Дуев
РАСЧЕТ РЕЖИМОВ С ПОЛНЫМ ИСПОЛЬЗОВАНИЕМ ИСХОДНЫХ РЕАГЕНТОВ В РЕЦИРКУЛЯЦИОННОЙ СИСТЕМЕ «РЕАКТОР-БЛОК РАЗДЕЛЕНИЯ»
Ключевые слова: реактор с рециклом, множественность стационарных состояний, устойчивость реактора с рециклом.
Рассматривается рециркуляционная система, состоящая из реактора и блока разделения. Проведен анализ стационарных состояний для реакции A + B -о C ^ P, проводимой в реакторе идеального смешения. Показано существование континуумов стационарных состояний и проведен расчет континуумов стационарных состояний на режимах с полным использованием исходных реагентов. Определены области стационарных значений концентраций на этих режимах и минимальные величины рециклов, необходимых для их функционирования.
Key words: reactor with recycle, multiplicity of steady states, stability of reactor with recycle.
The recycle system reactor - separation unit is considered. Analysis of the type of steady states for the reaction A + B о C ^ P takes place in the continuous stirred tank reactor is done. Calculation of continiums of the steady states for the regimes with the full using of basic reactants is given. The minimum of the values of the recycle for the regimes with the full using of basic reactants are defined.
Введение
При создании химических производств большой мощности особое значение приобретают вопросы интенсификации процессов, возможно более полного использования исходного сырья, минимального загрязнения окружающей среды непрореагировавшими исходными,
промежуточными и побочными продуктами химических процессов [1]. Экологические проблемы приобретают в последнее время особую значимость [2,3].
Эффективным способом решения проблемы минимизации отходов химического производства является рециркуляция непрореагировавших исходных веществ.[4]
Блок-схема рециркуляционной системы реактор - блок разделения представлена на рис. 1,
Rx*
F
G реактор блок разделения G
Рис. 1 - Блок схема рециркуляционной системы реактор - блок разделения
Здесь С - количество смеси, поступающее в систему в единицу времени, Р - количество рециркуляционной смеси в единицу времени, Р -количество смеси, поступающее в реактор в единицу времени, х - вектор концентраций в реакторе (со значком «0» - на входе в систему, со значком «*» - в рецикле, со значком «вых» - на на выходе системы).
2. Математическое моделирование режимов функционирования реактора идеального смешения с рециклом.
Пусть в политропическом реакторе идеального смешения, функционирующем в
рециркуляционной системе реактор-блок разделения
протекает последовательная реакция
А + В ^ С ^ Р.
Математическая модель политропического реактора идеального смешения в стационарном состоянии в рециркуляционной системе может быть представлена следующей системой уравнений.
Gx(0) - Vr, + Vr2 - Fx1 + Rx* = 0 Gx20) - Vr + Vr2 - Fx2 + Rx2 = 0 Vr — Vr2 — Vr3 — Fx3 + Rx3 = 0 Vr3 - Fx4 + Rx4 = 0 cpp(GT(0) + RT* - FT) + U(Tx - T) + V((Q1 - Q2)r1 + Q3r3) = 0
(1)
где х1>х2’Хз’Х4 - концентрации реагентов А, В, С и Р соответственно, Г (I = 1,2,3) - скорости
элементарных стадий реакции, 0| (I = 1,2,3) -
тепловые эффекты элементарных стадий реакции.
Полагаем, что концентрации измеряются в мольных долях, тогда концентрацию конечного продукта Р можно определить следующим образом:
x4 = 1 - х 1 - х2 - х3.
(2)
На режиме с полным использованием исходных и промежуточных реагентов А, В и С должны выполняться следующие условия [5]:
Rx| = Fxi, где i = 1,2,3.
(3)
Для организации этого режима необходимо также, чтобы исходные реагенты А и В подавались в систему в стехнометрическом соотношении [5]:
х(0) = х20) = 0,5.
Тогда система уравнений (1) рассматриваемом режиме запишется так:
на
— - Уіг + Уг2 = 0
2 1 2
- Уг., + Уг2 = 0
Уг, - Уг2 - Уг3 = 0 (4)
срр(СТ(0) + РТ* - РТ) + и(Тх - Т) +
У((01 - 02)г1 + Озг1) = 0
Пусть скорости элементарных стадий реакции имеют следующий вид:
г1 = [<1X1X2, г2 = к2х2, г3 = к3х3, где к1,к2,кз-
константы скоростей элементарных стадий реакции, определяющиеся по формуле Аррениуса:
ітах
К С
V Р 2Укэ ,
Г к с л 2 С Г к21
1Р 2Укэ , 2Ук2 V кэ у
і = 1,2
Вид континуума стационарных состояний (кривая аЬ) представлен на рис.2а.
При этом стационарное значение температуры на режиме определяется однозначно из уравнения теплового баланса:
1
Т
С
2
С рТ(0) + К рТ* + ит +
Рсрр + и V ср ср х
((<°1 - 02) + Оэ)
(11)
-Е.
____|_
к т
кі = Ае г , і = 1,2,3
(5)
Тогда, концентрация промежуточного
продукта С определяется так:
X, =■
С
2Ук
(6)
Минимальная величина рецикла, при
которой возможен рассматриваемый режим
определяется при равенстве стационарных значений концентраций А и В (х1 = Х2),
Х1 = х, =
Ук1
( С (
V V
„ к2
1 + -^-
кэ
/у
из условия (8):
(12)
Для определения концентраций исходных реагентов А и В имеется только одно уравнение, которое можно записать как:
Х2 =
С
2Укх
1 1
( \ 1 + ^2 кэ
(7)
Здесь концентрация реагента А - Х1 имеет
роль параметра, который может принимать любые значения в пределах интервала стационарных значений, границы которого определяются из условия существования режима с полным использованием исходных и промежуточных реагентов [5]:
к
X + X + X < — 1 2 э Р
(8)
2С
к
С(к2 + кэ) + ~С~
2Ук1кэ 2Укэ
1-
С(к2 + кэ) + С ~ 2Ук1кэ 2Укэ
(13)
Рассматривается реакция интересна еще и тем, что континуум стационарных состояний существует также и на режиме с полным использованием только исходных реагентов А и В, когда промежуточный продукт реакции С присутствует как в рецикле, так и на выходе системы. В этом случае концентрацию промежуточного продукта С в рецикле можно определить так:
Хэ = 1-р- (Х1 + Х2 )
(14)
При равенстве соотношения (8) с учетом формул (6) и (7) определяем минимальные и максимальные значения стационарных значений концентраций х1 , а следовательно, и х :
и следовательно, уравнение материального баланса реактора на этом режиме запишутся так:
С
- Ук1х1х2 + Ук2хэ +— = 0
X;
ітіп
К С
Л
Р 2Ук
К С
2Укэ
С
(
2Ук1
Л
(9)
С
- Ук^Х2 + Ук2Хэ + — = 0
Ук1Х1Х2 — Ук2Хэ — УкэХэ — Рхэ + К — Р(Х1 + Х2) = 0
(15)
Отсюда можно выразить концентрации х2
и хз через концентрацию-параметр Х1 :
1
э
2
2
Р
Vk2[ G + R ] +G (Vk3 + F)-Vk2Fx1
x2 =------------------^-=-)—2------- (16)
2 Vk^CV^ + F) + Vk2F
x3 =
G
— + R - F(x2 + х1) Vk3F
(17)
Так как концентрация Х1 может принимать любые значения, ограниченные условием
существования этого режима:
x1 + х2 <
R
F
(18)
то на этом режиме также существует континуум стационарных состояний, в котором не только концентрации х1 и х2 но и концентрация
промежуточного продукта С-Хэ может принимать любые стационарные значения из интервалов, ограниченных условием существования режима (18). Границы интервала стационарных значений
для концентраций Х1, Х2 и Хэ определяются из
решения системы уравнений (16), (17) и
следующего уравнения:
х.+х2 = ■
R
F
(19)
Вид континуума стационарных состояний на плоскости Х1Х2 представлен на рис.2б.
Минимальная величина рецикла на этом режиме ктіп2 определяется при равенстве
концентраций Х1 = Х2 из равенства (19). При расчетах показано, что К . > К . „.
А тіп1 тіп2
Рис. 2а - Вид континуума стационарных
состояний (кривая аЬ) на плоскости х^, при
R = 0.6, F = 1, i = 1,2
R „ = 0.4,
mini
0.05 < xi < 0.35,
Рис. 2б - Вид континуума стационарных
состояний (кривая a' b' ) на плоскости x^2, при R = 0.2, R . „ = 0.18, 0.06 < xi < 0.14, i = 1,2
min2 i
Заключение
Таким образом, показано, что для реакции A + B о C ^ P, проводимой в реакторе идеального смешения в рециркуляционной системе реактор-блок разделения на режимах с полным использованием исходных реагентов существует континуум (бесконечное множество) стационарных состояний. Причем в зависимости от величины рецикла в системе может существовать одномерный континуум (кривая на плоскости) на режиме с полным использованием исходных и промежуточных реагентов А, В и C и двухмерный континуум (поверхность в пространстве) на режиме с полным использовнием только исходных реагентов А и В.
Литература
1. Кафаров, В.В. Принципы создания безотходных химических производств / В. В. Кафаров // М.: Химия -1982 - 288с.
2. Дуев, С.И. Исследование режима с полным использованием исходных реагентов в рециркуляционной системе реактор-блок разделения / С.И. Дуев // Вестник КГТУ, 2010 - №10.
3. Дуев, С.И. Расчет стационарных состояний реактора в рециркуляционной системе реактор - блок разделения / Дуев, С.И. // Вестник КГТУ, 2012 - №16 - С.151-153.
4. Duev, S.I. Dinamic behaviour of the recycle system reactor-distillation column / S.I. Duev, A.I. Boyarinov // Proc.Int.conf. Distillation and absorption. Baden-Baden. Germany - 2002.
5. Бояринов, А.И. Множественность стационарных
состояний в системе: смеситель - реактор - узел разделения / А.И. Бояринов, С.И. Дуев // Теоретические основы химической технологии - 1980 -№6 - Т.14 -
С.903.
© С. И. Дуев - д.т.н., проф. каф. информатики и прикладной математики КНИТУ, douev@mail.ru.