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

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

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

вектор v = gradcp. Функция ср называется

потенциалом

скоростей.

Для потенциального поля уравнение (5.12)

приобретает

вид

-^-grad<p +

grad(-y-) = F — ^grad/> .

(8.5)

Будем опять иметь

в виду баротропную жидкость.

Положим

р (р)

тогда как нетрудно проверить,

1р grad/?=gradP

и уравнение (8.5) переписывается так:

F = g r a d ( - | j - + - £ + p ) .

(8.6)

Отсюда видно, что F есть потенциальный вектор, иными сло­ вами, безвихревое движение возможно только тогда, когда объем­ ные силы имеют потенциал, т. е. F = —grad V. Теперь (8.6) приво­ дится к виду

grad ( - £ - + - £ + ^ + / > ) = 0 ,

(8.7)

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

 

dt

V+P=W(f).

(8.8)

 

 

 

 

 

Здесь

—произвольная

функция. Выражение

(8.8) назы­

вается интегралом Коши.

 

 

 

 

Если жидкость несжимаема, то уравнение неразрывности при­

обретает вид

 

 

 

 

 

 

дх°- 1 ду2

1

дг*

-О,

(8.9)

 

 

 

и решение задачи о безвихревом движении сводится к отысканию одной функции ф, удовлетворяющей уравнению Лапласа (8.9), краевым и начальным условиям. Гидродинамическое давление опре­ делится из (8.8), а произвольная функция х¥ (t) может быть най­ дена, если задана зависимость р от t в одной точке поля. Если, кроме того, движение установившееся, а объемные силы сводятся к одной силе тяжести, то интеграл Коши принимает вид

z + ^ + - | - = c o n s t ,

(8.10)

внешне совпадающий с интегралом Бернулли (8.4). Однако между выражениями (8.4) и (8.10) есть существенная разница. Первое из них не предполагает, что движение безвихревое, и в общем слу-

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

вихревым.

 

 

 

Приведем

теперь без

доказательства

некоторые теоремы

•о вихрях.

 

 

 

Т е о р е м а

Л а г р а н ж а.

Если объемные

силы имеют потен­

циал и если в какой-то момент в некоторой части идеальной баро-

тропной жидкости нет

вихрей (rotv = 0), то их не было

раньше и

не будет позже в этой

части жидкости (подчеркнем, что

речь идет

об определенной массе жидкости, а не об определенной части пространства).

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

дующая

т е о р е м а

Г е л ь м г о л ь ц а .

Если

объемные

силы имеют потенциал, а жидкость идеаль­

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

дующие

моменты,

а

интенсивность любой вихревой

трубки

во все время движения

остается постоянной.

 

Эту теорему называют также теоремой о сохранении

вихревых

линий и

вихревых

трубок. Более общая теорема, не связанная

с идеальностью и баротропностью жидкости и потенциальностью объемных сил, принадлежит А. А. Фридману: необходимое и доста­ точное условие сохранения векторных линий и интенсивности век­

торных трубок вектора Е

(в нашем

случае £ = rotu) состоит в том,

что г е л ь м г о л ь ц и а н

этого

вектора

h e l m E = - ^

- ( E ,

V)v + E d i v v

должен быть равен нулю.

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

Возникновение вихрей связано с разными причинами, основные из них следующие:

1) небаротропность жидкости. Если в начальный момент изоба­ рические поверхности (поверхности равного давления) не совпа­ дают с изостерическими (поверхностями равной плотности), то в жидкости возникают вихри, тем более интенсивные, чем выше гра­ диенты плотности и давления и чем больше угол между векторами градиентов. Это главная причина образования вихрей в атмосфере и вызываемых ими циркуляционных течений;

4*

51

2) иепотеї-іциальпость действующих на жидкость сил, в первую очередь сил внутреннего трения (вязкости). Мы увидим в следую­ щих главах, что движение жидкости с внутренним трением вообще ие может быть потенциальным. Непотенциалыюсть вектора объем­ ных сил F, имеющая место в том случае, когда, кроме силы тяже­ сти, приходится учитывать еще силу Кориолиса, связанную с вра­ щением Земли, порождает вихри, даже если внутренним трением

вжидкости можно пренебречь;

3)индуцирование вихрей движущимися в жидкости телами, даже когда жидкость идеальная. Третья причина связана с внеш­ ней задачей гидромеханики, не рассматриваемой в этой книге.

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

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

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

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

