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

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

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

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

Добавлен: 21.10.2024

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

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

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

dW 2

tg б

dP2

 

p_______ H i.

/

D

dy

U ( \ — c)

D2 U ( 1 — c)

 

 

 

 

/2 [1 —

f (1 —

c)2]

(P*~Pi)-

 

 

 

 

 

g 'D 2

 

 

При выводе (2.6.8) сделано пока только одно приближе­ ние: множитель Vj при |3 в (2.6.5) выражен через Pj путем ис­ пользования геострофического соотношения. Если следовать обычному способу отфильтрования и вместо (2.6.7) подста­ вить в (2.6.5) соотношения геострофичности

 

 

 

1

dPj

 

(2.6.9)

 

 

и,- --------------

1

 

 

 

 

р/

ду

 

 

 

то получим систему уравнении:

 

 

 

d2P г _

Г k 2

,

Р _ _______ ^ i .

f

Pl =

dy2

[

^

£/(!+■

c)

 

U (1 -c c)

 

 

 

 

P

(Pi-P*),

 

 

 

 

 

g 'D i

 

 

(2.6. 10)

 

 

 

 

 

 

 

d2P 2

Й2

____ P_________ t £ i _

f _

 

dy2

 

 

/7(1—c)

D2

£/(1 — c)

 

 

 

 

 

72

(P,~Pl).

 

 

 

 

 

g 'A l

 

 

 

Квазигеострофическая модель такого вида с учетом гори­ зонтального изменения скорости основного течения была под­ робно исследована Педлоски (Pedlosky, 1964). Из сравнения (2.6.8) и (2.6.10) видно, что предположение квазигеострофич­ ности эквивалентно предположению AZ)/Dj<Cl и F = 0. По­ скольку нас в первую очередь интересует влияние негеострофических эффектов на поведение градиентно-вихревых волн, возвратимся к уравнению ,(2.6.8). Если предположить, что изменения глубины слоев ДD (см. рис. 19) на ширине течения

малы, по

сравнению с самими глубинами (т. е. AD/D<Cl)r

то вторые

члены в левых частях (2.6.8) будут также малы

по сравнению с первыми членами. В самом деле, порядок первых членов A/ у /2, а вторых

tg6

A Pj

AD

APj

Dj

l ~

Dj

l2

Таким образом, вторыми членами в (2.6.8) можно пренеб­ речь; и, кроме того, заменить Dj=Dj(y) постоянными сред­ ними значениями, что в очень большой степени упрощает ма­ тематическую сторону задачи. Используя эти упрощения и

Т'

99



формулу Маргулеса (2.6.1), получим вместо (2.6.8) систему уравнений с постоянными коэффициентами:

(PPj_

 

k2 —

g’Di

U

Px =

dy2

 

1 + c

 

 

 

 

 

g'Dx

 

■ k2

U2 (1 + c ) 2

(Pi-P*),

 

 

g'Dx

 

 

 

 

_L

(2.6 . 11)

 

 

 

 

 

d2P 2

k 2

g 'D 2

U

P*

dy2

 

 

 

 

 

 

 

Г P

 

— k2

U2 (1 — c)2

(P,~Pi).

. g'D2

 

g 'D 2

 

 

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

Pj (0) = Pj (l) =

0, или

dPj

__

dPj

=

0.

(2.6 . 12)

dy

у=о

dy y=t

 

 

 

 

 

Тогда собственные

функции системы

(2.6.11) будут

иметь

вид:

 

 

 

 

 

 

 

Pj = Aj sin К у, или Pj = Aj cos К у, где %п = - у к

Подставляя собственные функции в (2.6.11), получим си­ стему линейных однородных уравнений относительно Aj. При­ равнивая к нулю определитель этой системы, получим харак­

теристическое уравнение-для определения с:

(% + <*! - J ^ - )

+ С )» ]-

(2.6.13)

-----L [1 - cx2F2 (1 - с)2} [ т 2 + а2 -

+

+ i - [ l - a 2F 2( l - c ) 2]

100


Здесь введены обозначения:

 

k\i = fig'

-Щ-- Rf =

«; = P/klf,

Е,- =

U2/g' ~Y~,

ms = l 2n 4 / 4 .

В этих обозначениях

 

F = R~l = UVg' - f -

= (U2 - Uj)2/g'H

внутреннее число Фруда или обратная величина числа Ричардсона в двуслойной модели, H = 2D — полная глубина. Если положить Dl= D 2=D, т о индексы 1, 2 в (2.6.13) можно опустить. В этом случае, полагая Е-Я) (квазигеострофическое приближение) и /->оо ,(m-v0) — очень широкое течение, воз­ вращаясь к неподвижной системе координат, получим резуль­ тат Н. Филлипса (Phillips, 1951).

 

c = U ___ ±

+ а)

 

 

 

262 (1

 

 

 

±

 

(2.6.14)

 

2а (1 -р а)

 

 

 

Здесь

c = c r+ ic j= (®г+Шг)/й — размерная

комплексная

ско­

рость

в неподвижной системе координат,

сг — (фазовая

ско­

рость волн,, шг=Сг& — мнимая часть комплексной частоты — коэффициент возрастания волн.

Уравнение (2.6.13) в общем случае содержит и бароклинную неустойчивость и неустойчивость Гельмгольца, так как внутренние гравитационные волны не были отфильтрованы. Имея в виду, что в дальнейшем приложения модели будут от­ носиться к узким океаническим течениям, )|3-эффектом, как

показывают оценки, можно пренебречь.

В этом случае при

D{— D2= D и R = 0 уравнение (2.6.13) приводится к биквад­

ратному:

 

а (т + a) F с4 + [3 aF — (m -f a ) (m +

а + 1)] c2 -f

-f- [(m -pa)2{m + a) — a(m -j- a)F -f aE] = 0 . (2.6.15)

Если обозначить x — c2, то решение этого уравнения может быть записано в виде:

х =

За F + а) + а 1) ,

+ a) F

 

101


±[За F + а) -f- а -+- I)]3

 

[2а +

а) FY

 

 

_ f

+ а)2 - +

а) ,

аF [1 — (от + а)]

_ ^ g

16)

[

а (т + а) F

*

а + а) F

 

 

Очевидно,

отрицательные

или

комплексные

значения

х

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

В формуле (2.6.16) величина т (для первого собственно­

го значения М) равна я2Ц2ко= Lo/4l2

 

L0 = 2n!k0 = 2% j / V - f - -f~l-

I

Если задаться характерными параметрами Гольфстрима на участке от Флориды до м. Гаттерас (Stommel, 1963): (pi—р2)/р = 2 -10—3; H=2D=1000 м; /= 1 0 ~4 сек - 1; 1=60 км,

то Y g D /2 ^ 3 м/сек; L0=200 км. Тогда /га» 1. Поскольку ве­ личина сдвига скорости (1)2U\)j2 может достигать 1,8 м/сек, то параметр негеострофичности Г » 0,4. В районе Флоридско­ го пролива величина F может достигать единицы. Расчет по

формуле (2.6.16) (с учетом нижнего знака перед корнем) по­ казывает, что в этом случае для всех значений F < 1 и а < 1 бароклинной неустойчивости нет, т. е. соответствующие зна­ чения все действительны. Для рассматриваемой области из­ менения параметров подкоренное выражение в (2.6.16) всегда положительно. Более подробные расчеты (которые мы здесь не приводим) показывают, что с уменьшением значений т возникает бароклинная неустойчивость, причем в первую оче­ редь дестабилизируются длинные волны, соответствующие малым а. Однако следует иметь в виду, что геометрические свойства рассматриваемой модели таковы, что величина т

102