Исследование влияния рельефа дна на ветровую циркуляцию в двухслойном бассейне Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

УДК 551.465

И.И. Карпатович

Исследование влияния рельефа дна на ветровую циркуляцию в двухслойном бассейне

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

Ключевые слова: ветровая циркуляция, влияние рельефа дна.

Введение. Данное исследование мотивировано недавними наблюдениями глубинной циркуляции Черного моря с использованием всплывающих буев [1]. Скорость наблюдаемых глубинных течений в Черном море в среднем была равна 1 - 3 см/с, иногда она доходила до 7 см/с. В ходе эксперимента было запущено три буя - на глубине пикноклина 200 м, в промежуточном слое 750 м и в глубинном слое 1550 м. Данные, полученные с первого буя, указывают на активную роль мезомасштабных вихрей в динамике верхнего слоя моря. Скорость течения здесь достигала 15 см/с вдоль оси Основного Черноморского течения. Данные с буев на глубине промежуточного слоя 750 м и горизонте 1550 м также указывают на влияние мезомасштабных вихрей на течения, хотя скорость здесь в основном не превышала 5 см/с. Кроме того, характерная траектория буя на глубине 1550 м до захвата его мезомас-штабным вихрем проходила вдоль изолиний функции / / Н, что указывает на существенное влияние рельефа дна на придонные течения. Наблюдаемая циркуляция во всей толще моря была циклонической без изменения направления у дна, как это имеет место в океане. Таким образом, проведенный эксперимент дает новое представление о глубинной циркуляции в Черном море, т. к. ранее предполагалось, что течения в глубинной части этого бассейна должны быть очень слабыми. Слабая циркуляция в глубинных слоях не вызывала сомнения, т. к. Черное море имеет резко выраженный мелкий пикнок-лин [2]. Следовательно, если аппроксимировать стратификацию двухслойной моделью, то течение в нижнем слое должно быть намного слабее, чем в верхнем. Учет влияния рельефа дна на циркуляцию исследован в задачах применительно к океаническим масштабам [3] - [5] и типам рельефа дна. Менее исследованы течения в замкнутых бассейнах, где влияние рельефа может приводить к качественным изменениям структуры глубинной циркуляции. Отметим в этой связи результаты работы [6], в которой проведено аналитическое исследование влияния меридионального наклона дна на глубинные течения. Было показано, что такие наклоны могут приводить к существенной перестройке течений в поверхностных и глубинных слоях, а также значительно интенсифицировать глубинные течения в замкнутом бассейне. © И.И. Карпатович, 2011

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

73

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

Уравнения движения. Рассмотрим стационарную ветровую циркуляцию в двухслойном бассейне с глубиной H(x,у), используя простейшую двухслойную линейную модель:

- У у,

- g

У и1 = - g

Э и Э у

+ —1

Э х Э У

Э x

Э у +

+ A1

Э2 и1

+ л

Э ж1

"э7

Э2 у 1 Э г2

- fV2

У и 2

- g

- g

Э х Э у

Э и2 Э х

+

Э У2 Э у

ЭС2

- g1чг +

Э х

Э у + Э ж

А Э 2 и2

- gl

+ л

Э 2 2

Э 2 У2 2 Э 2 2

Э 2

0.

(1) (2)

(3)

(4)

(5)

(6)

Здесь и1 2 , у 2 , м12 - зональные, меридиональные и вертикальные составляющие скорости течения в верхнем и нижнем слоях; и С2 — уровень моря и отклонение поверхности раздела слоев; У = У0 + ¡5 У — параметр Корио-

лиса; g1 = g —Р; —р — разница плотности слоев; Л1 , Л2 — коэффициенты

Р

вертикального турбулентного обмена в верхнем и нижнем слоях. Граничные условия для уравнений (1) — (6) следующие: при г = 0

. Э щ тх . Э у Ту Л —1 = —х , А —1 = —У , ж = 0 ;

1 Э.

р

1 Э:

р

(7)

при г = И0

Э щ

Э у

и1 - u2), А1= -кЛ. У - У2 ) , Э г

74

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

0

э

г

. Э и2 . Э и, . Э у2 . Э у,

А = А, , А = А, , (8)

Э 2 Э 2 Э 2 Э 2

ж, = = 0 ;

при 2 = Н

. Э и2 — Э у2 , А,—2 = - ^2 и2 , А,—2 = - ^2 У2 , (9)

Э 2 Э 2

Э Н Э Н

