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

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

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

При выводе этого уравнения предполагалось, что t = const, а прямолинейная ось х, вообще говоря, не совпадает с дном. Од­ нако, если линия дна — плавная кривая, радиусы кривизны кото­ рой много больше глубины потока, то, как видно, по смыслу вы­

вода формулы

(24.1), можно принять за ось х линию дна

(т. е. кри­

вую), но тогда, конечно, уклон і будет функцией

х. Уравнение

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

(13.24) остается без

изменений:

 

 

 

 

 

 

 

 

 

Т

 

+ Т = ° -

 

 

 

 

 

<2 4 -2 >

Имея U = Q/F, напишем:

 

 

 

 

 

 

 

 

 

dt

 

F '

dt

F2

'

dt

'

 

dx

F

'

dx

F2

\

dx

dF

_ R

dH

 

I

dF

\

_

dF

dH

,

dF

R

dH

.

dF

dt

~~

dt

'

\

dx

 

 

dH

dx

 

dx

 

дх

^

dx '

Здесь (dF/дх)^— полная производная F по x, учитывающая из­ менение F как от непризматичности русла, так и от изменения Н с изменением х, а dF/dx есть та часть производной F по х, которая связана только с непризматичностью русла. С помощью этих соот­ ношений уравнения (24.1) и (24.2) приводятся к виду:

\

gF$ j дх

\ К2

 

gFZ

dx ) ^ 1

gF2

dx 1

 

 

 

1

dQ

 

 

(24.3)

 

 

 

gF

dt

 

 

 

 

 

 

 

 

 

5

T

+ - ^ = ° -

 

