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

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

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

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

<fi = UiX,

< р 2 = с / 2 х .

При волновом движении на эти потенциалы накладываются еще потенциалы прогрессивных волн. Для нижней жидкости, согласно (9.26), необходимо принять

4 l = UlX^rCl

ch k (z+h,)

sin (kx-at).

(9.36)

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

 

ср2== Ц2х -|- С, ch k (z - h2)

sin {kx - at),

(9.37)

чтобы удовлетворить краевому условию дф2ЛЗг = 0 при z = h%

Пусть уравнение поверхности раздела жидкостей

есть z = %{x, t).

г

 

 

 

1

- и2

t

X

^

Рис. 8. Волны на поверхности раздела двух жидкостей.

Скорость частицы нижней жидкости, расположенной у поверхно­ сти раздела, равна

dz

де.

dx

,

ді

dt

dx

dt

'

dt

Ho dx/dt = ui есть горизонтальная составляющая скорости час­ тицы, при малых волнах мало отличающаяся от Uit следовательно, с точностью до малых высших порядков

при z = £(x, t). Рассматривая малые волны, можно добиться, как

и в предыдущих случаях, выполнения

условия (9.38)

не при 2 = |,

а при 2 = 0. Учитывая это, полагая, как и выше

 

Ь{х, t)=acos(kx

— at),

(9.39)

и используя (9.36), получаем из (9.38)

 

 

C,JfeshAA, = a(o Ш , ) .

(9.40)

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

 

— C.2kshkh2=a(e-kU2).

(9.41)

Для

жидкости постоянной плотности I^P=-^-j, находящейся

в поле

силы тяжести {V=gz), при отсутствии внешних

перемен­

ных сил (F = const) из выражения (8.8)

находим

 

<Э<Р1

Piv i

і

і

^i = - P i - J j

§

Р і ^ - Ь const,

но

+ 2UlClk ch k (z+A,) cos (ta - e*) с точностью до малых второго порядка. Поэтому

( / 7 , ) . = £

= [С,р! (a — kUx) ch АЛ, — pig-fl] cos (Ах a*f) —^—|-const,

причем

в правой части

мы, как и выше, z = \ заменили на 2 = 0

(p2)z=zi

= \C2p2(? — kU2)ch

kh2 — p2ga]

cos {kx of)

