Режимы работы реактора: стационарный (установившийся) и нестационарный При

любой точке реактора во времени не происходит изменения концентрации и

температуры

При нестационарном режиме концентрация и температура в любой точке реактора меняются во времени:

Стационарный режим реализуется в реакторах непрерывного действия.
Нестационарный в периодических реакторах, и в непрерывных реакторах во время пуска или остановки реактора.

1

- уравнение теплового баланса

В общем виде уравнение материального баланса

применительно к реагенту можно представить в виде:
количество вещества, поступающего на вход в реактор в единицу времени – расход вещества в единицу времени (на химическую реакцию, а также выходящего из реактора) = накопление (убыль) вещества в реакторе в единицу времени

По аналогии уравнение теплового баланса в общем виде:
приход тепла в реактор в единицу времени – расход тепла в единицу времени = накопление (убыль) тепла в реакторе в единицу времени

Для реактора, работающего в стационарном режиме, когда концентрация и температура во времени не меняются, правая часть обоих уравнений равна 0.

2

по компоненту
Элементный баланс
3

и выходе из аппарата (составляется по всем компонентам системы)
4

– кипятильник, 4, 20 – роторно-плёночные испарители,
7, 8 –

реакторы, 9 – трубчатая печь, 10 – скруббер, 11 – воздушные конденсаторы, 12 – водяной конденсатор,
13 – отстойник, 14, 16, 18 – ректификационные колонны, 15, 17, 19 – конденсаторы.

5

продукту) составляется с учетом расходования или накопления вещества в результате

химической реакции:

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

7

по углероду для органических смесей
по азоту
по сере

по кислороду

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

8

химическая реакция:
А  В
где СА – концентрация реагента А

в реакционной смеси;
x,y,z – пространственные координаты;
Wx, Wy, Wz – составляющие скорости потока по осям х, у, z;
D – коэффициент диффузии вещества А;
rA – скорость химической реакции

Левая часть уравнения

характеризует изменение концентрации вещества

А во времени в элементарном объеме dV

Уравнение конвективного массообмена, дополненное слагаемым rA, которое
учитывает протекание химической реакции

9

отражает изменение концентрации в результате переноса вещества
путем молекулярной диффузии

слагаемое rA показывает изменение концентрации за счет химической реакции

10

или меньшем приближении можно описать одной из следующих моделей:
реактор идеального

смешения периодического действия РИС-П;
реактор идеального смешения непрерывного действия РИС-Н;
реактор идеального вытеснения непрерывного действия РИВ;
каскад реакторов идеального смешения непрерывного действия
К-РИВ-Н (ячеечная модель).

Уравнение реактора любого типа (модели) получают из уравнения материального баланса, составленного по одному из веществ, участвующих в химической реакции (обычно по одному из реагентов)

11

который периодически загружают исходные реагенты и также периодически выгружают продукты.
Режим

идеального смешения подразумевает, что в реакторе создается настолько интенсивное перемешивание, что в каждый момент времени концентрация реагентов одинакова по всему объему реактора и изменяется лишь во времени, по мере протекания химической реакции:

12

реагента по объему аппарата в данный момент времени (беря по

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

13

= CA,0 (1-αА)

- характеристическое уравнение РИС-П

16

nA  продукты
По закону действующих масс скорость этой реакции равна

произведению константы скорости и концентрации реагента в степени n:

Подставляя это выражение скорости реакции в характеристическое уравнение РИС-П
получаем:

Самостоятельно вывести уравнение для частного случая, когда порядок реакции равен 0

17

расчет времени пребывания
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

параметров и во времени, и по объему реактора – для

вывода характеристического уравнения используют уравнение материального баланса в виде:
NА приход - NA расход = NА ХимРеакции

левая часть уравнения
NА приход - NA расход = V (СА,0 – СА, вых),
где V – объемный расход сырьевого потока, л/мин;
СА – концентрация, моль/л

правая часть уравнения NА ХимРеакцц = rAVР,
где VР – рабочий объем реактора, л;
rA – скорость химической реакции, моль/(лмин)

23

nA  продукты
По закону действующих масс скорость этой реакции равна

произведению константы скорости и концентрации реагента в степени n:

Подставляя это выражение скорости реакции в характеристическое уравнение РИС-Н
получаем:

Самостоятельно вывести уравнение для частного случая, когда порядок реакции равен 0

25

аппарат, в котором отношение длины трубы к ее диаметру L/d

> 20

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

26

координате трехмерного пространства x)

Изменение скорости химической реакции по длине

реактора L

Изменение степени превращения реагента по длине реактора L

27

по длине реактора, а составляющие конвективного переноса по осям y

и z равны 0 имеем:

где W – линейная скорость движения реакционной смеси, L – длина реактора

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

28

в виде:
Как было отмечено выше, для реактора, работающего в стационарном

режиме,
концентрация в любой точке аппарата является величиной, постоянной во времени, то есть:

или

29

сечения равна
С учетом этого уравнение материального баланса преобразуется как:
или
После разделения

переменных и интегрирования получаем:

- характеристическое уравнение для РИВ, работающего в стационарном режиме

30

одинаковы, то есть время достижения заданной степени превращения реагента в

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

- характеристическое уравнение РИС-П

Сравнение РИС-П и РИВ

- характеристическое уравнение РИВ

31

в случае РИВ, но гидродинамическая обстановка в реакторах различная.
Проведем

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

32

действия (см. рисунок)

- эффективная константа скорости, л/(гмин), размерность массы в знаменателе соответствует

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

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Т)

= 149 С, сырье 100%-й АМСТ )
Вывод: в исследуемых условиях,

не зависимо от протекания химической реакции, концентрация водорода в жидкой фазе может быть приравнена к равновесной концентрации водорода в системе водород (газ) – растворитель (жидкость)

зависимости от логарифма концентрации водорода в жидкой фазе (mкат. =

75 г; t = 101С; Р = 13 МПа)
Вывод: порядок по водороду  1

зависимости от логарифма концентрации катализатора в жидкой фазе (Скат. =

125 – 250 г/л; t = 101С; Р = 13 МПа)

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

к виду:

изотермических условиях, когда переменной величиной в серии является только температура

реакции, начальная концентрация реагента, общее давление, концентрация катализатора в каждом эксперименте одинаковы.
Для каждого эксперимента, выполненного при отдельной температуре определяется численное значение скорости реакции.
Исходя из вида кинетического уравнения по значению скорости химической реакции определяется численное значение эффективной константы скорости для каждого значения температуры реакции, для нашего примера это:

4. После выполнения п. 3 мы получаем дискретную зависимость
kэфф = f(температура)
5. Преобразуем зависимость kэфф = f(температура) в координатах
ln kэфф = f(1000/T)

= 80 155 С; Р = 2,0 МПа)
Зависимость должна быть

линейна, то есть аппроксимироваться уравнением
y = a+bx
Угловой коэффициент данной зависимости b равен величине «- ЕА/R»
коэффициент а натуральному логарифму предэкспоненциального множителя kэфф,0

в ходе выполнения эксперимента.
В нашем примере это:
диапазон температур 80 –

155 С,
диапазон давлений 1 – 3 МПа,
диапазон концентраций катализатора 125 – 250 г/л

То есть мы можем вполне корректно рассчитать скорость реакции
для следующих рабочих условий:
t =123 C
P = 2,25 МПа
Скат = 200 г/л

Но не корректна экстраполяция за пределы варьирования, например расчет скорости для следующих рабочих условий:
t =123 C
P = 2,25 МПа
Скат = 280 г/л
и т.п.