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

книги / Теория химических реакторов введение в основной курс

..pdf
Скачиваний:
6
Добавлен:
12.11.2023
Размер:
5.34 Mб
Скачать
– адиабатное повышение температуры ре-

2.

CAi 1

CAi

CAi

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

цикле.

 

 

 

 

3.

Ci 1

ki 1

CAi 1

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

втекущем цикле.

4.Qрi 1 Ci 1 ( Hr ) – выделившееся (поглотившееся)

тепло в текущем цикле.

5. TAi 1 Qрi 1

N Cp

акционной массы в цикле.

6. Ti 1 Ti TAi 1 – теоретическая максимальная температу-

ра реакционной массы в цикле.

 

7.

T

 

mp Txi Tx b

 

– средняя температура хладагента

 

 

xi 1

 

mp b

 

 

 

 

 

врубашке охлаждения в данном цикле.

8.TTOi 1 Ti 1 Txi 1 – теоретическая максимальная разность

температур на стенке аппарата в цикле.

9. QTOi 1 K F TTOi 1 – количество теплоты, переданное через стенку от реакционной массы к теплоносителю.

10.

Txi

1 mp

b Cp.x

– среднее повышение темпера-

 

 

 

 

QTOi 1

 

туры теплоносителя.

 

 

11.

Txi 1

Txi 1 Txi 1 – средняя температура хладагента в

рубашке охлаждения для (i+2)-го цикла.

12.Qpi 1 Qpi 1 Qтоi 1 – количество теплоты, оставшееся

вреакционной массе после теплообмена в цикле.

13.

TAi 1

Qpi 1

– фактическое повышение температуры

 

 

 

N Cp

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

161

14.

 

 

 

– расчетная температура в реакторе (ре-

Ti 1

Ti TAi 1

акционной массы) к (i+2)-му циклу расчета.

15.

 

 

 

 

TTOi 1 TTOi 1

1/2 TAi 1 Txi 1 – уточненное значе-

ние разности температур на стенке аппарата (сравнение с TTO – не более 1 % отклонения).

16. QTOi 1 K F TTOi 1 – уточненное значение количества тепла, переданного через стенку к хладагенту.

 

17.

Txi +1

 

 

 

 

 

 

 

 

mp b Cр.x .

 

 

 

 

 

 

QTOi 1

 

 

цикл

 

 

 

 

 

 

18.

 

1 Qрi 1

 

– уточненное количество

 

Внутренний

Qрi

QTOi 1

тепла, оставшееся в реакторе.

 

 

 

 

19.

TAi 1

 

Qpi 1

 

– уточненное изменение температу-

 

N Cp

 

 

 

 

 

 

 

 

ры в реакторе, с учетом уточненного теплообмена.

20. Ti 1 Ti TAi 1 – уточненная расчетная температура в реакторе к (i+2)-му циклу.

 

TTOi 1

 

 

– сравнение

21. TTOi 1

1/ 2 TAi 1

Txi 1

сTTOi 1 (не более 1 % отклонение) (цикл) к № 16.

4.Вычисления в Mathcad. Рассмотрим учебный пример численного расчета химической реакции n-го порядка, протекающей на катализаторе в реакторе, выполненном в виде трубы диаметром d = 200 мм, длиной L = 10 м. Катализатор нанесен на поверхность керамических шариков радиусом r = 2,5 мм, которые пол-

ностью заполняют объем реактора.

Примем удельный вес катализатора γk = 3000 кг/м3; средний удельный вес смеси жидких реагентов – 800 кг/м3. Cpk = 12;

Cpg = 8 – теплоемкостькатализатораижидкости, соответственно. Начальная концентрация Cab = 4 кмоль/м3; ∆Hr – тепловой

эффект реакции, Дж/моль; M – средняя молекулярная масса реагентов; S – полная поверхность катализатора; INT – количество

162

расчетных интервалов (секций), на которые делится длина химического реактора; n – порядок химической реакции.

При вычислении на ЭВМ предварительно определяются некоторые параметры: масса и теплоемкость катализатора в секции, масса и теплоемкость жидкости, заполняющей свободный объем в секции; длительность пребывания реагентов в секции при заданной производительности реактора Vсек (Vsek), а также параметры α и β, соответствующие долям распределения выделяющейся теплоты между жидкой и твердой фазами соответственно:

163

Программы и примеры расчета распределения концентрации исходного реагента и температуры по секциям реактора представлены на рис. 54.

Рис. 54. Программа расчета реактора по секциям

164

Программа расчета концентрации СА и температуры Т отличаются друг от друга только блоком вывода результатов. Как было указано ранее, любой параметр, участвующий в расчете внутри программы, может быть выведен в качестве результата за счет правильной организации блока вывода.

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