^—(-const.

Так как pi = pa при z = £, то получаем

 

С,р[ (а ,%£/,) chft/z,— p{ga=C2p2

(а — ,Ш2 ) ch kh2 — p2g"a. (9.42)

Исключая отсюда Сі и С 2 с помощью выражений

(9.40) и (9.41)

и выражая о через скорость волны

с и угловую частоту ft, т. е. по­

лагая a = kc,

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

скорости

волны с:

 

 

 

Pi (с -

U,f cth ААх + Р2 (с -

£/2 )2 cth kh,= P l ~ Р 2 £ .

(9.43)

Не выписывая громоздкую общую формулу для с, ограничимся случаем, когда Ы и /г2 настолько велики по сравнению с длиной волны к = 2я/к, что можно принять cth khj. = cth kh2 = 1. Тогда

c _ P i ^ i +

P2^2 +

]

/

g ( P i — Рг) X

Pip2 (ЕЛ — ^ г ) 2

/g

Pl +

P2

~

Г

2те ( P l + p2 )

(p, + p2 )2

' ^ ' /

Знак плюс при корне

соответствует волнам,

распространяю­

щимся в направлении

оси Ох, а знак минус — волнам, идущим на­

встречу оси Ох. Из формулы

(9.44)

видно, что с вещественно,

пока

(£Л - Щ2

<

-У— .

 

(9.45)

v

1

~> ^ -

2лріР2

4

>

Если это условие

не выполнено

и с становится комплексным

числом, то решение,

которое нами

рассматривалось, перестает

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

циальное движение с поверхностью раздела

невозможно, оно дол­

жно перейти в какой-то иной вид движения

(заметим, что за счет

действия не учтенных в нашем

выводе капиллярных

сил, этот пе­

реход происходит при меньших

X, чем дает

неравенство (9.45), но

он обязательно

происходит).

 

 

 

Установив

факт физической

неосуществимости

рассматривае­

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

суждений не можем сделать никаких выводов об осуществимых

формах движения: для этого требуется

нелинейная теория. К во­

просу об осуществимых

движениях, имеющему

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

значение, мы вернемся в главе V I .

 

 

 

 

 

Если принять £Л = £/2 = 0, получим из формулы

(9.44)

 

 

V

~ъг

Р 1 +

Р 2

-

 

(9-46)

При р2 = 0 эта формула,

как и следовало

ожидать, переходит

в формулу (9.29). Из формулы

(9.46)

видно,

что, так как плот­

ность воздуха в 770 раз меньше

плотности воды,

то для системы

вода—воздух

 

 

 

 

 

 

 

 

l / "

Pi — Р2 ^

і

1

 

769

 

 

У

Pi +

Р2

 

770

770

'

 

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

10.Внутренние волны в стратифицированной жидкости

Вморях, озерах и водохранилищах плотность воды в разных слоях различна из-за разных причин, главные из которых — изме­

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

Составим уравнения внутренних волн, считая жидкость идеаль­

ной. Обратимся к уравнению Эйлера (5.12) и уравнению

неразрыв­

ности (5.10)

 

p - ^ f + g r a d / 7 - p F ^ O , - i -- |P - + d i v v = 0 .

(10.1)

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

4 - = - g - + v - g r a d P = 0

(10.2)

и уравнение неразрывности приобретает вид (5.18)

 

 

 

 

div v = 0 .

 

 

(10.3)

В общем

случае

 

 

 

 

 

 

- F = g r a d

+

[Q; v ] ,

 

(10.4)

где V — потенциал силы тяжести; К — вектор внешних

сил на еди­

ницу

объема,

отличных от силы

тяжести; Q — вектор

угловой

ско­

рости

вращения Земли. Последний член выражения (10.4), будучи

умножен на

р, дает так называемую

силу Кориолиса. Эта

сила

О час

100

300

Рис. 9. Внутренние волны в

Гибралтарском проливе — изменения

солености воды 16—18/V 1969 г.

существенна для внутренних

длиннопериодных приливных волн

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

будем ее учитывать.

 

 

 

 

 

 

Пусть в равновесном состоянии,

т. е. при установившемся те­

чении, К = 0,

v ^ v * ,

р=р%, р = р*. Для такого

течения уравнение

Эйлера п уравнения

(10.2) и (10.3) дают:

 

 

рй : (у*,

grad)v^ +

2p* [м # ,

v*] +

grad/»* +

p,grad V = 0 ,

(10.5)

 

 

v*

• g r a d p ^ O ,

d i v v ^ = 0 .

(10.6)

Положим

теперь

 

 

 

 

 

 

 

v = v i , . +

Av, />==р*+Д/>, р = Р * + л Р .

 

где Av, Ар, Ар — возмущения. Очевидно, что

 

 

 

 

 

р*»ьр,

Р * » Д Р -

 

 

Кроме того, будем считать, что локальные производные До и Др по t много больше субстанциональных (конвективных) производ­ ных, т. е.

 

 

 

d Av

д Av

d Др

д Др

 

 

 

 

 

dt

~~~ dt

 

' dt

~

dt

'

 

 

Для Ду такое предположение вполне соответствует предполо­

жению о том, что член v2 /2 в уравнении

(9.5)

мал. Таким

образом,

уравнения (10.2)

и (10.3) могут быть линеаризованы

по отношению

к возмущениям

и записаны

с

учетом

уравнений

(10.5)

и (10.6)

в следующем

виде:

 

 

 

 

 

 

 

 

 

Р* ^

^

+ P^v*grad)Av + p,; (Av,

g r a d ) v , +

 

 

 

 

+ gradA,o-f Apgrad У + К = 0 ,

 

(10.7)

~^-

+ v*

gradAp + Av • grad P : i : =0 ,

d i v A v = 0 .

(10.8)

Направим ось z

вертикально вниз и поместим начало координат

на свободной поверхности в состоянии равновесия. Тогда

 

 

 

 

 

V=-gz.

 

 

 

 

(10.9)

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

отсутствует, т. е. У^ =

0, а плотность зависит

только от

глубины,

т. е. р^. = р # ( г ) .

Если v^=0 ,

то

Av = v

и соответственно

Аи = и,

Av = v, Aw=w.

Система

(10.7)

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

 

 

 

 

d v

' gradAp-Apgk-\-K=0,

 

 

 

(10.10)

 

 

г* dt

 

 

 

 

 

 

 

 

 

 

 

l^-JrWA±_=Qi

 

 

d i v v = 0 .

 

(10.11)

Чтобы свести эту систему к одному уравнению, умножим ска­

лярно (10.10) на ід/дх + ]дІду.

Получим

 

 

 

 

 

d (да до).

д*Ар .

d*Lp

дКх

,

 

дКУ _ _ п

Р * ^ { ^ + ^ } + - Ш Г + - о ^ + - д 1 Г + - д у - - 0 -

Подставим сюда из уравнения

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

 

 

 

да

.

dv

 

dw

 

 

 

 

~dx

г

dy

 

dz

 

 

 

и продифференцируем результат по t

и z.

Будем

иметь

d*w

dp*

d*w

. /

д*

,

д*

)

d^Lp

dt* dz*

dz

dt* dz

~ \

dx*

1

dy*

j

dt dz

 

d*

(

dKx

,

dl<y

 

 

 

 

 

WdTK-Tx-r-nr,-*-

 

 

 

 

п п 1 9 ч

( 1 0 Л 2 )

(10.13)

(10Л4>

Проекция уравнения (10.10) на ось z есть

Дифференцируя это уравнение по t и подставляя d/Ap/dt из (10.11), найдем

 

 

 

д'2

 

(

 

д-w

 

.

rfp*

. дК,

\

 

, 1 Л 1

 

Подставляя

результат

в

(10.12),

получаем

окончательно

 

 

 

0Л-2

^ , 2 - Г

^-9

J

dfl ~

&

\ UA.-2 ^

dy2

 

 

dzdfi

~~

 

 

1

(

/

 

 

аг

\

дКг

 

 

д2

I д К

г

дКу

\ \

 

 

 

~

н

\ \

дх'2 "г" ду2 У

^

 

<^<?г

d*

'

cly

/ J '

 

У1"-1')

где Г = —

^ *

. При K = 0 это уравнение

есть уравнение

собст-

 

 

РФ

az

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

венных колебаний, называемое уравнением Релея.

 

 

 

 

 

Уравнение Релея допускает периодические частные решения

 

 

 

 

 

 

w=W(x,

 

у,

г)еш.

 

 

 

 

 

(10.18)

 

Подставляя

это значение w в

 

(10.17)

при К = 0,

получаем

 

 

 

 

 

 

 

 

 

 

S a - ^ + r - S - J - o .

 

 

(.0.19)

 

Частные решения

уравнения

(10.19)

будем

искать в

виде (ме­

тод

Фурье)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W{x,

у,

z)=X(x)Y(y)Z{z).

 

 

 

 

 

(10.20)

 

Подставляя W из (10.20) в

 

(10.19)

и

деля затем

результат

подстановки на W, получаем

уравнение

 

 

 

 

 

 

 

 

 

 

X"

,

У"

 

 

ш2

 

/

2"

, -г,

Z

 

 

 

 

 

 

 

 

*

1

У

 

 

gl-

2 " ( ^ | ~ + Г 4 " ) = 0 -

 

 

(Ю.21)

 

Поскольку

первое

и

второе

слагаемые

могут

зависеть

только

от х и у соответственно,

а

остальные

члены

уравнения — только

от

z,

уравнение

(10.21)

может

быть удовлетворено

только, если

 

 

 

 

 

X"

2

 

,

 

У"

02

 

,

 

 

 

 

 

 

 

 

 

 

i a i = c o n s t , —77—=pi=const,

 

 

 

 

 

 

 

 

 

«2

 

(

Z"

^

Z

 

2 = c o n s t

 

 

 

(10.22)

 

 

 

 

 

 

 

 

 

 

r - | ^ ) = T

 

 

 

и если, кроме того, выполнено условие

5 З а к а з № 428

65

Для того чтобы X и Y были ограниченными,

at и Pi, должны

быть мнимыми числами,

т. е. ai = ia, fii = ifi и, следовательно, у 2 =

= — (а 2 +(3 2 ) . Теперь уравнения (10.22)

приобретают вид

Л"' +

а 2 ^ = 0 ,

К' +

р 2 К = 0 ,

 

Z" +

r Z ' + (g-Г -

to2)

v Z = 0 .

(10.23)

Здесь

 

 

 

 

 

 

 

o2 +

p2

 

(10.24)

 

 

 

 

 

— положительное число. Первые

два уравнения

[10.23) дают

XV =• С cos (ax - f a) cos фу + />),

(10.25)

где С, a, b — произвольные постоянные.

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

условиям. Считая

дно

водоема

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

и имея

в виду,

 

 

 

 

 

что на свободной поверхности внут­

-Z

 

 

 

 

ренние волны

практически

отсутст­

 

 

 

 

вуют, следует принять, что при 2 = 0

EN

 

 

 

 

 

 

и при z = H вертикальная составля­

 

 

 

 

 

ющая скорости w равна нулю. Из

//

 

 

 

этого следует

 

 

 

 

\\

 

\

 

Z ( 0 ) = Z ( 7 / ) = 0 .

(10.26)

\

 

Запишем

теперь

третье

уравне­

I і

 

\

*

ние (10.23) так:

 

 

 

'

z

/

'

\

 

 

 

 

 

 

 

 

 

 

 

7

 

/

 

 

• г - ( р . - § ) + * С * > * - о .

Л

 

 

 

 

 

<7( £ )= Р*(£ Г-ш2 ).

(10.27)

\\ / A

 

 

/

 

Если q(z)

> 0

при всех z

и если

Рис. 10. Внутренние волны при

при этом Z > 0 , то

 

 

различных значеннях

параметра v.

 

d

 

I

dZ

 

 

 

 

 

 

 

dz

 

V*

dz •)<о

 

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

при монотонно

растущей

плотности р* (z)

произ­

водная dZldz монотонно уменьшается. Если

же Z < 0 , то dZjdz

монотонно увеличивается. Вид функций Z(z)

для

случая

<7(г)>0

приводится на рис. 10. Подчеркнем, что далеко

не при всяком зна­

чении v уравнение

(10.27), как

и всякое

однородное уравнение,

имеет решения, удовлетворяющие краевым условиям. Такие реше­

ния С у щ е с т в у ю т Т О Л Ь К О При НеКОТОрЫХ Определенных V = V i , V2, . . . ,

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

вечающие первым трем

собственным числам, а пунктиром — част­

ные решения уравнения

(10.27), не удовлетворяющие краевому ус­

ловию на дне водоема.

 

Если q{z) < 0

при всех z и если при этом Z > 0 ,

то

 

 

 

 

 

 

а

Р*

dZ

> 0

 

dz

dz

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

производная

dZ/dz

монотонно увеличивается от

поверхности ко дну. Краевое условие на дне при этом не может быть удовлетворено. Если 2 < 0 , то dZ/dz монотонно убывает с глу­ биной. Частные решения Z] и Z%, соответствующие этим случаям, представлены на рис. 11. Можно найти линейную функцию реше­ ний Zi и Z2 ) которая сама будучи решением удовлетворяет вместе с тем краевому условию на дне. Но она не только не будет удов­

