Например, Бобцов

Определение оптимального количества змеевиков для трубчато-пластинчатых конденсаторов холодильных машин

УДК 621.565.
Определение оптимального количества змеевиков для трубчато-пластинчатых конденсаторов холодильных машин
Лопаткина Т.А.

trusovsa@gmail.com ОАО «ИЭМЗ «Купол», г. Ижевск

Предложено два метода расчета воздушных конденсаторов.

Показано, что разбивка конденсатора на зоны охлаждения, конденсации и

переохлаждения хладагента, а также учет влияния гидравлического

сопротивления на температуру конденсации позволяют определить

оптимальное количество змеевиков. Первый метод расчета, названный

комплексным,

показывает неэффективность конденсаторов,

гидравлическое сопротивление которых снижает температуру конденсации

более чем на 2К. Второй инженерный метод целесообразно применять, если

гидравлическое сопротивление конденсатора снижает температуру

конденсации не более чем на 2К.

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

Definition of optimum quantity of circuits for fin-and-coil condensers of refrigerators

LopatkinaT.A.

trusovsa@gmail.com Open joint-stock company «IEMZ «Kupol», Izhevsk

It is offered the exact and simplified methods of calculation of the vertical air condenser. It is shown that condenser breakdown on zones of cooling, condensation and subcooling of a refrigerant, and also the account of influence of hydraulic resistance on condensation temperature allow to define optimum quantity of circuits. The exact method of calculation shows an inefficiency of condensers, which hydraulics resistance exceeds 2K. The

simplified method of calculation is expedient for applying, if hydraulic resistance reduces temperature of condensation no more than on 2 С.
Keywords: the air condenser, cooling, condensation, subcooling,
hydraulic resistance, a circuit, a calculation error.
При проектировании трубчато-пластинчатых воздушных конденсаторов выделяют три зоны по агрегатному состоянию хладагента: зону охлаждения, в которой парообразный хладагент охлаждается до состояния насыщения; зону конденсации и зону переохлаждения, в которой жидкий насыщенный хладагент охлаждается и становится переохлажденным. Во всех зонах отдаваемое хладагентом тепло воспринимается воздухом.
Чаще всего рекомендуют пренебречь расчетом зон охлаждения и переохлаждения, предлагают определять требуемую площадь теплообмена по зависимостям для зоны конденсации и проектировать конденсатор с запасом площади в 10 – 20 % относительно расчетной 4, 5 . Недостаток такого подхода - отсутствие информации о фактическом размере зон охлаждения и переохлаждения. Кроме того, невозможно определить, как влияет гидравлическое сопротивление конденсатора на его эффективность.
В работах 1,6 предлагается рассчитывать каждую из трех зон конденсатора по отдельности. Однофазные зоны в методах 1, 6 рассчитываются под конкретное количество параллельных змеевиков. Сопоставляя результаты расчета нескольких конденсаторов, отличающихся количеством змеевиков, можно определить их оптимальное количество. Процедура расчета представляется весьма сложной.
Ниже описан метод расчета вертикального конденсатора (рис. 1) и показано влияние гидравлического сопротивления (количества змеевиков) на его эффективность. Исходные данные для расчета:
- параметры воздуха на входе в конденсатор- температура tв1 , С; относительная влажность в1 %; объемный расход Vв , м3/ч; барометрическое давление Рб , Па;
- параметры хладагента - марка; температуры t2 - на входе в конденсатор, tк - конденсации, t3 - на выходе из конденсатора или переохлаждение tа , определяемое как tа = tк - t3 С) ;
- тепловая нагрузка на конденсатор Q , кВт или расход хладагента Ga , кг/с;

- параметры конструкции - размеры фронтального сечения (высота и ширина оребренной части); диаметр и толщина стенки трубок; геометрия и материал пластин.
Требуется определить площадь теплообмена, запас площади, аэродинамическое и гидравлическое сопротивление.
АГ Рис. 1. Схема вертикального воздушного конденсатора. Зоны с разными состояниями хладагента: I – охлаждения, II – конденсации, III – переохлаждения. Количество змеевиков в конденсаторах, шт.: А – 4, Г – 8. Принято, что хладагент движется сверху вниз, а воздух горизонтально, проходя последовательно все ряды трубок конденсатора. Данное направление потоков обеспечивает одну и ту же температуру воздуха на входе в зоны с разными состояниями хладагента, в отличие от 1, 6 . Предварительно задается количество рядов труб, шаг пластин, количество змеевиков и рассчитывается наружная площадь теплообмена Fн , м2. По исходным данным вычисляется массовый расход воздуха Gв , кг/с, далее рассчитывается зона охлаждения газообразного хладагента. Тепловую нагрузку в зоне охлаждения Qохл , кВт определяют по формуле:
Qохл Gа iа i" , (1)
где Gа , кг/с – расход хладагента, заданный или вычисленный на основании исходных данных; ia и i , кВт – энтальпия перегретого хладагента на входе в

