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

tot_book

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

231

Распределения скорости и избыточной температуры в ламинарном пограничном слое хорошо аппроксимируют полиномы третьей степени:

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

y

 

 

2

 

 

y

 

 

3

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

;

 

w

x

= w

 

0

+ a

 

 

 

+ a

 

 

 

 

 

+ a

 

 

 

 

(2.137)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

1

 

 

 

 

 

 

2

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δl

 

 

 

δl

 

 

 

δl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

y

 

2

 

 

y

 

3

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

,

 

ϑ = ϑ

0

+b

 

 

+b

 

 

 

 

 

 

+b

 

 

 

 

 

(2.138)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

l

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где a0,...,а3 и b0,...,b3 — постоянные коэффициенты. Коэффициенты полиномов (2.137) и (2.138) можно определить

из граничных условий для скорости и температуры:

wx

 

 

 

 

= 0; wx

 

= w0;

wx

 

 

= 0;

2wx

= 0;

(2.139)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y=0

 

 

 

 

 

y=δl

y

 

y=δl

 

y

2

y=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2ϑ

 

 

 

ϑ

 

 

 

= 0;

ϑ

 

 

= ϑ0;

∂ϑ

 

 

= 0;

 

 

= 0

(2.140)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y=0

 

 

 

 

y=∆l

y

 

y=∆l

 

y

2

 

y=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(последние условия в равенствах (2.139) и (2.140) соответствуют первым двум уравнениям системы (2.134), поскольку

wy y0 = wx yδl 0 ). Подставляя соотношения (2.139) и (2.140) в

полиномы (2.137) и (2.138) соответственно, получим

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

y

 

3

 

 

 

 

 

 

 

3

 

 

 

 

 

1

 

;

 

w

x

= w

 

 

 

 

 

 

 

 

(2.141)

 

 

 

 

 

 

 

 

 

 

0

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

δl

 

 

δl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

y

 

 

1

 

y

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

ϑ = ϑ

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

(2.138)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

l

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Расчет коэффициента теплоотдачи а на основе последнего уравнения системы (2.134) приводит теперь к выражению

α =

3

λ f

.

(2.143)

2

l

 

 

 

232

Таким образом, для окончательного расчета а необходимо найти l, например, из уравнения (2.136), но туда входит еще и δl, (косвенно, через wx), поэтому придется решить также уравнение импульсов (2.135). Подставляя wx из уравнения (2.141) в уравнение

(2.135), найдем

δl =

280 νx ,

 

(2.144)

 

 

13

w0

 

 

а подставляя wx, δl, ϑ в уравнение (2.136), получим

 

δl

 

=

1

;

 

(2.145)

l

 

3

Prf

 

 

l =

280 νx

1

 

(2.146)

 

 

13

w

3 Pr

f

 

 

 

 

0

 

 

(чтобы получить последние два результата, мы допустили, что l ≤ δl и что l и δl одинаковым образом зависят от продольной координаты х). Подставим теперь l в формулу (2.143); получим

α = 0,335λ f

νx 3 Prf ,

(2.147)

 

 

 

w0

 

откуда легко получить уравнение подобия

 

Nu xf =

αx

= 0,335Re0xf,5 Pr0f ,33 .

(2.148)

 

 

λ f

 

 

Этот результат позволяет рассчитать местные значения α, в то время как для практических расчетов часто достаточно знать средние по длине пластины значения коэффициента теплоотдачи

 

 

 

1

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α =

 

αdx при L

< xcr1. Интегрирование α по х

приводит к

L

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

уравнению подобия

 

 

 

 

 

 

 

 

 

 

 

 

 

0,5

0,38

Pr

f

0,25

(2.149)

 

 

 

 

 

NuLf

 

 

 

 

 

 

 

 

= 0,67 ReLf

Prf

Pr

 

.

 

 

 

 

 

 

 

 

 

 

 

w

 

233

В уравнение (2.149)

дополнительно введен поправочный

 

Pr

f

0,25

 

