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

tot_book

.pdf
Скачиваний:
40
Добавлен:
12.03.2015
Размер:
4.69 Mб
Скачать

241

Все сказанное справедливо для стабилизированного движения, которое наблюдается лишь на некотором удалении от входа жидкости в трубу. Процесс стабилизации сводится к следующему.

При входе в трубу вблизи стенок появляется динамический пограничный слой. Если течение ламинарное (рис. 118, а), то толщина этого слоя δl плавно увеличивается; при δl = r0 слой заполняет все сечение канала, после чего течение стабилизируется. На начальном участке l0 течение носит ламинарный характер, но эпюра скоростей здесь еще не имеет традиционного вида (2.165). Если же течение турбулентное (рис. 1 18,6), то на участке l0 происходит смыкание турбулентных пограничных слоев толщиной δт, а затем развитая турбулентность наблюдается почти во всем сечении трубы; лишь у самых стенок остается тонкий вязкий подслой, в котором скорость падает до нуля.

Длина начального участка при ламинарном режиме l0 = 0,05dRedf (т. е. тем больше, чем больше Redf), при турбулентном

l0 =15d (не зависит от числа Redf ).

Описанная картина течения определяет особенности конвективного теплообмена в трубах и каналах.

Рис. 118.

242

При ламинарном течении теплота передается в радиальном направлении только теплопроводностью. Поскольку слои жидкости движутся с различными скоростями w(у) (см. рис. 118,а), вдоль потока теплота передается конвективным путем.

Если на входе в трубу (при х = 0) температура жидкости Тf, постоянна по сечению и в общем случае не равна температуре стенок Tw, то по мере движения вдоль трубы значение Тf меняется. Вначале теплообмен происходит только в тонком пристенном слое (А), затем в передачу теплоты вовлекается все большая часть потока (В, С, D, Е) (рис. 119). Картина сходна с формированием эпюры скоростей (рис. 118,а). Толщина теплового пограничного слоя l увеличивается, на удалении от входа l0h (англ, heat — тепловой) l = d/2, слои смыкаются, зависимость Tf(d) становится параболической (эпюра D). Обычно считают, что l0h = 0,05dRedfPrf; на практике ламинарный режим характерен для движения вязких жидкостей, таких как масло, дизельное топливо и т. д.

Рис. 119.

243

Вспомним, что коэффициент теплоотдачи определяется в общем случае по формуле (2.118). На участке х l0h градиент

температуры Tn n=0 убывает быстрее, чем разность Tf Tw, поэтому

местный коэффициент теплоотдачи снижается. При х > l0h обе величины меняются пропорционально, значение α сохраняется почти неизменным. Поскольку речь идет о местном, зависящем от координаты х, значении α, температуру Tf рассматривают как среднюю в соответствующем сечении F:

 

 

 

 

1

Tf (x,r)dF.

Tf =Tf

=

F

 

 

 

 

F

 

 

 

 

 

Средний коэффициент теплоотдачи на участке длиной х

α = Tf qTw ,

где q — среднее значение q(x) на длине участка; Tw — средняя температура стенок на той же длине.

В случае обогрева стенок паром, когда изменение Tf (x) невелико, либо принимают

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Tf 1

+Tf 2

 

 

 

 

 

Tf =

,

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в начале и в конце рассматриваемого

где Tf 1,Tf 2 — значения Tf

 

участка трубы, либо рассматривают среднелогарифмическую температуру:

Tf =Tw ± T1 −∆T2 , ln T1

T2

где T1 = Tf 1 Tw1,T2 = Tf 2 Tw2 — температурные напоры в

начале и конце участка; знак + выбирают, исходя из направления теплового потока.

При электрообогреве стенок (q(x) = const) всегда пользуются первой из приведенных выше формул.

244

На рис. 119 видно, что участок, на котором стабилизируется α, всегда больше, чем участок стабилизации α:

l0h > l0h.

Для прямых труб постоянного сечения обычно определяют именно среднее значение α.

2.8.6. Теплообмен на стабилизированном участке течения. Интеграл Лайона

Рассчитаем коэффициент теплоотдачи на участке стабилизированного теплообмена (α = const) при граничных условиях II рода qw = const. При массовом расходе Gm на элементе длины dx выполняется первое начало термодинамики:

 

 

 

 

 

 

 

 

