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

книги из ГПНТБ / Палагин, Э. Г. Основы гидромеханики учебное пособие для метеорологов

.pdf
Скачиваний:
13
Добавлен:
20.10.2023
Размер:
8.24 Mб
Скачать

их фильтрации, то в соответствии с конкретно намечаемой целью намечаются и пути фильтрации.

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

§ 2. МЕТОД МАЛЫХ ВОЗМУЩЕНИЙ. ПАРАМЕТРЫ ВОЛН

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

Vi— V i ~ р = р-\-р', р= р + р /. Т=Т-\-Т'.

(16.1)

Здесь величины со знаком тильда ( ~ ) относятся к среднему дви­ жению, а отмеченные штрихом соответствуют малым возмуще­ ниям. При этом:

4 '- « . 1,

4 -

I . X - ч< 1.

« 1•

116.2)

г

р

у

Т

 

Подставляя (16.1) в уравнения движения (13.1) будем иметь

di'i

Op,

-

Ог,-

.1

с)с.*

Ос-

Ov'i

 

Ot

' 01

■ 1

Ох,

J

Ох i

'

'Ц/

OXj

~

 

 

1

/ Ор

 

Ор'

\j

or

 

 

 

 

0 -L г/ 1 Ох,

 

ох,

 

 

Прежде всего. преобразуем член с давлением в правой Для этого учтем, 1что

1

1

1

\

1

 

\

 

 

г.

0 Д ‘

f

)•

"

: (

Н 4 - ]

 

ь

 

 

 

 

(Здесь использовано разложение в ряд функцпи вида

= 1—а' + л'2+ . .., где мы ограничились линейным приближением). 11оэтому

1

| ар

Op j __

1

Op

1

Op'

o'

Op

и + [>'

\rlKi

O.Kj

о

0 x ;

у

Ox j

!.r

Ox-,

 

 

 

/ _

J p .

 

 

 

 

 

 

 

О2

( ) X [

 

 

 

 

190

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

dv> .

г

dvj ,

do'

.

.

dV;

. - dv)

< > Р

 

dt

-

°J

r

 

4

 

dt

J «*,

 

 

dXj

 

dX;

С)Х;

 

 

1

dp'

,

 

?'

dP

£833

 

 

 

и

dxi

 

 

p2

dXj

 

 

 

 

 

 

 

Уравнения для основного движения запишутся как обычно

Ah.

<b'i

А р .

r)t

dXj

r)X

Вычитая это уравнение из предыдущего, получим уравнение для возмущений в виде

dvj

+Vj Ah.

dv'i

1 dp'

dp

(16.3)

-

vj dXj

dx,

dx,

dt

dXj

 

Если характеристики основного движения нам известны, то (16.3) представляет собой линейные уравнения относительно не­ известных v/, // р'. Мы добились этого отбрасыванием нелинейных членов второго порядка малости. Такого рода процедура носит название линеаризации, а сами уравнения являются линеаризиро­ ванными по отношению к исходным.

Аналогичным образом можно получить линеаризированные уравнения неразрывности (13.2) п притока тепла (13.3), которые соответственно принимают вид:

 

 

do'

,

/ V o

do v 'j

0;

16.4)

 

 

dt

dx;

 

 

 

dx,

 

 

or

A L

 

 

dT'

/

 

dp

dp'

Vj

 

-I- О

dt

dx,

 

dXj - (

dt

dx,-

vi dx,

 

 

 

 

 

dT'

 

 

(16,5)

 

 

 

 

dX ;

dx,

 

 

 

 

 

 

 

Поскольку

волновые

 

возмущения

являются

периодическими

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

 

v 'i —

e 0 i s i n ( k ! x i ~ - o i ) ,

 

р' — Р„sin (/е,-А-,. - 4),

 

 

16.6*)

 

р‘ — р0 sin (А’;А,. — at),

 

г — Т0sin (k , x , —ot).

 

.4

\

