- •Расчет на эвм кожухотрубного экзотермического реактора идеального вытеснения
- •Введение.
- •Цель работы.
- •I. Постановка задачи
- •II. Алгоритм расчета
- •III. Оформление отчета по лабораторной работе.
- •IV. Вопросы и задания для самостоятельной работы.
- •Список рекомендуемой литературы.
- •Расчет на эвм кожухотрубного экзотермического реактора идеального вытеснения
- •400131 Волгоград, пр. Ленина, 28
- •400131 Волгоград, ул. Советская, 35
ФЕДЕРАЛЬНОЕ АГЕНТСТВО ПО ОБРАЗОВАНИЮ
ВОЛГОГРАДСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ
КАФЕДРА «ПРОМЫШЛЕННАЯ ЭКОЛОГИЯ И
БЕЗОПАСНОСТЬ ЖИЗНЕДЕЯТЕЛЬНОСТИ»
Расчет на эвм кожухотрубного экзотермического реактора идеального вытеснения
Методические указания к лабораторной работе по курсу
«Математическое моделирование стационарных систем
в химической технологии»
РПК
«Политехник»
Волгоград
2006
УДК 536.629
Рецензент Г.В. Рябчук
Печатается по решению редакционно-издательского совета
Волгоградского государственного технического университета
Расчет на ЭВМ кожухотрубного экзотермического реактора идеального вытеснения. Методические указания к лабораторной работе. Сост. А.Б. Голованчиков, Н.А. Дулькина, В.А. Козловцев, А.А. Шагарова. Волгоград. гос. техн. ун-т. – Волгоград, 2006. – 20 с.
Приведен алгоритм расчета кожухотрубного экзотермического реактора идеального вытеснения. Дана таблица идентификаторов с исходными и справочными данными и расчетными параметрами.
Предназначены для студентов дневной и вечерней форм обучения по направлениям 550800 «Химическая технология и биотехнология» и 551800 «Технологические машины и оборудование»
Табл.4, Библиогр. 4 назв.
© Волгоградский государственный
технический университет, 2006
Введение.
Кожухотрубчатые экзотермические реакторы идеального вытеснения находят широкое применение в химической и нефтехимической промышленности.
Основной проблемой при их расчете и проектировании является определение диаметра и числа трубок в трубном пучке, обеспечивающих необходимый объем реакционной массы в трубках и поверхность теплопередачи [1, 2].
Обычно коэффициент теплопередачи выбирается по справочной и учебной литературе. В настоящей лабораторной работе он варьируется до тех пор, пока его выбранное значение не становится равным расчетному.
Цель работы.
Изучить алгоритм расчета кожухотрубного экзотермического реактора идеального вытеснения.
Научить студентов выбирать справочные данные по учебной и научной литературе.
Научить студентов работать на ЭВМ в диалоговом режиме, варьируя задаваемый коэффициент теплопередачи так, чтобы он соответствовал его расчетному значению.
Уметь иллюстрировать проводимые расчеты в виде схем и графиков.
Делать выводы по результатам расчетов на ЭВМ.
I. Постановка задачи
Рассчитать технологические параметры и геометрические размеры кожухотрубного экзотермического реактора идеального вытеснения по заданной производительности, степени конверсии, начальной концентрации основного реагирующего компонента и данной для него зависимости скорости реакции от концентрации (дифференциальной кинетической кривой рис.1), а также тепловому эффекту реакции.
Рис.1. Дифференциальная кинетическая кривая химической реакции. |
Таблица 1. Зависимость скорости реакции от концентрации для дифференциальной кинетической кривой (рис.1).
Са |
72 |
70 |
68 |
66 |
64 |
62 |
60 |
58 |
56 |
υr |
0,53 |
0,521 |
0,516 |
0,51 |
0,5 |
0,495 |
0,482 |
0,476 |
0,46 |
Са |
54 |
52 |
50 |
48 |
46 |
44 |
42 |
40 |
38 |
υr |
0,42 |
0,4 |
0,38 |
0,361 |
0,33 |
0,30 |
0,26 |
0,23 |
0,21 |
Са |
36 |
34 |
32 |
30 |
28 |
26 |
24 |
22 |
20 |
υr |
0,19 |
0,18 |
0,167 |
0,155 |
0,15 |
0,144 |
0,14 |
0,135 |
0,13 |
Са |
18 |
16 |
14 |
12 |
10 |
8 |
6 |
4 |
2 |
0 |
υr |
0,125 |
0,12 |
0,1 |
0,08 |
0,06 |
0,04 |
0,010 |
0,005 |
0,0025 |
0 |
* Примечание: последнее значение скорости реакции υr лучше записать в виде числа больше нуля, например 0,0001.
Таблица 2. Рассчитанная на ЭВМ зависимость концентрации от времени реакции (интегральная кинетическая зависимость).
Са |
71 |
69 |
67 |
65 |
63 |
61 |
59 |
57 |
55 |
τ |
3,77 |
7,61 |
11,49 |
15,41 |
19,41 |
23,45 |
27,6 |
31,8 |
36,15 |
Са |
53 |
51 |
49 |
47 |
45 |
43 |
41 |
39 |
37 |
τ |
40,9 |
45,9 |
51,17 |
56,7 |
62,8 |
69,5 |
77,2 |
85,8 |
95,4 |
Са |
35 |
33 |
31 |
29 |
27 |
25 |
23 |
21 |
19 |
τ |
115,9 |
117 |
129 |
141,9 |
155,2 |
169,1 |
183,4 |
198,2 |
213,6 |
Са |
17 |
15 |
13 |
11 |
9 |
7 |
5 |
3 |
1 |
τ |
229,6 |
246,3 |
266,3 |
291,3 |
324,6 |
374,6 |
574,6 |
974,6 |
1174,6 |
Рис.2. Зависимость концентрации реагирующего компонента А
от времени реакции (интегральная кинетическая кривая,
построенная по данным таблицы 2).
На интегральной кинетической кривой обозначена точка А, соответствующая рассчитанным значениям конечной концентрации Сак= 7,2 и среднему времени пребывания реакционной массы в реакторе τс = 369,7 с. |
Таблица 3. Исходные и справочные данные и расчетные параметры программы «RИV».
№ |
Наименование параметра |
Размер- ность |
Обозначение |
Величина |
|||
в лекциях |
в програм-ме |
||||||
1 |
2 |
3 |
4 |
5 |
6 |
||
|
Исходные данные |
|
|
|
|
||
1 |
Производительность |
м3/час |
qV |
qv |
1 |
||
2 |
Степень конверсии по основному реагирующему компоненту |
- |
А |
xa |
0,90 |
||
3 |
Массив скоростей реакции |
кг/(м3с) |
υr |
vr(i) |
см. табл.1 |
||
4 |
Равномерный интервал разбиения дифференциальной кинетической кривой по концентрации |
кг/м3 |
ΔСA |
dсa |
2 |
||
5 |
Число точек в массиве скоростей реакции |
- |
n |
n |
35 |
||
6 |
Начальная концентрация реагента А |
кг/м3 |
САО |
са0 |
72 |
||
|
Справочные данные |
|
|
|
|
||
1 |
Молекулярная масса компонента А |
кг/кмоль |
МА |
ма |
112 |
||
2 |
Удельная теплоемкость реакционной массы |
кДж/(кгК) |
сР |
ср |
3,25 |
||
3 |
Тепловой эффект экзотермической реакции |
кДж/кмоль |
qt |
qt |
89800 |
||
4 |
Плотность реакционной массы при рабочей температуре в реакторе |
кг/м3 |
ρ |
ro |
780 |
||
5 |
Вязкость реакционной массы при рабочей температуре в реакторе |
Пас |
μ |
υi |
0,00065 |
||
1 |
2 |
3 |
4 |
5 |
6 |
||
6 |
Температура реакционной массы: - на входе в реактор - на выходе из реактора |
С
|
tН tК |
tн tк |
20 30 |
||
7 |
Температура хладагента (воды): - на входе в реактор - на выходе из реактора |
С
|
tХН tХК |
tхн tхк |
5 20 |
||
8 |
Удельная теплоемкость хладагента |
кДж/(кгК) |
cХ |
сх |
4,19 |
||
9 |
Толщина стенки труб трубного пучка |
м |
δ |
dw |
0,002 |
||
10 |
Термическое сопротивление ржавчины, накипи, солевого камня |
|
r |
r |
0,0082 |
||
11 |
Теплопроводность реакционной массы |
Вт/(мК) |
λГ |
lg |
0,58 |
||
12 |
Теплопроводность хладагента (воды) |
- |
λХ |
lx |
0,68 |
||
13 |
Теплопроводность стальной стенки труб трубного пучка |
- |
λС |
lw |
46,2 |
||
14 |
Длина труб трубного пучка |
м |
l |
l |
6 |
||
15 |
Плотность хладагента при его средней температуре |
кг/м3 |
ρХ |
rox |
956 |
||
16 |
Динамическая вязкость хладагента при его средней температуре |
Пас |
μХ |
υix |
0,00078 |
||
17 |
Коэффициент вязкости реакционной массы |
1/K |
аг |
ag |
0,067 |
||
18 |
Коэффициент вязкости хладагента |
1/K |
aX |
ax |
0,058 |
||
19 |
Коэффициент объемного расширения реакционной массы |
1/К |
bг |
bg |
0,00125 |
||
20 |
Задаваемый коэффициент теплопередачи (30÷270) [4,стр 175] |
Вт/(м2К) |
Кt |
kt |
38,85 |
||
|
Расчетные параметры |
|
|
|
|
||
1 |
Конечная концентрация реагирующего компонента А |
кг/м3 |
САК |
сак |
7,2 |
||
2 |
Массив средних концентраций реагирующего компонента А на интегральной кинетической кривой |
кг/м3 |
СА |
сa(i) |
см. табл.2 |
||
3 |
Массив времени на интегральной кинетической кривой |
сек |
τ |
t(i) |
см. табл.2 |
||
4 |
Среднее время пребывания реакционной массы в реакторе |
сек |
τС |
tc |
369,6 |
||
5 |
Тепловая мощность реактора |
кВт |
Q |
q |
7,387 |
||
6 |
Расход хладагента (воды) |
кг/час |
GX |
gx |
423,13 |
||
7 |
Средняя движущая сила теплопередачи от реакционной массы к хладагенту при прямотоке |
К |
ΔtC |
dtc |
12,33 |
||
8 |
Необходимая поверхность теплопередачи при заданном коэффициенте теплопередачи |
м2 |
F |
f |
15,42 |
||
9 |
Объем труб трубного пучка реактора |
м3 |
V |
υ |
0,1026 |
||
10 |
Расчетный внутренний диаметр труб трубного пучка |
м |
dв |
db |
0,0266 |
||
11 |
Теоретическое число труб трубного пучка |
- |
nT |
nt |
30,7 |
||
12 |
Диаметр кожуха (корпуса реактора) |
м |
DК |
dk |
0,453 |
||
13 |
Площадь межтрубного пространства, в котором движется хладагент |
м2 |
SМ |
sm |
0,138 |
||
14 |
Эквивалентный диаметр межтрубного пространства |
м |
DЭМ |
dem |
0,126 |
||
15 |
Число Рейнольдса для хладагента |
- |
ReX |
rex |
137,6 |
||
16 |
Средняя температура хладагента |
C |
tX |
tx |
12,5 |
||
17 |
Средняя температура реакционной массы |
C |
tг |
tg |
24,8 |
||
18 |
Средняя скорость реакционной массы в трубах трубного пучка |
м/с |
υC |
υc |
0,016 |
||
1 |
2 |
3 |
4 |
5 |
6 |
||
19 |
Число Рейнольдса для реакционной массы |
- |
Reг |
reg |
518,8 |
||
20 |
Число Прандтля для реакционной массы в трубах |
- |
Prг |
prg |
3,64 |
||
21 |
Число Прандтля для хладагента в межтрубном пространстве |
- |
PrX |
prx |
4,81 |
||
22 |
Отношение чисел Прандтля реакционной массы и стенки в степени 0,25 |
- |
|
pr0 |
0,959 |
||
23 |
Число Грасгофа для реакционной массы |
- |
Grг |
grg |
8,35*105 |
||
24 |
Число Нуссельта для реакционной массы |
- |
Nuг |
nug |
8,74 |
||
25 |
Коэффициент теплоотдачи для реакционной массы |
Вт/(м2К) |
г |
alg |
190,4 |
||
26 |
Удельная тепловая мощность для реакционной массы |
Вт/м2 |
qг |
qg |
477 |
||
27 |
Разность температур в ржавчине, стенке, накипи и солевом камне (процессов теплопроводности) |
К |
ΔtW |
dtw |
0,02 |
||
28 |
Температура ржавчины, солевого камня или накипи на стенке со стороны реакционной массы |
К |
X |
x |
22,33 |
||
29 |
Температура ржавчины, солевого камня или накипи на внешней стенке трубы со стороны хладагента |
К |
tW |
tw |
22,31 |
||
30 |
Число Нуссельта для хладагента в межтрубном пространстве |
- |
NuX |
nux |
9,14 |
||
31 |
Коэффициент теплоотдачи для хладагента |
Вт/(м2К) |
αX |
alx |
49,1 |
||
32 |
Удельная тепловая мощность для хладагента |
Вт/м2 |
qX |
qx |
481 |
||
33 |
Средняя удельная тепловая мощность |
Вт/м2 |
qC |
qc |
479 |
||
34 |
Разница удельных тепловых мощностей |
- |
U |
u |
-4,42 |
||
35 |
Расчетная поверхность теплопередачи |
м2 |
FP |
fr |
15,41 |
||
36 |
Расчетный коэффициент теплопередачи |
Вт/(м2К) |
Кtp |
KtR |
38,86 |
||
Рис. 3. Схема кожухотрубного реактора с прямоточным движением реакционной массы в трубах трубного пучка и хладагента: 1 – трубки трубного пучка; 2 – трубные решетки; 3 – корпус; 4 – крышки корпуса.
Вид А |
Рис. 4. Схема изменения средней температуры реакционной массы и хладагента поперек стенки трубки: 1 – ржавчина и накипь на внутренней стенке трубки; 2 – стенка трубки толщиной δ; 3 – ржавчина и солевой камень на внешней стенки трубки; 4 – тепловые пограничные слои (ограничены пунктирными линиями).
|
Рис. 5. Схема изменения температуры реакционной массы (1) и хладагента (2) по длине трубок трубного пучка (режим прямотока). |
