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

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

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

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

Добавлен: 21.10.2024

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

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

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

частиц, равномерно выпадающих из поверхностных слоев океана.

В метеорологии задача обтекания неровностей подсти­ лающей поверхности бароклинным потоком была достаточ­ но подробно исследована многими авторами. При рассмот­ рении океанографического варианта этой задачи мы будем следовать в основном работам автора (Тареев, 1964, 19656), в которых применен метод А. А. Дородницына (1938, 1950), но изменим соответствующим образом постановку задачи. Расположим начало координат на свободной поверхности океана и направим ось г вертикально вниз. Пусть форма неровностей дна и все движение не зависят от у и постоян­ ная скорость невозмущенного течения и направлена вдоль оси х. Считая морскую воду неоднородной несжимаемой жидкостью и пренебрегая эффектом вязкости, который мо­ жет быть существенным только в тонком придонном погра­ ничном слое (и в данном случае не представляет интереса), запишем уравнения задачи (считая движение квазистационарным и пренебрегая параметром Кориолиса)

, ди . ди \ др

 

р

и

ow

 

W-

 

 

др

 

gp-

(3.4Л)

 

~д7

 

 

 

дг

 

 

 

 

 

дг

 

 

 

 

 

 

 

 

 

- ^ + * 1

= 0,

 

 

 

 

 

 

 

 

 

дх

дг

 

 

 

 

 

 

 

 

и

ds

 

ds

 

г.

р =

,

ч

 

 

 

 

— +

w—

= 0 ,

р (s).

 

 

 

 

 

дх

 

дг

 

 

 

 

 

 

 

Здесь

и, w — составляющие

скорости

вдоль

осей х, z ;

р — плотность;

s — удельная

энтропия

морской

воды;

р — давление; g

— ускорение силы тяжести.

 

 

Пусть теперь:

 

 

 

 

 

 

 

 

 

 

 

 

и =

и +

и'

(х, z),

р = р0 (2) +

Р' (х, г),

 

 

 

w =

w' (х,

z),

s =

s0 (г) -f- s' (х, z).

 

(3.4.2)

 

 

 

 

Р =

Ро (2) +

р' (х, z).

 

 

 

 

Здесь

величины

со

штрихами — возмущения

основного

течения,

вызванные

обтеканием

малой

неровности

дна

| = 1(х),

причем | ( —оо)=0,

так что при х->----оо эти возму­

щения обращаются в нуль. Учитывая уравнение гидроста­ тики из (3.4.1) и (3.4.2), пренебрегая величинами второго порядка малости, получим:

154


 

dU'

 

dp’

p0u

dw'

dp'

gP

Ро^ dx

 

дх

дх

dz

 

d£_

>dsq

 

du'

dw'

(3.4.3)

u

 

= 0.

dx

W

dz

 

dx

dz

Считая в соответствии с вышесказанным, что изменения плотности связаны с йзменениями энтропии, но не давления, запишем:

Р =

Фо

ps •

(3.4.4)

dsn

 

 

 

Из уравнения (3.4.3) с учетом (3.4.4) получим уравнение для w'\

d2w' .

d2w |

1

dp0

dw' , klw' = 0.

(3.4.5)

dx2

dz2

p0

dz

dz

 

Здесь обозначено:

*g = _ I _ Y J L \ J*l .

Po^2 \ ds0 ) p dz

Учитывая, что относительные изменения плотности в океане порядка 10~2— 10~3, и сравнивая порядок членов в уравнении (3.4.5), найдем, что величина третьего члена в (3.4.5) по крайней мере на два порядка меньше остальных и им можно пренебречь1. Таким образом, с высокой сте­ пенью точности уравнение (3.4.5) можно переписать в виде:

d2w' .

d W i , 2 /

п

/о л c\

— — +

— — + k0w

= 0 .

(3.4.6)

дх2

ду2

 

 

Используя известные термодинамические соотношения с учетом уравнения гидростатики

I

Фо \

Т0

/

гф 0

\

ds0

\

dso 1 р

ср

\

dT

/ р

dz

перепишем выражение для kl:

dTa

Ро“2 L V дТ0 ) p dz

cp

- +

_g_

