Учебное пособие 800429
.pdf
|
|
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
Составим математическую модель системы, рассматривая каждый объект в отдельности друг от друга:
|
dТ1 (t) |
|
1 |
|
Твх (t) Т1 (t) ; |
|||||||||||||
|
|
|
dt |
|
|
|
1 |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
Т1 (0) Т0 ; |
|
|
|
|
|
||||||||||||
|
dТ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