С. И. Дуев
РАСЧЕТ СТАЦИОНАРНЫХ СОСТОЯНИЙ РЕАКТОРА В РЕЦИРКУЛЯЦИОННОЙ СИСТЕМЕ «РЕАКТОР-БЛОК РАЗДЕЛЕНИЯ»
Ключевые слова: реактор с рециклом, множественность стационарных состояний, устойчивость реакторных систем.
Рассматривается рециркуляционная система, состоящая из реактора и блока разделения. Проведен анализ стационарных состояний для реакции A + B ^ C + D, протекающей в реакторе идеального вытеснения в системе реактор - блок разделения. Показано, что на режиме с полным использованием исходных расчетов A и B существует континуум стационарных состояний. Определена минимальная величина режима, при которой возможен режим с полным использованием исходных реагентов A и B.
Keywords: reactor with recycle, multiplicity of steady states, stability of reactor system.
The recycle system reactor - separation unit is considered. Analysis of the steady states for the reaction A + B ^ C +
D takes place in the tubular reactor is done. It is proved that a continuum of the steady states exists on the regime with the full using of basic reactions A and B . The minimum of the value of the recycle for the regime with the full using of basic reactants is defined.
Введение
При создании химических производства большой мощности особое значение приобретают вопросы интенсификации процессов, возможно более полного использования исходного сырья, минимального загрязнения окружающей среды непрореагировавшими исходными,
промежуточными и побочными продуктами химических процессов [1]. Экологические проблемы приобретают в последнее время особую значимость [2,3].
Эффективным способом решения проблемы минимизации отходов химического производства является рециркуляция непрореагировавших исходных веществ.
Блок-схема рециркуляционной системы реактор-блок разделения представлена на рис. 1,
R*:i
Рис. 1 - Блок-схема рециркуляционной системы реактор-блок разделения
здесь, G- количество смеси, поступающими в систему в единицу времени, F - количество смеси поступающей в реактор в единицу времени, x -вектор концентрации в реакторе (со значком «0»- на входе в систему, со значком «*»- в рецикле, со значком «вых»-на выходе системы).
Предполагается, что блок разделения в общем случае может включать любые процессы разделения, такие как ректификация, абсорбация, экстрация и т.п. его функционирования определяется только заданием режимных
параметров процессов разделения, таких как флегмовое число, величина орошения в абсорберах и т.п., не регулируемых по составу получаемых продуктов. Поставленные условия соответствует принципу стабилизации потоков в системе
разделения. Между потоками существует однозначные соответствия: F = R + G. В этой системе, выделив целевой продукт в
технологическом операторе разделения, непрорегировавшие химические компоненты сырья через оператор смешения возвращают в оператор химического превращения. Типичными примерами данной рециркуляционной системы являются химико-технологические системы синтеза аммиака, синтезов этилового спирта из оксида углерода и водорода, этилового спирта, каталитической гидратации этилена в паровой фазе, производства ацетальдегида, гидратация ацетилена в жидкой фазе и другие.
Математическое моделирование реактора идеального вытеснения с рециклом
Пусть в изотермическом реакторе идеального вытеснения функционирующем в рециркуляционной системе протекает обратимая реакция второго порядка А + В ^ С + О.
Математическая модель реактора в стационарном состоянии может быть представлена следующей системой дифференциальных уравнений:
Р dx
д-^1 =-Г1(Х1,Х2) + Г2(Х3,Х4) (1)
^^Х2 =-Г1(Х1,Х2) + Г2(Х3,Х4)
Р ^3 , X , X
3 = Г1(Х1,Х2) + Г2(Х3,Х4)
о “I
р аХ
55 “О4 = + Г*2 (Х3, Х4 )
где -площадь поперечного сечения реактора, Х1, Х2, Х3,Х4 - концентрации реагентов А, В, С и й в реакторе. I - текущая длина реактора, Г^ = 1.2) - функции, определяющие скорости
прямой и обратной стадий реакции.
Для простоты анализа будем полагать, что Г,= к1Х1Х2 , Г2 = к2Х3Х4, где Х1(1 = 1.2) -
константы скоростей для прямой и обратной стадии реакции. Заметим, что для случая более сложной зависимости функции Г1= (Х1,Х2) и Г2= (Хз,Х4)
качественные результаты будут такими же. Граничные условия для реактора, функционирующего в рециркуляционной системе запишутся следующим образом:
РХ1(0) = 0х(0) + Рх* (2)
Рх2(0) = 0х20) + Рх2 РХз(О) = 0х3О) + Рх3 РХ4(0) = 0Х(4О) + РХ 4
где Х(0),Х|(1 = 1.4) - концентрация реагентов на
входе в систему и в рецикле.
Предположим, что системы разделения обладает достаточно высокой разделительной способностью для полного отделения непрореагировавших реагентов А и В от конечных продуктов реакции.
Рассмотрим режим, при котором в рецикле присутствуют три компонента реакции А, В и С. Причем, непрореагировавшие в реакторе исходные расчеты А и В полностью рециркулируют в реактор после стадии отделения их от конечных продуктов реакции, а один из конечных продуктов реакции С является распределяющимся, т. е. присутствует и в рецикле и на выходе системы.
Полагая, что концентрации компонентов во всех потоках измеряются в мольных долях для режима с полным использованием исходных реагентов А и В можно записать:
Рх* = РХ1(Ь) (3)
Рх2 = Рх2(Ь)
РХ3 = Р(1 — Х1 — Х2)
рх4 = 0
Тогда на режиме с полным использованием исходных реагентов А и В математическая модель реактора (1) и граничные условия (2) запишутся так:
р Сх1
Б СІІ
Р Сх2
Б СІ
Р Схз
Б СІ
= -к1х1х2 + к2х3х4 = -к1х1х2 + к2х3х4 = к1х1х2 - к2х3х4 0
х1(0) = 0х(0) + х^)
(4)
р 0 р
(5)
х2(0) = ^Т х20) + х2(Ь)
0
Хз(0) = рх30) + (Р — Х1(1_) — Х2(Ь))/Р
Х4 = 1 — Х1 — Х2 — Х3
Для функционирования этого режима необходимо, чтобы исходные реагенты А и В
подавались в систему в стехиометрическом соотношении, т.е. х(0) = х20). [4]
Поскольку первые два уравнения системы(4) и соответствующие или граничные условия (5) совпадают, то краевая задача(4), (5) будет иметь бесконечное множество (континуум) решений.
Поэтому, на рассматриваемом режиме будет существовать семейство стационарных состояний, в котором концентрации реагентов А и В будут принимать континуум стационарных значений. Результаты расчета краевой задачи (4), (5)
представлены на рис.2., где представлен вид семейства стационарных состояний на выходе
реактора в проекции на плоскость х1х2
Д.
Р
Рис. 2 - Вид семейства стационарных состояний (кривая аЬ) на выходе реактора на плоскости
х1х2
Стационарные концентрации исходных реагентов А и В на выходе реактора могут
приносить значения в
интервале хітіп < хі <хітах, і = 1.2 . Гранщы
интервала возможных стационарных значений
концентрации хі-хітіп, хітах ,і = 1.2 определяется
из условия существования режима с полным использованием исходных реагентов А и В [5]
х., + х2 <
р
г_ ^ (6)
Р + 0
Максимальные значения концентрации целевого продукта Э - Х4 достигается при равенстве концентрации А и В : Х1 = Х2. Как доказано в работе [7] это семейство стационарных состояний может в лучшем случае находиться на границе области устойчивости.
Минимальная величина рецикла, при
которой возможен режим с полным использованием исходных реагентов определяется при условии отсутствия конечных продуктов А и В в рецикле [6] Р(Х1 + Х2) = Р (7)
Так как Р = Р + 0, то из равенства (7) можно получить выражение для определения минимальной величины рецикла:
0(Х1 + Х2)
Р
1 — х1 — х2
(8)
В случае, если на режиме существует семейство стационарных состояний и концентрации Х^ = 1.2 могут принимать
различные значения в интервале
Х:т:п < Х: < Х: ,1 = 1.2 необходимо решить
задачу минимизации
= тіп С(х1 + х2)
Р
Хі Х2 1 — х1 — х2
(9)
Здесь варьируемыми переменными являются концентрации исходных реагентов А и В. Задача минимизации выражения (9) должны решаться с учетом системы уравнений для определения концентраций (4),(5). Для данной реакции она определится при равенстве концентраций Х1 = Х2 и равна:
= 20х1
Р
1- 2х1
(10)
Заключение
Таким образом, доказано существование континуума стационарных состояний на режиме с полным использованием исходных реагентов А и В, для реакции А+ В ^ О + й , проводимой в реакторе идеального вытеснения. Этот континуум стационарных состояний может быть устойчивым только на границе области устойчивости, поэтому для реализации рассматриваемого режима
необходима
регулирования.
система
автоматического
Литература
1. Кафаров, В.В. Принципы создания безотходных химических производств / В.В.Кафаров //М.: Химия -1982-288с.
2. Гафаров, А.Х. Мониторинг вредных выбросов при сжигании природного газа предприятий по выработке тепловой энергии в районах РТ // А.Х.Гафаров, Л.И.Лаптева // Вестник Казан. технол. ун-та.-2010 - №1-С/463-467.
3. Суркова, А.В. Аспекты экологических проблем
хромового дубления / А.В.Суркова,
И.Ш.Абдуллин,Б.Л.Журавлев, Ю.С.Парсанов // Вестник Казан. технол. ун-та.-20о9 - №6-С.91-94.
4. Бояринов, А.И. Множественность стационарных состояний в системе: смеситель - реактор - узел разделения /А.И.Бояринов, С.И.Дуев // Теоретические основы химической технологии - 1980 - №6 - Т.14 -С.903
5. Дуев, С. И Исследование режимов с полным использованием исходных и промежуточных реагентов в системе реактор-узел разделения / С.И.Дуев, А.И.Бояринцев, В.В. Кафаров // Системный анализ процессов химической технологии - М.,МХТИ-1979.-Вып.106.
6. Бояринов, А.И. Анализ стационарных состояний реактора идеального вытеснения с рециклом/ А.И.Бояринов, С.И.Дуев // Теор.основы хим.технол. -1988. - Т.22.-С402-404.
© С. И. Дуев - д-р тхн. наук, проф., каф. информатики и прикладной математики КНИТУ, douev@mail.ru.