Нелинейные волны и турбулентность в прибрежной зоне О. Коса Тузла

Термогидродинамика океана

УДК 551.465

В.З. Дыкман, В.А. Иванов, В.М. Кушнир

Нелинейные волны и турбулентность в прибрежной зоне о. Коса Тузла

Представлены результаты спектральной обработки измерений трех компонент пульсаций вектора скорости течения в диапазоне частот 0,1 ... 50 Гц для условий прибрежной зоны о. Коса Тузла (центральная часть Керченского пролива). Измерения выполнены на расстоянии 0,25 м от поверхности морского дна на глубине 1,5 м. Вычислено пространственное распределение характеристик поверхностных волн при помощи прикладного пакета SWAN, параметры которого подстроены таким образом, чтобы период волн соответствовал данным измерений. Анализ спектров пульсаций вектора скорости показал, что динамика указанной прибрежной зоны определяется нелинейными волнами конечной амплитуды (волны Стокса) и мелкомасштабной турбулентностью. Определены скорость придонного трения (динамическая скорость) - 0,018 м/с, турбулентное напряжение - 0,33 Н/м2 и коэффициент турбулентной вязкости (диффузии импульса) - 1,89-10-3 м2/с.

Ключевые слова: спектры пульсаций вектора скорости, нелинейные волны Стокса, мелкомасштабная турбулентность.

Введение

Динамические процессы в прибрежной зоне моря на глубинах до 15 ... 20 м относятся к числу наименее изученных и наиболее актуальных для практических приложений, что связано со значительным дефицитом экспериментальных исследований динамики этой зоны, с одной стороны, и с активизацией в ней литодинамических эффектов размыва, переноса и переотложения донного материала - с другой. Наглядным примером проявления таких лито-динамических эффектов является прибрежная зона о. Коса Тузла, расположенного в центральной части Керченского пролива.

Хорошо известно, что Керченский пролив играет исключительную роль в экономике Украины и России как важнейшая транспортная магистраль между портами Азовского и Черного морей с выходом в Мировой океан, зона рыболовства, рекреации, перспективный район для добычи энергетических ресурсов и строительных материалов. По этим и другим причинам проблема изменчивости рельефа дна под действием волн и течений приобрела здесь высокую актуальность, которая резко возросла после сооружения в 2003 г. дамбы от Таманского п-ова в направлении украинского о. Коса Тузла и искусственного углубления протока между оконечностью дамбы и юго-восточной частью этого острова. Эти изменения рельефа дна в значительной степени повлияли на естественное поле течений, следствием чего была интенсификация процессов размыва о. Коса Тузла. Так, на протяжении 2005 г. - первой половины 2006 г. площадь острова из-за отсутствия работ по его укреплению уменьшилась на 30 га. Огромные объемы донного материала, смываемые с © В.З. Дыкман, В.А. Иванов, В.М. Кушнир, 2012

ISSN 0233-7584. Мор. гидрофиз. журн., 2012, № 4

3

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

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

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

С целью получения данных прямых измерений характеристик динамики придонной области в прибрежной зоне моря в МГИ НАН Украины был создан комплекс аппаратуры «Донная станция», в состав которого входят трех-компонентные измерители пульсаций и среднего значения вектора скорости течения, ловушки взвешенного донного материала с дистанционным считыванием показаний, лазерный прозрачномер для определения концентрации взвеси и среднего размера частиц [3, 4]. Данные таких измерений комплекси-руются результатами численных расчетов параметров поверхностных волн на основе прикладного пакета SWAN. Динамика придонного слоя в значительной степени определяется этими процессами, которые, по существу, формируют придонные вихреволновые структуры и турбулентность. Цель настоящей работы - анализ данных «Донной станции» в комплексе с модельными расчетами поверхностных волн для определения основных механизмов и параметров придонной динамики, а также величин механических напряжений на поверхности дна, которые вызывают горизонтальный перенос и формирование вертикальных профилей донного материала.

Материалы и методы

Измерения пульсаций скорости течения были выполнены трехкомпо-нентным электромагнитным измерителем, установленным неподвижно в специальной раме на расстоянии 0,25 м от поверхности морского дна. На этой раме смонтирована аппаратура комплекса «Донная станция» (рис. 1), который был установлен в юго-восточной части о. Коса Тузла на глубине 1,5 м и на расстоянии ~10 м от береговой линии.

4

ISSN 0233-7584. Мор. гидрофиз. журн., 2012, № 4

Р и с. 1. Комплекс «Донная станция», расположенный на юго-восточной оконечности о. Коса Тузла: 1 - датчик пульсаций трех компонент вектора скорости течения u, v, w, температуры T и электропроводности С; 2 - измерители виброускорений по трем осям

