Prezentatsia_lektsia_Reaktora_2_Materialny_balans

Математическая модель химического реактора включает два уравнения:
- уравнение материального баланса
- уравнение теплового баланса
В общем виде уравнение материального баланса применительно к реагенту можно представить в виде:
количество вещества, поступающего на вход в реактор в единицу времени – расход вещества в единицу времени (на химическую реакцию, а также выходящего из реактора) = накопление (убыль) вещества в реакторе в единицу времени
По аналогии уравнение теплового баланса в общем виде:
приход тепла в реактор в единицу времени – расход тепла в единицу времени = накопление (убыль) тепла в реакторе в единицу времени
Для реактора, работающего в стационарном режиме, когда концентрация и температура во времени не меняются, правая часть обоих уравнений равна 0.
2

Уравнение материального баланса химического реактора
Варианты:
Общий (брутто-баланс)
Частный баланс по компоненту Элементный баланс
3

Общий материальный баланс реактора подразумевает равенство массовых потоков на входе и выходе из аппарата (составляется по всем компонентам системы)
4

1, 5, 6 – подогреватели, 2- колонна переиспарения 1-ФЭТ, 3 – кипятильник, 4, 20 – роторно-плёночные испарители,
7, 8 – реакторы, 9 – трубчатая печь, 10 – скруббер, 11 – воздушные конденсаторы, 12 – водяной конденсатор,
13 – отстойник, 14, 16, 18 – ректификационные колонны, 15, 17, 19 – конденсаторы.
5

6

Частный материальный баланс реактора по компоненту реакционной системы (реагенту или продукту) составляется с учетом расходования или накопления вещества в результате химической реакции:
wA,i – концентрация вещества А в i-ом потоке (размерность концентрации ???)
ΔGA, хр – масса вещества А расходующаяся (или образующаяся) в результате
химической реакции в единицу времени, кг/ч
7

Элементный баланс может быть составлен по любому элементу, например:
по углероду для органических смесей
по азоту
по сере
по кислороду
Как и для других видов материального баланса, количество элемента, по которому составляется баланс, должно быть одинаковым на входе и выходе из системы.
8

Уравнение материального баланса для элементарного объема аппарата, в котором протекает химическая реакция:
А → В
где СА – концентрация реагента А в реакционной смеси;
x,y,z – пространственные координаты;
Wx, Wy, Wz – составляющие скорости потока по осям х, у, z;
D – коэффициент диффузии вещества А;
rA – скорость химической реакции
Левая часть уравнения
характеризует изменение концентрации вещества
А во времени в элементарном объеме dV
Уравнение конвективного массообмена, дополненное слагаемым rA, которое
учитывает протекание химической реакции
9

слагаемое
отражает изменение концентрации реагента вследствие конвективного переноса
слагаемое
отражает изменение концентрации в результате переноса вещества
путем молекулярной диффузии
слагаемое rA показывает изменение концентрации за счет химической реакции
10

Гидродинамические модели реакторов. Вывод характеристических уравнений.
Любой химический реактор в большем или меньшем приближении можно описать одной из следующих моделей:
реактор идеального смешения периодического действия РИС-П;
реактор идеального смешения непрерывного действия РИС-Н;
реактор идеального вытеснения непрерывного действия РИВ;
каскад реакторов идеального смешения непрерывного действия
К-РИВ-Н (ячеечная модель).
Уравнение реактора любого типа (модели) получают из уравнения материального баланса, составленного по одному из веществ, участвующих в химической реакции (обычно по одному из реагентов)
11

Реактор идеального смешения периодического действия РИС-П
Аппарат с мешалкой, в который периодически загружают исходные реагенты и также периодически выгружают продукты.
Режим идеального смешения подразумевает, что в реакторе создается настолько интенсивное перемешивание, что в каждый момент времени концентрация реагентов одинакова по всему объему реактора и изменяется лишь во времени, по мере протекания химической реакции:
12

Графические зависимости:
1. Изменение концентрации реагента во времени СA=f(τ)
2. Изменение концентрации реагента по объему аппарата в данный момент времени (беря по абсциссе одну из координат трехмерного пространства x, y, z) СA=f(y)
3. Изменение скорости химической реакции по объему аппарата в данный момент времени (беря по абсциссе одну из координат трехмерного пространства x, y, z) rA=f(y)
4. Изменение степени превращения реагента по объему аппарата в данный момент времени (беря по абсциссе одну из координат трехмерного пространства x, y, z) αA=f(y)
13

14

15

При Vреакци.смеси = const CA = CA,0 (1-αА)
- характеристическое уравнение РИС-П
16

Предположим, что в реакторе протекает простая необратимая реакция «n»-го порядка
nA → продукты
По закону действующих масс скорость этой реакции равна произведению константы скорости и концентрации реагента в степени n:
Подставляя это выражение скорости реакции в характеристическое уравнение РИС-П
получаем:
Самостоятельно вывести уравнение для частного случая, когда порядок реакции равен 0
17

