Научная статья на тему 'Исследование свободных колебаний жидкости в ограниченном бассейне, представляющем приближенную модель Азовского моря'

Исследование свободных колебаний жидкости в ограниченном бассейне, представляющем приближенную модель Азовского моря Текст научной статьи по специальности «Физика»

CC BY
111
25
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
Ключевые слова
СЕЙШИ / СЕЙШЕООБРАЗНЫЕ КОЛЕБАНИЯ / СВОБОДНЫЕ ВОЛНЫ / ДЛИННЫЕ ВОЛНЫ / ТЕЧЕНИЯ / АЗОВСКОЕ МОРЕ / АНАЛИТИЧЕСКИЕ РЕШЕНИЯ / SEICHES / SEICHE-LIKE OSCILLATIONS / FREE WAVES / LONG WAVES / FLOWS / SEA OF AZOV / ANALYTIC SOLUTION

Аннотация научной статьи по физике, автор научной работы — Манилюк Ю.В., Черкесов Л.В.

На основе аналитического решения системы уравнений длинных волн рассмотрены особенности сейшевых и сейшеообразных колебаний (комбинация мод сейш) в модельном бассейне, имеющем характерные размеры Азовского моря. Установлено, что колебания уровня моря при сейшеобразных колебаниях принципиально отличаются от колебаний уровня при сейшевых колебаниях. Сейшеобразные колебания имеют одну узловую линию, которая поворачивается против часовой стрелки с течением времени. Максимальные возвышения уровня в отдельных областях акватории водоема могут приближаться к максимально возможным значениям, равным сумме амплитуд взаимодействующих волн и имеющим место в начальный момент времени. Максимальное значение модуля скорости волновых течений во всех рассмотренных случаях существенно меньше максимально возможного значения. Направление векторов волновых течений постоянно изменяется со временем, поворачиваясь по часовой стрелке.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Investigation of free liquid oscillations in a bounded basin representing an approximate model of the Azov Sea

Based on the analytic solution of long wave equations, the features of seiches and seiche-like oscillations (combination of seiche modes) in the model basin with the dimensions characteristic of the Azov Sea are considered. It is found that the sea level fluctuations induced by the seiche-like oscillations are fundamentally different from those resulting from the seiche oscillations. The seiche-like fluctuations have a single nodal line turning counter-clockwise with time. The level maximum elevations in certain areas of the basin can approached their maximum possible values which equal the sum of the interacting waves’ amplitudes and occur at the initial moment of time. In all the cases under consideration the maximum value of the wave flows’ velocity module is significantly smaller than the maximum possible one. The direction of the wave flows’ vectors is constantly changing with time turning clockwise.

Текст научной работы на тему «Исследование свободных колебаний жидкости в ограниченном бассейне, представляющем приближенную модель Азовского моря»

УДК 532.59 + 551.466

Исследование свободных колебаний жидкости в ограниченном бассейне, представляющем приближенную модель Азовского моря

© 2016 Ю.В. Манилюк, Л.В. Черкесов

Морской гидрофизический институт РАН, Севастополь, Россия E-mail: uvmsev@yandex.ru

Поступила в редакцию 20.10.2015 г. После доработки 26.11.2015 г.

На основе аналитического решения системы уравнений длинных волн рассмотрены особенности сейшевых и сейшеообразных колебаний (комбинация мод сейш) в модельном бассейне, имеющем характерные размеры Азовского моря.

Установлено, что колебания уровня моря при сейшеобразных колебаниях принципиально отличаются от колебаний уровня при сейшевых колебаниях. Сейшеобразные колебания имеют одну узловую линию, которая поворачивается против часовой стрелки с течением времени. Максимальные возвышения уровня в отдельных областях акватории водоема могут приближаться к максимально возможным значениям, равным сумме амплитуд взаимодействующих волн и имеющим место в начальный момент времени. Максимальное значение модуля скорости волновых течений во всех рассмотренных случаях существенно меньше максимально возможного значения. Направление векторов волновых течений постоянно изменяется со временем, поворачиваясь по часовой стрелке.

