ЭНЕРГЕТИЧЕСКОЕ, МЕТАЛЛУРГИЧЕСКОЕ И ХИМИЧЕСКОЕ МАШИНОСТРОЕНИЕ
39
УДК 621.565.83
Математическая модель для расчета теплового баланса радиатора холодильной системы, использующей эффективное излучение в космическое пространство
Канд. техн. наук А. П. ЦОЙ1, А. С. ГРАНОВСКИЙ2
1teniz@bk.ru, 2granovskiy.a@ramЫer.m
Алматинский технологический университет Республика Казахстан, 050012, г. Алматы, ул. Толе би, 100
В работе представлена методика расчета теплового баланса пластинчатого радиатора для холодильной системы, использующей эффективное излучение в космическое пространство, полученная путем аналитического решения уравнений теплового баланса. Представлены результаты моделирования процессов теплообмена при различных значениях скорости ветра, расхода хладоносителя, теплопроводности излучающей пластины радиатора. Расчеты показывают, что при малой толщине излучающей пластины радиатора, ее теплопроводность не оказывает значительное влияние на общую холодопроизводительность. Поэтому нет необходимости изготавливать ее из материалов с высокой теплопроводностью. Когда температура хладоносителя выше температуры воздуха, увеличение скорости ветра вызывает повышение холодопроизводительности за счет увеличения конвективной теплоотдачи. Существует некоторое максимальное значение скорости хладоносителя, после достижения которого, холодопроизводительность перестает увеличиваться. Представленные решения могут быть использованы при проектировании конкретных моделей радиаторов, а также при автоматизации процесса управления холодильной системой, использующей эффективное излучение в космическое пространство.
Ключевые слова: радиатор, эффективное излучение, ночное радиационное охлаждение, нетрадиционные способы охлаждения, холодильная техника, естественное охлаждение.
Mathematical model of heat balance calculation for the radiator of refrigeration system with effective radiation into space
Ph. D. A. P. TSOY1, A. S. GRANOVSKIY2
1teniz@bk.ru, 2granovskiy.a@rambler.ru
Almaty Technological University, Republic of Kazakhstan, 050012, Almaty, Tole bi str., 100
The method of calculating heat balance of the plate radiator for refrigeration system with effective radiation into space is presented. The method has been obtained by analytical solution of the heat balance equations. The results simulating heat transfer processes at different values of wind speed, coolant flow and thermal conductivity of the radiating fins are presented. Calculation shows that if the radiating fins of the radiator are thin its thermal conductivity does not have a significant impact on the overall cooling capacity. Therefore, there is no need to use high thermal conductivity materials. When the temperature of coolant is higher than the temperature of the air, increase in wind speed causes an increase in the cooling capacity due to increase of the convective heat transfer. There is a maximum value for the coolant speed at which the cooling capacities cease to increase. The solutions presented can be used in the design of specific models of radiators as well as in the automated management of the refrigeration system with effective radiation into space.
Keywords: radiator, effective radiation, nocturnal radiative cooling, non-traditional ways of cooling, refrigerating engineering, natural cooling.
Введение
Основным элементом системы охлаждения, использующей эффективное излучение, является радиатор [1, 2]. От его конструктивных характеристик в конечном итоге зависит эффективность работы всей системы.
Существует множество различных конструкций радиаторов. Ряд авторов предлагает использовать открытые радиаторы, в которых жидкость хладоноси-тель непосредственно контактирует с атмосферным
воздухом и охлаждается одновременно за счет радиационного и испарительного охлаждения [3]. Другие авторы предлагают использовать радиаторы, в которых хладоноситель проходит по закрытым каналам, образованным трубками или пластинами [4-6]. Радиаторы последней конструкции можно использовать для получения сравнительно низких температур в условиях холодного климата, когда во всей системе циркулирует антифриз.
Экспериментальные испытания проведены лишь для нескольких возможных конструкций радиаторов в очень ограниченных диапазонах температур [7].
Для выбора оптимальных параметров радиатора необходима математическая модель, описывающая изменение его холодопроизводительности в зависимости от множества факторов.
Известны исследования, в которых приводятся простейшие способы расчета холодопроизводительнос-ти радиатора пластинчатой конструкции [8, 9]. Однако, в данных работах вопрос о влиянии различных конструктивных особенностей на холодопроизводитель-ность радиатора изучен не достаточно подробно. В частности крайне важной задачей является подбор оптимальных по техническим и экономическим показателям материалов радиатора [10]. Неправильный выбор может неоправданно увеличить стоимость всей конструкции.
В данной же работе разрабатывается математическая модель, которая наиболее точно отражала бы все особенности конструкции радиатора и работала бы в широком диапазоне температур, как окружающей среды, так и поступающего в радиатор хладоносителя.
Как и в работах, ссылки на которые приведены выше, разрабатываемая модель описывает процессы теплообмена в радиаторе пластинчатой (плоской) конструкции. Выбор данной конструкции не случаен. Во-первых, данная конструкция имеет наименьшее сопротивление теплопередаче от хладоносителя через излучающую пластину окружающей среде. Поэтому, ожидается, что ее холо-допроизводительность будет весьма высокой. Во-вторых, подобная конструкция проста в изготовлении, и ее применение является наиболее оправданным с экономической точки зрения. В-третьих, провести расчет теплового баланса для данной конструкции сравнительно просто. В дальнейшем возможна разработка математических моделей для радиаторов более сложной конструкции.
При помощи разрабатываемой модели требуется определить влияние следующих параметров на холодопро-изводительно сть:
— расход хладоносителя;
— скорость ветра;
//////
Ти ав"в
Ы
Ы
— конструктивные особенности радиатора (материал излучающей поверхности).
При этом математическая модель должна давать точные результаты в диапазоне температур от -40 до 50 оС.
Объект и метод исследования
Исследуемая модель радиатора (рис. 1) имеет следующую конструкцию. Хладоноситель протекает по каналу образованному пластинами 2 и 4. Верхняя пластина 2 покрыта слоем селективного покрытия 1. Данный слой необходим для увеличения излучательной способности этой поверхности.
Ниже пластины 4 расположен слой теплоизоляции 5, снижающий теплоприток от окружающей среды к хла-доносителю. Под слоем теплоизоляции расположен корпус 6. Остекление радиатора отсутствует.
Излучающий слой 1 обращен к ночному небу. При работе радиатора между этим слоем и окружающей наблюдается радиационный теплообмен двпрад, за счет которого в основном происходит охлаждение хладоносителя в радиаторе. Расчет qвпрад в данном случае предлагается вести, как при процессе теплопередачи радиацией между двумя параллельными плоскостями. При этом слой атмосферного воздуха заменяется условной плоскостью с температурой Т . Величину Т можно рассчитать по известным методикам [9, 11].
Излучающий слой 1 также взаимодействует с атмосферным воздухом, в результате чего возникает конвективный тепловой поток qвп .
* конв
У нижней стенки корпуса наблюдается только конвективный теплообмен qнп .
* конв
В общем виде уравнение теплового баланса радиатора имеет вид:
е = е + е = F ^ + q) =
^ х ^ вп ^ нп ч-*вп 1 ни'
= F + qвп + qнп ), (1)
41 кони 1 пал 1 конв'" 4 '
т. е. теплота, подводимая с хладоносителем, уходит в окружающую среду через верхнюю излучающую пластину и через нижний слой теплоизоляции.
При составлении математической модели были приняты следующие допущения:
— температура излучающей поверхности во всех точках одинакова;
— конденсация влаги на поверхности радиатора отсутствует.
Для определения удельного теплового потока через верхнюю пластину радиатора qвп составим следующую систему уравнений:
=ах (тх ~тз );
Яви =
Тз - тх
к1/Х1+к2/Х2 = <с Т - Т )+)<<)Т 4 -Тн4).
(2)
Рис. 1. Схема для расчета теплового баланса радиатора
В данной системе три уравнения и три независимых переменных q , Т и Т . Решение системы осложнено на-
А -'вп х н
личием уравнения 4-ой степени.
Для решения системы выразим Т1 из первых двух уравнений как функцию от Т и q . Подставив получен-
ную для Т формулу в 3-е уравнение системы, получим многочлен 4-ой степени вида:
А q4 + В q3 + С q2 + D q + Е = 0,
1 вп 1 вп 1 вп ^кп ■
где A = 5с04; B = - 45сГ 04;C = 65сГ2 02;
£> =
^.ВП ВП 7 _,В1
ав tfB hl а
ч
а я
^-48о213 0 - 1
L =
BC D
- +—;
B
8A3 2A2 A '
,, 3B4 B2C BD E
M =-- +----- + —;
256A 16 A 4 A A
K 2
N =---M ;
12
K3 KM LL
P =--+---;
108 3 8
P P2 N
S =—+J—+—;
2 V 4 27
U = 3ÎS ;
W = --K + U-№, если U = 0; 6
5 -N
W = —K + U +-, если U Ф 0;
6 3U
X = 4k + 2W ;
B X -J-C3K + 2W + 2L / X )
q =---1----.
Чш 4A 2
Чнп = Œ -тау
1
1
D T =
ßx
(3)
Gx Cpx
(6)
1 2
Е = авп Т — авп Т — 5сГ + 5сГ4;
в х в в н X7
0 = — - Ь. _ ^2.
а х ^ 1
Уравнение (3) может быть решено при помощи метода Феррари.
* = - ^+С-;
8А2 А
(4)
Коэффициенты теплоотдачи воздуха ав и хладоносите-ля ах могут быть найдены по стандартным методикам [12].
При проведении расчетов важно учитывать изменение физических свойств атмосферного воздуха и хла-доносителя в зависимости от температуры. Для этого удобно использовать эмпирические формулы. К примеру, для водного раствора пропиленгликоля динамическая и кинематическая вязкость, плотность, температура застывания и теплоемкость могут быть рассчитаны по формулам, приведенным в работе [13].
Для проведения расчетов по предлагаемой методике необходимы данные о погодных условиях, таких как температура воздуха, его влажность, степень облачности. На основе этих данных может быть рассчитана температура ночного неба Тн.
Результаты расчетов и их обсуждение
Для проведения численных расчетов зададимся следующими параметрами радиатора.
Толщина отдельных слоев конструкции: И = 0,0005 м; И = 0,002 м; И, = 0,01 м; И = 0,002 м; И = 0,05 м; И = 0,002 м.
2 7 7 3 7 7 4 7 7 5 7 7 6 7
Длина канала, по которому движется хладоноситель I = 1 м; его ширина а = 1 м.
Коэффициенты теплопроводности слоев: = 0,23 Вт/(мК) для слоя масляной краски;
Х2 = 47 Вт/(м-К) и Х4 = 47 Вт/ (м-К) для стальных пластин;
Х5 = 0,025 Вт /(м-К) для слоя пенополиуретана; = 0,5 Вт/(мК) для слоя пластика.
Относительная излучательная способность слоя масляной краски 5 = 0,93 [13].
Хладоноситель — водный раствор пропиленгликоля с массовой концентрацией £ = 0,5.
Влияние теплопроводности излучающей поверхности на холодопроизводительность. Были проведены расчеты холодопроизводительности радиатора при различных значениях теплопроводности излучающей пластины в диапазоне от 0,025 до 400 Вт/(мК) при температурах ^ = 7 оС, ^ = 20 оС, ^ = 25 оС. Результат расчетов представлен на рис. 2.
Уравнение (4) представляет собой точное решение удельного теплового потока через верхнюю пластину радиатора.
Удельный тепловой поток через нижний слой теплоизоляции q может быть определен по формуле:
q, Вт/м2
(5)
При заданных значениях Т, Т, Тн и известных конструкционных параметрах радиатора по формулам (4) и (5) могут быть рассчитаны тепловые потоки qBn и qm через радиатор. Далее по формуле (1) может быть рассчитана общая холодопроизводительность радиатора.
Понижение температуры хладоносителя на выходе из радиатора:
100 90 80 70 60 50 40 30 20 10 о
к
Вт/(мК)
10
20
30
40
50
Рис. 2. Тепловые потоки в зависимости от теплопроводности излучающей пластины Х2 при температуре хладоносителя t = 25 оС: 1 — Q; 2 — авп 3 — авп ; 4 — анп
х х-х' 1 рад' 1 конв 1 конв
4
х
При увеличении коэффициента теплопроводности от 0 до 2 Вт/(м-К) холодопроизводительность радиатора резко вырастает. Дальнейшее увеличение теплопроводности не оказывает какого-либо влияния на холодопроизводительность.
Влияние скорости ветра на холодопроизводительность
Были проведены расчеты при изменении скорости ветра от 0 до 10 м/с и температуре хладоносителя t = 25 оС и tx = 18 оС (рис. 3, 4). В обоих случаях были приняты значения t = 7 оС и t = 20 оС.
н в
Теплоприток через нижний слой теплоизоляции крайне мал и не оказывает существенного влияния на холодопроизводительность.
В случае, когда t > t (рис. 3), увеличение скорости ветра приводит к увеличению интенсивности конвективной теплоотдачи q™ от излучающей поверхности к воздуху, что в свою очередь вызывает рост холодопро-изводительно сти.
Когда ^ < tE (рис. 4), тепло от воздуха передается излучающей поверхности, о чём говорит отрицательное значение q™ . При этом с ростом скорости ветра холодопроизводительность снижается. В заданных условиях величина конвективного теплопритока становится равной величине радиационного теплового потока при скорости ветра ув > 5 м/с.
На обоих графиках можно наблюдать незначительное уменьшение величины q"1^ с увеличением скорости ветра. Это изменение вызвано незначительным пониже-
нием температуры излучающей поверхности при увеличении коэффициента теплоотдачи авп .
Из рис. 2, 3 и 4 видно, что во всех рассмотренных случаях величина теплопритока через нижнюю пластину радиатора q™^ крайне мала. Следовательно, можно рассмотреть возможность снижения толщины слоя теплоизоляции.
Влияние скорости движения хладоносителя
Были проведены расчеты холодопроизводительнос-ти радиатора при температуре хладоносителя ^ = 25 оС и tx = 18 оС (рис. 5, 6), а также при tH =7 оС и tB =20 оС.
В обоих случаях при увеличении скорости движения хладоносителя v от 0 до 0,05 м/с (G = 0...0,5 кг/с) наблюдается резкий рост холодопроизводительности. Дальнейшее увеличение скорости вплоть до 0,5 м/с (в 10 раз) повышает холодопроизводительность лишь на 1.2 Вт. При скорости vx = 0,05 м/с хладоноситель изменяет свою температуру не более чем на AT = 0,49 0С.
Выводы
В данной работе представлена математическая модель для расчета теплового баланса пластинчатого радиатора холодильной системы, использующей эффективное излучение в космическое пространство. В результате анализа полученной модели важно отметить следующие моменты.
При малой толщине верхней излучающей пластины радиатора, ее коэффициент теплопроводности не оказывает значительное влияние на общую холодопроизво-дительность. Следовательно, излучающую поверхность
q, Вт/М2
е„.вт
В> м/с
Рис. 3. "Тепловые потоки в зависимости от скорости ветра при температуре хладоносителя t = 252 "С: 1 — Q ; 2 — дв
3 — qm ; 4 — q"'
1 коне' 1
рад'
60 40 20 0 -20 -40 -60
q, Вт/м2
\ ----------1 ----------1 ----------1
\
\
\ ~ "" " Ч
v„m/c
Рис. 4. "Тепловые потоки в зависимости от скорости ветра при температуре хладоносителя tх = 18 °С: 1 — Qx; 2 — д
З — дЗ ; 4 — днп
1 конв' 1 конв
рад'
97 96 95 94 93 92 91 90 89
v„ м/с
0,1
0,2
0,3
0,4
Рис. 5. Холодопроизводительность радиатора Qx в зависимости от скорости движения хладоносителя vx при tx = 25 оС
&, Вт
0,1 0,2 0,3 0,4
Рис. 66. .Холодопроизводительность радиатора Q в зависимости от скорости движения хладоносителя vx при t = 18 оЛ
V.M/C
можно изготавливать из наиболее дешевых материалов с низкой теплопроводностью, таких как сталь или полимерные материалы. Применения алюминия или меди не даст значительного увеличения холодопроизводитель-ности, но значительно увеличит стоимость конструкции, что нежелательно.
Если температура хладоносителя близка к температуре воздуха (|ДТ|=0.. .5 К) и скорость ветра мала (ув < 1 м/с), радиационный теплообмен по абсолютной величине превышает конвективный. При увеличении скорости движения хладоносителя v от 0 до 0,05 м/с наблюдается резкий рост холодопроизводительности. Дальнейшее увеличение скорости вплоть до 0,5 м/с практически не влияет на холо-допроизводительность.
Разработанная математическая модель может быть использована при проектировании радиаторов для конкретных условий эксплуатации, а также в системе автоматического управления холодильной системой, использующей эффективное излучение в космическое пространство.
Аналогичные исследования могут быть проведены для радиатора, в котором в качестве хладоносителя используется воздух [14].
Список условных обозначений
h — толщина слоя, м; AT — перепад температуры, К; A, B, C, D, E — коэффициенты уравнения; F — площадь, м2; G — массовый расход, кг/с; ß — тепловой поток, Вт; Т — температура, К; a — ширина, м; cp — изобарная массовая теплоемкость, Дж/(кгК); l — длина, м; q — удельный тепловой поток, Вт/м2; t — температура, C; K, L, M, N, P, S, U, W, X — коэффициенты для решения уравнения 4-ой степени методом Фер-рари; 5 — относительная излучательная способность; X — теплопроводность, Вт/(мК); £ — массовая доля; с — постоянная Стефана-Больцмана, 5,67-10-8 Вт/(м2-К4); ф — относительная влажность, %; v — скорость, м/с.
Индексы
в — воздух; вп — верхняя излучающая пластина радиатора; к — канал, по которому движется хладоноситель; конв — конвективный; н — небо; нп — нижняя поверхность корпуса радиатора; п — водяной пар; рад — радиационный; х — хладоноситель.
Список литературы (References)
1. Цой А. П., Грановский А. С., Бараненко А. В., Эгяит А. Я. Анализ холодильных систем, использующих тепловое излучения земной поверхности в космическое пространство // Вестник Алматинского технологического университета. 2013. № 1. С. 57-59. [Tsoy A. P., Granovsky A. S., Baranen-ko A. V., Eglit A. Ya. The analysis of the refrigerating systems
using thermal radiations of an earth surface in a space. Vestnik Almatinskogo tekhnologicheskogo universiteta. 2013. No 1. p. 57-59. (in Russian)]
2. Цой А. П., Бараненко А. В., Грановский А. С. Моделирование и математическая программа для расчета величины эффективного излучения // Вестник Международной академии холода. 2014. № 1. С. 7-10. [Tsoy A. P., Baranenko A. V., Granovsky A. S. Simulation and a mathematical program for calculating the value of effective radiation. Vestnik Mezhdun-arodnoi akademii kholoda. 2014. No 1. p. 7-10. (in Russian)]
3. Ali A. H. H. Passive cooling of water at night in uninsulated open tank in hot arid areas. Energy Convers. Manag. 2007. Vol. 48, No 1. P. 93-100.
4. Pat. 2011/062649 A1 USA, IPC F24J 2/44. Coaxial tube solar heater with nighttime cooling/Teoh S. T. (USA) — Fill. 23.11.2010; publ. 26.05.2011.
5. Golaka A. R. T., Exell R. H. B. Night radiative cooling and underground water storage in a hot humid climate: a preliminary investigation. Proc. 2nd Reg. Conf. Energy Technol. Towar. a Clean Environ. Phuket, 2003. Vol. 012. P. 10-16.
6. Pat. 6,357,512 B1 USA, IPC F25B 29/00; F04D 13/18; F24J 2/44. Passive heating and cooling system/Baer S. C., Mingenbach W. — fill. 26.06.2000; publ. 19.03.2002.
7. Potentials of night sky radiation to save water and energy in the state of New Mexico: report. Governor Richardson's water innovation fund; Mark Chalom, Bristol Stickney, Kate Snider. — New Mexico, 2006. — 100 p. — PSC #05-341-1000-0035.
8. Dobson R. T. Thermal Modeling Of A Night Sky Radiation Cooling System. Journal of Energy in Southern Africa. 2005. Vol. 16. No 2. P. 20-31.
9. Shuo Zhang. Cooling performance of nocturnal radiative cooling combined with microencapsulated phase change material (MPCM) slurry storage/Shuo Zhang, Jianlei Niu. Energy and Buildings. 2012. No 54. P. 122-130.
10. Цой А. П. Анализ эффективности холодильного радиатора с инфракрасным спектром излучения в космическое пространство // Казахстан-Холод 2014: Сборник докладов международной научно-технической конференции (27 февраля 2014 г.) — Алматы: АТУ, 2014. — С. 114-117. [Tsoy A. P. The analysis of efficiency of the refrigerating heat sink with an infrared spectrum of radiation in a space. Collection of reports of the international scientific and technical conference «Kazakhstan Cold 2014». 2014, p. 114-117. (in Russian)]
11. D. G. Leo Samuel, S. M. Shiva Nagendra, M. P. Maiya. Passive alternatives to mechanical air conditioning of building: A review. Building and Environment. 2013. No 66. P. 54-64.
12. Мухачев Г. А., Щукин В. К. Термодинамика и теплопередача: Учеб. для авиац. вузов. — 3-е изд., перераб. — М.: Высш. шк., 1991. — 480 с. [Mukhachev G. A., Shchukin V. K. Thermodynamics and heat transfer: Textbook. — Moscow, 1991. — 480 p. (in Russian)]
13. Thermophysical properties of brines: Models. — Zurich: M. Conde Engineering, 2011. — 11 p.
14. Hollick J. Nocturnal Radiation Cooling Tests. Energy Procedia. 2012. Vol. 30. P. 930-936.