1. Для конкретного уравнения скорости химической
реакции строим графическую зависимость:
Графический расчет времени пребывания
2. Вычисляем S - площадь под кривой, «отсекаемую» по значениям абсцисс,
соответствующим начальной и конечной степеням превращения реагента А
(в приведенном примере на момент начала реакции αA,0 = 0).
3. Рассчитываем время пребывания, необходимое для достижении степени
превращения αA,к:
18
На рисунке степень превращения обозначена хАо

Реактор идеального смешения непрерывного действия РИС-Н
Аппарат с мешалкой, в который непрерывно подают реагенты и также непрерывно выводят из него продукты
Вход
СА,ВЫХ
Принимается, что подаваемое в реактор сырье с концентрацией реагента СА,0 мгновенно смешивается с реакционной массой, находящейся в реакторе и по всему объему реактора устанавливается одинаковая концентрация реагента, равная его концентрации на выходе СА,ВЫХ
СА,ВЫХ
СА,ВЫХ
СА,ВЫХ
СА,0 – концентрация в гипотетической точке пространства реактора с координатами x0, y0, z0
19

Графические зависимости:
Изменение концентрации реагента по объему аппарата для времени пребывания (реакции) τ (беря по абсциссе одну из координат трехмерного пространства x, y, z) СA=f(y)
Изменение концентрации реагента по объему аппарата для нескольких значений времени пребывания (реакции) τ1, τ2, τ3 (беря по абсциссе одну из координат трехмерного пространства x, y, z) СA=f(y)
(Учесть, что величина перепада между начальной и конечной концентрацией зависит при прочих равных условиях от величины скорости химической реакции и от времени пребывания реакционной смеси в реакторе. Чем выше скорость и больше время, тем больше скачок концентрации)
3. Изменение скорости химической реакции по объему аппарата (беря по абсциссе одну из координат трехмерного пространства x, y, z) rA=f(y)
4. Изменение степени превращения реагента по объему аппарата (беря по абсциссе одну из координат трехмерного пространства x, y, z) αA=f(y)
20

Как соотносятся τ1 τ2 τ3 ?
21

22

Поскольку для РИС-Н (работающего в стационарном режиме) характерно отсутствие градиента параметров и во времени, и по объему реактора – для вывода характеристического уравнения используют уравнение материального баланса в виде:
NА приход - NA расход = NА ХимРеакции
левая часть уравнения
NА приход - NA расход = V⋅ (СА,0 – СА, вых),
где V – объемный расход сырьевого потока, л/мин;
СА – концентрация, моль/л
правая часть уравнения NА ХимРеакцц = rA⋅VР,
где VР – рабочий объем реактора, л;
rA – скорость химической реакции, моль/(л⋅мин)
23

- характеристическое уравнение для РИС-Н
24

Предположим, что в реакторе протекает простая необратимая реакция «n»-го порядка nA → продукты
По закону действующих масс скорость этой реакции равна произведению константы скорости и концентрации реагента в степени n:
Подставляя это выражение скорости реакции в характеристическое уравнение РИС-Н
получаем:
Самостоятельно вывести уравнение для частного случая, когда порядок реакции равен 0
25

Реактор идеального вытеснения непрерывного действия РИВ
К режиму РИВ близок трубчатый аппарат, в котором отношение длины трубы к ее диаметру L/d > 20
На вход реактора непрерывно подаются реагенты, которые превращаются в продукты реакции по мере перемещения их по длине реактора.
Любая частица потока движется прямолинейно и равномерно, не происходит ни продольного, ни поперечного перемешивания, то есть время пребывания каждой частицы в реакторе одинаково. Таким образом, модель РИВ предполагает движение реакционной смеси только в одном направлении – по х – координате, то есть по длине реактора L.
26

Графические зависимости:
Изменение концентрации реагента по длине реактора L (абсцисса соответствует координате трехмерного пространства x)
Изменение скорости химической реакции по длине реактора L
Изменение степени превращения реагента по длине реактора L
27

Так как изменение концентрации происходит только в одном направлении - по длине реактора, а составляющие конвективного переноса по осям y и z равны 0 имеем:
где W – линейная скорость движения реакционной смеси, L – длина реактора
Важно: с учетом малой величины скорости молекулярной диффузии по сравнению со скоростью конвективного переноса, диффузионная составляющая дифференциального уравнения материального баланса может быть приравнена к нулю:
28

С учетом вышеизложенного получаем для реактора РИВ
уравнение материального баланса в виде:
Как было отмечено выше, для реактора, работающего в стационарном режиме,
концентрация в любой точке аппарата является величиной, постоянной во времени, то есть:
или
29

Линейная скорость движения реакционной смеси в реакторе с постоянной площадью сечения равна
С учетом этого уравнение материального баланса преобразуется как:
или
После разделения переменных и интегрирования получаем:
- характеристическое уравнение для РИВ, работающего в стационарном режиме
30

Как видно, характеристические уравнения материального баланса для РИС-П и РИВ одинаковы, то есть время достижения заданной степени превращения реагента в этих реакторах одно и то же.
Но в РИС-П полное время процесса складывается не только из рабочего времени, рассчитываемого по характеристическому уравнению, но и из времени на вспомогательные операции - загрузку реагентов, выхода на рабочую температуру, охлаждение реакционной массы по окончании реакции, выгрузку реакционной массы, подготовку реактора к следующему синтезу
В РИВ вспомогательных операций нет, поэтому интенсивность РИВ выше, чем у РИС-П.
- характеристическое уравнение РИС-П
Сравнение РИС-П и РИВ
- характеристическое уравнение РИВ
31

