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

Совершенствование методики проектирования систем холодоснабжения открытых спортивных сооружений сезонной эксплуатации

Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»

№ 2, 2014

УДК 621. 59/56
Совершенствование методики проектирования систем холодоснабжения открытых спортивных сооружений сезонной эксплуатации
Проф. Эглит А. Я. fil24725@yandex.ru Киссер К.В. kisser90@list.ru
Университет ИТМО 921002, Санкт-Петербург, ул. Ломоносова, 9
В настоящий момент калорический расчет открытого ледового поля (катка) ведется по энциклопедическому справочнику «Холодильная техника», в свою очередь возникает ряд вопросов касающихся относительно коэффициента теплоотдачи, скорости ветра, температуры воздуха, влажности воздуха и коэффициента обратного величине испарения. В данной статье проведен анализ влияния климатических параметров для г.Алматы (Республика Казахстан) таких как скорость ветра, солнечная радиация, влажность воздуха, температура воздуха, которые непосредственно влияют на суммарный тепловой поток, предлагается альтернативная методика расчета теплопритоков для открытых искусственных катков сезонной эксплуатации, а также рассмотрен расчет тепловой нагрузки на холодильное оборудование. Данные полученные в ходе исследования могут быть использованы для теплового расчета открытого ледового поля проектными организациями.
Ключевые слова: тепловой расчет, суммарный теплоприток, открытый искусственный каток.
Improving methods of designing cooling systems outdoor sports facilities seasonal operation
Prof. Eglit A. Ya. fil24725@yandex.ru Kisser K.V. kisser90@list.ru
University ITMO 191002, Russia, St. Petersburg, Lomonosov str., 9
Currently caloric calculation of open ice field (rink) is performed based on Encyclopedic Handbook "Refrigeration", in turn, these are number of questions raised concerning relative heat transfer coefficient, wind speed, air temperature, air humidity, and an inverse magnitude evaporation coefficient. In the given article an analyzes of the influence of climatic parameters for the Almaty city (Kazakhstan) such as wind speed, solar radiation, air humidity, air temperature, which directly affect the total heat flow is performed, an alternate design procedure of heat influxes for artificial open ice rinks with seasonal operation is offered as well as that a thermal load calculation for refrigeration equipment is considered. The data received during research can be used for thermal calculation of open ice field by design organizations.
Keywords: thermal calculation, the total heat gain, open artificial ice rink.
В настоящее время открытые искусственные катки рассчитывают по энциклопедическому справочнику «Холодильная техника» [3], но при расчете

Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»

№ 2, 2014

теплопритоков по данному справочнику возникает ряд вопросов о коэффициенте теплоотдачи, скорости ветра, температуре воздуха, влажности и коэффициенте обратном величине испарения.
Согласно [3], расчет теплопритоков принято рассчитывать по следующим формулам:
1) из окружающего воздуха (без учета конденсации влаги): 2)

где: α=6–8 ккал/м2·час – коэффициент теплоотдачи ко льду от воздуха при
скорости его не выше 2 м/с [3]; tв – температура воздуха над поверхностью льда, ºС; tл – температура льда на поверхности катка, ºС. Температура воздуха над поверхностью льда рассчитывается по формуле:

С учетом естественной конвекции воздуха над ледовым полем конвективный коэффициент теплопередачи αконв может быть также рассчитан по следующей формуле [1] :
Коэффициент теплоотдачи ко льду от воздуха α , и в каком диапазоне скоростей, то можно предположить, что расчет коэффициент теплоотдачи ко льду от воздуха мог быть рассчитан по следующей формуле [4]:

где: υ – средняя скорость движения воздуха, м. Так же для открытых катков с учетом вынужденной конвекции, вызванной ветром над ледовым полем, коэффициент теплоотдачи ко льду от воздуха рассчитывается по формуле [1], [4]:

где: в – коэффициент теплопроводности воздуха, Вт/м·К; в – кинематическая вязкость воздуха, м2/с;

Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»

№ 2, 2014

в – средняя скорость ветра рассчитываемого периода, м/с; l – линейный размер охлаждающей плиты в направлении движения воздуха (принимается длина ледяного поля), м; Согласно [9] в существующих методиках предлагаются подобные эмпирические формулы для расчета конвективного коэффициента теплоотдачи, которые можно привести к следующим зависимостям:

– для свободной конвекции:
αконв = А· Δt n,
где А и n – безразмерные коэффициенты, зависящие от параметров воздуха в зоне ледового поля;
– для вынужденной конвекции:
αконв = A·V0,8·x-0,2
где х – линейный размер поверхности, м. Для того чтобы увидеть как влияет скорость ветра и температура окружающего на тепловой поток ко льду был построен график (рис.1), где видно что с увеличением температуры воздуха и скорости ветра над катком теплоприток заметно возрастает.

