Генерация поверхностных и внутренних волн в ограниченном бассейне перемещающимся барическим фронтом
Рассмотрена плоская задача о генерации поверхностных и внутренних волн в ограниченном вращающемся бассейне переменной глубины при прохождении над ним фронта атмосферного давления. Жидкость предполагалась двухслойной. Система нелинейных уравнений длинных волн решалась численно методом конечных разнос...
Gespeichert in:
| Veröffentlicht in: | Морской гидрофизический журнал |
|---|---|
| Datum: | 2009 |
| Hauptverfasser: | , |
| Format: | Artikel |
| Sprache: | Russian |
| Veröffentlicht: |
Морський гідрофізичний інститут НАН України
2009
|
| Schlagworte: | |
| Online Zugang: | https://nasplib.isofts.kiev.ua/handle/123456789/105051 |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| Zitieren: | Генерация поверхностных и внутренних волн в ограниченном бассейне перемещающимся барическим фронтом / С.Ф. Доценко, Н.А. Миклашевская // Морской гидрофизический журнал. — 2009. — № 3. — С. 3-18. — Бібліогр.: 27 назв. — рос. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraine| id |
nasplib_isofts_kiev_ua-123456789-105051 |
|---|---|
| record_format |
dspace |
| spelling |
Доценко, С.Ф. Миклашевская, Н.А. 2016-08-05T18:38:29Z 2016-08-05T18:38:29Z 2009 Генерация поверхностных и внутренних волн в ограниченном бассейне перемещающимся барическим фронтом / С.Ф. Доценко, Н.А. Миклашевская // Морской гидрофизический журнал. — 2009. — № 3. — С. 3-18. — Бібліогр.: 27 назв. — рос. 0233-7584 https://nasplib.isofts.kiev.ua/handle/123456789/105051 551.466.86 Рассмотрена плоская задача о генерации поверхностных и внутренних волн в ограниченном вращающемся бассейне переменной глубины при прохождении над ним фронта атмосферного давления. Жидкость предполагалась двухслойной. Система нелинейных уравнений длинных волн решалась численно методом конечных разностей для распределения глубины, соответствующего зональному сечению Черноморской котловины. Показано, что при прохождении над бассейном барического фронта генерируются баротропные и бароклинные колебания жидкости. Интенсивность возмущений зависит от скорости переноса и ширины фронта. Существует скорость движения фронта, при которой генерация внутренних волн происходит наиболее эффективно. После выхода фронта за пределы водоема образуется пакет внутренних волн, распространяющийся с отражением от одной боковой границы бассейна к другой. Plane problem on generation of surface and internal waves in a bounded rotating basin of variale depth by an atmospheric pressure front moving above the basin is considered. The fluid is assumed two-layer. The nonlinear system of long-wave equations is solved numerically by the method of finite differences for the depth distribution corresponding to the zonal section of the Black Sea basin. It is shown that barotropic and baroclinic fluid oscillations are generated as a result of a baric front translation above the basin. Intensity of fluid disturbances depends on translation velocity and the front width. There is a velocity of front translation at which the internal waves are generated most effectively. After the front leaves the basin a packet of internal waves which propagates reflecting from one lateral boundary to the other is formed. ru Морський гідрофізичний інститут НАН України Морской гидрофизический журнал Термогидродинамика океана Генерация поверхностных и внутренних волн в ограниченном бассейне перемещающимся барическим фронтом Article published earlier |
| institution |
Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| collection |
DSpace DC |
| title |
Генерация поверхностных и внутренних волн в ограниченном бассейне перемещающимся барическим фронтом |
| spellingShingle |
Генерация поверхностных и внутренних волн в ограниченном бассейне перемещающимся барическим фронтом Доценко, С.Ф. Миклашевская, Н.А. Термогидродинамика океана |
| title_short |
Генерация поверхностных и внутренних волн в ограниченном бассейне перемещающимся барическим фронтом |
| title_full |
Генерация поверхностных и внутренних волн в ограниченном бассейне перемещающимся барическим фронтом |
| title_fullStr |
Генерация поверхностных и внутренних волн в ограниченном бассейне перемещающимся барическим фронтом |
| title_full_unstemmed |
Генерация поверхностных и внутренних волн в ограниченном бассейне перемещающимся барическим фронтом |
| title_sort |
генерация поверхностных и внутренних волн в ограниченном бассейне перемещающимся барическим фронтом |
| author |
Доценко, С.Ф. Миклашевская, Н.А. |
| author_facet |
Доценко, С.Ф. Миклашевская, Н.А. |
| topic |
Термогидродинамика океана |
| topic_facet |
Термогидродинамика океана |
| publishDate |
2009 |
| language |
Russian |
| container_title |
Морской гидрофизический журнал |
| publisher |
Морський гідрофізичний інститут НАН України |
| format |
Article |
| description |
Рассмотрена плоская задача о генерации поверхностных и внутренних волн в ограниченном вращающемся бассейне переменной глубины при прохождении над ним фронта атмосферного давления. Жидкость предполагалась двухслойной. Система нелинейных уравнений длинных волн решалась численно методом конечных разностей для распределения глубины, соответствующего зональному сечению Черноморской котловины. Показано, что при прохождении над бассейном барического фронта генерируются баротропные и бароклинные колебания жидкости. Интенсивность возмущений зависит от скорости переноса и ширины фронта. Существует скорость движения фронта, при которой генерация внутренних волн происходит наиболее эффективно. После выхода фронта за пределы водоема образуется пакет внутренних волн, распространяющийся с отражением от одной боковой границы бассейна к другой.
Plane problem on generation of surface and internal waves in a bounded rotating basin of variale depth by an atmospheric pressure front moving above the basin is considered. The fluid is assumed two-layer. The nonlinear system of long-wave equations is solved numerically by the method of finite differences for the depth distribution corresponding to the zonal section of the Black Sea basin. It is shown that barotropic and baroclinic fluid oscillations are generated as a result of a baric front translation above the basin. Intensity of fluid disturbances depends on translation velocity and the front width. There is a velocity of front translation at which the internal waves are generated most effectively. After the front leaves the basin a packet of internal waves which propagates reflecting from one lateral boundary to the other is formed.
|
| issn |
0233-7584 |
| url |
https://nasplib.isofts.kiev.ua/handle/123456789/105051 |
| citation_txt |
Генерация поверхностных и внутренних волн в ограниченном бассейне перемещающимся барическим фронтом / С.Ф. Доценко, Н.А. Миклашевская // Морской гидрофизический журнал. — 2009. — № 3. — С. 3-18. — Бібліогр.: 27 назв. — рос. |
| work_keys_str_mv |
AT docenkosf generaciâpoverhnostnyhivnutrennihvolnvograničennombasseineperemeŝaûŝimsâbaričeskimfrontom AT miklaševskaâna generaciâpoverhnostnyhivnutrennihvolnvograničennombasseineperemeŝaûŝimsâbaričeskimfrontom |
| first_indexed |
2025-11-25T21:07:22Z |
| last_indexed |
2025-11-25T21:07:22Z |
| _version_ |
1850550070485712896 |
| fulltext |
ISSN 0233-7584. Мор. гидрофиз. журн., 2009, № 3
3
Термогидродинамика океана
УДК 551.466.86
С.Ф. Доценко, Н.А. Миклашевская
Генерация поверхностных и внутренних волн
в ограниченном бассейне перемещающимся
барическим фронтом
Рассмотрена плоская задача о генерации поверхностных и внутренних волн в ограничен-
ном вращающемся бассейне переменной глубины при прохождении над ним фронта атмо-
сферного давления. Жидкость предполагалась двухслойной. Система нелинейных уравнений
длинных волн решалась численно методом конечных разностей для распределения глубины,
соответствующего зональному сечению Черноморской котловины. Показано, что при прохож-
дении над бассейном барического фронта генерируются баротропные и бароклинные колеба-
ния жидкости. Интенсивность возмущений зависит от скорости переноса и ширины фронта.
Существует скорость движения фронта, при которой генерация внутренних волн происходит
наиболее эффективно. После выхода фронта за пределы водоема образуется пакет внутренних
волн, распространяющийся с отражением от одной боковой границы бассейна к другой.
Сейши представляют собой свободные стоячие волны в замкнутых или
полузамкнутых бассейнах [1 – 3]. Они могут вызываться различными внеш-
ними воздействиями, в число которых входят [4]: неравномерное распреде-
ление атмосферного давления над бассейном и его изменение со временем;
действие ветровых напряжений на поверхность бассейна; резкие повышения
или понижения уровня жидкости в какой-либо части водоема за счет притока
или оттока вод, выпадения дождя или града; сейсмотектонические подвижки
дна, если зона деформаций земной коры захватывает участок дна бассейна
[5]; внутренние волны, в частности приливные, отражающиеся от кромки
шельфа [6].
Сейши вносят заметный вклад в изменчивость гидродинамических по-
лей, действуя продолжительное время, достигающее 80 – 90 % безледостав-
ного периода в таких озерах как Байкал и Балхаш [4]. Сейшевые колебания
неизменно сопровождаются течениями, измерение которых позволяет иден-
тифицировать сейши и изучать их свойства, как это сделано, например, в ра-
ботах [4, 7]. Сейши влияют на перенос вещества и перемешивание в морской
среде, на перераспределение химических и биологических веществ. Значи-
тельная роль в этом отводится плотностной стратификации. Поэтому теоре-
тическое и экспериментальное изучение сейш в озерах, фиордах, морских и
океанских бассейнах с учетом бароклинности среды – важное направление
современных океанологических исследований.
Горизонтальная структура и параметры длинных внутренних волн в
замкнутых бассейнах (бароклинных сейш) зависят от геометрии бассейна и
плотностной (температурной) стратификации. Хотя первое теоретико-
С.Ф. Доценко, Н.А. Миклашевская, 2009
ISSN 0233-7584. Мор. гидрофиз. журн., 2009, № 3
4
экспериментальное исследование структуры внутренних сейш и связанных с
ними течений в озерах было выполнено около ста лет назад [8], внутренние
сейши изучены значительно менее детально, чем баротропные. Во многих
теоретических работах анализ бароклинных сейш проводился в рамках двух-
слойных, иногда трехслойных, моделей морской среды. Обзор исследований
внутренних сейш содержится в работе [9].
Генерация низших мод баротропных сейш при перемещении зон атмо-
сферных возмущений над ограниченными и полуограниченными бассейнами
подтверждена как данными наблюдений [10 – 14], так и результатами чис-
ленного моделирования в рамках теории длинных волн [15 – 19]. Возмож-
ность возбуждения внутренних сейш при атмосферных воздействиях на вер-
тикально стратифицированную жидкость в ограниченных и полуограничен-
ных бассейнах подтверждают результаты инструментальных измерений ко-
лебаний температуры и скорости течения во внутренних слоях бассейна [20 –
22] и колебаний самой поверхности моря [23]. Вклад внутренних сейш в из-
менчивость полей скорости и температуры в бассейне нередко является до-
минирующим [22].
Математические модели, ориентированные на анализ генерации внут-
ренних сейш атмосферными воздействиями в бассейнах переменной глуби-
ны, разработаны недостаточно. В равной мере это относится и к моделирова-
нию процессов генерации бароклинных сейш движущимися барическими
фронтами. Ниже исследуется плоская задача о возбуждении поверхностных и
внутренних сейш при прохождении фронта атмосферного давления над огра-
ниченным бассейном переменной глубины.
Математическая модель. Рассматривается плоская задача о генерации
волн в бассейне, занимающем область 0 ≤ x ≤ l, –H(x) < z < 0 (рис. 1), где x –
горизонтальная, z – вертикальная координата, отсчитываемая вверх от невоз-
мущенного положения свободной поверхности жидкости z = 0; l – ширина
бассейна; H = H(x) > 0 – переменная глубина бассейна в невозмущенном со-
стоянии. Вторая горизонтальная ось y направлена перпендикулярно плоско-
сти Oxz. Бассейн заполнен двухслойной жидкостью. В невозмущенном со-
стоянии толщина верхнего слоя равна h1 и постоянна, плотность жидкости
равна 1. Толщина нижнего слоя h2(x) = H(x) – h1 зависит от горизонтальной
координаты x, плотность 2 > 1.
Р и с. 1. Схема задачи
ISSN 0233-7584. Мор. гидрофиз. журн., 2009, № 3
5
Генератором волн является перемещающийся в положительном направ-
лении оси х с постоянной скоростью С над поверхностью бассейна фронт ат-
мосферного давления (рис. 1). Распределение атмосферного давления попе-
рек фронта задается в виде
CtxgFppp ),(10 ,
где p0 > 0 – среднее атмосферное давление, предполагаемое постоянным; p –
амплитуда возмущений атмосферного давления (в метрах водяного столба); g
– ускорение свободного падения; F() – безразмерная функция, удовлетво-
ряющая следующим условиям на бесконечности: F() → 0 при → + ∞,
F() → 1 при → – ∞. Таким образом, под барическим фронтом понимается
область гладкого перехода от среднего значения p0 к некоторому другому
постоянному значению p0 + p1g атмосферного давления. Характерная про-
тяженность зоны существенного горизонтального изменения атмосферного
давления равна W и выступает в роли ширины барического фронта.
Гидростатическое смещение свободной поверхности жидкости, вызван-
ное действием атмосферного давления, определяется по формуле, известной
в океанологии как приближение обратного барометра [1, 24]:
)(
1
0
pF
g
pp
z
. (1)
В рамках нелинейной теории длинных волн, предполагающей распреде-
ление давления в жидкости гидростатическим, движение двухслойной жид-
кости в бассейне переменной глубины с учетом вращения Земли описывается
в вертикальной плоскости Oxz системой из шести уравнений:
x
gfv
x
u
u
t
u
)
~
( 1
1
1
1
1
, 01
1
1
1
fu
x
v
u
t
v
, (2)
0
])[( 121121
x
uh
t
; (3)
x
g
x
gfv
x
u
u
t
u
21
2
2
2
2 )
~
(
, 02
2
2
2
fu
x
v
u
t
v
, (4)
0
])[( 2222
x
uh
t
(5)
с начальными условиями
)0( 02,12,12,1 tvu , (6)
ISSN 0233-7584. Мор. гидрофиз. журн., 2009, № 3
6
означающими, что жидкость в бассейне в начальный момент времени была
невозмущена.
Здесь uj(x, t), vj(x, t) (j = 1, 2) – проекции на оси x и y соответственно ос-
редненной по глубине горизонтальной скорости течения в верхнем (j = 1) и
нижнем (j = 2) слоях; 1(x,t), 2(x,t) – смещения свободной поверхности жид-
кости и границы раздела слоев от невозмущенных горизонтальных положе-
ний соответственно; )(
~
CtxpF ; = 1/2 < 1 и = 1 – >0 – безразмер-
ные параметры, характеризующие вертикальную плотностную стратифика-
цию; f – постоянный параметр Кориолиса. Заметим, что в начальный момент
времени t = 0 фронтальная зона располагается вне поверхности бассейна.
На боковых границах бассейна x = 0 и x = l, являющихся вертикальными
твердыми стенками, задаются условия
u1,2(0, t) = 0, u1,2(l, t) = 0, (7)
означающие непротекание жидкости через боковые границы в каждом слое.
Численный алгоритм расчета колебаний неоднородной жидкости.
Для численного решения начально-краевой задачи (2) – (7) на отрезке
0 x l при t ≥ 0 использована явно-неявная конечно-разностная схема [25].
Скорости u1,2 в момент времени t = tn (n = 0, 1, …) вычислялись в узлах сетки
x = xi (i = 0, …,N) с равномерным шагом = l/N (xi = i), а y-проекции скоро-
стей течения в слоях v1,2 и смещения границ 1,2 – в средних точках x = xi – /2
(i = 1, …, N) отрезков [xi–1, xi].
Для удобства записи разностного аналога задачи (2) – (7) временно обо-
значим неизвестные гидродинамические поля следующим образом: 1,
u = u1, v = v1 и Z 2, U = u2, V = v2. Введем величины
ntn , 0 ),,}(,{, ,000 n
N
n
Nni
n
i
n
i
nn UutxUuUuUu ,
ni
n
i
n
i
n
i
n
i txZVvZVv ,2/},,,{,,, (i = 1, ..., N ),
где – постоянный шаг интегрирования по времени. Толщины слоев в невоз-
мущенном состоянии зададим в узлах расчетной сетки: h1i = h1(xi), h2i = h2(xi).
Разностные аппроксимации уравнений (2) по явной схеме имеют вид
n
i
n
i
n
i
n
i
n
i
n
i
n
i
n
i
n
i
n
i gvfuuuuu
~~
5,0 112/111
1
, (8)
n
i
n
i
n
i
n
i
n
i
n
i ufvvuvv 2/1112/1
1 5,0
, (9)
где )(5,0 12/1
n
i
n
i
n
i uuu , )(5,0 12/1
n
i
n
i
n
i vvv , = /, i = 1, …, N – 1.
Разностный аналог уравнения неразрывности (3) для верхнего слоя запи-
сывается в неявном виде:
ISSN 0233-7584. Мор. гидрофиз. журн., 2009, № 3
7
n
i
n
i
n
i
n
i
n
i
n
i
n
i
n
i DuDuZZ 11
1
11
111
(i = 1, …, N), (10)
где )(5,0 1111
n
i
n
i
n
i
n
ii
n
i ZZhD – динамическая толщина верхнего слоя
в узле сетки x = xi.
Разностные аппроксимации уравнений (4) по явной схеме имеют вид
, ~~
5,0
111
2/111
1
n
i
n
i
n
i
n
i
n
i
n
i
n
i
n
i
n
i
n
i
n
i
n
i
ZZgg
VfUUUUU
(11)
n
i
n
i
n
i
n
i
n
i
n
i UfVVUVV 2/1112/1
1 5,0
, (12)
здесь )(5,0 12/1
n
i
n
i
n
i UUU , )(5,0 12/1
n
i
n
i
n
i VVV , i = 1,…, N – 1.
Уравнение неразрывности (5) для нижнего слоя в разностном виде запи-
сывается следующим неявным образом:
n
i
n
i
n
i
n
i
n
i
n
i DUDUZZ 12
1
12
11
(i = 1, …, N), (13)
где )(5,0 122
n
i
n
ii
n
i ZZhD – динамическая толщина нижнего слоя при x = xi.
В формулах (10) и (13) выполнена квазилинеаризация: использованы
значения горизонтальных скоростей на временном шаге n + 1, а смещения
свободной поверхности и границы раздела слоев – на предыдущем шаге по
времени. Аналог начальных условий (7) для системы разностных уравнений
(8) – (13) записывается в виде
000 ii Uu (i = 0, ..., N), 00000 iiii ZVv (i = 1, ..., N). (14)
Таким образом, исходная дифференциальная задача (2) – (7) заменена
конечно-разностной (8) – (14). Для численного определения гидродинамиче-
ских полей u1,2, v1,2 и 1,2 необходимо, начиная с n = 0, выполнять при перехо-
де от временного слоя n к временному слою n + 1 такую последовательность
вычислений полей по переменной x: (14) [(8), (9), (11), (12)] (13) (10).
Результаты численного анализа. Численный анализ процесса генера-
ции поверхностных и внутренних волн движущимся фронтом проводился для
распределения глубины бассейна H(x), показанного на рис. 2. Рельеф дна со-
ответствует зональному сечению Черного моря 42,66 с.ш. Ширина бассейна
l = 894 км. В вычислительных экспериментах задавались толщина верхнего
слоя h1 = 50 м и параметры плотностной стратификации 1 = 1012 кгм
–3
,
2 = 1017 кгм
–3
. В этом случае ≈ 0,0049. В численных экспериментах скорость
переноса барического фронта изменялась в интервале 0,5 мc
–1
< C < 10 мc
–1
,
ширина фронта – в диапазоне 50 км < W < 500 км, а перепад атмосферного
давления при пересечении фронта был задан равным p = 0,3 м.
ISSN 0233-7584. Мор. гидрофиз. журн., 2009, № 3
8
Р и с. 2. Распределение глубины бассейна, соответствующее зональному сечению Черномор-
ской котловины (42,66 с. ш.)
Функция F() в (1), описывающая гладкое распределение атмосферного
давления поперек фронта, задавалась по формуле
)0( 0 ),( 1 ),0 ( )2/(sin 2 FWFWWF . (15)
В двухслойной жидкости постоянной глубины могут существовать два
качественно различных типа свободных колебаний, описанных, например, в
работе [26]. Первый класс колебаний обусловлен смещениями свободной по-
верхности жидкости, и влияние на него неоднородности среды мало. Гори-
зонтальные скорости в слоях вдоль любой из осей х или y направлены в одну
сторону и близки по величине. Колебания свободной поверхности и границы
раздела слоев синфазны, отношение амплитуд колебаний свободной поверх-
ности и границы раздела слоев 21211 // hHAA . Эта мода колебаний является
баротропной, скорость распространения соответствующих ей волн
gHc 1 . Второй тип волн обусловлен неоднородностью жидкости (барок-
линная мода). Для него амплитуда колебаний границы раздела слоев значи-
тельно больше амплитуды соответствующих колебаний свободной поверхно-
сти жидкости, а именно, HhAA // 22221 . Горизонтальные скорости в слоях
направлены противоположно, различны по абсолютной величине, но полные
потоки жидкости в слоях близки по абсолютной величине, т. е. в линейном
приближении 02211 huhu и 02211 hvhv . Колебания свободной поверх-
ности и границы раздела слоев противофазны. Скорость распространения
внутренних волн в двухслойной жидкости Hhhgc /212 .
На рис. 3 показана форма свободной поверхности и границы раздела сло-
ев при различных положениях барического фронта над бассейном. Смещения
свободной поверхности и скачка плотности повторяют по форме распределе-
ние атмосферного давления, перемещаются вместе с фронтом и практически
совпадают. При скорости C = 1 мс
–1
отклик двухслойной жидкости на внеш-
нее воздействие является баротропным. Смещения свободной поверхности и
ISSN 0233-7584. Мор. гидрофиз. журн., 2009, № 3
9
скачка плотности существенно отличаются от гидростатического (1), что свя-
зано с ограниченностью бассейна и требованием сохранения массы каждого
из слоев:
l
dxtx
0
1 0),( ,
l
dxtx
0
2 0),( .
Р и с. 3. Вертикальные смещения поверхности жидкости (штриховая линия) и границы раздела
слоев (сплошная линия) при различных положениях xф = Ct передней границы барического
фронта: а – xф = l/3; б – xф = 2l/3; в – xф = l (ширина фронта W = 100 км, скорость его движения
C = 1 мс–1)
ISSN 0233-7584. Мор. гидрофиз. журн., 2009, № 3
10
При увеличении скорости переноса атмосферного фронта над поверхно-
стью бассейна динамический отклик двухслойной жидкости в бассейне пре-
терпевает качественные изменения. Такой случай показан на рис. 4.
Р и с. 4. Вертикальные смещения поверхности жидкости (штриховая линия) и границы раздела
слоев (сплошная линия) при различных положениях xф = Ct передней границы барического
фронта: а – xф = l/3; б – xф = 2l/3; в – xф = l; г – xф = l + W (ширина фронта W = 100 км, скорость
его переноса C = 3 мс–1)
Очевидно, что смещения свободной поверхности и скачка плотности
формируются из баротропной и бароклинной мод. Смещения свободной по-
верхности и границы раздела слоев близки по форме благодаря определяю-
щему вкладу баротропной моды. Бароклинность проявляется главным обра-
зом в волновых возмущениях скачка плотности. Эти возмущения представ-
ляют собой волновой след непосредственно за движущимся фронтом шири-
ны W изменения атмосферного давления от одного значения к другому
ISSN 0233-7584. Мор. гидрофиз. журн., 2009, № 3
11
(см. рис. 1). В процессе пересечения фронтом бассейна область, охваченная
бароклинными колебаниями, расширяется. После выхода фронта за пределы
бассейна (рис. 4, г) в нем продолжаются остаточные (бароклинные) колеба-
ния жидкости.
Скорость распространения внутренних волн в центральной части бассей-
на, где его глубина постоянна, с2 = 1,53 мc
–1
. Возмущения двухслойной жид-
кости, показанные на рис. 3 и 4, соответствуют перемещениям атмосферного
фронта со скоростями, меньшими (рис. 3) и большими (рис. 4) значения с2.
Изменение характера бароклинных движений жидкости при увеличении ско-
рости переноса барической аномалии можно объяснить переходом от докри-
тических значений скорости переноса фронта C < c2 к сверхкритическим зна-
чениям C > c2, когда в длинноволновом приближении генерируется незату-
хающий при удалении от генератора волн (зоны ненулевых градиентов атмо-
сферного давления ширины W, рис. 1) бароклинный волновой след, соответ-
ствующий низшей моде внутренних волн [27].
При выходе барического фронта через правую боковую границу за пре-
делы бассейна происходит отражение системы внутреннего волнового следа
от правой границы бассейна. В результате этого образуется пакет внутренних
волн, распространяющийся от правой к левой границе бассейна, а после от-
ражения от левой границы – в обратном направлении. Характер бароклинных
колебаний жидкости в бассейне иллюстрирует рис. 5.
Р и с. 5. Колебания границы раздела слоев в различных точках бассейна x при движении со
скоростью С = 2 м·с–1 барического фронта шириной W = 50 км: а – x = 0; б – x = l/2; в – x = l
ISSN 0233-7584. Мор. гидрофиз. журн., 2009, № 3
12
Время пересечения бассейна фронтом оценивается как (l + W)/C = 264 ч.
Такие внутренние волны вряд ли можно отнести к традиционным бароклин-
ным сейшам, поскольку они не сопровождаются синхронным движением
всей границы раздела слоев в бассейне, а многочисленные узловые линии,
характерные для сейш высокого порядка, совершают существенные переме-
щения в пространстве вместе с волновым пакетом.
Форма развитого пакета внутренних волн, возникающего при отражении
волнового следа от правой боковой границы при выходе барического фронта
за пределы бассейна, показана на рис. 5, в. После отражения пакет пересекает
бассейн и отражается от левой боковой границы, что иллюстрирует рис. 5, а.
На временной зависимости смещений скачка плотности в средней точке бас-
сейна (рис. 5, б) можно выделить три интервала усиления бароклинных коле-
баний жидкости. Хронологически они соответствуют, во-первых, прохожде-
нию через точку наблюдения волнового следа за барическим фронтом, во-
вторых, прохождению отраженного от правой границы пакета внутренних
волн, в-третьих, прохождению пакета внутренних волн после его отражения
от левой границы бассейна. Далее процесс повторяется. Все эти события раз-
делены во времени интервалами слабых колебаний скачка плотности.
Рассмотрим горизонтальную структуру гидродинамических полей в вол-
новом пакете, образовавшемся после прохождения над бассейном барическо-
го фронта. Типичный случай показан на рис. 6.
Р и с. 6. Горизонтальное распределение смещений границ 1,2 (а), горизонтальных проекций
скорости u1,2 (б) и v1,2 (в) в момент времени t = 502,1 ч (кривые 1 соответствуют j = 1, кривые 2
– j = 2; ширина фронта W = 50 км, скорость С = 2 м·с–1)
ISSN 0233-7584. Мор. гидрофиз. журн., 2009, № 3
13
Представленные на рис. 6, а распределения смещений границ указывают
на то, что в волновом пакете преобладает бароклинная составляющая и его
можно классифицировать как пакет внутренних волн. Так, колебания свобод-
ной поверхности и скачка плотности происходят в противофазе, направления
горизонтальных скоростей течения в слоях на одной вертикали противопо-
ложны (рис. 6, б, в). Наконец, амплитуда колебаний границы раздела слоев
значительно превосходит амплитуду колебаний свободной поверхности жид-
кости (рис. 6, а). Отметим интенсификацию течения в нижнем слое у левой
боковой границы бассейна (рис. 6, б, в), хотя основное тело пакета находится
в средней части бассейна. Такая особенность поля скорости связана с суще-
ствованием у этой границы ярко выраженной шельфовой зоны (см. рис. 2),
где толщина нижнего слоя меньше толщины верхнего.
Для анализа зависимости характеристик поля внутренних волн от пара-
метров внешнего воздействия выполнены расчеты колебаний двухслойной
жидкости в бассейне для различных значений скорости переноса и ширины
атмосферного фронта. Как отмечалось выше, интенсивность бароклинных
колебаний жидкости в бассейне после ухода барического фронта существен-
но зависит от скорости его переноса. Можно ожидать наиболее эффективной
генерации внутренних волн при движении фронта со скоростью С, близкой к
скорости распространения свободных длинных внутренних волн в централь-
ной части бассейна (с2 = 1,53 мc
–1
).
Обозначим через |),(|max 2,1
,
2,1 txA
tx
максимальные смещения свобод-
ной поверхности и границы раздела слоев на временном интервале от момен-
та выхода фронта за пределы бассейна до момента времени t = 600 ч. Зависи-
мости A1,2 от скорости переноса фронта С представлены на рис. 7.
Р и с. 7. Зависимость максимальных смещений свободной поверхности |),(|max 1
,
1 txA
tx
(кривая 1) и границы раздела слоев |),(|max 2
,
2 txA
tx
(кривая 2) от скорости С движения
фронта шириной W = 50 км после выхода фронта за пределы бассейна
Вклад баротропного компонента волнового поля (A1) относительно мал и
увеличивается (немонотонно) с ростом скорости фронта. Что касается вклада
бароклинного компонента в колебаниях скачка плотности (A2), то он наибо-
лее значителен при скорости С ≈ 1,8 мc
–1
, которая достаточно близка к ско-
ISSN 0233-7584. Мор. гидрофиз. журн., 2009, № 3
14
рости распространения длинных внутренних волн в центральной части бас-
сейна. Эффективность генерации внутренних волн при С < с2 быстро ослабе-
вает с уменьшением скорости переноса фронта.
При фиксированном перепаде p атмосферного давления поперек фронта
увеличение ширины фронта W приводит к уменьшению градиента барическо-
го поля (по модулю), который в уравнениях (2) и (4) выступает в роли гори-
зонтальной массовой силы. Таким образом, увеличение ширины фронта при
постоянном значении p должно приводить к уменьшению интенсивности
внешнего воздействия, а в результате – к ослаблению эффективности генера-
ции волн в замкнутом бассейне. Этот вывод подтверждает рис. 8.
Р и с. 8. Возмущения свободной поверхности (штриховая линия) и границы раздела слоев
(сплошная линия) для барического фронта различной ширины: а – W = 50 км; б – W = 100 км;
в – W = 500 км (момент времени t = 83 ч, скорость фронта C = 2 мс–1)
На всех рисунках фронт различной ширины располагается в средней час-
ти бассейна. Горизонтальная структура смещений скачка плотности форми-
руется из двух составляющих. Первая из них представляет собой некоторое
плавное среднее распределение Hh /21
*
2 , соответствующее баротроп-
ной реакции двухслойной жидкости на барическое воздействие. Вторая (от-
клонение 2 от *
2 ) – внутренние волны, влияние которых на форму свобод-
ISSN 0233-7584. Мор. гидрофиз. журн., 2009, № 3
15
ной поверхности малό. Внутренние волны наиболее значительны для узкого
фронта (рис. 8, а). Широкий фронт (рис. 8, в) – малоэффективный генератор
бароклинных волн в ограниченном бассейне. Для всех трех случаев, пред-
ставленных на рис. 8, наблюдаются значительные количественные отклонения
смещений поверхности жидкости от приближения обратного барометра (1).
В ограниченном бассейне постоянной глубины период бароклинных
стоячих колебаний жидкости находится по формуле [26]
2
2
21
22
2
f
Hl
hghm
Tm
, (16)
где m = 1, 2, … – номер (число узловых линий границы раздела слоев) барок-
линной моды. Как следует из (16), вне зависимости от номера моды колеба-
ний, уменьшение толщины верхнего слоя при h1 < h2 и H = const и относи-
тельного перепада плотности ε или увеличение ширины бассейна l вызывают
рост периодов колебаний бароклинных сейш.
Р и с. 9. Колебания скачка плотности 2 в средней точке бассейна x = l/2 при фиксированной
плотности верхнего слоя жидкости 1 = 1012 кгм–3, но различных перепадах плотности
12 между слоями: 1 – = 2,5 кгм–3; 2 – = 5 кгм–3; 3 – = 10 кгм–3 (ширина
фронта W = 50 км, скорость фронта С = 2 м·с–1)
Рис. 9 подтверждает зависимость периода остаточных бароклинных ко-
лебаний жидкости в бассейне от интенсивности скачка плотности: период
внутренних волн убывает с ростом перепада плотности ∆ρ между слоями.
Заключение. В рамках нелинейной теории длинных волн рассмотрена
плоская задача о генерации поверхностных и внутренних волн при переме-
щении над ограниченным вращающимся бассейном переменной глубины
фронта атмосферного давления. Барический фронт моделируется зоной мо-
нотонного увеличения/уменьшения атмосферного давления от одного посто-
янного значения к другому. Задача решалась численно методом конечных
разностей по явно-неявной схеме. Расчеты выполнены для бассейна, распре-
ISSN 0233-7584. Мор. гидрофиз. журн., 2009, № 3
16
деление глубины которого соответствует зональному сечению Черноморской
котловины. Глубина верхнего слоя принята постоянной.
Показано, что при прохождении над бассейном барического фронта ге-
нерируются поверхностные и внутренние волны. Их параметры существенно
зависят как от геометрии бассейна, так и от характеристик внешнего воздей-
ствия. При медленном перемещении фронта (C ≈ 1 мс
–1
) смещения свобод-
ной поверхности и скачка плотности повторяют по форме распределение ат-
мосферного давления, перемещаются вместе с фронтом и практически совпа-
дают, если толщина верхнего слоя относительно невелика. Подобный отклик
двухслойной жидкости на внешнее воздействие является баротропным. Сме-
щения свободной поверхности существенно отличаются от гидростатическо-
го, чему в силу сохранения массы жидкости способствует в немалой степени
ограниченность бассейна.
Что касается колебаний скачка плотности после выхода фронта за преде-
лы бассейна, то наиболее интенсивные внутренние волны генерируются при
скоростях движения фронта, близких к скорости распространения длинных
внутренних волн в центральной части бассейна. Баротропный компонент
усиливается по мере приближения скорости переноса фронта к резонансной
скорости возбуждения поверхностных волн gHC , которая значительно
превосходит возможные скорости переноса атмосферных возмущений в ре-
альных условиях.
Волновая реакция двухслойной жидкости на движущийся фронт зависит
от его ширины. Бароклинные возмущения среды наиболее значительны для
узкого фронта, когда градиент атмосферного давления велик. Широкий
фронт – малоэффективный генератор бароклинных возмущений жидкости в
бассейне. Во многих случаях наблюдаются значительные количественные
отклонения смещений поверхности жидкости от приближения обратного ба-
рометра.
После выхода фронта за пределы бассейна возникающие внутренние
волны представляют собой волновой пакет, распространяющийся от одной
боковой границы к другой. Такие внутренние волны вряд ли можно рассмат-
ривать как традиционные бароклинные сейши, поскольку они не сопровож-
даются синхронными движениями всей массы жидкости в бассейне, а много-
численные узловые линии совершают существенные горизонтальные пере-
мещения вместе с волновым пакетом.
Проведенные расчеты создаваемых движущимся фронтом возмущений
свободной поверхности жидкости и границы раздела слоев, а также течений в
слоях, не выявили существенного влияния нелинейности на бароклинную
динамику жидкости в бассейне рассматриваемой геометрии. Слабым нели-
нейным искажениям подвержены внутренние волны наибольшей амплитуды
в центральных частях бароклинных волновых пакетов. Как показано в работе
[18], нелинейные эффекты в поверхностных сейшах наиболее ярко выражены
на мелководье, в частности при перемещении барических фронтов над аква-
торией Азовского моря.
ISSN 0233-7584. Мор. гидрофиз. журн., 2009, № 3
17
СПИСОК ЛИТЕРАТУРЫ
1. Праудмэн Дж. Динамическая океанография. – М.: Иностранная литература, 1957. – 418 с.
2. Ламб Г. Гидродинамика. – М.; Л.: Гостехиздат, 1947. – 928 с.
3. Miles W.J. Harbor seiching // Annual Rev. Fluid. Mech. – 1974. – 6. – P. 17 – 33.
4. Судольский А.С. Динамические явления в водоемах. – Л.: Гидрометеоиздат, 1991. – 263 с.
5. Ichinose G.A., Anderson J.G., Satake K. et al. The potential hazard from tsunami and seiche
waves generated by large earthquakes within Lake Tahoe, California-Nevada // Geophys. Res.
Letters. – 2000. – 27, № 8. – P. 1203 – 1206.
6. Chapman D.C., Giese G.S. A model for the generation of coastal seiches by deep-sea internal
waves // J. Phys. Oceanogr. – 1990. – 20, № 9. – P. 1459 – 1467.
7. Abraham E.R.C. Seiche modes of Wellington Harbour, New Zealand // New Zealand J. Mar. and
Freshwater Res. – 1997. – 31, № 2. – P. 191 – 200.
8. Wedderburn E.M. Temperature observations in Loch Earn, with a further contribution to the
hydrodynamical theory of the temperature seiches // Trans. Roy. Soc. Edinburgh. – 1912. – 48. –
P. 629 – 695.
9. Mortimer C.H. Long internal waves in lakes: Review of a century of research // Univ. Wiscon-
sin-Milwaukee, Center for Great Lakes studies. – Spec. Rep. – 1993. – № 42. – 177 p.
10. Gomes D., Monserrat S., Tintore J. Pressure-forced seiches of large amplitude in inlets of the
Balearic Islands // J. Geophys. Res. – 1993. – 98, № C8. – P. 14437 – 14445.
11. Garcies M., Gomes D., Monserrat S. Pressure-forced seiches of large amplitude in inlets of the
Balearic Islands. 2. Observational study // Ibid. – 1996. – 101, № C3. – P. 6453 – 6467.
12. Rabinovich A., Monserrat S. Generation of meteorological tsunami (large amplitude seiches)
near the Balearic and Kuril Islands // Natural Hazards. – 1998. – 18. – P. 27 – 55.
13. de Jong M.P.C., Holthuijen L.H., Battjes J.A. Generation of seiches by cold fronts over the
southern North Sea // J. Geophys. Res. – 2003. – 108, № С4. – P. 14 – 19.
14. Vilibić I., Domijan N., Orlić M. et al. Resonant coupling of a traveling air pressure disturbance
with the east Adriatic coastal waters // Ibid. – 2004. – 109, № C10001,
doi:10.1029/2004JC002279.
15. Коновалов А.В., Черкесов Л.В. Генерация длинных нелинейных волн в замкнутом бассей-
не движущимися возмущениями атмосферного давления // Изв. РАН. Физика атмосферы и
океана. – 1995. – 31, № 5. – С. 713 – 718.
16. Еремеев В.Н., Коновалов А.В., Черкесов Л.В. Моделирование длинных баротропных волн в
Черном море, вызываемых движущимися барическими возмущениями // Океанология. –
1996. – 36, № 2. – C. 191 – 196.
17. Еремеев В.Н., Коновалов А.В., Манилюк Ю.В., Черкесов Л.В. Моделирование длинных
волн в Азовском море, вызываемых прохождением циклонов // Там же. – 2000. – 40, № 5.
– C. 658 – 665.
18. Доценко С.Ф., Миклашевская Н.А. Генерация сейш при перемещении барических образо-
ваний // Морской гидрофизический журнал. – 2007. – № 6. – С. 3 15.
19. Доценко С.Ф., Миклашевская Н.А. Генерация сейш в ограниченных бассейнах переме-
щающимися барическими фронтами // Там же. – 2008. – № 2. – С. 3 – 18.
20. Horn W., Mortimer C.H., Schwab D.J. Wind-induced internal seiches in Lake Zurich observed
and modeled // Limnol. Oceanogr. – 1986. – 31, № 6. – P. 1232 – 1254.
21. Lemmin U., Mortimer C.H. Tests of an extension to internal seiches of Defant’s procedure for
determination of surface seiche characteristics in real lakes // Ibid. – 1986. – 31, № 6. –
P. 1207 – 1231.
ISSN 0233-7584. Мор. гидрофиз. журн., 2009, № 3
18
22. Arneborg L., Liljebladh B. The internal seiches in Gullmar Fjord. Part I: Dynamics // J. Phys.
Oceanogr. – 2001. – 31, № 9. – P. 2549 – 2566.
23. Sirkes Z. Surface manifestations of internal oscillations in a highly saline lake (the Dead Sea) //
Limnol. Oceanogr. – 1987. – 32, № 1. – P. 76 – 82.
24. Лаппо С.С. Среднемасштабные динамические процессы океана, возбуждаемые атмосфе-
рой. М.: Наука, 1979. 181 с.
25. Вольцингер Н.Е., Пясковский Р.В. Основные океанологические задачи теории мелкой во-
ды. – Л.: Гидрометеоиздат, 1968. – 300 с.
26. Ле Блон П., Майсек Л. Волны в океане. Т. 1. – М.: Мир, 1981. – 480 с.
27. Доценко С.Ф. Генерация длинных внутренних волн в океане движущейся областью давле-
ний // Морской гидрофизический журнал. – 1990. – № 3. – С. 3 – 9.
Морской гидрофизический институт НАН Украины, Материал поступил
Севастополь в редакцию 20.02.08
ABSTRACT Plane problem on generation of surface and internal waves in a bounded rotating basin
of variale depth by an atmospheric pressure front moving above the basin is considered. The fluid is
assumed two-layer. The nonlinear system of long-wave equations is solved numerically by the meth-
od of finite differences for the depth distribution corresponding to the zonal section of the Black Sea
basin. It is shown that barotropic and baroclinic fluid oscillations are generated as a result of a baric
front translation above the basin. Intensity of fluid disturbances depends on translation velocity and
the front width. There is a velocity of front translation at which the internal waves are generated most
effectively. After the front leaves the basin a packet of internal waves which propagates reflecting
from one lateral boundary to the other is formed.
|