множитель М.А. Михеева

 

 

 

. Смысл его станет ясен, если

Pr

 

 

 

 

 

 

w

 

вспомнить, что для капельных жидкостей число Прандтля Рr снижается с увеличением Тf,. Отсюда следует, что при Tf > Tw

Prf

0,25

 

 

 

 

 

 

 

Prf

 

0,25

 

 

 

 

 

 

<1,

а

при

Tf <

Tw

 

 

 

 

>1,

т. е. множитель

Pr

 

Pr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

 

 

 

 

 

 

 

 

w

 

 

Pr

f

0,25

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

учитывает

 

связь

коэффициента

теплоотдачи с

Pr

 

 

 

 

 

 

 

w

 

 

 

 

 

 

 

 

 

 

 

 

 

направлением теплового потока.

 

 

 

 

 

 

 

 

 

 

 

Pr

f

0,25

 

 

 

 

 

 

 

 

 

 

Поправка

 

 

 

необходима

при

расчете коэффициента

 

 

 

Pr

 

 

 

 

 

 

 

 

 

 

 

 

 

w

 

 

 

 

 

 

 

 

теплоотдачи для капельно-жидких теплоносителей; для газов число Рr меняется слабо и при умеренных перепадах

Prf 0,25

температуры =1.

Prw

2.8.4. Теплообмен в турбулентном пограничном слое на плоской поверхности

Турбулентное течение, как отмечалось, характеризуется пульсациями скорости, температуры, давления и т. д. Если использовать достаточно точный измеритель скорости в пограничном слое, который сможет фиксировать мгновенные ее значения wx и wy, то при продольном обтекании плоской поверхности мы получили бы примерно следующие графики (рис. 114).

234

Рис. 114.

Скорость wx в любой момент времени можно представить в

виде

 

wx = wx + wx ,

(2.150)

где

wx — среднее значение скорости wx, a

wx — продольные

пульсации.

 

 

Поперечная составляющая скорости

 

 

wy = wy + wy ,

(2.151)

но

wy = 0, поэтому wy = wy , где wy — поперечная пульсация

скорости wy. При отсутствии значительных продольных градиентов температуры составляющая скорости wx не оказывает заметного влияния на теплообмен, а вот составляющая wy в пограничном

слое, где поперечные перепады температуры и скорости значительны, играет существенную роль. Если представить, что в турбулентном пограничном слое макроскопическая, но малая частица 1 переместилась в положение 2 (рис. 115), то условие неразрывности потребует, чтобы переход был скомпенсирован: частица 2 должна перейти в положение 1. Переходы 1–2 и 2–1 обеспечат сохранение сплошности потока, и уравнение неразрывности не нарушится. Однако частица 1 имела скорость wxl, а перешла в область частицы 2, где средняя продольная скорость равна wx2, и наоборот. Подобных обменов происходит много,

235

поэтому “сверху вниз” идет массовый поток ту, равный массовому потоку частиц “снизу вверх”. Кажется, что такой обмен должен приводить к выравниванию скорости, однако в действительности профиль средней скорости сохраняет устойчивость.

Рис. 115.

Это говорит о том, что выравниванию скоростей препятствует какая-то сила, а именно сила трения между слоями. Таким образом, на основе теоремы сохранения количества движения мы имеем право записать

ST

= my (wx1 wx2 ).

(2.152)

 

&

 

Здесь ST — турбулентное касательное напряжение

трения;

m&y — массовый поток, связанный с турбулентными пульсациями. Но этот же массовый поток m&y переносит энтальпию: частица 1

несет поток энтальпии cpf T1m&y , а частица 2 — cpf T2m&y . Разность

потоков энтальпии определяет плотность теплового потока в направлении, перпендикулярном поверхности:

qt = cpf m&y (T1

T2 ).

(2.153)

236

Исключая m&y из уравнений (2.152) и (2.153), получим

qt =

ST (T1

T2 )cpf

.

(2.154)

 

 

 

wx1 wx2

 