Глава III

Н Е К О Т О Р Ы Е Д В И Ж Е Н И Я И Д Е А Л Ь Н О Й Ж И Д К О С Т И

9. Простейшие волновые движения жидкости постоянной плотности

Чтобы выяснить одну из важнейших причин возникновения волн, рассмотрим действие мгновенных сил на жидкость. Для этого на­ пишем уравнение Эйлера (5.12) в таком виде

 

- 3 T = F — \ - & * * Р -

 

(9-І)

 

Положим

 

 

 

 

 

 

 

F = F 0 +

F»,

р=Ро+Р*,

 

 

 

где

и р%— составляющие

вектора

F и скаляра

р, отличные

от

нуля только в течение некоторого интервала времени от

^ = 0

до

t=x,

но принимающие в течение

этого интервала

очень

большие

значения, так что пределы

 

 

 

 

 

 

 

lim \F%dt=l,

 

lim

\p^dt—v:

 

 

 

конечны. Эти пределы называются соответственно импульсом мас­ совых сил и импульсом сил давления. Интегрируя уравнение (9.1) по t от ^ = 0 до t = x, переходя к пределу при т—>-0 и замечая, что для конечных Fo и ро соответствующие пределы равны нулю, полу­ чаем:

v - v 0 = I - g r a d ( - ^ ) .

