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

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

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

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

Добавлен: 21.10.2024

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

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

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

условие неустойчивости — изменение знака градиента абсо­ лютного вихря zv= § d2uldy2.

В зависимости от характера профиля скорости и (у) бэтаэффект может способствовать либо сохранению, либо изме­ нению знака zy в области течения.

Рис.

29.

 

Значения

|т г-|

для

Рис. 30.

Значения

|tj[

для

Р = 0,1.

Кривые 1 и

3 соответст­

{$ = 0,1.

Нумерация

кривых

та

вуют

левой

и правой

полуплоско­

же,

что и на рис. 29

 

сти,

бэта-эффект

учитывается:

 

 

 

 

.В = 0,05.

Кривая 2 вычислена при

 

 

 

 

отсутствии

бэта-эффекта

(В = 0)

 

 

 

 

Можно заметить, что в окрестности .нейтральной кривой Кочина (малые числа Ричардсона или очень короткие волны) влияние бэта-эффекта не существенно и сказывается главным образом в появлении отличных от нуля (но малых, порядка 10~2) тг (на нейтральной кривой Кочина при отсутствии бэтаэффекта имеет место «принцип изменения устойчивости», т. е.

Tj = Tr = 0) .

Значения % для геострофической системы (2.8.17) были также рассчитаны для ряда значений а, {$ (см. Приложение).

В рассмотренной области изменения

параметров

значения

|Tj | , и особенно |тг|, заметно занижены по сравнению с точ­

ными значениями из системы (2.8.14),

что связано,

очевидно,

с потерей .первых производных при использовании стандарт­

ного геострофического приближения, приводящего

к

(2.8.17).

 

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

122


Если принять, что основным механизмом меандров Гольф­ стрима является бароклинная неустойчивость (а на это опре­ деленно указывают результаты наблюдений Д. Ханзена, при­ веденные в § 2.6), то на участке от Флориды до м. Гаттерас в связи с малой глубиной, большим сдвигом скорости число Ричардсона достаточно мало невозможно, течение находится вблизи нижней части нейтральной кривой на рис. 26. Ниже м. Гаттерас увеличение глубины и уменьшение сдвига приво­ дит к увеличению числа Ричардсона, что, согласно рис. 26, соответствует резкому усилению неустойчивости и, в конечном счете, интенсивному образованию меандров. Такой характер движения в целом соответствует данным наблюдений в Гольф­ стриме. В численном решении использовался вариант метода

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

ных q= (dPi/dr), dP2/dr\) система (2.8.14) переходит в систему 4-х уравнений первого порядка. Для отхода от особых точек

т) = ± 1 (отгонка) используется подстановка q = w{r\)P. Реше­ ние уравнения (нелинейного) для матрицы w(r\) в окрестно­ сти особых точек ищется в виде степенных рядов с учетом характера особенностей. Первый член ряда определяется с

учетом граничных условий: w(—Ц )=ш 0, w(-\-l)— w0. После отхода от граничных особенностей уравнение для ®(ri) оказы­ вается неудобным для интегрирования, так как оно может со­ держать особенности внутри области интегрирования (напри-

мер, там, где Р = 0), поэтому вводится другое, эквивалентное первому уравнение для некоторой функции (матрицы) ф(г]), не содержащее отмеченных недостатков. Прогонка слева и справа и условия склейки в нуле приводят (после отделения вещественных и мнимых частей) к условию D(а, (3, В, тг, т,) = = 0, где D — некоторый определитель восьмого порядка.' (В процессе вычислений методом покоординатной минимиза­ ции минимизировался квадрат модуля этого определителя с точностью |Дт| = 10-310-5.)

§ 2.9. Градиентно-вихревые волны на материковом склоне океана

Уже было отмечено (§ 2.1), что с точки зрения принципа сохранения потенциального вихря эффект влияния топографии дна в баротропном океане качественно аналогичен |3-эффекту. Для того, чтобы выделить эффект из­ менения глубины в чистом виде, будет рассмотрен баротропный океан |(Тареев, 1971). Обобщение на двуслойный случай не представляет принципиальных затруднений. Для квазидву-

123


Т а б л и ц а собственных значений

1. Случай В = 0. Система (2.8.14). Собственные значения, принадлежащие «второму собственному решению» , отмечены индексом (II).

Значения т для гравитационной волны, приведенные для примера, отмечены индексом (Г. В.).

а

Э

К 1

К-1

а

Р

I V I

 

6

0

0,000

0,505

9

0,2

0,530

0,388

6

0

0,000

0,000 (II)

10

0,2

0,548.

0,375

7

0

0,006

0,494

10

0,2

0,002

0,275 (11)

8

0

0,111

0,500

11

0,2

0,563

0,359

10

0

0,172

0,505

11

0,2

0,005

0,281 (II)

11

0

0,189

0,511

12

0,2

0,575

0,345

12

0

0,203

0,516

12

0,2

0,003

0,278 (II)

6

0,1

0,358

0,586

5

0,3

0,437

0,375

7

0,1

0,402

0,564

7

0,3

0,516

0,302

8

0,1

0,437

0,541

9

0,3

0,531

0,256

8

0,1

0,006

0,353 (II)

10

0,3

0,544

0,234 •

9

0,1

0,467

0,525

11

0,3

0,553

0,214

9

0,1

0,008

0,403 (II)

10

0,1

0,494

0,511

12

0,3

0,561

0,194

10

0,1

0,007

0,430(11)

6

0,4

0,459

0,205

11

0,1

0,512

0,498

7

0,4

0,478

0,156

И

0,1

0,007

0,447 (II)

8

0,4

0,494

0,109

12

0,1

0,534

0,486

9

0,4

0,503

0.555

12

0,1

0,008

0,451 (II)

10

0,4

0,453

0,000

6

0,2

0,445

0,450

11

0,4

0,422

0,000

7

0,2

0,480

0,430

12

0,4

0,406

0,006

8

0,2

0,508

0,406

9

0,5

0,497

0,000

6

0,5

0,414

0,203

6

0,7

0,148

0,000

6

0,6

0,230

0,003

6

0,8

0,089

0,000

6

0,9

0,039

0,011

 

6

1,0

0,001

0,000

100

0,02

0,052

0,491

 

9

0,4

0,496

0,000 (Г. В.)

9

0,6

0,205

0,000

 

 

 

 

 

 

I2. Случай В = 0 ,05. Система

(2!.8 .14).

Значения хг для левой

 

и правой полуплоскости берутся с соответствующим знаком.

а

Р

\%г\

K I

а

Р

K I

l Til

2

0

—0,017

0,000

 

10

0,1

- 0 ,6 0 3

0,550

2

0

0,010

0,000

 

10

0,1

0,417

0,477

4

0

-0,067

0,475

 

10

0,1

—0,078

0,420(11)

6

0

—0,125

0,525

 

12

0,1

—0,628

0,519

8

0

—0,217

0,548

 

12

0.1

- 0 ,0 3 3

0,442(11)

10

0

—0,280

0,586

 

2

0,2

—0,063

0,000

12

0

—0,325

0,614

 

4

0,2

—0,439

0,539

12

0

0,172

0,395

 

4

0,2

0,280

0,486

2

0,1

0,048

0,000

 

6

0,2

—0,516

0,469

'4

0,1

—0,445

0,700

 

8

0,2

—0,562

0.417

4

0,1

0,142

0,516

 

10

0,2

—0,592

0,378

6

0,1

—0,519

0,638

 

12

0,2

—0,612

0,345

6

0,1

0,250

0,516

 

2,004

0,3

—0,008

0,000

8

0,1

—0,567

0,589

(И)

4

0,3

—0,453

0,420

8

0,1

—0,064

0,345

6

0,3

—0,512

0,338

6

0,3

—0,383

0,000

(II)

8

0,3

- 0 ,5 5 0

0,278

4

—0,4

—0,452

0,312

 

6

0,4

—0,491

0,203

8

0,4

—0,516

0,105

 

10

0,4

—0,414

0,000

2

0,5

—0,272

0,000

 

4

0,5

—0,441

0,195

6

0,5

—0,436

0,031

 

8

0,5

—0,330

0,000

10

0,5

,—0,302

0,000

 

100

0,02

- 0 ,3 9 4

0,519

124


Продэлжгние таблицы

 

3.

Случай В =

0,05. Геострофическая система (2.8.17.)

 

а

Р

IM

1т/ 1

а

3

\ ь \

1 \-|

3

0,1

—0,135

0,348

4

0,1

—0,136

0,434

5

0,1

—0,137

0,465

6

0,1

—0,145

0,469

7

0,1

—0,168

0,452

8

0,1

- 0 ,2 0 8

0,447

9

0,1

—0,238

0,456

10

0,1

—0,256

0,463

11

0,1

—0,272

0,466

12

0,1

—0,283

0,466

3

- 0,3

—0,041

0,210

4

0,3

—0,038

0,280

5

0,3

—0,034

0,284

6

0,3

—0,031

0,262

7

0,3

—0,003

0,227

8

0,3

- 0 ,0 0 3

0,188

9

0,3

—0,002

0,187

10

0,3

—0,002

0,110

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

< 2 ' 9 Л )

Здесь

 

 

d/dt =

+ и ----h v

I = ди/дх ди/ду,

от

дх

ду

/ — параметр Кориолиса, Я — глубина. Если считать /=const,

d /d t= 0,

то

при

из (2.9.1) следует результат Экмана

v\/H = 0,

т.

е.

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

стицы движутся приближенно вдоль изобат. Если пренебречь кривизной изобат, то можно считать Н = Н ( у ) , если ось х направить вдоль изобаты. Тогда для малых неустановивших-

ся отклонений от невозмущенного течения U— const,

направ­

ленного вдоль оси х,

получим

(при

 

£<Cf=const)

из

(2.9.1):

 

dj =

/

дН

 

 

 

(2.9.2)

 

-----V.

 

 

 

dt

Н

ду

 

 

 

 

Здесь и далее —d- - = ——

—р Ы

,

Если ——

дН

= const

dt

dt

 

дх

 

Н

ду

 

и движение не зависит от у

(t,=dv/dx), полагая о-7-ехрik(x

ct),

где k, С волновое число и фазовая скорость в направ­

лении

х, из (2)

получим дисперсионное соотношение, анало­

гичное

формуле

Россби:

c = U — (---------- —

W&2. Если

^

 

 

\

н

ду

}/ -

изобаты (и невозмущенное течение) направлены зонально, то, учитывая изменения / в приближении |3-плоскости, вместо пре­

дыдущей формулы получим: c==U— (р

дН

---- ^ -----г~)В от-

крытом

океане Н ~ 5

 

ду

км, ДЯ~1 км, Аг/ — 1000 км, что при

 

f

дН/ду * 10-13 CGS,

т. е. величину того

f ~ 10-4

сек-1 дает н

125


же порядка, что и |3. Таким образом, планетарные |3-волны должны существенно видоизменяться под влиянием изменений глубины. В районах'больших уклонов дна, например на ма­ териковом склоне, градиент «топогенной завихренности» (по

терминологии Экмана)