Э х Э у

Здесь Тх и Ту - составляющие тангенциального напряжения ветра; к, , к2 -

коэффициенты трения между слоями и о дно моря; и,, и2, у,, у2 - средние скорости течений в верхнем и нижнем слоях, определяемые выражениями

Ао п Н

_ __, "о _ __1

( ^ У1) = К I( U1, У1) ё2 , ( ^ У2) = н _ К I( У2) ё2 •

На поверхности моря и границе раздела слоев полагаем вертикальную скорость равной нулю, что допустимо ввиду линейности задачи и малости и по сравнению с толщинами слоев. Интегрируя уравнения (1) - (6) по вертикали в пределах каждого слоя, с учетом граничных условий (7) - (9) получим:

_/V = _ 8+ -Тт _ и2), (10)

Э х р К0 К0

/и, = -8+ Т - гУ,-У2), (11)

Э у рК К

Э и Э у, 1 + —1 = 0, (12)

Э х Э у

_ __к

Э х 1 Э х Н - К0 2 Н - К0

-/У2 =-8 ^ - 81 ^и2 + -к^( и, - и2), (13)

/и2 =-8 ^ - 81 ^-т^- у 2 + 77^ У, - У 2), (14)

Э у Э у Н - К0 Н - К

0

Э __э _

— и2 (Н-К0) + — уЛН-К,) = 0. (15)

Э х Э у

Стандартный вид уравнений неразрывности (12) и (15) позволяет ввести интегральные функции тока в верхнем и нижнем слоях по формулам

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

75

7 дщ — дщ V К0 = , щ К = , (16)

д х д у

Н-Ко) = Щ2 ( Н-Ко) = -дщ2. (17)

д х д у

Исключая перекрестным дифференцированием наклоны уровня и ¿Т2 в формулах (10), (11) и (13), (14), с учетом (16), (17) получим:

рЩ =^ - Кг АУ + А" У - 7ЙЛ2 М". Н), (18)

д х р К0 Н - К0 (Н -К0)2

(19)

ЬдУ2 -^1 (У2, Н)=-к+кт АУ2 + к щ +

д х Н - К0 Н - К0 К0

+ 2(к| + ^ 71 (У2, Н)--к-71(у1, Н).

(Н - К0)2 ' (Н - К,) К0 ^ ;

„ ^ дгу дтх т, , ч да дЬ да дЬ

Здесь го1^ т =—----- , I (а, Ь) =-----,

д х д у д х д у д у д х

д а д Ь д а д Ь

71 ( а, Ь ) =--+--.

д х д х д у д у

На границах бассейна поставим интегральное условие непротекания, которое сводится к равенству нулю функций тока в верхнем и нижнем слоях на границе области интегрирования. Будем интегрировать численно систему уравнений (18), (19) в прямоугольном бассейне (размером х0 по долготе и

у0 по широте) с котловинным профилем рельефа дна, качественно моделирующим Черное море.

Полагаем х0 = 1000 км, у0 = 300 км, К0 = 150 м, Н0 = 2 км,

Н = Н 0 + А Н 0(1-7-^42»-)71 - 7У-У//2)-),r, = 0, г, =г„соз р,

V4 у0/4 у у0

г0 = 0,008 Н/м2, к1 = 10-4 м/с, к2 = 2 -10-4 м/с, /0 = 10-4 с -1,

Ь = 1,6 -10-11 м-1 - с-1 , р = 103 кг - м-3 . В качестве масштаба длины выберем Ь = 20 км, близкий к радиусу деформации первой бароклинной моды Россби, масштаб скорости течений и = 1 см/с. Переходя в системе уравнений (18), (19) к безразмерным переменным, получим:

Ь0 ^= Т го! г т'- Г1 ( А^-АЩ /(1 + а Н') +

д х

(20)

а — /1(У ,Н')),

(1 + а Н')

76

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

А У -Ь+а1'Ун) =

Э х 1 + а Н = - (Г +12) АУ2/(1 + а Н) + /12 А У ■

'2

Г12 а

ТГ/1(у1, н') + 2Г+;,)а /1(У2, Н),

(1 + а Н) ' (1 + а Н')2

безразмерные параметры здесь следующие:

Т = т°т

о рКи2' _ к1

Т

А=Аи

а =

А Н

( Н о - Ко)

к 2т

К и'

(Но -Ко)и

к1 т

12

(Но -Ко)и

А

/о т а

и

а

АТ

/о

Основным параметром, определяющим режим циркуляции в нижнем слое бассейна, является а.

Метод интегрирования системы уравнений (2()), (21). Для интегрирования системы уравнений (2о), (21) введем в прямоугольной области равномерную сетку по координатам х и у с шагом К:

р(х, у) ®р(х,,у}), х, = 1к, у} = /К, К = хо/N=уо/М, о<,<N, о</<М. В каждом из слоев определим функции вихря:

¿1 = А у, (22)

г2 = А у2. (23)

Дифференциальные операторы заменим следующими конечно-разностными аппроксимациями, обеспечивающими второй порядок точности по пространственным координатам:

Эр Э х

Р 1 -р

Ар ®

Р

2 К

.+1, /

Эр

ЭУ

/+1 /-1 р -р

I ,-л,-1, / I г'л, , / + 1 I г'л, , /-1 + р +р + р

2 К -4 р^;

К

На границе прямоугольной области функция тока полагается равной нулю. Определим следующий двухшаговый итерационный процесс. На первом шаге из уравнений (2о), (21) находим функции вихря г1 и г2 во внутренних точках сетки. На втором шаге из уравнений Пуассона (22), (23) находим

/SSN 0233-7584. Мор. гидрофиз. журн., 2011, № 6

77

новое приближенное значение функции тока в каждом из слоев по методу верхней релаксации:

щ =щ +т(Щк+и +ук-1,у + щ,у+1 +щу-1 -4 Щ; -41 К2), к = 1,2.

Новые значения функций тока Щ 2 используем для вычисления новых приближений 2 из уравнений (20), (21) и т. д. Параметры релаксации ¡112 подбираются эмпирически так, чтобы была сходимость итераций. Вычисления заканчивались после выполнения в каждом узле сетки условия | у / у - 1 | < е, где е - малый параметр.

Система уравнений (20), (21) интегрировалась вышеуказанным методом при различных значениях параметра АН0: АН0 = 0 (плоское дно), АН0 = 100 м, АН0 = 150 м (рис. 1 - 3) с шагом 2 км на сетке 500 х 150 узлов. Параметры ¡и1 и ¡12 полагались равными 0,2, е = 10-7 .

Обсуждение результатов расчета скорости течений. Проведение трех численных экспериментов дает представление о влиянии рельефа дна на характер циркуляции в замкнутом бассейне. При плоском рельефе дна имеет место ярко выраженный пограничный слой у западного берега моря (рис. 1). Скорость течения наиболее интенсивна в западном пограничном слое (в верхнем слое ~ 51 см/с, в нижнем ~28 см/с), в юго-западной, северо-западной, южной и северной частях бассейна (в верхнем слое ~ 8 см/с, в нижнем ~1 см/с), менее интенсивна в центральной (в верхнем слое ~1 см/с, в нижнем ~ 0,2 см/с) и восточной (в верхнем слое ~ 3 см/с, в нижнем ~ 0,1 см/с) областях. При учете рельефа дна в обоих слоях происходит изменение режима циркуляции, причем в нижнем слое скорости течения становятся близкими к изолиниям функции / /Н (рис. 2, 3). При значении параметра А Н0 = 100 м скорость течения в западном погранслое ослабевает в верхней части моря до ~39 см/с, в нижней части до ~ 20 см/с; под воздействием рельефа дна она усиливается в южной части бассейна (в верхнем слое ~ 9 см/с, в нижнем ~4 см/с), в юго-западной части (в верхнем слое ~ 11 см/с, в нижнем ~5 см/с), в центральной (в верхнем слое ~ 2 см/с, в нижнем ~ 1 см/с) и в северной части (до ~2 см/с - только в нижнем слое). При АН0 = 150 м скорость течения в западном погранслое ослабевает (в верхней части моря до ~ 37 см/с, в нижнем слое до ~17 см/с); под воздействием рельефа дна она усиливается в южной части (в верхнем слое ~ 10 см/с, в нижнем ~ 4 см/с), в юго-западной части (в верхнем слое ~ 12 см/с, в нижнем ~ 5 см/с), в центральной части (в верхнем слое ~2 см/с, в нижнем ~1 см/с) и в северной части (до ~ 2 см/с - только в нижнем слое). На течения

78

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