(9.2>

Здесь Voскорость в момент £ = 0 до начала действия

мгновен­

ных сил.

 

Наибольший интерес представляет случай 1=0, когда мгновен­

ные силы есть только силы давления. Тогда

 

v - v 0 = - g r a d ( - ^ - ) .

(9.3>

Из этого выражения можно сделать следующие выводы.

1. Если до приложения импульса давлений движение

жидкости

было безвихревым, т. е. было vo = grad<p, то после приложения им­

пульса оно останется

безвихревым с потенциалом q> — л/р + Су

где С — произвольная

постоянная.

2. Безвихревое движение может возникнуть из состояния покоя (УО = 0) в результате приложения импульса давления. Но никаким импульсом давления нельзя вызвать вихревого движения покоя­ щейся жидкости.

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

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

дх* ^ ду2 Т" дгч — и -

Если р = const, то из уравнения (8.8) получаем

Будем предполагать, что движение совершается в поле силы тя­ жести без участия других массовых сил. Направим ось z верти­ кально вверх, тогда V=gz. Если движение достаточно медленное, то членом v2[2 в выражении (9.5) можно пренебречь и записать его так:

Р _

d?i

•gz,

(9.6)

р

dt

 

 

где

 

 

 

Так как уравнение (9.4)

сохраняет силу, если заменить

в нем <р

на фі, то индекс единица при ф в формуле (9.9)

можно отбросить

и написать ее в виде

 

f — Я - - * * -

<9-7>

Нормальная составляющая скорости на твердых поверхностях, ограничивающих жидкость, равна нулю, следовательно, на этих по­ верхностях имеет место следующее краевое условие:

0-8)

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

•&+gz=0.

(9.9)

Пусть уравнение свободной поверхности есть z — \(x, у , t). Крае­ вое условие (9.9) можно записать так:

г д<? (х, у, z,

t)

 

dt

- Ц ( , . , . о + # = 0 -

( 9 Л 0 >

Плоскость Оху совместим

со свободной поверхностью

жидкости

в положении равновесия. Мы ограничимся волнами, высота кото­

рых очень мала по сравнению с длиной,

тогда в

(9.10)

можна

брать значение dcp/dt

при 2 = 0, а не при z = £, т. е.

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

ді

 

 

1

(х,

у, 0,

t)

 

 

(9.11>

 

 

dt

 

 

g

 

д&

 

 

 

 

 

 

 

 

 

 

 

 

но так как

 

 

dx

 

 

rfy

 

dz

 

 

 

 

 

 

 

 

W —

 

 

 

 

 

и=—гг,

'

- V = —тг ,

—rr,

 

 

 

 

 

dt

 

dt '

 

dt

'

 

 

то на свободной

поверхности

 

 

 

 

 

 

dl

dS ,

dS

dx

, .55

dy

dt

.

d5

, „, й

/ Q 1 Г Л

®=1Г=-эт+-1Г-^гт-^г--^г=^г+и -йГ+г , ^г- (9Л2>

Так как высота волн мала по сравнению с их длиной, то плос­ кость, касательная к свободной поверхности, мало отличается or горизонтальной плоскости и, следовательно, дЦдх и дЦду в (9.12) пренебрежимо малы. Значит на свободной поверхности

dz

 

dt

 

и краевое условие (9.11) можно переписать так:

 

дер

1

д 2 ?

(9.13)

~~дТ~

g~ ' ~d~W

 

при z = 0.

 

 

 

Имея в виду волновые движения,

мы не будем формулировать

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

нения (9.4)

при краевых условиях

(9.8)

и (9.13). Для

этого по­

ложим

cp=cos (аг!-(-є) Ф (л:,

у,

z).

(9.14)

 

Подставляя значение <р в (9.4),

(9.8)

и

(9.13), видим, что функ­

ция Ф должна удовлетворять уравнению

 

 

 

 

^ " + ^ +

^ =

0

 

0М5>

и краевым

условиям

 

 

 

 

 

- Т и - *

 

 

(9.16)

на твердых границах и

(9-17)

ПрИ 2 = 0.

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

скорости

(9.14) есть

<p=COS(a* + E)i]>(jc, z).

 

(9.18)

 

 

 

Уравнение (9.15)

записывается так:

 

 

 

 

 

<?2ф

,

<Э2ф

О,

 

(9.19)

 

 

дх*

~

дг*

 

 

 

 

 

 

а краевые условия

(9.16) и (9.17) остаются

без

изменения.

Будем рассматривать водоем постоянной глубины h. Непосред­

ственной

подстановкой можно

проверить,

что

уравнение (9.19)

имеет частное решение вида

 

 

 

 

 

 

 

Ф(л-, z ) = P

(2) sin £ (* -- /)),

(9.20)

•если Р(г)

удовлетворяет дифференциальному уравнению

 

 

Р"

 

(z)-k2P(z)=0.

 

 

Из этого уравнения

 

 

 

 

 

 

 

P(z)=Clek2~\-C2e-'iZ.

 

 

 

Согласно краевому условию

(9.16), должно быть дФ/дг = 0 при

z — —h, т. е.

C,Ae* A - C 2 te - * f t =0 . Это условие удовлетворяется, если принять

т. е.

 

P{z)=Cohk{z+h).

(9.21)

Из условия же (9.17) получаем

 

e2=gk th kit

 

Следовательно, в данном случае

 

ср=С ch к (z+A)sin к {х - т)) cos (ot+e).

(9.22)

Интегрируя уравнение (9.11), найдем уравнение свободной по­ верхности

Здесь ар (л:)—произвольная функция. Но, согласно сделанному выбору системы координат, гр (х) = 0 и, следовательно,

е = - - у

t ] L o = а s i n / г - ^ s i n ^ + є > -

<9 -2 3 >

где

 

 

 

a = - ^ - c h A A

(9.24)

— амплитуда волны. Різ (9.23) видно, что в точках с абсциссами, удовлетворяющими условию кхц=пп (п — любое целое число или нуль), т. е.

хЪ '

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

Такие точки называются узлами. В

точках, в которых kx — т) =

= я я + я/2, амплитуды вертикальных

перемещений равны а, т. е.

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

стоячих волн. Длина

такой волны K = 2n/k,

а период Т = 2л/а.

Проекции скорости частиц на координатные оси равны:

a==~^-==-^-=Ck

ch k (z 4-/г) cos k (x — -q) cos^-J-s),

w—-^~—~^~

= Ck sh k (z-\-h)sink(x

y\) cos(at-\-z).

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

 

 

 

Ck

 

 

 

 

 

 

х — х0=—— ch k(zQ-\-h) cos k(x0

т}) sin (at-{-є),

 

 

 

 

 

Ck

 

 

 

 

 

Отсюда z — z0

 

sh

k (z0 - j - h) sink (x0

7j) sin (at-f- є).

 

 

 

 

z

z°- =

 

thk(z0+h)tgk(x0-rl),

 

 

 

 

 

••*o

 

 

 

 

 

т. е. траектории

частиц

есть прямые линии. При этом

в

узлах

sin£(xo — т)) = 0 и

частицы

совершают

колебательные

движения

точно в

горизонтальном

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

в пучностях

же,

где

cos£(xo — т)) =0 ,