конденсатор и энтальпия насыщенного пара хладагента при температуре конденсации tк соответственно.
Предварительно задают массовый расход воздуха Gв охл , приходящийся на зону охлаждения; рассчитывают температуру воздуха на выходе из зоны, среднелогарифмическую разность температур tохл , С, средние температуры воздуха и хладагента в зоне и параметры этих сред при средней температуре в зоне. Вычислив скорость воздуха в живом сечении и скорость хладагента при их средней температуре, находят по известным зависимостям коэффициенты теплоотдачи для воздуха в и для газообразного хладагента охл , Вт/(м2 С), а затем и коэффициент теплопередачи kохл , Вт/(м2 С). Требуемая для зоны охлаждения наружная площадь теплообмена Fн охл , м2 , определяется по формуле:

Fн охл

1000Qохл .
kохл tохл

(2)

Через эту площадь проходит массовый расход воздуха Gв охл пров , кг/с, значение которого используется для проверки предварительно принятого расхода воздуха Gв охл :

Gв охл пров



Fн охл Fн

.

(3)

Используя метод последовательных приближений, необходимо добиться, чтобы проверочный расход воздуха Gв охл пров отличался от принятого Gв охл не более чем на 0,1 %. Как только этого условие выполнится, наружная площадь теплообмена Fн охл для зоны охлаждения хладагента и массовый расход воздуха Gв охл будут в первом приближении определены. Рассчитывается гидравлическое сопротивление хладагента в зоне
охлаждения Ра охл , МПа, как для однофазного потока по зависимостям,
рекомендованным в 3 .

Гидравлическое сопротивление Ра охл вызывает снижение температуры конденсации tк , что сказывается на работе конденсатора и должно быть учтено в расчете. Уточняются наружная площадь теплообмена
Fн охл и массовый расход воздуха Gв охл следующим образом. Зная Ра охл , вычисляют температуру конденсации на выходе из зоны охлаждения, tк охл-к, С: это температура насыщения при давлении (Рк - Ра охл), где Рк , МПа – давление насыщенного хладагента при температуре tк . Вычисляют

энтальпию i охл-к , кДж/кг, насыщенного пара хладагента при температуре tк охл-к и переопределяют тепловую нагрузку в зоне охлаждения Qохл , кВт по формуле:

Qохл Gа iа iо"хл к .

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

Рассчитывают зону конденсации. Начальное значение тепловой нагрузки в зоне конденсации Qк , кВт определяется по формуле:

Qк Gа iо"хл к iо' хл к ,

(4)

где i , кДж/кг – энтальпия насыщенного жидкого хладагента при температуре конденсации tк охл-к , С.
Задают предварительно массовый расход воздуха Gв к , приходящийся на зону конденсации. Вычислив в зоне конденсации коэффициент теплопередачи kк , Вт/(м2 С) и среднелогарифмическую разность температур tк , С, в которой для температуры конденсации принято значение tк охл-к , определяют требуемую наружную площадь теплообмена для зоны конденсации Fн к , м2:

Fн к

1000 Qк kк tк

(5)

и проверочное значение массового расхода воздуха Gв к пров , кг/с:

Gв к пров

Fн к Gв Fн

.

(6)

Расчет площади теплообмена для зоны завершается при выполнении условия (Gв к - Gв к пров) 100 / Gв к 0,1. Гидравлическое сопротивление хладагента в

зоне конденсации Ра к , МПа, определяют по зависимостям из 5 , а коэффициент теплоотдачи от хладагента рассчитывают по формуле для
расслоенного режима из 5 .

Зная гидравлическое сопротивление Ра к , вычисляют температуру насыщенного жидкого хладагента на выходе из зоны конденсации, tк к-по , С
при давлении (Рк - Ра охл - Ра к) , МПа и уточняют наружную площадь теплообмена Fн к и массовый расход воздуха Gв к с учетом того, что на выходе из зоны температура насыщенного жидкого хладагента снизилась до
значения tк к-по . Прежде всего вычисляют энтальпию i к-по , кДж/кг, насыщенного жидкого хладагента при температуре tк к-по и переопределяют тепловую нагрузку в зоне конденсации Qк , кВт:

Qк Gа iо"хл к iк' по .

В среднелогарифмической разности температур tк учитывают температуру конденсации tк охл-к на входе в зону и tк к-по на выходе из зоны. Все теплофизические параметры хладагента берут при средней температуре конденсации в зоне. После теплового расчета зоны конденсации с
переопределенными величинами Qк , tк вновь рассчитывают гидравлическое сопротивление хладагента Ра к и сопоставляют его с полученным в предыдущем цикле. Если новое значение Ра к отличается от предыдущего более, чем на 0,1 % - выполняют новый цикл, основываясь на новом
значении Ра к . Расчет зоны завершают, когда значения гидравлического сопротивления в текущем и предыдущем циклах не будут различаться менее чем на 0,1 %.

Расчет зоны переохлаждения жидкого хладагента подобен расчету зоны охлаждения. Определяют тепловую нагрузку в зоне переохлаждения Qпо , кВт:

Qп о Gа iк' по iп о ,

(7)

где iп-о , кДж/кг – энтальпия переохлажденного хладагента на выходе из конденсатора.

Задав предварительно массовый расход воздуха Gв п-о , приходящийся на зону переохлаждения, определяют температуру воздуха на выходе из
зоны, среднелогарифмическую разность температур tп-о , С, средние температуры воздуха и хладагента в зоне и их параметры при средней

температуре в зоне. При вычислении среднелогарифмической разности
температур tп-о и средней температуры хладагента в зоне используют значение температуры конденсации tк к-по. Находят коэффициент теплопередачи kп-о , Вт/(м2 С) и определяют наружную площадь теплообмена Fн охл , м2 , требуемую для зоны переохлаждения:

Fн п о

1000Qп о .
kп о tп о

(8)

Расчет площади теплообмена зоны переохлаждения заканчивают, когда
будет выполнено условие (Gв п-о - Gв п-о пров) 100 / Gв п-о 0,1, где проверочное значение массового расхода воздуха для зоны переохлаждения Gв п-о пров , кг/с определено по формуле:

Gв п

о пров

Fн п Gв Fн

о

.

(9)

Рассчитывают гидравлическое сопротивление хладагента в зоне
переохлаждения Рп-о , МПа, по зависимостям из 3 . Влиянием Рп-о на результат теплового расчета зоны пренебрегают ввиду его малости.

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

Запас Fто Fн 100 ,
Fто

(10)

где Fто , Fн – наружная площадь теплообмена, имеющаяся согласно предварительно принятой конструкции конденсатора, и требуемая наружная
площадь теплообмена соответственно.

Описанным методом рассчитаны шесть конденсаторов, отличающихся количеством змеевиков, (см. таблицу 1 и рисунок 2). Конденсаторы А и Г имеют 4 и 8 змеевиков соответственно (рис. 1), конденсатор Е с 16-ю змеевиками конструктивно подобен конденсатору Г. Аппараты Б, В и Д с 5-ю, 6-ю и 12-ю змеевиками - условные; они рассчитаны для полноты представления зависимости результатов расчета от количества змеевиков.

Принято, что конденсаторы продуваются воздухом с начальной температурой tв1 = 32 С и относительной влажностью в1 = 40 % в количестве 10000 м3/ч. Требуется отвести тепло от хладагента R134a с параметрами Gа = 0,1378 кг/с, t2 = 63,2 С , tк = 47 С, tа = 5 . Барометрическое давление 100000 Па. Геометрические размеры конденсаторов: высота и ширина оребренной части 800 мм, диаметр медных трубок 9,52 мм, толщина трубок 0,35 мм, ребра алюминиевые толщиной 0,15 мм, шаг ребер 2,5 мм, форма - согласно 2 .
По мере увеличения количества змеевиков уменьшается требуемая площадь теплообмена для зоны конденсации Fн к (см. таблицу 1 и рис. 2). Причина - увеличение температурного напора по мере снижения гидравлического сопротивления в зонах охлаждения и конденсации и падения температуры конденсации Dtк охл и Dtк к.. Требуемая площадь теплообмена для зоны охлаждения Fн охл растет, т.к. снижается скорость хладагента и коэффициент теплопередачи kохл .

Таблица 1

Наименование параметра

Значение параметра для конденсатора АБВГДЕ

1. Тепловая нагрузка, кВт:

- в зоне охлаждения Qохл

2,60 2,59 2,58 2,58 2,58 2,58

- в зоне конденсации Qк

22,82 21,93 21,67 21,48 21,39 21,37

