Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Химические реакторы.-1

.pdf
Скачиваний:
4
Добавлен:
20.11.2023
Размер:
2.16 Mб
Скачать

Значение τ по этому уравнению является результатом усреднения времени пребывания в реакторе различных элементов потока. В этой связи возникает вопрос: обязателен ли учёт неравномерности времени пребывания элементов потока в реакторе по уравнению (6.2)?

Практика показывает, что при одном и том же среднем времени пребывания различные степени неравномерности приводят к разным результатам. Так, элементы потока с малым временем пребывания выносят из реактора много непрореагировавшего вещества, что снижает степень превращения. С другой стороны, в элементах, долго находившихся в реакторе, реакция протекает более глубоко. Но, учитывая, что скорость химической реакции снижается со временем, увеличение времени пребывания слабо сказывается на химическом превращении. Проигрыш за счёт быстро «проскочивших» через реактор элементов потока не компенсируется за счёт элементов, задержавшихся в нём.

На практике удобнее использовать безразмерное время пребывания Θ и безразмерную функцию распределения с(Θ):

 

 

t

;

 

 

(6.3)

 

 

 

 

 

 

 

 

τ

 

 

 

с( )

сi (t)

,

(6.4)

 

 

 

 

 

 

 

 

с0

 

 

 

 

 

 

 

 

i

 

 

где с0

– начальная концентрация индикатора в потоке,

 

i

 

 

 

 

 

 

 

 

 

 

с0

 

g

.

 

(6.5)

 

 

 

 

i

VR

 

 

 

 

 

 

Безразмерная функция распределения связана со средним временем пребывания и экспериментальным значением функции

отклика сi (t) следующей зависимостью:

 

с( ) с(t),

(6.6)

где выражение

 

 

с(t)

сi (t)

(6.7)

 

 

 

 

сi (t) dt

 

 

0

 

задает нормированную кривую отклика.

91

6.2. Обработка результатов экспериментального исследования структуры потоков

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

Расчётная формула метода Симпсона имеет следующий вид:

b

x f1 4 f

 

 

 

f x dx

2 2 f3 4 f

4 2 f5 4 fn 1 fn ,

(6.8)

a

3

 

 

 

 

 

 

 

где x – интервал между соседними значениями времени (интервал интегрирования); f1, f2, f3, …, fn – последовательные значения подынтегральных функций.

Использование метода Симпсона предполагает, что интервал интегрирования разбит на чётное количество отрезков (число точек при этом будет нечётным!), поэтому при обработке экспериментальных данных на с-кривой следует брать соответствующее количество точек.

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

Пример обработки с-кривой. Пусть в ходе опытов были получены следующие данные:

t, с

0

10

20

30

40

50

60

70

80

90

100

ci, усл. ед.

0

0

1

4

11

10

7

4

2

0,5

0

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

t

0

10

20

30

40

50

60

70

80

90

100

ci

0

0

1

4

11

10

7

4

2

0,5

0

t·ci

0

0

20

120

440

500

420

280

160

45

0

С учетом формул (6.2) и (6.8) находим среднее время пребывания потока в реакторе:

92

10 0 4 0 2 20 4 120 2 440 4 500 2 420 4 280 2 160 4 45

τ 3 ;

103 0 4 0 2 1 4 4 2 11 4 10 2 7 4 4 2 2 4 0,5 2 0

τ5860116 50,5 с.

6.3.Модели реальных реакторов

Втеории химических реакторов разработаны различные модели, позволяющие учесть неидеальность потока. Наибольшее распространение получили ячеечная и диффузионная модели.

Основой ячеечной модели является представление об идеальном перемешивании в пределах ячеек, расположенных последовательно, и отсутствии перемешивания между ячейками.

Параметром, характеризующим модель, служит число ячеек n. Функция распределения ячеечной модели имеет вид:

с Θ

nn

exp n Θ .

 

 

Θn 1

(6.9)

 

 

n 1 !

 

 

При n = 1 ячеечная модель переходит в модель идеального смешения, а при n → ∞ – в модель идеального вытеснения (рис. 19).

Рис. 19. Функция распределения c(Θ) ячеечной модели в зависимости от значения параметра n

при импульсном возмущении

93

