Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
РИВ экзо.doc
Скачиваний:
0
Добавлен:
01.07.2025
Размер:
875.01 Кб
Скачать

ФЕДЕРАЛЬНОЕ АГЕНТСТВО ПО ОБРАЗОВАНИЮ

ВОЛГОГРАДСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ

КАФЕДРА «ПРОМЫШЛЕННАЯ ЭКОЛОГИЯ И

БЕЗОПАСНОСТЬ ЖИЗНЕДЕЯТЕЛЬНОСТИ»

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

Методические указания к лабораторной работе по курсу

«Математическое моделирование стационарных систем

в химической технологии»

РПК

«Политехник»

Волгоград

2006

УДК 536.629

Рецензент Г.В. Рябчук

Печатается по решению редакционно-издательского совета

Волгоградского государственного технического университета

Расчет на ЭВМ кожухотрубного экзотермического реактора идеального вытеснения. Методические указания к лабораторной работе. Сост. А.Б. Голованчиков, Н.А. Дулькина, В.А. Козловцев, А.А. Шагарова. Волгоград. гос. техн. ун-т. – Волгоград, 2006. – 20 с.

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

Предназначены для студентов дневной и вечерней форм обучения по направлениям 550800 «Химическая технология и биотехнология» и 551800 «Технологические машины и оборудование»

Табл.4, Библиогр. 4 назв.

© Волгоградский государственный

технический университет, 2006

Введение.

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

Основной проблемой при их расчете и проектировании является определение диаметра и числа трубок в трубном пучке, обеспечивающих необходимый объем реакционной массы в трубках и поверхность теплопередачи [1, 2].

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

Цель работы.

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

  2. Научить студентов выбирать справочные данные по учебной и научной литературе.

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

  4. Уметь иллюстрировать проводимые расчеты в виде схем и графиков.

  5. Делать выводы по результатам расчетов на ЭВМ.

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К

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) по длине трубок трубного пучка (режим прямотока).

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]