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

книги / Общая физика. Гидродинамика и теплообмен

.pdf
Скачиваний:
2
Добавлен:
12.11.2023
Размер:
2.01 Mб
Скачать

Интеграл последнего уравнения

δ2

=

140

 

ν

x .

(6.18)

 

 

 

2

 

13 v0

 

Запишем в окончательном виде выражение длязависимости δ(x) :

δ =

280

 

ν

x 4,64

ν

x .

(6.19)

13 v0

 

 

 

v0

 

Формуле (6.19) можно придать безразмерный вид

δ =

4,64

=

4,64

,

(6.20)

 

 

x

v0 x

 

Rex

 

ν

где Rex − число Рейнольдса для течения в пограничном слое на рас-

стоянии х от края пластинки.

Далее посчитаем напряжение вязкого трения на поверхности пластинки:

(τтр )

 

=

3

 

µv0

=

3

 

µv0 13v0

= 0,323

µ2v03

= 0,323

µρv03

. (6.21)

 

 

 

δ

 

 

νx

x

 

c

2

 

2

 

280νx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сравним выражение (6.21) с известной [3] «формулой Блазиуса», полученной иным способом решения уравнений пограничного слоя для плоской пластины:

(τтр )

c

= 0,332 µρv03 .

(6.22)

 

x

 

 

 

 

Видим, что совпадение хорошее, так как расхождение результатов (небольшое различие в численном коэффициенте) непревышает 3 %.

Вместо напряжения вязкого трения можно ввести так называемый местный коэффициент сопротивления

 

 

=

(τтр )

 

C

 

c

,

(6.23)

f

 

 

 

ρv02 2

 

 

 

 

 

 

71

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

Подставим в (6.23) выражение (6.21) и получим

C f

= 0,646

ν

=

0,646

.

(6.24)

v0 x

 

 

 

 

Rex

 

Взадаче Блазиуса численный коэффициент имеет значение 0,664. Наконец, сосчитаем полное (по обеим сторонам пластинки) со-

противление трения пластинки длиной l (вдоль оси 0х)

 

 

Wтр

= 2l (τтр )cdx = 0,646

µρv03 l

dx

.

 

 

 

 

 

0

0

x

Так как l

dx

= 2

l , то

 

 

 

 

 

 

 

0

x

 

 

 

 

 

 

 

W = 1,292 µρlv3 .

(6.25)

 

 

 

тр

0

 

 

Физический смысл последней величины – это сила, приходящаяся на единицу длины (вдоль оси 0z) двусторонней пластины. Размерность Wтр – Н/м.

6.5. Тепловой пограничный слой

Аналогично гидродинамическому пограничному слою Г.Н. Кружилиным было введено понятие теплового пограничного слоя – слоя жидкости у твердой стенки, в пределах которого температура жидкости изменяется от температуры стенки Tc до температуры жидкости

T0 вдали от стенки.

 

 

 

 

Обозначим символом k тол-

 

щину

теплового

пограничного

 

слоя (рис. 6.2). В общем случае k

 

и δ не совпадают, но являются

 

величинами одного порядка ма-

 

лости

k = O(δ) .

Считаем, что

 

выполняется условие k << l , где

Рис. 6.2

l – длина пластины.

72

 

 

 

Пусть температура набегающего со скоростью v0 потока постоянна и равна T0 . Все изменение температуры жидкости сосредоточено в сравнительно тонком слое, прилегающем к поверхности тела. Внутри этого слоя ( y < k) Ty 0 , на границе и вне слоя ( y k) Ty = 0 и T = T0 . Кроме того, можно показать (см. подразд. 6.1), что выполняется условие

2T << 2T , x2 y2

которое означает, что теплопроводность вдоль слоя мала по сравнению с теплопроводностью поперек слоя. По аналогии с (6.3), после оценки порядка малости величин, стационарное уравнение энергии можно записать в виде

v

x

T

+ v