Обычно для определения параметров моделей (идентификации) используют дисперсию времени пребывания элементов потока в аппарате σΘ2 .

Дисперсия, или разброс случайной величины около среднего значения при непрерывном распределении для конечного числа эквидистантных точек (равные интервалы отбора проб), определяется по формуле:

 

 

σΘ2 ( 1)2 c( ) d .

(6.10)

0

 

Уравнение связи параметра ячеечной модели n с безразмерной дисперсией функции отклика модели на импульсное возмущение имеет вид:

n

1

.

(6.11)

 

 

2

 

 

 

 

В основе диффузионной модели лежит допущение от том, что при движении элементов потока отклонение от идеального вытеснения осуществляется за счёт механизма, аналогичного диффузионному. Это позволяет учитывать перемешивание с помощью коэффициента продольного перемешивания, или продольной диффузии DL (рис. 20). Причём коэффициент продольного перемешивания не тождественен коэффициентам турбулентной DT или молекулярной D диффузии, так как он учитывает отклонение от идеального вытеснения за счёт всех возможных механизмов в сумме (конвективного, молекулярного и турбулентного).

Рис. 20. Схема диффузионной модели

94

В настоящее время нет точных аналитических методик определения коэффициента продольного перемешивания DL. Его можно определить экспериментально с достаточной степенью точности через дисперсию времени пребывания элементов потока в аппарате σΘ2 :

σΘ2

2

Pe 1 e Pe ,

(6.12)

2

 

Pe

 

где Pe диффузионный критерий Пекле,

 

 

Pe

wl

.

(6.13)

 

 

 

 

 

DL

 

При значениях критерия Pe более 10 упрощённо можно принять:

σΘ2 Pe2 . (6.14)

При Pe = 0 (DL ) поток соответствует идеальному смеше-

нию. При Pe поток движется по схеме идеального вытеснения

(рис. 21).

Рис. 21. Функция распределения с(Θ) диффузионной модели в зависимости от значения критерия Пекле при импульсном возмущении

Учёт продольного перемешивания в реальном реакторе достигается введением в правую часть уравнения (5.21) дополнительного члена с коэффициентом продольного перемешивания:

95

с

с

 

2с

 

 

A w

A

D

A

r .

(6.15)

l

l2

τ

L

A

 

Уравнение (6.15) – это однопараметрическая диффузионная модель реального реактора.

Степень отклонения реального реактора от идеального может

быть выражена через соотношение их объёмов Vр , необходимых

Vид

для достижения одинаковой степени превращения ХА, в виде зависимости (рис. 22):

 

 

 

Vр

 

 

wl

 

 

 

 

 

 

 

 

 

 

 

f

 

, X A .

 

 

 

(6.16)

 

 

Vид

DL

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из рис. 22 следует, что

 

 

 

 

 

 

 

 

уменьшение комплекса

 

wl

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

DL

 

 

 

 

 

 

 

 

т.е. возрастание степени от-

 

 

 

 

 

 

 

 

клонения

реального

реактора

 

 

 

 

 

 

 

 

от идеального, влечёт за собой

 

 

 

 

 

 

 

 

необходимость

увеличения

 

 

 

 

 

 

 

 

объёма

реального

реактора

 

 

 

 

 

 

 

 

для достижения заданной сте-

Рис. 22. Зависимость

 

 

 

 

пени превращения, и эта тен-

 

 

 

 

денция растёт с увеличением

 

Vр

 

 

 

 

 

 

отношения

от степени

 

 

 

значения XA.

 

 

 

 

 

Vид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для учёта не только про-

 

 

 

 

 

 

 

 

превращения при различных

wl

 

дольного, но и

радиального

 

 

 

DL

 

перемешивания

в

реакторе

 

 

 

 

 

 

 

 

 

 

 

 

 

 

поступают так же,

как и при

выводе однопараметрической модели, но в этом случае необходимо ввести второй член в уравнение (5.21), отражающий перемешивание в радиальном направлении. При этом уравнение принимает вид:

96

c

c

 

2c

 

2c

 

A w

A

D

A

D

A

r ,

(6.17)

l

l2

 

τ

L

R R2 A

 

где DR – коэффициент радиального перемешивания; R – радиус ре-

актора.