= πdα(Tw Tf )dx = qwπddx,

(2.167)

 

 

 

 

cpf GмdTf

где d = 2r0 — диаметр трубы, поэтому

 

 

 

 

 

 

 

 

q

 

πd

 

 

 

 

4q

 

πd

 

2q

 

 

 

 

dTf

=

w

=

 

 

w

=

w

= const.

(2.168)

 

dx

cpf Gм

 

c

G πd 2w

cpf ρf r0w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pf м

 

 

 

 

 

 

Таким образом, при qw = const средняя температура Tf увеличивается линейно. Определим среднюю энтальпию h и среднюю температуру Tf в любом сечении трубы F0. Поскольку и скорость, и энтальпия распределены по сечению F0 неравномерно, средняя энтальпия

 

 

 

ρf whdF

 

 

 

 

=

F0

.

(2.169)

h

 

ρf wdF

 

 

 

 

 

 

 

 

F0

 

 

Равенство (2.169) для несжимаемой жидкости при h = cpfTf и cpf = const переходит в выражение для средней температуры

 

Tf wdF Tf wdF

 

 

Tf =

F0

=

F0

.

(2.170)

wdF

wF0

 

 

 

 

F0

245

Для круглой трубы dF = 2πrdr (здесь r — текущий радиус) и равенство (2.170) принимает вид

 

 

 

 

 

 

r0

w r dr

1

(2.171)

 

 

 

Tf = T

= 2TWRdR,

 

 

 

w

 

r

 

r

 

 

 

0

0

0

0

 

где W =

w

; R =

r

.

 

 

 

 

 

 

 

w

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

Интеграл (2.171) возьмем по частям, полагая Т = u, WRdR = dv:

 

 

 

 

1

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Tf

 

 

 

 

 

 

 

(2.172)

= 2 Tw WRdR

 

WRdR dT .

 

 

 

 

0

 

0

0

 

 

 

 

 

 

 

 

 

Рассмотрим первый интеграл правой части уравнения (2.172):

 

r0

 

r0

 

 

 

1

wrdr

 

2πwrdr

 

1

 

WRdR =

0

=

0

=

.

w0r02

2πw0r02

2

0

 

 

 

Отсюда следует, что

 

 

 

 

 

 

1

 

1

 

(2.173)

T

f

w

2

 

 

 

 

=T

 

 

 

WRdR dT.

 

 

 

 

 

 

0

0

 

 

Для того чтобы определить член dT, обратимся к уравнению энергии (2.123), которое в цилиндрических координатах и при осевой симметрии потока примет вид

c

 

 

ρ

 

T

+ w

T

 

+ w

T

=

 

 

 

 

 

pf

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

τ

 

 

 

 

r

 

 

 

 

 

 

=

1

 

(λ f T )r

 

T

 

+

(λ f

T )

T

,

 

 

 

 

 

 

(2.174)

 

 

r r

 

 

 

 

r

 

x

 

 

x

где wr — радиальная составляющая скорости.

 

 

В

стационарном

 

режиме

T

= 0,

а

на участке

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

τ

 

 

стабилизированного течения wx = 0, и, как следует из уравнения неразрывности, wxr = 0, т. е. wr = const. Однако на стенках (при r =

246

r0) w = 0, значит, в любой точке сечения wr = 0. Особо следует сказать о λT турбулентной теплопроводности. Рассматривая теплообмен в турбулентном пограничном слое, мы установили, что турбулентные пульсации переносят (дополнительно к теплопроводности) значительную теплоту поперек основного течения. Влияние пульсационных процессов учитывает закон, аналогичный закону Фурье, поэтому совместный перенос теплоты молекулярной и турбулентной теплопроводностью можно оценить по соотношению

q = −(λ f T )gradT.

 

(2.175)

Отметим, что слагаемые в правой части уравнения (2.174)

неравноценны. Как показал еще Л. Прандтль,

2T

= 0, поэтому

 

x2

 

последним слагаемым можно пренебречь, после чего уравнение примет вид

c

pf

ρ

f

w

T

=

1

(λ

f

T

)r

T

.

(2.176)

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r r

 

 

 

 

r

 

На стабилизированном участке течения (и теплообмена) не только средняя, но и местная температура должна изменяться линейно, поэтому в уравнении (2.176) производную в левой части можно заменить на константу из уравнения (2.168), при этом

