УДК 532.59.032
А.А. Абрашкин, Ю.П. Бодунова
ПРОСТРАНСТВЕННЫЕ СТОЯЧИЕ ВОЛНЫ НА ПОВЕРХНОСТИ
ВЯЗКОЙ ЖИДКОСТИ
Нижегородский государственный педагогический университет
Построено решение для стоячих пространственных линейных гравитационно -капиллярных волн на поверхности бесконечно глубокой вязкой жидкости. Используются лагранжевые координаты. Дан качественный анализ дисперсионного уравнения. Движение жидких частиц изучены в длинноволновом приближении. Найдены выражения для компонент завихренности.
Ключевые слова: волны на воде, вязкая жидкость, завихренность, переменные Лагранжа.
Изучение влияния вязкости на распространение поверхностных волн традиционно ограничивается рассмотрением плоских (двумерных) волн [1, 2]. Эффекты, связанные с учетом поперечной модуляции профиля волны, до настоящего времени не исследовались.
Отличительная особенность исследования эволюции поверхностных волн в вязкой жидкости состоит в том, что толщина пограничного слоя вблизи свободной поверхности может быть мала по сравнению с амплитудой волны. Вследствие этого ставить поверхностные граничные условия на уровне у = 0 (у - вертикальная декартовая координата) не представляется удовлетворительным. Лучшим вариантом является выбор системы отсчета, движущейся с волной, и использование ортогональной криволинейной координатной системы, в которой свободная поверхность является координатной поверхностью. При эйлеровом описании течения ее вид неизвестен, поэтому при изучении двумерных линейных волн граничные условия формулируются на плоской поверхности у = 0, а при рассмотрении квадратичных эффектов на синусоидальной поверхности, соответствующем профилю линейной волны [2, 4, 5]. Анализ эффектов в третьем приближении по малой крутизне волны предполагает знание профиля во втором порядке теории возмущений и т.д. Изучение пространственных волн при таком подходе представляется нереальным ввиду значительного объема требуемых вычислений.
Все эти трудности, однако, отсутствуют при лагранжевом описании волнового движения. Форма свободной поверхности при таком подходе задается условием Ь = 0, где Ь -вертикальная лагранжевая координата. Оно не зависит ни от порядка теории возмущений, ни от размерности задачи. В рамках лагранжевого подхода в работах [6, 7] выполнен ряд расчетов для двумерных поверхностных волн. В частности, выяснено влияние на дрейфовое течение в жидкости движения воздуха, вращения Земли [6] и упругих пленок на свободной поверхности [7]. Ниже построено и проанализировано решение для линейных пространственных стоячих волн в бесконечно глубокой жидкости.
Постановка задачи
Рассмотрим пространственную стоячую гравитационно-капиллярную волну в бесконечно глубокой, вязкой жидкости. Движение будем описывать в лагранжевых координатах а, Ь, с . Система уравнений гидродинамики в этих переменных имеет следующий вид [8]:
© Абрашкин А.А., Бодунова Ю.П., 2011.
х, = -1 [у, г, р]++у§¥, г, [у, г, х, ]] + [г, х, [г, х, х, ]+ [х, г, [х, г, х1 ]} р
Уи = 8 -1 [I, х, р]+у§у, г ,[у, г , у ] + [I, х, [г, х, у, ]]+[х, у, [х, у, х, ]} (2) р
1ц = -1 [х, у, р]+v{[у, г, [у, г, г, Ц+ [г, х, [г, х, г, Ц+ [х, у, [х, у, ^ ]} р
х(а, Ь, о,,), у (а, Ь, о,,), г (а, Ь, о,,) - координаты траектории жидкой частицы; , - время, квадратные скобки обозначают операцию взятия якобиана по переменным а, Ь, о; р - плотность жидкости; р - давление; V - вязкость; 8 - ускорение свободного падения (ось у направлена вертикально). Уравнение непрерывности (1) отражает то обстоятельство, что в начальный момент координаты жидких частиц равны соответственно а, Ь, о .
Систему уравнений (1), (2) следует дополнить граничными условиями. В случае свободных волн на бесконечно глубокой воде это требование непротекания на дне - у = 0 при Ь = -ж, и условие отсутствия вязких напряжений на свободной поверхности, которое можно записать так:
Tknk = (~Ро )ni; b = 0;
- xo(YaZc -YcZa) + yo(ZaXc -XaZc) + ^o(XaYc -XcYa) n = I 2 2 0 2 ; (3)
V (YaZc - YcZa ) + (ZaXc - XaZc ) + (XaYc - XcYa )
KJaa (1 + Yc) - 2YaYcYca + Ycc (1 + Ya2)
(1+ у а +У2Г '
где Т ^ - тензор вязких напряжений; р0 - постоянное внешнее давление; а - коэффициент поверхностного натяжения, п - внешняя нормаль к свободной поверхности у = у (а, 0, о) ; К - ее средняя кривизна. Выражения для компонент тензора Т^ получаются, если их представление в эйлеровых координатах записать в якобианной форме, затем перейти к новым переменным (по правилу "перемножения" якобианов) и учесть уравнение (1):
Тхх =- Р + ^р[х, ,у, г ];
Туу =- р - 2vр[уt, х, г ];
Т22 =- р - 2vр[It ,у, х]; Тху = Тух =vр(у, ,у, 1]-[х, , х, г]);
Тж = Тгх = vP\Zt ,у, 1]-[х, ,у, х]); Туг = Тгу = vрP[Уt, х,у]-[1,, х, г ]).
Представим все входящие в уравнения движения функции в виде ряда по малому параметру крутизны волны е = кЛ, к - волновое число, Л - амплитуда волны:
х = а + + й(е2); у = Ь + ел! + 0(е2); 2 = о + е^1 + 0(е ); р = р0-р8Ь + ер + 0(е ).
Подстановка этих соотношений в систему (2),(3) даст искомые уравнения для неизвестных функций.
Линейные волны
Система уравнений в первом порядке по малой крутизне волны имеет вид:
Ъа +Л16 +Cic = 0;
Slff = P1a - ЯЛи + vA^1t;
-1 л (4)
л и, =-р Р1ь- g^1b +vA^f;
Cut =-Р~1Рю - № +vAC1,,
где A - лапласиан по переменным Лагранжа. Граничные условия для нее на свободной поверхности ( b = 0 ) запишутся следующим образом:
Лш = 0;
Л1,с +с1,ь = 0; (5)
- Р1 + 2vprh,b = а(Л1аа + Л1СС)-
Будем искать решение методом разделения переменных, полагая
^ = Re A(b)ent sin kacos mc; Л = Re B(b)ent cos ka cos mc; C1 = Re C(b)ent cos ka sin mc; p1 = Re H (b)ent cos ka cos mc.
(6)
A/F -
Функции А, В, С, Н и постоянная п считаются комплексными. Проводя несложные алгебраические выкладки, придем к уравнению для функции А :
,2 2\ п 1 (у 2 22 2 п V „ к + т )+ — А +\к + т \к + т +— А = 0.
у} \ у)
Полагая А = ехр 1Ь, получим для I биквадратное уравнение, решением которого будут соотношения
¡2 = к2 + т1 = М2; ¡1 = к2 + т2 + - = И2.
у
Волновые возмущения должны спадать по мере погружения вглубь (при Ь ^ -ж ), поэтому функцию А следует выбрать в таком виде:
А = авМЬ + рвИЬ; М,ЯеN > 0. (7)
Выражения для функций В,С, Н находятся из системы (4) и выглядят так:
О М МЬ . О ИЬ п т МЬ . -г, ИЬ
В =--ав + 5в ; С =—ав + ке ;
к к (8)
Н = Я(¿М + п2)авМЬ - pgдвNЬ, к
а, р, 5, Я - комплексные постоянные, удовлетворяющие условию
кр + N5 + тЯ = 0, (9)
которое является следствием уравнения непрерывности (1). Для определения конкретного
вида этих постоянных подставим выражения (7), (8) в граничные условия (5). В итоге получим три уравнения для определения четырех постоянных:
2Ма + Жр- кб = 0;
—
2М — а + Ж - тб = 0; к
1 [п2 + 2vnM2 + 8М + р-1аМ3 ]а - [8 + 2vnN + р-1аМ2 ]з = 0. С учетом соотношения (9) условие совместности этой системы запишется так:
(п + 2vM2)2 + (8М + р-1аМ3) = ^2М3 N. (10)
Введем обозначения
ш2 = 8М + р_1аМ3; vM- = 6; п + 2vM2 = , (11)
ш
с их учетом выражение (10) можно свести к следующему уравнению
(52 +1)2 = 1663(5 -6). (12)
Оно совпадает с уравнением для плоских поверхностных волн в вязкой жидкости, только роль волнового числа в двумерной волне теперь играет величина полного волнового числа
М = 4к2 + т , где к - продольное волновое число, а — - поперечное (пространственный период модуляции). Значение — = 0 отвечает плоской волне. Выражение для параметра ю совпадает с частотой гравитационно-капиллярных волн, распространяющихся по поверхности идеальной жидкости.
Постоянные, входящие в решение для линейных волн, определяются равенствами
2МЖ _ 2М3 _ 2—Ш р =--о-га; б =-;-^-а; К =--^-^-а. (13)
М2 + Ж2 к(М + Ж2) к(М + N )
Соотношение (6)-(8), (12), (13) дают полное решение задачи о стоячих пространственных линейных волнах на поверхности вязкой жидкости. Величина а, входящая в выражения для остальных констант, задает амплитуду волны В0, которая
В = \В\1 Л= Яе- пМа
0 = ВЬ=0 = V к (М2 + N2). В качестве параметра крутизны волны е выступает величина е = МВ0 .
Анализ дисперсионного уравнения
Конкретные значения постоянных (13) находятся по известному решению дисперсионного уравнения (12). Оно дает возможность по заданным значениям вязкости жидкости, частоты ш и волнового числа М, определяющих величину параметра 6, найти величину числа 5 , а значит, и выражение для п , входящее в решение (6).
Проведем качественный анализ уравнения (12). Графиком функции, задаваемой его левой частью, служит вогнутая кривая с минимумом, равным 1 при 5 = 0. Левая же часть уравнения определяет прямую. Она пересекает ось ординат ниже нуля и имеет положительный угол наклона, тангенс которого равен 1663. Две эти линии либо пересекаются (при достаточно большом 6 ), и тогда уравнение (12) имеет два действительных корня, либо не пересекаются: в этом случае дисперсионное уравнение не имеет действительных корней. Оба
действительных корня, как следует из соотношений (10), (11), удовлетворяют условию Re N > 0 (смотри (7)).
Проверим теперь, удовлетворяют ли этому условию комплексные корни уравнения (12). Пусть ^ = у + ¿у2, где у, у2 - действительные числа. При условии, что у2 ^ 0, уравнение (12) эквивалентно следующей системе:
(у2 -у22 +1)2 = 4^2 +1663(у -6); (14)
^(у2 -у22 +1) = 463. (15)
Знак величины Яе N совпадает со знаком выражения Яе(у2 +1)= (у2 - у\ +1), поэтому, исходя из равенства (14) заключаем, что условию Яе N > 0 соответствуют только корни с положительной реальной частью. С другой стороны, правая часть выражения (14), очевидно, всегда больше нуля. Но это означает, что величина у удовлетворяет условиям
у1 <-2
6"'
1 +
2
Л
1
1 + 7 62
У
у1 > 2
6
С
1+% -1 62
У
При наличии только комплексных корней два из них имеют отрицательную реальную часть, а другая пара комплексно сопряженных решений (12) имеет положительную реальную часть. В случае существования пары действительных корней (для них у2 = 0, у > 0) два других комплексных корня соответствуют у < 0, а значит, не удовлетворяют условию Яе N > 0.
Длинноволновое приближение
Будем рассматривать волновые возмущения, для которых выполняется соотношение
М-1 »у/ 2у/ю = А-1. (16)
Величина А характеризует толщину (очень тонкого) пограничного слоя, а само неравенство (16) отражает тот факт, что длина волны существенно превышает этот характерный масштаб. Из соотношения (16) следует, что 6 << 1. При выполнении этого условия в уравнении (12) можно пренебречь правой частью, и его решением будут значения у = ± / или в размерных переменных
п = -2уМ2 ± /ю . (17)
Отсюда заключаем, что величина декремента равна Яе п = -2уМ 2, а частота колебаний равна частоте гравитационно-капиллярных волн в идеальной жидкости ю . Знак перед частотой определяет фазу колебаний и может быть любым; выберем для определенности плюс. С учетом неравенства (16) выражения (13) перепишутся так
_ 2ЫШ _ 2/уМ3 п 2ышШ ,, ,,Л ,10Л
Р«-а; 8«--а; Я «-а; N «А(1 + /). (18)
ю юк юк
Из последнего равенства следует, что вторые слагаемые в выражениях (7), (8) для функций
А, В, С, Н существенны только в поверхностном пограничном слое толщины А-1. Величина
Яе а будет равна кМ 1В0. Выражения (6)-(8) с учетом соотношений (16)-(18) описывают
пространственную слабозатухающую стоячую поверхностную волну.
Формулы для проекций вектора завихренности в лагранжевых переменных имеют вид
= [7,, 7, X] + [7, ,7, X] = е(С - ) + 0(82); = [7,, 7,7] + [X,, X,7] = 8(^с - С) + 0(82); = [7, ,7, 7] + [X,, X, 7] = 8Ка - ^ ) + 0(82).
Подставляя в них представления для траекторий жидких частиц (6), получим для компонент завихренности 0Х1, 071, первого порядка теории возмущений следующие представления:
2
2
у
у
2
2
2mMn2a
Q - = - Re-т-^-exp(Nb + nt) cos ka sin mc; Q , = 0;
vk (M 2 + N2) y1
_ „ 2Mn2a ч • ,
£2 - = Re-т-s-exp (Nb + nt)sinka cos mc.
v(M 2 + N2)
В длинноволновом приближении эти выражения примут вид
Qxl = 2mM< Reia exp[ A(1 + i)b + iwt]cos ka sin mc; k
Qzl = 2M< Re ia exp [ A(1 + i)b + i< t] sin ka cos mc.
Случай m = 0 соответствует плоской волне [9]. Полагая A = 0, v = 0, получим решение для пространственной потенциальной стоячей волны [10].
Библиографический список
1. Ламб, Г. Гидродинамика / Г. Ламб. - М.: Гостехиздат, 1947.
2. Филипс, О.М. Динамика верхнего слоя океана / О.М. Филипс. - Л.: Гидрометеоиздат, 1980.
3. Stokes, G.G. On the theory of oscillatory waves // Trans. Camb. Phil., 1847. V. 8. P 441-455; reprinted in Math. Phys. P. 1. 314-326.
4. Longuet-Higgins, M.S. Mass transport in water waves // Phil. Trans. R. Soc. Lond., 1953. A 245. 535-581.
5. Longuet-Higgins, M.S. Mass transport in the boundary layer at a free oscillating surface // J. Fluid Mech., 1960. 8. P. 293-306.
6. Weber, J.E. Effect of the air on the drift velocity of water waves / J.E. Weber, E. Forland // J. Fluid Mech., 1990. 218. Р. 619-640.
7. Weber, J.E. Mass transport induced by surface waves in viscous rotating fluid // Free surface flows with viscosity. Computational Mechanics Publications: Southampton UK and Boston USA, 1997.
8. Абрашкин, А.А. Вихревая динамика в лагранжевом описании / А.А. Абрашкин, Е.И. Якубович. - М.: Физматлит, 2006. - 175 с.
9. Abrashkin, A.A. Nonlinear Gravitational Waves on the Surface of Viscous Fluid: Lagrange Approach / A.A. Abrashkin, Yu.P. Bodunova // Physics of Wave Phenomena. 2010. V. 18. No. 4. Р.251-255.
10. Сретенский, Л.Н. Теория волновых движений жидкости / Л.Н. Сретенский. - М.: Наука, 1977. - 815 с.
Дата поступления в редакцию 29.04.2011
A.A. Abrashkin, Yu.P. Bodunova
SPATIAL STANDING WAVES ON THE SURFACE OF VISCOUS FLUID
Analytical solution for spatial standing linear gravitationally-capillary waves on the surface of viscous fluid is determined. Lagrange coordinates are used. The qualitative analysis of dispersion equation is done. The motion of fluid particles is studied in long-wave approximation. The components of vorticity are founded.
Key words: water waves, viscous fluid, vorticity, Lagrange variables.