Научная статья на тему 'ВЕРТИКАЛЬНЫЙ ПЕРЕНОС ИМПУЛЬСА ВНУТРЕННИМИ ВОЛНАМИ В ЗАПАДНОЙ ЧАСТИ СРЕДИЗЕМНОГО МОРЯ'

ВЕРТИКАЛЬНЫЙ ПЕРЕНОС ИМПУЛЬСА ВНУТРЕННИМИ ВОЛНАМИ В ЗАПАДНОЙ ЧАСТИ СРЕДИЗЕМНОГО МОРЯ Текст научной статьи по специальности «Физика»

CC BY
30
8
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
Ключевые слова
ВНУТРЕННИЕ ВОЛНЫ / ВОЛНОВОЙ ПОТОК ИМПУЛЬСА / СТОКСОВ ДРЕЙФ

Аннотация научной статьи по физике, автор научной работы — Слепышев А.А., Носова А.В.

Цель. Исследование вертикального переноса импульса внутренними волнами при учете вращения Земли и сдвигового течения в западной части Средиземного моря, изучение влияния турбулентной вязкости и диффузии на указанные волновые потоки и стоксов дрейф - цель данной работы. Методы и результаты. Для решения исходной системы нелинейных уравнений гидродинамики используется слабонелинейный подход. В первом порядке малости по амплитуде волны решается краевая задача для амплитуды вертикальной скорости, во втором порядке по амплитуде волны исследуются нелинейные эффекты: стоксов дрейф и вертикальные волновые потоки импульса. Указанная краевая задача решается двумя способами: методом возмущений, применявшимся ранее, и численным методом по неявной схеме Адамса третьего порядка точности. Метод возмущений предполагает разложение решения и частоты волны в ряд по малому параметру, пропорциональному горизонтальной турбулентной вязкости. В работе сравниваются результаты, полученные методом возмущений и численным методом. Показано совпадение результатов по расчету дисперсионных кривых обоими методами. Однако для декремента затухания волны метод возмущений дает завышенные значения, для второй моды бóльшие, чем для первой. Вертикальные волновые потоки импульса отличны от нуля, и метод возмущений дает завышенные значения для потока uw. Вертикальный волновой поток импульса vw практически не зависит от турбулентной вязкости и диффузии, и оба метода дают для него идентичные результаты. Компонента скорости стоксова дрейфа вдоль направления распространения волны также нечувствительна к турбулентной вязкости и диффузии, в то время как поперечная составляющая равна нулю при отсутствии турбулентной вязкости и диффузии. Выводы. Метод возмущений дает завышенные значения декрементов затухания волны, волнового потока импульса uw и поперечной составляющей скорости стоксова дрейфа. На волновой поток импульса vw и продольную составляющую скорости стоксова дрейфа горизонтальная турбулентная вязкость и диффузия влияния практически не оказывают.

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

VERTICAL TRANSFER OF MOMENTUM BY INTERNAL WAVES IN THE WESTERN PART OF THE MEDITERRANEAN SEA

Urpose. The work is aimed at studying both the momentum vertical transfer by internal waves with the regard for the Earth rotation and the shear current in the western part of the Mediterranean Sea, and the influence of turbulent viscosity and diffusion upon the indicated wave fluxes and the Stokes drift. Methods and Results. To solve the initial system of the hydrodynamics nonlinear equations, a weakly nonlinear approach was used. In the first order of smallness in the wave amplitude, the boundary problem for the vertical velocity amplitude was solved; in the second order in the wave amplitude, the nonlinear effects, namely the Stokes drift and the vertical wave momentum fluxes, were investigated. The indicated boundary problem was solved in two ways: by the perturbation method applied earlier and by the numerical one by the implicit Adams scheme of the third order of accuracy. The perturbation method assumes expansion of the solution and the wave frequency in a series by a small parameter proportional to the horizontal turbulent viscosity. The results obtained by the perturbation and numerical methods were compared. Coincidence of the results of calculating the dispersion curves by both methods is shown. However, for the wave damping decrement, the perturbation method yields the overestimated values, at that for the second mode the values are higher than those for the first one. The vertical wave momentum fluxes are nonzero, and the perturbation method yields the overestimated values for the flux uw. The vertical wave momentum flux vw is practically independent of turbulent viscosity and diffusion, and both methods give the identical results for it. The velocity component of the Stokes drift along the wave propagation direction is also insensitive to turbulent viscosity and diffusion, whereas the transverse component equals zero in the absence of turbulent viscosity and diffusion. Conclusions. The perturbation method provides the overestimated values of the wave damping decrements, the wave momentum flux uw and the transverse component of the Stokes drift velocity. The horizontal turbulent viscosity and diffusion exert practically no effect upon the wave momentum flux vw and the longitudinal component of the Stokes drift velocity.