Ключевые слова: сейши, сейшеообразные колебания, свободные волны, длинные волны, течения, Азовское море, аналитические решения.

Введение. Сейши являются одним из наиболее распространенных явлений в ограниченных водоемах и вносят значительный вклад в формирование их гидрологических режимов[1, 2]. Это обусловлено тем, что причин, вызывающих их, много, и для их генерации не требуется большого количества энергии [2, 3].

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

Теоретическое исследование сейш в реальных водоемах затруднено тем, что водоемы имеют сложную геометрическую форму и переменную глубину. В силу этого получить аналитические решения уравнений гидродинамики для них пока не удается. Поэтому приходится решать задачу численно. Для этого используются два подхода. При первом моделируется воздействие на поверхность водоема возмущающей силы, после прекращения действия которой исследуются колебания жидкости в водоеме. Возмущения могут быть весьма разнообразными в виде: циклона [4]; движущегося барического фронта [5]; волнопродуктора на жидкой границе бассейна, генерирующего волны на резонансных частотах [6]; поля ветра [7]. При этом неизвестно, из чего складывается результирующее колебание поверхности моря. Предполагается, что оно представляет собой суперпозицию нескольких мод сейш. При втором

16

МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 2 2016

подходе численно решается задача на собственные значения. Здесь также используется широкий набор разных методов на основе: метода конечных разностей [8]; метода конечных элементов [9]; ЛГ-алгоритма конформного отображения [10]. Так, для Азовского моря на основе метода конечных элементов получено решение задачи на собственные значения в статье [9]. В результате рассчитаны периоды и определена пространственная структура мод сейш.

Аналитические решения позволяют получить более полную картину явления и исследовать его детальнее. В работах [11 - 13] для некоторых бассейнов простой формы приведены формулы для расчета периодов сейш и отклонений свободной поверхности. Из бассейнов, для которых можно найти аналитическое решение, к Азовскому морю по своей форме наиболее близок прямоугольный бассейн постоянной глубины.

Целью данной работы является нахождение в линейном приближении аналитического решения задачи о сейшах для прямоугольного бассейна постоянной глубины, которое позволяет рассчитать периоды мод сейшевых колебаний, отклонения свободной поверхности и скорости волновых течений. На основе полученного решения исследуются сейшевые и сейшеообразные колебания. В работе [14] показано, что после прохождения над Азово-Черноморским бассейном барического фронта могут возникать низшие моды баротропных сейш. Поэтому нами рассмотрены сейшеобразные колебания, представляющие собой суперпозицию первых двух мод сейш.

Постановка задачи и получение аналитического решения. Рассмотрим свободные колебания однородной жидкости в замкнутом бассейне. Предположим, что волны длинные и атмосферное давление над акваторией бассейна постоянное. В таком случае движение жидкости описывается системой уравнений [15]

и = -, V = -^^^ , С = -(ки)х - к)у , (!)

где использованы следующие обозначения: х, у - декартовы координаты, и, V - проекции вектора скорости движения жидкости на осях х, у соответственно, С, - отклонение свободной поверхности, к - глубина бассейна, ^ -время, g - ускорение свободного падения.

На твердых боковых стенках бассейна должно выполняться условие непротекания

V = 0. (2)

Здесь Уп - нормальная к боковой стенке бассейна составляющая скорости.

Рассмотрим прямоугольный бассейн, имеющий постоянную глубину к. Длина бассейна равна а , ширина - Ь . С учетом этого граничные условия (2) примут вид

и(0, у) = 0, и (а, у) = 0; v(х,0) = 0, v(х, Ь) = 0. (3)

Будем искать решение системы (1) в виде функций, периодических по времени и удовлетворяющих граничным условиям (3):

МОРЖОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 2 2016

17