/

гдН

5

>Р, поэтому в дальней­

 

 

Нду

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

1 дР

_

(2-9.3)

fU(y) = - — — =

ду

1 =

рУ

 

 

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

Линеаризированные относительно (2.9.3) - уравнения дви­ жения:

ди

тг

ди

,

др_

 

dt

U

------ Ь v

дх

 

 

дх

dy

 

dv

I/

 

+ fu=------

ЁЕ.

(2.9.4)

~dt

дх

 

 

 

 

Р

ду

 

дают уравнение вихря (при /=& const):

 

ди

dv_ = 0

(2.9.5)

dt

U yyv + (f — U y>( дх

ду

 

(значок у внизу’'означает дифференцирование по у). Линеа­

ризированное

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

можно записать в

виде:

 

 

 

 

 

 

ди

dv

д(Н + |)

 

_ d | _

 

~дГ + ~ду Ня + ®

ду

■ V dt

 

 

д(Н + 1)

1

dp

(2.9.6)

 

 

ду

"ip

dT

 

 

 

Здесь l(x, у,

t) — отклонение

от'невозмущенного уровня и

p — gpl

(р= const — плотность,

р —отклонение от невозму­

щенного давления Р ).

 

 

(гироскопических)

Для

выделения градиентно-вихревых

 

волн в чистом виде целесообразно ввести квазигеобтрофическое приближение, фильтрующее гравитационные волны:

126