в северо-западной и восточной частях бассейна рельеф дна практически не влияет. Таким образом, с учетом рельефа дна течение усиливается в южной, северной (для нижнего слоя), юго-западной и центральной частях бассейна и ослабевает в западном погранслое. Вертикальная структура течений описывается функцией В, определяемой как отношение модулей вектора скорости течения в верхнем и нижнем слоях. Для значительной части бассейна функция В имеет структуру, близкую к баротропной. Бароклинная структура проявляется в западном погранслое, восточной части моря, а также в областях минимума функции у2 (рис. 1 - 3).

У

Vi

B

Р и с. 1. Функция тока в верхнем и нижнем слоях и амплитудная структура при Д H0 = 0 ISSN 0233-7584. Мор. гидрофиз. журн., 2011, № 6 79

У

0 5 10 15 20 25 30 35 40 45 50

У

0 5 10 15 20 25 30 35 40 45 50

B

15

10-

5-

0 5 10 15 20 25 30 35 40 45 50

Р и с. 2. Функция тока в верхнем и нижнем слоях и амплитудная структура при Д H0 = 100 м

80

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

1

0

15

10-

5-

00 5 10 15 20 25 30 35 40 45 50

15- -

10-

5-

0- -

0 5 10 15 20 25 30 35 40 45 50

B

15- -

10-

5-

0- -0 5 10 15 20 25 30 35 40 45 50

Р и с. 3. Функция тока в верхнем и нижнем слоях и амплитудная структура при Д H0 = 150 м

Выводы. Рассмотренная нами простейшая двухслойная модель без учета горизонтальной вязкости и нелинейных слагаемых дает качественное представление о влиянии рельефа дна на характер ветровой циркуляции для замкнутого двухслойного бассейна. В целом интегральная циркуляция в нижнем и верхнем слоях под влиянием рельефа дна усиливается (рис. 1 - 3), причем скорости течения в нижнем слое близки к изолиниям функции f / H. Течение у западного погранслоя несколько ослабевает; в центральной, южной и юго-западной частях бассейна оно усиливается как в верхнем, так и в нижнем слое. Последний факт подтверждается в эксперименте со всплывающими буями.

ISSN 0233-7584. Мор. гидрофиз. журн., 2011, № 6 81

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

1. Korotaev G, Oguz T, Riser S. Intermediate and deep currents of Black Sea obtained from autonomous profiling Floats // Deep-Sea Res. - 2006. - 2, № 53. - P. 1901 - 1910.

2. Korotaev G, Saenko O, Koblinsky C. Satellite altimetry observation of the Black Sea level // J. Geophys. Res. - 2001. - 106, № C1. - P. 917 - 933.

3. Саркисян А.С., Иванов В.Ф. Совместный эффект бароклинности и рельефа дна, как важный фактор в динамике морских течений // Изв. АН СССР. ФАО. - 1971. - 7, № 2. -С. 173 - 188.

4. Зырянов В.Н. Теория установившихся течений. - Л.: Гидрометеоиздат, 1985. - 248 с.

5. Коротаев Г.К., Шапиро Н.Б. К вопросу о влиянии рельефа дна на океаническую циркуляцию (двухслойная модель) // Морские гидрофизические исследования. - Севастополь: МГИ АН УССР, 1971. - № 5. - C. 55 - 63.

6. Коротаев Г.К. Интенсификация глубинной циркуляции в мезомасштабном бассейне под влиянием рельефа дна // Морской гидрофизический журнал. - 2005. - № 2. - C. 3 - 10.

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

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

После доработки 24.06.10

АНОТАЦ1Я Дослщжуеться глибинна циркулящя у замкнутому басейш з нерiвним дном. Ура-хування стратифкацп розглянуте на прикладi двошарово! модел^ яка яюсно описуе термодинамку невеликих морiв i водоймищ. Вплив рельефу дна виявляеться методом порiвняння розрахованих швидкостей течп в нижньому шарi для басейну з плоским дном i басейну з параболiчним рельефом. Показана штенсифкащя течп в нижньому шарi для замкнутого басейну при облку рельефу дна.

Ключовi слова: в^рова циркулящя, вплив рельефу дна.

ABSTRACT Deep-water circulation in a closed basin with uneven bottom is studied. Consideration of stratification is examined on the example of the two-layer model qualitatively describing thermodynamics of small seas and water basins. The topography effect is revealed by the method of comparison of the calculated flow velocities in the lower layer of the basin with flat bottom with those in the basin with parabolic relief. At consideration of bottom topography, the flow intensification in the lower layer of the closed basin is shown.

Keywords: wind circulation, effect of bottom relief.

82

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