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

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

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

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

Добавлен: 21.10.2024

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

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

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

u = ---------dp/dy,

v = ------ dp/dx

в уравнение вихря. Тогда

р/

и

 

Р/

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

из (2.9.5), (2.9.6)

(2.9.3)

 

 

 

f i f - U y )

'dp'

 

 

 

 

 

g(H+l)

dt

 

 

\ U (f — Uy) f

r r

(f~Uy)

 

(2.9.7)

gH

 

 

yy

H

 

 

Если включить g

в определение Я с учетом. (2.9.3) и поло­

жить p — p(y)expik(xct),

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

амплитуды р(у)

примет вид:

 

 

 

 

 

 

 

W ( y ) - c ]

d2p

k2p 1М.:иу>

p

+

 

dy3

 

gH н

+

[

-(t/g

-f)-

— Um]p = 0

(2.9.8)

 

н

 

dy

yyY

 

 

(здесь опущена черточка над р). На границах

течения у\ и

г/г нормальная компонента скорости должна обращаться в

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

дает:

 

Р Ш = Р Ш = 0-

 

(2.9.9)

Если вне

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

области

(z/b уг)

£/=const и

Я = const, то,

как следует из

(2.9.8), за

пределами \у\, у2)

амплитуда волновых возмущений экспоненциально затухает: р = ех р (± % ) где k = ]/& а + &о , ko=--f2/gH, а знак выби­

рается, чтобы обеспечить затухание при у-^+оо. Это приводит к другому виду граничных условий: •

± kp при у — ух или у — уг.

(2.9.10)

Границы могут находиться и на бесконечности. Уравнение (2.9.8) при условиях (2.9.9) или (2.9.10) может, вообще го­ воря, иметь решения с комплексными собственными значения­ ми C= Cr-\-iCi, которые приводят к неустойчивости. Условия, необходимые для неустойчивости, легко, доказать, следуя обычной методике. Обозначая выражение в последних квад­ ратных скобках (2.9.8) через К (у), полагая p — pr+ipi,

(U—c) _1 = ai+ ta2, где а \= (Я—с2) | U—с\~2; a2=Ci\U—c\-2,

отделяя вещественную часть (2.9.8) от мнимой, получим

d2Pr

И2 +

n f - V y )

__aiK{y)

pr— a1K{y)pi = 0 ,

[dif

 

gH

 

 

d2Pi

k2

f i f - U y )

UiK (y)

Pi ~j- a%K (y) Pr = 0.

dy2

 

gH

 

 

127


Отсюда следует (после умножения первого уравнения на * Pi, второго на рти вычитания):

j а $ ( у ) \ p f d y - „ : 0.

■ (2.9.11)

у.

 

Интегрируя

(2.9.11)

в пределах у\, у2, при любом виде гра­

ничных условий >(2.9.9) или (2.9.10) получим:

- j —(Pi

dy

-

Рг- ^ ~ ) = - а 2К(у) (рг2+ рЬ- (2.9.12)

dy \

 

dy ]

Следовательно, если с^Ф0 (неустойчивость), суммарный гра-

диент завихренности К(у) = U - U y )

дН ■— и,

должен ме­

Н

ду

уу

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

На материковом склоне А Н /Н ~\ и при U~ 1 м/еек, Дг/~100 км (оценки, характерные для Гольфстрима) первый член в выражении для К {у) будет на порядок больше кривиз­ ны поперечного профиля скорости.

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

Конечно, само по себе изменение глубины без поперечного

изменения скорости течения не может служить причиной неус­ тойчивости. В этом смысле полученный критерий не являет­ ся достаточным. В частности, при U = const, после умножения (2.9.8) на р— ргipi, интегрирования по частям и использо­ вании граничного условия (2.9.9) получим:

(U— с) J j |p '|2 + (& +

|Р I2} d y = j К(у)\р I2dy.

У,

У1

Откуда

 

у2

У2

О | [\Р' |2+

\p\2\}dy= ^ К(у)\ p\2dy, (2.9.13)

Ul

 

128


Уг

( U - c r) ==

_____

i K ( y ) \P |2 dy

(2.9.14)

y\_____________________

 

 

У2

 

 

j {lpl3+ (Й3+ _£г)|р|2}^

Hi