Текст научной работы на тему «ВЕРТИКАЛЬНЫЙ ПЕРЕНОС ИМПУЛЬСА ВНУТРЕННИМИ ВОЛНАМИ В ЗАПАДНОЙ ЧАСТИ СРЕДИЗЕМНОГО МОРЯ»

УДК 551.466.8

БО!: 10.22449/0233-7584-2022-4-358-371

Вертикальный перенос импульса внутренними волнами в западной части Средиземного моря

А. А. Слепышев А. В. Носова

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

Поступила в редакцию 10.03.2022; одобрена после рецензирования 13.04.2022; принята к публикации 04.05.2022.

Аннотация

Цель. Исследование вертикального переноса импульса внутренними волнами при учете вращения Земли и сдвигового течения в западной части Средиземного моря, изучение влияния турбулентной вязкости и диффузии на указанные волновые потоки и стоксов дрейф - цель данной работы.

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

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

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

Ключевые слова: внутренние волны, волновой поток импульса, стоксов дрейф

Благодарности: работа выполнена в рамках государственного задания по теме № 0555-20210004.

Для цитирования: Слепышев А. А., Носова А. В. Вертикальный перенос импульса внутренними волнами в западной части Средиземного моря // Морской гидрофизический журнал. 2022. Т. 38, № 4. С. 358-371. (1о1:10.22449/0233-7584-2022-4-358-371

© Слепышев А. А., Носова А. В., 2022

Vertical Transfer of Momentum by Internal Waves in the Western Part of the Mediterranean Sea

A. A. Slepyshev A. V. Nosova

Marine Hydrophysical Institute, Russian Academy of Sciences, Sevastopol, Russia

H slep55@mail.ru

Abstract

Purpose. The work is aimed at studying both the momentum vertical transfer by internal waves with the regard for the Earth rotation and the shear current in the western part of the Mediterranean Sea, and the influence of turbulent viscosity and diffusion upon the indicated wave fluxes and the Stokes drift. Methods and Results. To solve the initial system of the hydrodynamics nonlinear equations, a weakly nonlinear approach was used. In the first order of smallness in the wave amplitude, the boundary problem for the vertical velocity amplitude was solved; in the second order in the wave amplitude, the nonlinear effects, namely the Stokes drift and the vertical wave momentum fluxes, were investigated. The indicated boundary problem was solved in two ways: by the perturbation method applied earlier and by the numerical one by the implicit Adams scheme of the third order of accuracy. The perturbation method assumes expansion of the solution and the wave frequency in a series by a small parameter proportional to the horizontal turbulent viscosity. The results obtained by the perturbation and numerical methods were compared. Coincidence of the results of calculating the dispersion curves by both methods is shown. However, for the wave damping decrement, the perturbation method yields the overestimated values, at that for the second mode the values are higher than those for the first one. The vertical wave momentum fluxes are nonzero, and the perturbation method yields the overestimated values for the flux uw . The vertical wave momentum flux vw is practically independent of turbulent viscosity and diffusion, and both methods give the identical results for it. The velocity component of the Stokes drift along the wave propagation direction is also insensitive to turbulent viscosity and diffusion, whereas the transverse component equals zero in the absence of turbulent viscosity and diffusion.

Conclusions. The perturbation method provides the overestimated values of the wave damping decrements, the wave momentum flux uw and the transverse component of the Stokes drift velocity. The horizontal turbulent viscosity and diffusion exert practically no effect upon the wave momentum flux vw and the longitudinal component of the Stokes drift velocity.

Keywords: internal waves, wave momentum flux, Stokes drift

Acknowledgements: the study was carried out within the framework of the state assignment on theme No. 0555-2021-0004.

For citation: Slepyshev, A.A. and Nosova, A.V., 2022. Vertical Transfer of Momentum by Internal Waves in the Western Part of the Mediterranean Sea. Physical Oceanography, 29(4), pp. 334-346. doi:10.22449/1573-160X-2022-4-334-346

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

