147
Нахождение сейш сводится к решению задачи (4.1), (4.2) на
собственные частоты ...) ,2 ,1(
j
j
колебаний жидкости в бас-
сейне (частоты сейш) и соответствующие им собственные функции
Z = Z
j
(x, y) (распределения смещений свободной поверхности жид-
кости). В работе [2] доказано, что задача (4.1) – (4.3) имеет дискрет-
ный спектр, а система собственных функций является полной.
Только в небольшом числе случаев (бассейны модельной
формы) задача может быть решена аналитически [1, 3]. Нахождение
частот и распределений смещений уровня в сейшах для бассейнов
близкой к реальной геометрии
требует применения численных ме-
тодов. Сейши в Черном и Азовском морях могут быть найдены
только численно.
В работах [4 – 6] для расчета первых десяти сейш в Черном и
Азовском морях применен метод конечных элементов, который по-
зволяет с высокой степенью точности аппроксимировать береговую
черту. Для решения задачи о сейшах в Черном море использована
сетка, включающая 448 линейных конечных элементов с 260 узлами.
При нахождении сейш в Азовском море использована конечно-
элементная сетка из 284 треугольных элементов со 180 узлами.
Применение метода конечных элементов позволило свести краевую
задачу (4.1), (4.3) к обобщенной матричной задаче на собственные
значения, решение которой найдено с помощью преобразования Ха-
усхолдера и TQL-алгоритма. В результате были
найдены периоды
первых десяти баротропных сейш и соответствующие им распреде-
ления смещений свободной поверхности жидкости. Частично эти
результаты представлены на рис. 4.1 и 4.2. Аналогичный метод при-
менен в работе [7] для расчета сейш в Черном море.
В работе [8] для расчета сейш в Азовском море с учетом вра-
щения Земли и квадратичного донного трения
применен явный по
времени метод конечных разностей на разнесенных для смещений
уровня и проекций скорости течения прямоугольных сетках с шага-
ми по обеим осям 1,5 км (161
×250 узлов). В этом случае граница
расчетной области представляет собой ступенчатую функцию с па-
раллельными координатным осям звеньями и поэтому на границе
выполняется равенство нулю одной из проекций горизонтальной
скорости. Реализован резонансный метод выделения сейш. Он опи-
рается на исходную систему уравнений длинных волн. Суть метода
следующая. На входе в Керченский
залив задавались гармонические