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

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

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

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

Добавлен: 21.10.2024

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

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

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

или

 

 

 

 

 

 

 

 

 

 

 

 

Я ,=

Зо* +

-±-<г„.

 

 

0.3.11)

Подставляя

это значение'

Rc и

(1,3.10),

найдем,

что

значе-

ния о2

не

зависят от Qy:

 

3

 

Нетриви-

 

а3---- —a2s2 = 0.

альный корень

этого

уравнения

з

s2 (s =

ля, п — 1,2,

= —

3 ...). Полагая

л=1

(х= я),

из (1.3.11) при

Qy = 0 опять

получим для Rc рэлеевское значение. Следовательно,

гори-

Рис.

2.

Два

вида

зависимости

критического

числа

Рэлея

R c для

предельных

случаев

высоких

R Z(Q)

и

низких

R V{Q)

 

широт.

По

осям

координат

нанесены

десятичные

логарифмы

R c и Q.

На

вспомогательной вертикаль­

ной

оси

нанесены

безразмерные

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

размеры

 

конвек­

ционных

ячеек

 

(X/h).

 

Линии

XZ(Q),

k y (Q)

показывают

зави­

симость

этих размеров от пара­

 

 

 

метра

Q

 

 

зонтальная составляющая угловой скорости вращения Земли не влияет на величину циркуляционных ячеек, но критическое число Рэлея линейно растет вместе с Qy.

~ -Таким образом, влияние горизонтальной и вертикальной составляющих параметра Кориолиса приводит к увеличению критического числа Рэлея (хотя и по разным законам), но влияние fz и fy на горизонтальные размеры циркуляционных

ячеек существенно различно.

расчетов кривых

На

рис.

2

показаны

результаты

Rc = Ry{Q)

и Rc=Rz(Q)

в логарифмическом масштабе, а

также

кривые

X=XZ(Q) и

X=A,?/(Q) в

«полулогарифмиче­

ском масштабе», т. е. по оси ординат для X — кривых отло­ жены значения безразмерной длины волны X/h в зависимости от логарифма Q. Поскольку горизонтальная компонента па­ раметра Кориолиса не влияет на величину циркуляционных ячеек, кривая l = A,j/(Q) есть просто прямая, соответствующая постоянному рэлеевскому значению X/h—4.

Следует заметить, что один частный случай рассмотрен­ ной нами задачи, а именно случай высоких широт, был рас­ смотрен Л. С. Гандиным (1947) применительно к мезо-

22


метеорологическим явлениям и С. Чандрасекхарой (Chand­ rasekhar, 1953), применительно к астрофизическим пробле­ мам (конвекция в атмосферах звезд). В следующей работе (1954) С. Чандрасекхар (Chandrasekhar) рассмотрел также стабилизирующее действие магнитного поля нй конвекцию для частного случая, когда векторы силы тяжести, угловой скорости вращения и напряженности магнитного поля парал­ лельны. Следуя Чандрасекхару, легко показать, что ввиду малой электропроводности морской воды и малых значений напряженности земного магнитного поля влияние магнитного

поля на океаническую конвекцию практически

неощутимо.

Вскоре после работ

автора появилось

исследование Куо

(Кио, 1961), ссылка

дается на русский перевод (Куо, 1967),

в котором, однако,

также рассмотрен

только

предельный

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

Возвращаясь к нашей океанографической задаче, оценим реальную величину параметра Q, характеризующего враще­ ние Земли. Пусть характерная толщина конвективного слоя имеет порядок 104 см. Характерное значение молекулярной кинематической вязкости морской воды v~10~2 см2/сек. Так

как

/~ 1 0 -4 сек-1, то получим Q= /2/i4/v2«

1012. При

таких

больших значениях Q формулы (Т.3.7) и (1.3.11)

для

расче­

та

критических значений конвекции могут

быть

заменены

асимптотическими выражениями. Тогда для случая высоких

широт

(с учетом, что Rc = 3o2c) получим

Rc

3 ^я2

/з .

а для