летворять краевому

условию Z(0) = 0 , но и будет

иметь

экстремум

на поверхности и, следовательно, будет

 

-Z

 

 

 

отвечать не внутренним, а поверхностным

 

 

 

 

волнам.

 

 

 

\

\\

 

 

Таким образом, критерий существова-

Z

 

 

ния внутренних волн (при горизонталь­

 

 

 

ном дне) есть o(z)>0, т. е.

 

 

1

 

 

 

<VgT

=N{z),

(10.28)

 

 

 

 

 

\\

 

или иначе период Т внутренних волн дол­

 

 

жен удовлетворять

условию

 

 

 

\1

\

 

 

 

т>

(10.29)

Н

 

1

 

 

 

 

N(z)

 

 

 

частные

 

 

 

 

Рис. П. Некоторые

 

 

 

 

 

решения

уравнения

внут­

Величина

N(z)

в теории

внутренних

 

ренних

волн.

 

волн

называется частотой Вяйсяля.

 

 

 

 

 

Из

(10.28)

вытекает, что

в однородной жидкости,

характери­

зующейся тем, что р = const и N = 0, внутренние волны невозможны. Мы рассмотрели случаи, когда критерий выполняется при всехг

или не выполняется при всех z.

Если же

критерий (10.28) выпол­