Измерения выполнялись 22 - 25 августа 2008 г. и состояли из 32 серий продолжительностью 7 ... 24 мин. Преобладал устойчивый северо-восточный ветер со средней скоростью 7 м/с. Измерительные оси датчика скорости ox и oy (соответственно ы- и у-компоненты скорости течения) были установлены в горизонтальной плоскости, причем ось ох направлена вдоль меридиана на север, а третья измерительная ось oz (w -компонента) ориентирована вертикально вверх. Частота измерений 100 Гц, данные в виде физических величин пульсаций скорости течения (м/с) вводились непосредственно в память персонального компьютера. Так как характерный период поверхностных волн был равен ~2 с, продолжительность указанных серий обеспечила высокую статистическую достоверность спектральных характеристик основных волновых возмущений и индуцированной турбулентности в диапазоне частот 0,1 ... 50 Гц. Расчеты спектральных характеристик выполнены по хорошо из-

0233-7584. Мор. гидрофиз. журн., 2012, № 4

5

вестным алгоритмам путем Фурье-преобразования корреляционных функций [5]. Для обработки из общих массивов измеренных данных отбирались ряды по 10000 отсчетов, длина корреляционных функций составляла 600 отсчетов, что обеспечивало спектральное окно в указанном частотном диапазоне. Число степеней свободы при оценке всех спектров равно 33, что соответствует неравенству - 0,15 < 1о§(/ £) < 0,18, где Зв - действительное значение, S -выборочная оценка спектра. По оси ординат это совпадает с диапазоном выборочной изменчивости оценок спектров, равной 33% одного порядка (рис. 2 - 4).

Р и с. 2. Спектры горизонтальных и, V составляющих пульсаций скорости течения, квадрата когерентности Я2 и фазового сдвига Е

Основные особенности спектральных характеристик составляющих вектора скорости течения состоят в следующем (рис. 2, 3): пик спектральной плотности приходится на период 2 с, что характерно для поверхностных волн на небольших глубинах вблизи берега; резкое уменьшение когерентности соответствует частотам более 0,8 ... 1 Гц; фазовый сдвиг между вертикальной и горизонтальной составляющими вектора скорости в области высокой когерентности (0,4 ... 0,8 Гц) изменяется в пределах 71 ... 76°, что близко к квадратурной зависимости между горизонтальной и вертикальной волновыми скоростями; для диапазона частот от 0,8 . 1 до 50 Гц наблюдаются резкие, неупорядоченные колебания когерентности и фазового сдвига, характерные для мелкомасштабных турбулентных флуктуаций. Для этого диапазона частот спектры составляющих скорости изменяются как Би v f-3.

6

0233-7584. Мор. гидрофиз. журн., 2012, № 4

Р и с. 3. Спектры горизонтальной u и вертикальной ч составляющих пульсаций скорости течения, квадрата когерентности Я2 и фазового сдвига Е

Р и с. 4. Спектры собственных значений тензора спектральной плотности Х1, Х2 и угла Ф между осью Х\ и плоскостью меридиана

0233-7584. Мор. гидрофиз. журн, 2012, № 4

7

В общем случае спектральные характеристики пульсаций вектора скорости описываются тензором спектра Бу, который можно представить следующим образом [6]:

Sv =

V V V

ии ' иу? ин

V V V

уи' уу' ун

V V V

^ ЛАП1 ^ ^ ЛЛП У '

ну э нн

где ^и , ^, - автоспектры; Яиу = Сиу - 12иу - взаимный спектр; Сиу - ко-

спектр; Quу - квадратурный спектр.

Симметричная и кососимметричная части тензора Sv равны соответственно

V С С ии' иу^ ин О-Юи* -Юин

С = С V С иу? уу' ун II Ъиу А-Юун

С С V

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

Sv =

Suu , Сиу ^и Сиу + Юиу, ^

Sc =

V С

ии иу Сиу, Sуv

Sк =

0, ¡Яиу - Юу ,0

Фазовая траектория движения частиц воды на частоте / представляет собой замкнутую фигуру (окружность, эллипс или линию возвратно-поступательных колебаний), оси 1 и 1 которой являются собственными значениями симметричной части тензора Sv, а их значения инвариантны к ориентации системы отсчета. Пространственная ориентация большой оси 1 относительно оси ох определяется по соотношению

2С„„

р = агс^

причем

