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

книги из ГПНТБ / Дорфман Л.А. Гидродинамическое сопротивление и теплоотдача вращающихся тел

.pdf
Скачиваний:
18
Добавлен:
29.10.2023
Размер:
10.79 Mб
Скачать

ГЛАВА I

ОСНОВНЫЕ УРАВНЕНИЯ

§ 1. Дифференциальные уравнения движения

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

Вдальнейшем по необходимости будут введены соответ­

ствующие усложнения в уравнениях, связанные с условиями сжимаемости и интенсивного подогрева.

Для осесимметричных движений жидкости наиболее есте­ ственно представить уравнения движения в цилиндрической

системе координат г, ср,

 

г.

 

При этом

ввиду осевой

симме­

трии производные

встречающихся в

уравнениях величин по <р

будут

 

равны

нулю.

 

 

 

 

 

 

 

 

 

 

Поэтому в случае несжимаемой вязкой жидкости при

постоянстве

физических

 

констант

р

р,

к уравнения

Навье-

Стокса

осесимметричного

 

движения

без учета массовых сил

примут вид1):

 

 

 

 

 

 

 

 

 

 

dvr

!

 

 

dvr

r

dvr

 

 

_

1

др

(

 

dt

'

 

V>'

dr

""Г”dz

 

r

 

 

ftir

1

 

 

 

 

 

 

 

. (d2vr

.

1 dvr

 

vr . d2vr

 

 

 

 

 

 

 

 

V

 

 

'

7 ~dr

T2'^^

(71)

dvv

 

 

dva ,

dva

 

v„,vr

 

 

 

 

 

 

 

 

 

 

dt

+ Vr

----- 1

~5--- - —

 

 

 

 

 

 

1

r

dr

1

2

dz

1

 

r

 

 

 

 

 

 

!)

Лойцянский

Л.

 

Г.,

Аэродинамика пограничного слоя,

Гостехиздат, 1941, стр. 52.

10

 

ОСНОВНЫЕ

УРАВНЕНИЯ

[ГЛ. I

дщ

,

dvs

_

1

др

 

 

dt

r dr

дг

~

Р

дг

 

(1-1)

 

 

 

 

 

1

дг

\

 

 

 

 

 

г

д& Г

 

 

В этих уравнениях vr,

 

и

vz

представляют собой соот­

ветственно радиальную, тангенциальную (окружную) и осевую

составляющие

вектора скорости, р — давление,

v — кинема­

тическую вязкость.

 

К системе

уравнений Навье-Стокса (1.1)

присоединим

еще уравнение неразрывности, которое для несжимаемой жидкости примет вид:

dvr

 

vr .

dvz _

(1-2)

dr

'

г -т" дг

 

Решение систем (1.1),

(1.2)

требует задания начальных

и граничных условий. Начальные условия определяются зна­

чениями

искомых величин в начальный момент

времени

t~t0,

они имеют смысл для неустановившихся движений

жидкости.

 

Граничные условия определяют поле скоростей

на гра­

ницах,

а также давление на свободной поверхности жидкости

и на поверхностях раздела. В частности, важное значение имеет такое граничное условие, как условие прилипания

жидкости к твердой стенке, т. е. отсутствие относительной скорости скольжения жидкости по поверхности твердого тела.

Для компонентов касательных напряжений в случае осе­ симметричного течения будем иметь:

§ 2]

УРАВНЕНИЯ РЕЙНОЛЬДСА

11

 

§ 2. Уравнения Рейнольдса

для осредненного турбулентного движения

Известно,

что с увеличением скорости

поток жидкости

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

ших скоростях — турбулентным.

Явления, связанные с потерей устойчивости ламинарного движения, будут рассмотрены в дальнейшем на конкретных примерах потоков вокруг вращающихся осесимметричных тел. Здесь мы лишь заметим, что в случае развитого турбу­ лентного течения уравнения (1.1), (1.2) остаются в силе только для мгновенных действительных скоростей. Однако ввиду нерегулярности мгновенных скоростей и случайности характера пульсаций уравнения (1.1), (1.2) невозможно про­ интегрировать. Поэтому Рейнольдс предложил их преобра­ зовать в уравнения осредненного турбулентного движения. Для этого действительная скорость представляется в виде суммы двух слагаемых: осредненной скорости с составляю­

щими vr,

vr

vz

и

скорости

турбулентных пульсаций

с составляющими v', v',

vr:

 

 

 

 

 

 

 

г

 

z

 

 

 

 

vr = vr-j-v'r,

v,f = v.fy-uf,

 

(1.4)

Сделаем

подстановку

этих соотношений в уравнения (1.1),

(1.2) и произведем

их

осреднение во времени (обозначается

чертой сверху), полагая,

что

операция

осреднения

обладает

следующими свойствами:

 

 

 

 

 

 

 

т'^0’

 

 

 

4г=4Ь ■

(1-5)

где $ — одно

из

переменных

г,

ср, г,

t; ср, ф— любые две

функции этих

переменных.

 

 

 

 

Тогда получим

следующую

систему

уравнений осреднен­

ного осесимметричного движения вязкой несжимаемой жидко­ сти— уравнения Рейнольдса — в цилиндрической системе

12

 

 

 

 

 

ОСНОВНЫЕ

УРАВНЕНИЯ

[ГЛ. I

координат:

 

 

 

 

 

 

 

 

 

 

 

д^г I

- <4 ,

- д”г

 

 

 

 

 

 

1

dt 1

г

dr

1

г

 

dz

 

г

 

 

 

 

 

__

1 др ,

(d2vr

,

1 dvr

 

vr

,

d2vr\

 

p

dr

 

 

 

dr2

'

r

dr

 

r2

 

dz2 )

dvv

-

dvv

 

-

 

dv„

 

vrv„

 

 

 

 

~dT^!~Vr~dr

1“Vz ~дГ + ~~Г~

 

 

 

 

 

 

 

d‘-v,f

(

1

dVy

v,f

'

d2v^

 

 

 

 

dr2

'

r

dr

r2

 

dz2 /

(1-6)

 

 

 

 

 

 

0

-----

 

d

____

v v

 

 

 

 

dr \ r v/

 

dz \ ? г/ — 2 rr <p ,

 

 

1

dp

\

 

\ dr2

I

1

dr

I

d2ve \

 

 

p

dz

'

 

'

r

 

dz2

/

dvr ^Vr'dvz—r.

 

 

 

 

 

 

 

dr

r

 

dz

~

 

 

 

 

 

 

 

Сопоставляя систему уравнений (1.6) с уравнениями (1.1),

замечаем,

что

к

вязким

членам

уравнений

(1.1) добавляются

члены,

соответствующие

турбулентным

напряжениям трения,

возникающим благодаря турбулентным пульсациям скорости.

§ 3. Уравнение баланса энергии

Не останавливаясь на выводе1), приведем уравнение баланса энергии для осесимметричного потока несжимаемой жидкости:

+Ч2(»’+ 2Ш+2(»3+

,

(д^

v^2 ,

(dv^2 '

(dvs

 

 

\ dr

г )

dz )

\ dr + dz ] J ’

(1 -7>

1)Лойцянский Л. Г., Аэродинамика пограничного

слоя,

Гостехиздат,

1941,

стр. 53.

 

 

 

§ 3]

УРАВНЕНИЕ БАЛАНСА ЭНЕРГИИ

13

где

Т — температура, ср — теплоемкость при

постоянном

давлении. Член в фигурных скобках соответствует диссипации (рассеиванию) энергии в результате трения.

Вслучае турбулентного режима можно из уравнения (1.7) получить соответствующее осредненное уравнение баланса энергии.

Без учета диссипации оно примет вид:

Врассматриваемом случае постоянства физических кон­ стант р, ср, р, X уравнение баланса энергии можно решить,

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

ГЛАВА II

ЛАМИНАРНОЕ ТЕЧЕНИЕ ВОКРУГ ВРАЩАЮЩЕГОСЯ

ВСВОБОДНОМ ПРОСТРАНСТВЕ ДИСКА

§4. Ламинарное движение, создаваемое

вращающимся диском.

Точное решение уравнений Навье-Стокса

 

Пусть

бесконечная

плоская

пластинка

z = 0 вращается

в

вязкой жидкости вокруг

оси

г = 0 с постоянной угловой

скоростью

(О.

 

 

 

 

 

 

 

 

 

Рассматриваем движение жидкости в полупространстве

2^>0. Граничные условия задачи будут иметь вид:

 

 

■пг = 0,

 

и. = 0 при

2 = 0,

}

(2-1)

 

 

 

 

 

 

 

 

 

 

Vr — Q,

Uf = 0

 

При

2 = 00.

 

 

Скорость vz

при

2 = оо

не

может быть равна нулю, так

как диск

действует

в

качестве центробежного вентилятора

и

вызывает подсос,

 

создавая

отрицательное

значение уг

в

бесконечности

и

радиальное

движение

от центра,

в осо­

бенности вблизи

диска.

 

 

 

 

 

 

Оценим

вначале

толщину S

слоя жидкости,

увлекаемого

диском вследствие трения1). Пусть направление, в котором вдоль диска скользит поток и которое параллельно каса­ тельному напряжению на стенке тот, образует с окружным направлением угол &. Тогда радиальная составляющая каса­ тельного напряжения равна тот sin &. Она уравновешивается центробежной силой отбрасываемого пластинкой потока, т. е. пропорциональна ргш28. С другой стороны, окружная соста­ вляющая касательного напряжения totcos& пропорциональна

Ч Прандтль Л., Гидроаэромеханика, перев. с нем. ИЛ, 1949, стр. 436.

§ 4]

ТОЧНОЕ РЕШЕНИЕ

УРАВНЕНИЙ

НАВЬЕ-СТОКСА

15

градиенту

окружной составляющей

скорости около стенки,

т. е.

(— — знак пропорциональности). Исключая каса­

тельное напряжение тст из

соотношений

 

 

тст sin & — pro)23,

