Файл: Тареев, Б. А. Динамика бароклинных возмущений в океане.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-3—10-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