11,2 = 0,5[^и + Sуу ^ .

На рис. 4 представлены вычисленные спектры собственных значений (1,12) симметричной части тензора, а также фазового угла между ох и

0233-7584. Мор. гидрофиз. журн., 2012, № 4

8

большей осью эллипса Л1. Видно, что основное возмущение (период 2 с) проявляется в виде эллипса, большая ось которого направлена с северо-востока на юго-запад под углом 10,5° к плоскости меридиана, что соответствует среднему направлению ветра. Большая ось эллипса равна 0,05 (м/с)2/Гц, малая ось 0,02 (м/с)2/Гц. В области высоких частот, начиная с 0,8...1,0 Гц, фазовый угол изменяется хаотично, что указывает на неупорядоченные, некогерентные пульсационные движения, характерные для турбулентности.

К числу инвариантных параметров тензора спектра скорости относят индикатор вращения D = -2iQuv, характеризующий направление и интенсивность вращательных движений. Для основных возмущений с периодом 2 с действительное значение индикатора вращения по абсолютной величине равно 0,022 (м/с)2/Гц, что соответствует вращению вектора скорости по часовой стрелке, если смотреть на измерительные оси сверху вниз.

Оценки параметров поверхностных волн были выполнены по модели SWAN. Поле скорости ветра вычислено на основе натурных данных, усвоенных в модели MM-5, расчеты проводились методом вложенных сеток [7]. Основные результаты численных расчетов показаны на рис. 5. Параметр модели SWAN (скорость диссипации энергии обрушения волн) был подобран таким образом, чтобы вычисленный период поверхностных волн Tw на восточной оконечности о. Коса Тузла был близок к 2 с, что соответствует данным измерений пульсаций скорости течения на «Донной станции». При этом длина волны Л равна ~ 5,8 м, высота волн с обеспеченностью 3% составляет 0,6 ... 0,65 м, средняя высота и амплитуда волн равны 0,44 и 0,22 м. Отношения глубины моря к длине волны h / Л и средней высоты волны к ее длине H/Л равны 0,2586 и 0,07586 соответственно. Из диаграммы применимости различных моделей поверхностных волн (рис. 6) видно, что полученные значения h/Л и H/Л относятся к области нелинейных волн Стокса, профиль которых описывается соотношением

h

n=1

где Fn(h/Л) - параметр профиля волны, причем F1 = a , F2 = a2F22 + a4F24,

F3 = a3F33 + a3F35, F4 = a4F44 , F5 = a5F55, a - параметр высоты волны, который определяется из решения нелинейного уравнения

0,5kH = a + a3F33 + a5

Для рассматриваемых условий значение параметра a получено методом последовательных приближений и равно 0,232.

ISSN 0233-7584. Мор. гидрофиз. журн., 2012, № 4

9

Р и с. 5. Результаты расчетов поверхностных волн: а - направление распространения и высоты волн 3-процентной обеспеченности (м) по данным измерений 23.08.2008 г. в 15:00; б - длины волн (м); в - периоды волн (с); г - скорость и направление ветра. Параметры волн показаны на цветных шкалах

Я /А 0.12

0.08 -

0.04 -

2

1 ♦

3

у

о 0.2 о а о.б д.

Р и с. 6. Диаграмма применимости модели кноидальных волн (1), нелинейных волн Стокса (2) и линейных волн малой амплитуды (3). ИИ — отношение высоты волны к ее длине, НИ — отношение глубины к длине волны [8, 9]. Значком показано положение параметров поверхностных волн в прибрежной зоне о. Коса Тузла

10 ISSN 0233-7584. Мор. гидрофиз. журн., 2012, № 4

Горизонтальная ух и вертикальная уу составляющие скорости волнового течения вычисляются по следующим соотношениям [8, 9]:

5 со$Ь

со

= Т А °п~

к

$,тЬ

со8[п(кх - «)],

« Сп ^*У)81П[п(кх-«)],

к 8тЬ

где G1 = aG11 + а^13 + а5015, G2 = 2(а2G22 + a4G24), G3 = 3(a3G33 + a5G35), G4 = 4a4G44, G5 = 5a5G55 . Зависимости ¥ц (к /1) = /(к /1) для волн первых

пяти порядков приведены в табл. 1.

Т а б л и ц а 1

Зависимости параметров формы волны Стокса от г = к /1 при достоверности аппроксимации Я2

Зависимость ^ = /

Диапазон

Я1

^22 = 0,2 г-^24 = 0,505г ^33 = 0,077 г ^35 = 0,469г

-0,51 -1,87 -1,26

^44 = 0,09г Г55 = 0,062г-

-1,70

0,10 < г < 0,6 0,15 < г < 0,6 0,10 < г < 0,6 0,15 < г < 0,6 0,10 < г < 0,6 0,10 < г < 0,6

0,854 0,950 0,835 0,877 0,822 0,822

х

п=1

у

У

Параметры скорости волны зависят от г = к /1, соотношения для их расчета приведены в табл. 2. Вертикальная координата У в соотношениях (8) отсчитывается от поверхности дна.

Дисперсионное соотношение для нелинейных волн Стокса записывается

как

« = gk

где С1 = /1(к /1); С2 = /2(к /1) - параметры частоты волны и давления из [8, 9]. Использование для анализа модели нелинейных волн Стокса в данном случае обеспечивает наиболее точное воспроизведение реальных характеристик поверхностных волн.

Как уже отмечалось, высокочастотная часть спектра пульсаций скорости течения (начиная с частоты 0,8 ... 1 Гц) имеет структуру, характерную для мелкомасштабной турбулентности. Известно [10], что колмогоровские микромасштабы длины, скорости и времени при этом равны соответственно

0233-7584. Мор. гидрофиз. журн., 2012, № 4

11

¡к = К0'75^) , ^ =

где е0 - скорость диссипации энергии турбулентности:

= 2к|к 25У (к )дк,