горизонтального волнового числа:-

 

 

 

 

] //2 -J- тг

^ я2

'j .

 

 

 

Для низких широт

безразмерное

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

число

сохраняет рэлеевское

значение: l2 + m2 — s2jY 2,

а для

критического числа

Рэлея получаем Rc^ Q yl2.

Эти асимпто­

тические формулы

наглядно показывают

фундаментальное

различие влияния горизонтальной и вертикальной составля­ ющей угловой скорости вращения Земли на критические па­ раметры конвекции.

Оценим теперь фактическое значение числа Рэлея. Пусть

h = 104, у = а££=П0“1, 90—6i = 0,l град,

/е= 10-3, v= 10-2

CGS.

Тогда, согласно определению (1.3.6),

получим R = 1015,

что

намного превышает соответствующие

значения Rc, которые

могут быть получены из асимптотических формул. Представ­ ляется естественным при оценке критических параметров кон­ векции использовать молекулярные значения диссипативных коэффициентов k, v, как если бы движение начиналось из состояния покоя. Однако, когда фактическое число Рэлея

23


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

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

нографии

значения

vTyp6 = /Стурб = 102 сж2/се/с, получим

/? = 106. Но тогда мы должны использовать

то же значение

vTyp6 и в формуле для

Q, что дает Q= 104. Используя рис. 2,

найдем, что R соответствующее этому значению Q будет все-

таки меньше /?= 106, и конвекция все равно

должна иметь

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

§ 1.4. Сравнение теоретических результатов с данными наблюдений (некоторые замечания)

В изложенной выше теории уравнения для возмущений и граничные условия однородны, как это обычно бывает в линейных задачах гидродинамической устойчивости. Поэтому система решений вида (1.3.3) опреде­ лена с точностью до общего всем решениям постоянного амплитудного множителя. Чтобы определить этот множи­ тель, а вместе с ним и абсолютные значения всех неизвест­ ных величин, достаточно измерить значение одной из зави­ симых переменных в какой-либо точке в выбранной системе координат. Помимо определения критических параметров конвекции существенный интерес представляет определение абсолютных значений вертикальной компоненты скорости конвекции w, а также величины отклонений потенциальной температуры Ф от равновесного значения 0(z)= pz + ©o1. Только эта последняя величина может быть непосредственно найдена из данных по измерениям температуры в глубоко­ водных впадинах. Поскольку в решениях (1.3.3) физический смысл имеют только действительные части, для результи­ рующей температуры в размерных величинах имеем:

т = е (z) й(*, у, z) = pz + е0+

-f (0i — e0)fl cos — х cos — у sin — z.

(1.4.1)

п

п

п

 

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

1 В этом параграфе для удобства письма перейдем к новой систе­

ме обозначений: величины без штрихов будут обозначать размерные переменные.

24


производятся чаще, чем через 1000 м. Анализ немногочислен­ ных материалов наблюдений показывает, что наибольшая толщина конвективных слоев имеет место в Филиппинской и

Бугенвильской

впадинах,

где

на

 

 

 

 

 

толщину

слоя

 

приходится

по

 

 

 

 

 

пять наблюденных значений темпе­

 

 

 

 

 

ратуры.

Это позволяет с достаточ­

 

 

 

 

 

ной точностью

определить измене­

 

 

 

 

 

ние с глубиной отклонения изме­

 

 

 

 

 

ренной температуры от равновесно­

 

 

 

 

 

го значения 0(z). Поэтому в качест­

 

 

 

 

 

ве характерного

примера

мы

ис­

 

 

 

 

 

пользуем данные Б. Шульца, приве­

 

 

 

 

 

денные в табл. 1. На рис.

3 показа­

 

 

 

 

 

ны рассчитанные

по этим данным

 

 

 

 

 

потенциальные температуры относи­

 

 

 

 

 

тельно поверхности океана для об­

 

 

 

 

 

ластей отрицательного сверхадиаба­

 

 

 

 

 

тического

температурного

градиен­

 

 

 

 

 

та. В Филиппинской впадине наблю­

Рис.

3.

