CC BY
3
ИНЖЕНЕРНЫЕ СИСТЕМЫ В СТРОИТЕЛЬСТВЕ
НАУЧНАЯ СТАТЬЯ / RESEARCH PAPER УДК 697.9
DOI: 10.22227/1997-0935.2023.11.1802-1812
Разработка методики и программы расчета параметров воздуха при кондиционировании и вентиляции
крытых катков
Сергей Владимирович Чуйкин
Воронежский государственный технический университет (ВГТУ); г. Воронеж, Россия
АННОТАЦИЯ
Введение. В помещениях с искусственным льдом, к которым относятся крытые катки и ледовые арены, для предотвращения размягчения льда, образования тумана и выпадения конденсата конвективная составляющая теплообмена должна расходоваться на охлаждение воздуха зоны ледового поля до нормируемых значений. От точности расчета микроклиматических параметров зависит не только комфорт находящихся на поле людей, но и качественные характеристики льда. Актуальна разработка методики и программы расчета указанных параметров при проектировании систем вентиляции и кондиционирования воздуха.
Материалы и методы. Для повышения точности расчета предложена новая методика аналитического определения параметров состояния воздуха в зоне ледового поля крытых катков и ледовых арен. Данная методика отличается от существующих применением итерационного метода расчета при выборе оптимального соотношения расходов воздуха на первой и второй ступенях рециркуляции, что позволяет исключить дополнительное охлаждение притока в холодный период года.
Результаты. Разработанная методика реализована при составлении алгоритма программы расчета искомых параметров и учитывает современные данные численного моделирования температурных полей в обслуживаемой зоне, скорректированный диапазон величины градиента температуры по высоте рассматриваемого помещения и коэффициента массообмена, определяемого с помощью теории тройной аналогии тепломассообменных процессов. Выводы. Предложенные методика и программа расчета дают возможность повысить точность определения темпе-£ ф ратуры, влагосодержания, парциального давления и энтальпии при проектировании систем вентиляции и кондицио-
нирования помещений с искусственным льдом, что способствует сокращению избыточных энергетических затрат при обработке воздуха в центральном кондиционере.
КЛЮЧЕВЫЕ СЛОВА: вентиляция, кондиционирование, /^-диаграмма, процессы обработки воздуха, вентиляция ео щ крытых катков, параметры состояния, автоматизированное проектирование
1
Благодарности. Автор выражает благодарность редколлегии и рецензентам за профессиональный и ответственный О и подход к работе.
>
"(D ф
О -£=
О ££
со >
<° 2;
С\1 5
<л
(Л
.Е о cl"
ДЛЯ ЦИТИРОВАНИЯ: Чуйкин С.В. Разработка методики и программы расчета параметров воздуха при кондиционировании и вентиляции крытых катков // Вестник МГСУ. 2023. Т. 18. Вып. 11. С. 1802-1812. DOI: 10.22227/1997-0935.2023.11.1802-1812
с о
О ф Автор, ответственный за переписку: Сергеи Владимирович Чуикин, ser.chu@mail.ru.
Development of methodology and programme for calculation of air parameters for air conditioning and ventilation of indoor
skating rinks
Sergey V. Chuykin
S a
o E Voronezh State Technical University (VSTU); Voronezh, Russian Federation § ° -
^ ABSTRACT
Z £
co g Introduction. In premises with artificial ice, which include indoor ice rinks and ice arenas, in order to prevent ice softening, fog
— formation and condensate precipitation, the convective component of heat exchange should be spent on cooling the air of the ice
^ • field area to standardized values. Thus, the accuracy of the calculation of microclimatic parameters depends not only on the comfort
O jj of people on the field, but also on the qualitative characteristics of the ice. The presented article is devoted to the development of
O a methodology and a programme for calculating these parameters when designing ventilation and air conditioning systems.
^ S Materials and methods. To improve the calculation accuracy, a new method of analytical determination of air condition
S parameters in the ice field area of indoor ice rinks and ice arenas is proposed. This method differs from the existing ones
X
с
_ by application of iterative calculation method when choosing the optimal ratio of air flow rates at the first and second stages
jj jj of recirculation, which makes it possible to exclude additional cooling of the inflow during the cold season.
U > Results. The developed method is implemented in the algorithm of the programme for calculation of the required parameters
and takes into account the modern data of numerical modelling of temperature fields in the served area, the corrected range
1802 © С.В. Чуйкин, 2023
Распространяется на основании Creative Commons Attribution Non-Commercial (CC BY-NC)
of the temperature gradient value by the height of the room under consideration and the mass transfer coefficient determined using the theory of triple analogy of heat and mass transfer processes.
Conclusions. The proposed methodology and calculation programme make it possible to increase the accuracy of determining temperature, moisture content, partial pressure and enthalpy determination when designing ventilation and air conditioning systems for rooms with artificial ice, which contributed to the reduction of excessive energy costs when processing air in the central air conditioner.
KEYWORDS: ventilation, air conditioning, Id diagram, air treatment processes, indoor ice rink ventilation, condition parameters, computer-aided design
Acknowledgements. The author expresses gratitude to the editorial board and reviewers for their professional and responsible approach to the work.
FOR CITATION: Chuykin S.V. Development of methodology and programme for calculation of air parameters for air conditioning and ventilation of indoor skating rinks. Vestnik MGSU [Monthly Journal on Construction and Architecture]. 2023; 18(11):1802-1812. DOI: 10.22227/1997-0935.2023.11.1802-1812 (rus.).
Corresponding author: Sergey V. Chuykin, ser.chu@mail.ru.
ВВЕДЕНИЕ
Помещения с искусственным льдом, к которым относятся крытые катки и ледовые арены, в зависимости от вида проводимых спортивных и культурно-массовых мероприятий могут иметь различные архитектурно-строительные и планировочные решения, влияющие на подход к организации и проектированию систем жизнеобеспечения. Однако определяющим фактором при создании требуемого микроклимата рассматриваемых зданий является наличие обширной охлаждающей поверхности. В последние годы с помощью современных вычислительных средств был проведен ряд теоретических и экспериментальных исследований формирования тепловлажностного и воздушного режимов данных объектов [1-8], что позволило повысить точность получаемых результатов. Тем не менее, ввиду сложности и трудоемкости составления физико-математических моделей в современных программных средах, по-прежнему наибольшим распространением пользуется графоаналитическая методика, изложенная в трудах О.Я. Кокорина [9, 10]. Поскольку испарение и конденсация влаги при взаимодействии ледового поля с окружающим воздухом протекают по границам кристаллов [11, 12], что объясняется меньшей термодинамической устойчивостью границ по сравнению с основной массой кристалла, необходимо стремиться к заливке относительно монолитного слоя льда с крупной кристаллической структурой.
Ключевые проблемы и задачи при формировании более качественного ледового покрытия рассматриваются в работах [11-15], отмечается, что главным условием при теплотехническом расчете холодильного оборудования для заливки льда служит компенсация им тепловых притоков к поверхности льда от окружающей среды, величина которых определяется интенсивностью радиационно-конвективного теплообмена.
Таким образом, большое значение приобретает точность расчета параметров микроклимата, в значительной степени влияющих на тепловой баланс обслуживаемой зоны помещения и его поверхностей. Вышесказанное показывает возросшую акту-
альность разработки новых методов расчета температуры, влажности, энтальпии и других параметров состояния влажного воздуха, что и является основной целью данной статьи.
МАТЕРИАЛЫ И МЕТОДЫ
Рассмотрим более подробно существующую методику определения полного теплового потока к поверхности льда при ее радиационно-конвектив-ном теплообмене с окружающим воздухом. Принято считать, что теплопритоки можно найти с помощью равенства [15-17]:
От.л.пов От.кон ^ От.луч ^ От,в ^ ОТ,в ^ °г.лвд' (1)
где Q — конвективный приток теплоты от воздуха к поверхности льда, Вт; Qтлуч — приток лучистой теплоты от ограждающих конструкций и поверхностей к поверхности льда, Вт; Отесв — приток лучистой теплоты от естественного освещения, Вт; Огхв — теплопритоки к поверхности льда от осветительных приборов, Вт; Q — теплопритоки
7 ^ т. люд А
от людей на льду и в прилегающих зонах, Вт.
Величина конвективных притоков определятся по уравнению Ньютона - Рихмана, которое для рассматриваемого объекта принимает вид [9, 15-17]:
0г.
í (( - к \
(2)
где Fл — площадь ледового поля, м2; акон — ко- Ц эффициент теплоотдачи ледовой поверхности, Вт/(м2град); ^ — температура воздуха в зоне ледового поля, °С; (л — температура поверхности льда, °С.
Коэффициент теплоотдачи может определяться по эмпирическим зависимостям: • методика АБНМЕ (США):
■■ 3,41 + 3,55 • vë,
(3)
где vл — скорость воздуха у поверхности ледового поля, м/с;
• методика из справочника «Различные области применения холода» А.В. Быкова при вынужденном конвективном теплообмене для открытых катков:
< П i i
4
G Г
5 2
0 со
n СО
1 <
< -ь J со
U -
r i
П о
<3 o <
o i
O n
со со
м со
0
1
СО СП о о
cn
• ) ® 8
л '
00 п
1 Т
s У с о (D Ж
Ы 10
о о
10 10
U W
1803
(О (О
N N
О О
N N
¡г ш
U 3 > (Л С И
m «о
. т-
со щ
?!
Ф Ф
О ig
о о со со
I
о со сч
(Л (Л
.Е о с
Ю о
о Е
fe ° СП ^ т- ^
Е
от °
* А £ w
■S
II
О tn
«1 = 0,037 •
w0'8 • l
(4)
где Хв — коэффициент теплопроводности воздуха, Вт/(м-°С); vb — кинематический коэффициент вязкости воздуха, м2/с; wb — среднемесячная скорость ветра в расчетный период, м/с; l — линейный размер поля в направлении движения ветра, м;
• методика из справочника «Различные области применения холода» А.В. Быкова для крытых катков:
■Ъ
(5)
где 4 — коэффициент влаговыпадения, учитывающий выделение тепла при конденсации влаги на те-плообменной поверхности.
Приток лучистой теплоты от потолка к поверхности льда принято рассчитывать по формуле:
0г.луч Fr Члуч ßo
(6)
где qлуч — плотность лучистого теплового потока к поверхности льда от строительных конструкций и ограждающих поверхностей, Вт/м2; Ротр — коэффициент отражения (зависит от материала поверхности или конструкции).
Плотность лучистого притока теплоты в общем случае определяется по известной зависимости:
q ё
т V (тл4
100
_2
100
■ф
(7)
где сп — приведенный коэффициент излучения, Вт/(м2К4); Т1 и Т2 — температуры взаимодействующих поверхностей, К; ф — коэффициент облучения, учитывает взаимное расположение поверхностей теплообмена.
Приток лучистой теплоты от естественного освещения рассчитывают по формуле:
FOK ^ок Ва
(8)
0тисв = 0,6 ■ FR ■ qocB,
(9)
Необходимый объемный расход подаваемого воздуха определяется из условий равномерного заполнения зоны ледового поля приточным воздухом и создания предпосылок для понижения температуры притока до требуемого значения у поверхности льда. Подвижность воздуха у поверхности льда ограничивается величиной 0,3 м/с, в противном случае в зонах превышения этого показателя качество льда ухудшается. Добиться выполнения этого требования возможно лишь при омывании ледового поля обратными потоками, индуцируемыми направленными струями. Расчет воздухора-спределения при этом принято называть расчетом по обратному потоку, а за нормируемую скорость в обратном потоке принимается максимально допустимая скорость в обслуживаемой зоне.
Исходными данными при расчете микроклиматических параметров воздуха в обслуживаемой зоне помещений с искусственным льдом являются: расход притока (определенный из вышеупомянутых условий); температура и влажность наружного воздуха и воздуха в обслуживаемой зоне, которые зависят от вида проводимых мероприятий и климатической характеристики района строительства.
Расчет всех термодинамических параметров внутреннего воздуха производится путем последовательного решения уравнений (11)-(19), зная любые два из параметров можно определить остальные:
I = 1,005 ■ t + (2500 +1,8 ■ ^ ■ d ■ 10~3, (11)
где I — энтальпия воздуха, кДж/кг; t — температура воздуха, °С; d — влагосодержание воздуха, г/кг.
d = d ' ■ 103 = 622 ■
Ри
[P.-P« У
(12)
где ^ — влагосодержание воздуха, кг/кг; рп — парциальное давление водяного пара в воздухе, Па; Рб — барометрическое давление, Па (при отсутствии данный показатель принимается равным Р6 = 101 325 Па);
где Fок — площадь светового проема; qок — удельный тепловой поток через остекление, Вт/м2; В — по-
г ' ' ок
правочный коэффициент влияния вида остекления.
Теплопритоки к поверхности льда от осветительных приборов можно найти с помощью приближенной зависимости:
d = d' ■ 103 = 622 ■
Ф- Рш
(( -Ф-РшУ
Ф = ^■100
(13)
(14)
где q — удельные тепловыделения от световых приборов, Вт/м2.
Теплопритоки от людей на льду и в прилегающих зонах устанавливают по формуле:
Ст.лоод = П ■ 4лоод, (10)
где п — количество спортсменов и судей; qлюд — удельные тепловыделения от людей, Вт.
где рн — парциальное давление насыщенного водяного пара (давление насыщения), Па; рп — парциальное давление водяного пара в воздухе, Па. При температуре больше 0 °С:
„241,2+t
Pi = 0,6112 ■ е2 При температуре от 0 °С и менее:
22,489'?
рн = 0,6112 ■ e 272,88+i,
(15)
н
1804
где рн — парциальное давление насыщенного водяного пара или давление насыщения, Па; t — температура воздуха, °С.
353 -1,32 ■ 10-3 ■рп
(17)
где рп — парциальное давление водяного пара в воздухе, Па; Т — температура воздуха, К;
св.в. = св + сп ■ й ■= 1,8 + 1,005 ■ й ■ 10-3, (18)
где св — средняя удельная теплоемкость сухого воздуха, кДж/(кгград); сп — средняя удельная теплоемкость водяного пара, кДж/(кгтрад); d — влагосо-держание воздуха, г/кг.
Температура точки росы определяется по формуле:
237,7 ■ Г 17,27 ■ ^ + 1пГ^-1Г
1237 + ^ 1юо ))\
17,27 - Г 17,27 ■ + 1пГ фв 11
1237 + tв 1100 ))\
^ = Р(РПов - Рв ),
(20)
и массоотдачи. Поскольку концентрация воды в воздухе у ледового покрытия незначительно изменяется по нормали к межфазной поверхности, для рассматриваемой задачи представляется допустимым применить случай малой интенсивности мас-сообмена, который характеризуется почти полным отсутствием влияния поперечного потока массы на течение смеси влажного воздуха. Подробно тройная аналогия для процессов тепломассообмена рассматривалась в работе [20], в которой отмечалось, что для рассматриваемой задачи вынужденной конвекции критериальные уравнения примут вид:
Ки = / (Яе,Рг); Ки Л = / (Яе,РГд),
(21)
(19)
где tв — температура воздуха, °С; фв — относительная влажность воздуха, %.
Как отмечалось в работах [15-18], в существующих методиках, основывающихся на графоаналитическом построении процессов обработки воздуха на М-диаграмме, расчет начинают с нанесения точек Н и Вл, параметры которых соответствуют параметрам наружного и внутреннего воздуха. Однако считается, что приток, подаваемый в направлении ледового поля сосредоточенными струями, охлаждается при постоянном влагосодержании. Исследования, приведенные в трудах [6, 14, 19], показывают, что это не соответствует действительности. В этой связи для учета влияния ассимиляции избытков влаги на процесс изменения параметров приточного воздуха в обслуживаемой зоне допустимо применять построение луча процесса с угловым коэффициентом, связывающим начальные и конечные значения параметров воздуха. Для этого необходимо определить тепловлагопоступления к внутреннему воздуху. Первая величина устанавливается по уравнению (2), для расчета второй составляющей в инженерной практике принято использовать аналогичную зависимость:
где Re, Рг и Ргд — соответственно критерий Рей-нольдса, критерий Прандтля и диффузионный критерий Прандтля.
Омывание поверхности ледового поля вентиляционным воздухом представляет собой частный случай омывания плоской пластины, для которого можно записать:
Ш = 0,66Яе1/2Рг1/3,
(22)
где Ми — средний по характерному размеру критерий Нуссельта.
Для записи критериального уравнения массо-обмена необходимо тепловой критерий Прандтля заменить на диффузионный:
Шв = 0,66Яе1/2 Рг,1/3 ,
(23)
В формулах (22) и (23) критерии подобия определяются следующим образом:
^ = К = ^;Яе = ^, рЛ X ув
(24)
где I — характерный размер (длина поверхности), м; D — коэффициент диффузии, м2/с; 1 — коэффициент теплопроводности газовой смеси, Вт/(мК); V — скорость газовой смеси, м/с; V — кинемати-
вв
ческий коэффициент вязкости газовой смеси, м2/с.
В рассматриваемой задаче под газовой смесью понимается смесь воздуха с частицами водяного пара, где воздух считается газом постоянного состава. В этом случае коэффициент диффузии можно определить с помощью зависимости [20]:
где в — коэффициент массоотдачи, кг/(ч-м2-Па); рпов — парциальное давление воздуха при температуре жидкости, Па; рв — парциальное давление воздуха в обслуживаемой зоне, Па; Е — площадь поверхности льда, м2.
Соотношение (20) аналогично закону Ньютона - Рихмана, следовательно, процессы конвективного тепло- и массообмена описываются аналогичными дифференциальными уравнениями и, согласно тройной аналогии, можно определить количественную связь между коэффициентами тепло-
D = Dr
/'уЛ
т
\уо у
Л. р'
(25)
где D0 — коэффициент диффузии при условиях Т0 = 273 К, р0 = 0,101 МПа (D0 = 0,216 • 104 м2/с); Т и р — соответственно температура и давление смеси для произвольных условий.
Таким образом, согласно работе [20], при обтекании плоской пластины воздушной смесью допустимо с достаточной точность считать, что коэффи-
< п
I*
На
Ч
о Г м з
о со
§ С/3
< -ь о СО
и -
Г I
<<3 о <
§ л
СП '
со со
КЗ со
0
1
СП СП о о
• л ® 8
л '
ов оо I г
(Л п
л у с о Ф Ж
Ы 10
о о
10 10
ы ы
1805
(О (О
N N
О О
N N
> (Л
с и to «о
во ^
о н
циенты тепло- и массообмена связаны между собой соотношением:
Р = ^,
ср
(26)
где ср — теплоемкость влажного воздуха, Дж/(кгК), если считать, что смесь влажного воздуха имеет постоянный состав, то ее теплоемкость можно определить по формуле:
- - Ср1Свп + Ср2Св,
(27)
<и <и
где свп и св — соответственно массовые концентрации водяного пара и воздуха, кДж/(кгК).
Таким образом, зная величину коэффициента массообмена, можно найти объем влагопоступления к воздуху над поверхностью льда.
Поскольку для рассматриваемых объектов рекомендуется предусматривать рециркуляцию воздуха, что снижает энергетические затраты при его обработке [9, 16-18], большой интерес вызывает задача определения параметров воздуха после смешивания. В настоящий момент существует три способа организации воздухообмена зоны ледового поля по схемам «сверху-вверх», «сверху-вниз» и смешанной схеме, подробно они описаны в трудах [7, 9, 15-17].
Применительно к данной задаче нас интересует вопрос размещения воздухозаборных устройств для рециркуляции, которая может осуществляться как из верхней, так и из нижней зон, а также совместно. От этого зависят не только параметры смеси, но и дальнейшие расходы холода и теплоты на обработку приточного воздуха. Наиболее трудоемкими вычисления будут для смешанной схемы воздухо-
обмена, это связано с необходимостью определения соотношения объемов рециркуляционного воздуха, отбираемого совместно из верхней и нижней зон. При этом представляется целесообразным стремиться привести влагосодержание воздуха после второй ступени рециркуляции к значению, равному или максимально близкому к влагосодержанию приточного воздуха. В этом случае при равенстве названых параметров не требуется дополнительное осушение или увлажнение воздуха перед его подогревом в калорифере. Если добиться равенства влагосодержания проблематично, следует выбрать такое соотношение рециркуляционных расходов, при котором разница влагосодержания будет минимальной.
РЕЗУЛЬТАТЫ ИССЛЕДОВАНИЯ
Для учета всех упомянутых факторов была разработана методика аналитического расчета параметров состояния воздуха при вентиляции и кондиционировании обслуживаемой зоны над ледовым полем, структурная блок-схема которой приведена на рис. 1.
На основании приведенной выше методики разработана программа расчета параметров состояния воздуха в зоне ледового поля. Реализация программной части выполнялась на высокоуровневом свободном объектно-ориентированном языке программирования Python. Данное приложение имеет модульную архитектуру, которая является сложной в плане разработки, но позволяет развивать программный продукт, добавляя новые функциональные возможности и корректируя уже реализованные, не изменяя архитектуры всего программного продукта, а затрагивая лишь его отдельные компоненты — независимые модули. Пример работы программы показан на рис. 2.
о ё
о о со со
о со гм
от от
.Е о cl"
^ с ю О
S «
о Е
fe ° СП ^ т- ^
от от
О
Э
(Л (9
s *
I. Определение притоков теплоты к поверхности льда / I. Determination of heat inflows to the ice surface
--------------------------------------------------------------
a
Рис. 1. Структурная блок-схема методики расчета параметров состояния воздуха: a — определение притоков теплоты к поверхности льда; b — определение параметров состояния воздуха; c — определение расходов теплоты и холода в центральном кондиционере Fig. 1. Structural block diagram of the methodology for calculation of air condition parameters: a — determination of heat inflows to the ice surface; b — determination of air condition parameters; c — determination of heat and cold consumption in the central air conditioner
1806
Разработка методики и программы расчета параметров воздуха при кондиционировании _ . вло . ол 0
С. 1802—1812
и вентиляции крытых катков
2.1. Влагопоступления от людей 2.1. Moisture intake from people
____________Л___________________________
2. Определение параметров состояния воздуха при кондиционировании зоны ледового поля 2. Determination of air condition parameters during conditioning of the ice field area
2.2. Минимальный расход наружного воздуха 2.2. Minimum outdoor air flow rate
2.3. Расчет термодинамических параметров влажного воздуха / 2.3. Calculation of thermodynamic parameters of humid air
2.3.1. Параметры наружного воздуха (исходные данные: P6, f или p"j Последовательность решения уравнений: (15) или (16) [в зависимости от температуры], (14), (12), (11), (17), (18), (19) 2.3.1. Outdoor air parameters (initial data: Pb, ts, f или pn j Sequence of solution of equations: (15) or (16) [depending on temperature], (14), (12), (11), (17), (18), (19)
2.3.2. Параметры воздуха в обслуживаемой зоне (исходные данные: Pg, f, ф) Последовательность решения уравнений: (15) или (16) [в зависимости от температуры], (14), (12), (11), (17), (18), (19) 2.3.2. Air parameters in the served area (initial data: P., t, ф") Sequence of solution of equations: (15) or (16) [depending on temperature], (14), (12), (11), (17), (18), (19)
2.3.3. Параметры приточного воздуха (исходные данные: L^, P6, dn = dB) Последовательность решения уравнений: t" = ts + [3,6 • (Qt.KOh - £т.люд)/L ' P" ' CP~\> далее уравнения (11), (15) или (16), (12), (14), (17), (18), (19) / 2.3.3. Supply air parameters (input data: Ln,Pb, da = da) The sequence of solving the equations: tn = tB + |^3.6 • Qhcon - Qhpeop~)/L" ' P" ' c"~], then equations (11), (15) or (16), (12), (14), (17), (18), (19)
2.3.4. Параметры удаляемого воздуха (исходные данные: tB, h, grad t, P6)
'W + а„н((н - p)cn
Последовательность решения уравнении: f = tE + gradt(h - 2), dy = dn + -
далее уравнения (11), (15) или (16), (12), (14), (17), (18), (19) 2.3.4. Removal air parameters (initial data: ta, h, grad t, P^)
The sequence of solving the equation: ty — ta + grad t[h — 2), dy — dn H--
further equations (11), (15) or (16), (12), (14), (17), (18), (19)
L ■pn
W + Oon (p" - p"a \/1сПр La ■ Pa
2.3.5. Параметры смеси наружного и рециркуляционного воздуха (исходные данные: Ц"", I", d", P j 2.3.5. Mixture parameters of outdoor and return air (initial data:Lm, Iu, d11, РЛ
Схема воздухообмена «сверху-вверх» / Air exchange scheme "from top to top"
+ Да / Yes
Последовательность решения уравнений:
_(( - ¡ГУ'+ЦТ I, _( - + ifd,
' in ' ' in
Sequence of solving equations
(in - ЦГ'У + ЦТ1, (Ln - цу + ifd,
' in ' ' in
(11), (15) или (16), (12), (14), (17), (18) и (19)
(11), (15) or (16), (12), (14), (17), (18), and (19)
Схема воздухообмена «сверху-ввниз» / Top-to-bottom air exchange scheme
± Да / Yes
Последовательность решения уравнений:
I' _
( - ЦГ 'У' + ЦТ I" ' ( - Ц? " У + LTd"
L,
■, (11), (15) или (16), (12), (14), (17), (18) и (19)
Sequence of solving equations:
I' _
(in - iry■+с I' ' _ (in - СУ+Cd■
in
in
■, (11), (15) or (16), (12), (14), (17), (18), and (19)
Расчетные зависимости (смешанная схема): на 1-й ступени рециркуляции Ь 2 задается интервал 1 м3/ч
I, _
((- ir2 y ц?У+LTI" d ((-1,2 y Ц?У+iTd"
■, (11), (15) или (16), (12), (14), (17), (18) и (19
(Ln - Lrl) ' " (I, - Lp2
Calculated dependencies (mixed scheme): At the 1st stage of recirculation L 2 is given the interval 1 m3/h
^ _ Qln - Lrl)- Lf )/a + LfIs d _((Q - Lrl)- Ef )da + Efr
(11), (15) or (16), (12), (14), (17), (18) and (19
2-я ступень рециркуляции
L , Iy + (L - L , Ул L ,dy + (L - L , Ул
-^ ' , da ---, (11), (15) или (16), (12), (14), (17), (18) и (19)
2nd recirculation stage
L21' + (( - Lp2 У _ Lp2d' + (L - Lp2 У
Ln
Ln
(11), (15) or (16), (12), (14), (17), (18) and (19)
d= dn или / or dc ^ dn
Да / Yes
1 ■ ■ 1
2.3.6. Параметры воздуха при его подаче в калорифер для нагрева — точка О (исходные данные: d°) 2.3.6. Air parameters at its supply to the heater for heating — point O (initial data: d°)
Рис. 1. Структурная блок-схема методики расчета параметров состояния воздуха: a — определение притоков теплоты к поверхности льда; b — определение параметров состояния воздуха; c — определение расходов теплоты и холода в центральном кондиционере (продолжение)
Fig. 1. Structural block diagram of the methodology for calculation of air condition parameters: a — determination of heat inflows to the ice surface; b — determination of air condition parameters; c — determination of heat and cold consumption in the central air conditioner (continuation)
1807
< П
is
G Г
o CO n CO
J to
u -
r i
n о
< 3 o
П)
СЯ '
n 2
< 0
<
r 6 t ((
go
• ) ® 8
л *
OS n I T
s у с о <D X
hi 10 о о 10 10 u w
Нет / No
Нет / No
Нет / No
(О (О
N N
О О
N N
¡г ш
U 3 > (Л С И
to «о
. г
СО щ
?!
Ф Ф
О ё
о о со со
I
о со cn
от от
.£ о
cl"
с
Ю о
о Е
fe ° СП ^ т- ^
£
от °
* А
d° = d"
Да / Yes|
I' = Iе , d' = dc2, t' = tc2, pl = pl2,
ф' =фс 2 p» = p»2, С = pi, t; = tp2.
dc(dc2) > dn
Да / Yes + '
Расчетные зависимости (дополнительные данные ф0 = 90 %) Calculated dependenci s (additional lata ф0 = 90 %)
d"= 622•
( - Г
• ioo;
pi
pi > 611,15 Па / Pa
Да / Yes \
Pl = 0,6112 • e241,2"', I' = 1,005 • t' + (2500 +1,8 • t') ■ do • 10-3, p» =
cO, = 1,8 +1,005 • d' • 10-3, t' = 237,7 • [17,27 • f + li
.... . . ' 1 237 +t'
353 -1,32 • 10-3 • p°
/17,27 -[1727^ + JV
I 237 + Г 1100
22,489-•t" Pl = 0,6112 • e272№t', I" = 1,005 • t" + (2500 + 1,8 • t") d" • 10-3, pB = 353 - 1,32 • 10-3 • Pl,
T"
<B = 1,8 + 1,005 • d" • 10- , f = 237,7 •i17,27 • ^ + ln[ p [ 237 + t" 100))/ -f17,27 • ln[ [ 237 + t" 100))
Расчетные зависимости (дополнительные данные lo - Iе): (11), (12), (15) или (16), (14), (17), (18) и (19) Calculated dependencies (...): (11), (12), (15) or (16), (14), (17), (18) and (19)
2.4. Расчет напоров температуры и влагосодержания / 2.4. Calculation of temperature and humidity content
AtT = t - to, Atx = te - t o, Ad = dc - d0
3. Определение расхода теплоты и холода для обработки воздуха в центральном кондиционере 3. Determination of heat and cooling consumption for air handling in the central air conditioning unit
С
Рис. 1. Структурная блок-схема методики расчета параметров состояния воздуха: a — определение притоков теплоты к поверхности льда; b — определение параметров состояния воздуха; c — определение расходов теплоты и холода в центральном кондиционере (окончание)
Fig. 1. Structural block diagram of the methodology for calculation of air condition parameters: a — determination of heat inflows to the ice surface; b — determination of air condition parameters; c — determination of heat and cold consumption in the central air conditioner (ending)
Разработанная программа позволяет определять не только параметры состояния влажного воздуха при организации кондиционирования зоны ледового поля, но и рассчитывать затраты холода и теплоты на его обработку при различных схемах организации воздухообмена. Кроме того, для смешанной схемы осуществляется итерационный расчет оптимального соотношения объемов рециркуляционного воздуха, отбираемого из различных зон с отличающимися значениями температуры и влагосодержания.
Анализ приведенных на рис. 2, с результатов расчета показывает, что при воздухообмене в холодный
период года по смешанной схеме энергетические затраты на обработку приточного воздуха в центральном кондиционере значительно уменьшаются. Например, может не потребоваться дополнительное охлаждение подаваемого воздуха для его осушения, которого можно добиться путем выбора требуемого соотношения объемов рециркуляционного воздуха из верхней и рабочей зоны. Однако в теплый период года применение данной схемы может увеличивать затраты холода и теплоты по сравнению со схемой «сверху-вниз». По этой причине при круглогодичной эксплуатации помещений с искусственным льдом распределение
Нет / No
Нет / No
Нет / No
1808
< П
I*
4
G Г
5 3
b
Рис. 2. Пример работы программы: a — блок ввода исходных данных для теплового расчета; b — блок ввода исходных данных для расчета параметров состояния; c — блок вывода данных расчета
Fig. 2. Example of programme operation: a — input block of initial data for thermal calculation; b — input block of initial data for calculation of state parameters; c — calculation data output block
0 CO
§ со
1 < < -b J to
U -
r I
П о
< 3 o
oi
СЯ '
U S
§ 2
< 0 <
r 6
t °
• ) ® 8
л *
00 m
1 T
s У с о ® Ж
Ы 10
о о
10 10
U W
a
1809
Вывод данных расчета
(О (О
N N
О О
N N
К ш U 3 > (Л С И 2 "1 to «о
. т-
во щ
?!
Ф Ф
О ig
о о со со
I
о со сч
от от
.£ о
CL^
с
Ю о
S g
о Е
fe ° СП ^ т- ^
s
от °
* А
Точки ü о 1 и M ■а. и ь- _ И k 1Ч 5 О 9- и -5- « К а * г Ts Ь* S" TÈ Ь "г ь £ в" в а <У Qt, Вт Qx, Вт G, кг ч
Схема «сверху-вверх»
H -7.5 1.32 1.78 -3,15 89 -8.7 0.29
П 13.4 1,23 2,65 20,21 28 -4,2 0,430
вл 6,0 1.27 2.65 12.7 46 -4.2 0,430 2000 15000 92880 93790 63200
У 9.5 1.25 3.10 17.4 42 -2.4 0.503
С' 7.5 1.26 2,90 14.9 46 -2,9 0.480
о -2.9 1.30 2.65 3.59 90 -4.2 0.429
Схема «сверху-вниз»
H -7.5 1.32 1.78 -3,15 89 -8.7 0,29
П 13.4 1.23 2.65 20.21 28 -4.2 0.430
ВЛ 6,0 1.27 2.65 12.7 46 -4.2 0.430 2000 15000 92880 53400 2,04
У 6,0 1,27 2,65 12,7 46 -4,2 0,430
С' 4.4 1.27 2.55 10.8 49 -4.7 0.410
о 4.1 1.27 2.65 10.8 52 -4.2 0.43
Смешанная схема
H -7.5 1.32 1.78 -3,15 89 -8.7 0,29
П 13,4 1,23 2,65 20,21 28 -4,2 0,430
ВЛ 6,0 1.27 2.65 12.7 46 -4.2 0.430
У1 6,0 1.27 2.65 12.7 46 -4.2 0.430 2000 10500 4500 92880 46580
У2 1,25 3,10 17,4 42 -2,4 0,503
С 3,8 1,27 2,51 10,15 50 -4,8 0,410
С" 5.3 1.27 2.65 17.4 42 -4.1 0.430
о 5,3 1,27 2,65 17,4 42 -4,1 0,430
Рис. 2. Пример работы программы: a — блок ввода исходных данных для теплового расчета; b — блок ввода исходных данных для расчета параметров состояния; c — блок вывода данных расчета (окончание)
Fig. 2. Example of programme operation: a — input block of initial data for thermal calculation; b — input block of initial data for calculation of state parameters; c — calculation data output block (ending)
воздуха в теплый и холодный периоды года необходимо осуществлять по различным схемам.
ЗАКЛЮЧЕНИЕ И ОБСУЖДЕНИЕ
Разработанная аналитическая методика и программа расчета параметров состояния влажного воздуха дает возможность повысить точность расчета температуры, влагосодержания, парциального давления и энтальпии при проектировании систем вентиляции и кондиционирования помещений с искусственным льдом. Это способствует сокращению избыточных энергетических затрат при обработке
воздуха в центральном кондиционере. Отличительной особенностью предложенной методики является применение итерационного метода расчета при выборе оптимального соотношения расходов воздуха на первой и второй ступенях рециркуляции, что позволяет исключить при определенных условиях дополнительное охлаждение притока в холодный период года. Кроме того, использование тройной аналогии при рассмотрении конвективного тепломассообмена внутреннего воздуха с поверхностью льда позволило повысить точность расчета влагопо-ступлений к удаляемому воздуху за счет введения поправки к влагосодержанию.
СПИСОК ИСТОЧНИКОВ
1. Levy J.I., Lee K., Yanagisawa Y., Hutchinson P., Spengler J.D. Determinants of nitrogen dioxide concentrations in indoor ice skating rinks // American Public Health Association. 1998. Vol. 88. Issue 12. Pp. 1781-1786. DOI: 10.2105/ajph.88.12.1781
2. Yang C., Demokritou P., Chen Q., Spengler J. Ventilation and air quality in indoor ice skating arenas // ASHRAE Transactions. 2000. Vol. 106. Pp. 4405-4414.
3. Ivanov Y., Novikov N. Dual-channel digital control of energy consumption and air supply in micro-
climate systems of livestock premises // IOP Conference Series: Earth and Environmental Science. 2021. Vol. 937. Issue 3. P. 032002. DOI: 10.1088/17551315/937/3/032002
4. Русаков С.В. К расчету тепловых и влаж-ностных нагрузок ледовых катков. Нагрузка от радиационного переноса теплоты // Научный журнал НИУ ИТМО. Серия: Холодильная техника и кондиционирование. 2014. № 1. С. 10. EDN TOBGVX.
1810
c
5. KhaledM., Ali S., Jaber H., Faraj J., Murr R., Lemenand T. Heating/cooling fresh air using hot/cold exhaust air of heating, ventilating, and air conditioning systems // Energies. 2022. Vol. 15. Issue 5. P. 1877. DOI: 10.3390/en15051877
6. Анисимов С.М., Денисихина Д.М., Полу-шкин В.И. Решение задачи турбулентного переноса импульса, тепла, примеси в объеме «чаши» ледовой арены // Вестник гражданских инженеров. 2012. № 5 (34). С. 149-155. EDN PZVPRB.
7. Вишневский Е.П. Вентиляция и качество воздуха в крытых ледовых аренах // СОК. 2008. № 10. C. 34-39.
8. Денисихина Д.М. Исследование различных схем воздухораспределения ледовых арен // Известия высших учебных заведений. Строительство. 2015. № 5 (677). С. 38-48. EDN UBLRPF.
9. Кокорин О.Я. Инженерные системы помещений с искусственным льдом или снегом. М. : КУРС, 2022. 240 с.
10. Кокорин О.Я., Товарас Н.В. Снижение энергопотребления в системах кондиционирования воздуха для помещений с искусственными ледовыми полями // Холодильная техника. 2012. № 12. С. 12-15. EDN SEKLYR.
11. Дячек П.И., Ливанский Д.Г. Температурная неоднородность поверхности ледового поля крытых искусственных катков // Энергетика. Известия высших учебных заведений и энергетических объединений СНГ. 2010. № 2. С. 41-47.
12. Тарасова Е.В., Штым А.С. Изменение термодинамических параметров охлажденного воздуха при прямом контакте с естественным источником холода // Вестник гражданских инженеров. 2012. № 6 (35). С. 107-112. EDN QZATVX.
Поступила в редакцию 12 мая 2023 г. Принята в доработанном виде 15 мая 2023 г. Одобрена для публикации 2 сентября 2023 г.
13. Каутский А.А. Регулирование и контроль температуры льда искусственного катка // Холодильная техника. 1975. № 11. С. 51-52.
14. Пухкал В.А., Юстус Д.А. Воздухораспреде-ление в помещениях ледовых площадок со зрителями // Строительство уникальных зданий и сооружений. 2015. № 12 (39). С. 7-31. EDN VJPKGD.
15. Chuykin S.V., Glazkov S.S. Determination of the heat return coefficient of the ice surface for the mixed air distribution scheme // Scientific Herald of the Voronezh State University of Architecture and Civil Engineering. Construction and Architecture. 2013. Vol. 3 (19). Pp. 29-38. EDN QLHSSJ.
16. Chuykin S.V. Problems of developing methods of designing microclimate systems for premises with artificial // Russian Journal of Building Construction and Architecture. 2022. Vol. 4 (56). Pp. 18-28. DOI: 10.36622/VSTU.2022.56.4.002. EDN QJMDOJ.
17. Melkumov V.N., Chuykin S.V. A scheme and method of calculation for ventilation and air conditioning systems of ice arenas // Journal of Technology. 2017. Vol. 32. Issue 2. Pp. 139-146.
18. Melkumov V.N., Chuykin S.V. Organization of air distribution of covered multipurpose ice rinks // Scientific herald of the Voronezh state university of architecture and civil engineering. Construction and architecture. 2013. Vol. 3 (19). Pp. 17-28. EDN QLHSRZ.
19. Старкова Л.Г., Анисимова Е.Ю., Сорокина Н.Г., Дегтярь С.Л. Анализ схем вентиляции крытого катка с помощью числовой модели воздушных потоков // Вестник Южно-Уральского государственного университета. Серия: Строительство и архитектура. 2021. Т. 21. № 3. С. 42-50. DOI: 10.14529/ build210306
20. Цветков Ф.Ф. Тепломассообмен. М. : Изд-во МЭИ, 2005. 550 с.
Об авторе: Сергей Владимирович Чуйкин — кандидат технических наук, доцент кафедры строительной механики; Воронежский государственный технический университет (ВГТУ); 394006, г. Воронеж, ул. 20-летия Октября, д. 84; SPIN-код: 4528-7170, Scopus: 57216439704, ResearcherID: R-3144-2016, ORCID: 0000-0001-6040-323X; ser.chu@mail.ru.
REFERENCES
1. Levy J.I., Lee K., Yanagisawa Y., Hutchinson P., Spengler J.D. Determinants of nitrogen dioxide concentrations in indoor ice skating rinks. American Public Health Association. 1998; 88(12):1781-1786. DOI: 10.2105/ajph.88.12.1781
2. Yang C., Demokritou P., Chen Q., Spengler J. Ventilation and air quality in indoor ice skating arenas. ASHRAE Transactions. 2000; 106:4405-4414.
3. Ivanov Y., Novikov N. Dual-channel digital control of energy consumption and air supply in microclimate systems of livestock premises. IOP Conference Series: Earth and Environmental Science. 2021; 937(3):032002. DOI: 10.1088/1755-1315/937/3/032002
4. Rusakov S. To calculation of thermal and moist loadings ice skating rinks loading from radiant heat transfer. Scientific Journal of St. Petersburg State University of Low
< П iiï G Г
S 2
0 со n СО
1 <
< -ь J со
U -
r I
n о
< 3 o
n
со со
l\J со
0
1
cd СП о о
cn
• ) ® 8
л '
00 п
1 Т
s У с о (D *
Ы 10
о о
10 10
U W
1811
C.B. HyOKUH
W (0
N N
o o
N N
* o
U 3 > in E M
to co
. T-
AO (U
i!
<D <u
O 3
o o CD CD
I
o co CN
iO
Temperature and Food Technologies: Refrigeration and Air Conditioning. 2014; 1:10. EDN TOBGVX. (rus.).
5. Khaled M., Ali S., Jaber H., Faraj J., Murr R., Lemenand T. Heating/cooling fresh air using hot/cold exhaust air of heating, ventilating, and air conditioning systems. Energies. 2022; 15(5):1877. DOI: 10.3390/en15051877
6. Anisimov S.M., Denisikhina D.M., Polush-kin V.I. The solution of the problem of turbulent transfer of momentum, heat and scalar in the volume of the ice arena "bowl". Bulletin of Civil Engineers. 2012; 5(34):149-155. EDN PZVPRB. (rus.).
7. Vishnevsky E.P. Ventilation and air quality in indoor ice arenas. Plumbing, Heating, Air-Conditioning. 2008; 10:34-39. (rus.).
8. Denisikhina D.M. Research of different air distribution schemes in ice rink arena. News of higher educational institutions. Construction. 2015; 5(677):38-48. EDN UBLRPF. (rus.).
9. Kokorin O.Ya. Engineering systems of rooms with artificial ice and snow. Moscow, KURS Publ., 2022; 240. (rus.).
10. Kokorin O.Ya., Tovaras N.V. Reduction of energy consumption in air conditioning systems for rooms with artificial ice fields. Refrigerating equipment. 2012; 12:12-15. EDN SEKLYR. (rus.).
11. Diachek P.I., Livansky D.G. Temperature non-uniformity of ice-field surface at covered skating-rinks. ENERGETIKA. Proceedings of CIS higher education institutions and power engineering associations. 2010; 2:41-47. (rus.).
12. Tarasova E.V., Shtym A.S. Alteration of thermodynamic parameters of the cooled air in direct contact with a natural source of cold. Bulletin of Civil Engineers. 2012; 6:107-112. EDN QZATVX. (rus.).
Received May 12, 2023.
Adopted in revised form on May 15, 2023.
Approved for publication on September 2, 2023.
BioNoiEs: Sergey V. Chuykin — Candidate of Technical Sciences, Associate Professor of the Department of Construction Mechanics; Voronezh State Technical University (VSTU); 84 20-letiya Oktyabrya, Voronezh, 394006, Russian Federation; SPIN-code: 4528-7170, Scopus: 57216439704, ResearcherlD: R-3144-2016, ORCID: 0000-0001-6040-323X; ser.chu@mail.ru.
13. Kautskij A.A. Regulation and control of the ice temperature of an artificial ice rink. Refrigerating Equipment. 1975; 11:51-52. (rus.).
14. Pukhal V.A., Yustus D.A. Air distribution in spectator ice arenas. Construction of Unique Buildings and Structures. 2015; 12(39):7-31. EDN VJPKGD. (rus.).
15. Chuykin S.V., Glazkov S.S. Determination of the heat return coefficient of the ice surface for the mixed air distribution scheme. Scientific Herald of the Voronezh State University of Architecture and Civil Engineering. Construction and Architecture. 2013; 3(19):29-38. EDN QLHSSJ.
16. Chuykin S.V. Problems of developing methods of designing microclimate systems for premises with artificial. Russian Journal of Building Construction and Architecture. 2022; 4(56):18-28. DOI: 10.36622/ VSTU.2022.56.4.002. EDN QJMDOJ.
17. Melkumov V.N., Chuykin S.V. A scheme and method of calculation for ventilation and air conditioning systems of ice arenas. Journal of Technology. 2017; 32(2):139-146.
18. Melkumov V.N., Chuykin S.V. Organization of air distribution of covered multipurpose ice rinks. Scientific herald of the Voronezh state university of architecture and civil engineering. Construction and architecture. 2013; 3(19):17-28. EDN QLHSRZ.
19. Starkova L.G., Anisimova E.Yu., Soroki-na N.G., Deg-tiar S.L. Analysis of ventilation schemes for indoor ice skating rinks using a numerical air flow model. Bulletin of the South Ural State University. Ser. Construction Engineering and Architecture. 2021. 21(3):42-50. DOI: 10.14529/build210306 (rus.).
20. Tsvetkov F.F. Heat and mass exchange. Moscow, Publishing House of MEI, 2005; 550. (rus.).
.E o cl"
• c Ln o
s«
o E
fe ° CD ^
T- ^
■s
O (0
1812
