О решениях для внутренних захваченных волн с условиями прилипания на стенке

Автор: Смирнов Сергей Викторович

Журнал: Вычислительная механика сплошных сред @journal-icmm

Статья в выпуске: 3 т.1, 2008 года.

Бесплатный доступ

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

Еще

Короткий адрес: https://sciup.org/14320441

IDR: 14320441

Текст научной статьи О решениях для внутренних захваченных волн с условиями прилипания на стенке

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

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

работе [4] представлены результаты исследования совместного влияния вертикальной и горизонтальной вязкостей на волны Кельвина, длина которых превышает радиус деформации Россби.

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

2.    Основные уравнения

Рассмотрим несжимаемую вязкую жидкость, находящуюся в бассейне с плоским дном и одной прямой вертикальной стенкой. Будем пользоваться правой декартовой системой координат (х, y, z) с направленной вверх осью z и началом координат на невозмущенной поверхности. Пусть состояние покоя характеризуется плотностью р 0 (z) и давлением p 0 (z) . Запишем в приближениях Буссинеска и гидростатики систему основных уравнений теории океанских течений [5], линеаризованную относительно основного состояния покоя:

atu * - fv * = -р-1а хР * + ц(а xxu * + а yyu * )+а z (Ka zu *), a * + fu * = - р-1а yp * + ц (axxV * + a yyv *) + az (Ka zv *),

  • d    zp * = - g р*,

at р* + w *а z р0 = о, a w* + а u* + а v* = 0.(5)

zxy

Здесь и далее: t — время; и * , v * , w * — компоненты вектора скорости; Z — отклонение уровня жидкости от невозмущенного положения; p , р — отклонения давления и плотности от их невозмущенных значений; р * — характерное значение плотности; f — параметр Кориолиса, полагаемый постоянным; g — ускорение свободного падения; H — глубина бассейна, H > 0 при х 0; ц , K — коэффициенты горизонтальной и вертикальной турбулентной вязкости [5], полагаем ц > 0; N — частота Брента-Вяйсяля; d F — частная производная функции F по переменной 5 . Следует отметить [1], что приближение гидростатики может быть принято для крупномасштабных движений, в которых горизонтальный масштаб гораздо больше вертикального и становится несправедливым при рассмотрении тех движений, где необходим учет коротких внутренних волн. Запишем граничные условия на поверхности с учетом малости Z :

p*| = р*gz       w*| = дЛ-< z=0l

Зададим условия непротекания и прилипания на стенке и 1 .=0 = v 1 х=0 = 0(7)

и условия непротекания и свободного скольжения на дне w 1    =0,     a zu 1    = az v 1    =0.(8)

Iz=-H                   z lz=-H      z lz=-H- '

Оценка влияния трения о дно на внутренние захваченные волны выполнена в работе [3] на примере двухслойной жидкости. Показано, что в уравнениях движения для нижнего слоя трение о дно дает вклад в результирующий коэффициент трения. Отношение вклада трения о дно к вкладу трения между слоями равно rH 1 2/( r H 2), где H 1 — толщина верхнего слоя, r и г' — коэффициенты трения о дно и между слоями соответственно. Далее полагаем, что глубина H настолько превышает толщину термоклина, что влиянием трения о дно можно пренебречь. Исключая р из (3) при помощи (4), имеем

- 1 А * Л т 2 * р * 8 tz p = N w ,

где

N 2 = - g P- 1 d z Р о -

Известно [5], что структура уравнений (1), (2), (5), (9) позволяет искать решение задачи в виде:

и *= р»1 P (z) u (x, y, t),       v *= р*1 P (z) v (x, y, t);(11)

w*= - W(z)dt П(x, y, t);(12)

p * = P (z )П( x, y, t).(13)

Подставляя (11)-(13) в (5), запишем это уравнение как

81П   = 5 и + 5, v xy Р =        = -s2,                                                              (14) р,а zW где sn — постоянная разделения переменных. Рассмотрим простой случай [3], когда произведение KN2 является постоянной величиной. Подставив (11)-(13) в (1), (2), (9) с учетом (14) и краевых условий (6), (8), имеем где

