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

Учебное пособие 800429

.pdf
Скачиваний:
8
Добавлен:
01.05.2022
Размер:
2.55 Mб
Скачать

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

Таблица 2 Зависимость массы загрузки от времени режима осушки

t , час

 

1

 

 

3

 

 

5

 

x ,

0,0135

0,026

0,04

0,0135

0,026

0,04

0,0135

0,026

0,04

 

 

 

 

 

 

 

 

 

кг/м3

 

 

 

 

 

 

 

 

 

 

q ,

216

216

 

216

 

216

 

 

 

216

 

216

 

216

216

216

м3/сут

 

 

 

 

 

 

 

 

Мад ,

2,43

4,68

 

7,2

 

7,29

 

 

 

14,01

 

21,6

 

12,15

23,4

36

 

кг

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V , м3

0,362

0,698

 

1,075

 

1,08

 

 

 

2,091

 

3,22

 

1,813

3,49

5,37

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t , час

 

10

 

 

 

 

 

 

 

 

 

14

 

 

 

 

16

 

 

x ,

0,00016

 

0,00024

 

0,0004

 

0,00016

 

 

 

0,00024

 

0,0004

 

0,00016

 

0,00024

0,0004

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

кг/м3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q ,

216

 

216

 

216

 

216

 

 

 

216

 

216

 

216

 

216

216

м3/сут

 

 

 

 

 

 

 

 

 

 

Мад ,

0,428

 

0,428

 

0,72

 

0,4

 

 

 

0,599

 

1

 

0,461

 

0,685

1,152

 

кг

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V , м3

0,043

 

0,064

 

0,107

 

0,06

 

 

 

0,089

 

0,149

 

0,0688

 

0,102

0,172

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Проведём гидравлический расчет адсорбера.

 

 

 

 

 

Число Рейнольдса Re, определяют по формуле

 

 

 

 

 

 

 

 

Re

0.45

 

 

 

 

 

W1ср

dэ

,

 

 

 