К - молекулярная вязкость (к@ 1,0 10-6 м2/с), к - волновое число. Оценка е0 выполнена по известному соотношению А. Бенилова для развитого волнения и при кк > 1, что соответствует условиям проведения измерений [11]:

У™ Нч3г , 6я?Ч1

11 1

где у @ 1 - постоянная; УК - фазовая скорость; г - глубина установки датчика. Для рассматриваемых условий е0 = 3,17 10-5 м2/с3, что совпадает с диапазоном известных оценок этого параметра [11].

Т а б л и ц а 2

Зависимости параметров скорости волны Стокса от г = к /1 при достоверности аппроксимации Я2 в диапазоне 0,1 < г < 0,6

Зависимость О = / (г = к /1) Я2

о„ = 1 1,000

О13 = -369,23г4 + 614,11г3 - 384,57г2 + 105,01г -11,29 0,996

О15 = 3621,3г5 - 7535,6г4 + 6167,1г3 - 249г2 + 499,2г - 41 0,950

О22 = 0,0004г_4,02 0,982

О24 = -36640г6 + 84299г5 - 78843г4 + 38349г3 -10224г2 + +1414,6г - 78,54 0,999

О33 = -1085,2г5 + 2154,6г4 -1658,7г3 + 617,46г2 - 110,81г + 7,63 0,994

О35 = 0,0135г "2'637 0,934

= 187,5г5 - 393,02г4 + 323,52г3 -130,8г2 + 25,99г - 2,0292 0,999

О55 = -6626,3г6 +15178г5 - 14118г4 + 6820,9г3 -1803,9г2 + + 247,4г -13,73 0,999

0

Микромасштаб турбулентности ¡к равен 0,43 мм на частоте 5,6 Гц. Граница равновесного интервала турбулентности определена по положению безразмерного параметра к1к , соответствующего максимуму спектра диссипации турбулентной энергии. По данным [12], этот параметр приближенно ра-

12

0233-7584. Мор. гидрофиз. журн., 2012, № 4

вен 0,1, что совпадает с граничным масштабом равновесного интервала @ 2,5 см.

Результаты

По результатам расчетов волнового поля вблизи о. Коса Тузла по модели SWAN получены следующие данные: период волны Tw = 2 с, длина волны l = 5,8 м, средняя высота и амплитуда H = 0,44 м и aw = 0,22 м соответственно. Амплитуды горизонтальной Vx и вертикальной Vy составляющих волновой скорости на поверхности моря равны:

2pH cosh(kh) . 2pH

Vx =--= 0,74 м/с, Vy =-= 0,69 м/с.

x Tw sinh(kh) y Tw

Для расчета аналогичных параметров по модели нелинейных волн Сто-кса пятого порядка использован предварительно определенный параметр высоты волны а = 0,232 .

Амплитуды горизонтальной У8х и вертикальной У8у составляющих волновой скорости течения для первой гармоники волны Стокса вычислены на основе соотношений (8):

тл W ^ 4coshkh

Vsx =- (aG„ + a3Gi3 + a5^)—— = 0,69 м/с, k sinh kh

VSy = W(aG11 + a3G13 + a5G15) = 0,64 м/с.

k

Относительные расхождения значений амплитуды горизонтальной составляющей скорости по кинематическим соотношениям и по модели волн Стокса пятого порядка составляют 6,8 и 6,9% соответственно для горизонтального и вертикального компонентов скорости волнового течения.

Расчет указанных амплитуд составляющих скорости по модели волн Стокса пятого порядка на глубине 1,25 м выполнен по соотношениям (8). Результаты для первых трех гармоник таких волн приведены в табл. 3.

Т а б л и ц а 3

Параметры трех первых гармоник волн Стокса пятого порядка

Номер гармоники Период, с Частота, Гц VSx ,м/с VSy , М/С S, (м/с)2 sSy , (М/С)2

1 2,000 0,5 0,270 0,0714 0,0363 0,00253

2 1,000 1,0 4,53 -10-3 2,2 -10-3 1,01 -10-5 2,4 -10-6

3 0,667 1,5 1,59 -10-5 1,07 -10-5 1,25 -10-10 5,68 -10-11

Спектральные компоненты первых трех гармоник показаны на рис. 7. Видно, что первая гармоника совпадает с измеренным спектром пульсаций скорости, однако, начиная с частоты 0,8 ... 1 Гц, наблюдается значительное

ISSN 0233-7584. Мор. гидрофиз. журн., 2012, № 4

13

расхождение по экспериментальным данным между моделью волн Стокса и спектром скорости. В диапазоне частот 0,8 ... 1 - 10 ... 15 Гц спектр пульсаций пропорционален /_зд , что, как уже было отмечено, соответствует турбулентности.

Р и с. 7. Спектральные характеристики горизонтальных составляющих пульсаций скорости Би + Б„, первых трех гармоник нелинейных волн Стокса (красные кружки) и придонной турбулентности /_3,1

