CC BY
166
УДК 629.424.3:621.43.016
Я. К. СКЛИФУС, В. И. МОГИЛА (ВНУ им. В. Даля, Луганск)
ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТА ТЕПЛООТДАЧИ ПРИ КОНДЕНСАЦИИ ПАРА ВНУТРИ ТРУБ ТЕПЛОВОЗНЫХ РАДИАТОРНЫХ СЕКЦИЙ
В статье представлены расчет коэффициента теплоотдачи при конденсации пара, математическая модель распределения температуры в газовой и жидкой средах теплоносителя, модель формирования пленки конденсата на стенках труб.
Ключевые слова: конденсация пара внутри труб, коэффициент теплоотдачи, тепловой поток
В условиях развивающегося топливно-энергетического кризиса особое внимание на железнодорожном транспорте уделяется экономичности дизеля и его систем. Многолетний опыт судостроительных заводов доказал, что наиболее эффективными являются системы охлаждения (СО) дизелей, использующие фазовые переходы [1]. Уровень точности современного оборудования позволяет создавать подобные миниатюрные СО даже в масштабах легковых автомобилей, что с успехом доказано заводами BMW [2]. Очевидно, что СО с фазовыми переходами являются перспективным и могут быть применены в локо-мотивостроении [3]. При этом необходимо обеспечить замкнутый цикл, т. е. конденсацию пара, а для этого вполне реально использовать в качестве конденсаторов стандартные тепловозные радиаторные секции [4].
Количественной характеристикой теплопередачи в радиаторных секциях является коэффициент теплопередачи [5]:
K =
1
(-1+^). F+Л
а1 К F1 а 2
(1)
где К - коэффициент теплопередачи;
а1 - коэффициент теплоотдачи от жидкости к внутренней стенке; § - толщина стенки;
XМ - коэффициент теплопроводности стенки; , - площади соответственно внутренней и наружной поверхностей;
а2 - коэффициент теплоотдачи от наружной стенки к воздуху.
Известно, что а1 для воды составляет
4 650...6 400 Вт/(м2К), когда а2 для воздуха составляет 58. 175 Вт/(м2К). «Плоская» форма трубок радиатора является наиболее эффективной и экономичной с аэродинамической точки
зрения. При использовании радиаторных секций в качестве конденсаторов, необходимо рассчитывать а1 для конденсации пара внутри
«плоских» труб.
В данной работе рассмотрена задача теплообмена при пленочной конденсации насыщенного водяного пара в вертикальной трубе круглого и «плоского» (рис. 1) сечения. Вход пара сверху -направление движения пара совпадает с направлением силы тяжести и ускоряет движение пленки конденсата вследствие сил межфазного трения, что уменьшает толщину пленки и повышает коэффициент теплоотдачи. Режим процесса установившийся. Насыщенный пар с температурой ^ (равной температуре конденсации) движется со средней скоростью входа 0 в трубе радиуса Я . Стенка трубы имеет постоянную температуру < ts по всей длине.
Рис. 1. Сечение «плоской» трубки радиатора тепловоза: 01,02 - центры полуокружностей;
Я - радиусы полуокружностей; п - отношение радиуса к половине ширины трубы
Теплообмену внутри труб посвящены многие труды известных ученых, таких как В. Нус-сельт, В. П. Исаченко, Х. Хартманн, С. С. Кута-теладзе, Ф. Крейт, У. Блэк и другие. Однако в данных математических моделях рассмотрены лишь трубы круглого сечения, а также зачастую ряд теоретических расчетов заменен экс-
© Я. К. Склифус, В. И. Могила, 2012
периментальными коэффициентами. Представленные в данной статье вычисления базируются на трудах перечисленных ученых.
Построена математическая модель распределения теплового потока в двух средах: 1) газовая фаза (водяной пар в идеально-газовом состоянии), 2) жидкая фаза (пленка конденсата, вода). Формула (1) объединяет теплоотдачу а1 и а2 без расчета распределения теплового потока внутри стенки трубы. Очевидно, что температура теплоносителей является основным параметром, влияющим на величину теплоотдачи. Интерес представляет расчет теплопередачи на границе раздела фаз и в пристеночном слое конденсата. Характеристики слоя тепло-передающего вещества и температура на границе рассматривались как начальные условия для решения уравнения Фурье [6].
Основное дифференциальное уравнение распределения температур в теплоносителях имеет вид:
д 2Т 1 дТ д 2Т
дг2 г дг dZ
_ 1 дТ 2 _ D дЛ
-— + —2 _-—, (2)
д2u 1 ди д2u г0 ди
д р2 р д р д2 D дЛ
(3)
хт_ <г<г0
дZ |0, г > г0
при Z _ 0
при z _ 0
и и ^ 0 при -y/p2
Применив к дифференциальному уравнению преобразование Ханкеля по координате р, затем применив преобразование Лапласа по времени, с учетом начального условия и(^ = 0) = 0 , получили выражение, определяющее температуру поверхности при переходных процессах:
и _
J0(а р) Ji(a)
а
erf(—\[Dt )dа .
(4)
где а - переменная преобразования Ханкеля; 1 п
Jm (y) _ — jcos(mt - y sin(t))dt - интегральная n 0
формула Бесселя для действительного аргумента;
2 у
m _ 0,1,2,...; а erf(y) _-/=[exp(-t2)dt - функ-
л/я 0
ция ошибок.
Возвращаясь к исходным переменным, получили распределение температуры по координате y по времени (5):
где В - коэффициент температуропроводности; Т - температура теплоносителя; ^ - время;
X - осевая координата; г - радиальная координата. Далее, положив р = г/т0, г = Х/г0,и = X Т/дг0 , X - коэффициент теплопроводности; q - плотность теплового потока; г0 - половина линейного размера источника тепла (элементарного точечного параллелепипеда); уравнение (2) перепишем в виде:
Т _qTLi ТЫ) _ х J
Чг0 fJ0(Ур)J1(У)
y
erf
V г0
Решение исходного уравнения должно удовлетворять граничным условиям
и условию на бесконечности Т ^ 0 при
2 + X2 ^ да . Для переменной и эти условия записываются следующим образом:
ди = Г—1,0 <р<1 дг [0, р>1
■ да.
Интегралы определялись численно, методами Рунге-Кута и Симпсона. Функции Бесселя нулевого и первого порядка рассчитывались по интегральным формулам [7].
В результате численного решения уравнения Фурье были получены значения температур в каждой точке парового потока и пленки конденсата, и, соответственно, градиенты температур, необходимые для построения математической модели местного коэффициента теплоотдачи по длины трубы.
Плотность теплового потока пропорциональна градиенту температуры согласно закону Фурье [6]:
q = —Х-8гаа(Т), (6)
где q - вектор плотности теплового потока;
X - коэффициент теплопроводности.
Для решения второй задачи (распределения теплового потока в жидкой фазе и определение коэффициента теплоотдачи от жидкости к стенке) необходимо также рассчитать толщину пленки конденсата.
Расчет толщины пленки проводился исходя из баланса массовых расходов для элементарного кольца слоя конденсата с толщиной . Если (рис. 2) обозначить проникающую в элементарное кольцо массу конденсата тк1, выте-
со
0
кающую массу уже жидкого конденсата тк2 и конденсирующуюся в свободной поверхности массу тк3, то [8]:
тк1 + тк3 = тк2,
У
(7)
у=§( ^) У=0
ткх = 2прк | (Я - у) • Мк (у)—у, (8)
(9) (10)
тк2 = тк1 + —т—1 —2
тк3 = 2 п———— (Я - 5—, г • 5
Для участка (а) массовые расходы определяются формулами (8), (9), (10).
Для участка (б) массовые расходы рассчитывались по следующим формулам:
§( ^ )
ткА = рк | (Я - у) • м— (у— • 4(п -1)Я , (11)
(12)
тк2 = тк1 + —т—1 ёг
тк3 = (4(п - 1)Я— , (13)
г • 5
Из формул (8), (9), (10) для круглого сечения трубы следует отношение [8]:
—к -А— Я - 5 г Р— 5
У=§( г )
| (Я - у) м— (У)—у , (14)
У=0
Рис. 2. Вертикальное осевое сечение трубки с указанием баланса массовых расходов:
тк - массовые расходы; Я - радиус; —г - элементарное кольцо слоя конденсата; 5 - толщина пленки конденсата
Для отдельных участвующих массовых расходов для круглого сечения (см. рис. 2) исходя
из формуЛ т— = Рк • М • ^сеч ; тк3 • г = Ч • Sпов = =—7А— ^ "пов (где "сеч - пл°щадь сечения; ^пов
§( г)
- площадь поверхности конденсации; г -удельная теплота испарения; 5 - толщина пленки конденсата) рассчитаны отношения [8]:
Для «плоского» сечения трубы, для участка (а) использовалась формула (14), для (б) получена формула:
——•А— —
5( г )
— • | (Я - у)м— (у)—у , (15)
г Р—5 0
Еще неизвестное распределение скорости мк (у) было представлено в виде кубической параболы [8]:
м— ( у) = а + ьУ + сУ 2 + —У3,
(16)
Коэффициенты а, Ь, с, — определяются граничными условиями [8]:
где рк - плотность конденсата; у - расстояние от стенки; мк (у) - скорость конденсата; г - удельная теплота испарения. Для «плоского» сечения трубы элементарную площадку разобьем на две составляющие (см. рис. 1):
(а) две полуокружности
х1 е (-пЯ; - (п -1)Я) и ((п -1)Я; пЯ);
(б) прямоугольник х1 е (-(п - 1)Я; +(п - 1)Я), у1е (-Я; Я).
у = 0 у = 0
у = 0
у = 5
м— (у) = 0,
' дм >
ду
V ^ /0
П—
(17)
(18)
'дЧ ^
д2 у
1
0
'дм ^
ду
= —\—г-8Рк , (19)
Пк V )
Пк
(20)
\ ^ л
где тМ - напряжение на границе конденсата и стенки,
т5 - напряжение на поверхности пленки
вследствие трения пара, —р
—г
ния пара:
падение давления в направлении тече-
р— м
2Я 2
5
2
где С = 0,184 Яе-0 2;
Яе^ - критерий Рейнольдса для пара.
Согласно равенствам mD0 = • wd0 • £сеч и mD = mD0 - mk1 (где mD0 и mD - массовые расходы пара в начальном и текущем сечении) из формул (8), (9), (10) для круглого сечения следует отношение [8]:
2 Pk х
R2 w
_do _
пг! ~-т---т-
d (Я - 5)2 (R - 5)2 Pd
у=5( г)
х | (R - у^ (у)ф
у=0
(22)
^ 0 R
Pk
5( г )
| (R - у)Wk (у)Оу , (23)
ному количеству теплоты:
г • тк3
- Ы • ckmk2
^ Ы
(25)
где ск - теплоемкость конденсата; ^ - площадь поверхности трубы. Средний по длине трубы коэффициент теплоотдачи был определен согласно формуле:
а\ср ="
(26)
где г - длина трубы;
к - шаг интегрирования.
Все вычисления проводились в автоматическом режиме в специально написанных программах.
По результатам расчетов было проведено сравнение а1ср (рис. 3), полученных для круглых и «плоских» труб, с усредненными результатами [9], полученными для круглых труб по расчетным зависимостям других ученых: Хартмана Х. [8], Исаченко В. П. [10], Уонга Х. [11], Михеева М. А. и Михеевой И. М. [12].
Для «плоского» сечения трубы, для участка (а) использовалась формула (22), для (б) получена формула
1,5
(Я - 5) Pd (Я - 5) 0 В итоге для круглой трубы была получена система уравнений (14) и (22), а для «плоской» трубы - система из четырех уравнений (14), (15), (22), (23); в обоих случаях wk (у) определялось уравнением (16). Далее для определения толщины пленки конденсата 5(г) были численно решены полученные системы уравнений с начальными условиями (17)...(21) а также (24) г = 0 5(0) = 0
г = 0 Wd (0) = Wd 0, (24)
Решение уравнений (14), (15) проводилось методом Рунге-Кутта; интегралы (22) и (23) были решены методом Симпсона. Далее для заданного шага интегрирования была получена таблица рассчитанных значений толщины пленки конденсата. После аппроксимации полученных значений кубической параболой (и сплайном для точности) - была получена модель образования пленки конденсата по длине трубы для круглого и «плоского» сечения.
После вычисления массовых расходов, распределения температуры в средах и толщины пленки конденсата, был определен местный коэффициент теплоотдачи а1 согласно передан-
Ц5
аЖ5, Вт/(п2К)
/ УУ' у у 2 Г
1
3
к'е-10''
Рис. 3. График зависимости среднего коэффициента теплоотдачи а1ср при конденсации пара внутри
труб от критерия Рейнольдса пара Я^ :
1 - усредненные результаты для круглых труб, полученными по расчетным зависимостям других ученых (искусственно наложены на результаты программ для наглядности); 2 - результаты вычисления для трубы круглого сечения; 3 - результаты вычисления для трубы «плоского» сечения
В расчетах было принято значение внутреннего диаметра трубы d = 0,01 м (10 мм). Температура конденсации ts = 100 °С, температура стенки tw = 98 °С. Значение массового расхода пара 00 на входе в трубу было задано в диапазоне 1-10"5...0,05 кг/с, что соответствует в данном случае значению числа Рейнольдса Я^ =106,5...5,3-105 для пара. Представленные на графике (см. рис. 3) результаты расчетов для трубы «плоского» сечения соответствуют трубе с аналогичным Я^, но со значениями Я = 4 мм и п = 10 . При этих условиях (Я > 3,8 мм) [13] отсутствует перераспределение пленки конденсата по поверхности трубы под воздействием сил поверхностного натяжения. При этом значения а1 , полученные для круглых труб, должны совпадать с а1 для «плоских», что явно прослеживается на графике.
В результате применения принципов конфигурационного управления процессом теплоотдачи
удалось получить зависимость коэффициента теплоотдачи от входных параметров для различных сечений трубы. Численное решение уравнения Фурье позволяет моделировать тепловой поток аппроксимирующими многочленами адекватно реальному. Численные результаты расчета теплоотдачи для трубы круглого и «плоского» сечения, полученные нами, совпадают с результатами других исследователей с достаточно малой интегральной погрешностью. Это подтверждает правильность выбранного подхода.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Селиверстов, В. М. Утилизация тепла в судовых дизельных установках [Текст]/ В. М. Селиверстов. - Л. : Судостроение, 1973. - 256 с.
2. Синтезгаз. Паровой двигатель на BMW 5-ой серии [Электрон. ресурс] - 2011. - Режим доступа: http://sintezgaz.org.ua/energonovosti /173/parovoi-dvigatel-na-bmw-5 -i-serii
3. Могила, В. И. The prospects of increasing the effectiveness оf the cooling device of a diesel locomotive [Текст] / В. И. Могила, Я. К. Склифус // TEKA Commission of Motorization and Power industry in agriculture. - Lublin, 2010. - Т. X C. - С. 198-203.
4. Могила, В. И. The cooling devise of locomotive with vaporizing coolant [Текст] /
B. И. Могила, Н. И. Горбунов, Я. К. Склифус // TEKA Commission of Motorization and Power industry in agriculture. - Lublin, 2011. - Т. XI A. -
C. 169-176.
5. Камаев, А. А. Конструкция, расчет и проектирование локомотивов [Текст]: учебник для вузов /
Я. К СКЛ1ФУС, В. I. МОГИЛА
А. А. Камаев, Н. Г. Апанович, В. А. Камаев и др. / под ред. А. А. Камаева. - М. : Машиностроение, 1981. - 351 с.
6. Куликов, Ю. А. Системы охлаждения силовых установок тепловозов [Текст] / Ю. А. Куликов. - М., Машиностроение, 1988. - 280 с.
7. Корн, Г. Справочник по математике (для научных работников и инженеров) [Текст] / Г. Корн, Т. Корн. - М., 1974. - 832 с.
8. Hartmann, H. Wärmeübergang bei der Kondensation strömender Sattdämpfein senkrechten Rohren [Текст] / H. Hartmann // «Chemie - Ingeneur -Technic». - 1961. Bd. 33, volume 5. - 521 с.
9. Могила, В. И. Анализ погрешности расчетов коэффициента теплоотдачи при конденсации пара в трубах [Текст] / В. И. Могила, Я. К. Склифус // Вюник Схвдноукраш. нац. ун-ту iм. В. Даля. -Луганск, 2011. - Ч. 1. - № 4 (158). - С. 136-142.
10. Исаченко, В. П. Теплообмен при конденсации [Текст] / В. П. Исаченко. - М. : Энергия, 1977. -240 с.
11. Уонг, Х. Основные формулы и данные по теплообмену для инженеров [Текст] : Справочник / Х. Уонг ; пер. с англ.- М. : Атомиздат, 1979. -216 с.
12. Михеев М. А., Основы теплопередачи [Текст]/ М. А. Михеев, И. М. Михеева. Изд. 2-е. - М., «Энергия», 1977. - 344 с.
13. Химическая энциклопедия [Текст] : в 5 т. : т. 2 : Даффа-Меди / редкол. : И. Л. Кнунянц (гл. ред.) и др. - М. : Сов. энцикл., 1990. - 671 с.
Поступила в редколлегию 20.03.2012. Принята к печати 20.04.2012.
ВИЗНАЧЕННЯ КОЕФЩКНТА ТЕПЛОВ1ДДАЧ1 ПРИ КОНДЕНСАЦП ПАРИ ВСЕРЕДИН1 ТРУБ ТЕПЛОВОЗНИХ РАД1АТОРНИХ СЕКЦ1Й
У статп представлеш розрахунок коефщента тепловiддачi при конденсацп пари, математична модель розпод^ температури в газовш i рщкш середовищах теплоноая, модель формування плiвки конденсату на стшках труб.
Ключовi слова: конденсацiя пари усередиш труб, коефщент тепловщдач^ тепловий попк
Y. K. SKLIFUS, V. I. MOGILA
DETERMINATION OF HEAT TRANSFER COEFFICIENT WITH VAPOR CONDENSATION INSIDE THE TUBES DIESEL'S RADIATOR SECTIONS
The article presents the calculation of heat transfer coefficient during condensation of steam, the mathematical model of temperature distribution in the gas and liquid phases of the coolant and the model of the formation of the condensate film on the walls of the tubes.
Keywords: condensation of steam inside the pipes, the heat transfer coefficient, heat flux
