УДК 66.10 - 503.4.001.57
С. И. Дуев
НАХОЖДЕНИЕ МИНИМАЛЬНОЙ ВЕЛИЧИНЫ РЕЦИКЛА ДЛЯ РЕЖИМА С ПОЛНЫМ ИСПОЛЬЗОВАНИЕМ ИСХОДНЫХ РЕАГЕНТОВ
Ключевые слова: реактор с рециклом, множественность стационарных состояний, оптимизация реактора с рециклом.
Рассматривается рециркуляционная система, состоящая из реактора и блока разделения. С использованием метода неопределенных множителей Лагранжа, определена минимальная величина рецикла, при которой возможен режим с полным использованием исходных и промежуточных реагентов. Рассмотрены примеры поиска минимальной величины рецикла при проведении различных реакций.
Key words: reactor with recycle, multiplicity of steady states, optimal value of the recycle.
A recycle system "reactor - separation unit" is considered. Minimum value of the recycle on the operational mode with complete using of basic and intermediate reactants is found by method of undefined multipliers of Lagrange. The examples of the finding of minimum value of the recycle are given.
Введение
При создании химических производств большой мощности особое значение приобретают вопросы интенсификации процессов, возможно более полного использования исходного сырья, минимального загрязнения окружающей среды непрореагировавшими исходными и
промежуточными продуктами химических процессов [1].
Эффективными способами решения проблемы минимизации отходов химического производства является рециркуляция
непрореагировавших исходных веществ [2,3].
При функционировании рецир-куляционной системы реактор - блок разделения большой интерес вызывает режим, при котором достигается полное использование исходных реагентов.
Блок-схема рециркуляционной системы «реактор-блок разделения» представлена на рис.1.
Х®)_т® —►
1 г R
1 x,T 2
F
X™,T"">
Рис. 1 - Блок-схема рециркуляционной системы: реактор - блок разделения: 1 - реактор, 2 - блок разделения
При реализации режима с полным использованием исходных и промежуточных реагентов возникает задача нахождения минимальной величины рецикла, при которой этот режим может существовать.
Минимальная величина рецикла определяется при условии существования режима когда конечные продукты реакции отсутствуют в рецикле [3]:
Р ¿х| = Я (1)
¡=1
где 1 - число исходных и промежуточных реагентов.
Так как Р=Я+0, то из равенства (1) можно получить выражение для определения минимальной величины рецикла:
R = -
¡=i
1-¿*.
(2)
Однако, если на режиме существует континуум стационарных состояний и
концентрации Х^ = 1,...1 могут принимать различные значения, необходимо решить задачу минимизации [4]:
Rmi„ = min-
GX x.
¡=i
'1-EX
(3)
Варьируемыми переменными являются концентрации исходных и промежуточных
реагентов Х^Х2,...Х|. Задача минимизации (3) должна решаться с учетом уравнений для определения концентрации Х1 ,...Х|, которые можно интерпретировать как ограничения типа равенств. Эта задача может быть решена методом неопределенных множителей Лагранжа/5/.
Решим эту задачу, когда в системе функционирует реактор идеального смешения. В стационарном состоянии на режиме с полным использованием исходных и промежуточных реагентов должны выполняться следующие условия:
1.Полный возврат в реактор непрореагировавших исходных и промежуточных реагентов:
РХ^ ЯХ^ = 1,...| (4)
где 1 - число исходных и промежуточных реагентов.
2.Исходные реагенты должен подаваться в систему в стехиометрическом соотношении:
X(0) =--A*J M
(5)
где и - вектор столбец размерности равной числу
элементарных стадий реакций, А * - подматрица матрицы стехиометрических коэффициентов А,
X.....X
=1
X.T
G
состоящие из 1 строк, соответствующих исходным и промежуточным реагентам реакции, а скаляр М вычисляется следующим образом:
М = -1А * й (6)
где 1 - вектор - строка размерности 1, равной числу элементарных стадий реакции.
Тогда, в стационарном состоянии систему уравнений материального баланса реактора идеального смешения на режиме с полным использованием исходных и промежуточных реагентов для 1 исходных и промежуточных реагентов можно записать так [3]:
- О — А*(Уг-—й) = 0 У М '
(7)
где Г - вектор скоростей элементарных стадий реакций, V - объем реактора.
Так как минимизация выражения (3) эквивалентна минимизации суммы концентраций
Х Х| , то функцию Лагранжа [5]
можно записать в
1=1
следующем виде:
I в р
О,
ф(х,...ха,.л)=хх + емл^г - м,) (8)
¡=1 ¡=1 ¡=1 м
где, 3 - множители Лагранжа.
Значения концентраций Х1 ,...Х|, при которых величина рецикла принимает минимальное значение, определяется из следующей системы уравнений:
Э^(х1,...х 1, л1,...л з) _
ЭХ|
Э^(х1,...х 1, Л1,...Л 8)
ЭЛ;
= 0,1 = 1,...1
= 0,, = 1 ,...Б
(9)
Дифференцирование дает следующую систему уравнений для нахождения
Х1, —Х|, Л1,...Л3:
^ ^ ЭХ:
Х Л,Х А^ = 1,к = 1,-1, (10)
¡=1 ¡=1 дХк р О
ХА,(, Мй,)=0,1=1,Lз. -
Подставив найденные из системы уравнений (10) концентрации исходных и
промежуточных реагентов Х1,...Х| в формулу (2) определим минимальное значение рецикла, при котором возможен рассматриваемый режим.
Примеры определения минимальной величины рецикла в системе
1. Реакция А^В, проводимая в изотермическом реакторе идеального вытеснения.
Математическая модель этого реактора в стационарном состоянии может быть записана в виде
—^ = -кх1 в д£ 1
где в - площадь поперечного сечения реактора; к -константа скорости реакции; Х1 - концентрация продукта А, измеряемая в мольных долях; 1 -текущая длина реактора; 0 < 1 < _, Ь - длина реактора.
в
Обозначив а = —, мы имеем следующее решение уравнения (11):
х1(1) = х1 (0)-е-ак1, (12)
где Х1(0) - концентрация продукта А на входе в реактор.
Так как концентрации измеряются в мольных долях, то концентрация продукта В в реакторе определяется так:
х2 (1 ) = 1-х1(0)-е
ак1
(13)
На режиме с полным использованием исходного продукта А уравнение (11) должно быть дополнено следующим краевым условием:
—х1(0) = Ох/0) + —хД_),
(14)
(0)
где Х1 - концентрация компонента А на входе в рециркуляционную систему.
Решение уравнения (11) с краевым условием (14) даёт следующий результат:
П ааИ
1У' —1-еак1
О-ехр-
в-к -1 Я+О
(Я+О)1-ехр-
в-к -1 Я+О
(15)
Концентрация продукта В определяется из соотношения (13). Таким образом, на рассматриваемом режиме существует единственное стационарное состояние.
Запишем условия существования режима:
хдд<
Я
Я + О
(16)
Минимальная величина Ят:п определяется
как решение следующего уравнения: ~ ( в - к -_ \ О - ехр[ -Я+О ] =
Б - к - Я + в
(Я + О)[ 1 - ехр[ -
Я + в
(17)
Так, при значениях параметров, равных О = 4, в = 1, _ = 10, к = 1, концентрация исходного продукта А на выходе реактора равна Х1 (_) = 0,1 07, а минимальная величина рецикла -Ятш = 0,481.
2. Реакция А^В^С в политропическом реакторе идеального смешения.
Математическое модель реактора идеального смешения в стационарном состоянии на
Г
режиме с полным использованием исходного реагента А и промежуточного реагента В записывается так:
- Vr1 + в = 0
(18)
+ Угх - Уг2 = 0 ОсррТ(0) + я- СррТв - —Срр - Т + . V(г1AH1 + Г2 ЛИ2) + и(Т - Тх) = 0
Для простоты исследования положем, что скорости элементарных стадий реакции
определяются соотношениями: Г, = к^, Г2 — к 2 х 2,
Отсюда определяем выражения для нахождения концентраций А, В:
О
Х1 =
Ук1 О Ук,
(19)
в которых константа скорости к| (i = 1,2) зависит от температуры в соответствии с формулой Арренниуса:
k¡ = А^ ЯгТ,1 = 1,2, (20)
где А; - предъэкспоненциальный множитель, Е| -
энергия активации, ¡ - ой стадии реакции.
Единственное стационарное значение температуры на режиме определяется из следующей формулой:
Т=-
1
(О<ррТ°> +иТ+ясррТЯ +Олн+лн)) (21)
<^р+и
Минимальная величина рецикла, при которой возможен режим с полным использованием реагентов А и В, определяется из равенства(2):
Я =
Ят1п
О2 ( V
V к1
1 1
—+—
к
л
2 у
О V V
л
(22)
11
—+—
Л к2 У
Совместное решение уравнений (22) и (21) с учетом (20) даёт минимальное значение рецикла и температуру в реакторе, при которой это значение достигается.
3. Реакция А+В^2С в политропическом ректоре идеального смешения.
Материальный баланс рактива идеального смешения для исходных реагентов А и В можно записать следующим образом:
- Vr + Ох
(0)
Vг + Ох
(0)
—х1 + Ях* = 0 —х2 + Ях*2 = 0
(23)
На режиме с полным использованием исходных реагентов А и В должны выполнятся условия (4):
—x¡ = Ях*, I = 1,2 (24)
Учитывая, что на этом режиме исходные реагенты должны подаваться в систему в
стехеометрическом соотношении (5) (х(0) = Х 20)) и концентрации реагентов измеряется в мольных долях функцию Лагранжа (4) можно записать так:
(25)
ф(х1,х2, Л) = х1 + х2 + л[ Vг - О2
Пусть скорость реакции определяется простейшей формулой: Г = к - Х1 - Х2
После взятия производных система уравнений для нахождения неизвестных Х^Х2, Л запишется так:
^ = 1 + Л - V - кх2 = 0
Эх
Эф
Эх,
= 1 + Л - V - кх1 = 0
(26)
дф г, . О _ — = V -кх, - х2--= 0
дЛ
2
Отсюда видно, что минимальная величина рецикла достигается при равенстве концентраций:
Х1 — Х2 —
О
2Vk
(27)
Подставляя найденные значения Х1 и Х2 в (2) определим выражение для минимальной величины рецикла:
О - 2 -
Я = -
Ят1п
О
2Vk
(28)
1-2.
О
2Vk
Если реактор политропический, то константа скорости реакции будет зависит от температуры по формуле Аррениуса:
ЕЯг
к(Т) = А - е Т (29)
Уравнение теплового баланса для реактора идеального смешения записывается следующим образом:
ОсррТ0) + ^дНЦх2 +и(ТХ -Т)--—СррТ+ЯсррГ = 0
С учетом последнего уравнения системы (26) температура в реакторе определяется так:
(30)
Т =-
1
—Ср + и
(Сзр(ОТ(0) + ЯТ*) +
(31)
+ иТХ + О - х1(0)(лН))
Для нахождения Ят|п уравнение (31) должно решаться совместно с уравнением (28).
Подставляя выражение для Ят|п из (28) в (31) получаем уравнение для нахождения
2
температуры, при которой достигается минимальное значение рецикла в системе.
Т =^^(сррОТ0 + УТх +
РСр + и
О •
+ Ох1(0 )^(ЛИ}+сррЯТ*--1
2О 2Ук (Т)
(32)
1-2
1
О
Ук(Т)
Затем при этом значение температуры по формуле (28) вычисляется минимальная величина рецикла.
Заключение
Таким образом, с использованием неопределенных множителей Лагранжа найдена минимальная величина рецикла для режима с полным использованием исходных реагентов на котором существует континуум стационарных состояний.
Обозначения
А - предъэкспоненциальный множитель ср- удельная теплоёмкость, Дж/.К Ерэнергия активации 1-ой стадии реакции, Дж/.моль Р- поток смеси, поступающей в реактор, м3/с О-поток смеси, поступающей в систему, м3/с -AH¡-тепловой эффект 1-ой стадии реакции, Дж/.моль кгконстанта скорости 1-ой стадии реакции П-скорость 1-ой элементарной стадии рекции
Я-поток смеси в рецикле, м3/с ^-универсальная газовая постоянная, Дж/.моль, Б - площадь поперечного сечения реактора, м2 Т-температура смеси, ТХ- температура хладоагента, 1- время,
и коэффициент теплопередачи, Вт/м2*К У-объём реактора, м3
Х - вектор концентраций Х¡ -концентрация 1-ого реагента, мол.доли. ^-плотность смеси, кг/м3 Индексы: (0)- вход
¡ - номер реагента
- рецикл
Литература
1. Кафаров В.В. Принципы создания безотходных химических производств В.В. Кафаров// М.: Химия, 1982. 288с.
2. Дуев С.И. Расчет стационарных состояний реактора в рециркуляционой системе реактор-блока разделения/ С .И. Дуев// Вестник КГТУ, 2012 № 16 с .151 -153.
3. Бояринов А.И. Множественность стационарных состояний в системе: система - реактор - узел разделения/ А.И. Бояринов, С.И. Дуев// Теоретические основы химической технологий, 1980. Т.14. №6. С.903.
4. Дуев С.И. Расчет режимов с полным использованием исходных реагентов в рециркуляционной системе реактор-блок разделения. /С.И. Дуев.// Вестник КГТУ, 2013, №16, с. 167-169.
5. Бояринов А.И. Методы оптимизации в химии и химической технологии / А.И. Бояринов, В.В. Кафаров// М. Химия, 1975, 576с.
© С. И. Дуев - д-р техн. наук, проф. каф. информатики и прикладной математики КНИТУ, douev@mai1.ru.
*