Определение скорости придонного трения и* (динамическая скорость на дне) выполнено путем суммирования коспектров Си № (/) пульсаций скорости и, ^ в указанном диапазоне частот равновесного интервала придонной турбулентности 0,9 < / < 5,55 Гц:

и* —

(/ )Д /

/т

где интервал дискретности 0,01 с; длина корреляционной функции 600 точек; Д/ —1/12 Гц; и* — 0,018 м/с, что является типичным значением для аналогичных условий. Турбулентное напряжение т( на поверхности дна равно ри*2 = 0,33 Н/м2.

Представляет интерес сопоставление этой оценки на основе данных прямых измерений с результатами известных модельных оценок для волнового придонного пограничного слоя [13]. Известны различные модели типа (и*/ито) — ^(рг, / —1,2...), где р - параметры придонного пограничного слоя, в том числе величины, характеризующие неровность морского дна;

14

1ББМ 0233-7584. Мор. гидрофиз. журн., 2012, № 4

ито - амплитуда волновой скорости вблизи поверхности дна. Рассмотрим наиболее известные модели такого типа. 1. В работе [14]

и„

= 0,5 ),

= ехр[-5,977 + 5,213(а„/ён)-0,194 при 1,47 <а„/^ < 300,

= 0,32 при 1,47 > ак /< ,

где ён = 4ё - масштаб неровности дна или эквивалент шероховатости Нику-радзе; ё - средний размер фракций донного материала; ак = Ткито /2р -амплитуда горизонтальных волновых колебаний вблизи поверхности дна. В рассматриваемом случае ито = 0,26 м/с, а№ = 0,083 м, и* = 0,0165 м/с.

2. В работе [15]

к

и.

а

1и

где к = 0,41 - постоянная Кармана; и* = 0,0184 м/с. 3. В работе [16]

и* г\ 1 п л/Н \-0,182

и

= 0,104(—)-

где и* = 0,0084 м/с. 4. В работе [17]

к

и„

Р _2_ ^ + 5)

4 V Яе ё

где Яе =

Р

Н

5. В работе [18]

и* = 0,023 м/с.

0,5 к

и

[кег2(2^) + кш2^)]0,5

где V = ён /1 = 4ё /1; I = ю*/ю - масштаб высоты неровностей дна; кег( г), ке1( г) - функции Кельвина нулевого порядка. После аппроксимации

#0(г) = [кег2(г) + ке12(г)]0,5 = Аг, А = 0,5355, и = 0,2594, 5,5 < г < 0,55,

0233-7584. Мор. гидрофиз. журн., 2012, № 4

15

и

и

*

н

и

*

и

*