8tu - fv = -sn2 8xn + ц (axxu + 8yyu) - du,5 tv + fu = -s „2 8 y n + ц (a xxv + 5 yyv)- dv,

d,n + d и + 8, v = 0,

txy

Й d P P н 8Z ' ~ = 0, zN 2 s 2

н = 0, zH

dz P + P

N 2 z

= 0, z=0

d = KN 2/ s 2 n , П = n / s 2

Совокупность собственных функций Pn(z) задачи (18), (19) обладает свойством полноты [5]. Введем безразмерные переменные согласно соотношениям:

T = ft,    X = x /Ln,     Y = yILn;(22)

U(X, Y, T) = u(x, y, t)/ en,   V(X, Y, T) = v(x, y, t)/ en,    h(X, Y, T) = n(x, y, t);(23)

Ln = f .(24)

Здесь Ln — радиус деформации Россби внутренней моды с номером n . Перепишем в безразмерных переменных уравнения «мелкой воды» (15)-(17) и краевые условия (7)

5 TU - V = -5 Xh + E (5 XXU + 5 YYU)- DU,(25)

5 TV + U = -5 Yh + E (5 XXV + 5 YYV)- DV,(26)

dTh id XU id YV = 0,(27)

U |x=0=0, V |x=0=0,(28)

где

E = -£-, D = d.(29)

fLnf

3. Захваченные волны

Обозначим через c, X и l = 2 nd безразмерную фазовую скорость, длину волны и волновое число соответственно; расстояние, на которое смещается гребень волны за время, пока амплитуда уменьшается в e раз, есть

Re ( c ) Im ( c ) l

Пусть { U m , V m , h m } , гЛе

IT = IT A ikmn^^ + il ( Y cT )

Um Ume          , h => A ikmX+il(Y-cT)

nm me e           ,

Im( c )<0,   Im( k m ) < 0,

V = V A e‘imX + lI ( Y cT ) m m                 ,

являются решениями системы уравнений (25)-(27). Решения типа волн, захваченных вращением у стенки, будем искать в виде:

( U , V , h ) = ( U 1 + U 2, V 1 + V , , h 1 + h 2 ) ,            Im( k 1 )>Im( k 2).

Для решений вида (34) краевые условия (28) выполняются при

U A + U 2a=0,      V 1 A + V 2a = 0.

Перепишем уравнения (25)-(27), подставив выражения для U , V и h из (31) и (32):