y

T = a

2T .

(6.26)

 

x

 

y

 

y2

 

Используя уравнение

Фурье

qy

= −λ T ,

правую часть (6.26)

 

 

 

 

 

 

 

y

 

можно преобразовать следующим образом:

 

2T

=

λ ∂

T

= −

1

a

 

 

 

 

 

 

 

 

y

2

 

 

ρc

 

 

 

ρc y

y

 

qy , y

и подставить в исходное уравнение (6.26). В итоге получим стационарное уравнение ламинарного теплового пограничного слоя

vx

T

+ vy

T

= −

1

 

qy

.

(6.27)

x

y

ρc

y

 

 

 

 

 

Аналогично (6.9) можно вывести интегральное уравнение теплового пограничного слоя

d

k vx (T0

T )dy = −

qc

,

(6.28)

dx

 

0

 

ρc

 

где qc = qy y=0 .

Уравнение (6.28) впервые было получено Г.Н. Кружилиным.

73

6.6. Аппроксимация профиля температуры

Примем температуру поверхности тела постоянной: Tc = const (не зависит от х). Для удобства температуру жидкости будем отсчитывать от значения Tc . Введем так называемые температурные напоры:

Θ = T Tc ;

Θ0 = T0 Tc .

(6.29)

Кривую температурного профиля можно описать полиномом третьей степени

Θ = a + by + cy2 + dy3 .

(6.30)

Значения коэффициентов a,b,c,d находятся из следующих граничных условий:

1)y = 0; Θ = 0 , так как T = Tc . Из этого граничного условия и (6.30) получим коэффициент a = 0;

2)y = 0; 2Θy2 = 0 . Это условие следует из граничного усло-

вия прилипания для скорости: vx = 0 ; vy = 0 при y = 0 − и из уравне-

ний (6.26), (6.29). Тогда c = 0;

3)y = k; Θ = Θ0 – на внешней границе пограничного слоя;

4)y = k; ∂Θy = 0 . Это условие означает, что профиль темпера-

туры на границе слоя практически касается нормали к стенке. В результате, аналогично (6.14), получим

Θ

=

3

 

y

1

 

y 3

 

 

 

 

 

 

 

 

 

.

(6.31)

Θ0

 

 

 

 

 

2

k

 

2

k

 

6.7. Толщина теплового пограничного слоя

Так как k ≤ δ , то интегрирование (6.28) в пределах от y = 0 до y = k является интегрированием в пределах как теплового, так и гидродинамического пограничных слоев.

Интеграл левой части уравнения (6.28) с использованием (6.31), (6.14) примет следующий вид:

74

k vx (T0 T )dy =k vx (Θ0 − Θ)dy =

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

3

y

 

 

1

 

 

y 3

 

 

 

 

3 y

 

1

y 3

 

 

 

 

 

 

= v0Θ0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

+

 

 

 

 

 

 

dy

=

 

 

 

 

 

δ

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

2

 

 

 

 

 

δ

 

 

 

 

 

2 k

 

2 k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 k2

 

 

1 k 4

 

 

 

9 k3

 

 

 

3 k5

 

3 k5

 

 

 

1 k7

 

= v0

Θ0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

δ

 

 

8

 

δ

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

20 k

3

δ

28 k

3

δ

3

 

4

 

 

 

 

 

 

 

 

12 δk 20 kδ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 k 2

 

 

 

 

3 k 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= v0Θ0

δ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20

δ

 

 

 

 

280

 

δ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как kδ < 1, то второй член выражения мал по сравнению с первым, поэтому им можно пренебречь. В итоге можно записать

k

3

 

 

k 2

 

vx (T0 T )dy =

 

v0

Θ0

δ

 

.

(6.32)

20

 

0

 

 

 

δ

 

Преобразуем правую часть уравнения (6.28):

 

qc

 

λ ∂Θ

 

 

 

3

 

3y

2

 

 

3

 