Преобразуем программу расчета таким образом, чтобы все вычисления касались только этой секции. Для этого в операторе расчета ∆С заменим текущую концентрацию реагента CA на его исходное значение CAb. В расчете принято, что каждая новая порция сырья находится в первой секции ∆τ с, после чего она вытесняется в следующую секцию. Величина m – количество поданных порций сырья в первую секцию реактора (рис. 55). Вычисления показывают, что после 30-й порции расчетная температура в секции достигает астрономической величины. Причиной этого является невозможно большая величина падения концентрации реагента ∆С в первой секции, равная примерно 8000 кмоль/м3 при его фактической начальной концентрации4 кмоль/м.

В связи с этим в программу должны быть внесены ограничения по изменению концентрации. На рис. 56 во внутреннем цикле программы добавлен оператор «if», ограничивающий величину ∆С в первой секции. Результаты расчета представлены на рис. 57.

Равновесная температура в секции около 780 K, достигаемая примерно начиная с 50-й порции сырья, может оказаться неприемлемой с точки зрения технологии процесса. В связи с этим необходимо ввести в систему теплообменник, предназначенный для съема избыточного тепла химической реакции в первой секции. На рис. 57 представлена программа расчета секции со встроенным или наружным теплообменником. Величина F1 (м2) соответствует величине теплообменной поверхности, приходящейся на первую секцию; K – коэффициент теплопередачи (рис. 58).

165

Рис. 55. Программа расчета первой секции реактора

166

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

167

Рис. 57. Программа расчета реактора со встроенным теплообменником

168

Рис. 58. Влияние величины поверхности теплообмена F1 на технологические параметры работы химического реактора:

индекс 1: K = 0,6; F1 = 0,1; индекс 2: K = 0,6; F1 = 0,2

Реактор с регулятором [11]. В реакторе периодического действия протекает экзотермическая реакция конденсации А+В→С. Процесс проводится в аппарате, снабженном перемешивающим устройством и рубашкой охлаждения. Регулирования температуры осуществляется за счет изменения подачи реагента В, при этом компонент А заранее помещается в реактор в необходимом количестве. Расход охлаждающего агента через рубашку охлаждения постоянный. Температура реакционной массы измеряется термометром сопротивления, сигнал от которого передается на регулирующее устройство, управляющееподачейреагентаВ.

Система уравнений материального баланса по каждому веществу:

 

d(V CA )

k CA CB V ;

 

 

(217)

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

d(V CB )

w

 

C

Bin

k C

A

C

B

V ;

(218)

 

 

dt

 

Bin

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

169

d(V CC )

 

k CA CB V ;

(219)

dt

 

 

 

 

 

d (V )

w

,

(220)

 

 

 

dt

 

Bin

 

 

 

 

 

 

 

где СА, СB, СС – концентрации веществ в реакционной зоне; СBin – концентрация подаваемого вещества В; V – объем реактора, м3; wВin – объемныйрасходвеществаВвовходномпотоке, м3/ч.

Уравнения теплового баланса:

– для зоны реакции

dT

 

w

 

K F

 

 

Qр

 

 

 

 

 

Bin

(T

T)

m

(T T

)

 

k C

C ;

(221)

dt

 

V C

C

 

V in

 

ст

 

A

B

 

 

 

 

 

 

p

 

 

p

 

 

 

– для рубашки охлаждения

dTст

wx

(T

T

)

Kтп F

(T T

),

(222)

 

 

dt

 

x,0

ст

 

V Cp.ст

ст

 

 

Vохл

 

 

 

 

 

где Tin, T, Tст – температуры входного потока, в реакторе и стенки реактора, К; wx – объемный расход хладагента, м3/ч; F – площадь поверхности охлаждения, м2; Vохл – объем охлаждающей поверхности; Kтп – коэффициент теплопередачи, кВт/м2/К; ρ – плотность массы в реакторе, кг/м3; Ср, Ср.ст – удельная теплоемкость реагирующей массы и стенки реактора, кДж/(кг·К); Qр – тепловой эффект реакции, кДж/моль.

Уравнение регулятора

dw

 

 

dT

 

Kp

 

 

 

Bin K

 

 

 

 

 

(T

T ),

(223)

 

dt

tизодр

dt

p

 

 

зад

 

 

где Kр – коэффициент пропорциональности регулятора, м3/(ч·К);

tизодр – времяизодрома, ч;

Программа расчета и графики изменения основных параметров реактора от длительности процесса приведены на рис. 59.

170

Соседние файлы в папке книги