частицы

колеблются вертикально. Амплитуда ко­

лебаний

есть

 

 

 

 

 

 

 

Ck ,

, ,ч

sh k (z0 + h)

В случае бесконечной глубины (/г-н»со) эта амплитуда будет aehz\ Отсюда видно, что волны с глубиной очень быстро затухают:

например,

на глубине, равной длине

волны, т.

е. при Zo =Х =

=—2n/k,

амплитуда будет а е _ 2 я , т. е.

примерно

в 500 раз меньше

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

Мы считали до сих пор, что водоем простирается в направлении оси х неограниченно. В этом случае длины волн могут быть лю­

быми. Если же он при

х = 0

и

x = L ограничен

вертикальными

стенками, нормальными

к оси х,

то, согласно условию

(9.16), дол­

жно быть дФ/дх — О при х = 0 и при x = L , или

 

 

COS£T) = 0,

COS k (L — v))=0.

 

 

Отсюда следует

 

 

 

 

 

=

\-iriK,

£(/,_ 7)) = -|-4-/даГ

 

 

и значит

 

 

 

 

 

где т — любое натуральное число или нуль.

 

 

Вернемся к случаю

водоема,

неограниченного

в

направлении

оси л-. Нетрудно видеть, что всем уравнениям и краевым условиям

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

(9.22), но и потенциал

 

c?=Cch k(z-\-h)cosk(x

— -q) sin (at-[-в).

(9.25)

В силу линейности уравнений и краевых условий разность потен­

циалов

(9.22) и

(9.25), т. е. функция

 

 

 

 

 

 

 

<P=Cchfc(z + A)sin(fce —о/)

 

 

(9.26)

(здесь

для простоты

принято

г| = е = 0 ) , также

есть частное

реше­

ние

рассматриваемой

задачи.

Уравнение свободной

поверхности

в этом случае будет, по аналогии с (9.23),

 

 

 

 

 

 

 

 

£ = a c o s ( / b c — at).

 

 

(9.27)

 

Постоянные

а

и

С по-прежнему связаны

формулой

(9.24).

Из

(9.27) видно,

что профиль

волны, имеющий

вид

косинусоиды,

перемещается со скоростью

 

 

 

 

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

1) если глубина жидкости очень велика по сравнению с длиной волны, то по соотношению

и m t h г - = 1

имеем

 

 

 

 

 

' - / - g

r ;

(9-29)

2) если глубина жидкости очень мала по сравнению с длиной

волны, то

 

 

 

 

,,

 

2кк

2тсЛ.

 

t h

~

>

и тогда

 

 

 

 

 

c--=Vgh,

 

(9.30)

т. е. скорость перемещения волн не зависит от их длины. Скорости и траектории жидких частиц находятся, как и в пре­

дыдущем случае, по формулам:

 

 

 

 

u=^=-^=Ck

ch k{z+h)

cos (kx -

at),

(9.31)

•w=-^-=-^

=

Ckshk{z^-h)sin{kx-at),

 

(9.32)

 

ch k(z0+h)

cos (kx0-at),

 

(9.33)

Ck

 

 

 

 

 

z — zQ =—^— sh k (z0-f-/z)sin (kx0 — at).

 

(9.34)

Исключая отсюда

время

и выражая

С через

амплитуду волн

а, получаем уравнение

эллипса

 

 

 

 

( ^ ) 2 + ( ^ Н '

 

(9-35>

где

 

 

 

 

 

 

 

А = а - ch k 0 +

К)

 

 

 

 

 

sh kh

 

 

 

 

 

 

sh k (z0 +

h)

 

 

 

 

 

sh kh

 

 

 

 

Таким образом, частицы двигаются

по эллипсам с полуосями А

и В. Для волн, распространяющихся в положительную

сторону

оси Ох, частицы вращаются

по эллипсам

по часовой стрелке. Как

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

Обратимся к более общему случаю прогрессивных

волн — вол­

нам на поверхности раздела двух жидкостей (рис. 8).

Плоскость

Оху совместим с поверхностью раздела жидкостей в состоянии равновесия. Пусть движение совершается между двумя параллель­ ными плоскостями z=hi и z = hi, причем жидкости текут парал­ лельно оси Ох со скоростями Ui и ІІ2- Очевидно, что плотность Р2 верхней жидкости должна быть не больше плотности pi нижней

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