Mathematical modeling of heat and mass transfer processes of air masses in premises with an ice arena

Ivan A. Burkov; Бурков Иван Александрович; Michail  A. Kolosov; Колосов Михаил Анатольевич; Anton A. Zharov; Жаров Антон Андреевич; Leonid М. Kolishkin; Колышкин Леонид Михайлович; Nikolay  A. Andreev; Андреев Николай Алексеевич

doi:10.17816/RF322770

Математическое моделирование процессов тепломассообмена воздушных масс в помещениях с ледовой ареной

Авторы: Бурков И.А.¹, Колосов М.А.¹, Жаров А.А.¹, Колышкин Л.М.¹, Андреев Н.А.¹
Учреждения:
1. Московский государственный технический университет им. Н.Э. Баумана (научно исследовательский университет)
Выпуск: Том 110, № 2 (2021)
Страницы: 77-84
Раздел: Оригинальные исследования
URL: https://freezetech.ru/0023-124X/article/view/322770
DOI: https://doi.org/10.17816/RF322770
ID: 322770

Цитировать

Полный текст

Аннотация
Полный текст
Об авторах
Список литературы
Дополнительные файлы
Статистика

Аннотация

Обоснование: Системы вентиляции и кондиционирования ледовых арен должны обеспечивать комфортное пребывание людей в помещении и на льду, а также защиту поверхности льда от конденсации влаги. Использование метода конечных элементов при проектировании систем вентиляции позволяет на этапе расчета определить температурно-влажностный режим во всем объеме вентилируемого помещения и, при необходимости, исправить параметры вентиляции с минимальными затратами на этапе проектирования.

Цель – описание методики численного моделирования тепломассообмена в помещении с учетом влажного воздуха и лучистого теплообмена на примере небольшой тренировочной арены.

Методы: для моделирования было использовано программное обеспечение ANSYS CFX, использующее при расчетах метод конечных объемов. Объектом исследования являлся цифровой двойник малой тренировочной ледовой арены, позволяющий учитывать теплофизические процессы, происходящие в помещении. Исследование проводилось для установившегося (стационарного) режима тепломассообмена. Оценка результатов моделирования проводилась по контурам температур и влажности воздуха в наиболее характерных секущих плоскостях. Ввиду отсутствия возможности прямого определения влажности воздуха средствами программного комплекса, было приведено описание способа определения влажности по эмпирическому уравнению.

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

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

Ключевые слова

ледовая арена, температурные поля, влажный воздух, HVAC, ANSYS CFX

Полный текст

ВВЕДЕНИЕ

В мире и России сегодня проектируется и возводится множество спортивных залов с ледовыми аренами различных типов и размеров [1-3]. Системы кондиционирования и вентиляции (Heating Venting Air Conditioning, HVAC) подобных залов имеют много особенностей конструкции и индивидуальных условий эксплуатации, среди них: наличие искусственного ледового поля; жёсткие требования по параметрам льда и воздуха в зоне пребывания спортсменов и зрителей; сложность регулирования влажности в подобных помещениях, что часто приводит к образованию конденсата на внутренних поверхностях и т.п. Всё это ставит множество вопросов перед проектировщиками подобных объектов, ответы на которые можно получить с помощью моделирования. Особое внимание необходимо уделять защите ледовой арены от конденсации влаги, скапливающейся под потолком и над полем. Повышение влажности воздуха под потолком обусловлено несколькими факторами: лучистым теплообменом между потолком и льдом, внешним охлаждением потолка, некорректной работой вентиляционной системы. Практика показывает, что подобные проблемы, наносящие вред качеству ледовой поверхности, достаточно часто возникают на новых ледовых аренах, и это требует дорогой модернизации проекта. Кроме того, необходимо выполнить строгое зонирование пространства ледовых арен: трибуны для зрителей должны находиться в теплой зоне, где температура воздуха поддерживается в диапазоне +19…+22°С, а спортсмены должны находиться в зоне прохладного воздуха с температурой около +10°С (значение данной температуры зависит от вида спорта, которым занимаются спортсмены). Очевидно, что в любом случае недопустимо даже кратковременное образование тумана под потолком и у поверхности льда. Оптимальными величинами относительной влажности около поверхности льда являются 70…80%.