j

(2.2)

 

 

п

У

j

 

Jot COS it — U -g-

 

и принимая, что угол О

не зависит от радиуса, что

под­

тверждается наблюдениями,

мы получим:

 

 

8

 

 

 

(2-3)

следовательно, толщина увлекаемого диском слоя одинакова

на всех радиусах.

Отсюда следует, что касательное напряжение

тот — ргш28 — ргш j/vio,

а момент сопротивления диска радиуса R, который пропор­

ционален

произведению

касательного напряжения,

площади

и плеча,

будет

~ р/?4ш "j/w.

(2.4)

 

Л1 ~

Из этих рассуждений следует, что для интегрирования уравнений (1.1), (1.2) целесообразно ввести вместо z без­

размерное расстояние от стенки С.—

т. е. принять

v

(2-5)

Очевидно, что составляющие vr и vv пропорциональны гш, где множители пропорциональности есть функции С, а соста­

вляющая

вследствие уравнения неразрывности должна быть

пропорциональна 8w ='|/Ча>.

 

 

Таким образом, целесообразно сделать следующую замену

переменных *):

 

 

 

•ur = r^F^\

^ = rt»G(Q,

vs = y^H^), |

(2 6)

 

 

р = — pvuP (С).

/

 

!) К а г m a n Т h.,

