CC BY
33УДК 621.561
С. В. Здитовецкая, кандидат технических наук, ассистент (БГТУ);
В. И. Володин, доктор технических наук, профессор (БГТУ)
КОМПЛЕКСНЫЙ МЕТОД РАСЧЕТА ПАРОКОМПРЕССИОННЫХ ТРАНСФОРМАТОРОВ ТЕПЛА
Разработан метод совместного численного анализа параметров цикла и теплообменного оборудования контура парокомпрессионных трансформаторов тепла, учитывающий нестационарный режим работы и необратимые потери в аппаратах и трубопроводах контура. Метод реализован в виде пакета программ и может использоваться при проектировании или выборе трансформатора тепла с учетом хладагента и действительного оборудования, входящего в его состав.
The method of the joint numerical analysis of loop variables and the heatexchange equipment of a contour of vapor-compression transformers of heat, considering non-stationary operating mode and irreversible losses in devices and contour pipelines is developed. The method is realized in the form of the software package and can be used at design or a choice of the transformer of heat taking into account a coolant and the valid equipment entering into its structure.
Введение. В условиях широкого применения в различных отраслях промышленности парокомпрессионных трансформаторов тепла, технические системы нового поколения должны иметь высокую надежность, длительный срок эксплуатации, невысокую металлоемкость. Трансформаторы тепла включают широкую номенклатуру устройств и в зависимости от области применения классифицируются на холодильные машины, тепловые насосы и комбинированные холодильно-нагревательные устройства.
В процессе внедрения или проектирования трансформатора тепла встает проблема выбора энергоэффективного устройства в целом или его комплектующих, в первую очередь, тепло-обменных аппаратов и компрессора. Сделать такой выбор достаточно сложно, потому что на эффективность влияют различные внутренние и внешние факторы, которые определяются условиями эксплуатации.
Наиболее успешно данную проблему можно решить, используя расчетные методы анализа и проектирования, которые позволяют провести исследование системы с учетом заданных требований. Анализ показал, что предлагаемые модели численного расчета не позволяют совместно рассчитывать параметры цикла и теп-лообменные аппараты контура, учитывать нестационарный режим работы и потери давления в элементах контура [1-7]. Таким образом, разработка комплексного метода анализа, который включает взаимосвязанный расчет параметров цикла и теплообменников контура па-рокомпрессионных трансформаторов тепла, учитывает нестационарный режим, работу компрессора и потери давления в аппаратах и трубопроводах контура, является актуальной и представляет научный интерес.
Основная часть. Парокомпрессионный трансформатор тепла представляет собой сложную техническую систему, включающую как основное, так и вспомогательное оборудование (рис. 1). К основному оборудованию относят испаритель, компрессор, конденсатор и терморегули-рующий вентиль.
Рис. 1. Схема трансформатора тепла: 1 - испаритель; 2 - регенератор;
3 - компрессор; 4 - конденсатор;
5 - охладитель;
6 - терморегулирующий вентиль
Разработанный метод расчета представляет собой взаимосвязанные блоки для расчета элементов оборудования и контура трансформатора тепла. На рис. 2 представлена укрупненная схема комплексного расчета, которая показывает возможные направления исследования.
Расчетная схема образует замкнутую систему и позволяет осуществлять, как совместный расчет параметров цикла и теплообменников, входящих в контур, так и проектирование отдельных элементов оборудования.
Теплообменные аппараты технологического контура
Расчет параметров компрессора
Проектный оптимизационный расчет теплообменных аппаратов
Поверочный расчет теплообменников
X
Рис. 2. Схема комплексного метода расчета трансформатора тепла
Комплексный метод расчета включает следующие пункты:
• определение параметров цикла при обратимом и необратимом сжатии рабочего вещества в компрессоре;
• расчет эксплуатационных режимов работы, т. е. совместный расчет параметров цикла и теплообменных аппаратов контура заданной конструкции в нестационарном или стационарном режиме с учетом действительной работы компрессора и потерь давления в элементах контура;
• комплексное расчетное проектирование, т. е. взаимосвязанный расчет параметров цикла с оптимизацией теплообменного оборудования контура с учетом действительной работы компрессора и потерь давления в элементах контура;
• индивидуальный проектный или поверочный расчет теплообменников трансформатора тепла и технологических контуров, связанных с источником и потребителем теплоты.
Принципиальным отличием представленного метода является возможность анализа работы трансформаторов тепла в нестационарном режиме, учет потерь давления в элементах контура и возможность унифицированного оптимизационного расчета теплообменных аппаратов различного конструктивного исполнения и назначения.
Теплообменное оборудование может иметь различное конструктивное исполнение, что потребовало разработки унифицированного подхода к методике его численного расчета, которая учитывает и нестационарный режим.
Основные теплообменные аппараты (испаритель и конденсатор) рассчитываются по двум участкам: с фазовым переходом и для однофазного потока. Каждый из участков описывается отдельными системами уравнений теплопередачи и теплового баланса. В этом случае общая поверхность теплообмена
где F\ - поверхность участка кипения (конденсации), м2; Fn - поверхность участка однофазного потока перегретого пара, м2.
При анализе работы трансформатора тепла в нерасчетных режимах граница между участками изменяется при F = const. Это позволяет получить адекватные результаты.
Нестационарный режим работы трансформатора тепла определяется работой испарителя и изменяющимися внешними условиями. Теп-лоперенос в испарителе описывается системой уравнений:
dQ = kAtF; dx
dQ = nGi [cp, 1 ((,1- *1,2) + hfg (1- л) ]; (2)
dQ
= П2М2S, 2 '
dt
F = F + Fn,
(1)
где Q - теплота, Дж; т - время, с; к - коэффициент теплопередачи, Вт/(м2К); Дt - средний температурный напор, оС; п - эффективность с учетом потерь при нагреве или охлаждении среды; G - массовый расход, кг/с; ср - теплоемкость, Дж/(кгК); t - температура, оС; к^ - теплота фазового перехода, Дж/кг; х - степень сухости пара; М - масса охлаждаемой среды, кг. Индексы: 1 - хладагент; 2 - охлаждаемый теплоноситель в испарителе или нагреваемый теплоноситель в конденсаторе; 1,1 - вход; 1,2 - выход.
При численном интегрировании принимается, что для достаточно малого интервала времени йт режим работы квазистационарный. Для исследуемых случаев минимальный интервал времени йт устанавливался в результате вычислительного эксперимента таким образом, чтобы конечный результат существенно не изменялся. На каждом временном интервале квазистационарной работы рассчитываются также конденсатор и вспомогательные аппараты.
Система уравнений (1)-(2) дополняется расчетом коэффициента теплопередачи и сред-нелогарифмического температурного напора с учетом поправки на схему течения рабочих сред [4]. Коэффициенты теплоотдачи для однофазных потоков и при конденсации рассчитываются по соотношениям, приведенным в [3-5], при кипении хладагентов - по данным [7, 8]. При проектном расчете трансформатора тепла замкнутая система уравнений решается относительно поверхности теплообмена, а при поверочном - находится тепловой поток как функция параметров источника и потребителя тепла. Одновременно с тепловым расчетом теплообменников определяются потери давления при течении потоков рабочих сред, обусловленные сопротивлением трения и местными сопротивлениями. Эксплуатационные характеристики компрессора определяются на основании функциональных зависимостей коэффициента подачи, индикаторного и электромеханического КПД от степени сжатия.
Метод позволяет проводить также оптимизационные расчеты теплообменников. В зависимости от задачи проектирования в качестве показателя оптимальности выступают масса аппарата, тепловой поток или другой требуемый показатель. Явные ограничения на функцию цели определяют диапазон изменения конструктивных параметров, неявные - предельные значения потерь давления, мощности на прокачку теплоносителя, коэффициента ореб-рения для ребристых поверхностей теплообмена. Минимум функции цели находится итеративным методом сеток с переменным шагом.
Разработанный метод расчета реализован в виде пакета программ ИБЛТТЯ [9, 10]. Пакет позволяет рассчитывать параметры работы трансформаторов тепла, выбирать наиболее оптимальные конструктивные решения при проектировании и выборе теплообменного оборудования.
Сформирована база данных с параметрами компрессорного оборудования и распространенных хладагентов Я22, Я134а, И218, Я290, Я600а, которая при необходимости может рас-
ширяться. Метод расчета применим в интервале температур хладагента от -40 до 150оС соответственно в испарителе и конденсаторе.
Достоверность разработанного метода анализа подтверждена сравнением с опытными данными в процессе проведения вычислительного эксперимента [11].
Проводилось сравнение данных для холо-дильно-нагревательной установки, полученных с помощью разработанного пакета программ, с данными, полученными с помощью программы Refrigeration Utilities [12], которая позволяет рассчитывать только параметры цикла и не учитывает работу оборудования, входящего в состав контура. Результаты сравнения приведены в таблице. Холодильно-нагревательная установка работает в нестационарном режиме, который определяется работой испарителя. Аппарат выполнен в виде непроточного теплообменника объемного типа. В испаритель загружается охлаждаемое молоко объемом 1000 л с начальной температурой 35оС. Для интенсификации процесса оно перемешивается мешалкой. Теплота, отводимая от молока, с повышением потенциала в конденсаторе используется для нагрева воды, применяемой на технологические нужды. В исследуемой установке конденсатор представляет собой однозаходный змеевиковый теплообменник с вытеснителем. Температура конденсации поддерживалась на уровне 50оС. Температура воды на входе в конденсатор составляла 5оС. При работе установки должно выполняться предельное требование по времени охлаждения молока - не более 3 ч до температуры 4оС. В качестве рабочего вещества использовался хладагент R22.
Для сравнения выбраны два квазистационарных режима, соответствующие 10 и 130 мин работы установки, которые приведены в таблице. Результаты сравнения параметров, соответствующие выбранным интервалам времени, согласуются удовлетворительно. Максимальное отклонение результатов по холодопроизво-дительности составляет 1,8%, а по теплопроиз-водительности - 1,9%.
Сравнение параметров цикла холодильно-нагревательной установки, полученных с помощью HEATTR и Refrigeration Utilities
Параметр Время работы установки
10 мин 130 мин
HEATTR Refrigeration Utilities HEATTR Refrigeration Utilities
Удельная холодопроизводительность, кДж/кг 152,0 153,7 145,2 147,8
Удельная теплопроизводительность, кДж/кг 189,20 191,02 203,39 207,38
Удельная работа, кДж/кг 37,82 37,28 58,21 59,54
Степень сжатия 2,81 2,79 4,58 4,56
Коэффициент преобразования 5,08 4,12 3,49 2,48
Заключение. Разработан комплексный метод, реализованный в виде пакета программ, для совместного анализа параметров цикла, теплообменных аппаратов и элементов обвязки контура парокомпрессионных трансформаторов тепла с учетом нестационарного режима работы и необратимых потерь.
Данный метод анализа позволяет проектировать трансформаторы тепла и прогнозировать действительно достигаемые ими параметры с учетом конкретного оборудования в номинальном и нерасчетных режимах работы.
Литература
1. Быков, А. В. Холодильные машины и тепловые насосы / А. В. Быков, И. М. Калнинь, А. С. Крузе. - М.: Агропромиздат, 1988. - 286 с.
2. Бубялис, Э. Процессы энергопереноса в тепловых насосах / Э. Бубялис, В. Макаряви-чус; под ред. А. Жукаускаса. - Вильнюс: Мокс-лас, 1990. - 186 с.
3. Шавра, В. М. Основы холодильной техники и технологии пищевых отраслей промышленности /В. М. Шавра. - М.: ДеЛи принт, 2002. - 126 с.
4. Бажан, П. И. Справочник по теплообменным аппаратам / П. И. Бажан, Г. Е. Кавенец, В. М. Селиверстов. - М.: Машиностроение, 1989. - 367 с.
5. Гопин, С. Р. Воздушные конденсаторы малых холодильных машин / С. Р. Гопин, В. М. Шав-ра. - М.: Агропромиздат, 1987. - 149 с.
6. Теплообменные аппараты холодильных установок / Г. Н. Данилова [и др.]; под общ.
ред. Г. Н. Даниловой. - Л.: Машиностроение, 1986. - 302 с.
7. Интенсификация теплообмена в испарителях холодильных машин / А. А. Гоголин [и др.]; под общ. ред. А. А. Гоголина. - М.: Легкая и пищевая пром-сть, 1982. - 223 с.
8. Баттерворс, Д. Теплопередача в двухфазном потоке / Д. Баттерворс, Г. Хьюитт; под ред. Д. Баттерворса. - М.: Энергия, 1980. - 328 с.
9. Здитовецкая, С. В. Пакет прикладных программ для комплексного анализа компрессионных тепловых насосов / С. В. Здитовецкая,
B. И. Володин // Энергетика: изв. высш. учеб. заведений и энерг. об-ний СНГ. - 2009. - № 5. -
C. 85-90.
10. Пакет программ для расчета параметров парокомпрессионных трансформаторов тепла: свидетельство № 402 / В. И. Володин, С. В. Здитовецкая; Белорус. гос. технол. ун-т. -№ С20120010; заяв. 06.02.12; опубл. 01.03.12. -Минск, 2012. - С. 1.
11. Володин, В. И. Тепловой расчет моло-коохладительной установки при нестационарном режиме работы / В. И. Володин, С. В. Здитовецкая // Heat and Mass Transfer: Proc. of V Int. Forum, Minsk, 24-28 May 2004, Vol. 10 [Электронный ресурс]. - Минск, 2004. - 1 электрон. опт. диск (CD-ROM): зв., цв. - 6 с.
12. Интернет-портал Technical University of Denmark [Электронный ресурс]. - Copenhagen, 2011. - Режим доступа: http://www.et.dtu.dt/Cool. -Дата доступа: 01.09.2011.
Поступила 20.02.2013