- в зоне переохлаждения Qп-о

1,01 1,04 1,05 1,05 1,05 1,05

2. Скорость воздуха в живом сечении, 8,0 8,0 8,0 8,0 8,0 8,0 м/с

3. Скорость хладагента, м/с:

- в зоне охлаждения

8,8 7,0 5,8 4,4 2,9 2,2

- в зоне переохлаждения

0,43 0,35 0,30 0,22 0,15 0,11

4. Коэффициент теплопередачи, Вт/(м2 С):

39,3 35,6 32,6 28,1 22,4 18,8

- в зоне охлаждения kохл

62,5 61,8 61,6 61,5 61,4 61,4

- в зоне конденсации kк

43,8 40,2 37,2 32,5 26,3 22,3

- в зоне переохлаждения kп-о

5. Требуемая наружная площадь

теплообмена, м2

3,6 3,9 4,2 4,8 5,8 6,8

- для зоны охлаждения Fн охл

45,1 37,0 35,2 34,0 33,4 33,3

- для зоны конденсации Fн к

5,5 3,3 3,1 3,2 3,8 4,3

- для зоны переохлаждения Fн п-о

6. Падение температуры конденсации, С 0,46 0,27 0,17 0,08 0,03 0,02 - в зоне охлаждения фреона, Dtк охл 6,52 2,42 1,25 0,45 0,05 0,00 - в зоне конденсации фреона, Dtк к

7. Запас площади теплообмена, %

-33,5 -9,0 -4,8 -3,5 -5,9 -9,3

8. Доля требуемой площади 16,8 16,4 17,3 19,1 22,2 25,0 теплообмена для охлаждения и переохлаждения в общей требуемой площади теплообмена, %

F, м2

4,5 4
3,5 3
2,5 2
1,5 1
0,5 0

АБ

ВГД

конденсатора

Е Вариант

Рис. 2. Требуемая площадь теплообмена для конденсаторов с разным количеством змеевиков.

F – наружная площадь теплообмена, м2, в зонах: конденсации; - переохлаждения.

- охлаждения; -

Количество змеевиков в конденсаторах, шт.: А – 4, Б – 5, В – 6, Г – 8, Д – 12, Е – 16.

На величину требуемой площади теплообмена для зоны переохлаждения Fн п-о влияет как увеличение температурного напора в предыдущих зонах, так и снижение коэффициента теплопередачи kп-о . В результате действия разнонаправленных факторов наименьшую требуемую площадь теплообмена имеют конденсаторы В и Г с оптимальным числом змеевиков, при котором падение температуры конденсации в двух первых зонах (Dtк охл + Dtк к ) = (0,53 – 1,42) С и скорость хладагента в зоне охлаждения (4,4 – 5,8) м/с.
Если число змеевиков выше оптимального, то однофазные зоны занимают неоправданно большую площадь. Так в конденсаторе Е на них приходится 25 % от общей требуемой площади при доле тепловой нагрузки 14,5 % (таблица 1).
Второй способ увеличения эффективности конденсатора - сокращение змеевиков в однофазных зонах. Благодаря сокращению змеевиков преимущество высокой скорости хладагента в однофазных зонах сочетается с преимуществом низкой скорости в зоне конденсации: увеличиваются коэффициенты теплопередачи в зонах охлаждения и переохлаждения kохл , kп-о и одновременно обеспечивается наибольший температурный напор в зонах конденсации и переохлаждения. В результате уменьшается требуемая площадь теплообмена. Например, за счет сокращения количества змеевиков вдвое в однофазных зонах требуемая площадь теплообмена конденсатора Г (8 змеевиков) может быть снижена на 2 %, а теплообменника Д (12 змеевиков) - на 5,4 %.

Рассмотренный метод расчета трудоемок, поэтому возникает вопрос об инженерном, упрощенном методе, его возможной погрешности и границах применимости.
В упрощенном методе исключены пересчеты требуемой площади теплообмена для зон охлаждения и конденсации, обусловленные снижением температуры конденсации в них. Принято допущение, что температура конденсации в каждой зоне постоянна и соответствует значению на входе в зону. Снижение температуры конденсации в предыдущих зонах учитывается. Таким образом, в зоне охлаждения температура конденсации tк , в зоне конденсации tк охл-к , в зоне переохлаждения tк к-по . Как следствие, в расчете принято скачкообразное изменение энтальпии при переходе от зоны к зоне: i" i"охл-к , i'охл-к i'к-по . Тепловая нагрузка в зонах определяется по формулам 1, 4 и 7. Требуемая площадь теплообмена, как и прежде, определяется по формулам 2 и 3, 5 и 6, 8 и 9 для зоны охлаждения, конденсации и переохлаждения соответственно.
Рисунки 3 и 4 показывают, как снижается точность упрощенного метода расчета по сравнению с комплексным для конденсаторов А – Е.
Упрощенный метод расчета занижает падение температуры конденсации в двух первых зонах (Dtк охл + Dtк к) , С. Занижение тем значительнее, чем выше гидравлическое сопротивление конденсатора. Так (Dtк охл + Dtк к) = 7 С занижено упрощенным расчетом в 1,4 раза, а (Dtк охл + Dtк к) = 0,5 С только на 1,5 % (см. рисунок 3, конденсаторы А и Г). Подобным образом занижается требуемая площадь теплообмена. Погрешность ее расчета для конденсаторов А – Е