Уравнение (2.154) записано для двух произвольных точек 1 и 2. Если его применить для крайних точек турбулентного пограничного слоя (на верхней границе wx1 = w0 , T1 = Tf , , a на

нижней wx2 = wv — скорости на границе вязкого подслоя толщиной δv, T2 = Tv — температуре на границе вязкого подслоя), то получим

q

=

ST (Tf

Tv )cpf

.

(2.155)

 

 

t

 

w0

wv

 

 

 

 

В пределах вязкого подслоя δv, где скорость и температура изменяются линейно, плотность теплового потока и напряжение трения остаются постоянными, следовательно, в формуле (2.155) qт и Sт можно заменить на qw и Sw. Эти величины связаны со скоростью и температурой на границе вязкого подслоя:

 

 

 

λ f

(T

T );

 

 

 

w

q

w

=

 

S

 

= µ

v

,

δv

 

 

 

 

w

v

 

w

 

δv

поэтому

qw = λ f (TwµTv )Sw .

wv

Исключая Tv из равенств (2.155) и (2.156), получим

α =

 

 

Swcpf

 

 

.

w

1

+

wv

(Pr

f

1)

w

0

 

 

 

 

 

 

 

0

 

 

 

(2.156)

(2.157)

Соотношение (2.157) называют уравнением аналогии Рейнольдса: речь идет об аналогии между коэффициентом теплоотдачи α и напряжением трения Sw.

Остановимся на этом подробнее. Вспомним, что для газов число Pr 1, поэтому близко к 1 и выражение в квадратных скобках в формуле (2.157). Если теперь воспользоваться

 

 

 

237

 

 

 

 

 

 

 

 

 

 

 

 

ρ

f

w2

 

 

обычным соотношением

Sw = C f

 

 

0

,

(где Сf — коэффициент

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

трения), то формула (2.157) примет вид

 

 

St f =

 

α

 

=

 

C f

 

 

 

 

 

 

 

 

.

 

(2.158)

ρ

w c

 

 

 

2

 

 

 

f 0

pf

 

 

 

 

 

 

 

Обозначение St принято для числа Стантона (Stanton number), которое легко выразить через другие критерии подобия:

 

St f

=

Nu fx

,

 

(2.158)

 

Re fx Prf

 

 

 

 

 

 

 

откуда, при Pr 1,

 

 

 

 

 

 

 

Nu fx =

C f

Re fx Prf

C f

Re fx .

(2.159)

 

 

2

 

 

 

 

2

 

 

Формула (2.159) выражает суть аналогии Рейнольдса наиболее наглядно. Для того чтобы получить расчетные значения Nuf, а затем и α, необходимо использовать данные гидродинамики для коэффициента трения Cf (например, в виде уравнения Cf =

0,0592 Refx–0,2). При этом для газов

 

 

 

 

 

 

 

Nu fx = 0,0296Re0fx,8 .

 

 

 

(2.160)

Для капельных жидкостей соответствующая формула имеет

более сложную структуру:

 

 

 

 

 

 

 

 

 

 

 

 

 

Pr

f

0,25

 

0,8

 

0,4

 

 

 

(2.161)

Nu fx = 0,0296Re fx

Prf

Pr

 

.

 

 

 

 

 

 

w

 

Pr

f

0,25

 

 

 

 

 

 

 

 

играет здесь ту же роль,

 

Pr

 

 

Поправочный множитель

 

 

 

 

 

 

w

 

 

 

 

что и в уравнении для ламинарной части пограничного слоя (2.149). Посмотрим, как изменяется коэффициент теплоотдачи α вдоль поверхности. Согласно формуле (2.148), в ламинарной части слоя αlam ~ х–0,5, а в турбулентном пограничном слое αturb ~ х–0,2 (см. формулу (2.160)). Зависимость αtrans(x) для переходной области на

рис. 116 показана штриховыми линиями.

238

Рис. 116.

Вразд. 2.8.3 мы вывели формулу для расчета а при условии L

