Вычислительные технологии
Том 16, № 1, 2011
Моделирование умеренно длинных нелинейных пространственных волн на границе двух
>к
жидких потоков в горизонтальном канале*
Д. Г. Архипов1,2, Н. С. Сафарова1, Г. А. Хабахпашев1,2 1 Институт, теплофизики СО РАН, Новосибирск, Россия 2 Новосибирский государственный университет, Россия e-mail: theory@itp.nsc.ru
Проведено моделирование слабонелинейных планарных волн, которые бегут одновременно в разные стороны по границе раздела двух неглубоких слоев жидкостей, текущих с постоянными скоростями. Рассчитаны взаимодействия как протяженных, так и уединенных возмущений малой, но конечной амплитуды. Показано, что может наблюдаться интенсификация волн при их движении против потока и ослабление при распространении в спутном направлении.
Ключевые слова: взаимодействие возмущений, граница раздела, двухслойная жидкость, длинные волны, моделирование, неглубокие слои, нелинейные волны, солитон, стационарный поток, трехмерные возмущения, эволюция.
Изучению длинных волн на свободных границах слоев жидкостей посвящено много работ (см., например, литературу, цитируемую в монографиях [1-8]). В последнее время наибольший интерес исследователей вызывают плапарпые нелинейно-дисперсионные модели. В работе [9] предложен новый подход для моделирования возмущений поверхности первоначально покоящейся неглубокой жидкости, использующий основные достоинства известных моделей и в то же время свободный от некоторых их недостатков. Данный подход заключается в выводе одного нелинейно-дисперсионного уравнения для существенно трехмерного возмущения поверхности жидкости и нескольких простых линейных вспомогательных уравнений, определяющих векторы скорости частиц среды, которые входят лишь в члены второго порядка малости главного уравнения. В статье [10] этот метод обобщен на случай внутренних волн в водоемах со скачком плотности.
Цель настоящей работы — применение нового подхода к описанию эволюции умеренно длинных слабонелинейных трехмерных возмущений границы раздела двух слоев жидкостей, ограниченных горизонтальными крышкой и дном, при наличии стационарных течений и исследование влияния установившихся потоков на трансформацию волн с помощью численных решений ряда характерных гидродинамических задач.
* Работа выполнена при финансовой поддержке гранта Правительства России для государственной поддержки научных исследований, проводимых под руководством ведущих ученых в российских вузах, № 11.G34.31.0035.
1. Постановка задачи и вывод системы модельных уравнений
Будем предполагать, что, во-первых, горизонтальные крышка и дно канала неподвижны и недеформируемы, во-вторых, несжимаемые и невязкие жидкости не смешиваются, в-третьих, внутри каждого слоя неизменны плотность жидкости р и вектор скорости стационарного потока и, в-четвертых, "длина волны" Л существенно больше, а амплитуда возмущения ца значительно меньше равновесных глубин слоев кг (кг/Л ^ Па/кг ~ е, оде е — малый параметр, индекс I = 1, 2 — соответственно для параметров верхней и нижней жидкости), в-пятых, капиллярные эффекты невелики (модифицированное число Бонда Во = (р2 — р\)дЫк2/а > 1, здесь д — ускорение свободного падения, а — поверхностное натяжение), и наконец, в-шестых, течения в слоях потенциальны. Тогда исходную систему уравнений гидродинамики для каждой жидкости можно записать в виде
Биг _ /и2 рД
+ ^ + * = (2)
г рг дг
о»
где V = (д/дх,д/ду), Б/Бгг = д/дг + (иг • V) х и у — горизонтальные координаты (рис. 1), и — горизонтальная компонента вектора возмущения скорости жидкости, т — его вертикальная составляющая, г — вертикальная координата, г — время, р — давление в слое.
Поставим стандартные кинематические и динамическое краевые условия (условия непротекания жидкости через неподвижные крышку и дно, кинематические условия на границе раздела и условие капиллярного скачка давления на этой поверхности):
= 0 при г = — (—1)гкг(х, у),
= Бп/Бгг + ин • Vn, Рн = Р2г + а V2ц при г = ц(г, х, у),
где ц — отклонение границы раздела, индексом г помечены значения величин, относящиеся к этой поверхности.
Теперь преобразуем систему (1)—(3). Вначале найдем дивергенцию горизонтальных уравнений движения (1). Тогда в первых слагаемых левых частей этих уравнений дивергенции скоростей V • иг согласно уравнениям неразрывности можно заменить
Рис. 1. Волны на границе раздела и профили скоростей стационарного и возмущенного течений в неглубокой двухслойной жидкости между неподвижными крышкой и дном
на —dwi/dz. Затем проинтегрируем полученные уравнения и уравнения (3) по z от -h2 до п ПРИ l = 2 и от п до hi при l = 1, В результате получим следующие выражения:
Dit dz 2
- -V (uf) - -V2Pl
Jll dz 4- V. I
Di t
Dit dz 2 pi
Dn
dz = 0, (4)
-i + v- (Ы [г/+ (-l)%])=0. (5)
Здесь и далее угловыми скобками обозначены значения величин, усредненные по глубинам слоев. При сделанных предположениях (длинные волны малой амплитуды) профили вертикальных скоростей линейны по г:
ад = 1№(1 + (-1)'±) = (1 + (-1)-±)@ + и11^Ч). (6)
Если законы сохранения массы в слоях (5) и формулы для вертикальных скоростей жидкостей (6) линеаризовать, то их можно записать в виде
^ = (7)
Эти формулы понадобятся для замен в членах второго порядка малости основного уравнения.
Первые члены интегралов в (4) можно легко трансформировать с помощью правила Лейбница о дифференцировании по параметру, кинематических краевых условий на границе слоев, упрощенных формул для вертикальных скоростей (8) и уравнений (7):
П D dw^ Dwli Dn dwl -dz =
J Dit dz Dit Dit dz
(-i)l+1 hi
D2n D r. _ .
z=v Dit2 Dit
(mï=0 -- ж • ■" -- •• <-»+• <u<»2>-
При сделанных допущениях во вторых членах интегралов (4) можно заменить верхние пределы на 0 и вынести производные за знаки интегралов:
/ iv2(u f)dz= J iv2(u2)dz = (-l)^v2(u2).
(-i)l+1hi (-i)l+1hi
Наконец, для членов второго порядка малости воспользуемся также приближенными равенствами uii = (ui) и (uf) = (ui)2,
Для нахождения давлений в слоях проинтегрируем уравнения (2) по z от z до п-Используя динамическое краевое условие на границе раздела, имеем
п
Pl Pli . , ч . [ DwKl - = — + g{r] - z) + / -jr-dz. pi pi J Dit
Учитывая, что вертикальные скорости входят лишь в члены второго порядка малости, подставим в эти формулы линеаризованные профили вертикальной скорости (8):
Рг _ Рч , .л „ Л , ( л\1 г \
Рг Рг V 2Ни Б^2
С помощью данных выражений находим третьи члены интегралов в уравнениях (4):
1
п
Pi J
(-l)l+1hi
V2Pldz = [(-!)% + V] V2 + дЛ - f V20.
Еще раз заметим, что при вычислении интегралов в уравнениях (4) поправочные члены третьего и более высоких порядков малости, как всегда, опущены,
В итоге проведенные преобразования линейных, нелинейных и дисперсионных членов позволяют переписать уравнения (4) в следующем виде:
77 + 7ГГ • ^ - (-DWW) - (-1 )^v2№-
Dit2 Dit pi
-(-1 Ygh^n - VV2 + g^ - f V20 = 0. (9)
Чтобы исключить величины лапласианов давлений на границе раздела слоев из линейных членов уравнений (9), умножим эти уравнения на h2/p2 при l = 1 и нa hi/pi при l = 2, Сложив модифицированные таким образом уравнения, получаем
h2 D2n hi D2'n (h2 Duii hi Du2i\ , hih2—2 ( , \2 , \2\
— 7T7Ï + -7T7Ï +--FTT +--FTT ' +-V (Pi<ui> "ft u2 )-
P2 Dit2 Pi D2t2 P2 Dit Pi D2t PiP2
+ ghihv (— ~ —) V2r] - —V2(h2pu + hip2i)~
PiP2 P2 Pi PiP2
(h2 hA 2 hih2 2 /h2 D2r] hx D2r] \ \p2 pi ) 3 \P2 Dit2 Pi D2t2J
Величины Dun/Dit входят лишь в член второго порядка малости и поэтому с использованием линеаризованных уравнений (1) могут быть заменены на —Vpii/pi. Эти выражения можно сгруппировать с последним членом во второй строке уравнения (10), содержащим рц. В итоге получаем член —V[n ■ (h2Vpii + hiVp2i)]/(pip2), имеющий второй порядок малости.
Следовательно, Vp^ могут быть найдены из первого приближения, но тогда pii = pi, а в уравнениях (9) нужно опустить нелинейные и дисперсионные члены:
^-(-1У^2Рг-(-1УдЫ^ = 0. (11)
Dlt2 Pl
l =2 l =1
ближенную формулу для лапласиана давления на границе раздела слоев
V2pi = V2pf — pip2gHV2n/X,
где р/ — часть давления рг, определяемая скоростями стационарных течении
= ^ (ад - + [(и2 • V)2 - (и! • V)2],) , (12)
Н = Ь1 + Ь2, х = Р1Ь2 + р2Используя эти выражения и умножая все члены (10) на Р1 Р2/х-, основное уравнение можно переписать в следующем виде:
- с2^ + 2(и12 + и21) • ^ + - - (Р2(и2)2 - Р1(и1)2) -
• - - ■ - Ъ^п + = 0. (13)
Здесь с0 = д(р2 — р^ЬЬ/х, и12 = р^И/х И21 = Р2^и2/х V2 = [р1 Ь^И • V)2 + Р2Ь1(И2 • V)2]/х « = с0(р2Ь2 - Р1Ь2)/(Ь1к2х), в = Ь1Ь2(р1 Ь1 + р2Ь2)/(3х), ^з = Ь1Ь2(Р1Ь1И1 + Р2Ь2И2)/(3х^ ^ = ьЬ2[р1 Ь1 (И1 • V)2 + Р2Ь2(И2 • V)2]/(3х).
В уравнении (13) возмущения скоростей жидкостей остались только в квадратичных членах, поэтому для их определения достаточно простых уравнений первого приближения, т.е. линеаризованных уравнений горизонтального движения (1), Если для этих скоростей ввести потенциалы <рг, то для них справедливы уравнения Пуассона
V2?! = М = V**. (14)
Ь Пг г
Таким образом, полученная система модельных уравнений состоит из одного главного уравнения (13) и нескольких простых линейных вспомогательных уравнений (12) и (14), При отсутствии установившихся течений в обоих слоях система уравнений (13) и (14) согласуется с системой, полученной в [10],
2. Дисперсионный анализ и солитонные решения модели
В общем случае дисперсионное соотношение (зависимость циклической частоты ш от волнового числа к) для основного модельного уравнения (13) даже после его линеаризации весьма громоздко. Поэтому рассмотрим простую ситуацию, когда стационарный поток, направленный вдоль координаты х, имеется лишь в одном из слоев (И2 = (и2= (и2,0) и И1 = 0), а плоские возмущения также распространяются только вдоль этой координаты (все производные по координате у равны нулю). Подставляя в (13) решение в виде п = Па{ехр[г(кх — ш£)] + ехр[—1(кх — ш£)]}, находим дисперсионное соотношение в следующей форме:
,^21 + М2 , ш = к-——— ±
1 +рк2
, Ли21 + М2)2 + с2 - и2и21 - (1 + 13к2)к2(112и21Щ/3 - акМх) , ч
±к-гт^-• (15)
Заметим, что произведение и2и21, а также параметры в и х больше нуля. Тогда из зависимости (15) следует, что при наличии стационарного относительного движения
слоев жидкости решение может содержать растущие моды в коротковолновой области (аналог неустойчивости Кельвина—Гельмгольца). При отсутствии потоков жидкости уравнение (13) является корректным вариантом уравнения Буссинеска. При учете установившихся сдвиговых течений слоев в основном уравнении модели возникает ряд дополнительных членов. Один из них имеет четвертую производную по пространственным переменным, и при пренебрежении капиллярностью (а = 0) данное уравнение становится аналогом некорректного варианта уравнения Буссинеска. Как известно (см., например, [11]), некорректность таких уравнений связана с положительным знаком при старшей производной по горизонтальным координатам. При учете поверхностного натяжения вклад этого члена уменьшается, а при малых скоростях потоков жидкостей может быть полностью нейтрализован.
На рис. 2 показаны нейтральные кривые, вычисленные при значениях параметров задачи, соответствующих условиям проведения лабораторных опытов [12, 13] по распространению одиночных возмущений границы раздела керосина и воды. Слева от кривых — области устойчивости, справа — неустойчивости. На рис. 3 представлены зависимости обезразмеренной мнимой части фазовой скорости с = и/к от обезразмеренного волнового числа, рассчитанные при двух значениях скорости стационарного течения нижнего слоя (одно из них совпадает с тем, при котором выполнялся эксперимент [13]). Как и следовало ожидать, при учете поверхностного натяжения области устойчивости увеличиваются, что особенно заметно при малых скоростях потока жидкости.
В начале раздела 1 мы предположили, что исследуемые волны являются длинными. Поэтому может показаться, что рассмотрение нейтральных и дисперсионных кривых при кН > 5 не имеет смысла. Однако хорошо известно (см., например, [1]), что при численных расчетах различные эффекты перехода к конечным разностям и округления вводят малые высокочастотные осцилляции даже в том случае, когда решаемая задача удовлетворяет условию длинноволновости.
Модельная система уравнений имеет плоские установившиеся бегущие решения типа кноидальных волн. В частности, для уединенных возмущений, распространяющихся в направлении роста координаты х, получаем формулу
2
/
20 40 60 80 кН
Рис. 2. Кривые нейтральной устойчивости при р1 /р2 = 0.8, а = 0 (кривые 1 л 2) и а = 34 мН/м (крив ые 1' и 2!\ а также Н = 31.5 мм, Н1 /Н2 = 2.7 (кривые 1 и 1') и Н = 60 мм, Н1 /Н2 = 1.26 (кривые 2 и 2')
Рис. 3. Зависимости мнимой части фазовой скорости от волнового числа при Н = 60 мм, Ь1 /Н2 = 1.26, а = 0 (кривые 1 и 2) и а = 34 мН/м (крив ые 1' и 2'), а так же и2/с0 = 0.5 (кривые 1 л 1') и и2/с^ = 1 (кривые 2 и 2')
(16)
Здесь скорость волны С находится из решения квадратного уравнения
а
а
Чс20-и1и12-и2и21)-2г]а - +
--1----
2 /12 /ц 2Х2
и2 и21 СД и 12 р1 Р2Н
+
(С2 - С2) =о,
ее ширина Ь8 = 2^ В/А, где А = II2 - с\ - 2173(и12 + [/21) + С/1 С/12 + С/2С/21, В = /3112 -
акЬ/х - [2и8(Ь2и12 + ^21) + Ь2С/1 и 12 + Ь2СЦ^/З.
Отметим, что без потери общности скорость течения в одном из слоев или С/2) можно положить равной нулю, перейдя в соответствующую инерциальную движущуюся систему отсчета, В дальнейшем будем считать, что, как и в лабораторных экспериментах [12, 13] по распространению нелинейных уединенных плоских волн на сдвиговом потоке, и2 = 0,
3. Методика численного решения системы модельных уравнений
Пусть Л«, Лхх, Луу — операторы разностной аппроксимации вторых производных по соответствующим переменным со вторым порядком па трехточетном шаблоне, а Л:, Лх, Лу — аппроксимации первых производных центральными разностями. Тогда для уравнения (13) будем иметь трехслойную неявную конечно-разностную схему
Лйп - (Лхх + Луу)п + 2Л4[(С12 + СС21 )Лх + (У12 + Уи)Лу]п+ + [(грС2 + тОЛхх + 2(грС1^1 + тн С У2)Лх Лу + (трУ2 + тн ^2)Луу]п/х--а(Лхх + Луу)п2/2 + Ь1(Лхх + Луу)[тр(п21 + у2) - и2 - у2,}-
-[Лx(nЛxPf) + Лу(nЛyPf)]/х - вЛц(Лхх + Луу)п--2Ь1 [(трт^ С + С2)Лх + (тртн У1 + У2)Лу ](Лхх + Луу)Л:п/(3х)--Ь1Ь-2 [^ С2 + С/2 )Лхх + + С2У2)ЛхЛу + (^тУ + У>2)Луу ]х
х(Лхх + Луу)п/(3х) + Ь,1Ь,2(Лхх + Луу)(Лхх + Луу)п/Во = о,
где тр = р1 /р2, ть = Ь1 /Ь2 и все величины обезразмерены с помощью характерных параметров задачи: линейного размера Н и фазовой скорости очень длинных линейных волн в покоящейся жидкости с0, Указанная схема имеет второй порядок аппроксимации по всем трем независимым переменным. Решение итоговой системы алгебраических уравнений, полученной в результате разностной аппроксимации уравнения (13) по вышеприведенной схеме, проводилось с помощью итерационного метода Зейделя, который обладает хорошей сходимостью,
В общем случае исследование задачи на устойчивость чрезвычайно громоздко. Поэтому рассмотрим две простые, но характерные частные ситуации. Если стационарные течения жидкостей отсутствуют в обоих слоях, то условие устойчивости имеем в виде
(Д£)2 <
4в + (Д/)2
(Д/)2
(Дх)2(Ду )2
1 + паа + 4ЬЬ/[(Д/)2Во]'
(Дх)2 + (Ду )2'
Здесь ДЬ, Дх, Ду — шаги разностной сетки по времени и горизонтальным координатам. При р1 = 0 (слой однородной жидкости) указанное условие согласуется с полученным в [14] условием для подобного уравнения (отличаются только формы учета нелинейности и капиллярности).
Другим важным частным случаем является ситуация, когда стационарный поток, направленный вдоль координаты х, есть лишь в одном из слоев (например, нижнем) и плоские возмущения распространяются только вдоль этой координаты. Тогда находим
ш? <_4/3 + (Аж)2_
10. 1\ 9 тт , У*грЩ , АЫЬъ (Щ X Л1 '
1 + - х Г 2 1аЫ 2 (Лж)2 Во/_
Приведенное неравенство позволяет выявить влияние установившегося течения на условие устойчивости.
Входящие в уравнение (13) величины градиента давления на границе раздела слоев и векторов скоростей жидкостей определялись из простых линейных вспомогательных уравнений (12) и (14), Эти уравнения Пуассона решались с помощью метода быстрого преобразования Фурье по обеим горизонтальным координатам.
Ниже рассмотрены задачи, для которых в качестве исходного возмущения берутся либо стационарно бегущие плоские уединенные волны, либо первоначально покоящиеся возмущения. Следовательно, задание частных производных по времени на первых двух шагах по времени осуществляется очевидным образом.
На краях расчетной области по координатам х и у ставятся условия отражения от вертикальных стенок, заключающиеся в задании производных от отклонения границы раздела п по нормалям к стенкам равными нулю. Эти условия были реализованы следующим образом: значения искомой функции п в краевых и соседних с ними узлах брались равными ее значениям в ближайших внутренних узлах,
4. Результаты расчетов по распространению возмущений
Описанные в предыдущем разделе методика расчетов и разностная схема были протестированы с помощью аналитического решения (16), Как и следовало ожидать, возмущения двигались с постоянной скоростью и с неизменной формой. Так, при распространении солитонов па расстояния порядка 500^ их амплитуда изменялась менее чем на 1,5%, а закон сохранения массы выполнялся с точностью 10_6, Численные решения были устойчивы и сходились со вторым порядком по всем независимым переменным (времени и горизонтальным координатам). Подчеркнем, что солитонное решение (16) чрезвычайно чувствительно к ошибкам вычислений, так как его профиль формируется только при балансе вкладов нелинейных и дисперсионных членов, имеющих второй порядок малости.
На рис, 4 и 5 показаны профили границы раздела (п* = п/Ь2) при "лобовом" столкновении двух плоских уединенных волн (т, е, волны бегут как в направлении роста коорди-
х
но и различной полярности, В первом случае в момент наибольшего взаимодействия солитонов (при Ь = 24.2) максимальная величина отклонения границы раздела немного больше суммы амплитуд исходных волн. Во втором случае в аналогичный момент
Рис. 4. Трансформация двух плоских солитонов одной полярности П > 0) при их "лобовом" столкновении (значения параметров расчета: р1 /р2 = 0.8, Н1 /Н2 = 2.7, и2 = 0.5, Во = 11)
Рис. 5. Эволюция двух плоских солитонов разной полярности (пв < 0 и пв > 0) при их "лобовом" столкновении (значения параметров расчета: р1 /р2 = 0.8, Н1 /Н2 = 1.26, и2 = 0.5, Во= 51)
времени (при Ь = 49.4) наблюдается только отрицательное возмущение, поскольку соли-тон "впадина" шире солитона "бугор", а амплитуды начальных волн одинаковы. В дальнейшем уединенные возмущения возвращаются к своим исходным профилям.
В соответствии с формулами (15) и (16) скорости и ширины солитонов больше для спутных возмущений, чем при движении в противоположную сторону. Вычисления проведены для тех же значений параметров задачи, при которых были выполнены эксперименты [12, 13] по распространению волн против потоков.
х х
случаях. Следовательно, на характерные горизонтальные размеры волн приходилось не менее 50 узлов сетки.
Рис. 6. Трансформация первоначально аксиально-симметричного покоящегося колоколообраз-ного возмущения при р\/р2 = 0.8, /Н2 = 2.7, и2 = -0.5, Во = 11: а — £ = 0; 5 — £ = 80; в — сечения границы раздела при у = 0 в четыре момента времени
Рис. 7. Эволюция первоначально покоящегося крестообразного возмущения при р\/р2 = 0.8, Ь,1 /Н2 = 2.7, и2 = -0.3, ^2 = -0.4, Во = 11: о — Ь = 0; б - Ь = 20; в - Ь = 40
Рис. 8. Столкновение четырех квазиплоских уединенных возмущений при р1 /р2 = 0.8, Н1 /Н2 = 1.26, и2 = -0.5, Уч. = -0.1, Во = 51: а - £ = 0; б - £ = 49.4; в - £ = 100
Отметим, что для плоских внутренних волн, бегущих только в одну сторону в произвольно стратифицированном сдвиговом потоке уравнение Кортевега — де Вриза было впервые получено в работе [15]. Расширенное уравнение Кортевега — де Вриза для аналогичных возмущений сначала было выведено в [16] при слабой стратификации, а затем без использования этого приближения Буссинеска в работе [17]. Однако в последней статье проиллюстрированы лишь случаи, когда и > 0, а и2 = 0.
На рис. 6 представлена эволюция уединенного в пространстве возмущения, которое в начальный момент времени покоилось, а его форма была осесимметричпой и коло-колообразной: п(х,у) = По ехр[(ж2 + у2)/100]. Как и ожидалось, с некоторого момента времени за фронтом основной волны граница раздела опускается ниже равновесного уровня, а затем медленно стремится к этому уровню. Чтобы влияние стационарного течения было более очевидным, масштаб вертикальной оси на рис. 6, б" по сравнению с масштабом на рис. 6, о изменен в пять раз. Профили возмущения в наиболее интересном сечении показаны на рис. 6, е. Что касается значений безразмерных шагов по времени и горизонтальным координатам, то они в данном случае были равны 0.25.
На рис. 7 приведена трансформация крестообразного возмущения, которое первоначально также покоилось. Оно представляло собой суперпозицию четырех плоских уединенных решений: одно решение для волны, распространяющейся в направлении х
нии, и еще два аналогичных решения — по координате у, т. е. п = Пах^сЪ2(ж/Ьзх+) + Пах^сЪ2(ж/Ьзх-) + Пау^с:Ъ2(у/Ьзу+) + пауЗесЪ.2(у/Ьзу-), где индексы + и — отвечают случаям из > 0 и из < 0 соответственно (рис. 7, о). Хорошо видно, что образуются четыре квазиплоские волны, разбегающиеся в четыре стороны. На пересечениях этих волн имеются четыре пика возмущения границы раздела. Спутные возмущения превращаются в солптоны со следующими за ними волновыми пакетами. Возмущения, распространяющиеся против стационарных течений, преобразуются вначале в треугольные с крутым передним и растянутым задним фронтами (рис. 7, б), затем из них выделяются цепочки солитопов спадающей амплитуды (рис. 7, е). В данных расчетах ДЬ = Дх = Ду = 0.125. Амплитуды возмущений снижаются в спутных направлениях и возрастают в противоположных.
На рис. 8 показано взаимодействие четырех плоских уединенных волн, попарно движущихся навстречу друг другу. При этом возмущения, распространяющиеся вдоль ко-
ух рн ¡личную. В данном случае значения безразмерных шагов по времени и горизонтальным координатам были равны 0.2.
При решении всех пространственных задач расчетная область имела размеры 1024 х 1024
ких десятков узлов разностной сетки.
Заключение
В работе рассмотрена динамика умеренно длинных пространственных волн малой, но конечной амплитуды, одновременно распространяющихся под любыми углами друг к другу по границе раздела сдвигового потока двух жидкостей в горизонтальном канале. Анализ таких задач с помощью стандартных систем уравнений "мелкой воды" (например, уравнений движения (1) и законов сохранения массы (5)) чрезвычайно
затруднителен, так как векторы горизонтальных скоростей сред содержатся не только в линейных, но и в нелинейных членах всех уравнений.
Предложена комбинированная модель, которая состоит из одного основного нелинейного эволюционного уравнения для возмущения границы раздела и нескольких простых вспомогательных линейных уравнений. Последние необходимы для нахождения вклада в градиент давления на границе раздела слоев, связанного со стационарными течениями, и векторов горизонтальных скоростей жидкостей, входящих лишь в члены второго порядка малости главного уравнения.
Численные эксперименты подтвердили пригодность представленного подхода для решения ряда характерных плоских и пространственных задач динамики нелинейных волн. Использование полученной модели позволило исследовать как трансформации одиночных возмущений, так и различные их взаимодействия (столкновения и рн ¡бегания). Определено влияние скорости установившихся потоков на эволюцию протяженных и уединенных в пространстве возмущений. Продемонстрировано, что по границе раздела сдвигового течения способны одновременно распространяться плоские или квазиплоские уединенные волны разной полярности. Показано, что может наблюдаться ослабление возмущений, движущихся в спутном направлении, и, наоборот, интенсификация волн, бегущих против потока.
Авторы выражают благодарность I I.I I. Гешеву, В.К). Ляпидевскому, 11.11. Макаренко, E.H. Пелиновекому, К).Л. Степанянцу, Г.С, Хакимзянову, О.К). Цвелодубу и Л,Б, Чубарову за ценные советы и полезные обсуждения ряда вопросов.
Список литературы
[1] Уизем Дж. Линейные и нелинейные волны. М.: Мир, 1976.
[2] Овсянников Л.В., Макаренко Н.И., Налимов В.И. и др. Нелинейные проблемы теории поверхностных и внутренних волн. Новосибирск: Наука, 1985.
[3] Шокин Ю.И., Чуваров Л.В., Марчук А п. г.. Симонов К.В. Вычислительный эксперимент в проблеме цунами. Новосибирск: Наука, 1989.
[4] Пелиновский E.H. Гидродинамика волн цунами. Н. Новгород: Изд-во НПФ РАН, 1996.
[5] Ляпидевский В.Ю., Тешуков В.М. Математические модели распространения длинных волн в неоднородной жидкости. Новосибирск: Изд-во СО РАН, 2000.
[6] Хакимзянов Г.С., Шокин Ю.И., Барахнин В.В., Шокина Н.Ю. Численное моделирование течений с поверхностными волнами. Новосибирск: Изд-во СО РАН, 2001.
[7] Островский Л.А., Потапов А.И. Введение в теорию модулированных волн. М.: Физ-матлит, 2003.
[8] Шамин Р.В. Вычислительные эксперименты в моделировании поверхностных волн в океане. М.: Наука, 2008.
[9] Архипов Д.г., Хавахпашев Г.А. Новый подход к описанию пространственных нелинейных волн в диспергирующих средах // Докл. АН. 2006. Т. 409, № 4. С. 476-480.
[10] Архипов Д.Г., Сафарова Н.С., Хавахпашев Г.А. Моделирование нелинейных пространственных внутренних волн в морях и океанах со скачком плотности и пологим дном // Фундамент, и прикл. гидрофизика. 2009. Т. 2, № 2. С. 67-76.
[11] Daripa Р., Ни a W. A numerical study of an ill-posed Boussinesq equation arising in water waves and nonlinear lattices: Filtering and regularization techniques // Appl. Math. Comput. 1999. Vol. 101. P. 159-207.
[12] Букреев В.И., Гаврилов Н.В. Экспериментальное исследование уединенных внутренних волн в двухслойной жидкости // Там же. 1983. Т. 24, № 5. С. 51-56.
[13] Гаврилов Н.В. Уединенные внутренние волны большой амплитуды в двухслойной жидкости // ПМТФ. 1986. Т. 27, № 5. С. 49-54.
[14] Литвиненко А.А., Хавахпашев Г.А. Численное моделирование нелинейных достаточно длинных двумерных волн на воде в бассейнах с пологим дном // Вычисл. технологии. 1999. Т. 4, № 3. С. 95-105.
[15] Benney D.J. Long nonlinear waves in fluid flows //J. Math, and Phvs. 1966. Vol. 45. P. 52-63.
[16] Пелиновский E.H., полухина O.E., Лэмв К. Нелинейные внутренние волны в океане, стратифицированном по плотности и течению // Океанология. 2000. Т. 40, № 6. С. 805-815.
[17] Пелиновский Е.Н., Полухина О.Е. Уравнения Кортевега — де Вриза высшего порядка для внутренних волн в стратифицированных сдвиговых потоках // Изв. Акад. инж. наук. 2000. Т. 1. С. 117-133.
Поступила в редакцию 20 апреля 2010 г., с доработки — 22 ноября 2010 г.