икт (х У,1) = икт (у) 81п(кпха-) вш актГ, укт(х,у,0 = Укт(х)81п(тяуЬ~ )81п

<кт , (4)

Скт (X, У, ') = Скт (X, У)С08акт( ,

где сгкт - частота колебаний жидкости, к = 0,1, 2,..., ; т = 0,1, 2,...,. Подставив (4) в систему (1), получим

= -<^йкт (у)81п П, (5)

дх g а

д£кт (х, у) акт_ тпу

т =--^Укт (х) 81п. (6)

ду g Ь

Интегрируя уравнения (5) по х, а (6) по у, имеем

Скт (х, у) = ^аЩт (^008^ , (7)

gkл а

^(х,у) = ^кт(х)008. (8)

gmл Ь

Будем искать йкт (у) и Укт (х) в виде

йкт (у) = и0кт COS(mЛyЬ~1), ^кт (х) = ^0кт 008(кпха- ) . (9)

Здесь и0кт , У0кт - амплитуды составляющих скорости.

Приравнивая правые части уравнений (7), (8) друг к другу и заменяя йкт (у) и Укт (х) на выражения из (9), находим соотношение между амплитудами составляющих скорости

ат

У0кт =ТТ и0кт . (10)

Ьк

Из третьего уравнения системы (1), учитывая (4), (8), получим

2 х 2 Л

<У2та кпх тпу . _ ^ , __ к + та

~и0кт 008-] 81П <кт^ = -Ит0кт

1 и/ст т

gkл а Ь

ча кЬ2 ,

кпх тпу .

008-008-81П .

аЬ

Отсюда находим выражения для вычисления частот сейшевых колебаний жидкости в бассейне прямоугольной формы постоянной глубины:

<кт =^^2 + (11)

2 аЬ

Т кт = ~!=г , 2 2 • (12)

Л^ л/22Ь2 + т Ь2

18

МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 2 2016

С учетом (7), (9) формула (4) примет вид

^ / .ч crkma knx mny

Zkm ( X У,t) =—^U0km COS-COS—— COS ckmt . (13)

gkn a b

Выражение (13) можно записать так:

„ , . „ knx mny /л ..

Zkm (X У,t) = Z0 km COS-COS^~COSCkmt , (14)

a b

где Z0km - амплитуда отклонения свободной поверхности.

Формула (14) позволяет найти профиль свободной поверхности жидкости для заданных значений k, m, Z0km и геометрических параметров бассейна

a, b .

Используя (13) и (14), получим

Ckma

gkn

Z0km = U0km . (15)

С учетом выражения (11) соотношение (15) примет вид

ki+m

k g \ a2 + b2

Z0km = , U0kmj I 2 + , 2 . (16)

Отсюда имеем формулу для вычисления значений и0кт по заданному значеНИЮ Сокт :

-Л [к кЬ (17)

Щокт 2 2' ( ) » " л/к Ь + та

Подставляя выражение (17) в (10), получим формулу для вычисления у0кт по заданному значению ^0кт:

V -г [й та (18)

0кт 1,2,2 ^ 2 2' 1)

« п л/к Ь + т а

Используя (4), (9), (17), (18), выпишем формулы для расчета составляющих волновых скоростей:

, ч „ к кЬ . кпх тпу .

икт (X У, 1) - ^0кт д Т" I 2 2 2 2 вт-СО^ —— 81П , (19)

2Ь2 + т2а2 а Ь

/ ч ^ 18 та кпх . тпу . .....

vkm(XУ, 1) -^кт,Т I 2 2 2 -вш—— это/ • (20)

2Ь2 + та2 а Ь

Таким образом, для прямоугольного бассейна постоянной глубины получено решение системы уравнений длинных волн, позволяющее найти периоды мод сейшевых колебаний по формуле (12) и, задавая начальные значения амплитуд отклонений свободной поверхности £"0кт, рассчитать на основе со-

МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 2 2016

19

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

отношения (14) отклонения свободной поверхности и с помощью выражений (19), (20) волновые скорости соответствующих мод.

Сейшевые колебания. Анализ формул (14), (19), (20) позволяет разделить сейши на три вида. Вид определяется значениями индексов к, т. У сейш первого вида к Ф 0, т - 0, а и, V, ^ не зависят от координаты у . У сейш второго вида к - 0, т Ф 0, а и, V, £ не зависят от координаты х. Третий вид сейш имеет к Ф 0, т Ф 0, при этом и, V, £ зависят и от координаты х, и от координаты у.

Количество узловых линий уровня у каждой моды сейш равно сумме индексов к + т . Число узловых линий, параллельных оси у, равно к , а параллельных оси х, равно т .

Рассмотрим особенности волновых движений, имеющих место для всех указанных выше видов сейш. Расчеты проведены для прямоугольного бассейна длиной 450 км, шириной 250 км, глубиной 10 м, который является приближенной моделью Азовского моря, имеющей его характерные размеры.

В таблице приведены для первых восьми мод сейшевых колебаний значения периодов сейш и максимальные значения модуля скорости волновых

течений, рассчитанные по формуле |Укт| и2кт + vlm для амплитуды отклонения свободной поверхности, равной 0,3 м.

Периоды сейш и максимальные значения модуля скорости волновых течений

Номер моды Период, ч к т тЗХ, м/с

1 25,25 1 0 0,29

2 14,03 0 1 0,29

3 12,63 2 0 0,29

4 12,26 1 1 0,26

5 9,39 2 1 0,22

6 7,01 0 2 0,29

7 6,76 1 2 0,28

8 6,30 2 2 0,26

Как видно из таблицы, наибольший период, равный 25,25 ч, имеет продольная одноузловая сейша. Периоды остальных мод значительно меньше: 14,03 ч (одноузловая поперечная сейша), 12,63 ч (двухузловая продольная сейша), 12,26 ч (двухузловая продольно-поперечная сейша). Периоды сейш третьей (12,63 ч) и четвертой (12,26 ч) мод отличаются друг от друга незначительно. Это указывает на то, что не всегда можно однозначно идентифицировать моды сейшевых колебаний по длительности периодов. Максимальные скорости волновых течений одно- и двухузловых сейш первого ( к Ф 0, т - 0 ) и второго (к - 0, т Ф 0) видов равны 0,29 м/с. Максимальные скорости волновых течений сейш третьего вида ( к Ф 0, т Ф 0 ) всегда меньше, чем сейш первого и второго видов. Из всех рассмотренных здесь мод сейш третьего

20

МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 2 2016

вида наибольшую скорость волновых течений, равную 0,28 м/с, имеет седьмая мода (к = 1, т = 2), а наименьшую, равную 0,22 м/с, - пятая мода (к = 2, т = 1). Максимальные скорости волновых течений сейш ((к = 1, т = 1) и (к = 2, т = 2 )), имеющих одинаковое количество узловых линий, перпендикулярных как оси х, так и оси у, равны друг другу и составляют 0,26 м/с.

Рассмотрим подробнее структуру колебаний уровня. У сейш первого и второго видов все изолинии уровня параллельны. При этом у сейш первого вида изолинии параллельны оси у, а у сейш второго вида - оси х. В силу этого области, в которых отклонения уровня максимальны, расположены в окрестностях отрезков, параллельных одной из осей координат. Две из этих областей обязательно совпадают с боковой границей бассейна. У сейш третьего вида параллельны осям координат только узловые линии, остальные изолинии представляют собой кривые, расположенные в прямоугольниках, образованных пересекающимися узловыми линиями. Размеры областей, в которых отклонения уровня максимальны, значительно меньше, чем у сейш первого и второго видов. Для каждого из трех видов сейш, рассмотренных нами, в окрестности всех угловых точек бассейна существуют области, в которых отклонения уровня максимальны.

Сейшевые колебания вызывают течения, охватывающие весь объем жидкости в бассейне. Течение жидкости для сейш первого вида происходит параллельно оси х, а для сейш второго вида - параллельно оси у . Для сейш третьего вида картина волновых течений сложнее, чем для сейш первого и второго видов. Эти течения происходят под различными углами к осям координат.

Для сейш первого и второго видов области наибольших скоростей течений располагаются в окрестности узловых линий. Для сейш третьего вида изолиния скорости (для каждой моды она имеет свою величину) разделяет акваторию бассейна на ячейки. Выделяются ячейки, в которых скорость минимальная, и ячейки, в которых течения наиболее интенсивные. Расположение зон интенсивных волновых течений существенно зависит от номера моды. Имеются ячейки с постоянной скоростью, составляющей около 50% от максимального значения. Для каждой моды сейш выделяются ячейки, примыкающие к углам бассейна, в которых течения имеют минимальную интенсивность.

Для всех трех видов сейш встречаются как сходящиеся, так и расходящиеся потоки. Направления волновых течений изменяются на противоположные каждые полпериода соответствующей моды. Модуль скорости первые полпериода возрастает от нуля до максимального значения, а затем следующие полпериода убывает до нуля.

Сейшеобразные колебания. Рассмотрим взаимодействие двух мод сей-шевых колебаний на примере суперпозиции и .

На рис. 1, а приведены максимальные для всего бассейна значения отклонений уровня £о1 + за интервал времени 3т10 с шагом &, где т10- период первой моды сейшевых колебаний, & = т10 /20 . При этом амплитуды отклонений уровня взаимодействующих мод в начальный момент времени равны 0,3 м.

При ? = 0 отклонение уровня максимально и составляет 0,6 м. С течением времени значение максимума убывает, достигает величины 0,24 м через МОРЖОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 2 2016 21

3Д? (3,78 ч), затем возрастает до 0,39 м к моменту времени 6Д? (7,56 ч). После этого снова убывает до 0,32 м при ? = 7Д? (8,82 ч), затем возрастает до 0,58 м при ? = 10Д? (12,6 ч). Аналогичные колебания уровня происходят в следующие временные интервалы (т10, 2т10 ] и (2т10, 3т10 ]. Таким образом, за интервал времени, равный т10, максимум уровня совершает два колебания. В этом сказывается влияние более короткопериодной моды (т01/т10 = 0,55). Расчеты показали, что при взаимодействии двух мод максимальная амплитуда результирующей волны может быть близка к сумме амплитуд взаимодействующих мод и в данном случае достигала 98% от возможного максимума в момент времени, равный 20Д? (25,2 ч), и 88% от возможного максимума в момент времени, равный 40 Д? (50,4 ч). Отметим также, что в течение расчетного интервала времени максимум амплитуды результирующей волны был больше 10 см, за исключением узкого интервала в окрестности момента времени 45Д? (56,7 ч), в который максимум равен нулю. При этом значения максимального возвышения и максимального понижения уровня равны по модулю.

к

я =

и1

6_ ^ 0.2'

Т

20

—Г

40

Д( = Т10/20 а

Т

60

П

80

40

Д( = Т10/20 б

Рис. 1. Зависимость от времени максимумов отклонений уровня свободной поверхности и модуля скорости волновых течений: а - ^01 + С10, б - |У01 + У101

Теперь рассмотрим, как влияет взаимодействие двух мод на структуру колебаний уровня жидкости в бассейне. На рис. 2 показаны изолинии уровня результирующего колебания + £10 для различных моментов времени в течение интервала, равного периоду т10, для случая, когда амплитуды отклонений уровня взаимодействующих мод в начальный момент времени составляют ^001 = С010 = 0,3 м . В момент времени ? = 0 (рис. 2, а) ^тах = 0,6 м, при

* = Т10 /2 (рис. 2, б) Сшх = 0,29 м и при * = Т10 (рис. 2, в) Спах = 0,58 м.

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

22 МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 2 2016

0.6

0.5

0.4

0.4

0.3

0.1

20

60

80

поворачиваясь против часовой стрелки вокруг центра бассейна. Через интервал времени, равный периоду продольной одноузловой сейши г10, узловая

линия совершает поворот на 90°. Аналогичным образом перемещаются максимальные отклонения уровня. При этом в каждый момент времени минимальное и максимальное отклонения равны по модулю.

в

Рис. 2. Изолинии уровня результирующего колебания жидкости в бассейне £01 + Сю ■ а -/ = 0, б - / = т10/ 2 , в - / = т10

Также были рассмотрены случаи, когда £001 >С010, С001 < 4>010. Установлено, что разница амплитуд не влияет на характер движения. Изменяется только начальное положение этой узловой линии. Если £001 > Со10, то оно

МОРЖОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 2 2016

23

сдвинуто против часовой стрелки относительно положения, когда £001 = С010, а если С001 < С010, то оно сдвинуто по часовой стрелке.

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

На рис. 1, б приведены максимальные значения модуля скорости волновых течений |У01 + У10| в интервале времени продолжительностью 3 т10. Амплитуды отклонений уровня для обеих мод в начальный момент времени равны 0,3 м. Результаты расчетов показали, что максимальное значение модуля скорости за расчетный период равно 0,42 м/с, что составляет 72% от суммы максимальных амплитуд скоростей обеих мод и достигается в моменты времени 26 Дt (32,76 ч) и 48 Дt (60,48 ч). Также в момент времени t = 8Дt (10,08 ч) имеет место локальный максимум, равный 0,38 м/с (66% от максимально возможного значения). За интервал времени, равный т10, максимум модуля скорости совершает два колебания. При этом он на интервале от 0 до 3 т10 больше 10 см/с за исключением небольшого отрезка в окрестности момента времени t = 0 , в который скорость волновых течений равна нулю. Таким образом, в отличие от отклонения уровня за расчетный период, равный 3т10, модуль скорости не превышает 72% от максимально возможного значения и обращается в нуль один раз.

\ у,

/ / у

y, км ^ * м *

200^ 4 к £ ■/+■/*■/</</„

. АШВ

100-

У, км

200-:

150- ;

100- -

50-:

0-^ 0

а

-т!-г-

50 100 150 200

250 б

300 350 400 x, км

Рис. 3. Векторные поля скорости волновых течений У01 + У10 ■ а -1 = 0,3т10 , б -1 = 0,5т10 24 МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 2 2016

0

На рис. 3 для моментов времени t = 0,3г10 и t = 0,5г10 приведены векторные поля скорости волновых течений V01 + V10, которые иллюстрируют факт

того, что в результате взаимодействия двух мод сейш образуется система волновых течений, принципиально отличающаяся от системы волновых течений, имеющих место для одиночных мод сейш. В случае рассмотренных сейшеобразных колебаний направление течений постоянно изменяется со временем, поворачиваясь по часовой стрелке. При этом встречные и расходящиеся потоки отсутствуют.

Выводы.

1. В результате взаимодействия нескольких мод сейш получается непериодическое колебание, имеющее одну узловую линию, которая поворачивается вокруг центра бассейна с течением времени против часовой стрелки.

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

3. Максимальная величина модуля скорости волновых течений во всех рассмотренных случаях меньше максимально возможного значения и не превышает 72% от него.

4. Волновые течения для сейшеобразных колебаний принципиально отличаются от волновых течений для отдельных сейш. При сейшевых колебаниях течения имеют направление, которое является постоянным в продолжение полупериода, а затем изменяется на противоположное. Также для сейш характерно наличие встречных и расходящихся потоков. В случае сейшеоб-разных колебаний направление волновых течений постоянно изменяется со временем, при этом векторы скорости поворачиваются по часовой стрелке. Встречные и расходящиеся потоки отсутствуют.

СПИСОК ЛИТЕРАТУРЫ

1. Океанографическая энциклопедия. - Л.: Гидрометеоиздат, 1980. - 304 с.

2. Дарвин Д. Приливы и родственные им явления в солнечной системе. - М.: Наука, 1961. - 252 с.

3. Лабзовский Н.А. Непериодические колебания уровня моря. - Л.: Гидрометеоиздат, 1971. - 237 с.

4. Еремеев В.Н., Коновалов А.В., Манилюк Ю.В., Черкесов Л.В. Моделирование длинных волн в Азовском море, вызываемых прохождением циклонов // Океанология. - 2000. -40, № 5. - С. 658 - 665.

5. Доценко С.Ф., Миклашевская Н.А. Генерация сейш при перемещении барических образований // Морской гидрофизический журнал. - 2007. - № 6. - С. 3 - 13.

6. Матишов Г.Г., Инжебейкин Ю.И. Численное исследование сейшевых колебаний уровня Азовского моря // Океанология. - 2009. - 49, № 4. - С. 485 - 493.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

7. Черкесов Л.В., Шульга Т.Я. Численный анализ пространственной структуры свободных колебаний в Азовском море // Экология, экономика, информатика. - Сб. статей в 3-х томах. Т. 2. Системный анализ и моделирование экономических и экологических систем. - Ростов-на-Дону: Изд-во Южного федерального ун-та, 2015. - С. 689 - 699.

МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 2 2016

25

8. Иванов В.А., Янковский А.Е. Длинноволновые движения в Черном море. - Киев: Науко-ва думка, 1992. - 112 с.

9. Манилюк Ю.В., Черкесов Л.В. Математическое моделирование сейшевых колебаний в Азовском море на основе метода конечных элементов // Морской гидрофизический журнал. - 1994. - № 5. - С. 3 - 8.

10. Рабинович Б.И., Левянт А.С. Численное решение задачи расчета сейш на основе Rr-алгоритма конформного отображения // Природные катастрофы и стихийные бедствия в Дальневосточном регионе. Т. 2. - Владивосток: Изд-во ДВО АН СССР, 1990. -С. 328 - 342.

11. ЛамбГ. Гидродинамика. - М.: Гостехиздат, 1947. - 948 с.

12. Кочин Н.Е., Кибель И.А., Розе Н.В. Теоретическая гидромеханика. Т. 1. - М.: Гостехиздат, 1956. - 560 с.

13. Rabinovich A.B. Seiches and Harbor Oscillations (Chapter 9) // Handbook of Coastal and Ocean Engineering / Ed. Y.C. Kim. - Singapoure: World Scientific Publ., 2009. - P. 193 -236.

14. Доценко С.Ф., Миклашевская Н.А. Генерация сейш в ограниченных бассейнах перемещающимися барическими фронтами // Морской гидрофизический журнал. - 2008. -№ 2. - С. 3 - 19.

15. Сретенский Л.Н. Теория волновых движений жидкости. - М.: Наука, 1977. - 815 с.

Investigation of free liquid oscillations in a bounded basin representing an approximate model of the Azov Sea

Yu.V. Maniliuk, L.V. Cherkesov

Marine Hydrophysical Institute, Russian Academy of Sciences, Sevastopol, Russia e-mail: uvmsev@yandex.ru

Based on the analytic solution of long wave equations, the features of seiches and seiche-like oscillations (combination of seiche modes) in the model basin with the dimensions characteristic of the Azov Sea are considered.

It is found that the sea level fluctuations induced by the seiche-like oscillations are fundamentally different from those resulting from the seiche oscillations. The seiche-like fluctuations have a single nodal line turning counter-clockwise with time. The level maximum elevations in certain areas of the basin can approach their maximum possible values which equal the sum of the interacting waves' amplitudes and occur at the initial moment of time. In all the cases under consideration the maximum value of the wave flows' velocity module is significantly smaller than the maximum possible one. The direction of the wave flows' vectors is constantly changing with time turning clockwise.

Keywords: seiches, seiche-like oscillations, free waves, long waves, flows, Sea of Azov, analytic solution.

26

МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 2 2016

i Надоели баннеры? Вы всегда можете отключить рекламу.