Методика расчета тепловой и аэродинамической характеристик калориферов из биметаллических труб со спиральными накатными и навитыми алюминиевыми ребрами Текст научной статьи по специальности «Строительство и архитектура»

УДК 536.24+697.1

Методика расчета тепловой

и

и аэродинамическои характеристик калориферов из биметаллических труб со спиральными накатными и навитыми алюминиевыми ребрами

Канд. техн. наук А.Н.БЕССОННЫЙ ООО “НИИхиммаш”, д-р техн. наук, профессор В.Б.КУНТЫШ, Т.П.КАЛЕЙЧИК Белорусский государственный технологический университет

The most general methodic of heat transfer calculation for industrial bimetallic heaters consisting of pipes with spiral aluminium ribs is described. Methodic is based on distinct accounting of thermal resistances to heat transfer from heating warming medium to heated transverse air flow by frequently tested experimental dependences. The heat transfer of pipes and its aerodynamical resistance is suggested to calculate by generalized criteria I equations. These equations embrace wide interval of ribs and pipe parameters, arrangemental characteristics of bundle regime of air motion.

В последнее время получили широкое применение стандартизированные поверхностные рекуперативные воздухонагреватели — калориферы из биметаллических ребристых труб (БРТ). Они применяются в установках вентиляции и кондиционирования воздуха, в отопительно-вентиляционных системах зданий, для воздушного отопления производственных помещений, в сушильных установках различных технологических процессов, в системах утилизации тепла удаляемого воздуха промышленных цехов и общественных зданий, в системах предварительного подогрева первичного воздуха энергетических котлов.

Промышленные серийные калориферы конкретного ряда состоят из собранных в шахматные пучки БРТ одного типоразмера с постоянными значениями поперечного 5, и продольного S2 шагов размещения труб в решетках при числе поперечных (относительно направления воздуха) рядов z = 2, 3 и 4. Различные значения площади поверхности нагрева калорифера достигаются применением в поперечном ряду труб различной длины /0 при постоянном их количестве. В отдельных конструктивных исполнениях (калориферы значительной тепловой мощности) площадь поверхности нагрева изменяется путем увеличения числа труб в поперечном ряду при l0 = const. Следовательно, перемен-

ной является площадь поперечного сечения калорифера /,.

БРТ представляет собой несущую круглую трубу из углеродистой, нержавеющей стали, латуни, меди со спиральными накатными алюминиевыми ребрами. Перспективно применение БРТ со спиральными навитыми из алюминиевой ленты ребрами, для которых характерны в 1,5...2 раза меньший расход алюминия, в 2,5...3 раза более низкие затраты электроэнергии на оребрение трубы [7]. Кроме того, в технологическом процессе оребрения не нужна смазочно-охлаждающая жидкость, загрязняющая окружающую среду.

В применяющихся методиках расчета и подбора стандартизированных калориферов из БРТ коэффициент теплопередачи к [Вт/(м2-К)] принимают либо по табличным данным [9] в зависимости от массовой скорости воздуха м>р [кг/(м2-с)] (м> - скорость воздуха; р - его плотность), либо по скорости греющей воды V, или рассчитывают по эмпирическим формулам функционального вида [1,2, 8], которые получены по результатам испытаний головных образцов:

£=/(и>р, V); (1)

*=/(м;р,/0); (2)

к =f(wp,f^), (3)

где 10 - длина трубы.

Потери давления воздуха на пучке труб калорифера вычисляют по уравнению Ар =/И>). (4)

Формулы (1) - (4) имеют частный характер: область их применения ограничена конкретными значениями геометрических параметров ребер и трубы, их материалом, конструкцией ребра, формой поперечного сечения трубы, углом атаки воздушного потока компоновочными характеристиками труб в пучке (5р 52, г). Распространение формул на иные типоразмеры БРТ, характеристики пучка и модульные ряды калориферов невозможно.