<xcrl. Если же длина поверхности настолько велика, что L > xcr2, то

значение α находят путем вычисления трех интегралов:

 

1

xcr1

L

xcr 2

 

 

 

 

 

 

 

 

(2.162)

α =

 

 

αlamdx +

αturbdx + αtransdx ,

L

 

0

xcr 2

xcr1

 

 

 

 

 

 

 

где αlam, αturb, αtrans — значения коэффициента теплоотдачи для ламинарного, турбулентного и переходного участков течения соответственно.

Последний интеграл мало влияет на α, так как xcrl достаточно близко к xcr2, поэтому выражение (2.162) сводят к более простому:

 

1

xcr1

L

 

 

 

 

 

 

 

(2.163)

α =

 

 

αlamdx +

αturbdx ,

L

 

0

xcr 2

 

 

 

 

 

 

 

где αlam и αturb определяют по формулам (2.148) и (2.161) соответственно.

Подставим их значения, проинтегрируем и приведем равенство (2.162) к безразмерному виду; получим

 

 

 

 

Pr

f

0,25

(2.164

0,5 0,33

0,8

0,5 0,4

 

 

 

NuLf = [0,67 Recr1 Prf

+ 0,037(ReL

Recr1 )Prf

]

Pr

 

.

)

 

 

 

 

 

w

239

В качестве критического принимают обычно Recrl Если ReL >> Recrl , то формула (2.164) упрощается:

 

 

 

Pr

f

0,25

0,8

0,4

 

 

 

NuLf = 0,037 ReL

Prf

 

 

 

.

Pr

 

 

 

 

 

w

= 5·105.

(2.165)

При расчетах конвективного теплообмена газовых потоков с

 

Pr

f

0,25

 

плоской поверхностью множитель

 

 

 

в формулах (2.164) и

Pr

 

 

 

 

 

 

w

 

(2.165) можно опустить.

Более строгий вывод уравнений подобия для конвективного теплообмена в ламинарном и турбулентном пограничных слоях требует численного интегрирования системы уравнений (2.134).

2.8.5. Теплообмен при вынужденной конвекции в трубах и каналах

Мы рассмотрели конвективный теплообмен для случая, когда жидкость ограничена твердой поверхностью только с одной стороны; такой процесс происходит, например, при течении по достаточно глубокому открытому каналу, при обтекании тела, движущегося в весьма большом объеме жидкости, и т. д. Если сечение потока ограничено со всех сторон, то характер движения и теплообмен со стенками заметно меняются. На начальном участке пограничный слой ведет себя так же, как при обтекании плоской поверхности, но затем он может достичь толщины, равной радиусу трубы или канала.

Затем уже не удастся разделить поток на невозмущенное ядро и динамический пограничный слой: силы трения проявятся по всему сечению. Поскольку трение связано с изменением скорости, следует в первую очередь установить, как меняется эпюра скоростей в зависимости от числа Рейнольдса.

В трубах режим движения зависит от числа Рейнольдса

Redf = wνd ; в качестве характерного размера выбирают диаметр трубы, а среднюю скорость задают из соотношения

240

w = GFv ,

где Gv — объемный расход жидкости; F — площадь поперечного сечения трубы (для круглых труб диаметром d F = πd2/4).

Напомним, что при Redf < 2,3·103 течение носит ламинарный характер, при Redf > 104 является турбулентным, а область 2,3·103 < Redf < 104 характеризует переходный режим течения.

При ламинарном течении жидкости по трубе диаметром 2r0 скорость w(y) на расстоянии у от оси трубы определяется

соотношением Пуазейля:

w( y) = w0 1ry 2 ,0

где w0 — скорость на оси трубы.

При турбулентном течении профиль скоростей по форме близок к трапеции, параболическое соотношение (2.166) не выполняется. Обычно используют зависимость

w = f (Redf ); w0

при Redf < 2,3·103 принимают w = 0,5 ; при больших значениях Redf w0

пользуются экспериментальными данными (рис. 117).

Рис. 117.

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