Математическое исследование испарительного охлаждения теплообменного аппарата с пористыми ребрами Текст научной статьи по специальности «Механика и машиностроение»

УДК 621.573.001.2

Математическое исследование испарительного охлаждения теплообменного аппарата с пористыми ребрами

Канд. техн. наук, чл.-кор. МАХ В.В. КИРЕЕВ Ангарская государственная техническая академия

In the calculation of the evaporative cooling apparatus the method was used in which during integration of the equation for calculation of quantity of heat obtained by the air from film by convection and evaporation, the formula of Newton-Rihman and Dalton were applied. The obtained linear differential equations of first order can be used for determination of areas of refrigeration of apparatuses with porous fins at different temperatures of water film and parameters of atmospheric air or for evaluation of refrigerating effect at the pre-determined surface.

Теплообменники с ребристыми поверхностями применяются в тех случаях, когда теплообмен происходит между двумя теплоносителями с большим и малым коэффициентами теплоотдачи. Как правило, это система жидкость - воздух.

В результате патентного поиска и теоретических исследований выявлено, что применение каплеобразных труб в воздушных теплообменных аппаратах вместо труб круглого сечения позволяет увеличить эффективность теплообмена.

При впрыске воды в поток нагнетаемого воздуха часть влаги образует пленку на поверхности теплообмена, в результате испарения температура стенки понижается [1,3].

Другая часть мелкодисперсной влаги испаряется в потоке воздуха, который насыщается парами воды, дви-

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

Теплообменный аппарат (см. рисунок) состоит из корпуса 1, впускного коллектора 2, труб 4, выпускного коллектора 3, сетчатых пластин 5, перегородок коллекторов 7, межтрубного пространства б, форсунок 8.

Трубы в поперечном сечении выполнены в виде вытянутой капли, что позволяет увеличить поверхность теплообмена и уменьшить турбулизацию воздуха. Для улучшения теплообмена грубы расположены в шахмат-

Теплообменный аппарат (а) и схема движения охлаждающей среды (б)

ном порядке, причем выпуклой стороной они крепятся к сетчатой пластине, а противоположной стороной - к последующей сетчатой пластине 5.

Теплообменный аппарат работает следующим образом.

, Охлаждаемая среда поступает в впускной коллектор 2, затем под давлением поступает вверх по трубам 4 в выпускной коллектор 3. Благодаря перегородке 7 выпускного коллектора 3 среда направляется вниз по другому ряду труб 4 в впускной коллектор 2, затем перегородкой 7 впускного коллектора 2 направляется вверх по следующему ряду труб 4 в выпускной коллектор 3 и выводится из теплообменника.

Перед подачей воздуха в межтрубное пространство 6 теплообменного аппарата в струю воздуха нагнетается мелкодисперсный распыл воды форсунками 8. При впрыске мелкодисперсного распыла воды в поток нагнетаемого воздуха влага испаряется на поверхности теплообменного аппарата, понижая температуру стенки.

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

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

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

Для получения дифференциального уравнения изменения температуры ребра по высоте теплообменника и установления связи его с параметрами, определяю-

щими это изменение, выделим элементарный слой поверхности ребра единичной ширины с высотой йх. Для этого объема составим уравнение теплового баланса. Плотность теплового потока, проходящего через тепловыделяющую поверхность в пределах выделенного элемента qc, запишем в виде [2]: с1£)с _ Э?„

с1хЛ

а =--------— — ос я ——

1 рН'О и>

(1)

где 0С - количество теплоты; р№ — плотность воды; с - теплоемкость воды;

- среднеобъемная температура ребра по толщине в пределах выделенного объема.