300 q к2он7д0., Вт/м2
240 210 180 150 120
90 60 30 0

υ= 8 м/с

υ= 7 м/с

υ= 6 м/с

υυυυυυυυυ=========

5 м/с 344,,55м/ммс//сс 3 м/с 2,5 м/с 2 м/с 1,5 м/с 1 м/с

υ= 0,5 м/с

8 9 10 11 12 13 14 15 16 17 18 19 20 tв, ºС

Рис.1 Теплоприток от окружающего воздуха qk=f(tв,v)

3) от грунта или воздушного пространства с нижней стороны катка:

где: k=0,4–0,5 ккал/м2·час·ºС – коэффициент теплопередачи изоляции основания катка [1],[2],[3];
t – температура грунта или воздуха под основанием катка, ºС;

Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»

№ 2, 2014

tк – температура песка, бетона или эвтектика (в аккумуляционных плитах), равная средним температурам кипения холодильного агента, циркулирующего рассола или замерзания эвтектического раствора.
4) от солнечной радиации с учетом 50% тепла, отражаемого от поверхности льда:

где: qR = 550 ккал/м2·час — расчетное напряжение солнечной радиации для летнего периода при падении лучей на горизонтальную поверхность.
Согласно [5], [6], [7] для летнего периода qR ≈ 500 Вт/м2 (Казахстан, г.Алматы) для 1500 ч, а по данным измерений для конец марта-начало апреля 1989г. [8], qR ≈ 730 Вт/м2 для 1500 ч. На период с 1300 ÷1400 для летнего периода по данным [5], [6], [7] qR ≈ 720 Вт/м2, а для весеннего периода qR ≈ 900 Вт/м2 [8]. Соответственно можно сделать вывод о том, что данные в СНиП и СП являются не совсем точными(см.рис.2,3).

q, Вт/м
900 2 800 700 600 500 400 300 200 100
0
4 5 6 7 8 9 10 11 12 13 14 15 16 17 1τ8, ч19ас прямая солнечная радиация

q, Вт/м 1000 2

900

800
qпад
700

600 qпог
500

400

300
qотр
200

100

8 9 10 11 12
Падающая

с13ол1н4 е1ч5н1а6я17р1а8ди19ац20иτ2я,1ч2а2с

Рис.2 Солнечная радиация, поступающая на горизонтальную
поверхность

Рис.3 Солнечная радиация (конец марта– начало апреля 1989г)

5) от конденсации и замерзания влаги:

Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»

№ 2, 2014

где: δ = 25—30 кг/м2·час — коэффициент, обратный величине коэффициента испарения;
х — влагосодержание воздуха с учетом влажности его над поверхностью льда, кг/кг;
х " — влагосодержание воздуха в пограничном слое при средней температуре поверхности льда, к г / к г ;
r = 680 ккал/кг — теплота конденсации водяных паров с учетом замерзания влаги на поверхности льда.
Влагосодержание воздуха с учетом влажности его над поверхностью льда:

где: cра – удельная теплоемкость сухого воздуха при постоянном давлении при 0ºС (1 кДж/(кг·ºС);
cрv – удельная теплоемкость паров воды при постоянном давлении при 0ºС (1,805 кДж/(кг·ºС);
hg – удельная энтальпия насыщенных водяных паров при 0ºС (2501 кДж/кг); h – удельная энтальпия воздуха, кДж/кг; t– температура наружного воздуха, ºС.
Аналогично были построены графики для теплопритока от конденсации и замерзания влаги, из которых видно влияние коэффициента обратному коэффициенту испарения и влажности воздуха в диапазоне температуры воздуха от 8÷20°С.

Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»

№ 2, 2014

q (конд.+зам) , Вт/м2

145 135 125 115 105
95 85 75 65 55 45 35 25 15
8

δ=30 кг/м2·час δ=28 кг/м2·час δ=27 кг/м2·час δ=25 кг/м2·час
tв, ºС 12 16 20

Рис.4 Теплоприток от конденсации и замерзания влаги qkонд+зам= f(tв, φ,δ) при φ=60 %

q (конд.+зам) 250 , Вт/м2
225

φ=90 %

200 φ=75 % 175

150 φ=60 % 125

100 φ=45 %
75

50

25
0 8

tв, ºС 12 16 20

Рис.5 Теплоприток от конденсации и замерзания влаги qkонд+зам= f(tв, φ,δ) при δ =30 кг/м2·час

Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»

№ 2, 2014

Рассмотрим расчет тепловой нагрузки на холодильное оборудование открытого катка для города Алматы (см. табл.1,2,3,4,5).

Таблица 1

Теплоприток из окружающего воздуха (без учета конденсации влаги)

Октябрь

Температура льда tл, ⁰С Температура наружнего воздуха tн, ⁰С
Скорость ветра, м/с
Коэффициент теплоотдачи ко льду α , Вт/м2·⁰С
Удельный теплоприток из окружающего воздуха qконв ,Вт/м2

-3 26 0 2,625
76,374

Октябрь
-3 22 0

Октябрь (ЭС ХТ) -3 18 1

2,529

8,127

63,485 168,026

Таблица 2

Теплоприток от солнечной радиации с учетом 50% тепла, отражаемого от

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

Октябрь

Октябрь

Октябрь (ЭС ХТ)

Максимальная солнечная qR,Вт/м2
Удельный теплоприток от радиации qр, Вт/м2

радиация солнечной

216,318 108,159

216,318 108,159

216,318 108,159

Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»

№ 2, 2014

Таблица 3 Теплоприток от конденсации и замерзания влаги из воздуха

Температура наружного воздуха tн, ⁰С Относительная влажность наружного

воздуха, %

Относительная влажность льда, %

Влагосодержание

воздуха

в

пограничном слое при средней

температуре поверхности льда х", кг/кг

Влагосодержание воздуха с учетом

влажности его над поверхностью льда х, кг/кг

Температура льда, 0С

Теплота конденсации водянных паров с

учетом замерзания влаги на

поверхности льда r, ккал/кг

Коэффициент обратный величине

испарения δ, кг/м2·час

Удельный теплоприток от конденсации

и замерзания влаги из воздуха qконд+зам ,Вт/м2

Октябрь 26 49 100 3
10,164 -3 680
25
141,647

Октябрь 22 49 100 3
10,729 -3 680
25
152,801

Октябрь (ЭС ХТ)
18 64 100 3
9,998 -3 680
25
138,366

Суммарный коэффициент теплопритоков на 1 м2 площади катка

Октябрь Октябрь

Температура наружнего воздуха tн, ⁰С

Удельный

теплоприток

из

окружающего воздуха qконв,Вт/м2 Удельный теплоприток от солнечной

радиации qр, Вт/м2

Удельный теплоприток от конденсации

и замерзания влаги из воздуха qконд+зам ,Вт/м2

∑q, Вт/м2

26 76,374 108,159
141,647 326,18

22 63,485 108,159
152,801 324,45

Таблица 4
Октябрь (ЭС ХТ)
18 168,026 108,159
138,366 414,55

Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»

№ 2, 2014

Результаты исследований по г. Алматы

Таблица 5

Показатель

Методика

Предлагаемая методика

Составляющие тепловой нагрузки на лед:

Расчетная температура наружного воздуха, 0С

18

22 26

Тепловая нагрузка на лед, Вт/(м2·К)

414,55 324,45 326,18

Параметры работы холодильной установки при tл= -30С,хладагент- R404А:

Эффективный теплоотдачи, Вт/(м2·К)

коэффициент

20

13 12

Относительная температура трубной решетки

0,6022

0,6838 0,6474

Температура хладоносителя, 0С

-17 -15 -19

Температура испарения, 0С

-22 -20 -24

Температура конденсации, 0С

33 37 41

По данным результатов видно, что температура кипения повышается на 2К и температура конденсации повышается на 4К. Тепловая нагрузка на лед для расчетной температуры воздуха 18°С превышает допустимый суммарный теплоприток при открытом типе в условиях мягкой зимы без морозов 200— 300 ккал/м2·час практически в 1,5 раза. Из этого следует, что запуск открытого катка в октябре месяце не целесообразен, в связи с высокими энергозатратами. Если рассматривать рассчитанный суммарный теплоприток по предложенной нами методике для октября, то начало работы открытого катка, возможно.
Вывод: Результаты исследования предлагаются для использования проектными организациями.
Список литературы
1. Рекомендации по проектированию инженерного оборудования искусственных катков (хладотехническая часть). ЛенЗНИИЭП Госгражданстроя СССР. Л., 1972.
2. Лихтенштейн Э. Л. Исследование температурного поля искусственного катка. Автореферат диссертации на соискание степени канд. техн. наук. Новосибирск, 1970.

Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»

№ 2, 2014

3. Холодильная техника: Энциклопедический справочник. Госторгиздат, 1961. С. 378–380.
4. Различные области применения холода/Под ред. А.В. Быкова.-М.: Агропромиздат, 1985.-272с.
5. СНиП РК 2.04-01-2010 «Строительная климатология», Изд.: Агентство по делам строительства и жилищно-коммунального хозяйства, Алматы 2011г.
6. СП 131.13330.2012 Строительная климатология, Москва 2012 7. Погода и климат. Погода в Алматы за 2011г.– http://pogoda.ru.net/climate/36870.htm 8. УДК 621.56/59.004.14 Исследование режимных параметров холодильной
установки катка «Медео» Цой А.П. 9. Богословский В. Н. Строительная теплофизика. Учебник для вузов. – М.: Высшая
школа, 1982.