aΘ0

 

 

 

 

 

 

 

 

 

 

 

=

 

 

= aΘ0

 

 

=

 

.

(6.33)

 

 

 

 

 

 

2k

3

 

 

 

ρc

ρc y

 

y=0

2k

 

 

y=0

 

2 k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Объединим (6.32) и (6.33) в интегральное уравнение теплового пограничного слоя (6.28) и получим

3

 

 

 

d

 

 

k

2

 

 

 

3 aΘ0

 

 

v0

Θ0

 

 

 

δ

 

 

 

 

=

 

 

 

 

 

 

,

20

 

 

 

δ

2

 

 

k

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

k 2

 

10 a

 

 

 

 

 

δ

 

 

 

=

 

 

 

 

 

.

 

(6.34)

 

 

 

 

 

 

 

 

 

 

 

 

dx

δ

 

 

 

v0

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Введем новую переменную γ − отношение толщин теплового и гидродинамического пограничных слоев:

γ =

k

< 1.

(6.35)

δ

 

 

 

 

 

 

75

Исходя из полной аналогии гидродинамического и теплового пограничных слоев (уравнения, граничные условия, распределения скорости и температуры), можно полагать, что δ и k зависят от координаты x одинаково, тогда их отношение (6.35) не будет зависеть от х, т.е. можно полагать, что γ = const. Перепишем уравнение (6.34) с использованием (6.35):

d

(δγ2 ) =

10

 

a

.

 

 

 

dx

 

v0 γδ

Продифференцируем левую часть последнего уравнения:

2γδ dγ + γ2 dδ = 10 a . dx dx v0 γδ

Первое слагаемое левой части уравнения обнуляется. Тогда остается

γ3δ

dδ

=

10a

.

(6.36)

 

 

 

dx v0

 

Используем ранее полученное выражение (6.17):

δdδ = 140ν dx 13v0

и получим для (6.36)

γ3 =

13a

=

13

 

; Pr = ν .

 

 

 

 

 

 

14ν

14Pr

 

 

a

 

В итоге получим

 

 

 

 

 

 

 

γ =

0,976

 

1

.

(6.37)

 

 

 

 

 

3 Pr

 

3 Pr

 

Рассмотрим некоторые частные случаи.

1. Для капельных жидкостей Pr 1, например, для воды при 20 ºС число Прандтля Pr = 7, тогда по формуле (6.37) получим γ ≈ 0,5 < 1, т.е.

условие (6.35) выполняется.

76

2. Для воздуха в интервале температур от 0 до 90 ºС число Прандтля Pr 0,7 , тогда γ ≈ 1,1. Получается, что k ≥ δ , хотя разница в тол-

щинахгидродинамического итепловогопограничных слоев невелика.

3. Для

жидких металлов

Pr << 1, например, для ртути

Pr = 2,72 102

при 20 ºС, тогда

γ ≈ 3,2 и изложенная выше теория

непригодна.

Подставим полученное ранее выражение (6.19) для δ в (6.37) и получим для толщины теплового пограничного слоя

k =

4,53

 

ν

x ,

(6.38)

 

 

 

3 Pr

 

v0

 

т.е. k зависит от x по закону k ~

x .

 

 

 

Еще одно соотношение для k можно вывести из ранее получен-

ной формулы (6.20) и (6.37):

 

 

 

 

 

k =

 

4,53 x

 

 

 

 

 

,

(6.39)

 

Re

3

 

 

 

Pr

 

 

 

 

x

 

 

 

где Rex – число Рейнольдса для течения в пограничном слое на расстоянии х от края пластинки.

6.8. Коэффициент теплоотдачи

Полученные соотношения позволяют оценить коэффициент теплоотдачи α с учетом ламинарного пограничного слоя на пластинке. Плотность теплового потока на стенке может быть определена из уравнения Фурье

qc

 

T

= −λ

.

 

 

y y=0

