CC BY
39
УДК 662.987
А. А. Зуев, А. А. Кишкин, М. И. Толстопятов, Д. А. Жуйков
ТЕПЛООТДАЧА ВРАЩАТЕЛЬНЫХ ТЕЧЕНИЙ В ТУРБОМАШИНАХ НА ОСНОВЕ ДВУХСЛОЙНОЙ МОДЕЛИ ТУРБУЛЕНТНОГО ПОГРАНИЧНОГО СЛОЯ*
Рассмотрен частный случай течения газового потока, индуцированного вращением диска. Записаны интегральные соотношения уравнения энергии пространственного пограничного слоя для вращательного течения при степенном законе распределения скорости. Получены зависимости для оценки локальных коэффициентов теплоотдачи вращательных течений с различным распределением скорости.
Ключевые слова: теплоотдача, турбулентное течение, уравнение энергии, закон теплообмена.
Течение газового потока для турбонасосных агрегатов (ТНА) носит весьма сложный характер. Поперечный градиент давления оказывает существенное влияние на поверхностное трение и теплообмен. Влияние неучтенных в расчете тепловых нагрузок на узлы и конструкцию может привести к упругим деформациям, а также к нерасчетным параметрам вязкости и плотности рабочего тела, что влечет за собой отклонение от расчетного режима работы агрегата и падение КПД.
Широкое распространение в практике получили методы расчетов по учету теплоотдачи основанные на обобщении экспериментальных данных. Подобные методы показывают удовлетворительные результаты при рассмотрении относительно простых течений. В более сложных условиях, а именно при реализации вращательных течений, зависимости, основанные на обобщении экспериментальных исследований, не применимы, это объясняется достаточно скудными экспериментальными данными в этой области, а также отсутствием точных аналитических решений.
Большие возможности для поиска более точных зависимостей дает теория пограничного слоя. Многочисленными исследованиями установлено, что в большинстве важных для практики тепло- и массооб-мена между поверхностью тела и жидкостью процессов основная часть изменения температуры и концентрации также происходит в области, прилегающей к поверхности тела [1].
Рассмотрим частный случай турбулентного течения (Яеш > 3 105) газового потока около вращающегося диска: распределение скорости (и) по толщине пограничного слоя (5) к скорости в основном потоке (и), апроксимируется степенной функцией
и
и
(1)
Учтем, что можно выделить несколько зон турбулентного пограничного слоя. В работе [2] опытные данные по замеру профиля скорости соответствуют линейной зависимости
У
(2)
которая получена в предположении, что касательные напряжения полностью определяются как вязкостные, а динамическая скорость в выражении (2), характеризующая собой напряжение трения,
и * = 'т'
(3)
может быть сопоставлена с потерей напора ДА прямолинейного равномерного потока в цилиндрической трубе по формуле
аь й
то = анр 41
(4)
где АН - потеря напора по длине I; й - диаметр трубы.
Очевидно, что подход определения характеристик через значения напряжения трения весьма трудоемок с точки зрения аналитических выкладок. Выразим динамическую скорость через параметры пограничного слоя. Потребуем, чтобы на границе вязкого подслоя (5л) степенной закон (1) сопрягался с линейным
законом (2), т. е. при у = 5л будем иметь
(и * )’.
У
■ = и
, 8 ' (5)
V V 8.
Учитывая, что на внешней границе вязкого подслоя
и X
V
где ал - постоянная. Можно записать
5 =ОХ
и
(6)
Тогда с учетом (5), (6)
*
и а л =
и
1
а лу 17 и*5,
После несложных преобразований выражение для квадрата динамической скорости будет иметь вид
(и ‘ )2 = Цт
14
5и
(7)
V
*Исследование выполнено при поддержке Министерства образования и науки Российской Федерации, соглашение 14.В37.21.1835.
Рассмотрим случай, при котором толщины температурного и динамического пограничных слоев будут одинаковы, т. е. при Pr = 1
Т - Т 0 Т s- Т 0
и
U
(8)
S
=ЛУ
О*Ф = 1
0
1 ( У і7 О
ґ
T - To
Л
1 -T - T
V ts ToУ
dy =
dy = 0,097 -S.
(9)
Закон теплообмена в виде критерия Стантона для степенного профиля распределения скорости имеет вид
St = ■
X д
PCU дУ
Л
V ts To У
С учетом (8) выражение для определения толщины потери энергии в пограничном слое запишется в виде
0,558-X-
aл6 -v3 -S'
(13)
tv У
Ср •и1/4 • р
После упрощения закон теплообмена для турбулентного течения газа имеет вид
St =
X
0,558
1
П гг0,25 1,5 0,75 с **0,25 '
pCvU , «-л v , ^ ,
(14)
t ф
Учитывая закон теплообмена
St =
qo
»|дТ
У=0
pCpU (Ts- To) pCpU (Ts- To I
PCpU
T - To V Ts - To У
y=o
определим производную температурного пограничного слоя на стенке:
дУ
Т - To
V Ts - To У
д Г и
У=0 дУ IU » У=0
(10)
д ( У і7
дУ V О У ly=0 7-О1/7 - У6/7 |У=0'
Рассмотрим двухслойную модель пограничного слоя, с учетом (1) и (2) в ламинарном подслое имеем
u (U* }2 У U '
vU
тогда
д Г и
дУ V U
y=o
д
дУ
(/ *\2 і (U I У
vU
У=0
(11)
vU |y=° дУ
T - To V TS - To У
= const.
Следует отметить, что для практического применения закона теплообмена (14) необходимо определить значение ал . Воспользуемся выражением (3):
p
(U * I2
и выражением для закона трения на плоской пластине при турбулентном обтекании для степенного профиля распределения скорости:
pU2
= 0,012 56
г uО** і0 25
v
V У
тогда выражение (3) примет вид
= -° = 0,012 56U: p
( ЛЯ** і 0,25
U
V
V У
(15)
Из условия подобия динамического и температурного пограничных слоев получаем значение производной
д
дУ
Л
V Ts To У
д ( и
дУ V U
И2
vU
окончательно имеем
дУ
т - To
V TS- To У
0,012 56-U 0,
0,75с**0,25 v , Stv
(16)
Примем во внимание выражение для динамической скорости (7) и толщины потери энергии температурного пограничного слоя (9), окончательно запишем
Видно, что выражение (12), полученное из условия смыкания ламинарного подслоя с турбулентным профилем при у = 8л, и (16), полученное из закона трения степенного турбулентного профиля, имеют сходную форму и отличаются лишь численным коэффициентом. Потребовав равенства производных, определим коэффициент ал:
0,558 - U
дУ
т - To
V TS - Т0 У
2
an -V
vU-О*Ф,
1/4
0,558-U -
1/4
0,75
a„2 -v
(12)
0,012 56-U
,,0,75 5? **0,25 v , -°ф ,
U-S'
tф
a.2 -v
(17)
S
0
X
0
д
Лучшее совпадение с полученным значением а л = 12,559 дает значение, полученное в работе [3], расхождение составляет 10 %; при ал = 5, полученном в [4], расхождение в 3 раза.
При подстановке (14) в интегральное соотношение уравнения энергии, полученное в [5], учтем что, для вращательного движения (линия тока - кольцевая линия) удобнее выполнить запись в цилиндрических
^ дR
координатах: ф = а;; у = R ; Hф = R ;
ду dR
= 1;
3HV 5
Hy = 1; --------= 0; — = 0, тогда будем иметь
дф дф
д Js
sdR ‘ф + ~r ‘ф =
0,558 1
Ф0
(■ + s2)
(18)
п т 7-0,25 1,5.0,75 о**0,25 „г (Т Т V
PCpU аЛ v , 8^ pCp (Ts- T0)
Выражение для определения коэффициента теплоотдачи имеет вид
а = р • Cp • U • St.
(19)
После преобразования (18) с учетом (19) можно записать выражения для определения коэффициента локальной теплоотдачи для вращательных течений газового потока с различным распределением скоростей. Для вращательного течения по закону свободного вихря, с распределением скорости UR = C = const имеем:
2,18 -p-Cp-U
4
Pr5
2 J s
A:
Чал • Re®y
(20)
Для вращательного течения, индуцированного вращением диска (закон твердого тела), с распределе-
U
нием скорости R = ю= const
а =
0,528-p-Cp-U
4
Pr5
2Js
- Re„
(21)
В (23) и (24) J - относительная характерная толщина, полученная в работе [6], е - тангенс угла скоса донной линии тока.
Таким образом, получены зависимости для оценки локального коэффициента теплоотдачи вращательных течений газового потока при различном распределении скоростей; выражения позволяют учитывать локальную температуру движущегося потока, распределение температуры в теплопроводящей стенке, элементарное количество тепла, переданного через единицу поверхности; вывод выражений осуществлялся с учетом подобия температурного и динамического пограничных слоев, что характерно для частного случая Рг = 1.
Библиографические ссылки
1. Кейс В. М. Конвективный тепло- и массообмен : пер. с англ. М. : Энергия, 1972. 448 с.
2. Войткунский Я. И., Фаддеев Ю. И., Федяевский К. К. Гидромеханика : учебник. 2-е изд., перераб. и доп. Л. : Судостроение, 1982. 456 с.
3. Емцев Б. Т. Техническая гидромеханика : учебник для вузов. 2-е изд. перераб. и доп. М. : Машиностроение, 1987. 440 с.
4. Романенко П. Н. Тепломассообмен и трение при градиентном течении жидкостей. М. : Энергия, 1971.
5. Интегральное соотношение уравнения энергии температурного пространственного пограничного слоя / А. А. Зуев [и др.] // Вестник Рыбинской гос. авиац. технол. акад. имени П. А. Соловьева. 2010. № 2 (17). С. 37-42.
6. Относительные характерные толщины динамического пограничного слоя при различных законах распределения скорости / А. А. Кишкин [и др.] // Вестник СибГАУ. 2009. № 3 (24). С. 129-133.
A. A. Zuev, A. A. Kishkin, M. I. Tolstopyatov, D. A. Zhuykov
ROTATIONAL FLOWS HEAT TRANSFER IN TURBOMACHINES BASED ON TWO-LAYER MODEL OF TURBULENT BOUNDARY LAYER
Individual case of flow of gas flow induced by rotation of disk is considered in the article. Integral relations of energy equation of dimensional boundary layer for flow with rotary power law velocity distribution is presented. Dependencies for estimation of local heat transfer coefficients of rotational flows with different velocity distribution are obtained.
Keywords: convective heat transfer, turbulent flow, energy equation, heat transfer law.
© Зуев А. А., Кишкин А. А., Толстопятов М. И., Жуйков Д. А., 2012