cpf ρ f w

 

 

 

 

2qw

 

=

1

 

(λ f

T )r

T

.

(2.177)

 

cpf ρ f r0w

 

 

 

 

 

 

 

 

 

 

 

r r

 

 

 

 

 

r

 

Вновь используем безразмерные скорость W и радиус R:

 

2q

 

 

 

r

 

 

 

 

 

 

 

 

λ

 

 

 

T

 

 

 

 

 

 

w 0

 

 

 

 

 

+

 

T

 

R

 

;

(2.178)

 

 

 

 

 

 

 

 

 

 

 

 

λ

 

 

 

 

WRdR = d 1

λ

 

 

R

 

 

f

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

проинтегрировав, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2q

w

r

 

R

 

 

 

 

λ

 

 

 

T

 

 

 

 

 

 

 

 

 

0

WRdR = 1

+

λ

T

 

R

R

.

 

 

(2.179)

 

 

 

 

 

 

 

 

 

 

 

λ

f

 

0

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение (2.179) разрешим относительно T , после этого уравнение (2.173) примет форму

 

 

 

 

 

 

 

 

247

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dR

 

 

 

 

 

4q

r

1

 

WRdR

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

Tf

= Tw

 

w 0

 

 

 

 

.

(2.180)

λ f

 

 

 

λ

 

 

 

 

 

 

 

0

 

+

 

R

 

 

 

 

 

 

 

 

1

λ

T

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Это равенство удобней представить в безразмерной форме:

 

 

 

 

 

 

 

 

R

 

 

2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dR

 

 

 

 

 

 

 

 

 

 

 

WRdR

 

 

 

 

 

 

 

 

 

0

 

 

 

 

(2.181)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Nudf

= 2

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

1+ λT

 

 

 

 

 

 

 

 

 

 

 

λ

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение (2.181) (интеграл Лайона) позволяет вычислить

 

df

=

αd , а затем

и

α

 

при

 

условии, что известен

профиль

Nu

 

 

 

 

 

λ f

 

 

 

 

 

 

 

 

 

 

скорости на участке стабилизированного течения. Формула справедлива как для ламинарного, так и для турбулентного течения. Рассмотрим оба варианта подробнее.

2.8.7. Теплообмен при ламинарном течении в трубах

При ламинарном режиме Redf < 2300, λT = 0, а профиль скорости на стабилизированном участке имеет вид (2.166) или, поскольку w0 = 2w ,

 

 

 

 

 

 

2

 

 

2(1

 

 

 

).

 

 

 

 

 

 

 

 

 

 

 

 

w

=

2w

 

r

 

 

 

=

R

2

(2.182)

 

 

 

1

r

 

 

, W

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

Подставим значение W в интеграл Лайона (2.181) при λT = 0, найдем

 

 

 

 

 

 

 

 

 

 

 

 

 

248

 

 

 

 

 

 

 

 

 

 

 

 

 

R

2(1R

2

 

 

2

 

1

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)RdR

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Nudf

=

2

 

 

 

 

 

 

 

dR

=

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

R

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

48

 

 

 

λ

f

 

 

 

= 8

2(R R3 )dR dR

 

=

4,36,

α = 4,36

 

 

.

 

 

11

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

d

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.183)

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким

 

 

 

образом,

 

коэффициент

теплоотдачи

на

стабилизированном участке течения и теплообмена автомоделен (т. е. независим) относительно числа Рейнольдса Redf. Однако формула (2.183) имеет ограниченное применение. Дело в том, что длина участка тепловой стабилизации при Redf < 2300 оказывается очень большой. Действительно, при qw = const l0h/d = 0,055 Redf Prf. Для высоковязких жидкостей, например для органических теплоносителей, у которых Рrf = 10…100, при Redf = 103 получаем l0h/d 5,5·(102…103), и при d = 10 мм = 10–2 м длина участка стабилизации составит l0h = 5,5·(1…10) м. Как правило, реальная длина теплообменных участков оказывается значительно меньшей, а это значит, что пограничные слои не успевают сомкнуться и число Нуссельта можно рассчитывать по формуле (2.148). Необходимо только учесть, что жидкость в ядре потока на входном участке трубы ускоряется, чего нет при обтекании плоской поверхности, поэтому формула для расчета местных коэффициентов теплоотдачи (при l << l0h) имеет вид