(

дро

 

dT0

 

\

дТ0

 

T0

dz

Po

р

Л- Л

(. Фо

 

>

0.

(3.4.7)

Pocp V dT

/ p.

 

 

 

Здесь

70(г )— равновесное

значение

абсолютной

темпера­

туры:

изменения

солености морской

воды не принимаются

во внимание;

ср — теплоемкость при

постоянном давлении.

В

качестве

граничного

условия для w' на дне имеет

обычное кинематическое

соотношение

w = и ----,

линеа-

 

 

 

 

 

 

dx

 

1 Заметим, что это приближение опять-таки

в точности

соответст­

вует приближению

Буссинеска,

неоднократно

использовавшемуся выше.

155


ризируя которое и принимая во внимание малость £ по сравнению с Н — глубиной океана, получим:

w' = и — при zzzH.

(3.4.8)

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

хю' = 0 при г = 0,

(3.4.9)

кроме того, w' (—002),

w' (+ 002) — ограничено.

(3.4.10)

Как известно, в основной толще вод океана, за исклю­ чением приповерхностного слоя главного термоклина тол­ щиной порядка 500—1000 м, температура меняется крайне медленно и это изменение можно считать линейным. Поэто­ му, полагая в (3.4.7) dT0/dz=const и заменяя р0 его средним

значением, получим, что kl = const.

Положим (отбрасывая в дальнейшем штрих у а/), что

w = - ^ 7 U ^ r + 4>(x,z),

(3.4.11)

лах

тогда, учитывая (3.4.6) и неоднородное граничное условие (3.4.8), получим для <р неоднородное уравнение:

(Э2ф

Э2ср

&оф = ■ U ■

d%

+ k20 J L

(3.4.12)

5л:2

 

Н

\ dx3

 

dx

 

с граничными условиями:

 

 

 

 

Ф (х, 0) = ф (х, Н) = ф(— оо, z) — 0,

 

(3.4.13)

 

 

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

 

Удовлетворяя граничным условиям по z, положим

 

 

 

О

 

 

 

 

 

ф =

5 ] Ф» (*) sin Y„z

 

 

 

(3.4.14)

 

 

7 1 = 1

 

 

 

 

Подставляя

(3.4.14) в (3.4.12)

и разлагая

правую часть

в ряд по sinyn^, получим для ф„(х) уравнение:

 

 

42<Р/г

+ (kt— yl) Ф„ = исп

+ kl

 

(3.4.15)

dx3

 

 

dx3

'dx

 

156


где

2

(3.4.16)

с„ = (— 1)" пп

Интеграл этого уравнения, удовлетворяющий

условиям

(3.4.13) по х, будет

X

Фя = р р = = 2 [ Sin [К&о — у2(х — £)] х

(3.4.17)

при kl — у«> 0,

при у» —&о > 0■

(3.4.18)

Пусть теперь

 

 

 

| 0

при

—оо < х < 0

(3.4.19)

( g

sinp*

при 0 < х < [о о .

 

Как покажут приведенные

ниже численные

оценки для

интересующих нас значений параметров, слагаемыми, кото­

рые

даст формула (3.4.18),

молено в

первом приближении

пренебречь,

так как они образуют * знакопеременный ряд,

члены которого убывают во всяком

случае не

медленнее,

чем

1In.

Решения

вида

(3.4.18)

появляются,

когда

Упко >0,

т. е. (как

будет

видно ниже)

при

достаточно

больших п.

Поэтому,

учитывая появление

сильно

убываю­

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

Подставляя (3.4.19) в (3.4.17), выполняя квадратуры и простые тригонометрические преобразования, получим

Учитывая (3.4.11). (3.4.14), (3.4.16), (3.4.19), получим

(3.4.21)

157


Причем N соответствует последнему значению п, при кото­

ром все еще выполняется неравенство

kt — у« >0.

Уместно написать еще одно решение для случая, когда

рельеф дна задан функцией:

 

 

 

 

 

 

 

 

 

 

0

при

 

х < 0

 

 

 

 

£0/2 при

х —0

(3.4.22)

 

 

 

1£

при

 

х > 0

 

Тогда, очевидно, £' (х) = £8),

где б(х)

— функция Ди­

рака. Используя основное свойство б-функции:

ъ

 

 

 

 

 

 

 

 

 

J / (X -

 

1) 6<л) dl = ( -

1)"/<"> (х),

а < х < Ь,

а также формулы (3.4.11), (3.4.17), (3.4.18), получим:

w = ы£0

— б (х) + У

- ^ п-

sin V

kl — ylx х

 

Н

n=i

V k 2 v2

 

 

 

 

 

г к0

Yo

 

 

X sin

-

 

c«Yre

■sin

 

Y„xsm y„2.

 

 

Yo

 

 

 

 

л=ЛГ+1

 

 

 

 

(3.4.23)

 

 

 

 

 

 

 

 

 

Эта формула имеет особенность при х=0, что объяснимо физически, так как в окрестности этой точки линейная тео­ рия неприменима. Первая сумма в (3.4.23) должна вычис- I

ляться д л ях > 0,

сходимость

бесконечной

суммы в

(3.4.23)

при |х |=^=0

обеспечивается

наличием

экспоненциальных

множителей.

Для

х > 0 можем вместо (3.4.23) написать:

 

N

с пУ2п

•sin Vm- Y«xsinY„2.

 

W =

 

(3.4.24)

 

 

П~\ V k i ­ l l

Теперь проведем численные оценки, основываясь на ре­ шении (3.4.21), достаточно близко соответствующем реаль­

ным условиям.

___ Отг

Пусть X = — = 2 км, т. е. 0 = 3,14-10—5 см-1. Как из-

Р

вестно, за исключением верхнего слоя главного термоклина (500—1000 м), температура в океане изменяется с глубиной очень слабо. Поэтому для оценки положим dT0fdz и возьмем согласно (3.4.7) минимальное значение k0:

kl = —-— (

) л Г И .

Рйи 2 \

дТ0 J p V ср

158