книги из ГПНТБ / Козырев, А. П. Теория тепловых и гидродинамических процессов в атомных энергетических установках учеб. пособие
.pdfЯа входном участке трубы образуется пристенная ооласть толщиной сгг , в которой сосредоточено все изме
нение температуры от величины Ьст до |
Эта об |
ласть называется тепловым пограничным слоем. Толщина
Рис. 6 .4 . Развитие профиля температуры на термичес ком начальном участке
слоя |
8 Т по мере продвижения |
потока вдоль канала рас |
|
тет |
и на расстоянии ■£ |
от |
входа происходит смыка |
ние тепловых слоев. Начиная с этого момента, в теплооб мене участвует вся жидкость, и происходит изменение температуры потока на оси канала.
Входной участок длиной ^ГС называется участком
тепловой стабилизации или начальным термическим участ ком. Стабилизация профиля температуры и скорости опре
деляется числом Рг = ■§- , поскольку кинематическая
вязкость 9 представляет собой коэффициент диффузии импульса, или скорости, а температуропроводность о - коэффициент диффузии тепла, или температуры. Коэффици ент диффузии определяется как скорость диффузии какойлибо субстанции в среде при градиенте потенциала, рав ном единице. При числе Рг = I тепло и импульс диффунди руют в жидкости с одинаковой скоростью, и толщина Рт определяется толщиной гидродинамического пограничного слоя <?Г . Если Р2 <с I , то молекулярная вязкость
меньше молекулярной теплопроводности, тепловые возмуще ния более интенсивно распространяются в ядро потока,
профиль |
температуры развивается быстрее профиля скорос |
||||
ти, при |
этом |
(?r |
, а ртс <вгс- При pz ^ |
I |
со |
ответственно |
с?г <с^. и |
-втс^ в г.с * Роль числа Рг |
в |
теп |
лообмене чрезвычайно велика. На рис. 6.5 показан диапа зон изменения чисел Рг для различных веществ.
Коэффициент теплоотдачи связан с градиентом темпе ратур уравнением теплообмена
210
|
|
Л t \dZ /г =г |
» |
(6 .9 ) |
|
|
'О |
|
|
где Л t = tcr- t |
- |
температурный |
напор; |
|
Z0 |
- |
радиус трубы. |
|
|
Ю'2 Ю ч Ю ° |
Ю* |
iO2 103 |
|
метлллы |
л&за/ |
|
Масло |
|
|
|
|
Рис. 6 .5 . Спектр чисел Прандтля |
для |
различных жид- |
|
|
костей |
|
|
Наибольший температурный градиент имеет место на |
входе, |
|||
где |
температура по всему сечению постоянна, а |
на |
стен |
|
ке |
скачок температуры равен tCT - |
. Поэтому |
на участ |
ке термической стабилизации локальный коэффициент тепло
отдачи |
о(_х уменмается |
и при термически стабилизиро |
||
ванном течении становится |
постоянной |
величиной.(рис .6.6 ). |
||
Средне интегральная по длине канала величина |
ос будет |
|||
больше_стабилизкрованного значения, |
и поэтому при рас |
|||
чете <=* |
необходимо вводить поправку |
на |
длину тру- |
Рис.6 .6 .Изменение коэффици- |
Рис. 6 .7 . Полностью раз |
|||||
вита теплоотдачи на вход- |
витый |
профиль темпера- |
||||
ном участке |
канала |
туры в |
трубе |
|||
Величина |
€ TZ |
зависит от |
числа |
/?е |
, |
коэффициента |
теплопроводности Л |
, входных условий, |
наличия гид |
||||
родинамической |
стабилизации |
и многих |
других |
211
факторов. При ламинарном течении жидкости |
гтс^ т |
ы , |
|
а при турбулентном |
~(10-г15)Ы и протяженность |
||
входного участка заметно меньше, чем при ламинарном |
|
||
режиме. |
|
|
|
При термически |
стабилизированном течении |
в трубе |
в |
отличие от полностью развитого неизменяемого по длине профиля скорости в сечении канала профиль температур изменяется по длине. Однако для некоторых способов обогрева, например при q = const или tCT = const , в трубе устанавливается профиль безразмерной темпера туры, который не изменяется по длине (рис. 6 .7 ). Обобщенный безразмерный профиль температуры выражается как
|
J _ ^СТ_ |
t |
|
- f |
|
|
|
|
|
^ст |
Ь |
где tж |
- средняя по |
энтальпии температура (темпе |
ратура смешения). Условие инвариантности профиля тем-
пературы |
Т |
по длине |
записывается в виде |
||||
|
|
д Т |
_ |
|
1 |
|
|
|
|
|
'-i |
II |
(б.Ю) |
||
|
|
дх |
|
|
|||
|
|
-дх( ( t - t |
|
||||
при |
этом |
при |
t = |
> |
T = i. |
|
в каналах на |
В |
зависимости от величины числа Re |
||||||
блюдается |
ламинарный, |
переходный или турбулентный ре |
жим течения. Факт перехода ламинарного режима в тур
булентный в гладких трубах при |
2000 был уста |
|||||
новлен Рейнольдсом |
в |
1883 г . |
Поздней экспериментально |
|||
было установлено, |
что |
Re кр |
зависит |
от |
степени |
воз |
мущения потока на входе. Б частности, |
в |
отдельных |
опы |
|||
тах путем тщательного |
устранения внешних |
возмущений |
потока удалось затянуть переход к турбулентному режи
му до Re^n^lO*. При |
Re > RekD в технических тру |
хла |
<2 |
212
бах и каналах устанавливается развитое турбулентное те чение. Ему соответствует турбулентный режим теплоот дачи. Режим течения при Re = 2 • 10s -f 10^ называет ся переходным. Закономерности теплоотдачи для всех трех режимов различны.
Универсальное уравнение теплоотдачи при
стабилизированном течении
Изучение теплоотдачи в настоящее время идет как по пути теоретических исследований, так и посредством экспериментирования. Все большее развитие получают аналитические методы исследования, позволяющие оценить влияние основных факторов на теплоотдачу, а в ряде за дач получить точные расчетные зависимости (см ., напри мер, [ 38, 42, 59 ] ) i Инженерный анализ многих при кладных задач конвективного теплообмена базируется на интегрировании уравнения энергии с привлечением гид родинамической аналогии Рейнольдса. Сущность такого анализа рассмотрим на примере вывода известного уни версального интегрального соотношения Лайона, позво ляющего получить целый ряд важных результатов в част ных случаях.
Пусть стационарный осесимметричный поток гидродина мически и термически стабилизирован. Если пренебречь зависимостью физических свойств жидкости от темпера туры и теплопроводностью в осевом направлении
( ^ / 'дх% = |
0) |
, то осредненное уравнение энергии |
в цилиндрических координатах принимает вид |
||
а |
д Ч |
i dt |
дг г г д?
213
Вводя коэффициент |
турбулентной |
теплопроводности д |
, |
||||||
уравнение энергии |
( б .II) можно |
записать |
в виде |
|
г |
||||
_д |
г(Л + Л т) |
dt |
|
_ |
dt |
(6. 12) |
|||
дг |
дг |
|
|
|
|||||
L |
|
|
|
|
|
|
|||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
(6 .IS ) |
|||
д г W |
А ) д г I |
а |
д х |
||||||
|
|
||||||||
В уравнениях |
(6 .12), (б .13) и далее |
знак |
осреднения |
во |
времени актуальных значений температуры и скорости для простоты опущен.
да |
получим |
|
( б .13) |
от |
0 до |
г |
. Тог- |
||
|
г |
|
|
|
|
||||
|
|
|
|
|
|
|
|||
|
|
|
Г |
и>г |
dt . |
|
(6 .14) |
||
|
А |
! дг |
а |
дх |
dz . |
|
|||
|
|
х |
|
|
|
||||
Поскольку |
<°dt |
, |
|
, |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
I |
1 7 * - Ь ' * ' |
|
|
|
|
|||
то |
после интегрирования |
(6.1*0 |
от |
г |
до |
г |
имеем |
||
|
t.ст- ‘ - у |
(1 + * т/тдг. dz. |
|
(6.15) |
|||||
Определим производную |
^ / дх |
|
• |
|
|
|
|
Если плотность теплового потока вдоль трубы постоянна
( |
~ const) , |
тогда |
из уравнения |
баланса тепла мож |
но |
написать |
|
|
|
|
dtж |
— |
= const . |
(6.16) |
|
d F |
Gc’.Р |
|
|
где t
р- поверхность нагрева трубы;
Ы Р = £fczodx. Тогда (б.1б) перепишется в виде
|
|
|
_ |
£ Чет |
2ct(tCT^ |
tm ) |
|
|
|
||||
ы х |
|
~ f>Cpmx %0 |
|
|
|
= const. (6.17) |
|||||||
|
J*cP |
^ |
Zo |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
||||
Если |
U |
= |
const ( |
термически |
стабилизированное |
||||||||
течение), |
то |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
tcr |
~ |
= |
Con&t |
|
|
|
|
|
|
||
и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dter |
dt«, |
|
|
|
|
|
|
|
|
||
|
с/х |
|
— — = const |
|
|
|
|
|
(6.18) |
||||
|
|
d x |
|
|
|
|
|
|
|
||||
Дифференцируя |
условие (б.Ю ) |
и решая |
его |
относительно |
|||||||||
'х |
|
|
|
получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
&L |
= d tCTf t - tm |
|
dt* |
t —/ |
|
|
||||||
|
|
|
|
l'C T |
L |
(6 .19) |
|||||||
|
dx |
|
d x ( tcT- ^ |
|
d x |
^ст |
|
|
|||||
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Последнее выражение с учетом (6.18) дает |
|
|
|
||||||||||
|
dt |
|
dtCT |
dt* |
2d(tCT- t J |
|
= const.(6.20) |
||||||
|
dx |
|
dx |
d x |
|
^ сР Ч с го |
|
||||||
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|||
Таким образом, |
при заданном граничном условии вели |
||||||||||||
чины tej-, |
|
|
, |
t(%) |
меняются вдоль |
оси |
канала |
||||||
линейно (рис. |
|
6.8) |
|
|
|
R = у |
|
W |
- |
||||
Вводя |
безразмерные величины |
, |
|||||||||||
и учитывая |
(6 .20), |
преобразуем |
зависимость |
(6 .15) к* |
|||||||||
виду |
|
|
|
|
|
|
|
|
|
|
|
|
|
215
t - t = |
Щ |
|
f t / w * |
U |
1в.я , |
CT |
J*c, |
R |
( ' * |
|
|
|
P |
|
|
||
ИЛИ |
!г J W R d R |
|
|
||
|
|
|
|||
7 = Nu f f f - A i dR |
|
( 6. 22) |
По определению средняя безразмерная температура равна единице и в то же время определяется интегралом
ZI
Т= — - Г и>гТс(г = 2[WRTdR = i . (б .23)
to г? J |
J |
Рис. 6 .8 . Распределение |
температур по длине трубы |
|||
|
|
при |
= Const |
|
Подставим (6.22) |
в (6.23)* е |
|
||
|
|
'гJWRdR |
|
|
- L . t f w x f f |
dRldR . |
(6.24) |
||
« и |
I |
( { |
( *> v ) * |
|
Преобразуя |
(6 .24), получим известный интеграл |
Лайона |
||
j |
d (JWRdR)2 |
(6.25) |
||
_ L |
= А |
11---------------d R . |
||
Nu |
} |
0 fKr/*,)R |
|
2 1 6
Уравнение (6.25) является универсальным, поскольку может быть использовано для анализа теплообмена в ламинарных и турбулентных потоках теплоносителей с различ
ными числами |
P z |
|
. Для расчета теплоотдачи по форму |
||
ле (б .25) |
необходимо знать закон распределения скоростей |
||||
и отношение |
; |
^ т / , |
, |
. |
|
_ |
|
/ л |
|
||
Выражение |
|
/ |
+ |
Ат/^ » входящее в интеграл, можно |
|
преобразовать |
следующим образом: |
Лт |
, |
Л |
9со/э |
Яг |
А |
/ -ь — |
= { + |
|
л |
|
|
Л |
|
|
|
|
|
|
|
|
|
|
|
где Pz |
= |
- |
турбулентное |
число |
Прандтля* по |
физическому смыслу определяет неподобие рассеивания теплосодержания и количества движения при турбулентных
пульсациях. Коэффициенты а |
и ^)г |
учитывают тур |
|
булентный перенос. Они аналогичны коэффициентам тем |
|||
пературопроводности и кинематической вязкости, имеют |
|||
ту же размерность [м ^/сек]и |
соответственно |
называются |
|
коэффициентами турбулентного |
переноса |
тепла |
и количест |
ва движения . Однако в отличие от « и ^ эти коэффициенты не являются физическими параметрами,
а |
зависят от числа |
Re |
, |
рг |
и расстояния от стен |
|
ки |
ч + . Число |
Ргт |
находится |
в пределах от 0,5 до |
||
2 и в расчетах, |
исходя |
из |
гидродинамической аналогии |
|||
процесса переноса тепла и количества движения, часто |
||||||
принимают Рг |
I . Коэффициент |
связан с полем |
скоростей и может быть найден дифференцированием эпю
ры скоростей, полученной экспериментально |
или задан |
|
ной определенным законом. |
Н и |
|
Согласно уравнению (6.25) число |
определяет |
|
ся не только гидродинамикой потока, |
но и числом Рг . |
|
Методы аналитического расчета теплоотдачи, |
основанные |
|
на использовании интеграла Лайона, |
интенсивно разви |
217
ваются и подробно рассмотрены в монографиях по теплооб мену [ 38, 42 ] . Частные случаи применения зави симости будут рассмотрены ниже.
Теплоотдача._при ламинарном режиме
При ламинарном режиме течения ('/?е$2000) отсутству ют турбулентные пульсации. Тогда при Л г = 0 интеграл Лайона упрощается и принимает вид
-i |
/ г* |
\2 Ы R |
|
|
- г ! [ vtRdR)T |
(6.27) |
|
|
|
|
Если течение гидродинамически и термически стабили зировано и неизотермичность потока слабо влияет на фи зические свойства жидкости, то параболический закон распределения скоростей в ламинарном потоке можно за писать в относительных величинах:
|
|
W |
= 2 (i - R 2). |
|
|
(6.28) |
|||
С учетом (6.28) можно рассчитать значение интеграла |
|||||||||
(6 .27): |
|
|
|
i |
я |
|
|
|
|
|
|
/ |
|
|
н_ |
|
|
||
|
|
ej[2j(l-R2) M R ] |
|
|
|||||
|
|
Nu |
|
|
|
||||
|
|
о |
О |
|
т |
да |
|
|
|
или |
|
|
|
|
|
|
|
|
|
|
|
|
Nu = |
-- |
= R,36 |
|
|
(6.29) |
|
|
|
|
|
|
|
|
|
|
|
Таким образом, |
коэффициент |
теплоотдачи |
оС |
при ламинар |
|||||
ном режиме |
для |
любых теплоносителей зависит только от |
|||||||
Д |
и |
d 0 |
и не |
зависит |
от и?, у 3 , |
ср |
, что |
объясняется |
чисто молекулярным процессом переноса теп |
|
ла. Решение |
(6.29) получено при условии |
const . |
218
При граничном условии |
± = const было получено, |
что |
|
И.и - |
3,66 |
|
(б.ЗО) |
Полученные теоретические решения (б .29) и (б .30) не учитывают влияния поля температур на физические свой ства потока и поле скоростей и поэтому совпадают с опытными данными только прималых температурных напо рах. При переходном и турбулентном режиме поле темпе ратур (неизотермичность потока) несущественно влияет на поле скоростей. Однако при больших градиентах тем ператур по сечению в ламинарном режиме теряется его основное свойство - параболический закон распределе ния скоростей. Кроме того, наличие неодинаковой тем пературы по сечению является причиной возникновения подъемных сил и свободной конвекции, которая наклады вается на вынужденное движение. Таким образом, в ус ловиях подвода или отвода тепла при любой ориентации трубы в поле тяжести всегда возникают вторичные тече ния, обусловленные разностью плотностей жидкости в по токе. О влиянии этих вторичных течений на гидродина мику ламинарного потока можно судить по соотношению сил вязкости и подъемных сил, т .е . по критерию G % или Ra . Существуют предельные числа Рэлея, ниже которых свободной конвекцией можно пренебречь и выше которых ее следует учитывать.
В зависимости от влияния свободной конвекции на теплоотдачу различают вязкостный и вязкостно-гравита ционный режим неизотермического движения. При вязкост ном режиме силы вязкости преобладают над подъемными силами. При этом режиме температурный напор в потоке сказывается только на деформации параболического про филя скорости из-за влияния температуры на вязкость.
219