= (^V^)1^n){Owyп/(и+1), u* = 0,041 м/с.

A

mo

6. В работе [19]

где и* = 0,0108 м/с.

^^ = 0,122

Um„ H d

Среднее значение динамической скорости по этим оценкам составляет < и* >= 0,0196 м/с, что примерно на 8,5% отличается от данных прямых измерений. Наиболее близкие значения к данным по экспериментальным спектрам дают модели из работ [14, 15], при этом относительное расхождение не превышает 3%.

Полученное значение скорости трения и* использовано для определения возможности применения модели переменного логарифмического пограничного слоя для рассматриваемых условий. При этом, учитывая, что характерное время формирования волнового ламинарного и турбулентного пограничных слоев (толщины 0,8 мм и 3,5 см) значительно (по крайней мере в 10 и более раз) меньше периода волны, орбитальную волновую скорость можно рассматривать как «медленное» внешнее течение. Известно [17], что профиль амплитуды средней скорости в таком слое равен

и

п

при и* = 0,018 м/с, у = 0,25 м, п = 1,0 • 10-6 м2/с, и(0,25) = 0,468 м/с. С другой стороны, по данным вычисленного экспериментального спектра горизонтальных составляющих скорости течения (рис. 7) дисперсия пульсаций скорости равна 0,096 (м/с)2, стандартное отклонение оу = 0,31 м/с, а амплитуда

скорости течения -/2 • 0,31 = 0,438 м/с, что достаточно хорошо соответствует модельной оценке.

Согласно [17], коэффициент турбулентной вязкости (диффузии импульса) в логарифмическом пограничном слое имеет вид

u

KT = ки*y ,

где к = 0,42; y - расстояние от поверхности дна. При и* = 0,018 м/с y = = 0,25 м, KT = 1,89 10-3 м2/с, или 18,9 см2/с, что соответствует оценкам аналогичных параметров, полученным различными методами [20].

Подвижность различных фракций донного материала существенно зависит от их размеров [7, 15]. Поэтому гранулометрический состав играет важную роль в динамике изменений рельефа дна под действием волн и течений. В настоящей работе использованы осредненные данные инженерных геологических изысканий в нескольких зонах Керченского пролива [7]. В табл. 4 16 ISSN 0233-7584. Мор. гидрофиз. журн., 2012, № 4

представлено распределение основных фракций донного материала по их размерам d, а также других параметров. В их числе безразмерный параметр D* = d(ggV"2)0'33; число Шильдса (Shields numbers) Y = u*2/gdg, g = (ps - pw)/ pw; критическое число Шильдса Ycr, которое определяется по функциональной зависимости Ycr = f (D) [7, 15]; критическая скорость трения u*cr = (Ycrgdf)0'5. Также определена безразмерная концентрация взвеси

на дне Cb =

Ps -Pw Pb -Pw

= 0,035(lg

0

u*D*

)3,75, где ps, кг/м3 - плотность взвеси в

воде, pw - плотность воды, pb @ 2600 кг/м3 - плотность донного грунта.

Т а б л и ц а 4

Процентное содержание фракций донного материала по их размерам и некоторые другие параметры

Средний размер d, мм Процентное содержание Безразмерный параметр D. Гидравлическая крупность взвешенных частиц WS, м/с Гидравлическая крупность частиц на поверхности дна №ь, м/с Показатель степени m = = Ws /kV* Число Шильдса Y-103 Критическое число Шильдса Ycr-103 Критическая скорость трения u*cr, м/с

3,5 12,7 67,95 0,18 0,26 23,78 5,89 44,0 1,11

1,5 25,8 29,12 0,14 0,17 18,66 13,8 34,5 0,82

0,75 12,5 14,56 0,087 0,096 11,5 27,5 30,6 0,62

0,375 13,8 7,28 0,046 0,052 6,11 55,0 39,0 0,39

0,175 31,2 3,40 0,015 0,017 2,0 118 70,6 0,20

0,075 4,1 1,46 0,0028 0,0036 0,39 275 164 0,082

Решение уравнения для концентрации взвеси на расстоянии г от поверхности дна в логарифмическом придонном пограничном слое можно записать следующим образом [7]:

C (z) = Cb (-)-Ws/ku* = Cb (-)-z0 Zr>

0

где г0 = dY /6 - минимальный размер неоднородностей на дне, равный для

всех фракций 3,44 • 10-6 м.

Безразмерная концентрация взвеси на расстоянии 0,25 м от поверхности дна для всех фракций донного материала, за исключением самых мелких частиц (d = 0,075 мм), является исчезающе малой величиной. Для указанных

мелких частиц С(0,25) = 8,44 10-6 или в размерном виде Ср = 13,5 г/м3. При средней скорости течения иг = 0,1 ... 0,15 м/с удельный перенос взвеси составляет 1,69 г/м2-с.

Амплитуда горизонтальной скорости волнового течения основной гармоники равна 0,28 м/с (табл. 3) при толщине слоя вертикальных колеба-

ISSN 0233-7584. Мор. гидрофиз. журн., 2012, № 4

17

m

ний 0,18 м. Интегральный перенос взвеси в этом слое вычисляется по соотношению

С, = Сь = )-т _ ^)-т ],

• 20 1 - т 20 20

где г1 = 0,16 м, г2 = 0,34 м. Удельный перенос мелкодисперсной взвеси в слое установки ловушки для рассматриваемых условий равен 1,54-10-6 или в размерном виде

2,95-10-4 кг/с •м (1,84-10" м /м-с). Учитывая, что размер отверстия Ц вертикальной ловушки равен 0,016 м, средняя скорость горизонтального переноса взвеси и = 0,12 м/с, скорость вертикального опускания частицы Ж, = 0,0028 м/с, объем взвеси QV, попадающей в ловушку из указанного слоя, составляет

QV = СЦ

или 1,48-10-6 м3 за 6 ч экспозиции ловушки. При площади ловушки 1,256-10-3 м2 толщина слоя взвеси за указанное время составит ~ 1 ... 1,2 мм.

Другой, независимый, метод оценки толщины ку слоя взвеси основан на соотношении [15]

дК =_Э^ ,

Э^ Эх

где q = дт ът(ю1 - кх) = 1,84-10-7 м3/с-м - удельный перенос взвеси в слое 0,16 - 0,34 м. При отрицательных значениях Эq / Эх толщина слоя взвеси увеличивается за счет осаждения частиц грунта, при положительных - уменьшается за счет увеличения интенсивности размыва. С учетом экранирующего действия ловушки при интенсификации размыва, т. е. в моменты времени, когда Эq / Эх > 0, использовано среднее отрицательное значение Эq / Эх .

к 0,5т"

Опуская вычисление интеграла-— - кх)Л, находим:

0,5Тм> 0

АИУ = АТ = 6,34 • 10-8 АТ ,

р

где АТ - время экспозиции. При АТ = 6 ч АИУ @ 1,4 мм, что совпадает с приведенной выше оценкой.

Как известно [7, 15], перенос донного материала осуществляется в двух формах: в виде взвеси, соотношение (24), и в виде движения частиц непосредственно на поверхности дна. Первая форма переноса рассмотрена выше, для второй формы использованы основные соотношения модели Ван Райна ^ап Щп). В соответствии с этой моделью удельный перенос донного материала зависит от соотношения эффективной и критической скоростей трения, т. е. от отношения и*^ /н*сг. В рассматриваемом случае н*^ для всех фрак-

18

