Файл: Тареев, Б. А. Динамика бароклинных возмущений в океане.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 21.10.2024

Просмотров: 129

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

мущенном

состоянии движение отсутствует,

имеет место

уравнение

гидростатики:

 

 

 

Фо

е Ро = °-

(3.1.1)

 

дг

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

ut fv =

1

--------г р*.

 

 

Ро

Щ+ fu —

1

Ру-

 

 

Ро

wt

gp

1

= ---------Рг,

 

Po

Ро

гидродинамики .относительно (3.1.1) с точностью до малых

(3.1.2)

их + Vy+ Wz = 0,

Р< + - -iBsdz -w = 0.

Здесь как

обычно и,

и — компоненты скорости по

гори­

зонтальным

осям х, у\

w — компонента

скорости

вдоль

оси г,

направленной вертикально вверх; t — время; f — пара­

метр

Кориолиса;

g — ускорение силы тяжести; р — возму­

щение

давления;

po = po(z)— потенциальная

плотность

в не­

возмущенном

состоянии;

р — возмущение

плотности.

Ниж­

ние индексы указывают на дифференцирование по соответ­

ствующей переменной.

Жидкость

считается неоднородной,

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

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

с изменениями давления в волне,

пренебрегается. Рассмат­

ривается только устойчивая стратификация, т. е. dpo/dz^.0. Из системы (3.1.2) удобно исключить все неизвестные, кро­ ме вертикальной компоненты скорости w, для которой наи­

более естественно

могут быть

сформулированы

граничные

условия. В результате такого

исключения

получаем

урав­

нение

 

 

 

 

 

 

 

A wtt +

N2Д2 w -f /2 wzz

 

 

JL

Фо

,

 

 

 

 

 

Po

dz

 

А =

д2

д2

А2

дх2

л .

 

(3.1.3)

 

дх2

ду1

 

ду2

 

Здесь N — частота Вяйсяля—Брента; А, А2 — трехмер­ ный и двумерный операторы Лапласа.

При выводе уравнения (3.1.3) было использовано так называемое приближение Буссинеска, т. е. пренебрегалось изменениями p0(z) в уравнениях движений по горизонталь­

136


ным осям. Это

приближение для океана

всегда справед­

ливо1.

формулировке

граничных

условий.

Если

Перейдем к

2 = —# = const — дно, то, очевидно,

 

 

 

 

 

о» (— //) = 0

 

 

 

(3.1.4)

Если глубина бесконечна, то

со(2 = оо,

х,

у,

t) =0.

Ниже

будет показано,

что если N достаточно быстро

стремится к

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

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

ственно из

уравнений

движения (3.1.2)

и из

условия

dpldt = 0 при

z=l{x, у,

0, где р = р0 + р,

£=£(к,

у, t)

отклонение свободной поверхности от невозмущенного поло­

жения

2 = 0. Линеаризируя это условие с учетом (3.1.1), по­

лучим

 

 

 

 

 

 

dp

__ _др_ .

dp0

dp

— gPoW’

 

dt

dt

dz

dt

или, исключая давление с помощью уравнений

движения (3.1.2),

(wtta+ f 2w)z g Д2и) = 0 при z = £

=^0.

(3.1.5)

Система (3.1.2) и уравнение (3.1.3) инвариантны относи­ тельно сдвигов по горизонтальным осям координат и оси времени, и, следовательно, для w (и других зависимых переменных) всегда существуют решения, имеющие вид плоских волн,

ш = ср ( г ) е Нк^ +куу- ш ) .

( 3 . 1 . 6 )

Здесь ф(г) — амплитудный множитель; к» — круговая

часто­

та. Любое решение может быть получено путем

интегриро­

вания (или суммирования

в случае

дискретного

спектра)

элементарных

волн

вида

(3.1.6)

в

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

волновых

1 Интересно

отметить,

что

если не

использовать приближение

Бус-

синеска и во всех уравнениях

считать р0

существенно

зависящей

от г,

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

(3.1.3) получим:

 

 

 

 

 

 

 

 

Д wtt + IV2 Д 2 w+ / 2 wzz

-------wztt 0 .

 

 

( 3 . 1 . 3 )

 

 

 

 

 

 

Я

 

 

 

 

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

 

при обычном

граничном

усло­

вии (3.1.8) при

//= оо и f2<co2 одним

из

решений

этого

уравнения

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

извест­

ное дисперсионное

соотношение <в2=gk,

не

зависящее

от

плотностной