Значения коэффициентов теплопередачи, вычисляемые по (1) - (3), включают составляющие термического сопротивления в процессе теплопередачи в скрытой форме, поэтому нельзя выявить влияние конкретного способа теплообмена на интенсификацию передачи тепла от греющего теплоносителя к нагреваемому воздуху. Представление теплопередающей характеристики калорифера формулами указанного вида крайне затрудняет научно обоснованный выбор направления исследований по повышению энергетической эффективности и оптимизации конструкторско-компоновочных характеристик пучка труб калорифера и геометрических параметров БРТ.

В настоящее время промышленность производит много различных типоразмеров БРТ. Однако методика расчета коэффициента теплопередачи калорифера, состоящего из труб, параметры которых отличаются от параметров труб серийных биметаллических калориферов, отсутствует. Тем самым априори исключается возможность проектирования и изготовления индивидуальных калориферов, наиболее полно удовлетворяющих требованиям конкретного потребителя. В ряде случаев возникает ситуация, при которой допустимая потеря давления воздуха в калорифере не может быть обеспечена серийным аппаратом из модульного ряда. Тогда необходим калорифер с иными значениями 5, или г, так как слабо влияет [7] на аэродинамическое сопротивление пучка.

Цель работы — создание универсальной методики расчета коэффициента теплопередачи и аэродинамических характеристик биметаллических калориферов из шахматных пучков БРТ с накатными и навитыми алюминиевыми ребрами. Методика яв-

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

Сущность методики заключается в том, что коэффициент теплопередачи БРТ К [Вт/(м2-К)], отнесенный к полной наружной площади поверхности оребрения трубы, вычисляется по аналитической формуле Я=(Д1+Я2+Лз+Я4+Л5)->, (5)

о _ 1 Ф^О

где —----------- приведенное термическое сопротив-

ление при теплоотдаче с внутренней стороны трубы, м2К/Вт;

*2 =

Ні 4 К

- приведенное термическое сопротив-

ление стенки несущей трубы, м2К/Вт;

В - V

л3 - Кк------- приведенное термическое контакт-

ное сопротивление трубы, м2-К/Вт;

Я - термическое контактное сопротивление, м2К/Вт;

и _ <К 5а

л4 —--------— - приведенное термическое сопротив-

(1 А

к а

ление стенки алюминиевой ребристой оболочки, м2-К/Вт;

И5 — — - приведенное термическое сопротивление а

при теплоотдаче от оребрения к воздуху, м2К/Вт; ф — коэффициент оребрения трубы;

а,, а — соответственно коэффициент теплоотдачи от греющего теплоносителя к внутренней поверхности трубы и приведенный коэффициент теплоотдачи от оребрения к воздуху, Вт/(м2-К); с1Г £?к, й? — соответственно внутренний диаметр трубы; диаметр контактной зоны, равный наружному диаметру несущей трубы; диаметр трубы по основанию ребер, м;

5 , 5а - толщина стенки несущей трубы и алюминиевой оболочки соответственно, м;

X , Аа -коэффициент теплопроводности стенки несущей трубы и ребристой оболочки соответственно, Вт/(м-К).

Расчет численных значений Л,, У?2, Л4 не вызывает затруднений. Коэффициент теплоотдачи от греющей воды, движущейся внутри прямых гладких труб кало-

а, = 0,021 —Ке°'8 Рті0’4,

(6)

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

4

где Я, - коэффициент теплопроводности воды, Вт/(м-К); Ке1 = V Д/у, — число Рейнольдса;

V, — коэффициент кинематической вязкости воды,

м2/с.

РГд - число Прандтля для воды.

Значения X,, V,, Рг1 принимаются по средней температуре воды в калорифере.

При использовании в калорифере конденсирующегося внутри прямых гладких труб водяного пара широко применяют следующую зависимость [7]:

а, =(3955 + 116м/п)3/1,21 /10, (7)

где IV — скорость пара при входе в трубу, м/с.

Совершенно очевидно, что аналитическая зависимость (5) имеет гораздо более общий характер, чем (1) - (3), однако достоверный расчет коэффициента теплопередачи возможен лишь при надежном расчете термических сопротивлений /?5 и /?к, то есть при наличии обобщенного критериального уравнения для коэффициента теплоотдачи воздуха и значения термического контактного сопротивления Лк в диапазоне температур эксплуатации БРТ - калориферов и материалов контактной зоны.

Анализ литературных данных показал, что для расчета конвективного коэффициента теплоотдачи воздуха, омывающего шахматные пучки из круглых труб со спиральными накатными и шайбовыми ребрами или спиральными навитыми ребрами из ленты Ь-образного поперечного сечения, наиболее обоснованно применение критериального уравнения [7], базой для которого явились экспериментальные исследования [4] более 100 пучков из труб с алюминиевыми ребрами различных геометрических и компоновочных параметров, включающих всю размерную гамму БРТ стандартных калориферов [3]:

ак =0,132—С С

5

0 43 1 С/Г т / і \ “о -0,54 Г лЛ

5, -<іп 5 5

V °) Ч / ч /