Сравнение РИС-Н и РИВ
В РИС-Н вспомогательных операций нет, как и в случае РИВ, но гидродинамическая обстановка в реакторах различная.
Проведем сравнение этих реакторов исходя из изменения концентрации реагента и скорости химической реакции:
в РИВ-Н наблюдается постепенное уменьшение концентрации реагента по длине реактора (см. соответствующий график), а в РИС-Н – резкое падение концентрации от начального до конечного значения (см. соответствующий график).
Поскольку скорость реакции прямо пропорциональна концентрации реагента, профили скорости имеют вид такой же, как и профили концентрации (см. соответствующие графики).
Вывод: поскольку в РИВ в среднем концентрация реагента и скорость реакции выше, чем в РИС-Н, интенсивность РИВ выше чем РИС-Н.
Обратить внимание: для реакции псевдонулевого порядка, когда скорость не зависит от концентрации, различие между реакторами РИВ и РИС-Н вырождается.
32

Реакция гидрирования боковой цепи альфа-метилстирола исследовалась в реакторе смешения периодического действия (см. рисунок)

Общий вид кинетического уравнения
в соответствии с законом действующих масс
где kэфф. - эффективная константа скорости, л/(г⋅мин), размерность массы в знаменателе соответствует массе катализатора;
Скатализатор - “концентрация“ катализатора ( отношение массы катализатора к объему жидкой фазы, г/л );
a, b - наблюдаемые порядки по олефину и водороду;
с - показатель степени при “концентрации” катализатора.

Кинетические кривые изменения концентрации реагента и продукта реакции
(mкат. = 75 г; t = 133 °С; P = 2,0 МПа)

Определение скорости реакции по линейным участкам кинетических кривых (те же данные, что и для предыдущего рисунка)
Исходя из вида кинетической кривой расходования олефина dC/dτ=const
до остаточной концентрации порядка 0,3 моль/л скорость реакции не зависит от концентрации олефина, то есть наблюдаемый порядок по реагенту равен 0

Растворимость водорода в АМСТ и ИПБ в зависимости от температуры; давление PΣ = 3,0 МПа
Полученные экспериментальные результаты по растворимости для температур 60 ÷ 150°С и давлений 4 ÷ 30 атм (0,4 ÷ 3,0 МПа) обобщены в виде температурных зависимостей коэффициента Н/, аналогичного коэффициенту Генри (Н), но без учета давления паров органических компонентов:
для альфа-метилстирола Н/ = 5,31*ехр(629,8/Т)
для изопропилбензола Н/ = 2,81*ехр(781,6/Т)

Растворимость водорода (дм3/кг) в присутствии и отсутствии
катализатора ( t = 149 °С, сырье 100%-й АМСТ )
Вывод: в исследуемых условиях, не зависимо от протекания химической реакции, концентрация водорода в жидкой фазе может быть приравнена к равновесной концентрации водорода в системе водород (газ) – растворитель (жидкость)

К определению порядка по водороду: логарифм скорости гидрирования альфа-метилстирола в зависимости от логарифма концентрации водорода в жидкой фазе (mкат. = 75 г; t = 101°С; Р = 1÷3 МПа)
Вывод: порядок по водороду ~ 1

К определению порядка по катализатору: логарифм скорости гидрирования альфа-метилстирола в зависимости от логарифма концентрации катализатора в жидкой фазе (Скат. = 125 – 250 г/л; t = 101°С; Р = 1÷3 МПа)

Итого:
порядок по олефину принимаем равным 0
порядок по водороду принимаем равным 1
порядок по катализатору принимаем равным 1
и преобразуем кинетическое уравнение к виду:

Как определяются параметры уравнения Аррениуса
kэфф,0 и ЕА
Проводится серия экспериментов в изотермических условиях, когда переменной величиной в серии является только температура реакции, начальная концентрация реагента, общее давление, концентрация катализатора в каждом эксперименте одинаковы.
Для каждого эксперимента, выполненного при отдельной температуре определяется численное значение скорости реакции.
Исходя из вида кинетического уравнения по значению скорости химической реакции определяется численное значение эффективной константы скорости для каждого значения температуры реакции, для нашего примера это:
4. После выполнения п. 3 мы получаем дискретную зависимость
kэфф = f(температура)
5. Преобразуем зависимость kэфф = f(температура) в координатах
ln kэфф = f(1000/T)

Температурная зависимость логарифма k эфф.
(Cкат. = 125 г/л; t = 80 ÷155 °С; РΣ = 2,0 МПа)
Зависимость должна быть линейна, то есть аппроксимироваться уравнением
y = a+b⋅x
Угловой коэффициент данной зависимости b равен величине «- ЕА/R»
коэффициент а натуральному логарифму предэкспоненциального множителя kэфф,0

Итоговое кинетическое уравнение скорости химической реакции: