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