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

книги из ГПНТБ / Картвелишвили Н.А. Потоки в недеформируемых руслах

.pdf
Скачиваний:
7
Добавлен:
25.10.2023
Размер:
12.62 Mб
Скачать

Нетрудно проверить, что

_ g r a d v 2 = ( v , V ) v + [ v , r o t v ] .

Используя эту формулу и полагая, что внешние силы имеют потенциал, т. е. что F = —grad У, а жидкость имеет постоянную плотность, приведем уравнение (11.1) к виду

- | f + g r a d r = [ v , r o t v ] .

(11.2)

Здесь

Для установившегося движения уравнение (11.2) имеет вид

g r a d r = j v , rotv] .

(11.3)

Так как векторное произведение двух векторов есть вектор, нор­

мальный к обоим векторам, то из

(11.3) следует,

что вектор

grad Г перпендикулярен к плоскости,

в которой лежат

векторы v

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

ния, т. е. траектория rotv, лежат

на поверхности

Г = const, т. е. Г==

= const вдоль линии тока; этот

вывод

был

получен

в §

8 иным

путем. Но с переходом от одной

траектории

к

другой

Г,

конечно,

в общем случае изменяется.

Чтобы

Г

было

постоянно

во всей

массе жидкости, как видно из

 

(11.3),

должно

быть

[v, rotv] = 0,

т. е. rotv = 0, что соответствует

безвихревому

движению, либо век­

торы v и rot v должны иметь одинаковое направление. Этот послед­ ний случай — единственное вихревое движение с одинаковым зна­ чением Г для всей массы жидкости. Такое движение называется однородно-винтовым. Это движение в прямоугольном канале впер­ вые было исследовано И. С. Громеко [20], значительное внимание уделено ему А. Я. Миловичем [41]. М. В. Потапов [47] обратил внимание на то, что однородно-винтовое движение характеризуется резко неравномерным распределением скоростей по сечению по­ тока, очень далеким от реально наблюдаемых распределений, ко­ торые в призматических руслах более или менее близки к равно­ мерным (поступательная скорость одинакова во всех точках попе­ речного сечения). Отступление от равномерного распределения есть следствие главным образом гидравлических сопротивлений, кото­ рые приводят к тому, что максимальная в сечении скорость превы­ шает среднюю только на 20—30%, а не в 2,5 раза, как в одно­ родно-винтовом движении. В связи с этим М. В. Потапов рассмат­ ривает иное, так называемое продольно-винтовое движение, с вих­ ревыми линиями, параллельными дну и стенкам канала, которое значительно ближе к движениям реальных жидкостей, чем однород­ но-винтовое.

Уравнение (11.3) в проекциях на координатные оси записыва­ ется так:

 

— 5Р=2( - іто - С<а),

(11.4)

 

 

(11.5)

 

 

(11.6)

Здесь s,

т], £— проекции на оси х, у, z компонентов

угловой

скорости —

rot v:

 

Присоединим

к написанным

уравнениям

р и с

^

Поперечная

также уравнение неразрывности

 

циркуляция

в

прямо-

ди

, до

, dw

п

/ 1 1 0 Ч

Уольном канале

(один

дх

ду

dz

=0.

(11.8)

 

винт,

т=п=1).

Рассмотрим прямоугольный канал с горизонтальным дном; на­

правим оси координат как показано на

рис. 12. Так как течение

не изменяется вдоль оси х,

то

 

 

 

 

дТ

ди

ди

dw

^

п л п .

дх

дх

дх

дх

=0,

(,11. У)

 

а так как вихревые линии

параллельны

оси х, то т] = £=0.

Урав­

нение (П.4) при этом удовлетворяется тождественно, уравнения (11.5) и (11.6) приобретают вид:

дТ

дТ

(11.10)

dv

dz

 

а из уравнений (11.7) и (11.9) следует

 

— = 0

да

=0.

 

dz

и'

 

Таким образом, и = const, т. е. продольная составляющая ско­

рости

одинакова во всех точках потока. Уравнение же неразрыв­

ности

(11.9) приобретает вид

 

dz -0.

 

 

ду

'

(11.11)

 

d v

 

d w

 

Продифференцируем

первое уравнение (11.10) по z,