неоднородности жидкости.

 

 

 

 

 

 

 

 

 

137


чисел (kx, ky). Подставляя (3.1.6) в (3.1.3), (3.1.4) и (3.1.5),

получим:

 

0,1 (д т ~ * !) ф _ / ! Л5Г +

4’л'!'г = 0

(3 IJ)

 

при Z — 0;

ф ~ >

----у = 0;

у (—

Н) — 0.

(3.1.8)

Тильда над k в

(3.1.8)

означает,

что если мы положим

k = 0,

то придем к теории длинных волн,

соответствующих

гидростатическому

приближению в третьем уравнении ис­

ходной

системы

(3.1.2). Поскольку уравнения системы

(3.1.2)

инвариантны относительно вращений вокруг

верти­

кальной оси 2 в (3.1.7), (3.1.8), во все последующие соот­ ношения kx и kv будут входить в виде симметричной комби-

нации k2 = kx +ky. Уравнения (3.1.7), (3.1.8) немедленно приводят к одному важному следствию, касающемуся коле­

баний

с инерционной

частотой /.

Полагая со-*-/,

вместо

(3.1.7)

получим

 

 

 

 

 

( N * - f 2)k2Ф =

0.

 

(3.1.9)

Так как в океане, как правило,

 

то в (3.1.9) либо

k2 — 0,

либо <р=0. Поскольку ищется

нетривиальное

реше­

ние, то k2= 0, а ф=ф(г)

— произвольная функция. Но тогда

из (3.1.6) и уравнения

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

системы (3.1.2) сле­

дует, что ду/дг=0, т. е. Ф — const. С учетом граничных усло­

вий (3.1.8) фн=0. Таким

же путем

получим,

что динамиче­

ское

давление р = 0.

Горизонтальные

компоненты скорости,

легко

определяемые

в

этом случае

из

(3.1.2), равны

и = и' (z)e~ift\

 

 

 

 

 

 

 

и = — iu' (z)

 

-£(»+—-)

 

 

 

= и’ (z) e

V

2

где и’ (z)

произвольная функция 2. Следовательно, инерционные коле­ бания с частотой f представляют .собой чисто горизонталь­ ное бездивергентное движение без градиентов давления, никак не связанное со стратификацией и архимедовыми си­ лами. («Длина волны» таких колебаний равна бесконечности

(k — О), определяемая формальным путем

фазовая скорость

cg= d(S)jdk также

бесконечно велика, однако групповая ско­

рость cg = da/dk,

очевидно, равна нулю, т.

е. никакого пере­

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

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

138


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

затухать на характерном расстоянии D ~ Y A J f , где Ае— коэффициент горизонтального турбулентного трения. В ра­

боте

(Ямпольский,

1961) эта величина интерпретирована

как

«длина волны

инерционных внутренних волн», однако

из сказанного выше следует, что инерционным колебаниям не может быть приписана определенная длина волны, а ве­ личина-' D является просто горизонтальным аналогом глу­ бины трения в смысле Зкмана. В связи с этим можно заме­ тить, что довольно распространенный термин «инерционные внутренние волны» едва ли можно признать удовлетвори­ тельным, так как, являясь предельным случаем внутренних волн, инерционные колебания не будут внутренними волна­ ми в строгом смысле слова. В частности, они не связаны ни с архимедовыми силами стратификации, характеризуемыми параметром N, ни с волновым переносом энергии. Выводы автора справедливы только для случая, когда изменения параметра Кориолиса с широтой не учитываются. Вблизи предельных широт этот эффект важен.

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

меняется

в

пределах толщины термоклина. В этом

случае

за счет

кинематического

условия на

дне

w = u(dH/dx) +

v(dH/dy)

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

ния,

и если

и и v колеблются с инерционной

частотой,

то,

как

следует

из последнего

уравнения

(3.1.2), на

этой

же

частоте будут наблюдаться колебания плотности (или тем­ пературы). По-видимому, такие колебания температуры с инерционной частотой, описанные в работе (Haurwitz, Stommel, Munk, 1959) !, были отмечены на береговом скло­ не одного из Бермудских островов. Удачным для наблюде­ ний обстоятельством следует считать то, что один из термо­

метров

находился на дне

на

глубине 50 м (глубина

сезон-

1 Работа

Гаурвитца,

Стоммела,

Манка

(1959), а также

работы

Краусса

В.

(1959)

и

Раттри

М.

(1956)

опубликованы в

сборнике

«Внутренние волны».

М.,

«Мир»,

1964.

 

 

 

139