Файл: Тареев, Б. А. Динамика бароклинных возмущений в океане.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 21.10.2024
Просмотров: 133
Скачиваний: 0
u = ---------dp/dy, |
v = ------ dp/dx |
в уравнение вихря. Тогда |
|||||
р/ |
и |
|
Р/ |
получим уравнение для р: |
|||
из (2.9.5), (2.9.6) |
(2.9.3) |
||||||
|
|
|
f i f - U y ) |
'dp' |
|
|
|
|
|
|
g(H+l) |
dt |
|
|
|
\ U (f — Uy) f |
r r |
(f~Uy) |
|
(2.9.7) |
|||
gH |
|
|
yy |
H |
|
|
|
Если включить g |
в определение Я с учетом. (2.9.3) и поло |
||||||
жить p — p(y)expik(x—ct), |
уравнение для |
амплитуды р(у) |
|||||
примет вид: |
|
|
|
|
|
|
|
W ( y ) - c ] |
d2p |
k2p — 1М.:иу> |
p |
+ |
|||
|
|||||||
dy3 |
|
gH н |
|||||
+ |
[ |
-(t/g |
-f)- — |
— Um]p = 0 |
(2.9.8) |
||
|
н |
|
dy |
yyY |
|
|
|
(здесь опущена черточка над р). На границах |
течения у\ и |
г/г нормальная компонента скорости должна обращаться в
нуль, что в квазигеострофическом приближении |
дает: |
|||
|
Р Ш = Р Ш = 0- |
|
(2.9.9) |
|
Если вне |
рассматриваемой |
области |
(z/b уг) |
£/=const и |
Я = const, то, |
как следует из |
(2.9.8), за |
пределами \у\, у2) |
амплитуда волновых возмущений экспоненциально затухает: р = ех р (± % ) где k = ]/& а + &о , ko=--f2/gH, а знак выби
рается, чтобы обеспечить затухание при у-^+оо. Это приводит к другому виду граничных условий: •
— ± kp при у — ух или у — уг. |
(2.9.10) |
Границы могут находиться и на бесконечности. Уравнение (2.9.8) при условиях (2.9.9) или (2.9.10) может, вообще го воря, иметь решения с комплексными собственными значения ми C= Cr-\-iCi, которые приводят к неустойчивости. Условия, необходимые для неустойчивости, легко, доказать, следуя обычной методике. Обозначая выражение в последних квад ратных скобках (2.9.8) через К (у), полагая p — pr+ipi,
(U—c) _1 = ai+ ta2, где а \= (Я—с2) | U—с\~2; a2=Ci\U—c\-2,
отделяя вещественную часть (2.9.8) от мнимой, получим
d2Pr |
И2 + |
n f - V y ) |
__aiK{y) |
pr— a1K{y)pi = 0 , |
|
[dif |
|
gH |
|
|
|
d2Pi |
k2 |
f i f - U y ) |
UiK (y) |
Pi ~j- a%K (y) Pr = 0. |
|
dy2 |
|||||
|
gH |
|
|
127
Отсюда следует (после умножения первого уравнения на * Pi, второго на рти вычитания):
j а $ ( у ) \ p f d y - „ : 0. |
■ (2.9.11) |
у. |
|
Интегрируя |
(2.9.11) |
в пределах у\, у2, при любом виде гра |
|
ничных условий >(2.9.9) или (2.9.10) получим: |
|||
- j —(Pi |
dy |
- |
Рг- ^ ~ ) = - а 2К(у) (рг2+ рЬ- (2.9.12) |
dy \ |
|
dy ] |
Следовательно, если с^Ф0 (неустойчивость), суммарный гра-
диент завихренности К(у) = U - U y ) |
дН ■— и, |
должен ме |
|
Н |
ду |
уу |
нять знак внутри области. Полученный простой критерий так называемой баротропной неустойчивости (т. е. неустойчиво сти, источником которой является кинетическая энергия не возмущенного движения) может иметь интересные приложе ния к анализу баротропной компоненты течений, проходящих вблизи материкового склона (например, начальный участок Гольфстрима, Г. Стоммел (1963)).
На материковом склоне А Н /Н ~\ и при U~ 1 м/еек, Дг/~100 км (оценки, характерные для Гольфстрима) первый член в выражении для К {у) будет на порядок больше кривиз ны поперечного профиля скорости.
Отсюда следует, что расчеты баротропной неустойчивости Гольфстрима, основанные на бездивергентных моделях (Гаурвитц и Пановски, 11950), недостаточны для определения фак тического характера движения.
Конечно, само по себе изменение глубины без поперечного
изменения скорости течения не может служить причиной неус тойчивости. В этом смысле полученный критерий не являет ся достаточным. В частности, при U = const, после умножения (2.9.8) на р— рг—ipi, интегрирования по частям и использо вании граничного условия (2.9.9) получим:
(U— с) J j |p '|2 + (& + |
|Р I2} d y = j К(у)\р I2dy. |
У, |
У1 |
Откуда |
|
у2 |
У2 |
О | [\Р' |2+ |
\p\2\}dy= ^ К(у)\ p\2dy, (2.9.13) |
Ul |
|
128
Уг
( U - c r) == |
_____ |
i K ( y ) \P |2 dy |
(2.9.14) |
y\_____________________ |
|
||
|
У2 |
|
|
j {lpl3+ (Й3+ _£г)|р|2}^
Hi
Следовательно, если U = const, то с*=0 и течение устойчиво. Это очевидно также из энергетических соображений, так как для нарастания волновых возмущений, как известно, необхо-
димо |
|
Уг |
%Uydy^>0, где т = —рш>. |
|
выполнение неравенства j |
||||
|
|
!Л |
ci=0, |
соотношение (2.9.14) |
Когда Я = const и, следовательно, |
||||
дает удобную |
возможность для приближенного вычисления |
|||
дисперсионных соотношений вариационным методом. |
||||
С |
другой |
стороны, из (2.9.8) |
при U—const следует, что |
|
если |
Я (у) — симметричная функция в |
интервале \у\, г/г) |
(симметричный относительно гребня подводный хребет или
ложбина), то |
К (у) = |
----- — дН/ду |
— антисимметрична. Сле- |
|
|
|
н |
|
и антисим |
довательно, решениями (8) будут симметричные |
||||
метричные в |
(уи г/2) |
функции: р(у) — ± р {—у). |
Тогда из |
|
(2.9.14) ввиду антисимметричности |
К(у) следует U—с=0, |
т. е. волны, распространяющиеся вдоль симметричных под водных хребтов, сводятся к установившимся отклонениям, движущимся со скоростью среды.
На материковом склоне условия существенно несиммет ричны, поэтому целесообразно подробнее рассмотреть пример достаточно простого, но характерного изменения глубины.
Когда Я меняется линейно, глубина по мере отдаления от берега стремится к бесконечности. В действительности, ма териковый склон имеет конечную ширину, поэтому ниже будет рассмотрен достаточно типичный случай
Я = Я0(1 —е*У). |
(2.9.15) |
Начало координат возьмем на береговом крае материко вого склона. Как видно из (2.9.15), в открытом океане (г/-*— оо) Я = Я 0. Для U—const (и, следовательно, веще ственных с) уравнение (2.9.8) перепишем в виде:
d?P , |
f._______ дН_ |
р |
= 0 |
(2.9.16) |
||
dy2 |
(с - U) Н |
ду |
||||
|
|
|
||||
или, обозначая |
%=y.f{U—с)-1, a2= f 2/gH, -из |
(2.9.15) и |
||||
(2.9.16) получим: |
|
|
|
|
|
|
d2p |
К |
а2 |
Р = |
0. |
(2.9.17) |
|
~dy* |
1- е кУ |
|||||
1-<*» |
|
|
|
9 Б. А. Тареев |
129 |
Поскольку нормальная к береговой черте компонента ско рости должна обращаться в нуль, граничные условия для
(2.9.17): р = 0 при у = 0, р — ограничено при г/->— оо (2.9.18).
Если ехр(—ху)_— х, то вместо (2.9.17) получим
хЦх — 1) ^ Е - + х ( х — 1) -£ -+ [* ,' — a'2x) — k'2(x— \)]p=0.
|
|
|
|
|
|
|
|
|
|
(2.9.19) |
Здесь а’2= а2х~2, k'2x~ 2, |
к' — кх~2. |
Граничные |
условия для |
|||||||
(2.9.19): |
|
|
|
|
|
|
|
|
|
|
|
|
р (1) = 0, |
р (оо) — ограничено. |
|
(2.9.20) |
|||||
При р = хтф (х), |
т = ± V k ’2+ к’ уравнение (2.9.19) |
приво |
||||||||
дится к виду: |
|
|
|
|
|
|
|
|
|
|
|
х(х — 1) ср" + |
(х — 1) (2т + 1) ф' + (Г — а'2) ф = 0.- |
(2.9.21) |
|||||||
Определяя параметры а, |
(3, у равенствами |
|
|
|||||||
|
a = m ± V k’2+ |
а'2, |
р = т ± У k’2+ а'2, |
|
||||||
|
у = |
2т + 1, (т = ± V k'2 + х'2 , |
|
(2.9.22) |
||||||
приводим уравнение (2.9.21) к гипергеометрическому:] |
|
|||||||||
|
х ( х — 1)ф" + [а -Ь Р + |
1) л: — у)] ф' + <*Рф = 0. |
(2.9.23) |
|||||||
Общее решение (2.9.23) в рассматриваемом интервале: |
|
|||||||||
|
Ф = хга F (а, |
1 + а — у, |
1 + а — р; г -1) + |
|
||||||
|
ч- х~&F (Р, |
1-гР — У, |
1-гР —а; х~]). |
(2.9.24) |
||||||
Здесь F ^ Л, В , С; — j — гипергеометрический ряд, |
1 < | х | <оо. |
|||||||||
Из |
(2.9.24) получаем: |
|
|
|
|
|
|
|
||
|
р = Clxm~aF ( А, |
В, С; — ) + с2хт~Р F (Л', В', |
С'; — ). |
|||||||
|
|
|
|
|
|
|
|
|
|
(2.9.25) |
Но |
независимо |
от |
знака |
т, |
если в |
выражениях 2.9.22 для |
||||
а, р выбрать верхние знаки, получим |
|
|
|
т — а — — V k '2 + а'2 < 0 и т — р = V k'2 + а’2 > 0 ,
поэтому условие ограниченности р (оо) приводит к С2= 0 , и решение для р с учетом использовавшихся выше обозначений:
130