УДК 621.576:614.8.084
ББК [31.392.3:31.368]:[35.112:52.81]
Ю. В. Шипулина, М. Ш. Каримов, М. Ф. Руденко
МЕТОДИКА ОПРЕДЕЛЕНИЯ ТЕПЛОВЫХ НАГРУЗОК НА РЕАКТОР ГЕНЕРАТОРА-АДСОРБЕРА ГЕЛИОЭНЕРГЕТИЧЕСКОЙ ХОЛОДИЛЬНОЙ УСТАНОВКИ
Yu. V. Shipulina, M. Sh. Karimov, M. F. Rudenko
METHODS OF DETERMINATION OF HEAT LOADS ON THE REACTOR OF THE GENERATOR-ADSORBER OF THE SOLAR-POWERED REFRIGERATING INSTALLATION
Рассматривается моделирование оптической части конструкции генератора-адсорбера. Приведена математическая модель расчета генератора-адсорбера, проанализировано распределение температурных полей и тепловых потоков в поперечном сечении в простой и оребренной конструкции генератора-адсорбера. Проанализированы различные конструктивные решения реактора, повышена эффективность его функционирования в гелиоэнергетической установке.
Ключевые слова: генератор-адсорбер, моделирование оптической части, зеркальные концентраторы, математическая модель, тепловые потоки, распределение температурных полей.
The modeling of the optical part of the generator-adsorber is considered. A mathematical model for calculating the generator-adsorber is presented and the distribution of temperature fields and heat flows in the cross section in a simple ribbed design of the generator-adsorber are analyzed. The different reactor designs are studied and the efficiency of their functioning in the solar-powered installation is increased.
Key words: generator-adsorber, modelling of the optical part, mirror concentrators, mathematical model, heat flows, distribution of temperature fields.
Гелиоэнергетические холодильные установки адсорбционного типа циклического автономного действия могут устанавливаться на морских буровых установках, стационарных островных и береговых объектах, в портах и терминалах и предназначаются для получения охлажденной питьевой воды, кондиционирования и охлаждения помещений, замораживания и хранения медикаментов и скоропортящейся продукции и других нужд.
В гелиоэнергетических адсорбционных холодильных установках генератор-адсорбер совмещенного типа [1-4] является одним из важных аппаратов. Основной элемент генератора-адсорбера - реакторы, установленные в гелиоприемное устройство типа «горячий ящик» между зеркальными концентраторами. В реакторы генератора-адсорбера засыпается адсорбент, например активированный уголь. Под действием солнечной энергии днем происходит обогрев наружных поверхностей корпуса реактора, теплота через стенки передается насыщенному аммиаком адсорбенту и с определенной температуры начинается выделение адсорбата из адсорбента, который в виде газообразной фазы сжижается в конденсаторе. Этот период времени работы гелиоэнергетической холодильной установки называется периодом регенерации. При охлаждении сухого адсорбента (активированного угля) в ночной период времени в реакторах начинается поглощение аммиака из испарителя - начинается период зарядки. Работа гелиоэнергетической холодильной установки протекает циклично - так же циклично, как изменяется температура окружающей среды.
Эффективность работы гелиоэнергетической холодильной установки циклического действия во многом зависит от работы реакторов генератора-адсорбера. Обзор отечественных и зарубежных литературных источников и исследование собственных конструкций гелиоэнер-гетических холодильных установок выявил недостатки в работе совмещенных генераторов-адсорберов, состоящих из гелиоприемного устройства типа «горячий ящик», концентрирующих зеркальных элементов, реакторов и некоторых других элементов. Это требует перспективной проработки конструкции именно этого аппарата. Основными недостатками аппарата являются: недостаточное количество теплоты (в отдельные дни) для полного обогрева всей массы насы-
щенного сорбента в реакторах при суточном освещении солнцем и, как следствие, низкие значения температуры в объеме насыпной массы сорбента; снижение эффективности работы вследствие значительного объема воздуха под прозрачной оболочкой гелиоприемника и большого количества реакторов в одной конструкции генератора-адсорбера; повышение стоимости аппарата в связи с увеличением массогабаритных характеристик.
Конструирование оптической части гелиоприемного устройства генератора-абсорбера гелиоэнергетической холодильной установки можно начать с элементов корпуса реактора и зеркал в аппарате.
Методики расчета и анализ эффективности оптических характеристик гелиоприемных устройств с плоскими зеркальными концентраторами на различные формы адаптирующих поверхностей были рассмотрены нами в [5-7]. На рис. 1 представлены соотношения оптического моделирования конструкции генератора-адсорбера. В основу моделирования оптической части положено два принципа: первый - угол падения солнечного луча на зеркальную поверхность полностью отражается, и угол падения равен углу отражения луча; второй - вся энергия солнечного луча (прямого и отраженного), падающая на адаптирующую поверхность, целиком поглощается. Гелиоприемное устройство генератора-адсорбера представляет собой удлиненную конструкцию (на рис. 1 представлен профиль конструкции), состоящую из плоских зеркальных поверхностей, имеющих определенный угол раскрытия 0, град, расположенных симметрично относительно центральной плоскости.
W
Рис. 1. Моделирование оптической части конструкции генератора-адсорбера гелиоэнергетической холодильной установки
В симметричной части зеркал устанавливается круглая трубка корпуса реактора с наружным радиусом r, мм. Ширина раскрытия зеркал W, мм. Полная длина конструкции аппарата
l, мм. Q и G, мм2, - поверхности зеркал, отражающих солнечные лучи. S, мм2, - поверхность трубки цилиндрического реактора, адаптирующая солнечное излучение. F, мм2, - поверхность металлической подложки, адаптирующей солнечное излучение. На основании анализа концентрической эффективности солнечного облучения круглой поверхности в плоских зеркалах при соотношении W/2R = 3,5 был выбран угол раскрытия зеркал 0, град. Длина зеркал Lз рассчитана таким образом, чтобы при перпендикулярном облучении концентраторов (вдоль оси симметрии аппарата) и при однократном отражении лучей от верхней до нижней кромки поверхностей G и Q солнечные лучи полностью попадали на внешнюю поверхность S реактора. При таких условиях и перпендикулярном освещении солнцем цилиндрическая конструкция корпуса реактора будет иметь трехкратное облучение энергией солнечной радиации. Кроме того, нижняя подложка DE с поверхностью F работает как прямое ребро и передает дополнительную теплоту на нижнюю часть реактора при хорошем контакте с корпусом. Тогда основные элементы конструкции аппарата можно выразить через соотношения радиуса реактора R и фиктивного угла раскрытия зеркал 0.
h1 = R-cos (0/2) /sin (0/2); (1)
L0 = R /sin (0/2); (2)
W/2 = 3,5-R; (3)
L3 = (W/2)/sin (0/2) - R /sin (0/2); (4)
H = (L0 + L3) /cos (0/2). (5)
На основании оптических параметров моделирования и баланса энергетических потоков разработана теплоэнергетическая модель конструкции генератора-адсорбера гелиоэнергетиче-ской холодильной установки, профиль конструкции которого представлен на рис. 2. Конструкция трапецеидальной формы имеет прочный металлический корпус 1 и изоляцию 2 с пяти сторон (торцевые стороны конструкции не показаны), зеркальные отражатели 3, внутренний стальной корпус реактора 4 с активированным углем 5, двухстекольное покрытие 6 со стороны облучения солнцем и металлическую подложку под реактором 7.
Рис. 2. Модель генератора-адсорбера гелиоэнергетической холодильной установки:
1 - корпус элемента гелиоприемника; 2 - изоляция корпуса; 3 - плоские зеркальные отражатели;
4 - наружный корпус реактора; 5 - адсорбент (активированный уголь); 6 - двухстекольное покрытие;
7 - металлическая подложка; 8 - крепежные болты
На корпус реактора, как видно на модели, подводится теплота в виде потока dq 1(х) - лучистой энергии на видимую поглощающую поверхность площадью лЯ7 и отраженной dq2(т) от зеркал О, Q на поглощающую поверхность площадью 2-пЯ7, в виде потока dq3(т) - от теплоты, подводимой за счет металлической плоской поверхности подложки.
й^1(т) [й^пр.к(т) - й^пр.к(т)Ат - dq1пр,(T)Бт], (6)
где dqпр.к(т) - доля полезной теплоты, проходящей в «горячий ящик» и идущей на обогрев корпуса реактора, Вт; Ат, Бт - соответственно коэффициенты поглощения и отражения солнечной энергии от трубки.
dq2(т) №пр.к(т) - й^пр.к(т)Ат - ^пр.«^^ (7)
где Со - коэффициент оптической концентрации на поглощающую поверхность, Со = (2-пЯ7 )/(2-л^Я^/).
dqз(т) = А • 5о• т• йТо(т) ґк(тЬ), (8)
где Ат, А0 - коэффициент теплопроводности материала трубы и металлической подложки соответственно, Вт/(м- К); 5т, 5о - толщина трубы и подложки соответственно; Я, г - внешний и внутрен-
ний радиусы трубы, м; То - температура у основания подложки, К; т = [2 • а/(А0 • 50)]°’5 - комплекс
теплофизических переменных; а - коэффициент теплоотдачи от металлической подложки конвекцией, Вт/(м2К); b = 2т - длина ребра, м; th (m b) = [e(mb) + e<mb')y[e(-mb') - e<mb)] - функция гиперболического тангенса, число е = 2,7183.
Теперь рассмотрим, как и из чего формируются полезная тепловая нагрузка на реактор генератора-адсорбера, равная
dq 'пол(т) = dq i(t) + + ^(т). (9)
Величина теплового потока, падающего на остекленную оболочку генератора-адсорбера, выполненного по типу конструкции «горячего ящика», выражается следующим уравнением:
йп.д = Z dq пад (т) =( ^(т) ■ C0s i + dId (т) ) W ■ h (10)
где - мощность полной падающей солнечной энергии на поверхность гелиоприемника за время облучения в течение светового дня, Вт; dq^ (т) - мощность падающей солнечной энергии за конкретный промежуток времени (т = 1 ч); l - длина реактора генератора-адсорбера, м; dIд(Т) - мощность потока прямого солнечного излучения, Вт/м2; dId(x) - мощность потока диффузионной солнечной радиации, Вт/м2; cos i - коэффициент перевода потока солнечной энергии на перпендикулярную к гелиоприемнику поверхность.
cosi = cos(90 - ф + ^)sin5 + sin(90 - ф + ц) cos5■ cosт , (11)
где i - угол падения солнечных лучей на перпендикулярную поверхность гелиоприемника, град; ф - угол широты местности, град; ц - угол наклона поверхности гелиоприемника к горизонту земли, град; 5 - угол склонения солнца (функция времени года, по астрономическому календарю), град; т - часовой угол движения солнца, град.
Тогда тепловой поток в единицу времени, проходящий через остекленную оболочку в «горячий ящик», будет выражаться уравнением
dqпр(т) Кпз Кпп [й^пад(т) - й^пад(т) А - й^пад(т) Б - 8ст со(Тст(т)/100) ] (12)
где Кпз, Кпп - соответственно коэффициент потери энергии при прохождении потока солнечной радиации от запыленности стекла и коэффициент потери энергии из-за двухстекольного покрытия в конструкции «горячего ящика»; A, Б - соответственно коэффициенты поглощения и отражения солнечной энергии от стекла; 8ст со (Гст(т)/100)4 - собственное излучение стекла по формуле Стефана - Больцмана, Вт; 8ст - степень черноты стекла; со - коэффициент излучения абсолютно черного тела; Гст(т) - абсолютная температура стекла в единицу времени, К.
Количество солнечной энергии, падающей на поверхность трубки корпуса реактора, будет складываться из следующих составляющих: энергия от потока на видимую часть трубки; энергия отраженного потока от плоских зеркальных концентраторов; энергия от нижней подложки под цилиндрическим реактором.
Следует учесть, что проходящая в «горячий ящик» солнечная энергия расходуется также на подогрев элементов конструкции генератора-адсорбера (металлические части, изоляционный материал конструкции «горячего ящика», зеркал и т. п.), воздуха в замкнутом объеме и потери через внешние стенки «горячего ящика» (как изолированные поверхности, так и верхнюю остекленную часть) за счет разности значений температуры нагретых элементов аппарата и температуры внешней окружающей среды, причем эти потери сильно зависят от внешних условий: от скорости обдувания внешней поверхности корпуса ветром и направления ветра, от толщины и качества изоляции, от толщины и количества стекол гелиоприемного устройства, от влажности.
Энергию на подогрев элементов конструкции генератора-адсорбера можно определить так:
dqi (т)= Yjmi'Cibti , (13)
где mt - масса соответствующего материала i-го элемента конструкции; ci - теплоемкость материала i-го элемента конструкции; Ati - разность значений температуры i-го элемента конструк-
ции за период между замерами и расчетами тепловых потоков.
Энергию на подогрев воздуха в замкнутом объеме гелиоприемника можно определить по следующей формуле:
dq^) = 0,9 (4 • V 8в • с0/Я-^ГСр.в(т)/100)4, (14)
где dq^^) - собственное поглощение энергии воздухом в единицу времени, Вт; V = {12,25 т2*
*[tg (0/2)] - лт2 - r2/[sin (0/2)-l] }-l - объем воздуха в замкнутом пространстве «горячего ящика»,
м3; Н = 3,5-r-tg (0/2) - толщина воздушного слоя, м; dTp.^) - средняя абсолютная температура воздуха внутри «горячего ящика» в единицу времени, К; 8в - степень поглощения (черноты) воздуха.
Определим энергию потерь из «горячего ящика» в окружающее пространство через изолированную поверхность и стекло. Условно считаем, что элемент конструкции гелиоприемника имеет верх площадью W-l, м2, дно площадьюМдна = 1-{[3,5тtg (0/2) - т/tg (0/2)] л [90 - (0/2)]/180}, м2; две боковые стенки общей площадью 2Ьз • l, м2; два торца площадью Ттор = 2-[(2т + 2-5ю)(Н - hj) + + (2,5 т + 5из) 2 (H- hj) + лт2 (1/2 - 0/360) - т2] l, м2.
Тогда
dq’пот(т) = К’ (МдНа + 2-Ьз-1 + ГТОр) [dT^Cr) - dT^Cr)]; (15)
dq”пот(т) = К”- W-l ^Гср.в(т) - dT^)], (16)
где dq'пот(т) - потери тепла через изолированные элементы конструкции генератора-адсорбера в единицу времени, Вт; dq’nm-CO - потери тепла через элемент остекления конструкции «горячего ящика» в единицу времени, Вт; K’, K” - соответственно коэффициенты теплопередачи через элементы изолированной конструкции и стекло, Вт/(м2К); dT^^) - изменение температуры наружного воздуха в единицу времени, К.
К’ = 1/[(1/ав) + (5юДиз) + (1/анар)],
К” = 1/[(1/ав) + [(2-5сДст)+(5вАв)] + (1/анар)],
где ав, анар - коэффициенты теплоотдачи от внутренней и наружной поверхности плоской конструкции «горячего ящика» соответственно, Вт/(м2К); 5из,5ст, 5в, Аиз, Аст, Ав - соответственно толщины, м, и коэффициенты теплопроводности, Вт/(мК) изоляции, стекла, воздуха между стеклами; ав = const; анар - зависит от условий внешней среды и ее можно рассчитать по уравнению [8].
При умеренном и сильном ветре
анар = 0,0296-^нар Яев0ЛРГв0,43ЧРГв/РГст)0,25/Г,
где Анар - коэффициент теплопроводности наружного воздуха, Вт/(м. К); Г - параметр конструкции «горячего ящика», обдуваемого ветром, м; Яев = (ювТ)^в - критерий Рейнольдса; шв^- скорость ветра, м/с; vв - коэффициент кинематической вязкости наружного воздуха, м2с; Ргв • = ув/ав -критерий Прандтля для воздуха при наружной температуре; Ргст - критерий Прандтля для воздуха при температуре наружной поверхности «горячего ящика»; ав - коэффициент температуропроводности воздуха, м2.с.
Тогда доля полезного тепла, проходящего в «горячий ящик» и идущего на обогрев корпуса реактора, будет равна:
й^пр.к(т) dqпp(т) - dqсв(т) - dq пот(т) - dq пот(т) - dqi. (17)
При высоких значений температуры наружной стенки трубы следует учитывать собственное излучение ее поверхности, равное:
dqСТ(т) = 2 • л • R 18т • со • dTТ(т)/100)4. (18)
Тогда суммарное полезное тепло, передаваемое насыщенному адсорбенту, будет равно:
dqпол(т) = dq1(т) + dq2(т) + dq2(т) - dqCTCt). (19)
Полное полезное тепло, передаваемое реактору генератора-адсорбера за световой день, с 6 часов утра до 18 часов вечера, будет равно:
т=18
Q^ = f dqПол (т). (20)
т=06
Для расчета полезных тепловых нагрузок, воспринимаемых реактором генератора-адсорбера гелиоэнергетической холодильной установки, на основании формул (1)-(20) была разработана расчетная программа в программе Mathcad.
Расчетные значения полезной тепловой нагрузки по разработанной программе были проверены экспериментально на собранном стенде, представляющем конструкцию генератора-адсорбера. Замеры солнечной радиации определялись по показаниям пиранометра и актинометра, температура поверхностей корпуса реактора, стекол и других элементов - хромель-копелевыми термопарами с цифровым вольтметром, параметры теплоносителей - термометрами. Через корпус конструкции пустого реактора пропускали водоглицериновый раствор.
Полезная тепловая нагрузка за контрольный промежуток времени определялась по следующей формуле:
пол(х) ^ж'£ж'рж(Твых — Твх)
где Gж - объемный расход жидкого теплоносителя, м3/с; сж - теплоемкость жидкого теплоносителя, Дж/(кг-К); рж - плотность жидкого теплоносителя, кг/м3; Твых, Твх-- температура жидкого теплоносителя на входе и выходе в корпус реактора, оС.
На рис. 3 представлены графики замеров прямой и диффузионной составляющих солнечной радиации на открытом полигоне в г. Астрахани в июне месяце и расчетные значения тепловых нагрузок по приведенным выше уравнениям на корпус реактора генератора-адсорбера ге-лиоэнергетической холодильной установки. Здесь же приведены значения полезных тепловых нагрузок по экспериментальным данным теплотехнических замеров.
t ч
Рис. 3. Графики тепловых нагрузок на реактор генератора-адсорбера в течение светового солнечного дня в июне месяце: 1 - параметры прямой солнечной радиации; 2 - параметры диффузионной солнечной радиации; 3 - расчетные значения тепловых нагрузок; 4 - экспериментальные значения тепловых нагрузок
На графиках видно удовлетворительное совпадение расчетных данных с экспериментальными. Полного совпадения не происходит из-за тепловой инерции конструкции гелиоприемно-го устройства, которая в данной методики не учитывалась. Однако для практических расчетов она большого значения не имеет.
СПИСОК ЛИТЕРАТУРЫ
1. Руденко М. Ф. Разработка и исследование эффективности адсорбционной гелиохолодильной установки // Вестн. Междунар. академии холода. - 2003. - № 1. - С. 34-37.
2. Buchter F., Dind P. H., Pons M. An experimental solar-powered absorptive refrigeration tested in Burkina-Faso // Int. J. Refrigeration. - 2003. - 26. - P. 79-86.
3. Thermal conductivities study of new types of compound adsorbents used in solar adsorption refrigeration / Z. Huizhong, Z. Min, L. Jing, Y. Guoding, Z. Zhijun // Energy Convers Manage. - 2009. - Vol. 50. - P. 1244-1248.
4. Boubakri A. A new conception of an adsorptive solar powered ice maker // Renewable Energy. - 2003. -Vol. 28. - P. 831-842.
5. Руденко М. Ф. Эффективность гелиоприемных устройств с концентраторами для систем тепло-и хладоснабжения. - Саратов: Лаборатория нетрадиционной энергетики ОЭП СНЦ РАН, 2001. - 64 с.
6. Ильин А. К., Руденко М. Ф., Коноплева Ю. В. Оценка эффективности концентраторов солнечной энергии // Изв. вузов. Машиностроение. - 2002. - № 9. - С. 33-36.
7. Руденко М. Ф., Туркпенбаева Б. Ж., Чивиленко Ю. В. Теоретические основы проектирования экологически безопасной гелиоэнергетической техники для производства теплоты. - Актау: Изд-во АктГУ им. Ш. Есенова, 2008. - 165 с.
8. Исаченко В. П., Сукомел А. С., Осипова В. А. Теплопередача: учеб. для вузов. - М.: Энергоиздат, 1981. - 416 с.
REFERENCES
1. Rudenko M. F. Razrabotka i issledovanie effektivnosti adsorbtsionnoi geliokholodil'noi ustanovki [Development and research of efficiency of adsorbing solar-powered installation]. VestnikMezhdunarodnoi akademii kholoda, 2003, no. 1, pp. 34-37.
2. Buchter F., Dind P. H., Pons M. An experimental solar-powered absorptive refrigeration tested in Burkina-Faso. Int. J. Refrigeration, 2003, 26, pp. 79-86.
3. Huizhong Z., Min Z., Jing L., Guoding Y., Zhijun Z. Thermal conductivities study of new types of compound adsorbents used in solar adsorption refrigeration. Energy Convers Manage, 2009, vol. 50, pp. 1244-1248.
4. Boubakri A. A new conception of an adsorptive solar powered ice maker. Renewable Energy, 2003, vol. 28, pp. 831-842.
5. Rudenko M. F. Effektivnost’ geliopriemnykh ustroistv s kontsentratorami dlia sistem teplo- i khlados-nabzheniia [Efficiency of solar receivers with concentrators for the systems of heating and cooling]. Saratov, La-boratoriia netraditsionnoi energetiki OEP SNTs RAN, 2001. 64 p.
6. Il'in A. K., Rudenko M. F., Konopl'va Iu. V. Otsenka effektivnosti kontsentratorov solnechnoi energii [Assessment of the efficiency of concentrators of solar energy]. Izvestiia vuzov. Mashinostroenie, 2002, no. 9, pp. 33-36.
7. Rudenko M. F., Turkpenbaeva B. Zh., Chivilenko Iu. V. Teoreticheskie osnovy proektirovaniia eko-logicheskoi bezopasnosti gelioenergeticheskoi tekhniki dlia proizvodstva teploty [Theoretical foundations of designing of ecologically safe solar-powered equipment for heat production]. Aktau, Izd-vo AktGU im. Sh. Esenova, 2008. 165 p.
8. Isachenko V. P., Sukomel A. S., Osipova V. A. Teploperedacha [Heat transfer]. Moscow, Energoizdat, 1981. 416 p.
Статья поступила в редакцию 12.10.2012
ИНФОРМАЦИЯ ОБ АВТОРАХ
Шипулина Юлия Викторовна - Астраханский государственный технический университет; канд. техн. наук; доцент кафедры «Безопасность жизнедеятельности и гидромеханика»; aleera78@mail.ru.
Shipulina Yulia Victorovna - Astrakhan State Technical University; Candidate of Technical Sciences; Assistant Professor of the Department "Life Security and Hydromechanics"; aleera78@mail.ru.
Каримов Марат Шайдулаевич - Астраханский государственный технический университет; соискатель кафедры «Безопасность жизнедеятельности и гидромеханика»; aleera78@mail.ru.
Karimov Marat Shaidulaevich - Astrakhan State Technical University; Postgraduate Student of the Department "Life Security and Hydromechanics"; aleera78@mail.ru.
Руденко Михаил Фёдорович - Астраханский государственный технический университет; д-р техн. наук, профессор; зав. кафедрой «Безопасность жизнедеятельности и гидромеханика»; mf.rudenko@mail.ru.
Rudenko Mikhail Fedorovich - Astrakhan State Technical University; Doctor of Technical Sciences, Professor; Head of the Department "Life Security and Hydromechanics"; mf.rudenko@mail.ru.