Nu xf

где Nu xf = αx

λ f

= 0,335Re0,5

Pr0,33

Prf

0,25

x 0,1

(2.184)

 

 

 

 

 

.

Pr

 

xf

f

 

 

d

 

 

 

 

w

 

 

 

 

; Rexf = wνx ; x — координата, отсчитанная от

входного сечения трубы.

Если l и l0h — величины одного порядка, то формула (2.184) оказывается непригодной, и необходимо использовать эмпирическое уравнение подобия

249

 

 

 

 

d 0,4

Pr0,33

Prf

0,25

 

 

 

(2.185)

Nu

df

=1,4 Re

df

 

 

 

 

ε

l

,

Pr

 

 

l

f

 

 

 

 

 

 

 

 

 

 

 

 

w

 

 

 

 

которое справедливо при Redf > 10 и l/d > 10 (здесь εl — коэффициент, учитывающий зависимость α от l/d и Redf),. Формулы (2.183)–(2.185) не учитывают, что плотность жидкости по сечению и длине трубы при теплообмене не остается постоянной. Это приводит к появлению подъемной силы, а значит, в число определяющих критериев должно войти число Рэлея Radf = Grdf·Ргf. Однако эксперименты показали, что подъемная сила заметно влияет на теплообмен лишь при Radf > Racr = 8·105. Режим теплообмена при Radf < Racr называют вязкостным, и формулы (2.183)-(2.185) относятся именно к этому режиму. При Radf > 8·105 теплообмен называют вязкостно-гравитационным, и здесь Nudf = f(Redf, Radf). Теплообмен при вязкостно-гравитационном режиме исследован недостаточно, для расчета коэффициентов теплоотдачи следует обращаться к специальной литературе.

2.8.8. Теплообмен при турбулентном течении в трубах

Профиль скорости в трубе на стабилизированном участке при Redf > Recr обычно описывают набором из трех уравнений, пригодных соответственно для вязкого подслоя, пристенного турбулентного слоя и турбулентного течения в ядре. Однако с достаточной точностью можно представить профиль скорости степенной зависимостью (закон одной седьмой):

w

 

(r

r)v

 

17

 

(2.186)

 

=8,74

0

 

 

,

v

 

 

 

ν

 

 

 

 

в которой в качестве масштаба скорости использована так называемая «динамическая» скорость (скорость трения)

v = Sw .

 

ρf

 

250

Формула (2.186) удовлетворительно согласуется с известным соотношением для коэффициента трения при турбулентном течении в трубах

ξ =

0,31

.

(2.187)

1

 

Redf4

 

Кроме профиля скорости при турбулентном течении в интеграл Лайона (2.181) необходимо подставить значение турбулентной теплопроводности λT. Обычно эффективную теплопроводность λT+λf представляют в виде λf (1+λT/λf), для чего необходимо определить отношение

 

λ

 

 

λT cpf

ρ f

 

a

 

aT νT ν f ρ f

 

 

T

=

 

 

 

 

=

T

=

 

 

 

cpf ρ f λ f

a f

a f νT ν f ρ f

 

λ f

 

 

 

где

a

=

 

λT

 

— турбулентная

 

 

 

 

 

T

 

 

cpf ρf

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

Prf

µ

T ,

(2.188)

Pr

µ

 

f

 

 

T

 

 

температуропроводность;

Pr = νT — турбулентный

аналог

числа

Прандтля; µT

T

aT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

турбулентная динамическая вязкость.

 

 

 

Обычно PrT

 

 

считают

постоянной величиной ( PrT =0,8).

Величину µT определяют, исходя из модели Л. Прандтля,

следующим образом:

 

 

 

 

 

 

 

 

 

 

µ

T

= ρ

f

χ2 (r

r)2

w

,

(2.189)

 

 

 

 

 

 

 

 

0

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

где χ = 0,4 — одна из универсальных констант турбулентности.

 

Подставляя значение w из равенства (2.186) в формулу (2.189)

и учитывая, что v

= w χ

8

, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

µT

= 0,25Redf (1R)87 Prf .

(2.190)

 

µ f

 

 

 

 

 

 

 

 

 

 

Интеграл Лайона (2.181) с учетом равенств (2.186) и (2.190) позволяет получить формулу

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]