( Напомним, что

V

/г,-х-, J.

191

В (16.6) величины, отмеченные индексом «О», представляют собой амплитуду соответствующей характеристики.

Из условия периодичности вытекает, что:

J

2 -

2 -

j '

(1 6, /)

1

^ i

где Т — период, т. е. время

повторения

одинаковых

состояний, а

— расстояние по различным осям между одинаково расположен­ ными па профиле волны точками, например между двумя, наибо­ лее удаленными от осп х точками А и В (рис. 16.1), лежащими на соседних гребнях (точка С лежит на подошве волны)

Z

рис. Ifl.l

Общая длина волны ?.= j /

/у, i — ~ — круговая частота,

А’.- именуются волновыми числами. Общий аргумент Лулу - о/ на­ зывается фазой. Форма волнового возмущения какого-либо эле­ мента в (16.6) может быть получена приравниванием ее постоян­ ной. Например, для давления, если обозначить // = /(.V|. ,\-2, .г3, О

/То sin (k;Xj - О/) —/ (Л'|, Л'2, Л'3, / ) = COITSt.

(16.8)

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

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

192

Составляющие фазовой скорости с,- легко найти через посред­

ство соотношения

 

 

 

 

 

 

df

df

dx1

Of

dx2

Of

dx о

0,

dt

+ dx,

dt

+ dx.

df

Ox.,

dt

получаемого дифференцпроваппем (16.8) по времени. Поскольку

dx;

СУТЬ производные от координат профиля по времени, то они находятся как результат дифференцирования неявной функ­

ции нескольких переменных, т. е.

 

 

 

 

 

dXj

dfjdt

 

 

 

 

dt

Of.dXj

 

 

 

 

С учетом равенства (16.8) будем иметь

 

 

 

dxi

*

 

Ч_

 

 

С;— ' dt

 

 

Т

 

 

НЛП

 

 

 

 

 

 

/

V

=

Т '

(16.9)

- T

V '

S

" -

 

 

 

Формула (16.9) утверждает, что фазовая скорость равна отно­ шению длины волны к периоду. Отметим, что в силу специфики за­ дач часто существует зависимость сг = сг (ки k2, k3), именуемая

дисперсионным уравнением.

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

Любую периодическую функцию можно разложить в ряд Фурье, представляющий собой сумму бесконечного числа гармо­ ник с различными амплитудами и частотами, и изучать процесс с этих позиций. Разумеется, практически рассматривается конеч­ ное число членов этого ряда, исходя при этом из заданной точ­ ности. Так или иначе, практически мы, как правило, имеем дело с совокупностью волн, так называемым волновым, пакетом. Тогда для любой функции / в (16.6) следует записать

ОС

СГтО •

(16.10)

/ошЛл (kimX

ш —I

 

 

Установим фазовую скорость группы волн, иначе именуемую

групповой скоростью.

13 Зак. 1 1 2

193

Вначале рассмотрим две волны, имеющие одинаковую ампли­ туду п близкие по своим значениям величины волновых чисел k t, к-,' н частот о, о', т. с.:

 

 

/ 1=/« sin (kjXi at) ,

 

 

 

 

h — fn sin

(A’/.vf --- a't).

 

 

Суммарный эффект

 

 

 

 

 

 

 

/ = / | + / 2= 2 /о cos

 

J 4 -V,-

 

 

t 1 sin >

 

 

А’; Ь А/

a

t

a

 

 

В с и л у близости к;

н к/,

а также сг п а'

множитель

 

 

¥ = 2 /о cos

A’i -

к{

о — а

б

(16.

 

 

,----- -V,

 

-

слабо меняется н его можно

рассматривать

как

переменную

амплитуду, значения которой меняются

в пределах от 0 до 2 /0 на

s k i —kj'

 

~

 

 

 

 

 

длине - ц-

---- ^----=

------вдоль каждои оси.

 

/

+*

 

«i

ас,•

 

 

 

 

Общая скорость перемещения обоих волн может быть охарак­ теризована как скорость точки с каким-либо постоянным значе­

нием F (А], х2, д'з, /)=const.

Отсюда,

как и

ранее, находим компо­

ненты фазовой скорости

 

 

 

 

 

 

 

dXj

_

 

dF/di _ л,

о/

(п

dt

~

 

dF/dXi

А>,- —/?,■' ‘

Если же имеем

непрерывную

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

групповая скорость должна

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

как

предел

 

 

Нш

Аз

 

 

( 16. 12)

 

 

 

 

 

 

 

А к.;

 

О

 

 

 

Очевидно, что групповая скорость отлична от фазовой скорости отдельной волны.

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

йс

f == V /ош-ехр | i(k\m Х|-|-А’2ш-V’2-f- А’зш А’з -- Лш( ) | .

(16.13)

т!

(Напомним, что согласно формуле Эйлера, еМ! = ао$ </+/ sin у, т. е. вещественная часть дает косинусоидальную волну, мнимая — си­ нусоидальную, i — комплексная единица).

194

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

§ 3. АКУСТИЧЕСКИЕ ВОЛНЫ

Упругие волны с диапазоном частот 20ч-20'103 герц, которые воспринимаются человеческим ухом, называются акустическими или звуковыми волнами. Существенно, что последовательные скач­ ки сжатия и разрежения, представляющие собой звуковую волну, происходят настолько быстро, что обмен теплом между областью волновых движении и средой происходить не успевает, и процесс протекает адиабатически. Таким образом, уравнение притока тепла имеет вид (9.21):

, , = v

;16.14)

Кроме того, силой тяжести в данном случае можно пренебречь, что, как мы убедимся далее, связано с большими числами Фруда,

1

характерными для этого процесса. Таким образом -рг- <£. 1.

Если источник звуковых колебаний находится в покоящейся среде, то

-

др

д р

dp

dp

■О,

(16.15)

'

дх1

дх.,

д.А',

дх2

 

 

а если игнорировать силу тяжести:

др

_

др

(16.15')

дх3

~

дх3

 

В этом случае (16.3) и 16.4) упрощаются и соответственно при­ нимают вид:

 

 

dv', - =

1

др'

(16.16)

 

 

dt

 

р

дх,

 

 

др'

 

dv\

= 0.

(16.17)

 

 

dt

+ р

- dXj

н

dv[

 

dv< '

 

 

Напомним, что —— — V

дХ;

.)•

 

 

 

dXj

/ i

 

 

Поскольку адиабатическое течение соответствует случаю баротропной жидкости, то можно записать

д р

d p

д о

d x i

d p

d x t '

13*

195

Логарифмируя и дифференцируя выражение (16.14), в свою очередь получим

dp __ у. do

Откуда

Р

 

 

 

 

 

 

 

 

dp_

•//;

 

(16.18)

 

do

 

 

 

 

 

 

Поэтому

 

 

 

 

 

_др_

у.р

до

 

 

дх,

о

д.V/

 

Подставив сюда разложение (16.1) с учетом (16.15) и (16.15'),

будем иметь

 

 

 

 

др’

 

- - : - ( л : Л (1

(!'/

дх,

1д. ■' \ дх,

~0jr7

 

 

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

др ’

у.р д о '

д о’

дх.

 

дХ;

дхi

Здесь введено обозначение ч2-

у.р .Из

(16.18) ясно, что у.р всегда

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

растает и наоборот, так что -^£->0. Это обстоятельство и отмече­

но тем, что величина а взята в квадрате.

можно

представить

Таким

образом, в итоге уравнения (16.16)

в виде

 

 

 

 

 

 

 

 

 

 

 

 

Д Д = ^ 4 1 Д £ 1 .

 

 

О

 

 

 

р

дх,

 

 

 

 

Пели продифференцировать каждое из уравнении

(16.1'J)

по х,-

и сложить результаты, то получим

 

 

 

 

 

дЧ\

,

г/“ Г2

, д2г'з

__

а 2 ( д2 [/

.

d‘V

, dV

 

dt дх1

 

dt Ox-,

dt дх-

 

дх?

 

дх:

дх:

 

В свою очередь,

дифференцируя

(16.17) но

/,

получаем

зави­

симость

 

 

 

dh:\

д~ г'-,

d2v'z

1 d V

dt дх,

dt дх.

dt дх.

О '~dF

196

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

3V

= я 2А(/

(Лр' — лапласиан).

 

(16.20)

—flf-—

 

Наиболее часто встречающаяся звуковая полна является сфе­

рически симметричной, так

что в

рассмотрение следует

принять

только одну координату г.

При этом (см. Б.

26)

AV .

2 до'

()г'

~ir

н уравнение (16.20) примет вид

 

 

 

 

 

 

 

c/V _

,

/ г>у

2

Ф'

 

 

(16.21)

dtf

~а

 

г)г

г

дг

 

 

 

 

 

 

Решение этого уравнения будем искать в виде

 

 

 

</=-

V

'

ехр | i (кп, г —- (7„, /) |

 

 

(16.22)

/лИ I

 

 

 

 

 

 

(■‘1»1 = const).

В данном случае амплитуда обратно пропорциональна радиусу, ибо, помещая источник в начало координат, мы там должны иметь

особенность первого порядка.

 

по / и г

и подставив

Выполнив все дифференцирования (16.22)

результаты в (16.21),

получим

 

 

 

со

д

ехр [i(k,„r

- гг,,,/)] (з=„

k‘f„a-) =

0.

V

т I

 

 

 

 

 

Последнее

выражеппс может

обратиться

в нуль, если только

 

 

=

 

 

(16.23)

С учетом (16.7) последнее равенство можно преобразовать и представить в виде

Из (16.9) следует, что фазовая скорость всех волновых гармо­ ник равна постоянной величине

ст = ± а .

(16.24)

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

197

Общее решение уравнения' (16.2!) получим если в (16.22) под­ ставим (16.23) и (16.24). Оно имеет вид

00

4

[i к’п, ( г - о/)] +ехр

\ikm(r-fn/)] |.

/ — V

— "г f exp

m l

1

 

 

Или, используя формулы Эйлера и отбирая члены с вещественной частью,

{•

X

“ [ cos А*„, (/* cii) -(-cos km

| .

^

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

Введенное нами понятие скорости звука относится к невозму-

.

Г 7Р

 

щепной среде, гак что а — |/

— . При наличии общего движения

под ней следует понимать величину а\— I

- L~. Это, так иазы-

 

I

Р

ваемая, местная скорость звука, представляющая собой скорость распространения колебания относительно среды в какой-либо точке. Из (16.18) следует, что скорость звука повышается с умень­ шением упругости, ибо при этом одному и тому лее изменению плотности соответствуют большие изменения давления. Так, ско­ рость звука в воде больше скорости звука в воздухе. В пределе, если среда абсолютно жесткая, а-* °с.

С учетом уравнения

состояния p = RpT скорость

звука в газе

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

 

 

 

 

 

 

 

 

a = yxRT.

 

 

(16.25)

Мы видим, что она зависит лишь от физических констант газа

и его

температуры в °К.

Для

воздуха, в

частности, и=1,4;

R — 287

м2■сек~2• град- 1

и

а =

20,1 уТ

м-сек~]. С

повышением

температуры скорость звука возрастает.

Поскольку

в атмосфере

в большинстве случаев

температура с высотой

уменьшается, то

скорость звука ь нижних слоях выше, чем в верхних. Поэтому зву­

ковые лучи,

идущие от наземного источника,

искривляются

с высотой.

 

 

 

 

Наконец,

оценим влияние силы тяжести на

распространение,

звука. Для

воздуха при 0°С а — 332 м/сек. Считая а характерной

скоростью,

 

2 -а.

характерным линейным разме­

а длину волны /.= ----

 

 

 

 

а- п

ром, можем число Фруда представить в виде Fr =

07. -2— . Даже

При a = 20 сек~', чему соответствует минимальное значение fr, оно

равно ~100 н, следовательно,

1- Таким образом пренебреже­

ние влиянием гравитации вполне оправдано.

§ 4. ГРАВИТАЦИОННЫЕ ВОЛНЫ

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

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

этого вытекает, что v2 — v2' = 0 и, кроме

того, все производные

— =0. Во-вторых, считаем, жидкость

несжимаемой, т. е

ОJCо

 

р— р— const, f/ = 0. Последнее условие означает, что в получаю­ щихся при .этом упрощенных уравнениях не находят отражения акустические эффекты, ибо звук может распространяться только в сжимаемой среде. Тем самым мы отфильтровываем звуковые волны и рассматриваем интересующий пас процесс в чистом виде, без включения не интересующих пас в данный момент звуковых колебаний.

В результате этих предположений уравнения (16.3) и (16.4) упрощаются п принимают соответственно вид:

(ь;х

.

dv,.

,

dv,.

 

dv'x

 

dv'x

1

dp\

 

сП

 

 

' '■*

dz

Vx

dx

bz

dz

p

dx ' (16.26)

dv'

- и

%

,

dv

+ ?,v

dv'z ,

 

dvi

1

dp'

(16.27)

dt

+ ' - I F

 

 

dz

p

dz

 

d.x

 

 

 

 

 

 

 

 

,

do'y

 

0.

 

 

(16.28)

 

 

 

 

 

dy

~

 

 

 

 

 

 

 

dx

 

 

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

свободной поверхностью. При этом будем считать, что два потока

с плотностями pi = pi и из—р2* движутся друг над другом, причем

* Индекс 1 относится к верхнему потоку, а 2 к нижнему.

199

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