0233-7584. Мор. гидрофиз. журн., 2012, № 4

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

Выводы

1. Спектральная обработка данных измерений пульсаций скорости течения в диапазоне частот 0,1 ... 50 Гц показала, что при средней скорости устойчивого ветра 7 м/с в прибрежной зоне о. Коса Тузла (глубина 1,5 м) формируется сложная динамика поля скорости течения. Ее основными составляющими являются нелинейные поверхностные волны Стокса и мелкомасштабная турбулентность, индуцированная этими волнами. Согласно расчетам параметров нелинейных волн, выполненным по модели пятого порядка, амплитуда первой гармоники волновой скорости по данным измерений соответствует модельным расчетам с относительным расхождением не более 7%. Амплитуды второй и третьей гармоник волновой скорости на расстоянии от дна 0,25 м значительно меньше амплитуд турбулентных пульсаций на этих же частотах.

2. В диапазоне частот от 0,8 ... 1,0 до 12 ... 15 Гц спектры пульсаций скорости изменяются пропорционально f ~3,1, что по характеру изменчивости спектров когерентности и фазы соответствует режиму мелкомасштабной турбулентности. Скорость диссипации турбулентной энергии при использовании известного соотношения А. Бенилова равна 3,17 10-5 м2/с3. Микромасштаб турбулентности и соответствующая граничная частота fr составляют 0,42 мм и 5,5 Гц.

3. По данным экспериментальных спектров пульсаций скорости течения определены динамическая скорость (скорость трения) на дне (0,018 м/с) и турбулентные напряжения (0,33 Н/м2), под действием которых происходит перемещение и взмучивание донного материала - формирование профилей взвеси. Обычно аналогичные оценки находят на основе полуэмпирических моделей. Выполненные расчеты для шести моделей показали, что указанные параметры имеют такой же порядок величины, как и вычисленные на основе экспериментальных спектров. Относительное отклонение среднего значения по результатам расчетов по шести моделям и полученной оценки не превышает 8,5%, для двух моделей - 3%.

4. При u* = 0,018 м/с, y = 0,25 м и n = 1,0 • 10-6 м2/с модельная оценка амплитуды скорости течения для логарифмического пограничного слоя u(0,25) = 0,468 м/с согласуется с аналогичной оценкой по данным экспериментального спектра горизонтальных пульсаций вектора скорости u = 0,438 м/с. Коэффициент турбулентной вязкости (диффузии импульса) равен 18,9 см2/с, что в целом совпадает с оценками аналогичного параметра, полученными различными методами.

5. Анализ транспорта донного материала при рассмотренных условиях показал, что подвижной является самая мелкая фракция со средним размером 0,075 мм. Под действием турбулентных напряжений эти частицы образуют профили взвеси, концентрация которой на расстоянии 0,25 м от поверхности

ISSN 0233-7584. Мор. гидрофиз. журн., 2012, № 4 19

дна равна 13,5 г/м3. Двумя независимыми методами определена толщина слоя взвеси в ловушках за период экспозиции 6 ч. Она равна 1 ... 1,5 мм, что соответствует фактическим данным.

Авторы выражают благодарность С.В. Федорову за выполненные расчеты характеристик поверхностных волн на основе прикладного пакета SWAN и О.И. Ефремову - за полезное обсуждение и рекомендации при подготовке статьи.

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

1. Спиридонова Е.О. Оценка влияния гидрометеорологических факторов на распределение донных осадков в Керченской бухте // Системы контроля окружающей среды. -Севастополь: МГИ НАН Украины, 2007. - С. 235 - 239.

2. Ломакин П.Д., Спиридонова Е.О., Чепыженко А.И. и др. Антропогенные и природные источники взвешенного вещества в водах Керченского пролива // Морской экологический журнал. - 2008. - 7, № 4. - С. 51 - 59.

3. Иванов В.А., Дыкман В.З., Ефремов О.И. и др. Современные методы и средства контроля морской среды / Под ред. В.А. Иванова. - Севастополь: ЭКОСИ-Гидрофизика, 2006. - 112 с.

4. Дыкман В.З., Маньковский В.И., Маньковская Е.В. Прозрачномер для определения параметров крупной взвеси в море методом измерения флуктуаций показателя ослабления направленного света // Экологическая безопасность прибрежной и шельфовой зон и комплексное использование ресурсов шельфа. - Севастополь: ЭКОСИ-Гидрофизика, 2008. - Вып. 16. - С. 236 - 242.

5. Коняев К.В. Спектральный анализ случайных океанологических полей. - Л.: Гидроме-теоиздат, 1981. - 207 с.

6. Белышев А.П., Клеванцов Ю.П., Рожков В.А. Вероятностный анализ морских течений. -Л.: Гидрометеоиздат, 1983. - 264 с.

7. Иванов В.А., Кушнир В.М., Федоров С.В. Динамика изменений рельефа дна в Керченском проливе по данным космических съемок оптическими сканерами // Экологическая безопасность прибрежной и шельфовой зон и комплексное использование ресурсов шельфа. - Севастополь: ЭКОСИ-Гидрофизика, 2010. - Вып. 22. - С. 127 - 155.

