Введение. Энергосбережение и защита окружающей среды от вредных выбросов является актуальной проблемой. Создание и поддержание требуемых параметров микроклимата обеспечивают инженерные системы, для работы которых требуются источники теплоснабжения, оказывающие существенное влияние на жизненный цикл строительного объекта, а также  являются наиболее затратной статьей коммунальных расходов Перспективным направлением энергосбережения является использование возобновляемых и нетрадиционных источников энергии, к которым относятся тепловая энергия недр земли, воздуха, воды, а также тепловые выбросы вытяжных  вентиляционных систем, технологических процессов промышленных предприятий [1–3], для утилизации  которых  применяют тепловые насосы.В зависимости от источника энергии различают геотермальные, водяные и воздушные теплонасосные установки. Системы с использованием теплонасосных установок имеют меньшие эксплуатационные расходы и выделяют меньше вредных веществ в атмосферу, по сравнению с традиционными источниками теплоснабжения [4–7]. Следовательно, разработка энергоэффективных конструкций теплонасосных установок, оптимизация их работы на различных этапах функционирования инженерных систем, является актуальной задачей управления жизненным циклом объектов строительства.Наиболее перспективными с точки зрения эксплуатации и изготовлении являются воздушные теплонасосные установки по сравнению с геотермальными [8–11]. Однако несмотря на преимущество воздушных тепловых насосов, их применение ограничено в регионах с холодным климатом, в том числе в России, в связи с образованием на наружной поверхности испарителя инея и ледяной корки. Это приводит к повышению сопротивления теплопередачи и снижению коэффициента преобразования теплоты. Существующие способы удаления наледи требуют дополнительных затрат энергии, что снижает эффективность работы установки. Наиболее эффективным является метод разрушения наледи с помощью механических колебаний [12, 13].В рамках программы Приоритет-2030, нами разработан испарительный блок воздушного теплового насоса, включающий теплообменник оригинальной формы, магнитострикционный излучатель для удаления наледи, разработан способ крепления сердечника излучателя к поверхности теплообменника, а также изготовлена теплонасосная установка и проведены ее испытания [14]. Целью настоящей работы является повышение теплопроизводительности и коэффициента трансформации теплоты воздушного теплового насоса на основании совершенствования системы удаления наледи и улучшения гидродинамических условий проведения тепло-массообменных процессов в испарительном блоке установки с системой MOVEBIT. Материалы и методы. Исследования проводили на экспериментальном стенде воздушного теплового насоса, включающего модернизированный теплообменник испарительного блока с системой удаления наледи MOVEBIT, представленный на рис. 1.  Рис. 1. Экспериментальный стенд с воздушным тепловым насосом, оснащенный системой MOVEBIT Холодопровод испарительного блока представлен в виде стержня диаметром 80 мм, для увеличения поверхности теплообмена имеет 12 пластин высотой 28 мм. Стержень теплообменника изготовлен из сплава алюминия методом экструзии. Поперечное сечение стержня составляет 14,24 см2 (0,001424 м2), длина стержня 1,5 м. Стержень теплообменника находится в трубе длиной 1,8 м с утеплителем. Внутренний диаметр трубы 100 мм. Соответственно площадь живого сечения для прохода воздуха составляет0,006426 м2. площадь одного стержня теплообменника составляет 0,984 м2. В дальнейших расчетах принимаем площадь теплообменной поверхности одного стержня, равной 1 м2. Испарительный блок лабораторной установки содержит 4 теплообменника типа «труба в трубе» площадью теплообмена 4 м2 соответственно (рис. 2). Рис. 2. Наружный блок теплонасосной установки с трубами и утеплителем ( один стержень показан без трубы) В испарительном блоке воздушного теплового насоса происходят процессы охлаждения воздуха и массообменные процессы с конденсацией паров влаги на поверхности теплообменника с последующим образованием инея и ледяной корки. Интенсивность процессов зависит от теплофизических свойств теплоносителя, гидродинамических условий обтекания теплообменной поверхности и формы теплообменника.  Для описания теплообменных процессов в модернизированной конструкции испарительного блока наиболее пригодной является критериальное уравнение Нуссельта (Nu) для теплообмена в кольцевом канале при критерии Прандля (Pr) в диапазоне 0,7 -100 и соотношением диаметров 1,2<(d2/d1)<1,4 [15]: ,        (1)здесь d2 – внутренний диаметр наружной трубы, м; d1 – эквивалентный диаметр внутреннего стержня, м; Reэкв – критерий Рейнольдса, рассчитанный по эквивалентному диаметру канала.Отношение температур среды и поверхности учитывается температурным коэффициентом:   ,               (2)где Tc – температура наружного воздуха, подаваемого в межтрубное пространство, К; Tж – температура поверхности внутреннего стержня, равная температуре фреонового хладагента К.                 Преимущество этого уравнения заключается в учете отношения диаметров наружной d2 и внутренней d1 трубы, температур, режима движения воздуха в межтрубном пространстве, оцененного по значениям критерия Рейнольдса (Re), рассчитанного по значениям эквивалентного диаметра межтрубного канала. Свободный коэффициент А в уравнении (1) можно рассматривать как критерий, оценивающий влияние структуры поверхности трубы на интенсивность тепломассообменных процессов.Основная часть. Исследования проводили при расходах наружного воздуха, подаваемого в испарительный блок от 400 м3/ч, до 1200 м3/ч, при этом скорость воздуха в межтрубном пространстве менялась от 4,2 м/с до 13 м/с. Исследования проводили при температурах среды от-26 °С до +10 °С.Обработка экспериментальных данных позволила установить численное значение коэффициента в уравнении (1), который составил 0,043, вместо 0,017 в исходном уравнении для гладкой поверхности внутренней трубы. Увеличение коэффициента объясняется изменением формы внутренней трубы и увеличением площади ее поверхности по сравнению с традиционной гладкой.В таблице 1 представлены экспериментальные и расчетные значения коэффициента теплоотдачи в испарительном блоке теплонасосной установки.Практическая значимость полученного уравнения заключается в возможности рассчитать количество теплоты, отбираемого в испарительном блоке, т.е. холодопроизводительность установки. определить холодильный и отопительный коэффициенты и, соответственно, рассчитать количество отданного теплонасосной установкой теплоты, т.е. теплопроизводительность.В испарительном блоке воздушного теплового насоса происходят сложные тепло-массообменные процессы, изучению которых посвящено много работ [16–18]. Процессы теплоотдачи от наружного воздуха к поверхности теплообменника сопровождаются образованием конденсата, так как температура теплообменной поверхности испарителя ниже точки росы подаваемого наружного воздуха. Поэтому в испарительном блоке образуется наледь, отрицательно влияющая на эффективность работы воздушной теплонасосной установки. Кинетика образования наледи имеет важное значение для регулирования тепломассообменных процессов теплонасосной установки и режима работы магнитострикцинного излучателя для удаления наледи. Таблица 1 Расчетные и экспериментальные значения коэффициента теплоотдачи в испарительномконтуре при различных температурах и скорости воздуха№п/пtср, оСtф, оСV, м/сρ, кг/м3λ, Вт/м Кν 105, м2/сRePrNuαрас, Вт/м2Кαэкс, Вт/м2К19,2-11,64,81,2490,02511,40658520,70439,0757,1564,6625,4-19,34,21,2660,02481,37452390,70535,5751,4651,5430,8-20,55,01,2880,02451,33664140,70742,0560,0065,034-4,0-25,64,41,3110,02411,29758150,70938,8754,6758,485-9,3-28,85,11,3370,02371,25469700,71145,0962,4065,416-17,1-30,25,21,3780,02311,19274750,71448,1965,0570,407-20,2-32,74,61,3950,02291,16867490,71544,4759,4365,698-26,4-32,44,01,4300,02241,12061190,71741,5654,4360,21910,0-6,712,81,2460,02511,413155320,70485,77125,76120,1105,0-13,413,01,2680,02481,371162550,70588,79128,30112,92110,2-16,412,51,2900,02441,331160950,70788,36125,86120,3712-5,0-19,713,01,3160,02401,289172890,70993,88131,63127,69 Как видно из таблицы, эффективность теплообмена в основном зависит от гидродинамических условий в межтрубном пространстве. На рис. 3 представлена зависимость экспериментальных и расчетных значений коэффициента теплоотдачи в межтрубном канале испарительного блока от критерия Рейнольдса, которая показала удовлетворительную сходимость. Рис. 3. Зависимость коэффициента теплоотдачи от числа Рейнольдса:1 – расчетные значения; 2 – экспериментальные значения   Кинетика процессов образования конденсата и инея зависит от температуры, относительной влажности, скорости воздуха, подаваемого в испарительный блок. Исследования показали, что наибольшее влияние оказывает влажность воздуха. Поэтому представляет интерес получение зависимости роста наледи в зависимости от относительной влажности воздуха.Влажность воздуха регулировали парогенератором, скорость движения воздуха в проходных каналах испарителя  регулировали расходом воздуха, который в блоке из четырех элементов составлял 400, 600 и 1200 м3/ч. Скорость воздуха устанавливалась в пределах: 5,1–5,5; 7,5–8,2, и 12–13 м/с соответственно. Исследования проводились  при температурах подаваемого воздуха  -5; -10; 0; 5; 10 °С. Относительную влажность поддерживали в диапазоне 40–100 %. Удельное количество влаги, (г влаги/кг сухого воздуха) конденсирующейся на поверхности теплообменника, рассчитывали с помощью i-d диаграммы по уравнению : ,                   (3)где d1 – влагосодержание воздуха на входе в испаритель (г /кг); d2 – влагосодержание воздуха при температуре газового трубопровода испарителя (г/кг).Значения влагосодержания определяли по i-d диаграмме при соответствующих параметрах воздуха, поступающего в испарительный блок, и температуре хдадагента в испарителе. Данные для расчетов взяты из протоколов проведенных испытаний. Анализ температурного поля в кольцевом сечении испарительного блока позволяет предположить, что вся сконденсировавшаяся влага превращается в иней плотностью 300 кг/м3. Количество образующего инея (W) зависит от массовой производительности установки по воздуху и может быть рассчитано по уравнению: ,                      (4)где G –массовый расход воздуха, кг/ч.Массовый расход воздуха находили по уравнению: ,                              (5)где L – объемный расход воздуха, м3/ч; ρ - плотность влажного воздуха кг / м3. .                         (6)Скорость образования инея рассчитывали по уравнению, мм/ч:  ,                         (7)где плотность инея ρи =300 кг/м3, F – поверхность теплообменника, 4 м2.Исследования показали, что для стабильной и эффективной работы воздушной теплонасосной установки сброс надели с поверхности теплообменника испарительного блока необходимо производить при толщине слоя h=3,0-5,0 мм.Расчет времени работы между циклами включения магнитострикционного излучателя рассчитывали по уравнению, ч:  .                              (8)Таким образом, можно рассчитать скорость образования наледи при заданном расходе, температуре и влажности наружного воздуха, и соответственно время между циклами работы до снятия наледи по уравнениям (3) – (8).В табл. 2 представлены результаты расчета скорости образования наледи при расходе воздуха L=600 м3/ч, G=600·1,31=786 кг/ч.Рассчитано время образования ледяной корки толщиной h=5мм, т.е. период работы установки до включения излучателя. Таблица 2 Результаты экспериментальных и расчетных значений кинетики роста наледи№ п/пW’t, оСφ, %W’W, кг/чv’, мм/чt, мин16,13-5606,044747,003,967525,320605,584385,883,658235,195605,124024,323,35368948,21-5808,046319,445,276057,40807,585957,884,9660,567,345807,125596,324,6664,38710,13-51009,667592,766,3271,2089,3201009,587529,886,2747,8599,1951009,127168,325,9750,25  В результате проведения планированного эксперимента было получено уравнение множественной регрессии, кг/кг:    ,   (9)Статистическая значимость уравнения проверена с помощью коэффициента детерминации и критерия Фишера. Как видно из таблицы, экспериментальные значения удельного количества конденсата (второй столбец) и расчетные значения (четвертый столбец) показывают близкие значения с расхождением 1–3 %Для массовой производительности500 кг/час по результатам эксперимента рассчитана скорость образования конденсата v’ в диапазоне изменения температуры Х1 ( t от -5 до +5  °С) и относительной влажности в диапазоне 60–100 % (Х2). Полученные данные описываются уравнением регрессии, мм/ч:  ,          (10)По результатам комплексных исследований получен патент на изобретение и разработана программа для ЭВМ. для автоматизации и мониторинга работы воздушного теплового насоса с системой MOVEBIT, включающие расчет теплопроизводительности, коэффициента трансформации, управление электрической частью излучателя в широком диапазоне температур [19, 20]. Полученные объекты интеллектуальной собственности подтверждают практическую значимость применения воздушных тепловых насосов с системой MOVEBIT в регионах с холодным климатом для создания с помощью инженерных систем требуемых параметров микроклимата в помещениях различного назначения.Выводы. В настоящей работе представлены результаты исследования влияния гидродинамических условий на тепло- и массообменные процессы в испарительном блоке теплонасосной установки. Получено уравнение зависимости критерия Нуссельта от критериев Рейнольдса, Прандтля, а также отношений температур наружного воздуха и фреонового холодопровода. Данное уравнение позволяет определить холодопроизводительность, теплопроизводительность и коэффициенты трансформации тепловой энергии в широком диапазоне температуры, расходов и скоростей движения наружного воздуха. На основании полученных данных разработаны программы для расчета теплопроизводительности и коэффициента трансформации теплоты для всех этапов жизненного цикла инженерных систем, включая: сбор, хранение, обработку и передачу данных. Получено уравнение кинетики роста наледи на поверхности теплообменника испарительного блока, позволяющее рассчитать время работы установки до включения магнитострикционного излучателя. Установлено, что удаление наледи целесообразно проводить при толщине корки не более 5 мм. Полученные результаты использованы для разработки программы работы электрической части излучателя.Таким образом, применение магнитострикционных излучателей в испарительном блоке воздушного теплового насоса позволяет использовать возобновляемые источники энергии в инженерных системах создания микроклимата помещений, с том числе в вентиляционных установках с рекуперацией теплоты, что в перспективе значительно расширит географию применения воздушных теплонасосных установок в том числе в регионах с холодным климатом.