Пакеты вычислительной гидродинамики (CFD) в сочетании с мощными вычислительными машинами сегодня позволяют с высокой точностью проанализировать и изучать особенности работы систем вентиляции в просторных залах [4–14]. Использование пакетов CFD позволит избежать ошибок в конструкции ледовых арен на стадии проектирования, позволит сократить стоимость эксплуатации арены и повысит комфорт людей. Вместе с тем, особенности залов с ледовыми аренами накладывают ограничения на возможности моделирования, и поэтому сегодня нет устоявшихся подходов к расчёту и проектированию подобных систем. В связи с этим, проведение отработки и верификации таких программ целесообразно проводить на основе реализованных проектов.

ПОСТАНОВКА ЗАДАЧИ

На примере небольшой проектируемой ледовой тренировочной арены (до 150 зрителей и 100 спортсменов во время групповой тренировки) необходимо провести расчет тепломассобменных процессов в помещении, определить температурные и скоростные поля течения воздуха, а также его влажность. В рассматриваемой арене имеют место следующие процессы: воздушный тепломассообмен в помещении от приточно-вытяжной вентиляции с заданным объемным расходом воздуха, тепловыделение от ламп освещения, людей и конвекторов, охлаждение воздуха ледовой ареной, лучистый теплообмен между потолком и ледовой поверхностью, влаговыделение от людей и регулярной заливки ледовой поверхности. Предполагается, что ледовая арена является не отдельно стоящим зданием, а расположена внутри торгового центра, поэтому не учитывается теплоприток от стен и потолка внутрь помещения.

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ

Математическое моделирование осуществлялось в ANSYS 18.1 CFX. Данное приложение позволяет рассчитывать движение жидкостей и газов в пространствах сложной формы с учетом турбулентности и гравитации [15–20]. Расчет проводился в стационарном режиме, процессы, возникающие при первичной наморозке льда, не рассматривались.

Расчетная геометрическая модель показана на рис. 1. Она состоит из одного объемного тела – текучей среды, который состоит из двух веществ – сухого воздуха и водяного пара. На рис. 1 показаны места размещения основных элементов технологических систем, участвующих в формировании параметров воздуха в зале.

Рис. 1. Общий вид геометрической модели зала ледовой арены. 1 – сопла приточной вентиляции; 2 – вытяжные отверстия, 3 – воздухонагреватели; 4 – батареи центрального отопления; 5 – осветительные приборы; 6 – игроки; 7 – поверхность льда; 8 – зрители.

Fig. 1. General layout of the ice arena geometric model. 1 – nozzles of supply ventilation; 2 – exhaust openings, 3 – air heaters; 4 – central heating batteries; 5 – lighting; 6 – players; 7 – ice surface; 8 – spectators.

Сетка конечных элементов состояла из 19.8 млн элементов в форме тетраэдра с размерами от 5 мм до 500 мм. Для более точного моделирования процесса теплоотдачи радиаторов, ламп, зрителей и хоккеистов, в областях воздуха, граничащих с данными источниками теплоты, были использованы призматические пограничные элементы. Общий вид сетки конечных элементов в области зрительских трибун показан на рис. 2.

Рис. 2. Сетка конечных элементов в зоне зрительских трибун. Показаны измельчение элементов в форме тетраэдра в зоне граничных условий и пограничные слои призматической формы, непосредственно примыкающие к граничному условию.

Fig. 2. Finite elements grid in the area of the viewing stands. Shredding of elements in the form of a tetrahedron in the zone of boundary conditions and the boundary layers of a prismatic shape directly adjacent to the boundary condition are shown.

Граничные условия, описывающие происходящие в зале ледовой арены процессы, приведены в таблице 1. Данные граничные условия характерны для режима работы в теплое время года, этим объясняется отсутствие теплопритоков от нагревателей в представленном расчете, т.к. в данном режиме работы они не используются. Летний режим работы был выбран для расчета как более нагруженный с точки зрения влаговыделения в помещении. Количество зрителей и хоккеистов выбрано максимально допустимым по проекту арены при их одновременном нахождении в помещении.

Таблица 1. Граничные условия, использованные в расчете. Для зрителей и спортсменов граничные условия приведены для одного человека, общее количество зрителей – 150 человек, спортсменов – 100 человек

Table 1. Boundary conditions used in the calculation. For spectators and athletes, the boundary conditions are given for one person, the total number of spectators is 150 people, athletes – 100 people

