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

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

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

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

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

Расчет на эвм экзотермического реактора идеального смешения

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

РПК

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

Волгоград

2006

УДК 536.629

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

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

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

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

Ил. 4. Табл.3. Библиогр.: 7 назв.

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

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

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

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

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

  1. Введение

Экзотермические реакторы идеального смешения широко применяются в химической и нефтехимической промышленности.

Основной проблемой при их расчете является определение поверхности рубашки и при необходимости поверхности змеевика, обеспечивающих отвод тепла реакции от реакционной массы к хладагенту [1-3].

  1. Цель работы

  1. Изучить алгоритм расчета экзотермического реактора идеального смешения.

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

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

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

  1. Постановка задачи

Рассчитать технологические параметры и геометрические размеры экзотермического реактора идеального смешения по заданной производительности, степени конверсии основного реагирующего компонента, тепловому эффекту реакции и зависимости концентрации основного реагирующего компонента от времени (интегральной кинетической кривой, рис.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Н

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

Среднее время пребывания в реакторе идеального вытеснения

с

τВ

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Э

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,61105

19

Число Прандтля реакционной массы

-

Pr

Pr

25,45

20

Объемный расход хладагента

м3

QW

qw

1,7210-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д)

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