Уравнение (6.17) – двухпараметрическая диффузионная модель реального реактора.

Выводы

Типовые модели идеальных реакторов (РИВ-Н, РИС-Н) предполагают идеальную структуру потоков и не всегда адекватны реальным процессам.

Различные траектории и скорости движения элементов потока в химическом реакторе обусловливают их различное время пребывания в нём, что оказывает большое влияние на химические превращения.

Степень отклонения реального реактора от идеального определяется экспериментально путём анализа функций отклика на стандартные возмущения потока с последующим синтезом математической модели реактора. Наибольшее распространение получили ячеечная и диффузионная модели.

Функция распределения ячеечной модели имеет вид:

с Θ

nn

exp n Θ ,

 

Θn 1

 

 

n 1 !

 

где n – число ячеек – основной параметр ячеечной модели.

При n = 1 ячеечная модель переходит в модель идеального смешения, при n → ∞ – в модель идеального вытеснения.

Число ячеек n связано с дисперсией 2 функции отклика на

импульсное возмущение уравнением

n 12 .

97

Основным параметром однопараметрической диффузионной модели является коэффициент продольного перемешивания DL в составе диффузионного критерия Пекле:

Pe wl .

DL

Критерий Пекле может быть определён следующим образом:

σΘ2 Pe22 Pe 1 e Pe .

При Pe = 0 (DL ) поток соответствует идеальному смешению.

При Pe поток движется по схеме идеального вытеснения.

В двухпараметрической диффузионной модели учитывается коэффициент радиального перемешивания DR .

Контрольные вопросы

1.Что понимают под структурой потоков в реакторе?

2.Для чего применяются знания о структуре потоков в реак-

торах?

3.В чём заключается суть методики исследования структуры потоков в реакторе?

4.Каковы типовые математические модели структуры потоков

вхимических реакторах?

5.Какие основные параметры характеризуют ячеечную и диффузионную модели?

98

7. ТЕПЛОВЫЕ РЕЖИМЫ РЕАКТОРОВ

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

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

7.1.Классификация реакторов

Взависимости от теплового режима реакторы подразделяют на три группы:

1) адиабатические;

2) изотермические;

3) политропические.

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

Визотермическом реакторе путём подвода или отвода тепла поддерживают постоянную температуру в течение всего процесса.

Вполитропическом реакторе температура непостоянна, при этом часть тепла подводится или отводится.

7.2. Уравнение теплового баланса реактора

Основой для расчёта реакторов с учётом теплового режима служит дифференциальное уравнение конвективного теплообмена,

вкоторое вводят дополнительные члены, учитывающие отвод тепла

врезультате теплообмена и тепло реакции:

99

 

 

T

ρCP

 

 

T

wy

T

wz

T

 

 

 

ρCP

t

wx

x

y

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2T

 

2T

 

2T

 

 

 

 

 

 

 

(7.1)

λ

x

2

y

2

z

2

 

F K T r H R ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где ρ – плотность реакционной массы; СР – удельная теплоёмкость реакционной массы; х, у, z – пространственные координаты; wx, wy, wz – составляющие скорости движения потока в направлении осей; λ – коэффициент теплопроводности реакционной массы; F' – удельная поверхность теплообмена (на единицу объёма реакционной массы); K – коэффициент теплопередачи; Т – разность температур реакционной массы и хладоагента,

Т T Tхл ,

Т – температура реакционной массы; Tхл – температура хладоаген-

та; HR – тепловой эффект реакции; r – скорость реакции. Уравнение (7.1) – общее уравнение теплового баланса реактора. Левая часть этого уравнения выражает скорость накопления те-

пла в любой точке, для которой составляется тепловой баланс. Группа членов правой части уравнения характеризует измене-

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

Примем, что изменение температуры происходит только в одном направлении по длине реактора l (или по оси x), тогда:

T

 

T

;

T

 

T

0.

(7.2)

x

 

l

 

y

 

z

 

 

Изменение температуры в реакторе за счёт теплопроводности обычно мало и его можно не учитывать, тогда:

 

2T

 

2T

 

2T

0.

 

λ

x

2

y

2

z

2

 

(7.3)

 

 

 

 

 

 

 

 

 

С учётом указанных упрощений уравнение теплового баланса (7.1) примет вид:

100