Laminate und turbulente Reibung,

ZAMM,

16

ЛАМИНАРНОЕ ТЕЧЕНИЕ ВОКРУГ ДИСКА

[ГЛ. II

Подстановка (2.6) приводит уравнения (1.1), (1.2) к виду (штрихами обозначены производные по С):

F2 —G2+F'H=FZ'

2FG + G'H=G",

(2.7)

НН' = Р'-^-Н", 2F + H' = Q,

а граничные условия будут:

F = 0,

0=1,

Н = 0 при

F = 0,

0 = 0

(2-8)

при

Заметим, что при такой замене вязкие члены в правой части уравнений (1.1) сводятся только к последним слагаемым, представляющим производные по направлению, нормальному

к плоскости

диска, так что уравнения движения тождественны,

уравнениям пограничного слоя.

 

 

Первые два уравнения (2.7) совместно с последним дают

возможность

определить F,

О,

И,

а третье — значения Р.

Если И— с

при С-> оо,

то можно формально построить!)

разложения по

степеням

для

F,

О и Н, удовлетворяю­

щие дифференциальным уравнениям и условиям на бесконеч­ ности. Нетрудно определить, что первыми членами разло­ жений будут

 

 

Ai + B2

е-2Л

Л(^+В2)

_3ес _

 

 