Ее'

0,73

(8)

где А, - коэффициент теплопроводности воздуха, Вт/(м-К);

/г, 5 -высота и шаг ребра;

С, =Дг)—поправочный коэффициент на число поперечных рядов в пучке, принимаемый равным для

ъ = 2, 3,4, 5, 6 и более —0,91; 0,98; 0,99; 0,995; 1,0; С¥ = Д\|/) — поправочный коэффициент на угол атаки \|/ потоком воздуха пучка БРТ, принимаемый равным для ці = 90°, 80°, 70°, 60° -1,0; 1,035; 1,037; 1,167.

т = 0,53 — 0,019ф - показатель степени;

5=а/(°’551 )2 + —диагональный шаг труб;

Не = \vsfv — число Рейнольдса для воздуха; н- — скорость воздуха в наименьшем проходном сечении пучка, м/с;

у — коэффициент кинематической вязкости воздуха, м2/с;

Графическая зависимость =Д\|/) [6] аппроксимирована выражением С, = -1664,83 + 132,41Л|#- - 4,36\|/2+ 0,076\у3 -

- 7,44ч/4+3,86у5. (9)

Физические параметры воздуха А..У принимают при его средней температуре в пучке калорифера.

Погрешность расчета ак по (8) не превышает ±15 %. Уравнение обобщает опытные данные в диапазонах: Ие = 150...85000; Щз - 2,88...16,1; Ш = 0,4...6,4; (5,-^)/(52-й0) = 0,46...2,61; <р = 5...22.

Связь приведенного коэффициента теплоотдачи а с конвективным коэффициентом теплоотдачи «к устанавливается выражением (10) [11]

а = ак[(Гр/Г)иЕц^Ргр/Ги (10)

где Р = Р'6 + /^—площадь поверхности ребер на 1м длины трубы;

/г,—площадь боковой поверхности ребер на 1 м длины трубы;

площадь поверхности торцов ребер на 1м длины трубы;

^гр- площадь поверхности трубы длиной 1м по основанию ребер, не занятая последними;

Р = + Р —полная площадь поверхности оребре-

ния 1м длины трубы;

ц — коэффициент формы поперечного сечения ребра. Для ребер прямоугольного и близкого к нему сечения ц = 1,0.

Е = ДРк,Шс1^ — коэффициент эффективности ребра, принимаемый по графикам, приведенным, например, в [7, 11]. Для алюминиевых ребер труб стандартных калориферов Е = 0,94.. .0,98.

<1 - наружный диаметр оребрения, м;

(3/г = /г ^2ак /(ХрА) — безразмерный комплекс;

Хр — коэффициент теплопроводности материала ребра, Вт/(мК);

Д — средняя толщина ребра, м;

\|/р- коэффициент неравномерности распределения коэффициента теплоотдачи по поверхности ребра

[И],

vj/p= 1 -0,058р/г. (11)

Значения Rk для биметаллических стальных труб и алюминиевых накатных ребер в интервале изменения средней температуры зоны контакта = 50...230 °С с погрешностью ±10 % вычисляют по [7, 4] как

Rk=0,22-10_3 + 2,5-10-6 (tK - 95). (12)

Можно принимать (при tK< 90 °С) для труб из углеродистых сталей /?к=1,83-104 м2-К/Вт, для нержавеющей стали — (3,2...3,7)ТО^ м2-К/Вт, для латуни ЛОМш-70-1 -

0,7-Ю"4 м2К/Вт. '

Для БРТ с навитыми алюминиевыми L-образными ребрами значения /?к приводятся в работе [5]. Такие трубы в настоящее время по различным причинам еще не нашли широкого применения.

Для расчета потерь давления воздуха Ар (Па) на шахматном пучке калорифера рекомендуется [7] применять обобщенное уравнение [11], которое по исследованиям [4] может быть распространено на БРТ со спиральными накатными и навитыми ребрами:

Ар = 2,7zCCJl/d)™ Re^25 pw2, (13)

где С, = f (z) — поправочный коэффициент на число

поперечных рядов, принимаемый равным для z -

1, 2, 3, 4, 5 и более - 1,23; 1,18; 1,08; 1,02; 1,0; для труб с навитыми ребрами С, = 1,0;

С¥ =Д\(/) — поправочный коэффициент на угол атаки воздухом труб пучка, принимаемый равным для \|/ = 90°, 80°, 70°, 60° соответственно 1,0; 1,079;

1,188; 1,60.

Графическая зависимость С¥=Д\|/) [6] нами аппроксимирована выражением

