CC BY
176
УДК 621.1.016.7
Л. В. Галимова, А. А. Камнев, О. О. Лазаренко, К. Гуиди Тоньон
МОДЕЛИРОВАНИЕ И ЭКСЕРГЕТИЧЕСКИЙ АНАЛИЗ ОДНОСТУПЕНЧАТОЙ АММИАЧНОЙ ЭКСПЕРИМЕНТАЛЬНОЙ ХОЛОДИЛЬНОЙ МАШИНЫ
Одним из направлений решения проблем энергосбережения является анализ работы технических систем, потребляющих и преобразующих энергию.
Эксергетический анализ любого исследуемого объекта дает после соответствующей обработки информацию не поверхностную, а показывающую скрытые в энергетическом балансе характеристики, необходимые для исследования и совершенствования объекта.
Задача эксергетического анализа - оценка степени термодинамического совершенства технической системы в целом, а также выявление тех стадий технического процесса, на которых сосредоточены основные потери эксергии (потенциальной работоспособности системы).
Оценка энергоэффективности агрегатов проектируемой, эксплуатируемой и модернизируемой холодильной установки осуществляется на основе определения обобщенного показателя энергетической эффективности (эксергетического КПД), значение которого показывает степень приближения установки к идеальной [1].
По существующим оценкам, при выборе процесса определяются 40 % потерь, при проектировании - еще 40 %, и, таким образом, примерно на 80 % потерь на этапе производства воздействовать уже невозможно. Поэтому эксергетический анализ имеет особенно большое значение на стадиях эскизной проработки и проектирования. Для стадии эксплуатации часто достаточны качественные эксергетические соображения [2].
В дальнейшем эксергетический анализ проводится для существующей экспериментальной холодильной машины.
При рассмотрении возможностей уменьшения потерь будут обсуждаться и результаты, которые можно получить путем изменения основного процесса или принципа действия системы.
Целью исследований являлся анализ цикла аммиачной экспериментальной холодильной машины в рабочих режимах.
Для достижения поставленной цели были решены следующие задачи:
- сняты параметры действующей холодильной машины;
- разработана объединенная методика теплового расчета, определения эксергетических потерь в элементах холодильной машины и в соответствии с ней произведены расчеты;
- разработаны блок-схема и программный код теплового расчета, определения эксерге-тических потерь в элементах машины и цикле в целом;
- проведен численный эксперимент и установлена адекватность программного кода;
- построены диаграммы потоков и потерь эксергии по элементам схемы [3];
- построены зависимости эксергетических потерь в элементах схемы от рабочих параметров;
- проанализированы полученные зависимости и даны качественные эксергетические соображения.
Методы исследования
Схема исследуемой экспериментальной холодильной установки представлена на рис. 1.
Лабораторная установка предназначена для обеспечения холодом систем непосредственного и рассольного охлаждения, обслуживающих ряд потребителей (охлаждаемая камера, льдогенератор, батарея рассольного охлаждения). Для общности выводов приведен анализ цикла холодильной машины.
Рабочими параметрами являются: давление конденсации, давление кипения, температура конца сжатия, температура перегрева холодильного агента, температура охлаждающей воды на входе и выходе из конденсатора, ее расход, температура окружающей среды.
Для теплового расчета по методике использовались параметры характерных точек цикла, взятые из программы Соо1Раск. Расчет проведен в диаграмме lgp-h.
Определение массового расхода холодильного агента вели по тепловому балансу конденсатора.
Рис. 1. Схема экспериментальной установки: 1 - компрессор; 2 - отделитель масла; 3 - конденсатор; 4 - регулирующий вентиль; 5 - испаритель; 6 - ротаметр;
7 - счётчик. В - подача воды; К - слив воды; М - выпуск масла из аппарата. Линия трактов: хладагента; — - — ■ воды.
Состояние потока: ■>
Измерение параметра: ^
пар;
температуры;
жидкость.
—■- давления
На рис. 2 представлены методика, результаты теплового расчета, определения эксергети-ческих потерь и КПД элементов и цикла холодильной машины в целом на примере одного из режимов с использованием программы Mathcad.
рк := 0.83 МПа Ыар:= 18
р0 := 0.23 МПа t2 := 96 С
^е := 8.4 tв1 := 16.5
tв2 := 18
2.4 3.6
Gв:
Gв = 0.667 Ю := -15.5 tk := 19
h1(t0, tпе) := h2(tk, t2)
кг
с
= 1500 1679.5
tлб:= 18 Lкон:= 2
м
м
DOкон= 0.6
алб := 0.01
S1(t0, ^е) := 6.048 S2(tk, t2) := 5.998
S3(tk) := 1.228
Рис. 2. Методика, результаты теплового расчета, определения эксергетических потерь и КПД элементов и цикла холодильной машины в целом на примере одного из режимов с использованием программы МаШсаё
111(10) := 1441.4 54(10*^ := 1259
^ := 264 4 h4(x) := h3(x)
I := h2(tk, 12) - h1(t0, 1пе) _
' ’ у 4 ’ ' Qкон:= &в4.18-(1в2- tв1)
I = 179.5 Qкон= 4.18
qk:= h2(tk, t2) - h3(tk) ДQкон:= -алб-3.14- DOконLкон ^к- 1лб
3
qk = 1.415 х 103 ДQкон= -0.038
^:= Qкон+ ДQкон Qk= 4 142
Qk
G :=---
^ G = 2.927 х 10-3
N := & (112^, t2) - М(10, tпе)) N = 0 525
0 := 111 (10)
Q0:= &- q0
q0 := h11(t0) - h3(tk) q0 = 1.177 х 103
Q0= 3.445
Q0 е :=----
: N Qk1:= ^ NQk1= 3.971
е = 6 557 п^:= 0.9
^- Qk^ -100
D :=----------!----- ^ nk := 0.83
Qk D = 4.141 1
tk = 19
N
Еп:=---------- Еп = 0.703 Т0:= 273 + tнaр
г^-пг^ Т0 = 291
Де21 := (h2(tk, t2) - М(10, 1пе)) - , t2) - Б1(10 , 1пе)) Де21 = 194.05
1
^:= -Де21 + (h2(tk, t2} - 1п1(10 , 1пе))
пМ^ Гk = 46.245 12) - 53(11
Де34:= h3(tk) - h3(tk) - T0(S3(tk) - 54(Ю, tk))
Де23:= h2(tk, 12) - h3(tk) - T0(S2(tk, 12) - S3(tk)) Де23 = 2703
Де34= 9.021
Де41 := h3(1k) - М(10 , 1пе) - , 1k) - S1(t0, 1п0))
Де41 = 157.999
Гдв := 0.1
1
^2^, 12) - Ь1 (10, 1пе))-——
Гдв = 24.029
Проверка
еобщ:= Гk + Де23 + Де34 + Де41 _ „
^ еобщ= 240.295
&
Dk:= Dk= 0.192 ^ Де23-&
Еп Dкон :=----- Dкон = 0.112
Еп
Врв - Де34-& Де4Ю
р Еп Врв = 0.038 Е^=—— Eq = 0.658
100 - Dk■ 100 , ,
пek:=--------------- (1 - Dk- Dкон- Dрв]- 100
90 пek = 0.897 I ^
1 перв:=---------------------------
Eq- 0.05 100 - Dk■ 100 - 11.2
пеисп:=--------:----
Eq пеисп= 0.924 перв = 0.945
Рис. 2. Окончание
При проведении расчетов было установлено, что в одном из режимов нарушено эксерге-тическое неравенство. Причиной является низкая температура воды, охлаждающей компрессор, и создание за счет этого опасного режима работы компрессора, связанного с возможной конденсацией пара на холодной поверхности цилиндра. В дальнейшем из анализа этот режим исключен. Результаты расчетов используются для установления адекватности программного кода.
На рис. 3 приведена блок-схема программы, которая построена на основе приведенной методики с использованием зависимостей необходимых термодинамических параметров, полученных при помощи программ Соо1Раск, Mathcad.
Рис. 3. Блок-схема программы
В соответствии с блок-схемой разработана программа на языке Visual Basic и проведен численный эксперимент для тех же рабочих режимов.
Результаты численного эксперимента для одного из режимов представлены в виде интерфейса программы на рис. 4.
^ Исходные Результаты расчета |
Рк, МПа 075 Ц сИ [1Ш55 Ы,кВт ^^^^^^ЦИ108,248
Р0. МПа 0203 |Щ ^Щ23 д0,и/кс
12, С 100 МЦь1Л [шмэ 00,кВт
1нар, С 17 Ь2\,и/к рЗтЗЗе ере
1пе, С М1, С 5.63 Тз М.Шкс И2,М/кс рбЙгЗ ОкЬкВт □ МЯ5И1
№2, С 6и,мЗ/с Шз 0775 И31,Ы/1< ЬЗ,М/кс |24^23 |^Пб8 Еп,кВт 51 ,Ы/кд 16,07561 |11,0347 '
пк 083 ЩЩ |21ТЭ1Э «2.Ы/кд
ndk П9 дк,Ы/кс р5^6 $3,к*1/кд 11,00127 •
1_кон,м 2 Б,мЗ/с |^и2073 Т0.К
00коп,м 06 Щк,кВт |?58233 е21,Шк
Рис. 4. Интерфейс программы
Определение адекватности разработанной программы было проведено путем сравнения результатов численного эксперимента и расчетов по методике.
Графически результаты сравнения представлены на рис. 5 в виде зависимости эксергети-ческой холодопроизводительности от температуры окружающей среды.
ш I"! й"" ■Едщ
--Щ---Щ---Щ — □ Ецр.
ллл
14,5
15,2
16,83
Рис. 5. Результаты сравнения значений эксергетической холодопроизводительности, полученных в результате численного эксперимента и расчетов по методике
Отклонение сравниваемых величин находится в пределах (4,3.. .6,7) %, что говорит о достаточной точности разработанной программы.
Результаты исследования и их обсуждение
На основе результатов численного эксперимента построены диаграммы потоков и потерь эксергии для четырех режимов работы холодильной машины.
Диаграммы потоков и потерь эксергии позволяют судить о степени участия каждого элемента холодильной машины в потере вводимой эксергии (мощности электродвигателя).
При анализе диаграмм потоков и потерь эксергии (рис. 6) было установлено, что в одном из опытов (^о = -18,9, Ъ,с = 14,5) значения потерь эксергии в компрессоре минимальны вследствие приближения процесса сжатия к адиабатическому. В связи с этим полученное значение потерь, равное 25 %, может быть рекомендовано как минимально возможное для данного типа компрессоров.
Е2
120 100 80 60 40 20 0
1 2 3 4 5 6 7 8
Е3
120 100 80 60 40 20 0
1 2 3 4 5 6 7 8
Е4
120 100 80 60 40 20 0
1 2 3 4 5 6 7 8
Рис. 6. Диаграммы потоков потерь эксергии для четырех режимов
Характер зависимости на рис. 7 соответствует данным, приведенным в [1] для АК-2АВ-20/10. Отклонение сравниваемых величин находится в пределах 10 %.
▲ Т|е(^) ^“Линейный (ґ0)
-20 -18 -16 -14 -12 -10
Рис. 7. Зависимость эксергетического КПД цикла от температуры кипения
Анализ зависимости на рис. 8 показал, что рекомендованные значения потерь в конденсаторе наблюдаются при очень малой разности температур воды на входе и выходе из конденсатора. Таким образом, действительные потери в конденсаторе не соответствуют рекомендован-
-------70
50
40
-------30
20
10
0
ному температурному напору по воде при ее одноразовом использовании. Если предположить, что характер зависимости не меняется, то при рекомендованных значениях разности температур потери в компрессоре составят от 40 до 50 %. Причиной расхождения в величинах потерь является наличие дополнительного термического сопротивления на теплообменной поверхности.
50
45
40
35
30
25
20
15
10
5
0
Акд(ЛГводы)
Логарифмический Dкд(ATводы)
0 1 2 3 4 5
Рис. 8. Зависимость потерь в конденсаторе от температурного напора охлаждающей воды
В подтверждение данного предположения была получена зависимость изменения конденсационного температурного напора от действительных потерь (рис. 9). Как видно из графика, при рекомендованном интервале значений конденсационного температурного напора потери в конденсаторе будут в пределах 30 %.
*■ АТкдфкд)
—Логарифмический
А^дфкд)
Рис. 9. Зависимость конденсационного температурного напора от потерь в конденсаторе
Потери в регулирующем вентиле, определяющие эксергетический КПД в пределах
0,97...0,98, соответствует данным, приведенным на рис. 10 [1].
^.в (Рк1Р00)
~1Л°.в (РфР0)
Рис. 10. Зависимость потерь в регулирующем вентиле от степени сжатия
Заключение
Эксергетический анализ экспериментальной одноступенчатой аммиачной холодильной машины при низкой температуре охлаждающей воды, обеспечивающей приближение к адиабатическому сжатию, позволил установить предельные значения потерь в поршневых компрессорах.
Анализ работы конденсатора позволил сделать заключение о загрязненности теплообменной поверхности.
Полученные значения потерь в элементах экспериментальной установки могут быть использованы при анализе действующих холодильных установок.
СПИСОК ЛИТЕРА ТУРЫ
1. Холодильные машины: Учеб. / Под ред. Л. С. Тимофеевского. - СПб.: Политехника, 1997. - 99 с.
2. Бродянский В. М., Фратшер В., Михалек К. Эксергетический метод и его приложения. - М.: Энерго-атомиздат, 1988. - 287 с.
3. Галимова Л. В. Лабораторный практикум по дисциплине «Теоретические основы энергоресурсосбережения в химической технологии». - Астрахань: Изд-во АГТУ, 2005. - 87 с.
Статья поступила в редакцию 31.01.2008
MODELLING AND EXERGY ANALYSIS OF THE ONE-STAGE AMMONIA EXPERIMENTAL REFRIGERATION MACHINE
L. V. Galimova, A. A. Kamnev, O. O. Lazarenko, C. Guidy Tonyon
The efficiency estimation of refrigeration machine units is determined on the basis of exergy efficiency definition which shows a degree of plant approach to an ideal. The purpose of matter is the analysis of the cycle of experimental refrigeration machine in operating conditions. For achievement of this aim the parameters of the working refrigeration machine have been taken, the incorporated technique of thermal design and exergy losses definition in elements of the refrigeration machine is developed and according to it calculations are carried out, the block diagram and the program code for thermal design, exergy losses definition in elements of the machine and in the cycle as a whole is developed, numerical experiment is carried out and adequacy of the program code is established, schemes of exergy flows and losses in elements of the circuit are constructed, the analysis of the received dependences is made and qualitative ex-ergy reasons are given. Exergy analysis of the experimental one-stage ammonia refrigeration machine in conditions of low temperature of the cooling water, providing approximation to adiabatic compression, allows to define limiting values of losses in piston compressors. The analysis of the condenser work helps to make the conclusion about impurity heat exchanging surface. The received values of losses in elements of experimental plant can be used in the analysis of real refrigeration machine
Key words: the exergy analysis of the one-stage ammonia experimental refrigeration machine, exergy losses, cooling capacity, schemes of exergy flows.
