CC BY
16
Анализ результатов наблюдений
и методы расчета гидрофизических полей океана
УДК 532.59
В.А. Иванов, Л.В. Черкесов, Т.Я. Шульга
Экстремальные отклонения уровня и скорости течений, вызываемые постоянным ветром в Азовском море
Методом численного моделирования исследовано влияние скорости постоянного ветра на течения и колебания уровня Азовского моря. С использованием трехмерной нелинейной численной модели установлены районы, где имеют место максимальные и минимальные величины скоростей поверхностных стационарных течений и наибольшие отклонения уровня в открытой части Азовского моря. На приведенных картах показаны пункты акватории, в которых достигаются экстремальные значения отклонений уровня и наибольшие скорости течений. Выполнен анализ влияния скорости постоянного ветра на время достижения экстремальных значений подъемов и понижений уровня и скоростей течений.
Ключевые слова: численное моделирование, о-координатная модель, сгонно-нагонные процессы, поверхностные течения, экстремальные отклонения уровня.
Введение
Изучение колебаний уровня моря связано с научными и практическими задачами, включая строительство береговых гидротехнических сооружений, создание морских навигационных карт, прогноз наводнений, определение движения берегов. При решении некоторых прикладных проблем возникает необходимость определения положения экстремальных отклонений уровня по всей акватории Азовского моря. Представляет интерес установление тех районов моря, в которых скорости стационарных ветровых течений будут наибольшими.
В ряде работ, посвященных изучению течений и волнения в Азовском море, исследовался вопрос о максимальных сгонах и нагонах, имеющих место на ряде береговых станций, а также о наибольших значениях скоростей возникающих течений. В статье [1] с использованием трехмерной нелинейной о-координатной математической модели решена задача о сгонно-нагон-ных колебаниях в Азовском море при воздействии постоянного ветра, в [2] аналогичные исследования проведены для полей нестационарного ветра, в [3] — для возмущений, вызванных прохождением циклонов.
В данной работе на основании анализа результатов численных экспериментов, проведенных с использованием указанной выше модели, установлены районы расположения наибольших изменений уровня моря и скоростей стационарных течений, вызванных постоянным ветром в открытой части Азовского моря.
© В.А. Иванов, Л.В. Черкесов, Т.Я. Шульга, 2011
24
0233-7584. Мор. гидрофиз. журн, 2011, № 2
Постановка задачи
Граничные и начальные условия. Введем систему координат, в которой ось х ориентирована на восток, у - на север, г - от невозмущенной поверхности вертикально вверх. Исходные уравнения представляют собой нелинейные уравнения движения однородной несжимаемой жидкости в приближении теории мелкой воды [4, 5]:
du 1 Эр Э . . Э . . Э ^ Эи
/У + = (Тп) + (Т12) + км —, (1)
dt р Эх Эх Эу Эг Эг
dv 1 Эр Э . . Э . . Э ^ ЭУ
— + /и + --£ = _ (Т>21) + ~ (^22) + ^ Км— , (2)
dt р Эу Эх Эу Эг Эг
+ ЕР = 0, (3)
Эг
Эи Эу дм ^
— + — + — = 0. (4)
Эх Эу Эг
Здесь все обозначения общепринятые.
Компоненты тензора турбулентных напряжений Рейнольдса Т^ в уравнениях (1), (2) определяются через градиенты осредненных скоростей по формулам
Эи . ( Эу Эи ^ „ . Эу
Т11 = 2 Ам -ч , Т12 = Т21
Эх
-+ -
Т22 = 2 Ам— . (5)
22 м
Эх Эу) Эу
Коэффициент горизонтальной турбулентной вязкости Ам вычисляется с помощью формулы Смагоринского [6]. Для параметризации коэффициента вертикальной турбулентной вязкости Км использовалась теория Мелло-ра - Ямады [7].
Граничные условия на свободной поверхности имеют вид
Э^ Э^ Э^
— + и— + у—, Эt Эх Эу
Мг=с = ^г~ + + , (6)
(Эи Эу^
КмI Э' Эг.
= (Т0х,Т0у). (7)
г=?
Здесь Г0х = СкЖх^ и Т0у = СкЖу^ - проекции касательных напряжений трения ветра, Жх, Жу - компоненты вектора скорости ветра W на высоте 10 м над уровнем моря. Для численных расчетов коэффициент поверхностного трения С4 выбран в соответствии с работой [8].
0233-7584. Мор. гидрофиз. журн., 2011, № 2
25
На боковых границах выполняется условие прилипания. На дне (г = — H у )) нормальная составляющая скорости равна нулю:
( dH dH ^
w + u--+ v-
ÖX dy ;
= 0. (8)
z=-H
Придонные касательные напряжения трения связаны со скоростью квадратичной зависимостью [4]:
K fdu dv^
Km l äz' dz.
t t), (9)
где Tlx = CbWu2 + v2 , t1y = CbWu2 + v2 . Здесь коэффициент донного трения C b находится по формуле Cb = k 2ln-2( z2/ z0), где z2 - шаг по вертикали в придонном слое; z0 = 0,003 м - параметр шероховатости, характеризующий гидродинамические свойства подстилающей донной поверхности.
В начальный момент времени t = 0 движение жидкости отсутствует, свободная поверхность горизонтальна:
u(x, y, z ,0) = 0, v(x, y, z,0) = 0, Z(x, y, z,0) = 0. (10)
Численная реализация модели и выбор параметров. В исходных уравнениях (1) - (4), граничных условиях (6) - (9) и начальных (10) осуществляется переход от координаты z к о-координате [4, 5]: x* = x, y* = y, о= (z - Z)/(H + Z), t* = t, где - 1 < s< 0.
В численных экспериментах используются равномерные шаги по горизонтальным координатам (Ax = Ay = 1350 м) и по о-координате. Количество расчетных уровней по вертикали равно 11. Уравнения решаются с шагом At =18 с для определения осредненных двумерных компонент скорости и уровня моря и с шагом 10At - для вычисления отклонений от найденных средних и вертикальной компонент скорости. Рельеф дна и конфигурация береговой линии Азовского моря в узлах сетки взяты с навигационных карт.
Установившиеся движения в Азовском море формируются при нормальных условиях погоды. Эти условия характеризуются тем, что скорость ветра, достигнув определенного распределения над морем, продолжительное время мало меняется. В численных экспериментах при моделировании стационарных течений поле ветра задается в виде однородного по пространству ветра, действующего над Азовским морем начиная с момента времени t = 0. Его скорость (Wst) первые 3 ч (0 < t < 3 ч) нарастает со временем по линейному закону, достигает своего наибольшего значения и далее не меняется. Выход течений на установившийся режим осуществляется под действием постоянного западного ветра, скорость которого изменяется в пределах 5 - 15 м/с (Wst1 = 5, Wst2 = 10, Wst3 = 15 м/с).
26 ISSN 0233-7584. Мор. гидрофиз. журн., 2011, № 2
z=-H
Момент времени выхода течений на установившейся режим (^ = ¿0) определяем из условия, что относительные изменения полной энергии, осреднен-ной по объему моря на смежных шагах по времени, не превосходят 1% [1].
Анализ результатов численных экспериментов
Анализ опубликованных метеорологических данных [9] позволил выделить для района Азовского моря характерные типы погоды. Антициклонический тип погоды возникает в теплый период (апрель - октябрь) и характеризуется ветрами со скоростью более 8 м/с, когда Азовское море оказывается в тыловой части антициклона, смещающегося с запада на восток. В это время преобладают восточные и северо-восточные ветры. Слабый циклонический тип погоды [10] возникает при прохождении над Азовским морем слабых атмосферных фронтов или периферий отдельных циклонов. Средняя повторяемость его 12% в холодный период и 22% в теплый. Этот тип погоды характеризуется слабыми и умеренными ветрами переменных направлений, облаками среднего и нижнего ярусов; иногда выпадают слабые осадки. Циклонический тип погоды возникает при прохождении над Черным и Азовским морями средиземноморских циклонов или циклонов, развивающихся на хорошо выраженных атмосферных фронтах севернее параллели 50° с. ш. Повторяемость его зимой и летом не превышает 15%. Для этого типа погоды характерны ветры со скоростью 6 - 14 м/с вначале восточных, а затем южных и западных направлений. Западные ветры обычно бывают порывистыми и шквалистыми. Моделирование проводилось для наиболее продолжительного слабого циклонического типа погоды.
На рис. 1 представлены пункты расположения максимальных отклонений уровня в открытой части Азовского моря, вызываемых действием постоянного ветра различных скоростей. Эти пункты на рисунках отмечались в течение всего расчетного периода (72 ч). Отсюда (рис. 1, а - в) видно, что расположения максимумов при различных скоростях действующего ветра имеют как общие черты (максимумы сосредоточены в восточной части акватории моря), так и различия. Отметим, что воздействие ветра больших скоростей приводит к смещению координат максимумов в сторону Таганрогского залива. Действие ветра со скоростью 15 м/с вызывает максимумы отклонений уровня и в западной части Азовского бассейна.
Р и с. 1. Положения максимумов отклонений уровня в Азовском море в различные моменты времени под действием западного ветра со скоростью 5 м/с (а); 10 м/с (б); 15 м/с (в)
0233-7584. Мор. гидрофиз. журн., 2011, № 2
27
На рис. 2 показаны положения минимальных отклонений уровня Азовского моря с дискретностью 1 ч, вызываемых действием постоянного ветра различных скоростей. Отсюда видно, что наибольшие понижения уровня в основном имеют место в западной части моря. При больших скоростях постоянного ветра минимумы отклонений уровня наблюдаются и в восточной части Азовского моря.
Р и с. 2. Положения минимумов отклонений уровня в Азовском море в различные моменты времени под действием западного ветра со скоростью 5 м/с (а); 10 м/с (б); 15 м/с (в)
В табл. 1 приведены максимальные и минимальные значения отклонений уровня моря, время их достижения и соответствующие им координаты в зависимости от скорости ветра в течение всего расчетного периода. Данные в табл. 1 являются обобщением результатов, представленных на рис. 1 и 2.
Т а б л и ц а 1
Максимальные и минимальные значения отклонений уровня моря (м), время их достижения (ч, мин) и соответствующие им координаты (км) для трех значений скорости ветра
Скорость Zmax, tmax, xmax, .Уmax, Zmin, xmin, .Утит
ветра, м/с м ч, мин км км м ч, мин км км
£ ii 5 0,20 17.10 366,51 219,61 -0,10 13.23 4,02 95,07
Wst2 = 10 0,69 18.46 367,86 218,27 -0,39 13.47 4,02 107,13
ii 5 1,65 19.48 343,69 224,96 -1,05 14.09 5,37 87,04
Из анализа приведенных данных видно, что координаты максимумов и минимумов С для рассмотренных скоростей ветра изменяются незначительно. Различия состоят в величинах экстремальных отклонений уровня и времени их достижения. Из сравнения значений максимумов уровня следует, что они возрастают в 3,5 и 8,3 раза для больших в 2 и 3 раза значений (0,2 м при 0,69 м при 'й2; 1,65 м при 'й3). Как видно из табл. 1, под действием ветра большей скорости достигаются и большие значения понижений уровня: -0,1 м при 'й1; -0,39 м при 'й2; -1,05 м при Анализ результатов расчетов показывает, что действие ветра со скоростями и (в 2 и 3 раза большими, чем приводит к увеличению Спт в 4 и 10 раз соответственно.
28
ISSN 0233-7584. Мор. гидрофиз. журн., 2011, № 2
Время достижения экстремальных отклонений уровня отмечается через 10 — 16 ч после начала действия ветра и зависит от его скорости. По данным, приведенным в табл. 1, при сравнении изменения величины скорости ветра и времени достижения наибольших и наименьших значений ^получено, что это время больше для больших скоростей ветра. Так, разница во времени достижения £тх при и составляет 2 ч 38 мин, а разница во времени достижения £шп при тех же скоростях ветра — 46 мин. Таким образом, стационарный ветер в 3 раза большей скорости дольше на 14% формирует максимальный подъем уровня и на 4% — минимальное понижение уровня моря.
Сравним величины максимальных и минимальных значений £в открытой части Азовского моря (табл. 1) с экстремальными значениями сгонов и нагонов, отмеченными в ряде пунктов его побережья [1]. Наибольшие нагоны при скоростях ветра 5, 10 и 15 м/с отмечаются в Таганрогском заливе и составляют 0,17; 0,60 и 1,42 м соответственно. При тех же скоростях ветра в Азовском море достигаются значения ¿тх, равные 0,20; 0,69 и 1,65 м. Таким образом, максимальные возвышения уровня в открытой части Азовского моря на 15, 13 и 14% превышают нагоны на побережье. Аналогично, сравнивая наибольшие значения сгонов (Геническ) и ¿шп, отметим, что при рассматриваемых скоростях ветра сгон и максимальное понижение уровня отличаются мало.
В настоящее время прогнозы течений практически не составляются, причиной тому является недостаточное количество наблюдений, которые позволили ли бы проверить и уточнить расчеты скоростей течений. Чтобы разрабатывать и использовать методы расчета и прогноза скоростей морских течений в оперативной практике, необходимо хотя бы в нескольких отдельных точках моря проводить регулярные наблюдения за скоростями течений.
На рис. 3, 4 показаны районы, где имеют место максимальные и нулевые величины скоростей поверхностных течений в открытой части Азовского моря, вызываемых действием постоянного ветра различных скоростей. Эти пункты на рисунках показаны с дискретностью 1 ч в течение всего расчетного периода (72 ч). Из рис. 3 видно, что при различных скоростях действующего ветра положение максимумов смещено в сторону Таганрогского залива. Как следует из рис. 4, положения нулевых значений скоростей течений сильно разбросаны по акватории.
Р и с. 3. Положения максимумов скоростей течений Азовского моря на горизонте 1 м в различные моменты времени под действием западного ветра со скоростью 5 м/с (а); 10 м/с (б); 15 м/с (в)
0233-7584. Мор. гидрофиз. журн., 2011, № 2
29
Р и с. 4. Положения минимумов скоростей течений Азовского моря на горизонте 1 м в различные моменты времени под действием западного ветра со скоростью 5 м/с (а); 10 м/с (б); 15 м/с (в)
Т а б л и ц а 2
Максимальные значения скоростей течений (м/с) на различных глубинах Азовского моря, время их достижения (ч, мин) и соответствующие им координаты (км) для трех значений скорости ветра
Скорость ветра, м/с
Горизонт, м
Um
max?
м/с
xm
км
ym
км
Wst1 = 5
1
3 5 10
0,16 0,13 0,12 0,08
8.00 8.02 8.22 9.44
236,29 236,39 237,01 235,34
174,08 173,44 173,31 173,00
1 0,49 8.35 236,29 172,74
Wst2 = 10 3 0,44 8.55 236,60 173,32
5 0,41 8.26 237,09 171,73
10 0,29 9.34 236,86 171,98
1 1,06 9.45 237,63 171,40
Wst3 = 15 3 0,97 10.47 237,99 172,28
5 0,91 11.19 237,76 171,88
10 0,64 12.45 238,46 171,50
tmax, ч, мин
В табл. 2 приведены максимальные значения скоростей течений на различных горизонтах, время их достижения и соответствующие им координаты в зависимости от скорости ветра в течение всего расчетного периода. Из анализа данных, представленных в этой таблице, видно, что при увеличении скорости постоянного ветра в 2 и 3 раза скорости течений в поверхностном слое достигают значений 0,16; 0,49 и 1,06 м/с (возрастают соответственно в 3 и 6,6 раза). Такое же соотношение между величинами скоростей течений для
30 ISSN 0233-7584. Мор. гидрофиз. журн., 2011, № 2
рассматриваемых Wst прослеживается и на других горизонтах. С глубиной происходит уменьшение скорости течений, ее наименьшее значение отмечается в придонном слое. Сравнивая координаты максимумов, отметим, что и с глубиной, и при разных значениях Wst их положение практически не меняется. Время достижения наибольших значений скоростей течений увеличивается (на 4 и 18%) под действием ветра больших скоростей (Wst2 и Wst3 соответственно). Оценивая разницу во времени достижения максимумов скорости, видим, что tmax под действием постоянного ветра рассматриваемых скоростей в поверхностном слое достигается на 18, 12 и 21% раньше, чем в придонном слое.
Выводы
В результате численного моделирования циркуляции вод Азовского моря, вызванной действием постоянного ветра различных скоростей, получены следующие выводы:
— установлены пункты в открытой части моря, в которых отмечаются максимальные значения отклонений уровня и наибольшие скорости течений;
— получены оценки влияния скорости действующего постоянного ветра на время достижения наибольших подъемов и понижений уровня и максимальных значений скоростей течений на различных горизонтах.
СПИСОК ЛИТЕРАТУРЫ
1. Иванов В.А., Фомин В.В., Черкесов Л.В., Шульга Т.Я. Исследование эволюции поля примеси в Азовском море при наличии стационарных течений // Доп. НАН Укра'ши. -2007. - № 7. - С. 116 - 120.
2. Иванов В.А., Фомин В.В., Черкесов Л.В., Шульга Т.Я. Исследование влияния течений, вызванных западным ветром, на эволюцию областей загрязнения в Азовском море // Там же. - 2007. - № 3. - С. 112 - 117.
3. Иванов В.А., Фомин В.В., Черкесов Л.В., Шульга Т.Я. Исследование сгонно-нагонных явлений в Азовском море, вызванных атмосферными возмущениями // Там же. - 2006. - № 11. - С. 109 - 113.
4. Blumberg A.F., Mellor G.L. A description of three dimensional coastal ocean circulation model // Three-Dimensional Coastal Ocean Models / Ed. N. Heaps. - Washington, D. C.: American Geophysical Union. - 1987. - 4. - P. 1 - 16.
5. Фомин В.В. Численная модель циркуляции вод Азовского моря // Научные труды Укр-НИГМИ. - 2002. - Вып. 249. - C. 246 - 255.
6. Smagorinsky J. General circulation experiments with primitive equations. I. The basic experiment // Mon. Wea. Rev. - 1963. - 91. - P. 99 - 164.
7. Mellor G.L., Yamada T. Development of a turbulence closure model for geophysical fluid problems // Rev. Geophys. Space Phys. - 1982. - 20. - P. 851 - 875.
8. Hsu S.A. A mechanism for the increase of wind stress coefficient with wind speed over water surface: A parametric model // J. Phys. Oceanogr. - 1986. - 16. - P. 144 - 150.
ISSN 0233-7584. Мор. гидрофиз. журн., 2011, № 2
31
9. Гидрометеорология и гидрохимия морей СССР. Т. V. Азовское море. - СПб.: Гидроме-теоиздат, 1991. - 234 с.
10. Лоция Азовского моря. - http://katamaran.rU/maps/azov/2.htm.
Морской гидрофизический институт НАН Украины,
Севастополь
E-mail: [email protected]
Материал поступил в редакцию 13.11.09 После доработки 27.04.10
АНОТАЦ1Я Методом чисельного моделювання дослщжено вплив швидкост постшного в^ру на течи та коливання рiвня Азовського моря. З використанням тривимiрноI нелшшно! чисель-но! моделi знайденi райони розташування максимальних та мiнiмальних значень швидкостей поверхневих стацiонарних течiй i найбшьших вiдхилень рiвня у вiдкритiй частиш Азовського моря. На наведених картах вiдмiченi пункти акваторiI, в яких досягаються екстремальнi зна-чення вщхилень рiвня та найбiльшi швидкостi течш. Зроблено аналiз впливу швидкостi постшного в^ру на час досягнення екстремальних значень тднять та опускань i швидкостей течiй.
Ключовi слова: чисельне моделювання, о-координатна модель, зпнно-напнш процеси, поверхневi течiI, екстремальнi вщхилення рiвня.
ABSTRACT Method of numerical simulation is used to study influence of constant wind velocity upon the currents and the sea level fluctuations in the Azov Sea. The three-dimensional non-linear model permits to define the regions where maximum and minimum values of surface stationary currents' velocities and locations of the largest level deviations in the open part of the Azov Sea take place. The represented charts show the points in the sea area where extremal values of the level deviations and the highest currents' velocities are achieved. Analyzed is the influence of constant wind velocity upon the time of achieving extremal values both of the level elevation and lowering, and the currents' velocities.
Keywords: numerical simulation, o-coordinate model, surge processes, surface currents, extremal level deviations.
32
ISSN 0233-7584. Мор. гидрофиз. журн., 2011, № 2