а второе

по у и вычтем второе

уравнение из первого. Учитывая

при этом

уравнение

(11.11), получаем

 

 

 

или

 

 

 

 

 

где

 

(W, gradS)=0,

 

 

(11.13)

 

W^Vrf+w1

 

 

(11.14)

 

 

 

 

— скорость поперечной

циркуляции.

 

 

 

Первое

уравнение

(11.7) и уравнения (11.11)

и (11.12)

обра­

зуют систему трех уравнений с тремя неизвестными v,

w,

|. Эта

система нелинейна, так как нелинейно уравнение

(11.12)

и в квад­

ратурах она в общем виде неразрешима. По М. В. Потапову, мы

ограничимся

частным случаем, в

котором

квадратура

возможна.

К этому случаю мы придем,

написав уравнение (11.12)

в виде

 

 

 

V

 

 

 

w

 

 

 

 

 

 

 

 

 

ді

 

~~

~~дТ~ •

 

 

 

 

 

 

 

 

дг

 

 

 

ду

 

 

 

 

 

 

Обозначая каждое из этих отношений через К, получаем

 

 

 

 

 

w =

- X - g - .

 

 

(11.15)

Частный

случай,

который

мы

имеем

в

виду,

есть

случай ?» =

= const>0.

Полагая

А, = 2/А2,

уравнения

(11.15)

запишем так:

 

 

дї

 

k2

 

di

 

k2

 

 

 

/111 r\

 

~W=~

W>

-dr=Tv-

 

 

 

 