сти. Таким образом, мелкомасштабная турбулентность обусловливает вертикальный обмен в океане [3]. Турбулентность порождается при обрушении поверхностных и внутренних волн [3, 4]. Гидродинамическая неустойчивость течений и внутренних волн также приводит к генерации турбулентности [2, 5, 6]. Нередко гидродинамическая неустойчивость течений порождает внутренние волны, которые тоже становятся неустойчивыми и генерируют более мелкомасштабные волны [7]. Следовательно, происходит каскадный перенос энергии в малые масштабы, вплоть до турбулентности. Обрушение внутренних волн отнюдь не типичный процесс в океане, гораздо чаще турбулентность поддерживается слабым сдвигом скорости течения во внутренней волне 1.

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

Турбулентная вязкость и диффузия позволяют описать воздействие мелкомасштабной турбулентности на внутренние волны, которые при учете турбулентной вязкости и диффузии затухают [8-10]. Вертикальные волновые потоки импульса при этом отличны от нуля [11, 12]. Однако эти потоки находились при решении краевой задачи для амплитуды вертикальной скорости методом возмущений путем разложения решения и частоты волны в ряды по малому параметру, пропорциональному горизонтальной турбулентной вязкости. В настоящей работе эта краевая задача решается численно по неявной схеме Адамса третьего порядка точности. Представляет интерес сравнить волновые потоки импульса, полученные численным методом и методом возмущений. Это же касается и скорости стоксова дрейфа частиц жидкости. Вертикальные волновые потоки импульса у инерционно-гравитационных внутренних волн отличны от нуля и при отсутствии турбулентной вязкости и диффузии при наличии течения, у которого компонента скорости, нормальная к направлению распространения волны, зависит от вертикальной координаты [13, 14]. В этой связи актуально исследовать влияние горизонтальной турбулентной вязкости и диффузии на указанные процессы.

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

1 Островский Л. А., Соустова И. А., Цимринг Л. Ш. Воздействие внутренних волн на мелкомасштабную турбулентность в океане. Н. Новгород, 1981. 15 с. (Препринт № 31 / ИПФ АН СССР). 360 МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ том 38 № 4 2022

Уравнения движения для волновых возмущений имеют вид:

ди , „ ди ди ди _ dUn

(и + и0)--ь V--ь м--/V + м-

д дх ду дг dz

дР „ д2и ^ д2и К-- + К-

0

р0(0)дх дх2 д

ду , „ д д д _ дР „д V ^д V ^

--ь (и + и0)--ъу--ь м--ь /и =---ь К—- + К—(2)

дt дх ду дг р0 (0)ду дх ду

дм т т „ дм дм дм --ь ( + и0)--ь V--ь м-=

дt дх ду дг

дР ^ д^м ^ д^м р

=--+ К-— + К--- g—-—

Р0(0)д2 дх2 д2 р0(0)

(3)

др ггчдр др др dр0 д2р _ _ д2р

— + (и + и0)— + V— + м— = + М—^ + М—у, (4)

дt дх ду дг dz дх ду

ди ду дм ^ —= 0, (5)

дх ду дг

где х, у, г - декартовые координаты, ось г направлена противоположно ускорению свободного падения g; и, V, м - проекции волновой скорости течения в указанной системе координат; р0 - невозмущенная средняя плотность воды; Р и р - волновые возмущения давления и плотности; / - параметр Кориолиса; К, М - коэффициенты горизонтальной турбулентной вязкости и диффузии, которые полагаются постоянными; и0 (г) - скорость среднего течения. Сделаем оценку горизонтального масштаба изменения средней плотности. Для

этого используем геострофические соотношения 2. Из них следует оценка для

р g

указанного масштаба Ьр =--0—¡- = —:-¡-. Профиль вертикального гра-

|др0 / Щ /|ди0 / дг| диента скорости среднего течения приведен на рис. 1, а. Максимальное значение модуля вертикального градиента скорости течения составляет 2,2 ■ 10-2 ■ 1/с, тогда Ьр - не менее 5,19 х 106 м, т. е. много больше длины волны. Таким образом, зависимостью средней плотности от горизонтальной координаты можно пренебречь.

В качестве граничных условий используем условия «твердой крышки» на поверхности и отсутствие тангенциальных напряжений 3 [9]:

„ ,^ дм дм „

м = 0, К-= 0, К— = 0, г = 0, (6)

дх ду

2 Каменкович В. М. Основы динамики океана. Ленинград : Гидрометеоиздат, 1973. С. 128.

3 Миропольский Ю. З. Динамика внутренних гравитационных волн в океане. Ленинград : Гидро-

метеоиздат, 1981. С. 30.

МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ том 38 № 4 2022 361

на дне - условие непротекания и также отсутствие тангенциальных напряжений 3 [9]:

dw 8w

w = 0, K— = 0 , K— = 0, z = -H, 8x 8y

(7)

где H - глубина моря.

Р и с. 1. Вертикальные профили градиента скорости течения (a), скорости течения (b), частоты Брента - Вяйсяля (c)

F i g. 1. Vertical profiles of the current velocity gradient (a), the current velocity (b) and the Brunt -Vaisala frequency (c)

У дна тангенциальные напряжения нулевые, так как коэффициентами вертикального обмена пренебрегаем по сравнению с коэффициентами горизонтального турбулентного обмена.

Линейное приближение. Решения линейного приближения ищем в следующем виде [11-14]:

u = u10 (z)A exp(i(kx - roí)) + c.c., v = v10 (z)A exp(i(kx - roí)) + c.c., wl = w10 (z)Aexp(i(kx - roí)) + c.c., P = po(z)Aexp(i(kx - roí)) + c.c., Pj = p10(z)Aexp(i(kx - roí)) + c.c., (8)

где A - амплитудный множитель; k - горизонтальное волновое число; ю - частота волны; c.c. - комплексно сопряженные слагаемые.

После подстановки (8) в систему (1) - (5) следует связь амплитудных функЦий ui0, P0, р10, vi0 с Wjq [П, 12]:

i dw,.

Ро(0) k

k dz

Q dwin dün

Q = ro-k• U0,

(9)

-win --

if 2

Pió =-

k dz dz 10 k(iQ-k2K) dz dPo „ _ if

dWio + ikKdWio

dz

w

dw,.

iQ-k M dz

k (iQ- k2K) dz

(10)

Функция w10 удовлетворяет следующему уравнению [11, 12]:

(Q + ik2K) ÍQ2 + 2iQk2K - f2 - k4K21 ^ - kf2 dUó ^

+

+k

2

d U

- kQ - ik3K (Q + ik2K)2 + kN''

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

(Q + ik2 K )2 (Q + ik 2M)

(11)

wió = 0

где N2 =--—--— - квадрат частоты Брента - Вяйсяля.

Ро(0) &

Граничные условия имеют вид

^о(О) = м>ю(-И) = 0.

(12)

Остальные граничные условия в (6), (7) выполняются автоматически.

Краевая задача (11), (12) в [11, 12] решалась методом возмущений, когда решение и частота волны разлагались в ряд по малому параметру, пропорциональному значению горизонтальной турбулентной вязкости. В настоящей работе эта задача решается численно по неявной схеме Адамса третьего порядка точности. Сравнение дисперсионных кривых и зависимости декремента затухания волны от частоты для первых двух мод приводится ниже.

Нелинейные эффекты. Связь между скоростью течения в представлении Эйлера и и Лагранжа иь с точностью до членов второго порядка по крутизне волны имеет следующий вид 4 [15]:

иь = и+(|иьахУ)и .

(13)

Данное интегральное уравнение решается методом итераций. С точностью до членов, квадратичных по амплитуде волны, после осреднения по периоду волны получается выражение для средней лагранжевой скорости 4:

иь = и +(| иахУ )и

(14)

где вектор и(ио, V ) - скорость среднего течения (V = 0); и - поле волновых

эйлеровых скоростей, черта сверху означает осреднение по периоду волны. Скорость стоксова дрейфа частиц жидкости представлена вторым слагаемым в (14) и определяется по следующей формуле 4 [15]:

I

и = (|и^У)и .

(15)

Компоненты скорости стоксова дрейфа вдоль и поперек направления распространения волны имеют вид [11, 14]

и =

А А

1 А.

ю dz

щ

dz

+ с.с.

(16)

. * г d V = А А--

5 11 1

ю dz