Граничное условие	Температура / Тепловыделение	Влаговыделение	Расход воздуха
Ледовая арена	–6°С	25 кг/ч
Освещение	60 Вт/м²
Зрители	132 Вт/чел.	40 г/ч/чел.
Хоккеисты	191 Вт/чел.	135 г/ч/чел.
Вентиляция	20°С		16 000 м³/час
Лучистый теплоприток от потолка ко льду	60 Вт/м²

В расчете была использована модель турбулентности k-epsilon, тип остатка – RMS, величина остатка решателя – 10^-5. Моделирование осуществлялось с распараллеливанием расчета на 8 ядер.

При расчете решатель ANSYS Solver итерационным методом решает уравнения неразрывности, сохранения энергии, турбулентности, теплопроводности. Относительная влажность воздуха определяется отношением давления водяного пара к равновесному давлению насыщенных паров при данной температуре. При этом, равновесное давление определяется по эмпирическому уравнению Тетенса в зависимости от температуры.

РЕЗУЛЬТАТЫ МОДЕЛИРОВАНИЯ

Был проведен стационарный расчет с граничными условиями, приведенными в таблице 1. В результате моделирования были получены распределения температур, относительной влажности и скорости воздуха во всем объеме зала ледовой арены. На рисунке 3 показаны данные параметры в среднем вертикальном сечении арены.

Как видно из рис. 3a, температура воздуха в области пребывания зрителей находится в комфортных для людей пределах – не ниже +19°С. От зрительских трибун к потолку наблюдается восходящий поток теплого воздуха. Наибольшим источником влаги в помещении является ледовая поверхность (рис. 3b). При нормальной эксплуатации ледовая поверхность каждый час заливается теплой водой (около +60°С), часть из которой испаряется и формирует приток влаги в помещении. Влаговыделение от зрителей малозаметно, т.к. происходит непосредственно под приточной вентиляцией, подающей относительно сухой воздух. Также из поля относительной влажности видно, что под потолком ледовой арены отсутствует зона повышенной (более 80%) влажности.

Рис. 3c показывает траекторию течения приточного воздуха. Как видно из данного рисунка, в области пребывания хоккеистов имеется постоянное движение воздуха с незначительной (до 0.2 м/с) скоростью, обеспечивающий комфортное пребывание людей в данной зоне.

Рис. 3. Результаты моделирования ледовой арены в среднем сечении арены. a – поле температур; b – поле относительной влажности; c – линии тока из приточного сопла, расположенного с противоположной стороны от зрителей.

Fig. 3. The results of modeling in the middle section of the ice arena. a – the temperature field; b – the field of relative humidity; c – flow lines from the supply nozzle located on the opposite side of the spectators.

ВЫВОДЫ

В результате моделирования параметров зала с ледовой ареной были получены поля температур, скоростей и относительной влажности воздуха во всем объеме. Расчетным путем было установлено, что выбранные характеристики работы системы вентиляции обеспечивают требуемые эксплуатационные параметры воздуха. Подтверждено, что при данных условиях достигается обеспечение комфортных параметров воздуха в зонах пребывания людей. Предлагаемая методика моделирования позволяет численно определить основные особенности процессов тепломассообмена в помещениях с нехарактерными значениями температур и влажности.

Учет лучистого теплообмена необходим для более точного моделирования распределения параметров в помещении. В данной работе моделирование лучистого теплообмена от потолка ко льду было осуществлено путем задания теплового потока, полученного экспериментальным путем. При этом, для различных ледовых арен он может отличаться в зависимости от размеров арены, условий эксплуатации и внешних факторов. Для более точного учета лучистого теплообмена возможно использовать метод Монте-Карло. Данный метод позволит рассчитывать лучистый теплообмен между всеми поверхностями зала ледовой арены и повысит точность моделирования параметров воздуха.

ДОПОЛНИТЕЛЬНО

Вклад авторов. Все авторы внесли существенный вклад в разработку концепции и подготовку статьи, прочли и одобрили финальную версию перед публикацией.

Конфликт интересов. Авторы заявляют об отсутствии конфликта интересов, связанного с подготовкой и публикацией статьи.

Источник финансирования. Авторы заявляют об отсутствии внешнего финансирования при проведении исследования и подготовке публикации.

ADDITIONAL INFORMATION

Authors’ contributions. All authors made a substantial contribution to the conceptual development and preparation of this article and read and approved the final version before publication.

Competing interests. The authors declare that they have no competing interests.

Funding source. This study was not supported by external sources of funding.