2с2

 

4с<

 

 

 

Д 2 _|_ Д2

 

 

 

 

--Л4^(17Л2+В2)е-4О:+ •••’

О = Ве-^

В(42-г№) с_3л, ,

 

12В

 

 

(2.9).

 

 

 

 

 

 

 

 

 

 

 

 

'

18c6

-!-•••-

Н = — с 4------ е-с-

2c3

 

 

 

1

с

 

 

т

6c5

°

288c7

(17Л2 +fi2) e

4e’ + • • ■

Остается определить

только

значения постоянных А, В и с.

9 Cochran W. Q., The flow due to rotating disc, Proc. Cambr. Philos. Soc., t. 30, № 3, 1934.

§ 41

 

ТОЧНОЕ

РЕШЕНИЕ УРАВНЕНИЙ

НАВЬЕ-СТОКСА

 

17

 

С другой стороны, можно получить

формальное разло­

жение

вблизи

С = 0,

удовлетворяющее

(2.7)

и

граничным

условиям

при £ = 0:

 

 

 

 

 

 

 

 

 

 

 

 

Г,

у

 

1

r2

 

1

,

уЗ

1 ,'iyi

 

1

у5

,

 

 

1

 

Т7

■ а.(у,

 

С

 

-g

 

 

bo-,

-gQ- п0,

 

 

 

 

 

 

 

(1

 

arfin

\

/

bn

I

апЬп

\

Г7

I

 

 

 

 

 

 

________re

I

______ ®_

0 0

 

 

 

 

 

 

360

 

90

/

’К 315

 

1260 /

 

 

 

 

о=1+*ос+4

 

 

 

(«<л - щ4 - 4 v—

(2.10)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f ао

|

^0

ув

|

(

ао

60

 

й0^0 \ р |

 

 

 

 

 

\ 90

'

45 / ’

'

\ 315

315

 

252 /

 

 

 

Н = -а^.

 

 

 

 

 

'О’* + W

 

 

 

 

 

 

 

 

где

а0,

^ — неизвестные постоянные.

находятся из условия

 

Значения постоянных А, В, с, а0, Ьо

непрерывного

сращивания значений

функций F, О и Н и

их

производных,

получаемых

из

разложений (2.9)

и

(2.10).

Если задаться

приближенными значениями

а0

и

Ьо,

то

с по­

мощью

(2.10)

находятся величины

F,

G и Н и их

произ­

водные для малых значений С Затем методом Адамса можно

продолжить

решения

дифференциальных

уравнений

(2.7)

к большим значениям С. Сравнением полученных значений

функций и их производных в

точках С= 1,9 и

2,5 со зна­

чениями, которые даются разложениями (2.9), Кокрэн1) чис­

ленным путем уточнил величины постоянных а0 и Ьо и нашел

значения с,

А и В. Его результаты

дают:

 

 

 

 

 

 

 

 

ао = О,51О; Ьо = — 0,616;

 

 

 

 

 

(2.И)

 

 

 

 

с = 0,886;

А — 0,934;

В = 1,208.

 

 

 

Значения функций F, О, И, Н', G' и Р,

вычисленные

Кокрэном, приводятся в таблице 1. Первые три функции

графически представлены на рис. 1.

Из графиков видно, что

расстояние от стенки, на котором окружная скорость течения

понижается до половины окружной скорости стенки, равно

§0,5 ~ j/v/d).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4) См.

сноску на

стр.

16.

ГОС.

ПУБЛИЧНА

 

I’-X

 

 

 

 

 

 

 

 

 

НАУЧН-ТЕХНЙЧЕСИАЯгУ«/|

2 Зак. 944. Л. А. Дорфман

БИБЛИОТЕКА СССР |

18

ЛАМИНАРНОЕ ТЕЧЕНИЕ ВОКРУГ ДИСКА

[ГЛ. II

Таблица 1

Значения функций, определяющих распределение скоростей и давления вблизи диска,

 

 

 

вращающегося в неподвижной жидкости

 

 

 

 

 

(по В.

Кокрэну)

 

 

 

_

1/ ш

F

О

—р

F'

-G'

C=z у

 

0

 

0

1,000

0

0

0,510

0,616

 

0,1

 

0,046

0,939

0,005

0,092

0,416

0,611

 

0,2

 

0,084

0,878

0,018

0,167

0,334

0,599

 

0,3

 

0,114

0,819

0,038

0,228

0,262

0,580

 

0,4

 

0,136

0,762

0,063

0,275

0,200

0,558

 

0,5

 

0,154

0,708

0,092

0,312

0,147

0,532

 

0,6

 

0,166

0,656

0.124

0,340

0,102

0,505

 

0,7

 

0,174

0,607

0,158

0,361

0,063

0,476

 

0,8

 

0,179

0,561

0,193

0,377

0,032

0,448

 

0,9

 

0,181

0,517

0,230

0,388

0,006

0,419

 

1,0

 

0,180

0,468

0,266

0,395

—0,016

0,391

 

1,1

 

0,177

0,439

0,301

0,400

—0,033

0,364

 

1,2

 

0,173

0,404

0,336

0,403

—0,046

0,338

 

1,3

 

0,168

0,371

0,371

0,405

—0,057

0,313

 

1,4

 

0,162

0,341

0,404

0,406

—0,064

0,290

 

1,5

 

0,156

0,313

0,435

0,406

—0,070

0,268

 

1,6

 

0,148

0,288

0,466

0,405

—0,073

0,247

 

1,7

 

0,141

0,264

0,495

0,404

—0,075

0,228

 

1,8

 

0,133

0,242

0,522

0,403

—0,076

0,210

 

1,9

 

0,126

0,222

0,548

0,402

—0,075

0,193

 

2,0

 

0,118

0,203

0,572

0,401

—0,074

0,177

 

2,1

 

0,111

0,186

0,596

0,399

—0,072

0,163

 

2,2

 

0,104

0,171

0,617

0,398

—0,070

0,150

 

2,3

 

0,097

0,156

0,637

0,397

—0,067

0,137

 

2,4

 

0,091

0,143

0,656

0,396

-0,065

0,126

 

2,5

 

0,084

0,131

0,674

0,395

-0,061

0,116

 

2,6

-

0,078

0,120

0,690

0,395

—0,058

0,106

 

2,8

0,068

0,101

0,721

0395

—0,052

0,089

 

3,0

 

0,058

0,083

0,746

0,395

—0,046

0,075

 

3,2

 

0,050

0,071

0,768

0,395

—0,040

0,063

 

3,4

 

0,042

0,059

0,786

0,394

—0,035

0,053

 

3,6

 

0,036

0,050

0,802

0,394

—0,030

0,044

 

3,8

 

0,031

0,042

0,815

0,393

—0,025

0,037

 

4,0

 

0,026

0,035

0,826

0,393

—0,022

0,031

 

4,2

 

0,022

0,029

0,836

0,393

—0,019

0,026

 

4,4

 

.0,018

0,024

0,844

0,393

—0,016

0,022

 

оо

 

0

0

0,886

0,393

0

0

Соседние файлы в папке книги из ГПНТБ