(2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

где

dэ

средний

диаметр

частицы, dэ

=

0,0036

м [2];

,

– пористость, , = 0,49 [2]; – кинематическая вязкость, м2/с.

11

 

 

 

 

 

2 Твых

 

Твх

,

 

 

 

 

 

 

(3)

 

 

f

 

 

 

 

 

 

Тср

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где Тср

– средняя эффективная температура газа в адсорбере, К;

Твх – температура на входе

 

 

в адсорбер, Твх

=

293 К;

 

Твых

температура на выходе из адсорбера, Твых

 

= 433 К.

 

 

 

Среднюю температуру в адсорбере Тср.ад ., К определяют по

формуле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т

 

 

Твх

Твых

.

 

 

 

 

 

 

 

(4)

 

 

ср

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Среднюю эффективную температуру газа в адсорбере Т , K

вычисляют по формуле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т

 

 

 

 

 

 

Твх Твых

.

 

 

 

 

 

 

 

(5)

 

 

ср.эф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тср

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Среднюю эффективную плотность газа в адсорбере

 

, кг/м3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ср.ад.

 

вычисляют по формуле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

273

 

,

 

 

 

 

 

 

 

(6)

 

 

ср.ад.

 

0

Tср

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

0

– плотность газа на входе адсорбера, кг/м3.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Плотность газа на входе адсорбера

 

 

 

0

, кг/м3

определяют по

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

формуле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

104

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

(7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R Tвх

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Принимая время режима

 

 

 

осушки

 

t

 

=

 

12

часов,

диаметр

адсорбера D = 144 мм.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определяют скорость газа на входе в адсорбер W1 , м/с по

формуле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W1

 

 

 

 

 

 

 

 

m R

Tвх

 

 

 

 

 

.

 

 

 

(8)

 

 

 

 

 

104 0,7854 0,1442

 

 

 

вх

Среднюю эффективную скорость газа в адсорбере W1ср , м/с вычисляют по формуле

12

W

W

Т

ср

W

Твх

Твых

.

(9)

 

 

 

 

 

1ср

1

Т

вх

1

2Т

вх

 

 

Подставляя значения в формулы (2) – (9) получаем следующее:

 

 

 

 

W1

 

0,00165 52,8 300

 

0,082 м/с;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 104 0,7854 0,1442

 

 

 

 

 

W

0,082 293 433 0,101 м/с;

 

 

 

 

 

 

 

1ср

 

 

 

 

 

 

 

2 293

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т

 

 

 

293 433

363

К;

 

 

 

 

 

2 104

1, 29

кг/м3;

ср.ад.

 

 

2

 

 

 

 

0

52,8 293

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ср.ад.

1,29

273

1,041

кг/м3;

f

 

 

2 293 433

0,771 ;

363

363

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Re

 

 

0.45

 

 

 

 

 

 

0,101 0, 0036

33, 25 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

1 0, 49

 

 

 

 

 

14,71 10 6

 

 

0, 49

 

 

 

 

 

 

 

 

 

75

 

 

 

 

15

 

 

1

3,515 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

33,25

 

 

 

 

 

 

 

 

 

 

 

33,25

 

 

 

 

 

30,6 3,515

1, 4614

0,771 436629,8 .

 

 

 

0,0036

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, перепад давления на адсорбере при режиме регенерации р , МПа может составить

 

 

 

 

 

 

 

 

 

 

W 2

 

 

 

 

 

 

 

 

 

ср

1ср

,

(10)

 

 

р

 

 

 

 

 

 

 

 

 

 

 

 

2 g

 

 

 

 

 

 

 

 

 

 

 

 

 

436629,8

1, 041 0,1012

 

236,32 кгс/м2 = 0,0024 МПа.

р

2 9,81

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На

 

режиме осушки

 

эта

 

величина

составляет примерно

0,08 МПа.

Литература

1.Коршак А.А. Основы нефтегазового дела / А.А. Коршак, А.М. Шаммазов. – Уфа: ДизайнПолиграфСервис, 2001. – 249 с.

2.Юрьев А.С. Справочник по расчётам гидравлических и вентиляционных систем / Юрьев А.С., Пирогов С.Ю. и др – С.-Пб: АНО НПО «Мир и семья», 2001. – 1154 с.

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

13

УДК 620.9

Р.А. Константинов – магистрант, В.Ю. Дубанин, канд. техн. наук, доц. – научный руководитель

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ТЕПЛОЭНЕРГЕТИЧЕСКОЙ СИСТЕМЫ, СОСТОЯЩЕЙ ИЗ ПОСЛЕДОВАТЕЛЬНО СОЕДИНЕННОЙ ГИДРОДИНАМИЧЕСКОЙ СТРУКТУРЫ

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

Ключевые слова: математическая модель, теплоэнергетическая система, теплоноситель, гидродинамическая структура

Теплоэнергетические системы, представляющие собой последовательное соединение зон перемешивания и зон вытеснения, получили широкое распространение в химической и нефтехимической промышленности. Ярким примером таких систем является реакторное оборудование химических производств. Как правило, в целях изучения процессов, протекающих в данных системах, используют их идеальные варианты – аппараты идеального вытеснения и идеального перемешивания [1]. Идеализация протекающих процессов позволяет исключить влияние менее значимых факторов, что в значительной степени облегчает процесс моделирования таких систем.

Рассмотренная система может быть представлена в виде схемы, показанной на рис. 1. Задача моделирования состоит в нахождении температуры теплоносителя на выходе из системы Tвых (t) и

определении влияния технологических параметров (объема V аппарата, расхода U теплоносителя, начальной температуры теплоносителя Tвх (t) и т.п.) на конечный результат.

При математическом моделировании системы, приведенной на рис. 1, будем рассматривать каждый ее объект в отдельности друг от друга. При этом изменение температуры в зоне идеального перемешивания представим в виде системы дифференциальных уравнений первого порядка [2].

14

Рис. 1. Теплоэнергетическая система, состоящая из последовательно соединенных зон идеального перемешивания и идеального вытеснения

Для первого аппарата:

 

 

dT1 (t1 )

 

 

1

T (t ) Т (t ) ;

 

 

 

 

 

 

 

dt1

 

 

 

 

вх 1

1

1

 

 

 

 

1

 

 

 

(1)

 

 

 

 

Т1 (0) T0 .

 

 

 

 

 

 

 

 

 

 

 

 

Для второго аппарата:

 

 

 

 

 

 

 

 

 

 

 

dT2 (t2 )

 

 

1

T (t ) Т

 

(t

 

) ;

 

 

 

 

 

2

2

 

dt2

 

1 1

 

 

 

 

 

2

 

 

 

(2)

 

 

 

 

Т2 (0) T1 (t1 ).

 

 

 

 

 

 

 

 

 

 

 

 

где 1 V1 / U,

2

V2 / U – время пребывания теплоносителя в

аппаратах 1 и 2 соответственно, с.

Система уравнений, описывающих изменения температуры в зоне идеального вытеснения, выглядит следующим образом

 

T3

(z,t3 )

u

T3

(z,t3 )

;

 

dt3

 

dz

 

 

 

 

 

 

 

 

 

Т3 (0) T0 ;

 

(3)

 

 

 

 

 

Т3 (0,t) T2 (t2 ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где u U / S3 – скорость движения теплоносителя в аппарате 3,

м/с.

При этом очевидно, что температура на выходе из системы будет равна изменению температуры в третьем аппарате по его длине, то есть

Tвых (t) T3 (l3 ,t),

15

где l3 V3 / S3 – длина аппарата 3.

Для решения систем 1-3 необходимо привести их к безразмерному виду. Для этого введем следующие поправочные

коэффициенты:

 

t / 1

безразмерное

время;

1 / 2

коэффициент,

учитывающий

изменение

времени

во втором

аппарате;

u 1

коэффициент, учитывающий

изменение

l3

скорости по сечению третьего аппарата.

Сучетом введенных безразмерных коэффициентов системы (1)-

(3)примут вид:

 

d 1 ( )

вх ( ) 1 ( ) ;

 

 

 

 

 

 

 

 

 

 

d

(4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 (0) 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d 2 ( )

1 ( ) 2

( ) ;

 

 

 

 

 

 

 

 

 

 

d

(5)

 

 

 

 

 

 

 

 

 

 

 

2 (0) 1 ( ).

 

 

 

 

 

 

 

 

 

3 (Z , )

 

3 (Z , )

;

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 (0) 0;

 

 

 

(6)

 

 

 

 

 

 

 

 

 

 

3

(0, ) 2 ( ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

Z z / (u 3 ) ;

1

(T1 T0 ) / T0 ;

2 (T2 T0 ) / T0 ;

3 (T3 T0 ) / T0 .

По определению передаточной функции [1] имеем:

вых

 

вх

 

 

 

W (Р) L(

) / L(

) ,

(7)

где W (Р) – передаточная

 

функция

 

всей системы;

L

преобразование Лапласа.

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

W (Р) W1 (Р) W2 (Р) W3 (Р),

(8)

где W1 (Р) , W2 (Р) , W3 (Р) – передаточные функции аппаратов 1,

2 и 3 соответственно. При этом имеем:

16

W (Р)

1

 

, W (Р)

 

, W (Р) е( Р/ ) .

 

 

 

 

 

 

1

Р 1

2

Р

3

 

 

 

Подставляя данные значения в (8), получим выражение передаточной функции для всей системы

W (Р)

е( Р/ )

(9)

.

(Р 1)(Р )

Выражая вых ( ) из (7) получим дифференциальное уравнение

для нахождения температуры теплоносителя на выходе из системы (в безразмерном виде)

 

вых

( ) L 1 W (P) L

( ) ,

 

 

вх

 

где L – обратное преобразование Лапласа.

Подставляя (9) в полученное выражение, придем к окончательному виду уравнения для нахождения искомой величины

 

 

 

е( Р/ )

 

 

 

 

вых

( ) L 1

 

 

L вх

( )

,

(10)

 

 

 

 

(Р 1)(Р )

 

 

 

 

Связь Tвх (t) с

вх ( ) и

Tвых (t) с вых ( )

выглядит следующим

образом:

вх ( )

вых ( )

В итоге имеем

Твх (t) Т0 ;

Твых (t) Т0 .

 

е( Р/ )

 

 

Твых (t) L 1

 

L Tвх (t) T0 T0 ,

(11)

 

 

(Р 1)(Р )

 

 

Зная Tвх (t) , T0 , U ,

V и S с помощью выражения (11) можно

легко найти Tвых (t) .

 

 

 

 

Так для системы, изображенной на

рис.1, при T0

= 20 С;

U = 0,015 м3/с; V = 0,1 м3; V = 0,2 м3;

V = 0,3 м3; S

= 0,3 м2 и

1

2

 

3

 

изменении температуры теплоносителя на входе в систему по закону Tвх (t) 40(1 e 0,1t ), используя уравнение (11), получим графическое изображение зависимости изменения температуры теплоносителя от

времени на выходе из системы Tвых (t) (рис. 2). Изменяя входные параметры системы, можно получить наилучшее их соотношение с требуемым выходным значением температуры теплоносителя.

17

Рис. 2. Зависимость температуры теплоносителя от времени на выходе из системы, приведенной на рис. 1

Результатом проведенного в данной работе математического моделирования теплоэнергетической системы с последовательным соединением элементов (рис. 1) стала возможность разработки алгоритма наглядной оптимизации теплоэнергетических параметров объектов системы, позволяющая достичь наилучших техникоэкономических показателей системы.

Литература

1.Бондарь А.Г. Математическое моделирование в химической технологии / А.Г. Бондарь. – Киев: «Вита школа», 1973. - 280 с.

2.Цветков Ф.Ф. Тепломассообмен: учебник для вузов / Ф.Ф. Цветков, Б.А. Григорьев. - М: Издательский дом МЭИ, 2011. - 562 с.

3.Берд Р. Явления переноса / Р. Берд, В. Стьюарт, Е. Лайтфут. - М.: «Химия», 1974. - 688 с.

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

18

УДК 620.9

С.А. Воробьев – магистрант, Ю.Н. Агапов, д-р техн. наук, проф. – научный руководитель

АНАЛИЗ ВЛИЯНИЯ ВЗАИМНОГО РАСПОЛОЖЕНИЯ АППАРАТОВ ИДЕАЛЬНОГО ПЕРЕМЕШИВАНИЯ ПРИ НАЛИЧИИ ЗАСТОЙНОЙ ЗОНЫ

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

Ключевые слова: математическая модель, теплоэнергетическая система, застойная зона, вычислительный эксперимент

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

Рассмотрим данный объект, который изображен на рис. 1.

Рис. 1. Теплоэнергетическая система

19

1. Составим математическую модель системы

1 (t)

 

 

 

1

 

Твх

(t)

Т1 (t) ;

 

 

 

 

 

 

 

 

 

 

 

 

dt

1

 

 

 

 

 

 

 

Т

 

(0) T

 

 

 

 

 

 

 

 

1

 

 

 

0

 

 

 

 

 

 

2 (t)

 

 

 

1

 

Т1

(t) (1 )Т2 (t) T3 (t) ;

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

2

 

 

 

 

3 (t)

 

 

 

 

Т2

 

 

Т3 (t)

 

 

 

 

 

 

 

 

 

 

(t)

 

 

dt

 

 

 

 

 

 

 

 

3

 

 

 

 

Т

2

(0) Т

3

(0) Т

0

.

 

 

 

 

 

 

 

 

 

 

 

2.Приведём систему к безразмерному виду

t1 , 1 Т1 ( ) Т0 , вх ( ) Твх Т0 ,

1 .

2

Запишем систему, применив введённые параметры

d 1 ( )

вх ( ) 1 ( ) ;

 

 

d

 

 

 

 

 

 

 

 

 

 

 

1 (0) 0;

 

d ( )

 

1 ( ) (1 ) 2 ( ) 3 ( ) ;

 

 

 

2

 

 

 

d

 

 

2 (0) 0;

 

 

d 3 ( )

 

 

1 2 ( ) 3 ( ) ;

 

 

 

 

 

 

d

 

3

 

 

 

(0) 0.

 

 

 

 

 

3

 

 

 

 

 

 

(1)

(2)

3. Решение системы:

Для аппарата идеального перемешивания без застойной зоны уравнение передаточной функции запишем в виде:

w1 (р)

 

,

(3)

 

р

 

 

 

а передаточная функция для аппарата идеального перемешивания с застойной зоной имеет вид:

20