(а + e к 2 - icAlU A - V A + ikh A =0     ( m = 1,2),

m            m m m m

UA +(a + sк2 -ic\lVA + ilhnA=0      (m = 1,2), m              m            m n, /, lchA - kUA - lVA = 0                  (m = 1,2), m m m m                                     \, /, где а = D1l + El,    e = E1l.

Исключим h m из (36)-(38), а из полученных уравнений исключим U A и V 2 A при помощи (35); в результате имеем:

( ( ( e k 1 2 + а ) c - ic 2 ) l 2 + ik^ ) U A + ( ik 1 - c ) lV 1 A =0, (40) ( ( ( e к 2 2 + а ) c - ic 2 ) l 2 + ik 2 2 ) U 1 A + ( ik 2 - c ) lV 1 A=0, (41) ( c + ik 1 ) U A + ( ( e k А + а ) c + i ( 1 - c 2 ) ) /V/ 0, (42) ( c + ik 2 ) U A + ( ( e k 2 2 ) c + i ( 1 - c 2 ) ) V/ 0. (43)

Для существования нетривиальных решений системы уравнений (40)-(43) необходимо задать три условия. Запишем сначала условия существования нетривиальных решений для четырех систем уравнений, состоящих из [(40), (41)], [(42), (43)], [(40), (43)] и [(41), (42)] соответственно; это дает

( 1 - il 2 c e ) k1 k 2 + ( i + l 2 c e ) c ( k1 + k 2 ) + ( c + i а ) cl 2 =0,

i e ck 1 k 2 + e c 2 ( k 1 + k 2 ) - c 2 - ic а + 1 = 0,

[( i -L / 2^c) 2 -I- (in 4-^4- il2p\^ — A\lr24-(n — f^272 2-|_T2]_i- 11 i l c e ) e ck 2 I i cx + c + il e ) c 1 1 k 1     ( cx ic ) e c l I k 1 + k 2 ) +

+ k1 k 2 + ic ( k 2 - k 1 ) - l 2 c 4 - 2 il 2 а c 3 + ( l 2 а 2 + 1 + l 2 ) c 2 + il 2 а c = 0,

4- 12r-e\pr-k 2 -I-     -I- с -I- il 2f)^ Й 2 -I- (n — iA^^2!2 ( b-2 4- b- 2\ 4- i (i +1 ce) eck 1 + (iex + c + il e) c 1 pc2 (^а ic) ec l (k1 k2 ) +

+ k 1 k 2 + ic ( k1 - k 2 ) - 1 2 c 4 - 2 il 2 а c 3 + ( l2 а 2 + 1 + 1 2 ) c 2 + il 2 а c = 0.

Первую пару условий преобразуем, умножив (44) на c 2 e , а (45) — на ( i + l 2 c e ) c , а затем взяв сумму и разность полученных соотношений; во второй паре вычтем (47) из (46). Имеем

2ec (1 - icel2)k 1 k2 + 212c3e + (i + 2ie12а)c2 - (а + el2) c - i = 0,(48)

2(i + cl e)ec (k1 + k2) — ic +(el +а)c + i = 0,(49)

(cl 2e + i + c а-ic2)(k1 + k2)-2 c = 0.(50)

Исключая k 1 + k 2 из (49), (50), получим уравнение для c:

(1 - 4г2l2)c4 + i2 (l2г - 2г + a) c3 - (l4г2 + 2 + 2l^a + a2) c2 - i2 (l2г + a) c +1 = 0; (51)

исключая k 1 или k 2 из (48), (49), получим уравнение для k:

2(i + cl2г)гc2k2 + (i - ic2 + (гl2 + a)c)k -

—2iгl c +(21 га +1) c + i (l г + a) c — c — 0,      k e{k 1, k2}, решения которого имеют вид k1,2

—i + ic — (l г + a) c ± ^JM। 2M2

где

M 1 — (ca + i +12cг - ic2) - 8( l2cг + i )гc3 (ica -1 + c2 + 212c2гa- 2il2c3 г + il2c8) ,(54)

M 2— 2( l2 c г + i )г c2,(55)

M 1 ^ 0,     M2 ^ 0.(56)

Уравнения (51) и (52) и есть искомые условия существования нетривиальных решений системы (40).

Можно показать, что решения уравнений (51) и (52) удовлетворяют уравнениям (2.7) и (3.4) работы [4]:

г k4 (1 - il2 c г) + k2 (a- ic + г l2 (1 - 2 ic (a- ic )))-

- ic + ( a- ic ) l 2 ( 1 - ic ( a- ic ) ) — 0,               k e { k 1 , k 2 } ,

1 + ic г ( k 1 + k 2 ) 2 ( 1 - il 2 c г ) — 0.

Подставив в (57) и (58) выражения для k 1 и k 2 из (53) и применив второе условие из (56), получим факторизованные уравнения, в которых один из сомножителей совпадает с левой частью уравнения (51).

Пусть H = 1 км, f = 10-4 с-1, K = 0,01 м2с-1. Выберем два значения коэффициента горизонтальной турбулентной вязкости из рассмотренных в работе [6]: Ц 1 — 103 м2с-1 и ц 2 — 102 м2с-1. Оценим порядок величин a и г и относительную роль вертикальной и горизонтальной составляющих турбулентной вязкости в случае внутренней волны с X * — 25 км при L — 50 км, N — 10-2 с-1. Из (22) следует, что l — 2 п L / X * , где l « 12. Подставляя приведенные выше числовые значения и р — 103 м2с-1 в (20), (29) и (39), имеем D « 4^10-4, E « 4-10-3, a « 5^10-2, г « 340-4. Аналогично, при ц — 102 м2с-1 имеем D « 4-10-4, E ~ 4-10-4, a ~ 5-10-3, г ~ 3-10-5. Поскольку a = D/l+El « El , при данных условиях влиянием вертикальной турбулентной вязкости можно пренебречь.

Решения уравнений (51) и (52), удовлетворяющие условиям (33), представлены на рисунке 1 для трех наборов параметров D и E . Графики функций Re(c), Л , Re(k) и Im(k),

Рис. 1. Графики функций, зависящих от λ : а — Re( c ); б Λ ; в — Re( k ); г — Im( k ) (сплошные линии — D = 4 10-4, E = 4 10-3; штриховые линии — D = 4 10-4, E = 4 10-4; пунктирные линии — D = 0,1; E = 0,01)

зависящих от λ , показаны на рисунках 1, а–г соответственно.

Отметим, что на рисунке 1 не показаны решения с Λ < 0,02. Решения для волн с длиной больше радиуса деформации Россби ( λ > 1) получены и проанализированы в работе [4]. Важными особенностями решений в области волн c λ < 1 являются существенное уменьшение масштаба убывания захваченных волн с уменьшением длины волны и наличие волн с отрицательной действительной компонентой фазовой скорости c . Здесь и далее решения с Re( c ) > 0, соответствующие волне Кельвина, помечены символом F , решения с Re( c ) < 0 — символом B .

На рисунке 2 представлены два решения уравнений (25)–(27) с граничными условиями (28), полученные при λ = 0,6; D = 4∙10-4; E = 4∙10-3. Решения нормированы на максимальные значения h max = 1; ( U B ) max 0,21; ( V B ) max 1,03; Λ B –3,4; ( U F ) max 0,19; ( V F ) max 0,95; Λ F 1,7. На рисунках 2, а и г показаны изолинии h B и h F с шагом 0,2; на рисунках 2, б и д U B и U F с шагом 0,05; на рис. 2, е V B и V F с шагом 0,2. Изолинии с положительными значениями показаны сплошными линиями, с отрицательными значениями — штриховыми линиями, с нулевым значением — пунктиром. Следует отметить наличие тонких пограничных слоев на рисунках для U и V и существенное различие масштабов убывания для решений F и B .

Рис. 2. Изолинии h , U и V при X = 0,6; D = 4 - 10-4; E = 4 ^ 10 -3

(сплошные линии — положительные значения; штриховые линии — отрицательные значения; пунктир — нулевые значения).

4.    Асимптотические решения

Полагая, что £ является малым параметром, будем искать приближенные решения в виде асимптотических рядов:

c = c(0) + £1/2c(1) + £c(2) + £3/2c(3) +... ,(59)

k1 = k^ + £1/2k1(1) +£k^ + £3/2k1(3) + ... ,(60)

k2 =£-1/2k20) + k21) +£1/2k22) +£k23) +... .(61)

Подставляя (59)–(61) в уравнения (48)–(50) и приравнивая к нулю коэффициенты при разных степенях £ , получим последовательность уравнений для определения c (0), k 1 (0), k 2 (0) и так далее. Положив £ = 0 в (59) или (60), получим уравнение

(0)2

c -1 + i a c = 0, решениями которого являются

(0) F

- i a , 2

- i a .

Необходимо отметить, что рассматривается та ветвь функции квадратного корня, которая отображает плоскость c разрезом вдоль отрицательной действительной полуоси на правую полуплоскость. Для k2(0) имеем квадратное уравнение k20)2 - ic(0) + a = 0

и для него выбираем такие решения уравнения (64), которые удовлетворяют условию (33):

k20)=k20) I ™    =-1 V2—a - i1 v2+a,

2f     2 c(0) = c?)     22

k2°) = k2°)I     ™ =1V2 — a — i1V2 + a.

2в     2 c(0)=c0 2             2

Остальные переменные можно выразить через с (0) и k2(0). Например, k - ic(0),

(0) i + a c (0)      k (1) = (a 4 - 4 a 2 + 2 ) ic (0) + 3 a-a 3

2 2 ic (0) -a ,    2          ( 4 -a 2 )( 1 - i a c (0) )

Асимптотические решения для волны Кельвина получены в работе [4] и проанализированы при X > 1. Для k 1 (1) запишем приближенное выражение, полагая a <1:

(1)        (0) (0) 2

Л-1     ,* i *v2 c It

B lB

2 1/3

E 1/3 ,

X в

2п

l B

« пE1/322/3.

Необходимо отметить важную роль второго слагаемого в приближенном решении k ( a )

1F для волны Кельвина, обеспечивающего существенное уменьшение масштаба убывания с уменьшением длины волны. Графики |Im k(X)| (Рис. 3) построены в логарифмических координатах при D = 4∙10-4 и E = 4∙10-3. Вертикальные штриховые линии пересекают горизонтальную ось в точках с координатами Х|Л = 0,1 » 0,12, Xв « 0,8 и X|s = 0,1 * 160. Исключены из рассмотрения короткие волны с длиной X < Х|Л = 0,1, которые быстро затухают. Решениям k1 и k2 уравнений (51) и (52) соответствует пунктир, приближенным решениям k(a) = k1(0) + s1Z2k(1) и k2a) =£-1/2k20) + k21) — сплошные линии. Отметим, что в области относительно коротких волн графики для k1(a) и k2(a) практически совпадают с графиками для k1 и k2.

Рис. 3. Функции |Im k ( λ )| для k 1 и k 2 (пунктир), k 1 (a) и k 2 (a) (сплошные линии)

5.    Заключение

В работе рассмотрены решения для внутренних волн, захваченных вращением Земли у вертикальной стенки в бассейне с плоским дном. Система уравнений для захваченных волн приведена к удобному для решения виду: уравнение четвертой степени (51) для фазовой скорости волны c и квадратное уравнение (52) для волновых чисел k 1 и k 2 . Обратные значения мнимых частей k 1 и k 2 характеризуют масштабы экспоненциального убывания для захваченной волны и для пограничного слоя соответственно. В области относительно коротких волн ( λ < 1) решения характеризуются следующими важными свойствами: 1) масштаб убывания волн Кельвина существенно уменьшается с уменьшением длины волны; 2) имеется второе решение в виде волны противоположного направления (Re( c ) < 0). Построены асимптотические решения для c, k 1 и k 2 , применимые для относительно коротких волн, и получена оценка (69) для λ B — верхней границы диапазона длин волн второго решения.

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

Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (проект № 09-01-98519) и ДВО РАН (проект № 06-III-A-03-072).

Список литературы О решениях для внутренних захваченных волн с условиями прилипания на стенке

  • Марчук Г.И., Дымников В.П., Залесный В.Б. Математические модели в геофизической гидродинамике и численные методы их реализации. -Л.: Гидрометеоиздат, 1987. -296с.
  • Ле Блон П., Майсек Л. Волны в океане. Пер. с англ. -М.: Мир, 1981. -В 2-х томах. -845с.
  • Martinsen E.A., Weber J.E. Frictional Influence on Internal Kelvin Waves//Tellus. 1981. -V. 33, N. 4. -P. 402-410.
  • Davey M.K., Hsieh W.W., Wajsowics R.S. The Free Kelvin Wave with Lateral and Vertical Viscosity//J. of Phys. Oceanogr, 1983. -V. 13, N. 12. -P. 2182-2191.
  • Каменкович В.М. Основы динамики океана. -Л.: Гидрометеоиздат, 1973. -240с.
  • Zalesny V.B. Numerical Simulation and Analysis of the Sensitivity of Large-Scale Ocean Dynamics//Russ. J. Numer. Anal. Math. Modelling. -1996. -V. 11, N. 5. -P. 421-443.
Статья научная