Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Козырев, А. П. Теория тепловых и гидродинамических процессов в атомных энергетических установках учеб. пособие

.pdf
Скачиваний:
13
Добавлен:
22.10.2023
Размер:
15.28 Mб
Скачать

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

нение температуры от величины Ьст до

Эта об­

ласть называется тепловым пограничным слоем. Толщина

Рис. 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

 

 

 

^ст

Ь

где

- средняя по

энтальпии температура (темпе­

ратура смешения). Условие инвариантности профиля тем-

пературы

Т

по длине

записывается в виде

 

 

д Т

_

 

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