С¥=-143,016 + 11,04\|/ - 0,323\(/2ч- 0,0046\|/3-

— 3,17 8\|/4 -ь 8,629\|/6, (14)

lld= 0,15... 12,5; у <90...60°.

Линейный определяющий размер, м

I ~ (Ftp/F)d0 + (Fp/F)^OJS5(d2-d20). (15)

Эквивалентный диаметр наименьшего проходного

сечения пучка, если таковым является фронтальное поперечное сечение (как во всех стандартных калориферах), рассчитывают по формуле:

= 2ВД - й0) - 2/гД]/(2/г + 5). (16)

Если наименьшим проходным сечением является диагональное, то в (16) вместо 5, подставляется 52.

Уравнение действительно в интервале Яе, = \vllv = = (2...180)Т03.

Физические свойства воздуха р,у принимаются по его средней температуре в пучке.

Погрешность вычисления Ар составляет ±20 %.

***

С целью проверки применимости изложенной методики нами выполнены сравнительные расчеты коэффициента теплопередачи и потерь давления воздуха при массовых скоростях воздуха во фронтальном сечении ыр -= 3 и 7 кг/(м2-с) для серийных паровых калориферов № 6 и 10 Костромского завода с трехрядным и четырехрядным шахматным пучком. Конструктивные размеры калориферов даны в [3, 8]. Средняя температура водяного пара в пучке принята 50 °С, давление сухого насыщенного воздуха 0,12 МПа. Пучок труб состоит из БРТ с накатными алюминиевыми ребрами размерами dXd0XhXsXA = 39,0 х 20,0 х 9,5 X 3,4 X 0,825 мм; ф = 9,5. Разбивка труб в решетках равносторонняя с шагами 5, = Б2'= 41,5 мм; 52 = 36 мм. Несущая труба с наружным диаметром dк = 16 мм и толщиной стенки 5ст =1,2 мм выполнена из углеродистой стали. Скорость пара на входе в трубы калориферов >уп= 4 м/с. Результаты расчетов сведены в таблицу.

Характе- ристика Массовая скорость wp, кг/(м2с) №6 № 10

Число поперечных рядов в калорифере г

3 4 3 4

К Вт/(м2-К) 3,0 49,3/49,9 45,0/50,4 46,8/49,9 43,8/50,4

7,0 69,5/79,3 67,4/80,1 66,0/79,3 65,5/80,1

Др, Па 3,0 45,3/68,9 57,7/90,2 45,3/68,9 57,7/90,2

7,0 213,0/277,7 244,5/363,6 213,0/277,7 244,5/363,6

В таблице значения коэффициента теплопередачи в числителе вычислены по формуле (3) с использованием данных испытаний промышленных головных образцов [2, 9] (далее обозначается кп), а в знаменателе — по обобщеному уравнению (8) (далее кр). Значения потерь давления воздуха в числителе получены по зависимости (4) из [2, 8] (далее Дрп), а в знаменателе - по предлагаемому уравнению (13) (далее Дрр).

Расхождение результатов расчета по разным формулам вычисляли в % следующим образом:

5*= [(*,-*)/*„] 100 (17)

иб/7 = [(Дрп-Дор)/Дрп]100. (18)

Расчетные значения коэффициента теплопередачи по обобщенной формуле превышают данные традиционной методики, базирующейся на испытаниях промышленных головных образцов, на 1,2...22,3 %. Результат является приемлемым, так как находится в поле накопительной погрешности исходных уравнений в аналитической формуле (5). И хотя ни в одном из справочных изданий по калориферам и результатам их заводских испытаний не указывается погрешность формул (1) - (4), если судить по имеющимся результатам промышленных испытаний головных образцов регенераторов [11] газотурбинных установок, тем не менее погрешность их не ниже ±15 %. Также хорошо известно, что .расчетные уравнения и соотношения, полученные при экспериментальном исследовании моделей образцов теплообменных устройств, и построенные на их базе обобщенные критериальные уравнения являются наиболее достоверными и надежными по сравнению с полученными на базе иных опытных данных. Сказанное не противоречит необходимости проведения испытаний натурных образцов аппаратов, так как подобные испытания позволяют установить влияние концевых эффектов поверхности теплообмена, неравномерности скорости воздуха перед пучком, турбулентности набегающего потока и других факторов, обусловленных масштабным фактором.

Расчетные потери давления воздуха на пучке превышают полученные по частному уравнению (4) на 48,7 — 52 %. Результат удовлетворительный с учетом того, что и частные, и обобщенные зависимости для потерь давления имеют бблыную погрешность, чем тепловые характеристики.

Таким образом, наиболее целесообразными для расчета теплоаэродинамических характеристик биметаллических калориферов из БРТ являются аналитическая зависимость (5) и уравнение (13). Вместе с этим остается актуальной задача дальнейшего обобщения опытов по аэродинамическому сопротивлению моделей пучков БРТ и построению критериального уравнения большей точности.

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

1. Богословский В.Н., Поз М.Я. Теплофизика аппаратов утилизации тепла систем отопления, вентиляции и кондиционирования воздуха. - М.: Стройиздат, 1988.

2. Ефимов А.Л., Косенков В.И., Яковлев И.В. Системы кондиционирования воздуха. - М.: Изд-во МЭИ, 2002.

3. Каталог выпускаемой продукции. Костромской калориферный завод,2003.

4. Кунтыш В.Б., Кузнецов Н.М. Тепловой и аэродинамический расчеты оребренных теплообменников воздушного охлаждения. - СПб.: Энергоатомиз-дат,1992.

5. Кунтыш В.Б., Пиир А.Э. Контактный теплообмен в биметаллических трубах со спирально-навитыми алюминиевыми ребрами Ь-образного поперечного сечения // XIII Школа-семинар молодых ученых и специалистов под рук. акад. РАН А.И.Леон-тьева “Физ. основы эксперимент, и мат. моделирования процессов газодинамики и тепломассообмена в энергетических установках”, 20 - 25 мая 2001г., Санкт-Петербург. Т. 2. - М.: Изд-во МЭИ, 2001.

6. Кунтыш В.Б., Федотова Л.М. Влияние угла атаки воздушного потока на теплообмен и сопротивление шахматного пучка оребренных труб // Известия вузов. Энергетика. 1983. № 4.

7. Основы расчета и проектирования теплообменников воздушного охлаждения. Справочник / А.Н.Бессонный, Г.А.Дрейцер, В.Б.Кунтыш и др.; Под общ. ред. В.Б.Кунтыша, А.Н.Бессонного. - СПб.: Недра,1996.

8. Расчет, проектирование и реконструкция лесосушильных камер / Е.С.Богданов, В.И.Мелехов, В.Б.Кунтыш и др. - М.: Экология, 1993.

9. Рымкевич А.А., Халамейзер М.Б. Управление системами кондиционирования воздуха. - М.: Машиностроение, 1977.

10. Теплообменные аппараты из профильных листов / В.М.Антуфьев, Е.К.Гусев, В.В.Ивахненко и др. - Л.: Энергия, 1972.

11. Юдин В.Ф. Теплообмен поперечнооребренных труб. -Л.: Машиностроение, 1982.