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

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

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

 

 

p

2

 

 

 

 

 

 

 

w2 3

(р)

 

 

3

 

 

 

 

 

,

(4)

 

 

 

 

 

2

 

 

 

 

p 1 p

 

2

 

 

2

 

 

 

 

 

 

3

 

 

3

 

 

Передаточная функция всей системы находится по формуле:

 

 

 

 

w p w

p w

p

выхL p

;

 

 

(5)

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

1

 

 

2 3

 

 

 

 

вхL

 

 

 

 

отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

выхL

p w p вхL

p ;

 

 

 

(6)

 

 

 

 

 

вых ( ) L 1 w(р) L вх ( ) .

 

 

 

(7)

По условию задачи:

U = 0,01 м3/с,

V V = 1 м3,

V

= 0,2 м3,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

3

 

T0 = 20 °C,

= 0,4.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Находим недостающие параметры ранее записанных уравнений

 

 

 

 

 

 

 

 

 

 

 

 

 

V

;

 

 

 

 

 

 

 

 

 

 

(8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

100 сек,

 

 

0, 2

 

20 сек.

 

 

 

 

 

1

2

 

 

3

 

 

 

 

 

 

 

 

0,01

 

 

 

 

 

 

 

0,01

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.

 

 

 

 

 

 

 

 

 

 

 

Найдем функцию температуры на входе

 

 

 

 

 

 

 

 

 

 

 

 

вх Tвх (t) T0 ;

 

 

 

 

 

 

(9)

 

 

40 1 exp( 0,1 t

20 40 1 exp( 0,1 t /

1

 

вх

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

20 20

 

2 1 exp( 10

1 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставим функцию температуры на входе в выражение

функции температуры на выходе:

 

 

 

 

 

 

 

 

 

 

 

 

 

вых ( ) L 1 w (р) L 20 2 1 exp( 10 1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С помощью программы Maple, получим график изменения температуры среды на выходе из системы во времени, изображённый на рис. 2.

21

Рис. 2. График зависимости температуры на выходе из системы во времени

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

Литература

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

2.Цветков Ф.Ф. Тепломассообмен / Ф. Ф. Цветков, Б.Л. Григорьев. –

М.: Изд. МЭИ, 2011. – 562 с.

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

22

УДК 620.9

В.И. Даншина – магистрант, А.В. Бараков, д-р техн. наук, проф. – научный руководитель

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

ИПЕРЕМЕШИВАНИЯ

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

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

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

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

Рассмотрим данную систему, которая изображена на рис. 1.

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

23

Составим математическую модель системы, рассматривая каждый объект в отдельности друг от друга:

 

1 (t)

 

1

 

Твх (t) Т1 (t) ;

 

 

 

dt

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т1 (0) Т0 ;

 

 

 

 

 

 

2 (t)

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Твх (t) Т1 (t) ;

 

 

 

 

 

 

 

2

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

(0) Т0 ;

 

 

 

 

 

 

Т2

 

 

 

 

(1)

 

Т (z

, t)

u3

Т (z , t)

;

 

 

3

 

3

 

 

 

 

3

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

z3

 

 

 

(z3 , t) Т0 ,

Т3 (0, t) Твх (t);

 

Т3

 

Т

 

 

Т

 

(l

, t).

 

 

 

 

вых

 

 

 

 

 

3

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где 1 V1

U1 , 2

V2

U2

 

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

аппаратах 1 и 2 соответственно, с; u US3 – скорость движения

теплоносителя в аппарате 3, м/с.

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

t , t l3

, (z, ) Т3 (z, ) Т

0 , ( ) Твх Т0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

u3

 

 

3

 

 

 

Т0

 

 

вх

Т0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

(1 ) U

,

u3 1

 

 

 

U V1 S3

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V2

 

 

 

l

S

 

(1 ) U V

 

 

 

 

 

 

 

2

 

 

 

 

 

 

3

 

 

3

 

 

3

 

 

 

 

 

 

 

 

 

 

 

U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

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

втором

аппарате;

 

 

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

учитывающий

изменение

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

С учетом введенных безразмерных переменных система 1 примет вид:

24

d 1

( )

вх

( ) 1 ( );

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

(0) 0;

 

 

 

 

 

d ( )

вх ( ) 2 ( ) ;

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

(0) 0;

 

 

 

(2)

2

 

 

 

 

3 (z, )

 

 

3

(z, )

 

 

 

;

 

 

 

z

 

 

 

 

 

 

 

 

(z, ) 0,

3 (0, ) вх ( );

3

вых

3 (l, ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Идеальная функция для аппаратов перемешивания имеет вид:

w1 (р)

1

 

, w

 

(р)

 

,

 

 

2

 

р 1

р

 

 

 

 

 

а для функции идеального вытеснения

w3 (р) exp( р) ,

тогда

w(р) (

1

 

 

 

) exp(

р ).

р 1

р

 

 

 

 

Подставим полученные выражения в уравнение температуры на выходе:

 

 

вых ( ) L 1 w(р) L вх ( ) .

 

 

 

По условию задачи: U = 0,01

м3/с,

V

= 0,5

м3, V

= 0,3 м3,

 

 

 

 

 

 

 

1

 

2

V

= 0,2 м3, S = 0,2 м2,

= 0,2.

 

 

 

 

 

 

 

3

3

 

 

 

 

 

 

 

 

 

 

 

Находим недостающие параметры:

 

 

 

 

 

 

1

0,5

 

62,5; 2

 

0,3

 

150;

 

 

 

 

 

 

 

 

 

 

 

0, 2

 

 

 

 

 

0,8 0,01

 

0,01

.

 

l3 1; 0, 42; 3,125.

Найдем функцию температуры на входе:

25

вх

 

 

 

 

 

40 1 exp( 0,1 t 20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 1

exp

 

0,1

t

1

 

 

 

 

 

 

 

 

 

 

 

20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

exp

 

0,1

 

 

t

 

 

 

1 2

1 exp

 

 

 

 

 

1,

 

1

 

t

 

 

0,1 t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V3

 

 

 

 

20 сек.

 

 

 

 

 

 

 

 

 

где t l3

V3 : U

0, 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u3

 

 

S3

 

S3

 

U

 

 

0,01

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставив полученное значение в выражение, получим:

вх

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

1 exp

 

2

1.

 

Подставим функцию температуры на входе в выражение

функции температуры на выходе:

 

 

 

 

 

 

1

w (р)

 

2

1

 

 

 

 

вых ( ) L

 

 

exp 2 1 .

 

 

 

 

 

 

 

 

 

 

 

На рис. 2 изображен график изменения температуры теплоносителя на выходе из системы, с помощью которого можно проанализировать математическую модель.

Рис. 2. График зависимости температуры на выходе от времени

26

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

Литература

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

2.Цветков Ф.Ф. Тепломассообмен / Ф. Ф. Цветков, Б.Л. Григорьев. –

М.: Изд. МЭИ, 2011. – 562 с.

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

УДК 621.181

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

СРАВНЕНИЕ ПРОИЗВОДИТЕЛЬНОСТИ ПАРОВОГО КОТЛА, РАБОТАЮЩЕГО НА ГАЗЕ ИЛИ МАЗУТЕ

В данной работе рассматривается котельная с котлоагрегатами малой производительности ДЕ-10-14ГМ, топливом которых служит топочный мазут марки 100. В настоящее время есть возможность перехода котельной на газообразное топливо. Это возможно вследствие того, что топочная камера котлов серии ДЕ предназначена как для сжигания мазута, так и для сжигания газообразного топлива

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

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

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

Производственные и отопительные котельные должные обеспечивать бесперебойное и качественное теплоснабжение предприятий и жилищно-коммунального сектора. Повышение

27

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

При разработке новой конструкции газомазутных парогенераторов серии ДЕ особое внимание было обращено на увеличение степени заводской готовности парогенераторов в условиях крупносерийного производства, снижения металлоемкости конструкции парогенераторов, приближение эксплуатационных показателей парогенераторов к расчетным. Газомазутный вертикально-водотрубный паровой котел с естественной циркуляцией ДЕ-10-14ГМ производительностью 10 тонн пара в час предназначен для выработки насыщенного или слабоперегретого пара, в системах отопления, вентиляции и горячего водоснабжения.

В паровых котлах ДЕ-10-14ГМ, предназначены для сжигания газа и мазута, топочная камера располагается справа от поверхности нагрева конвективных пучков (при виде с фронта) и имеет глубину, равную длине котла. Основными составными частями котла являются верхний и нижний барабаны, конвективный пучок, фронтовой и боковой экраны, образующие топочную камеру.

Сравнение производительности производится с помощью расчета. Исходными данными для данного расчета являются основные характеристики котлоагрегата ДЕ-10-14ГМ, а также составы топлив.

Основные характеристики котла

Наименование

Единица

Значение

п/п

измерения

газ

мазут

 

 

 

 

 

 

1

2

3

4

5

 

 

 

 

 

1

Паропроизводительность

т/час

10,0

10,0

 

 

 

 

 

2

Давление пара на выходе

кгс/см2

14

14

 

из котла

(ата)

 

 

 

3

Температура:

 

 

 

 

- насыщенного пара

°С

194

194

 

- питательной воды

°С

100

100

 

 

 

 

 

4

Тип горелок (число

-

ГМ-7(1)

ГМ-7(1)

 

горелок, шт.)

 

 

 

 

5

Номинальная тепловая

МВт

5,24(4.5)

5,24(4.5)

 

мощность горелки

(Гкал/час)

 

 

 

 

 

 

 

 

28

Продолжение таблицы

1

 

 

 

2

 

 

 

 

3

 

4

 

5

6

 

Объем топочной камеры

 

 

м3

 

17,4

 

17,4

7

 

Площадь поверхности:

 

 

м2

 

 

 

 

 

 

 

- нагрева радиацией

 

 

 

м2

 

39,02

 

39,02

 

 

 

- нагрева конвекцией

 

 

 

м2

 

110,0

 

110,0

 

 

 

- водяного экономайзера

 

 

-

 

236

 

236

8

 

Марка водяного

 

 

 

 

м3/час;

 

ЭП2-236

 

ЭП2-236

 

 

 

экономайзера

 

 

 

 

кг/час

 

 

 

 

9

 

Расход топлива

 

 

 

 

ккал/м3∙час

 

718

 

673

10

 

Теплоснабжение объема

 

 

ккал/кг∙час

 

361,6∙103

 

361,6∙103

 

 

 

топки

 

 

 

 

 

 

 

 

 

 

 

11

 

Температура газов на

 

 

 

°C

 

1114

 

1069

 

 

 

выходе из топки

 

 

 

 

 

 

 

 

 

12

 

Температура газа за

 

 

 

°C

 

273

 

310

 

 

 

котлом

 

 

 

 

 

 

 

 

 

 

 

13

 

Температура уходящих

 

 

°C

 

146

 

174

 

 

 

газов

 

 

 

 

 

 

 

 

 

 

 

14

 

Расчетный КПД

 

 

 

 

%

 

92,1

 

90,99

 

 

 

 

 

 

 

 

 

 

15

 

Газовое сопротивление

 

 

кгс/м2

 

195,7

 

195,7

 

 

 

котла

 

 

 

 

 

 

 

 

 

 

 

16

 

Расположение труб

 

 

 

-

 

коридорн.

коридорн.

 

 

 

 

 

 

 

 

 

 

 

17

 

Диаметр и толщина

 

 

 

мм

 

51 2,5

 

51 2,5

 

 

 

стенки труб экранов

 

 

 

 

 

 

 

 

18

 

Шаг труб боковых

 

 

 

мм

 

55

 

55

 

 

 

экранов

 

 

 

 

 

 

 

 

 

 

19

 

Поперечный шаг труб

 

 

 

мм

 

110

 

110

 

 

 

 

 

 

 

 

 

 

 

20

 

Продольный шаг труб

 

 

 

мм

 

110

 

110

 

 

 

 

 

 

 

 

 

 

 

Состав

топлива

мазут малосернистый: W p 3 %,

QP 0,3 %,

QP

40211 кДж/кг,

S p

 

0,3 %, Cp

84,65 %,

 

H P 11,7 %,

H

 

 

 

 

op+k

 

 

 

 

 

 

 

Ap 0,05 %.

 

 

 

 

 

 

 

 

 

 

 

 

 

Состав

 

топлива

газ

Минск-Москва:

CH4 92,8 % ,

C2 H6 2,8 % ,

C3H6 0,9 % ,

C4 H10 0,4 % ,

C5 H12 0,1 % ,

N

2

2,5 % ,

CO 0,5 % ,

Qр 36550 кДж/кг.

 

 

 

 

 

 

2

 

н

 

 

 

 

 

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

29

КПД котельной установки , % определяем по формуле

100 (q2 q3 q4 q5 q6 ) ,

где q2 – потери теплоты с уходящими газами, %; q3 – потеря теплоты от химической неполноты сгорания, %; q4 – потеря теплоты от механической неполноты сгорания, %; q5 – потеря теплоты от наружного охлаждения, %; q6 – потеря в виде

физической теплоты шлаков, %.

Были получены следующие данные: КПД котла, работающего на мазуте 90,5 % ; КПД котла, работающего на газе 91,9 % .

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

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

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

Литература

1.Тепловой расчет котельных агрегатов (Нормативный метод) / Под редакцией Н.В. Кузнецова, В.В. Митора. - М.: Энергия, 1973. - 295с.

2.Гусев Ю.Л. Основы проектирования котельных установок / Ю.Л. Гусев. - М.: Стройиздат, 1967. - 291с.

3.Родданик С.Л. Справочник по котельным малой производительности / С.Л. Родданик, Я.Б. околовский. - М.: Энергия, 1975. - 368с.

4.Роддатис К.Ф. Котельные установки: учеб. пособие / Роддатис К.Ф.

-М.: Энергия, 1977. - 432с.

5.Краснощеков Е.А. Задачник по теплопередаче / Краснощеков Е.А., Сукомел А.С. - М.: Энергия, 1980. - 232с.

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

30