( 2 4 ' 4 >

Из (24.3) видно, что уравнение невозмущенного установивше­ гося движения есть

Рассматривая малые возмущения, с малыми производными, по­ ложим # = # а + Д # , Q = Q0 +AQ и отбросим в (24.3) и (24.4) выс­ шие степени и произведения АН, AQ и их производных. Используя далее уравнение (24.5), получим следующие уравнения возмущен­ ного движения:

 

 

 

dH \

AQ

 

. 0о

CLQ

 

d АО

.

1

 

dbQ

і

 

 

 

 

 

Q

 

 

 

gF2

'

dx

^

gF

 

dt

 

 

, ( - , _

BQ2

\

dhH

j

n

o

\ 2

 

dK

t_lr>_B_

 

dF

dB \]

>[l

a gF*

)

dx

I

у

 

[ K3

"

dH

 

^ 3

I

F

'

dx

dxjV

 

 

-

^ .

( З

 

4

_

^ . ^

) ^

)

д я

_

0 ,

 

(24.6)

 

 

 

gF3

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Из уравнения

(24.7)

видно,

что существует

такая

функция

Ф(*. О, что

 

 

д Ф

 

!

 

д Ф

 

 

 

 

 

Подставляя AQ и Д Я из (24.8)

в уравнение

(24.6),

получаем

уравнение

 

 

 

 

 

 

 

 

 

 

 

 

 

-ЙГ+РЛх)

_

+ 2 а £ /

- ш г + Л

( л ) - 5 Г + Л ( ^ ) - 5 ^ - = 0 ,

 

(24.9)

в котором

 

 

 

 

 

 

 

 

 

 

 

 

 

 

^ W = - f

 

( ' - ^ ) . р 3 )=\^ 2 ^-1),

 

 

 

 

 

 

 

2BFy.

а

/

3

 

d)F

1

 

дВ

 

 

 

 

1(2

W2 \F

 

дх

В

 

дх

 

г / о в

і

дв \ дн . , ^

 

1 Ч

і

a s )

 

 

-

а ^ ( 3

 

в" • " Ж " ) Т + ( а Х

-

 

-дх~

\ •

 

 

a X и (.і определяются

по формулам

(21.8),

но теперь это, конечно,

будут не постоянные величины, а функции

х.

 

 

 

 

 

Пусть при х = 0 имеет место гармоническое

возмущение с угло­

вой частотой со, тогда решением уравнения

(24.9)

будет

 

 

 

 

 

Ф(х,

t)=<?(x)cosut

+

 

ty(x)slnu>t.

 

 

(24.10)

Из формул

(24.10) и (24.9) следует

 

 

 

 

 

 

 

 

 

-

ш2<? - j - со/7 ,ф+2u>a Щ' +р2<?' +/? 3 с?"=0,

 

 

 

 

- ш2ф -

со/7, <р + 2соа с/ср' +/? 2 ф' + / > 3 ф * = 0 ,

 

 

(^4.11)

где штрихи обозначают дифференцирование по х. Из (24.10)

видно,

что квадрат амплитуды

волны в сечении с абсциссой х

есть

ф2 +ір2

и что, следовательно, волна затухает при перемещении

вниз по

руслу, если

 

 

 

срср' +

ф ф ' < 0 .

 

 

 

 

 

 

(24.12)

 

 

 

 

 

 

 

 

 

 

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

нарастает,

т. е. когда

неравенство

(24.12)

 

заменяется

равенством

 

 

 

 

срср' +

ф ф ' = 0 .

 

 

 

 

 

 

(24.13)

Выясним, каким условиям должны удовлетворять функции ри р2, рз, U для того, чтобы имело место равенство (24.13). Из этого равенства следует

" V = P W .

(24.14)

где р(х)—некоторая неизвестная функция. Из (24.14)

получаем

<Р' + Ф' = - Р ( ? - « ! > ) , < Р ' - Ф ' = Р ( « Р + Ф ) .

(24.15)

Складывая уравнения (24.11), а затем вычитая второе уравне­ ние из первого и используя соотношения (24.15), получаем

(со2 +2шаUp +р3р2) (<р + ф) +

+/?2 р +/>зР') (<Р Ф ) = О ,

(ш/?, +/?2 р +/?3 р') (<р + Ф) — (<°2 + 2

ш « + АР2 ) ( ? — Ф ) = 0 .

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

сВ2 + 2свас/р+/73 р2=0, « У і + Л Р + Л Р = 0 .

Из первого уравнения (24.16)

а£/ + / а 2 £ / 2 _ р

/>з

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

Pi (*U2-p3)

\р>Рг-2а

[(р2 -р'г) и+Ръи] \-<гр\и'г+

 

 

 

'2

 

+Р2{р2-р'г){*и2-Рз)

+ *РзРъии'--^=0.

(24.17)

Это соотношение не зависит от частоты со и, следовательно, оно дает границу 5-устойчивости и неустойчивости по отношению к лю­ бым малым возмущениям (см. § 19). Нас, однако, интересует вы­ полнение условия (24.12), а не (24.13). Очевидно, что для того, чтобы получить критерий выполнения условия (24.12) в формуле (24.17) знак равенства следует заменить знаком больше ( > ) или знаком меньше ( < ) . Чтобы сделать выбор между этими двумя знаками, достаточно установить знак неравенства для одного ка­ кого-либо случая. Таким частным случаем может служить призма­ тическое русло, для которого левая часть (24.17) с точностью до положительного множителя равна

I*2 — 2 а ц + а — р - .

Сопоставляя этот результат с критерием (21.15), легко прийти к выводу, что критерий устойчивости для непризматического по­

тока получается

из

(24.17)

заменой

знака

равенства

знаком

меньше ( < ) :

 

 

 

 

 

 

А 2и2г)

1р3-

[{р2 -р'з)

U +

p3U')}~aplU'2+

 

 

 

 

 

'2

 

+Р2 (Р2-Рг)

(*U2-Рз)

+ ьРзРъии' -

< 0.

(24.18)

Следует заметить, что к непризматическим одномерным пото­ кам относятся не только потоки в непризматических руслах, но и неравномерные потоки в призматических руслах. В противополож­ ность призматическим потокам условие (24.18) для непризматиче-

172

О 20 60 100 150 200

300

WO Xfif

Рис. 28. Зоны устойчивости и неустойчивости течения в прямоугольном сужающемся лотке переменного уклона.

1 = 0,207 — 0.00128 .ї+0,000002559 Х-; 6 = 3 — 0,004 X; т=0.

О 20 60 100 150200

300

-400 хм

Рис. 29. Зоны устойчивости и неустойчивости в трапецеидальном сужа­ ющемся лотке переменного уклона.

/=0,207 — 0,00123 А"+0,000002559 Х'\ 6=3 — 0,004 х; т = 1.

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

Для иллюстрации на рис. 27—29 приведены результаты расче­ тов устойчивости по критерию (24.18) для трапецеидальных русел с переменной шириной по дну Ъ, переменным п постоянным укло­ ном і и постоянным заложением откосов т. Нанесенные на этих рисунках кривые есть кривые свободной поверхности при невозму­ щенном движении с разными расходами Q м3/с, рассчитанные чис­ ленным интегрированием уравнения (24.5). На каждой из этих кривых штриховкой отмечены участки, на которых не выполняется условие устойчивости (24.18). Коэффициент шероховатости п для всех трех русел составляет 0,02, корректив количества движения а принимался равным 1,08.

25.Самопроизвольное возникновение бегущих волн

вдвумерных потоках

Неустановившееся движение двумерного потока выражается уравнением (15.3), уравнением получающимся из (15.3) заменой индекса 1 индексом 2, а индекса 2 индексом 1, и уравнением (15.14). Учитывая результаты § 21, мы и здесь будем пренебрегать инерционными силами, связанными с компонентами скорости по координатному направлению хз, и запишем перечисленные урав­ нения так:

 

 

 

 

 

1

 

дН

\\

,

„ „ . . .

 

Щ

 

 

 

 

 

 

 

 

 

 

 

 

'

Н

 

дГ~

 

 

а

dL,UHU2 Ї

,

аa

 

OLdL%HU,U2tiU,U2

.

а

 

д ц

ни

1

L % _

 

дх,

1

L,L

2

 

дх

2

 

 

/.J

дх,

 

-*н[ч+чя

 

- £ ) ]

+

^

:

+

^

+ ~

- х

 

dL2L,HUl

"т"

а _

'

 

dI*HU,Ut

 

 

 

а_

'

д ц

HUz-

 

Л

дх2

L2L,

 

дх,

 

 

 

L*

дх2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L2L,

 

_oL, NUf=0i

 

 

 

 

 

(25.2)

 

 

 

 

 

дх2

 

 

 

 

 

 

 

 

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

будет Ui = U0, L'2 = 0, Я = Яо- В результате возмущения

будем иметь

[/і = £/о+Л£/,

U2=AU2, # = Я о + Л # . Однако

анализ

устойчивости

по отношению

к произвольному возмущению

вектора

U очень сло­

жен и громоздок, и мы ограничимся только такими возмущениями,

которые меняют только длину вектора

U, но не

его направление,

т. е. положим Д£/2 = 0. Тогда, согласно

формулам

(15.23),

Ji=J-=-D+T'!T-

 

( 2 5 - 4 )

Уравнение (25.2) для возмущенного движения вырождается при этом в уравнение, связывающее At/ и дАН/дх2, которое не пред­ ставляет интереса, уравнения же (25.1) и (25.3) приобретают вид:

г , / ,

.

1

дН \ і , . . 1

дни .

 

д-ЩШ1

-J-.J^LHIP-O,

(25.5)

дН

.

1

dL2HU

 

 

dt

1

ЦЦ

дху

= 0 .

(25.6)

 

 

 

 

Входящий в уравнение (25.5) уклон сопротивлений / зависит от величины w, определяемой по формулам (15.17) и (15.22) — (15.26). Компонент U3 в этих формулах учитывает влияние искривления струи в плоскостях главной кривизны координатной поверхности (хи х2) на гидравлические сопротивления. Нужно подчеркнуть, что учет этого влияния не есть учет влияния деформаций линий тока локальных усредненных скоростей, о котором говорилось в § 21: здесь имеется в виду только кривизна линий тока вектора U сред­ ней скорости по нормали к координатной поверхности. Поэтому введение в формулы для сопротивлений компонента U3 не может исправить некорректность, связанную с попыткой учесть силы инер­ ции, связанные с компонентом &з локальной скорости, которые здесь не учтены. Мы будем также кривизну дна потока считать пре­ небрежимо малой. Тогда

 

 

дЦ

(25.7)

 

 

дх2

 

 

 

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

 

 

 

, 1

1 ,

} U2.

(25.8)

•\

D ^

gLxL2 дх2

 

Введем

новую

переменную

q = LiHU

 

и

напишем

уравнения

(25.5) и (25.6) в следующем виде:

 

 

 

 

 

 

 

 

 

 

 

 

 

дН

 

 

 

К2

1

 

1

dq

 

 

 

 

 

 

 

 

 

 

 

gL2H

dt

 

 

 

 

 

 

 

 

 

 

 

а

 

 

dLi

q*

(25.9)

 

gHLJL2

'

<>хх \Li

Н

 

 

gL-tLr2

 

 

дХх

 

 

 

 

 

 

 

 

 

 

 

 

 

дН

_j_

1

 

J2_

 

 

 

 

 

(25.10)

 

 

 

 

 

dt

і

цц

 

 

 

 

 

 

 

 

 

 

dx,

 

 

 

 

 

 

Здесь К определяется

формулой

 

 

 

 

 

 

 

 

 

 

 

1

 

1

1

1

,

 

 

dL,

 

 

(25.11)

 

 

 

К2

 

 

 

 

 

 

 

dxc,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и является аналогом пропускной способности. Уравнение

(25.9) мо­

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

, .

1

/

 

.

a

 

q*

\

дН

 

 

 

 

 

 

 

 

 

 

 

dL2

 

1

 

 

 

 

 

 

 

 

К2

 

gL.Ll

 

д*і

 

 

 

 

 

 

 

 

 

2з<7

 

dq

.

1

 

 

 

 

 

(25.12)

 

 

 

 

 

 

 

dx.

 

 

 

 

dt '

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как всегда єз<0, то легко заметить,

что это уравнение —

аналогия уравнения

(24.3).

 

 

 

 

 

 

 

 

 

 

Запишем

уравнения

(25.12)

и

(25.10)

 

в

приращениях:

2(ч +

чГі+і

 

дН

\

Д?

 

 

2aq

 

d&q

 

dhq

 

дхх

I

q

 

gH-LxL'2

 

dx.

 

 

dt

 

 

 

 

 

 

 

 

о /

dK

 

 

За

 

дЦ

\

 

3

 

 

dH

Д Я

 

/<3

dH

 

gH'L.Ll

 

dx.

 

 

 

 

dx

H

 

 

 

 

 

 

gLi

 

Hi

dhHdxi

=

0,

 

(25.13)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а д я

 

 

 

dbq

=

0.

 

 

(25.14)

 

 

 

 

 

dt

 

L,L

 

dx.

 

 

 

 

 

 

 

 

 

1^-2

 

 

 

 

 

 

 

Здесь H и q относятся к невозмущенному движению, но ин­

декс нуль для простоты опущен.

 

 

 

 

 

 

 

 

 

Дальнейшее

решение

выполняется

в

точности так же, как

в § 24, оно приводит к условию

устойчивости

(24.18), но и с дру-

гимн значениями функции рь рг, рз, U, а именно:

 

 

2-?Я/.2

є /

1

4

<ЭЯ

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

д/<

За

 

(3Z-2

За?2

дН

 

 

<?Я

g H 3 ^

'

дхі

gLxL\Hb

 

ГЯ

1

дЦ

, 1

д і я

-з" ^ 4 я 3

 

 

 

 

 

дхі

 

 

Критерий (24.18) в данном случае указывает на возрастание или затухание возмущений при различных значениях хг, т. е. вдоль определенных этим значением линий тока. При этом устойчивость может иметь место на одних линиях тока и не иметь на других.

Проверка устойчивости по

критерию (24.18) в данном случае,

так же как и в предыдущем,

не представляет затруднений. Глав­

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

вглаве IX.

26.Вторичные течения в турбулентных потоках

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

тивлений,

которые моделью иде­

Рис. 30.

Схема

циркуляционных

по­

альной жидкости не учитываются.

токов в

трубах

треугольного и

пря­

Вторичные

же

циркуляционные

 

моугольного сечений.

 

потоки, которые

имеются в виду в данном случае, возникают без ка­

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

ким же, как, например, турбулентность. Эти потоки в

треугольной

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

трубах представлены

на рис. 30. Их

существова­

нием объясняют

[48, 57] специфическую

форму изотах

продольных

компонентов скорости: в трубах выпуклого сечения изотахи не об­ разуют выпуклых фигур, а имеют в разных точках кривизну раз­ ных знаков, как это видно на рис. 31. Аналогичная особенность изотах наблюдается и в открытых руслах (рис. 32). Следовательно,

12 З а к а з № 428

177

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

Рис. 31. Изотахи в трубе прямоугольного сечения.

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

б)

 

 

 

с"

 

 

 

 

 

 

1

 

С

 

 

 

 

 

 

 

Рис. 33. Изотахи и линии градиента

 

 

модуля

скорости (а)

и график коэф­

 

 

фициента шероховатости

(б)

потока

 

 

в бесконечно широком русле с раз­

Рис. 32. Изотахи

в открытом

личной

шероховатостью

на

смочен­

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

канале.

 

ном периметре.

 

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

мер по формуле Маннинга С=—где

l = h, коэффициент шеро-

ховатости слева от точки А будет ППІ, а справа от точки А — п = = №• Однако распространить этот расчет на окрестности точки А нельзя. Даже если учесть искривление изотах и линий градиента вблизи этой точки (рис. 19), предполагая, что поток призматиче­ ский, мы все равно придем к противоречию. Обобщенный гидрав­ лический радиус / вблизи точки А будет отличен от h, но при дви­ жении поперек потока он будет изменяться плавно, в то время как коэффициент шероховатости при переходе через точку А сделает скачок (рис. 33 б). В связи с этим при таком подсчете получится, что средняя скорость U.i: в отсеке между двумя бесконечно близ­ кими линиями градиента по формуле (16.10) также сделает скачок

при переходе через точку А, а такого скачка не может

быть. В дей­

ствительности

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

и справа от

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

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

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

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

Таким образом, поперечная циркуляция выравнивает касатель­ ные напряжения на смоченной поверхности. В связи с этим необхо­ димо заметить, что опыты с трубами треугольного и прямоуголь­ ного сечений, а также с круглыми трубами, имеющими продольный гребень на внутренней поверхности, показывают [57], что зависи­ мость коэффициента сопротивления от числа Рейнольдса, рассчи­ танного по гидравлическому радиусу, для этих труб в точности со­ впадает с соответствующей зависимостью для обычных круглых труб, в которых распределение касательных напряжений по смо­ ченной поверхности заведомо равномерное. Отсюда напрашивается вывод: те средние напряжения, которые рассматриваются при обоснованиии формул для гидравлических сопротивлений, построенных на понятии гидравлического радиуса, в действительности близки к фактическим (выравненным за счет поперечной циркуляции), так как иначе можно было бы ожидать различия в законах сопротив­ ления для труб с различными формами сечения. Согласно § 16, ис­ кривление линий тока за счет поперечной циркуляции должно

12*

179

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