Следовательно, если U = const, то с*=0 и течение устойчиво. Это очевидно также из энергетических соображений, так как для нарастания волновых возмущений, как известно, необхо-

димо

 

Уг

%Uydy^>0, где т = —рш>.

выполнение неравенства j

 

 

ci=0,

соотношение (2.9.14)

Когда Я = const и, следовательно,

дает удобную

возможность для приближенного вычисления

дисперсионных соотношений вариационным методом.

С

другой

стороны, из (2.9.8)

при U—const следует, что

если

Я (у) — симметричная функция в

интервале \у\, г/г)

(симметричный относительно гребня подводный хребет или

ложбина), то

К (у) =

----- — дН/ду

— антисимметрична. Сле-

 

 

н

 

и антисим­

довательно, решениями (8) будут симметричные

метричные в

(уи г/2)

функции: р(у) — ± р {—у).

Тогда из

(2.9.14) ввиду антисимметричности

К(у) следует U—с=0,

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

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

Когда Я меняется линейно, глубина по мере отдаления от берега стремится к бесконечности. В действительности, ма­ териковый склон имеет конечную ширину, поэтому ниже будет рассмотрен достаточно типичный случай

Я = Я0(1 —е*У).

(2.9.15)

Начало координат возьмем на береговом крае материко­ вого склона. Как видно из (2.9.15), в открытом океане (г/-*— оо) Я = Я 0. Для U—const (и, следовательно, веще­ ственных с) уравнение (2.9.8) перепишем в виде:

d?P ,

f._______ дН_

р

= 0

(2.9.16)

dy2

(с - U) Н

ду

 

 

 

или, обозначая

%=y.f{U—с)-1, a2= f 2/gH, -из

(2.9.15) и

(2.9.16) получим:

 

 

 

 

 

d2p

К

а2

Р =

0.

(2.9.17)

~dy*

1- е кУ

1-<*»

 

 

 

9 Б. А. Тареев

129



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

(2.9.17): р = 0 при у = 0, р — ограничено при г/->— оо (2.9.18).

Если ехр(—ху)_— х, то вместо (2.9.17) получим

хЦх — 1) ^ Е - + х ( х — 1) -+ [* ,' — a'2x) — k'2(x— \)]p=0.

 

 

 

 

 

 

 

 

 

 

(2.9.19)

Здесь а’2= а2х~2, k'2x~ 2,

к' — кх~2.

Граничные

условия для

(2.9.19):

 

 

 

 

 

 

 

 

 

 

 

р (1) = 0,

р (оо) — ограничено.

 

(2.9.20)

При р = хтф (х),

т = ± V k ’2+ к’ уравнение (2.9.19)

приво­

дится к виду:

 

 

 

 

 

 

 

 

 

 

х(х — 1) ср" +

1) (2т + 1) ф' + (Г — а'2) ф = 0.-

(2.9.21)

Определяя параметры а,

(3, у равенствами

 

 

 

a = m ± V k’2+

а'2,

р = т ± У k’2+ а'2,

 

 

у =

+ 1, = ± V k'2 + х'2 ,

 

(2.9.22)

приводим уравнение (2.9.21) к гипергеометрическому:]

 

 

х ( х — 1)ф" + [а -Ь Р +

1) л: — у)] ф' + <*Рф = 0.

(2.9.23)

Общее решение (2.9.23) в рассматриваемом интервале:

 

 

Ф = хга F (а,

1 + а — у,

1 + а — р; г -1) +

 

 

ч- х~&F (Р,

1-гР — У,

1-гР —а; х~]).

(2.9.24)

Здесь F ^ Л, В , С; — j — гипергеометрический ряд,

1 < | х | <оо.

Из

(2.9.24) получаем:

 

 

 

 

 

 

 

 

р = Clxm~aF ( А,

В, С; — ) + с2хт~Р F (Л', В',

С'; — ).

 

 

 

 

 

 

 

 

 

 

(2.9.25)

Но

независимо

от

знака

т,

если в

выражениях 2.9.22 для

а, р выбрать верхние знаки, получим

 

 

 

т а — — V k '2 + а'2 < 0 и т — р = V k'2 + а’2 > 0 ,

поэтому условие ограниченности р (оо) приводит к С2= 0 , и решение для р с учетом использовавшихся выше обозначений:

130