Теплоотдача от стенки к жидкости определяется граничным условием III рода – законом Ньютона−Рихмана:

qc = α(Tc T0 ) .

77

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

α = −

 

λ

T

 

 

 

 

.

(6.40)

 

 

Tc

T0

y y=0

 

Внаших новых обозначениях (6.29) выражение (6.40) примет вид

 

 

λ

∂Θ

 

λ

 

3 Θ0

 

3 λ

 

 

α =

 

 

y y=0 =

 

 

 

 

=

 

k .

 

 

Θ0

Θ0

 

2 k

2

 

Таким образом, можно записать, что α ~ 1 k ~ 1

x .

Перепишем (6.41) с учетом (6.39):

 

 

 

 

 

 

 

 

3λ

3

 

 

 

 

λ

 

3

 

α =

 

 

Rex

Pr

= 0,331 x

 

Rex

Pr .

2 4,53 x

 

Введем в рассмотрение безразмерное число Нуссельта

Nux = αλx .

(6.41)

(6.42)

(6.43)

Тогда (6.42) с учетом (6.43) после несложных преобразований можно записать в общепринятом виде

Nux = 0,33Re0,5x Pr1 3 .

(6.44)

Контрольные вопросы

1.Что такое пограничный слой?

2.Вывести уравнения Прандтля для пограничного слоя.

3.Вывести интегральное уравнение пограничного слоя.

4.Аппроксимация профиля скорости в приближении Польгаузена.

5.Получить выражение для толщины пограничного слоя на плоской пластине.

6.Дать определение местного коэффициента трения.

7.Что такое тепловой пограничный слой?

8.Аппроксимация профиля температуры.

9.Получить выражение длятолщинытеплового пограничного слоя.

10.Вывести критериальное уравнение для числа Нуссельта.

78

7.ТУРБУЛЕНТНОСТЬ

7.1.Гидродинамическая устойчивость

Имеются два основных режима течения жидкости: ламинарный и турбулентный. При ламинарном режиме частицы жидкости движутся упорядоченно, без перемешивания, слоисто (например, рассмотренное ранее течение Пуазейля в круглой трубе). При других условиях на упорядоченное движение жидких частиц накладывается хаотическое, пульсационное движение, приводящее к разрушению слоистой структуры и перемешиванию слоев. В этом случае направления и значения скоростей отдельных частиц беспрестанно меняются. Такой режим течения называется турбулентным.

Возникновение турбулентного режима связано с потерей устойчивости движения жидкости, которая наступает при определенных значениях гидродинамических параметров. Если в ламинарном потоке малые возмущения затухают и не приводят к изменению его структуры, то такой поток является устойчивым. Малые возмущения возникают, например, на входе жидкости в трубу из резервуара – «неплавность» входа в трубу.

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

Впервые научные исследования возникновения турбулентного движения проводил в 1883 г. английский физик О. Рейнольдс. Он изучал поведение подкрашенной струйки в потоке жидкости, текущей в круглой прозрачной цилиндрической трубе (рис. 7.1).

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

79

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

Рис. 7.1

Устойчивый ламинарный режим в круглой трубе существует при числах Рейнольдса Re Reкр , где Reкр 2300 . В этом случае в потоке

затухают любые внешние малые возмущения. При Re > Reкр в зависи-

мости от условий может существовать либо ламинарный, либо турбулентный режим. При Re >> Reкр реализуется развитаятурбулентность.

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

минарном потоке периодически появляются кратковременные очаги турбулентности, которые могут на отдельных участках заполнять все сечение потока (турбулентные пробки). Другими словами, при одном и том же значении Re в одном и том же сечении трубы может происходить чередование ламинарных и турбулентных режимов – так называемая перемежаемость режимов. Ее количественная характеристика – коэффициент перемежаемости γ − представляет собой долю времени ∆t существования турбулентного режима в данном сечении трубы в течение некоторого времени t: γ = ∆tt .

80