ФЕДЕРАЛЬНОЕ АГЕНТСТВО ПО ОБРАЗОВАНИЮ
ВОЛГОГРАДСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ
КАФЕДРА «ПРОМЫШЛЕННАЯ ЭКОЛОГИЯ И
БЕЗОПАСНОСТЬ ЖИЗНЕДЕЯТЕЛЬНОСТИ»
Расчет на эвм экзотермического реактора идеального смешения
Методические указания к лабораторной работе
РПК
«Политехник»
Волгоград
2006
УДК 536.629
Расчет на ЭВМ экзотермического реактора идеального смешения: метод. указ. к лабораторной работе / сост.: А. Б. Голованчиков, Н. А. Дулькина, В. А. Козловцев, А. А. Шагарова; Волгоград. гос. техн. ун-т. – Волгоград, 2006. – 18 с.
Приведен алгоритм расчета экзотермического реактора идеального смешения. Дана таблица идентификаторов с исходными и справочными данными и расчетными параметрами.
Работа выполняется в процессе изучения дисциплины «Математическое моделирование стационарных систем в химической технологии».
Предназначены для студентов дневной и вечерней форм обучения по направлениям 550800 «Химическая технология и биотехнология» и 551800 «Технологические машины и оборудование».
Ил. 4. Табл.3. Библиогр.: 7 назв.
Рецензент Г.В. Рябчук
Печатается по решению редакционно-издательского совета
Волгоградского государственного технического университета
© Волгоградский государственный
технический университет, 2006
Введение
Экзотермические реакторы идеального смешения широко применяются в химической и нефтехимической промышленности.
Основной проблемой при их расчете является определение поверхности рубашки и при необходимости поверхности змеевика, обеспечивающих отвод тепла реакции от реакционной массы к хладагенту [1-3].
Цель работы
Изучить алгоритм расчета экзотермического реактора идеального смешения.
Научить студентов выбирать справочные данные по учебной и научной литературе и выбирать стандартные реакторы по результатам расчетов.
Уметь иллюстрировать проводимые расчеты в виде схем и графиков.
Делать выводы по результатам расчетов на ЭВМ.
Постановка задачи
Рассчитать технологические параметры и геометрические размеры экзотермического реактора идеального смешения по заданной производительности, степени конверсии основного реагирующего компонента, тепловому эффекту реакции и зависимости концентрации основного реагирующего компонента от времени (интегральной кинетической кривой, рис.1) [4-6]. Исходные и необходимые справочные данные приведены в таблице 2.
Рис.1.
Интегральная кинетическая кривая
зависимости концентрации основного
исходного компонента А
от времени реакции.
Таблица 1. Исходные данные интегральной кинетической зависимости концентрации реагирующего компонента А от времени реакции (рис. 1), при рабочей температуре в реакторе.
i |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
, с |
0 |
4 |
8 |
12 |
16 |
20 |
24 |
28 |
33 |
Са, моль/м3 |
30 |
22 |
15,6 |
11,7 |
9,3 |
7,3 |
5,8 |
4,8 |
3,9 |
i |
10 |
11 |
12 |
13 |
14 |
15 |
16 |
17 |
18 |
, с |
38 |
44 |
50 |
60 |
72 |
88 |
113 |
144 |
187 |
Са, моль/м3 |
3,3 |
2,9 |
2,6 |
2,3 |
2,0 |
1,7 |
1,3 |
0,9 |
0,5 |
Таблица 2. Результаты расчетов дифференциальной кинетической зависимости скорости реакции от концентрации при рабочей температуре в реакторе.
Саc, моль/м3 |
26 |
18,8 |
13,65 |
10,5 |
8,3 |
6,55 |
5,3 |
4,35 |
3,6 |
υr, моль/м3с |
2 |
1,6 |
0,975 |
0,6 |
0,5 |
0,375 |
0,25 |
0,18 |
0,12 |
Саc, моль/м3 |
3,1 |
2,75 |
2,45 |
2,15 |
1,85 |
1,5 |
1,1 |
0,7 |
υr, моль/м3с |
0,0677 |
0,05 |
0,03 |
0,025 |
0,01875 |
0,016 |
0,0129 |
0,0093 |
Рис.
2. Дифференциальная кинетическая кривая
(1) зависимости скорости реакции от
концентрации реагирующего компонента
А
и рабочая линия (2) реактора идеального
смешения.
Таблица 3. Исходные и справочные данные и расчетные параметры программы «ИPИC» [4-6].
№ |
Наименование параметра |
Размер- ность |
Обозначение |
Вели-чина |
|
в лекциях |
в программе |
||||
1 |
2 |
3 |
4 |
5 |
6 |
|
Исходные данные |
|
|
|
|
1 |
Производительность по реакционной массе |
м3/час |
qυ |
qv |
1 |
2 |
Начальная концентрация основного реагирующего компонента А |
моль А/м3 |
САО |
са0 |
30 |
3 |
Степень конверсии |
- |
АК |
хак |
0,92 |
4 |
Массив концентраций компонента А в интегральной кинетической зависимости |
мольА/м3 |
СА |
сa(i) |
см. табл. 1 |
5 |
Массив времени реакции в интегральной кинетической зависимости |
с |
τ |
t(i) |
см. табл. 1 |
6 |
Число точек в вышеназванных массивах |
- |
n |
n |
18 |
7 |
Температура реакционной массы на входе в реактор |
С |
tН |
tн |
25 |
|
Справочные данные |
|
|
|
|
1 |
Рабочая температура реакционной массы в реакторе и на выходе из него |
С |
tK |
tk |
35 |
2 |
Тепловой эффект реакции по компоненту А |
кДж/моль А |
qt |
qt |
2755 |
3 |
Плотность реакционной массы при рабочей температуре |
кг/м3 |
ρ |
ro |
968 |
4 |
Вязкость реакционной массы при 0 С |
Пас |
μ0 |
υi0 |
0,0005 |
5 |
Температурный коэффициент вязкости реакционной массы |
К-1 |
β |
be |
0,062 |
6 |
Теплоемкость реакционной массы |
кДж/кгК |
сР |
ср |
3,920 |
7 |
Теплопроводность реакционной массы |
Вт/мК |
λ |
la |
0,62 |
8 |
Начальная температура хладагента (воды) на входе в рубашку |
С |
tWH |
twн |
10 |
9 |
Конечная температура хладагента (воды) на выходе из рубашки |
С |
tWK |
twk |
25 |
1 |
2 |
3 |
4 |
5 |
6 |
10 |
Теплоемкость хладагента (воды) |
кДж/кгК |
сW |
cw |
4,190 |
11 |
Вязкость хладагента (воды) при 0 С |
Пас |
μW |
υiw0 |
0,001 |
12 |
Коэффициент вязкости хладагента (воды) |
К-1 |
βW |
bew |
0,0057 |
13 |
Плотность хладагента (воды) при средней температуре в рубашке |
кг/м3 |
ρW |
row |
980 |
14 |
Теплопроводность хладагента (воды) при средней температуре |
Вт/мК |
λW |
law |
0,68 |
15 |
Ширина кольцевого сечения рубашки, в котором движется хладагент |
м |
δР |
dp |
0,006 |
16 |
Толщина стальной стенки реактора |
м |
δС |
dc |
0,006 |
17 |
Теплопроводность стальной стенки реактора |
Вт/мК |
λС |
lac |
46,2 |
18 |
Термическое сопротивление стенки рубашки (ржавчина, накипь, солевой камень) |
|
rC |
rc |
0,00046 |
19 |
Относительный диаметр мешалки в реакторе идеального смешения |
- |
dОТ |
dot |
0,3 |
20 |
Число оборотов промышленной мешалки |
об/с |
nM |
nm |
7 |
|
Расчетные параметры |
|
|
|
|
1 |
Конечная концентрация реагирующего компонента А |
моль/м3 |
САК |
сак |
2,4 |
2 |
Массив средних концентраций компонента А для дифференциальной кинетической зависимости |
моль/м3 |
САС |
сaс(i) |
см. табл.2 |
3 |
Массив средних скоростей реакции компонента А для дифференциальной кинетической зависимости |
моль/м3с |
υr |
υr(i) |
см. табл.2 |
4 |
Скорость химической реакции в реакторе идеального смешения, соответствующая конечной концентрации САК |
моль/м3с |
υrк |
υrk |
0,0267 |
5 |
Среднее время пребывания в реакторе идеального смешения |
с |
τCM |
tcm |
1035 |
6 |
Среднее время пребывания в реакторе идеального вытеснения |
с |
τВ |
tв |
56,67 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
Отношение объемов реактора идеального смешения и вытеснения |
- |
VOT |
υot |
18,26 |
8 |
Объем реакционной массы в реакторе идеального смешения |
м3 |
VP |
υP |
0,2875 |
9 |
Тепловая мощность реактора |
кВт |
Q |
q |
10,58 |
10 |
Массовый расход хладагента в рубашке |
кг/час |
GW |
gw |
605,9 |
11 |
Расчетный диаметр реактора |
м |
DА |
da |
0,66 |
12 |
Высота цилиндрической части реактора |
м |
HA |
ha |
0,66 |
13 |
Эквивалентный диаметр кольцевого зазора в рубашке |
м |
dЭ |
dе |
0,012 |
14 |
Диаметр пропеллерной мешалки |
м |
DM |
dm |
0,198 |
15 |
Площадь сечения кольцевого зазора рубашки |
м2 |
SM |
sm |
0,0171 |
16 |
Вязкость реакционной массы при температуре в реакторе |
Пас |
μ |
υi |
0,0004 |
17 |
Вязкость хладагента при средней температуре в рубашке |
Пас |
μW |
υiw |
0,0009 |
18 |
Число Рейнольдса для реакционной массы в реакторе с мешалкой |
- |
Re |
Re |
6,61105 |
19 |
Число Прандтля реакционной массы |
- |
Pr |
Pr |
25,45 |
20 |
Объемный расход хладагента |
м3/с |
QW |
qw |
1,7210-4 |
21 |
Скорость хладагента в рубашке |
м/с |
υX |
υx |
0,0134 |
22 |
Число Рейнольдса для хладагента в рубашке |
- |
ReX |
Rex |
175,8 |
23 |
Число Прандтля для хладагента в рубашке |
- |
PrX |
Prx |
9,54 |
24 |
Средняя движущая сила по температуре |
К |
ΔtC |
dtc |
16,37 |
25 |
Средняя температура хладагента |
С |
tXC |
txc |
18,6 |
26 |
Отношение чисел Прандтля для реакционной массы |
- |
|
Pr0 |
0,957 |
27 |
Число Нуссельта для реакционной массы |
- |
Nu |
Nu |
2477 |
28 |
Коэффициент теплоотдачи для реакционной массы |
Вт/м2К |
α |
al |
775,4 |
1 |
2 |
3 |
4 |
5 |
6 |
29 |
Отношение чисел Прандтля для хладагента |
- |
Pr0X |
Pr0x |
1,035 |
30 |
Число Нуссельта для хладагента |
- |
NuX |
Nux |
1554 |
31 |
Коэффициент теплоотдачи для хладагента |
Вт/м2К |
αХ |
alx |
880,1 |
32 |
Удельная тепловая мощность горячей реакционной массы |
Вт/м2 |
qГ |
qg |
5493 |
33 |
Удельная тепловая мощность хладагента |
Вт/м2 |
qX |
qx |
5426 |
34 |
Средняя удельная тепловая мощность |
Вт/м2 |
qC |
qc |
5458 |
35 |
Температура реакционной массы на стенке реактора |
С |
Х |
х |
28 |
36 |
Температура хладагента на стенке рубашки |
С |
tW |
tw |
24,80 |
37 |
Поверхность боковой цилиндрической стенки реактора |
м2 |
FC |
fc |
1,37 |
38 |
Поверхность эллиптического днища реактора |
м2 |
FД |
fd |
0,59 |
39 |
Общая поверхность стенки и днища реактора |
м2 |
FA |
fa |
1,96 |
40 |
Необходимая расчетная поверхность теплопередачи |
м2 |
F |
f |
1,95 |
41 |
Высота рубашки в реакторе |
м |
HP |
hp |
0,66 |
42 |
Коэффициент теплопередачи |
Вт/м2К |
Kt |
kt |
331,6 |
б) с рубашкой на
боковой стенке
(Hp<Ha)
Вид I
Рис. 4 Схема
изменения температуры реакционной
массы и средней температуры хладагента
поперек стенки рубашки.
1 – стальная
стенка;
2 – ржавчина со
стороны хладагента;
3 – накипь со
стороны реакционной массы;
4 – тепловые
пограничные слои.
Рис. 3 Схема реактора
а) с рубашкой на
боковой стенке и днище (максимально
возможная поверхность теплопередачи
Fc
+ Fд)