Распределение

дается повышение

потенциальной

потенциальных

темпера­

температуры с глубины 6000 м, где

тур

в Филиппинской

(/)

01 = 1,01°,

до придонного горизонта

и

Бугенвильской

(2)

измерений (9788 м), где

0о= 1,31°.

океанических

впадинах

в области

отрицатель­

Следовательно, равновесный темпе­

ного

сверхадиабатиче­

ратурный градиент

 

 

 

 

ского

градиента

 

 

р =; _Ё1—

— 0,8 • 10-6

град/см.

 

 

 

 

 

 

h

 

 

 

 

 

 

 

Отклонения температуры от равновесного распределения,

изображенного на рис. 3 прямыми линиями 0г0о, с большой степенью точности совпадают с полученным ранее теоретиче-

ским изменением, по синусоидальному закону

я

sin — г, где

O ^ z ^ /i, причем h — толщина конвективного

слоя, равная

3788 м для Филиппинской впадины и 3400 м для Бугенвиль­

ской впадины (рис. 3).

Для

определения

безразмерного

амплитудного множителя

O'

положим

в формуле

(1.4.1)

х = у = 0, z=h/2 (начало

координат на

дне,

так что

z=h/2

соответствует на рис. 3,1 расстоянию 1894 м от дна), и най­ дем из рис. 3, что соответствующее отклонение температуры

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

значения

при z = h/2

равно примерно

0,06°,

т. е. /01—0о/О= О,О6 град,

а так

как

(0О—0i) =0,3 град, то

0 = 0,2. Заметим,

что кривые на

рис.

3 вычислены при

пред­

положении

изменения температуры по

синусоидальному за­

кону (2.4.1)

при 0 = 0,2,

а точки на кривых соответствуют на­

блюденным

значениям

потенциальной

температуры. Теперь

25


можно определить величину вертикальной скорости. Исполь­ зуя (1.3.3), а также вторую формулу (1.3.9) при п = 0 и фор­ мулы (1.2.7), получим.

ш = ---- —ad cos — х cos — у sin— г.

(1-4.2)

h

h

h

h '

 

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

ратур

на расстоянии,

по горизонтали равном по крайней

мере

2яН/УI2+ т2,

так как данные единичной батометри-

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

части конвективной ячейки находится вертикаль

измерений

с координатами х = у 0.

Поэтому, полагая в формуле (1.4.2)

x=y = 0, z=hj2, можно

только заключить, что

максималь­

ная величина w будет не меньше по абсолютной величине, чем

wn = — oft.

(1.4.3)

0 h

 

Поскольку Филиппинская впадина расположена в низких широтах, для а можно принять рэлеевское значение сг=1,5л2.

Соответствующее число Рэлея -^ - /г4 =

10-3/г4. При

столь

kv

иметь, по-видимому,

большом h (3788 м) движение будет

турбулентный характер и величина

k

в формуле

(1.4.3)

должна иметь смысл коэффициента турбулентной теплопро­

водности

[k — v — Ю2- ^ - ) .

Тогда,

учитывая, что # = 0,2;

а=1,5 л2;

\

сек j

 

что tw0 = 0,8• 10—3

см/сек.

&—3,8-105 см, найдем,

Однако есть основания полагать,

что

значение k порядка

IQ2 CGS,

обычно

принимаемое

для

поверхностных

слоев

океана, для абиссальных областей слишком велико и оцен­ ки, следующие из теорий абиссальной термохалинной цир­ куляции, дают для величины k порядок 1 CGS. Если принять, что условия по горизонтали одинаковы и конвекция вызвана геотермическим притоком тепла, среднее значение которого известно (g=l,3-10-6 кал/см2сек), то учитывая, что равно­

весный градиент температуры

в наших условиях (3 = 0,8-

•10_6 град/см, можно оценить

«эффективный коэффициент

теплопроводности»:

 

х = 2L ~ 2CGS.

Р

Так как удельная теплоемкость Ср и плотность р равны в

единицах CGS примерно единице, то k = и =2см 2/сек,

срР

26