(Dtк охл+Dtк к), град. С
Поргешность, %

85
7
0
6
-5
5
4 -10
3 -15
2
-20
1
0 -25
-1 -30

АБ

ВГД

ЕЕ

Вариант конденсатора

АБВГ Вариант конденсатора

Д

Рис. 3. Падение температуры конденсации в зонах охлаждения и конденсации для конденсаторов с

Рис. 4. Погрешность расчета требуемой площади теплообмена.

разным количеством змеевиков.

- точный метод;

- точный метод;

упрощенный метод; -

-

упрощенный метод; варианты

- расчет при постоянной

конденсаторов – см. подписи к рис. температуре конденсации tк.во всех зонах.

2.

упрощенным методом (рисунок 4) вычислена по формуле:

Погрешность Fн упрощ Fн точн 100 , %
Fн точн
где Fн упрощ , Fн точн – требуемая наружная площадь теплообмена по упрощенному и по комплексному методам.
Упрощенный метод расчета можно использовать для поиска оптимального количества змеевиков конденсатора, если суммарное падение температуры конденсации в двух зонах (Dtк охл + Dtк к) не превышает 2 С. В этом случае обеспечивается погрешность расчета требуемой площади теплообмена не более 6 % (см. рис. 3 и 4). При (Dtк охл + Dtк к) более 2 С растущая погрешность упрощенного метода расчета (рис. 4) не позволяет обнаруживать снижение эффективности конденсатора в полной мере. Так для

конденсатора А точным методом по формуле 10 получен недостаток площади теплообмена 33,5 %, а упрощенным – недостаток 3,2 %.
Как правило, конденсаторы рассчитывают при постоянной температуре конденсации tк во всех зонах. На рис. 4 показана погрешность расчета конденсаторов А – Е при tк = const. Погрешность требуемой площади теплообмена в этом случае оказалась в 1,2 – 1,5 раза выше, чем по упрощенному методу, во всем рассмотренном диапазоне гидравлических сопротивлений.
Предложенные методы расчета вертикального воздушного конденсатора с разбивкой на зоны охлаждения, конденсации и переохлаждения хладагента, а также с учетом влияния гидравлического сопротивления на температуру конденсации позволяют определить оптимальное количество змеевиков. Упрощенный метод расчета целесообразно применять, если гидравлическое сопротивление снижает температуру конденсации не более чем на 2 С. Согласно расчетам, количество змеевиков оптимально, если суммарное гидравлическое сопротивление змеевика в зонах охлаждения и конденсации хладагента снижает температуру конденсации на (0,5 – 1,5) С.
Список литературы
1. Бялый Б.И. Тепломассообменное оборудование воздухообрабатывающих установок ООО «Веза» / М.: ООО «Инфорт», 2005. 278 с.
2. Емельянов А.Л., Кожевникова Е.В. Исследование коэффициентов теплоотдачи в воздушных трубчато-пластинчатых калориферах // Холодильная техника -2011, -№7, - С.15-19.
3. Идельчик И.Е. Справочник по гидравлическим сопротивлениям / под ред. М.О. Штейнберга. М.: Машиностроение, 1992. 672 с.
4. Маринюк Б.Т. Аппараты холодильных машин (теория и расчет) / М.: Энергоатомиздат, 1995. 160 с.
5. Теплообменные аппараты холодильных установок: 2-е изд., перераб. и доп. / под общ. ред. Г.Н. Даниловой. Л.: Машиностроение, Ленингр. отделение, 1986. 303 с.
6. Чумак И.Г., Чепурненко В.П., Чуклин С.Г. Холодильные установки: 2-е изд. перераб. и доп. / М.: Легкая и пищевая промышленность, 1981. 344 с.