8. Доусон Т. Проектирование сооружений морского шельфа. - Л.: Судостроение, 1986. -286 с.

9. Skjelbreia L., Hendrickson A. Fifth order gravity wave theory // Proc. VII Conference on Coastal Engineering, 1961. - Р. 184 - 196.

10. Филлипс О.М. Динамика верхнего слоя океана. - Л.: Гидрометеоиздат, 1980. - 319 с.

11. Монин А.С., Озмидов Р.В. Океанская турбулентность. - Л.: Гидрометеоиздат, 1981. -320 с.

12. Монин А. С., Яглом А.М. Статистическая гидромеханика. - М.: Наука, 1967. - 720 с.

13. Кушнир В.М. Структура волнового пограничного слоя над ровным и неровным дном // Морской гидрофизический журнал. - 2005. - № 6. - С. 188 - 202.

14. Bijker E. W. Longshore current computation // Proc. ASCE. J. Waterway. Harborsand Coastal Engineering Division, 1971.

15. Иванов В.А., МихиновА.Е. Прогноз динамики наносов в прибрежной зоне моря (Практические рекомендации и примеры расчетов). - Севастополь, 1991. - 50 с. (Препринт / АН УССР. МГИ).

16. Faria A.F., Thornton E.B., Stanton T.P. et al. Vertical profiles of longshore currents and related bed shear stress and bottom roughness // J. Geophys. Res. - 1998. - 103, № C2. -Р. 3217 - 3232.

20 ISSN 0233-7584. Мор. гидрофиз. журн., 2012, № 4

17. Прандтль Л. Гидроаэромеханика. - M.: Иностранная литература, 1949. - 520 с.

18. Grant W.D., Madsen O.S. Combined wave and current interaction with a rough bottom // J. Geophys. Res. - 1979. - 84, № C4. - Р. 1797 - 1808.

19. Ninomiya K., Tagaya K., Murase Y. A study on suction and securing of sit-on-bottom type offshore structure // OTC. - 1972. - Р. 873 - 886.

20. Кушнир В.М. Придонный пограничный слой в Черном море: экспериментальные данные, турбулентная диффузия, потоки // Океанология. - 2007. - 47, № 1. - С. 39 - 48.

Морской гидрофизический институт НАН Украины, Материал поступил

Севастополь в редакцию 21.02.11

E-mail: vdykman@alpha.mhi.iuf.net, После доработки 18.04.11 vaivanov@alpha.mhi. iuf. net, kushnirv@yahoo. com

АНОТАЦ1Я Представлеш результати спектральноТ обробки вишрювань трьох компонент пульсацш вектора швидкосп течГТ в дiапазонi частот 0,1.. .50 Гц для умов прибережноТ зони о. Коса Тузла (центральна частина КерченськоТ протоки). Вишрювання виконаш на вщсташ 0,25 м вщ поверхш морського дна на глибиш 1,5 м. Обчислено просторовий розподш характеристик поверхневих хвиль за допомогою прикладного пакету SWAN, параметри якого шдстро-еш так, щоб перюд хвиль ввдповвдав даним вишрювань. Аналiз спектрiв пульсацш вектора швидкосп показав, що динамiка вказаноТ прибережноТ зони визначаеться нелшшними хвилями кiнцевоТ амплгтуди (хвилi Стокса) i дрiбномасштабною турбулентшстю. Визначенi швидкiсть придонного тертя (динамiчна швидкiсть) - 0,018 м/с, турбулентна напруга - 0,33 Н/м2 i коефь цiент турбулентно!' в'язкостi (дифузп iмпульсу) - 1,89-10-3 м2/с.

Ключовi слова: спектри пульсацiй вектора швидкосп, нелiнiйнi хвилi Стокса, дрiбномас-штабна турбулентшсть.

ABSTRACT Results of spectral processing of measurements of three components of fluctuations of the current velocity vector in the frequency band 0.1 ... 50 Hz are represented for the conditions of the coastal zone of the island Kosa Tuzla (the central part of the Kerch strait). The measurements are carried out at 0.25 m distance from the bottom on the depth 1.5 m. Spatial distribution of the surface waves' characteristics are calculated due to application of the software SWAN whose parameters are adjusted in such a way that the waves' period corresponds to the measurement data. Analysis of the spectra of current vector fluctuations shows that dynamics of the indicated coastal zone is conditioned by nonlinear waves of final amplitude (the Stokes waves) and small-scale turbulence. Defined are the velocity of near-bottom friction (dynamic speed) - 0.018 m/s, turbulent stresses - 0.33 N/m2 and the turbulent viscosity coefficient (diffusion of momentum) - 1.89 • 10-3 m2/s.

Keywords: spectra of velocity vector fluctuations, nonlinear Stokes waves, small-scale turbulence.

ISSN 0233-7584. Мор. гидрофиз. журн., 2012, № 4

21