CC BY
94
УДК 662.997
М. Ф. Руденко, Б. Ж. Туркленбаева*
Астраханский государственный технический университет Актауский государственный университет им. Ш. Есенова Республика Казахстан
МОДЕЛИРОВАНИЕ ЭНЕРГЕТИЧЕСКОГО БАЛАНСА И РАСЧЕТ ГЕЛИОПРИЕМНЫХ УСТРОЙСТВ ДЛЯ ВОДОПОДОГРЕВА
Введение
Из существующих способов подогрева воды в районах с жарким климатом наиболее экономически целесообразным является солнечный нагрев в так называемых гелиоэнергетических устройствах типа «горячий ящик».
Ниже приведено описание аналитической модели энергетического баланса и предлагается методика расчета гелиоприемных устройств типа «горячий ящик» для проектирования солнцеиспользующих водонагревательных установок, которые отличаются от моделей, описанных в [1-3].
Построение модели теплофизических явлений
Для расчета тепловых нагрузок, воспринимаемых поверхностью гелиоприемных систем, обогреваемых солнцем, построим математическую модель физических явлений, происходящих в гелиоэнергетических устройствах типа «горячий ящик» (рис. 1).
Рис. l. Модель гелиоэнергетического устройства: 1 - корпус изолированный; 2 - ребра плоские;
3 - трубка для теплоносителя; 4 - остекленная поверхность
Удельный тепловой поток, падающий на остекленную оболочку элемента «горячего ящика» в данный момент времени, выражается следующим уравнением:
dg^(x) = [dI4(x)-cos i + dId(x)] F/l, (l)
где F = l-k-n - площадь солнечной энергии, падающей на поверхность гелиоприемника, м;
l - длина трубок, м; k - элемент ширины конструкции гелиоприемника, м; n - число элементов конструкции; оІд(т) - поток прямого солнечного излучения на l м2 поверхности, горизонтальной к земле (функция от времени), Вт/м2; dl^x) - поток диффузионной (отраженной от сферы) солнечной радиации на l м2 поверхности (функция от времени), Вт/м2; cos i - коэффициент перевода потока солнечной энергии на перпендикулярную к гелиоприемнику поверхность.
3l
cos i = cos(90 - ф + ^)-sin 5 + sin(90 - ф +ц) • cos 5 • cos т ,
где i - угол падения солнечных лучей на перпендикулярную поверхность гелиоприемника, град; ф - угол широты местности, град; ц - угол наклона поверхности гелиоприемника к горизонту земли, град; 5 - угол склонения солнца (функция времени года, по астрономическому календарю), град; т - часовой угол движения солнца, град.
Тогда тепловой поток в единицу времени, проходящий через элемент остекленной оболочки в «горячий ящик», будет выражаться уравнением
^^пр(т) Кп.з Кп.п [$^пад(т) - ^^пад(т) А- ^^пад(т) R - 8ст с0 (Тст(т)/100) ], (2)
где Кпз, Кпп - соответственно коэффициент потери энергии от запыленности стекла при прохождении потока солнечной радиации и коэффициент потери энергии за счет непрозрачных переплетов элементов конструкции «горячего ящика»; A, R - соответственно коэффициенты поглощения и отражения солнечной энергии от стекла; 8ст со (Тст(т)/100)4 - собственное излучение стекла по формуле Стефана - Больцмана, Вт; 8ст - степень черноты стекла; с0 - коэффициент излучения абсолютно черного тела; Тст(т) - абсолютная температура стекла в единицу времени, К.
Количество солнечной энергии, падающей на поверхность элемента конструкции, будет складываться из трех составляющих: энергии от потока на видимую часть трубки; энергии от потока на плоскую поверхность между трубками, которая является плоскими ребрами; энергии от конвективного теплообмена всей поверхности трубки с воздухом в «горячем ящике».
Следует учесть, что солнечная энергия, проходящая в «горячий ящик», расходуется также на подогрев воздуха в замкнутом объеме и потери через внешние стенки ящика (как изолированные, так и остекленную) за счет разности температур с внешней окружающей средой, причем эти потери сильно зависят от наличия паров воды в воздушном объеме ящика, от скорости обдувания внешней поверхности ящика ветром и направления ветра, от толщины и качества изоляции, от толщины и количества стекол гелиоприемного устройства.
Количество энергии на подогрев воздуха в замкнутом элементе объема гелиоприемника определим следующим образом: условно выделим элемент шириной k из конструкции «горячего ящика», тогда
^^св(т) = 0,9 • (4 • V- 8в • с{)/И)(ёТСрв(т)/100),4 (3)
где а?дсв(т) - собственное поглощение энергии воздухом в единицу времени, Вт; V = l• kh - объем воздуха в замкнутом пространстве элемента конструкции, м3; h - толщина воздушного слоя, м; Т в(т) - средняя абсолютная температура воздуха внутри «горячего ящика» в единицу времени, К; 8в - степень поглощения (черноты) воздуха.
Определим энергию потерь из «горячего ящика» в окружающее пространство через изолированную поверхность и стекло ящика. Условно считаем, что элемент конструкции гелиоприемника имеет верх и дно площадью l k, м2, каждое и четыре боковые стенки с поверхностями соответственно 2• kh и 2• l• h/n, м2. Тогда
^д'пот(т) = K’ • (l • k + 2 • l• h/n) • d^. в(т) - Тн(т)], (4)
dq”mT(T) = K”• l• k d[Tcv, в(т) - Тн(т)], (5)
где dq'поT(т) - потери тепла через изолированные элементы конструкции коллектора в единицу времени, Вт; dg’noT-Ci;) - потери тепла через элемент остекления «горячего ящика» в единицу времени, Вт; K’, K” - соответственно коэффициенты теплопередачи через элементы изолированной конструкции и стекло, Вт/(м2 • К); dTK(T) - изменение температуры наружного воздуха в единицу времени, К.
К’ = 1/[(1/ав) + (5юДиз) + (1/анар)],
К” = 1/[(1/ав) + (5СДСТ) + (1/анар)],
где ав, анар - соответственно коэффициенты теплоотдачи от внутренней и наружной поверхности плоской конструкции «горячего ящика», Вт/(м2 • К); ав = const - можно считать величиной постоянной; анар - зависит от условий внешней среды, и ее можно рассчитать по следующему уравнению [6]:
анар = 0,0296 • ^нар Яев0,8 • Ргв0,43 • (Ргв/РгСт)0,25/Г,
где ^нар - коэффициент теплопроводности наружного воздуха, Вт/(м • К); Г - параметр «горячего ящика», обдуваемого ветром, м; Яев = юв • РЬв - критерий Рейнольдса (юв • - скорость ветра, м/с; vв -коэффициент кинематической вязкости наружного воздуха, м2 • с); Ргв = vв/ав - критерий Пран-дтля для воздуха при наружной температуре; Ргст - критерий Прандтля для воздуха при температуре наружной поверхности «горячего ящика»; ав - коэффициент температуропроводности воздуха, м2 • с; 5из, 5ст - толщина соответственно изоляции и стекла, м; ^из, ^ст - коэффициенты теплопроводности соответственно изоляции и стекла, Вт/(м2 • К).
Тогда доля полезного тепла, проходящего в коллектор и идущего на обогрев трубок, будет равна:
^пр.кСО <^пр(т) - dqсв(т) - dq пот(т) - dq пот(т) .
На трубки, как отмечалось выше, подводится тепло в виде потока dq1(т) лучистой энергии на видимую часть поглощающей поверхности площадью п Я• l• п, в виде потока dq2(T) от теплоты, подводимой за счет оребрения плоской поверхности, и в виде доли потока dq3(т) от конвективного теплообмена от нагретого воздуха в «горячем ящике» к стенке трубок поверхностью
2 • п Я1п.
dql(т) = [dqпр. к(т) - dqщ, к(т) А - dq1пр. к(т)• Ят][(п Я• Iп)/(кI• п)], (7)
где - Ат, Ят - соответственно коэффициенты поглощения и отражения солнечной энергии от трубки.
dq2(т) = • 5о • т • dТо(x) • 1к(т • Ь)(к - п • Я), (8)
dqз(т) = {2 • п • г• Id(TCр. в(т) - Тж(т)) / [1/ав + 5Т • (^(г/Я))Дт + г/(аж \ • Я)}, (9)
где ^т, 5т, 5о, Я, г - соответственно коэффициент теплопроводности материала трубы и оребре-
ния, Вт/(м • К), толщины трубы и оребрения, внешний и внутренний радиусы трубы, м;
То, Тж, аж - соответственно температура у основания ребра, К; абсолютная средняя температура жидкости в трубе, К, и средний коэффициент теплоотдачи теплоносителя в трубе, Вт/(м2 • К); т = [2• ат/(^т • 5о)]0,5 - комплекс теплофизических переменных; Ь = к/2 - длина ребра, м; етЬ + е~(тЬ)
Ж(т• Ь) = ---_(тЬ) , где е = 2,7183 - функция гиперболического тангенса.
При высокой температуре наружной стенки трубы следует учитывать собственное излучение ее поверхности, равное:
dqCТ(т) = 2 • п Я • I• 8т• со• dТТ(т) / 100)4. (10)
Тогда суммарное полезное тепло, передаваемое жидкому теплоносителю, протекающему внутри трубы, будет равно:
dqпол(т) = dq1(x) + dq2(x) + dqз(т) - dqст(т). (11)
Полное полезное тепло, передаваемое теплоносителю за световой день, будет равно:
т =18
бпол = Е dqПOЛ СО, (12)
т = 6
где т - время сияния солнца в дневное время: 6 часов - рассвет, 18 часов - заход солнца.
При ламинарном течении жидкости внутри трубы (естественная циркуляция теплоносителя) коэффициент теплоотдачи аж можно найти по уравнению [4]
Маж = 0,15 • Яеж0,33 • Огж0,1 • Ргж0,43, (13)
где №ж = 2г • аж/^ж - критерий Нуссельта; Яеж = 2гтаж^ - критерий Рейнольдса для жидкого теплоносителя при скорости протекания юж; Огж = g • (2 • г)3рж АТ^ж - критерий Грасгофа; Ргж = vж/аж -критерий Прандтля, где vж, а*, рж - соответственно коэффициент кинематической вязкости, м2 • с; коэффициент температуропроводности, м2 • с; коэффициент объемного расширения жидкого теплоносителя в пределах изменения температур АТ; g - ускорение свободного падения, м/с2.
Проверка модели теплофизических явлений
Для расчета тепловых потоков, воспринимаемых коллектором гелиоприемника, рассчитанных на основании аналитической модели по формулам (1)-(13), была разработана расчетная программа в пакете МаШса^
На рис. 2 представлены расчетные зависимости тепловых потоков. Значения коэффициентов поглощения, отражения и пропускания выбирались из справочной литературы [5]. На этом рисунке представлены также экспериментальные значения полезных тепловых нагрузок, полученные при прохождении теплоносителя через солнечный коллектор (по данным [4]) и рассчитанные по следующей формуле:
Qпол. ^ж сж рж (Твых ТвхХ
3/
(14)
где Gж - объемный расход жидкого теплоносителя, м /с; сж - теплоемкость жидкого теплоносителя, Дж/(кг • К); рж - плотность жидкого теплоносителя, кг/м3; Твых, Твх - температура жидкого теплоносителя на выходе и входе в гелиоэнергетический коллектор, К.
Время, ч
Рис. 2. Зависимости тепловых потоков солнечной радиации в течение светового дня
Сдвиг расчетных и экспериментальных значений связан с тепловой инерцией экспериментального гелиоприемника, которая зависит от конструкции, массы, температуры нагрева гелиоприемной части и в первой половине дня оказывает отрицательное воздействие на нагрев теплоносителя, а в конце дня удлиняет время подогрева.
Подобрав графической интерполяцией коэффициент инерции гелиоприемного коллектора, позволяющего сдвинуть значения расчетных полезных тепловых нагрузок, полученные по формулам (1)-(13), по оси времени на экспериментальные значения полезных тепловых нагрузок, рассчитанные по формуле (14), можно получить хорошее совпадение.
Заключение
Таким образом, с достаточной степенью точности можно рекомендовать разработанную аналитическую модель и рассчитанную по ней методику для определения тепловых нагрузок на гелиоприемные коллекторы типа «горячий ящик», применяемые достаточно широко в водонагревательных установках для получения горячей воды. При этом достаточно знать значения солнечной радиации, время года и погодные условия.
СПИСОК ЛИТЕРАТУРЫ
1. Маэда Я., Аман Т., Танака Н. Характеристики плоских коллекторов солнечной энергии // Миядзаки дайгаку ногапугу кэнкю хококу. - 1980. - № 26. - С. 15-43. ВЦП НТЛ. Перевод № КГ-75273 (Киевская редакция).
2. Петухов Б. В. Метод расчета солнечных водонагревателей // Использование солнечной энергии. -М.: Изд-во АН СССР, 1957. - Сб. 1. - С. 177-201.
3. Руденко М. Ф., Черкасов В. И. Методика расчета гелионагревательных установок // Вестн. Астрахан. гос. техн. ун-та. Механика. - 1998. - С. 96-101.
4. Руденко М. Ф., Некрасов В. П., Черкасов В. И. Исследование эффективности высокотемпературных коллекторов // Вестн. Астрахан. гос. техн. ун-та. - 1996. - № 2/95. - С. 169-173.
5. Холодильная техника. Кондиционирование воздуха. Свойства веществ: Справ. / Под ред. С. Н. Богданова. - СПб.: СПбГАХПТ, 1999. - 320 с.
6. Исаченко В. П., Осипова В. А., Сукомел А. С. Теплопередача. - М.: Энергоиздат, 1981. - 418 с.
Статья поступила в редакцию 26.11.2007
MODELLING OF ENERGY BALANCE AND CALCULATION OF SOLAR POWERED DEVICES FOR WATER HEATING
M. F. Rudenko, B. Zh. Turklenbaeva
A physico-mathematical model of energy balance of solar powered devices as " a hot box" is developed for water heating by the energy of solar radiation. On the basis of this model the technique and analytical dependences for the calculation of solar powered devices is resulted. There is a calculated program developed in MathCad. The theoretical model is confirmed by natural experimental researches. It is possible to recommend the received results for the definition of heat loads on solar powered devices as "a hot box".