няется не при всех z от 0 до

Н,

а только

при некоторых, т. е. при

0 < Z i < Z < 2 2 < # , В н у т р е н н и е

ВОЛНЫ ВОЗМОЖНЫ В СЛОЄ Z i < Z < Z 2 , они

носят название нехарактерных внутренних волн.

Обратимся к важному частному случаю экспоненциального рас­

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

 

 

 

 

Р* (2) = = Р*(0)ег °г ,

Г ( 2 ) = -

dz

= r

0 = c o n s t .

(10.30)

 

 

 

 

 

Отметим прежде всего

очень важный

факт,

что весьма

малое

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

например, РцДЯ) = (1+є)р„.(0),

где

е — величина порядка

Ю - 4 .

Нетрудно подсчитать, что в этом

случае при глубине водоема

# =

= 10 м

частота Вяйсяля

есть

~ 1 0 - 2

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

(10.29)

будет Г > 6 2 8 с « 1 0

мин.

 

 

 

В рассматриваемом частном случае третье уравнение

(10.23)

приобретает вид

 

 

 

r + r 0 Z 4 - № - ( o 2 ) v Z = 0 ,

(10.31)

и при условии, что

 

 

 

4 ( £ Г 0 - о г Ь - Г 5 > 0 ,

 

имеет решение

 

 

 

Z = e

2

(С, coste-|-CsinXz),

 

в котором

 

 

 

Х=4-У"4 ( g - r 0 - u . 2 ) v - r g .

(10.32)

Из условия Z( 0 ) = 0

находим СІ = 0, а из условия Z ( # )

= 0 по­

лучаем

 

 

 

Ш=пк,

п=1, 2, . . .

(10.33)

Таким образом, решение, удовлетворяющее краевым условиям,

будет

 

 

Z(z)=Ce

2 sia-^ff-.

(10.34)

Из (10.32), (10.33) и (10.24) нетрудно найти следующее соот­ ношение:

Объединяя теперь формулы (10.18), (10.25) и (10.34), получаем следующее выражение для w.

w=De

2 sin

cos (a-x-\-a) cos (§y-\-b) cos (<&t-\- c).

(10.36)

Здесь D и с — новые

произвольные постоянные. Частное

реше­

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

ний (10.36) с разными сс и

Р,

удовлетворяющими соотношению

(10.35). Решение (10.36) дает

стоячие (целлюлярные) волны: на

линиях

 

 

тк -f- -jr

а

тъ-\- -= Ь

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

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

1.Факт существования внутренних волн, их принципиальное отличие от волн поверхностных (подчеркнем, что внутренние волны, так же как и поверхностные, могут вызываться непериодическими возмущениями).

2.Внутренние волны возможны только в стратифицированной

жидкости, при этом достаточно самого незначительного различия

вплотностях верхних и нижних слоев.

3.Теория, рассмотренная в данном параграфе, имеет в виду случай, когда поступательное движение жидкости отсутствует. Однако анализ волн на поверхности раздела двух жидкостей, вы­ полненный в предыдущем параграфе, показывает, что поступатель­ ное движение жидкости имеет принципиальное значение и спо­ собно приводить к качественно новым явлениям неустойчивости поверхности раздела. А так как скачок плотности на поверхности раздела есть предельный случай постепенного, но значительного изменения плотности на небольшом интервале изменения z, то сле­ дует ожидать, что учет течений приведет к существенным новым результатам.

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

Более подробное изложение теории внутренних волн и всего

круга проблем, связанных с этой теорией, имеется в монографии В. Краусса [34]..

11. Течение с поперечной циркуляцией

Рассмотрим уравнение (5.12)

J * f + ( v, v ) v = F _ _ L g r a d / ? .

( и л )

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