( l u 6 >

Дифференцируя первое уравнение

(11.16) по у,

а второе по z

и складывая, получаем,

учитывая

первое

уравнение

(11.7),

 

_дЧ_ ,_дЧ_

 

k2 (I dv

 

dw \

, 9 t

 

 

 

ду2

дг2

 

 

 

 

 

 

 

 

 

 

или

 

 

ду* ~'

 

' - ^ = 0 .

 

 

 

 

(11.17)

 

 

 

дг2

 

 

 

 

 

 

 

д Ч

 

д Ч

 

 

 

 

 

 

 

 

Таким образом, задача свелась к системе трех линейных урав­

нений (11.16) и (11.17).

 

 

 

 

 

 

 

 

 

 

 

Положим

% = YZ,

где Y зависит

только

от у,

a

Z

только от z.

Подставляя это значение | в уравнение

(11.17), получаем

 

 

 

- V

+

4 - + # = 0 .

 

 

 

(11.18)

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

 

 

 

 

 

 

 

 

 

 

 

 

 

^у-^-о2,

 

 

^

=

- р

2

,

 

 

(11.19)

где а

н (З — вещественные

постоянные,

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

условию

 

 

 

 

а2 +

р2=АЗ_

 

 

 

 

(11.20)

Из уравнений (11.19)

 

 

 

 

 

 

 

 

 

 

 

 

Y =

С[ cos ay -(- С2

sin ay,

 

 

 

 

Z = D ,

cos p z + D 2

sin pz,

 

 

где Сі, Co, Di, A>-—произвольные

постоянные. Найденное

решение

должно

удовлетворять очевидным

краевым

условиям:

и = 0 при

у = 0

и

у = Ь; до = 0

при

z = 0,

z = /i.

Подставляя

\=YZ

в (11.16), находим из этих условий

 

 

 

 

 

 

 

 

С , = Д = 0 ,

sina6=0,

 

sinp/z=0,

 

 

откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

* = • —

,

Р =

— .

 

 

 

где /п и п — произвольные

 

целые

числа,

положительные

или отри­

цательные. Следовательно,

 

 

 

 

 

 

 

 

 

 

 

i j = = C s i n ^ ^ s i n - ^ ,

 

C = O D 2 ,

 

(11.21)

 

 

2Сля .

 

/пяу

 

ляг

 

,., w n m

 

 

v = = ~ M ~ s m

 

Ъ C 0 S _ A '

 

( 1 L 2 2 >

 

 

™ = - П Р Ї - С 0 8

- Т - 8 Ш

Т -

 

О 1 - 2 3 )

 

 

£2

 

- 2 ( 4 + 4 ) -

 

 

(П.24)

Если вместо выражений

(П.19)

взять

 

 

 

 

 

 

Г II

 

 

"711

 

 

 

 

 

 

 

Г

 

'

 

Z

- = р

или

 

 

 

 

Y"

 

.,

 

Z"

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(условие (11.20) в обоих случаях может быть удовлетворено), или если рассматривать случай К= const < 0, т. е. k2 < 0, то функции Y и Z (или по крайней мере одна из них) будут выражаться не три­ гонометрическими, а гиперболическими функциями, и тогда крае­ вые условия будут удовлетворяться только при С = 0, т. е. при отсутствии вихрей. Поэтому указанные комбинации интереса не представляют.

Кинетическая энергия поперечной циркуляции на единицу длины потока составляет

А ь

.,2

Р ^ rfy dz =—cj— bit,

О О

где v0 — среднеквадратическая

поперечная

скорость. Подставляя

в последнее выражение формулы

(11.22) и (11.23) и выполняя ин­

тегрирование, находим

 

C2bh

 

а _

 

 

Р

2/г2

 

и, следовательно,

 

 

 

где е = Э/рЫг — удельная энергия

поперечной

циркуляции, т. е. ки­

нетическая энергия поперечной циркуляции, отнесенная к единице массы жидкости.

Уравнения линий тока

(2.2)

в данном случае будут

 

 

 

 

dx

 

dy

 

 

 

 

 

 

dz

 

 

 

- J J

 

и

ІСш

тку

nnz

 

 

2Ст%

/тису

,

пкг

\

• '

 

г ^ т — s i n

т — C O S — г —

 

 

— Г о ? — c

o s

 

г —

 

 

 

 

k2h

b

h

 

 

k2b

 

b

 

h

 

 

 

 

Последнее равенство можно записать в виде

 

 

 

 

 

 

т

тку .

гшг

,

,

п

 

тъу

nnz

,

 

 

или

так:

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, уравнение проекций линий тока

на

плоскость

yOz есть

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

i n ^ s i n ^

£ .

= A .

 

 

 

 

(П.27)

 

Здесь А — постоянная,

которая

может

иметь

значения

от

—1

до 1, при этом, если m и п положительны, то

 

 

 

 

 

 

А = 1 п р и у = ( 2 / + 4 ) А

 

г = ( 2 ^ + ^ ) 4 ,

 

 

 

 

Л = - 1 п р и у = ( 2 / + 4 ) ^ ,

г = ( 2 А + 4 - ) 4 ,

 

 

 

Л = - 1 П р и у - ( 2 / + 4 - ) А ,

z = ( 2 A + 4 - ) 4 ,

 

 

и при любых целых тип

 

 

 

 

 

 

 

 

 

 

 

 

 

А—0

при у—i

— ,

z=k

— .

 

 

 

 

 

В этих формулах i n k — произвольные положительные целые числа. В потоке имеется \тп\ винтов; на рис. 12—14 это показано соответственно для т=1, п=\\ т=2, я = 1; т = 3, п = 2. При изме­ нении знака т или п (но не т и п одновременно) направления

вращения этих винтов и направления вихревых линий, показанных прямыми стрелками, изменяются на противоположные.

Чтобы получить уравнение проекций линий тока на плоскость

хОу,

обратимся

к первому

равенству

(11.26), которое представим

в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ішу

 

 

nizz

 

 

dy

g

a№h

 

 

 

 

 

 

Sin—JT-COS

Л

=

 

в.dx

'

2ъпс '

 

 

 

(11.28)

Возводя

в квадрат

(11.27)

и

(11.28) и

складывая,

получаем

уравнение, не содержащее z,

 

 

 

 

 

 

 

 

 

 

 

sin-- тку

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos-

тку

] Л

—A* sin с?

 

 

 

 

 

 

 

 

 

в это уравнение,

получаем

 

 

 

 

Рис.

13.

Поперечная

цирку­

 

dx2--

 

£>2 rf?2

 

 

 

 

ляция в

прямоугольном

ка­

 

1 —

(1 —

.42)

Sin2 =р

 

 

 

нале

(два винта,

т—2,

п=

 

 

 

 

 

 

 

 

 

 

 

 

=

!)•

 

 

 

 

 

D

иЫгЧ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2т&\тп\С

-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-x0=D

 

Йф2

 

 

 

 

 

 

 

 

 

 

 

 

 

VI—К2

Sin2 <}/

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

/ С 2 = 1 — Л 2 ;

F — символ

эллипти­

 

 

 

 

 

 

 

ческого интеграла первого рода. Урав­

 

 

 

 

 

 

 

нение (1.28) равносильно

следующему

Рис.

14.

Поперечная

циркуля­

 

телу -Ksn(K,

 

 

 

 

 

 

ция

в

прямоугольном

канале

cos

X ~ D

X a

) ,

 

(П.29)

(шесть

винтов, т=3,

 

п=2).

где

sn символ

так

называемого

эллиптического

синуса

 

Якоби.

Это

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

сложную периодическую

по

х

кривую.

В заключение необходимо подчеркнуть, что мы рассмотрели только возможные течения с поперечной циркуляцией в прямо­ угольном канале. На вопрос о том, какие значения т и п будут иметь место в действительности, будет ли действительное циркуля­ ционное течение представляться формулами (11.22) и (11.23) или двойными рядами из членов вида (11.22) и (11.23), в которых сум­ мирование будет производиться по т и п и т. д., изложенная тео­ рия ответа не дает. Выяснение всех этих обстоятельств требует анализа конкретных причин возникновения циркуляции.

Глава IV

Г И Д Р А В Л И Ч Е С К А Я И Д Е А Л И З А Ц И Я Т У Р Б У Л Е Н Т Н О Г О Т Е Ч Е Н И Я

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

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

ществуют

и имеют первостепенное

значение для

ряда

областей

науки, в первую очередь

для гидрологии и гидротехники,

задачи,

в которых

пренебрегать

эффектом

сопротивлений

принципиально

невозможно, т. е. нужно рассматривать не идеальную, а реальную жидкость, описываемую уравнениями Навье—Стокса.

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

малых

скоростях потока,

точнее

при

малых

числах

Рейнольдса

Re = VL/v, где

V-—характерная

(например,

максимальная в сече­

нии)

скорость

потока,

L — характерный

линейный

поперечный

размер потока

(например, диаметр

трубы),

v — кинематический

коэффициент вязкости, фигурирующий в уравнениях Навье— Стокса.

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

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

Чтобы получить уравнения турбулентного потока, рассмотрим уравнение (5.14)

- I f - + ( v • grad) v = F - - f grad p + -j - div T.

(12.1)

Изучая несжимаемую жидкость, характеризующуюся тем, что div v = 0, можно написать

 

 

(v

• grad)v=( v • grad) v-f-v div v.

(12.2)

Но, с другой стороны, написав проекции

(12.2)

на координатные

оси Ox, Оу, Oz, можно убедиться, что

 

 

 

 

 

 

 

(v

• grad) v-f-v div v = d i v ( v v ) ,

 

(12.3)

где тензор

(vv)

есть диадное произведение

вектора

v на самого

себя. Теперь уравнение

 

(12.1)

можно записать так:

 

 

d v

' - d i v ( v v ) = F — - g r a d / 7 +

— div Т.

(12.4)

Рассматривая турбулентное течение, положим

 

 

 

 

 

 

 

Р=Р+Р',

v = v

+ v \

 

 

(12.5)

где р и v математические

ожидания

давления

и скорости; р' и

v' флуктуации

этих

функций. Математическое

ожидание можно

мыслить как результат усреднения

р и v по реализациям процесса.

Пусть

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рі=Рі(х,

У> z .

0.

 

 

 

 

 

 

 

 

 

 

Vi = V ; ( A ' , у , Z,

t)

 

 

 

 

— давление

и скорость

 

в момент t

в точке

(х, у,

z) при і-том вос­

произведении исследуемого

движения

(при

і-той

реализации).

В силу случайности

процесса

ріфр,, Vi=^Vj,

если

£ = / . По опреде­

лению

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п

 

 

 

п

 

 

 

Р=

l

i m

 

2

А>

v , =

lim 4"

2

v,.,

(12.6)

 

 

 

л->• оо

1 — \

 

я - » - с с

; = = 1

 

причем предел понимается, конечно, не в обычном, а в

вероятностном смысле. Очевидно, что

p'=v'

= {vv')

=

0,

( v ' v ' ) ^ O ,

dv

dv

 

- j

,

~~дТ~~д~Г

grad /7 = grad/>,

div /==div T,

 

divv =

divv,

dv

dv

dv

 

dv

 

dv

dv

dx

dx '

dy

 

dy

'

dz

dz

Подставив p и v яз

(12.5) в (12.4) и взяв

операцию

матема­

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

 

d t

, div(vv) =

F-grad/? + ^ d i v | 7 , - P

( v ' v ' ) ) .

(12.7)

Строго говоря, намследовало бы также учесть флуктуацию

плотности,

положив р = р + р', чего мы, однако,

делать не

будем,

так как в дальнейшем будет рассматриваться почти исключительно

однородная жидкость, для которой p = const. Учитывая,

что divv =

= 0,

и используя (12.3) и (12.2), можно записать (12.7)

в

виде

-

| f + (v • g r a d ) v = F - ^ g r a d ^ + - ^ d i v { r - P ( v 7 v T ) } .

(12.8)

Уравнение (12.7) носит название уравнения усредненного дви­ жения. Компоненты тензора вязких напряжений Т определяются формулами (6.13). Вследствие линейности этих формул компоненты

усредненного тензора Т выражаются через компоненты и, v, w усредненной скорости так же, как компоненты актуального тен­ зора Т через компоненты и, v, w актуальной скорости.

Актуальное турбулентное движение не может быть установив­ шимся. Когда говорят об установившемся турбулентном движе­ нии, имеют в виду движение усредненное, а именно движение, ха­ рактеризующееся тем, что р и v не зависят от t. В этом случае флуктуации давления р' и скорости v ' есть стационарные случай­ ные процессы, в силу эргодичности которых определения р И V по формулам (12.6) совпадают со следующими определениями:

 

т

 

 

т

 

р= 11т

[pdt,

v = lim -

L

f v dt.

(12.9)

T - V CO

J

T CD

~ I

J

 

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

Диада (v'v') в уравнении (12.8), являющаяся симметричным тензором

(v'v')=

(12.10)

называется тензором пульсации. Тензор —p(v'v') часто называют тензором турбулентных напряжений. Этот тензор добавляется к тензору вязких напряжений и по отношению к усредненному дви­ жению действительно играет роль тензора дополнительных на­ пряжений. Обозначим через Т% тензор

 

 

7-, = 7 - - p ( v V )

 

(12.11)

и запишем уравнение

(12.7) в виде

 

 

 

 

 

^ ' . dlv(vv) = F — 1р

- grad/>+' - £1 - dlv

7V

(12.12)

 

dt

 

 

 

 

 

Уравнение

(12.12)

формально совпадает с (12.4) и может

трак­

товаться не

только

как уравнение

усредненного

движения,

если

под р и v понимать усредненные давления р и скорость v, а под тензором Т* — тензор (12.11), но и как уравнение актуального дви­

жения, если под р, v,

понимать

актуальные

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

и тензор вязких напряжений.

 

 

Векторное уравнение

(12.12)

эквивалентно

трем скалярным.

В ортогональной прямолинейной системе координат эти три урав­

нения и уравнение неразрывности

divv записываются так:

да

да*

 

dav

 

daw

-X-±

dp

dt

dx

 

dy

 

dz

dx

 

 

p

 

 

 

dz-

 

dzr

, d-r

(12.13)

 

 

 

dx

 

dy

dz

 

 

 

 

 

dv

dva

 

dv*

 

dvw

 

 

dt

dx

 

dy

 

dz

 

dy

 

 

 

dz yx

 

d' yy

dz J

(12.14)

 

 

 

dx

 

dy

 

dw

dwu

,

dwu

 

dw*

=z

dp

dt

dx

 

dy

 

д.

dz

 

 

 

 

J_

 

 

dzzy

dz,

(12.15)

 

+ P

I

dx

 

dy

dz

 

 

da

.

dv

dw

0,

(12.16)

 

 

dx

1

dy

1 dz

где Xxx, txy, и т. д. — компоненты тензора Т%

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