(о* - гк2 К Д к dz

/dw1

* ^ 10

+ с.с.,

(17)

где А = Аехр(5ю • () , 5ю = 1т(ю) - мнимая часть частоты, декремент затухания волны.

Из (8) - (10) следуют выражения для вертикальных волновых потоков импульса ищ, ^ [11, 14]:

ищ =

т\А

* Л

— щ10 ^ =

dz dz

/аЧо

(о + гк2Кк dz

+ с.с.

(18) (19)

о

о

о

к

4 Дворянинов Г. С. Эффекты волн в пограничных слоях атмосферы и океана. Киев : Наукова думка, 1982. С. 25.

364 МОРЖОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ том 38 № 4 2022

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

Результаты расчетов. Гибралтарский пролив соединяет Средиземное море с Атлантическим океаном. Воды Атлантического океана, менее соленые, с высокой скоростью входят в Гибралтарский пролив, располагаясь в верхнем 100-метровом слое. Глубже находятся более соленые воды, которые более медленно вытекают из Средиземного моря. На выходе из Гибралтарского пролива в Средиземное море формируется слой с резким перепадом скорости течения и плотности воды, соответствующие вертикальные профили скорости течения и стратификации показаны на рис. 1, Ь, с [16]. Полусуточный прилив, набегая из Атлантического океана на Гибралтарский пролив, генерирует внутренние волны. В частности, в ходе натурного эксперимента были обнаружены мощные волновые пакеты низшей моды внутренних волн с периодом 14 мин [16]. С целью выяснения вертикальной структуры внутренних волн краевая задача (11), (12) решается численно по неявной схеме Адамса третьего порядка точности при К=2М. На рассматриваемых масштабах типичное значение коэффициента горизонтальной турбулентной диффузии составляет М=1 м2/с [17]. Для фиксированной частоты волны волновое число и декремент затухания волны находятся методом пристрелки. Результаты численных расчетов сравниваются с результатами решения краевой задачи (11), (12) методом возмущений [12], когда решение и частота волны ищутся в виде ряда по малому параметру, пропорциональному горизонтальной турбулентной вязкости.

Из рис. 2 следует, что дисперсионные кривые, построенные методом возмущений и путем численного решения краевой задачи (11), (12), фактически совпадают. На рис. 3 показана зависимость мнимой части частоты волны от действительной части частоты для первой (1) и второй (2) мод, рассчитанных по методу возмущений (кривые 1в, 2в) и путем численного решения краевой задачи (11), (12) (кривые 1ч, 2ч). Метод возмущений дает завышенные значения модуля мнимой части частоты волны, для второй моды - большие.

Нормирующий множитель А находится по известной величине максимальной амплитуды вертикальных смещений [12]. Вертикальные волновые потоки импульса ищ (формула (1 8)) для 1 4-минутных внутренних волн первой моды с амплитудой 16 м, полученные как при численном решении краевой задачи (11), (12), так и при применении метода возмущений, показаны на рис. 4. Метод возмущений дает завышенные значения вертикального потока импульса ищ .

Р и с. 2. Дисперсионные кривые первой (1) и второй (2) мод внутренних волн F i g. 2. Dispersion curves of the first (1) and second (2) modes of internal waves

Р и с. 3. Зависимость декремента затухания от частоты волны, полученная методом возмущений (кривые 1в, 2в) и численным решением краевой задачи (кривые 1ч, 2ч) F i g. 3. Dependence of the damping decrement upon the wave frequency obtained by the perturbation method (curves 1p, 2p) and by the numerical solution of the boundary problem (curves 1n, 2n)

Р и с. 4. Профили вертикального волнового Р и с. 5. Профили вертикального волнового

потока импульса vw

F i g. 5. Profiles of the vertical wave momentum flux vw

потока импульса uw , полученные методом возмущений (сплошная) и численным методом (штриховая)

F i g. 4. Profiles of the vertical wave flux of momentum uw obtained by the perturbation method (solid line) and by the numerical one (dashed line)

Вертикальный волновой поток импульса ^ (формула (19)) показан на рис. 5. Метод возмущений и численный метод дают идентичные результаты, более того, на волновой поток ^ практически не влияют горизонтальная турбулентная вязкость и диффузия.

Волновые вертикальные потоки импульса ищ, ^ и соответствующий турбулентный поток и'щ' сравниваются на рис. 6. Турбулентный поток импульса определяется по градиентной гипотезе: и'щ' = -К~~° . Коэффициент

dz

вертикального обмена в верхнем 150-метровом сильно стратифицированном слое обратно пропорционален частоте Брента - Вяйсяля: Кг = 8,4^ • 10-4 м2/с, здесь N соответствует частоте Брента - Вяйсяля (цикл/ч) [18]. Глубже, где стратификация слабая, коэффициент турбулентного обмена пропорционален N [19]. В целом волновые потоки импульса по модулю превосходят турбулентный. Волновой поток ^ по абсолютной величине превосходит поток им/.

Рис. 6. Профили волновых (сплошная и

штриховая) и турбулентного u W' (пунктир) вертикальных потоков импульса F i g. 6. Profiles of the wave (solid and dashed

lines) and turbulent u W' (dotted line) vertical momentum fluxes

Р и с. 7. Вертикальные профили поперечной составляющей скорости стоксова дрейфа, полученные методом возмущений (сплошная) и численным решением краевой задачи (11), (12) (штриховая) F i g. 7. Vertical profiles of transverse component of the Stokes drift velocity obtained by the perturbation method (solid line) and by the numerical solution of the boundary problem (11), (12) (dashed line)

Вертикальные профили поперечной к направлению распространения волны составляющей скорости стоксова дрейфа V, рассчитанные методом возмущений и численным методом, показаны на рис. 7. Метод возмущений

дает завышенные значения в верхнем 70-метровом слое. При отсутствии турбулентной вязкости и диффузии V = 0 .

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

Р и с. 8. Вертикальные профили продольной компоненты скорости стоксова дрейфа F i g. 8. Vertical profiles of the longitudinal component of the Stokes drift velocity

Р и с. 9. Профили вертикального волнового

потока импульса vw при увеличенных коэффициентах турбулентного обмена, полученные методом возмущений (черная кривая), численным методом (красная кривая) и в невязком случае (зеленая кривая) F i g. 9. Profiles of the vertical wave momentum flux vw at the increased turbulent exchange coefficients obtained by the perturbation method (black curve), by the numerical one (red curve) and in the inviscid case (green curve)

Представляет все-таки интерес исследовать влияние турбулентной вязкости и диффузии на волновой поток импульса и на скорость стоксова дрейфа и при больших коэффициентах турбулентного обмена. На рис. 9 показаны профили вертикального волнового потока импульса ™, рассчитанного по методу возмущений и численным методом при М = 50 м2/с, а также для невязкого случая при М = К = 0 . Метод возмущений дает сильные искажения вертикального потока импульса в верхнем 70-метровом слое. Волновой

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

На рис. 10 показаны вертикальные профили скорости стоксова дрейфа ы! 14-минутных внутренних волн первой моды, рассчитанные методом возмущений, численным методом и для невязкого случая при М = К = 0 . Метод возмущений дает заметное искажение скорости стоксова дрейфа в верхнем 25-метровом слое.

Р и с. 10. Вертикальные профили скорости стоксова дрейфа us при увеличенных коэффициентах турбулентного обмена, полученные методом возмущений (черная кривая), численным методом (красная кривая) и в невязком случае (зеленая кривая)

F i g. 10. Vertical profiles of the Stokes drift velocity us at the increased turbulent exchange coefficients obtained by the perturbation method (black curve), by the numerical one (red line) and in the inviscid case (green curve)

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

Выводы

1. Вертикальный волновой поток импульса uw при решении краевой задачи (11), (12) методом возмущений превышает поток uw при решении этой задачи численным методом по неявной схеме Адамса третьего порядка точности.

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

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

несколько увеличивает волновой поток импульса у^ и скорость стоксова дрейфа и5.

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

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

1. Булатов В. В., Владимиров Ю. В. Волны в стратифицированных средах. Москва : Наука, 2015. 735 с.

2. Коротаев Г. К., ПантелеевН. А. Экспериментальные исследования гидродинамической неустойчивости в океане // Океанология. 1977. Т. XVII, вып. 6. С. 941-953.

3. Wunsch C., Ferrari R. Vertical Mixing, Energy, and the General Circulation of the Ocean // Annual Review of Fluid Mechanics. 2004. Vol. 36. P. 281-314. doi: 10.1146/annurev.fluid.36.050802.122121

4. Самодуров А. С., Любицкий А. А., Пантелеев Н. А. Вклад опрокидывающихся внутренних волн в структурообразование, диссипацию энергии и вертикальную диффузию в океане // Морской гидрофизический журнал. 1994. № 3. С. 14-27.

5. Коротаев Г. К., Пантелеев Н. А. Гидродинамическая неустойчивость внутренних волн в океане при нестационарном сдвиге // Известия Академии наук СССР. Физика атмосферы и океана. 1977. Т. 13, № 10. С. 1044-1054.

6. Подымов О. И., Зацепин А. Г., Островский А. Г. Вертикальный турбулентный обмен в черноморском пикноклине и его связь с динамикой вод // Океанология. 2017. Т. 57, № 4. С. 546-559. https://doi.org/10.7868/S0030157417040049

7. Пантелеев Н. А., Щербаков А. Н., Щербакова Е. Н. Наблюдение гидродинамической неустойчивости в океане // Известия Академии наук СССР. Физика атмосферы и океана. 1989. Т. 25, № 6. С. 616-626.

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

8. Ле Блон П., Майсек Л. Волны в океане. Москва : Мир, 1981. Ч. 2. 363 с.

9. LeBlond P. H. On the damping of internal gravity waves in a continuously stratified ocean // Journal of Fluid Mechanics. 1966. Vol. 25, iss. 1. Р. 121-142. doi: 10.1017/S0022112066000089

10. Островский Л. А., Соустова И. А. Верхний перемешанный слой океана как сток энергии внутренних волн // Океанология. 1979. Т. XIX, вып. 6. С. 973-981.

11. Слепышев А. А. Вертикальный перенос импульса внутренними волнами при учете турбулентной вязкости и диффузии // Известия Российской академии наук. Физика атмосферы и океана. 2016. Т. 52, № 3. С. 342-350. doi:10.7868/S0002351516030111

12. Слепышев А. А., Носова А. В. Генерация вертикальной тонкой структуры внутренними волнами при учете турбулентной вязкости и диффузии // Морской гидрофизический журнал. 2020. Т. 36, № 1. С. 5-19. doi:10.22449/0233-7584-2020-1-5-19

13. Анкудинов Н. О., Слепышев А. А. Вертикальный перенос импульса внутренними волнами в двумерном потоке // Известия Российской академии наук. Механика жидкости и газа. 2021. № 3. С. 39-47. doi:10.31857/S0568528121030026

14. Воротников Д. И., Слепышев А. А. Вертикальные потоки импульса, обусловленные слабонелинейными внутренними волнами на шельфе // Известия Российской академии наук. Механика жидкости и газа. 2018. № 1. С. 23-35. doi:10.7868/S0568528118010036

15. Longuet-Higgins M. S. On the transport of mass by time-varying ocean currents // Deep Sea Research and Oceanographic Abstracts. 1969. Vol. 16, iss. 5. P. 431-447. https://doi.org/10.1016/0011 -7471(69)90031 -X

16. Watson G. Internal Waves in a Stratified Shear Flow: The Strait of Gibraltar // Journal of Physical Oceanography. 1994. Vol. 24, iss. 2. P. 509-517. https://doi.org/10.1175/1520-0485(1994)024<0509:IWIASS>2.0.m;2

17. Боуден К. Физическая океанография прибрежных вод. Москва : Мир, 1988. 324 с.

18. Интенсификация вертикального турбулентного обмена в районах сопряжения шельфа и континентального склона в Черном море / В. А. Иванов [и др.] // Доповщ Нацюнально! академи наук Украши. 2008. № 6. С. 108-112.

19. Самодуров А. С. Взаимодополняемость различных подходов для оценки интенсивности вертикального турбулентного обмена в естественных стратифицированных бассейнах // Морской гидрофизический журнал. 2016. № 6. С. 37-48. doi:10.22449/0233-7584-2016-6-37-48

Об авторах:

Слепышев Александр Алексеевич, ведущий научный сотрудник, отдел турбулентности, ФГБУН ФИЦ МГИ (Россия, 299011, г. Севастополь, ул. Капитанская, д. 2), доктор физико-математических наук, ResearcherlD: V-6948-2017, ORCID ID: 0000-0002-92597558, slep55@mail.ru

Носова Анна Владимировна, ведущий инженер-исследователь, отдел турбулентности, ФГБУН ФИЦ МГИ (Россия, 299011, г. Севастополь, ул. Капитанская, д. 2)

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