Тепло, полученное от элемента теплоотдающей поверхности (1x1, расходуется на повышение температуры воздуха, а также на увеличение его влагосодержания.

Конвективная составляющая qa может быть выражена следующим равенством [2]:

Яа =

= Р*сгВ8*

сіх-1 ' ° рв в Эх • 1

Массообменная составляющая (2):

Яа =

СІОр

дс1„

(2)

(3)

с1х-1 в Эх-1 где г - удельная теплота парообразования;

1В и й?в - соответственно среднеобъемные температура и влагосодержание по толщине потока в пределах выделенного объема.

Полагая, что процесс тепло- и массоотдачи в пределах выделенного объема является установившимся во времени, примем, что

Яс=Яа+Яр- (4)

После подстановки в выражение (4) уравнений (1), (2), (3) получим уравнение теплового баланса в дифференциальной форме:

дт

' РвСРв#в

Э/„

- + Р„£в>--

э</„

(5)

" н’ Эх ■ 1 р Эх-1 ' '~вов Эх Г

После соответствующих преобразований получим

8С = Р*8„СРЛ* = Рв^А • (6)

Уравнение (6) совместно с уравнениями движения, сплошности, энергии и массообмена представляют систему дифференциальных уравнений, которые описывают скоростные, температурные и влажностные поля потоков распиливаемой воды и воздуха.

Для решения этой системы необходимо принять граничные условия [2] для водной аэрозоли и потока воздуха: скорость пленки, образующейся на пористой вертикальной поверхности \¥пл = 0, скорость воздуха у ее поверхности равна скорости поверхности пленки. Градиент скорости на внешней границе пограничного слоя

^ = 0. йу

Касательные напряжения сопротивления на границе раздела фаз равны между собой [1]:

dW„.

dy

= -Hb

dW.

>=®ш

dy

>=8„,

(7)

где |%,ЦВ - коэффициенты динамической вязкости воды и воздуха соответственно.

На границе стекающей пленки с твердой поверхнос-

dt . тью —— = 0, на границе раздела фаз температура

*У

пленки равна температуре соприкасающегося с ней слоя воздуха ? = / . Градиент температуры воздуха

на внешней границе равен

dt --2- = 0

dy

В качестве граничных условий принимается, что парциальное давление водяного пара максимально при температуре поверхности пленки. На внешней границе пограничного слоя градиент парциальных давлений

Фп

: 0 ИЛИ

dm

0.

dy dy

Вода поступает на поверхность со скоростью Wnjl = = const и температурой ?пл = const.

Поток воздуха также поступает в канал с постоянными скоростью WB = const, температурой tB = const и парциальным давлением пара рп , = const.

Количество испарившейся воды в элементе охладителя незначительно, поэтому расход воды на входе и выходе элемента gw= const.

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

На основании принятых допущений баланс тепла, отданного пленкой и полученного потоком воздуха,

(8)

можно выразить следующим уравнением:

= Рв£в[Срв (*.2 ~ *.1) + Г«2 - 4,1 )]•

В уравнении (8) теплофизические параметры воды и воздуха определяются при средних температурах этих сред по длине потоков.

При расчете аппаратов испарительного охлаждения распространены методы, в которых при интегрировании уравнения (5) для расчета количества тепла, отнимаемого от пленки воздухом, конвекцией и испарением, используются формулы Ньютона-Рихмана и Дальтона:

=<*(*» “О. (9)

где а - коэффициент теплоотдачи.

Яц =РрГ(Рп,с-Рп,о)’ (Ю)

где |3 - коэффициент массоотдачи, отнесенный к разности парциальных давлений.

При совместном решении уравнений (2) и (9) получаем:

Эг , ч

Рв^®Срв Эх4 = - (И)

При совместном решении уравнений (3) и (10) получаем:

dd a Рв^вГॳ" р пс_Рпо

(12)

(13)

Подставляя выражения (11) и (12) в уравнение (5), получим:

К .

РкСрк8”дхЛ

= а(*к-1в) + $рг(рпс-Рпо).

Линейные дифференциальные уравнения первого порядка (11), (12) и (13) могут быть использованы для определения площадей охлаждения аппаратов с пористыми ребрами при различных температурах пленки воды и параметрах атмосферного воздуха или для оценки охлаждающего эффекта аппарата при заданной поверхности.

Список литературы

1. Барановский Н.В., Коваленко Л.М., ЯстребенецкийА.Р. Пластинчатые и спиральные теплообменники. - М.: Машиностроение, 1973.

2. Исаченко В.П., Осипова В.А., Сукомел А.С. Теплопередача. -М.: Энергоиздат, 1981.

3. Ульянов Б.А., Бадеников В.Я., Щелкунов Б.И., Патрушев К.Ю. Расчет теплообменных аппаратов. - Иркутск: ИГТУ, 2001.