CC BY
15
2021 Математика и механика № 69
УДК 531.38
DOI 10.17223/19988621/69/8
Вин Ко Ко, А.Н. Темнов
ТЕОРЕТИЧЕСКОЕ ИССЛЕДОВАНИЕ ЭФФЕКТОВ КОЛЕБАНИЙ ДВУХ НЕСМЕШИВАЮЩИХСЯ ЖИДКОСТЕЙ В ОГРАНИЧЕННОМ ОБЪЁМЕ
Исследованы нелинейные колебания двухслойной жидкости, полностью заполняющей ограниченный объём. Используя две основные несимметричные гармоники, возбуждаемые в двух взаимно перпендикулярных плоскостях, исследованы дифференциальные уравнения нелинейных колебаний поверхности раздела двухслойной жидкости. В результате построены области неустойчивости вынужденных колебаний двухслойной жидкости в цилиндрическом баке, а также области параметрического возбуждения для различных соотношений плотностей жидкостей. Для построения областей неус-тойчивостей при приближенном решении нелинейных дифференциальных уравнений использован метод Бубнова - Галеркина.
Ключевые слова: цилиндрический бак, двухслойная жидкость, нелинейные колебания, область неустойчивости, гидродинамические коэффициенты.
Нелинейная теория движения ограниченного объема жидкости со свободной поверхностью, а также с поверхностью раздела двух жидкостей, представляющая собой особый раздел механики, используется при решении ряда практических задач. Обеспечение устойчивого полета современных и перспективных летательных аппаратов ракетно-космической техники и достижение точности управления ими невозможно без тщательного описания динамических процессов, происходящих в сложной механической системе.
Литературу по волновым движениям жидкостей, состоящих из слоев разных плотностей и имеющих отношение к настоящей работе, условно можно представить двумя направлениями.
К первому отнесем работы, связанные с исследованием волновых движений жидкостей разных плотностей, занимающих открытую область пространства. Здесь прежде всего необходимо отметить основополагающие работы Сретенского Л.Н. [1], Ландау Л.Д. [2], в которых изложены основные сведения и методы исследования колебаний двух жидкостей. Из иностранных работ отметим фундаментальные статьи по колебаниям двухслойной жидкости Thorpe S.A., Camassa R., M. La Rocca [3-5].
Среди современных работ, связанных с рассматриваемой тематикой нелинейных колебаний двухслойной жидкости, следует отнести работы [6, 7], в которых проведено экспериментальное исследование профиля двумерно гравитационных волн и показано, что для таких волн имеют место вторичные циркуляционные течения, пронизывающие всю жидкость. В работе Порубова А.В. [8] установлено, что двухмерные нелинейные внутренние волны двухслойной жидкости могут быть описаны при помощи двухмерного обобщенного уравнения Гарднера. Басинским К.Ю. в работе [9] рассмотрена нелинейная задача о распространении волны на свободной поверхности вязкой жидкости в плоском случае. Закономер-
ности реализации гравитационного нелинейного волнового движения в двухслойной жидкости с конечной глубиной верхнего слоя исследованы в работе [10].
К второму направлению отнесем работы, в которых изучаются волновые движения жидкостей в ограниченном объеме неподвижного или подвижного твердого тела. Из последних работ отметим работы Калиниченко В. А., [11, 12], в которых рассмотрен эффект влияния верхнего слоя вязкой жидкости на колебания двухслойной жидкости в прямоугольном сосуде. В статье Мерзлякова А.В., Крюковой Е.А. [13] рассмотрена плоская задача о колебаниях в прямоугольном баке двухслойной жидкости, разделенной жёсткой горизонтальной проницаемой перегородкой, а в работах [14, 15] - подобная задача в круглом цилиндрическом баке с упругой перегородкой.
В работе Науменко В.В., Стрельниковой Е.А. [16] показано, что использование потенциала двойного слоя позволяет учесть эффект краевого вихревого шнура и дает полное соответствие математической модели реальному процессу в задаче о свободных колебаниях пологой оболочки в идеальной несжимаемой жидкости.
В работах [17-21] рассмотрены малые колебания двух- и трёх-слоистых жидкостей в полостях различной формы для случаев подвижного и неподвижного твердого тела. При проведении экспериментов с жидкостями, полностью заполняющими круглый цилиндрический бак, вблизи основного резонанса было замечено вращательное движение слоёв жидкостей, подобное движению свободной поверхности однородной жидкости.
Особенности линейных и нелинейных колебаний однородной жидкости, частично заполняющей полость подвижного и неподвижного твердого тела, рассмотрены в [22-25].
Целью данной статьи является исследование нелинейных колебаний поверхности раздела двухслойной жидкости и получение теоретической интерпретации наблюдаемого эффекта.
1. Постановка задачи
Рассмотрим осесимметричный сосуд произвольной формы, полностью заполненный двумя жидкостями. Введем систему координат Oxyz , с началом в точке O на невозмущённой поверхности раздела жидкостей (рис. 1). Жидкости плотности Pj и р2 предполагаются идеальными и несжимаемыми. Обозначим через h и h2 глубины каждого слоя жидкости при отсутствии вращательного движения поверхности раздела. Систему координат Oxyz расположим так, чтобы в невозмущенном положении механической системы тело - жидкости ось Ox была перпендикулярна невозмущенной поверхности раздела жидкостей Г0.
Смоченные поверхности полости обозначим через S(i)(i = 1,2), а возмущенную поверхность раздела жидкостей - через Г (см. рис. 1).
Уравнение возмущенной поверхности раздела можно представить в виде, разрешенном относительно координаты x:
z = x - f (y, z, t) = 0. (1)
В предположении отсутствия вихревого движения в каждой жидкости сформулируем задачу о нелинейных колебаниях поверхности раздела жидкостей, полностью заполняющих полость подвижного твердого тела, совершающего поступательное движение по закону: U (t) = S cos rat (см. рис. 4).
Рис. 1. Система координат и основные обозначения для тела с двухслойной жидкостью Fig. 1. Coordinate system and basic designations for a body with a two-layer fluid
С учетом допущений постановка задачи состоит из уравнения Лапласа, условий непротекания на смачиваемых поверхностях, а также кинематического и динамического условий на возмущенной поверхности раздела и имеет вид
V2Ф(1) = 0, в т , У2Ф(2) = 0, в т2;
dФ
(1)
Р 2-
dv
дФ(2)
~дГ'
-=0, на S1
dФ'
(2)
dv
-=0, на S2
dФ(1) dФ
(2)
dv
dv
на Г;
(2) (3)
дФ'
,(1) ^
-Рг
dt
+1 [р2 (VФ(2) )2 - Р1 (VФ(1))2 ] + (Р2 - Р1) U• г
= (Р1 - Р2 )ё •r на Г •
(4)
Г |г = Г |г + 1Х/ , V - внешняя нормаль к соответствующей границе области, занимаемой жидкостью.
При решении подставленной задачи будем предполагать, что квадрат отклонений ] возмущенной поверхности раздела от невозмущенной имеет порядок малости е. При вычислениях будем удерживать линейные, квадратичные, кубичные члены / , /2, /3, пренебрегая величиной е2.
Рассмотрим сначала случай свободных движений жидкостей в неподвижной полости. С этой целью представим потенциалы скоростей каждой жидкости в виде следующей суммы:
Ф(k) (х, y, z, t) = X a (t)B(k) (x, y, z), k = 1,2,
(5)
i=1
где Ф(к^ - потенциалы скоростей верхней и нижней жидкостей, В(к-1 - функции координат верхней и нижней жидкостей, а,- - обобщенные координаты волновых
S
2
движений жидкостей на поверхностях раздела /-й гармоники. Здесь и в дальнейшем суммирование по / производится по натуральному ряду чисел от единицы до бесконечности. Верхние индексы параметров (1) и (2) относятся соответственно к верхней и нижней жидкостям.
Предположим далее, что известна некоторая ортогональная на области Г0
система функций / (у, х), составляющих вместе с постоянной величиной полную систему функций. Отклонение поверхности раздела жидкостей разложим по системе функций / :
/ = С)/ (У, х). (6)
Подставив (5) в уравнения Лапласа (2) и в кинематические граничные условия (3), получим краевые задачи
V2В'к) = 0 , в т,
ёВ/,)
ё V
= 0, (к = 1,2),
СВ(г)
ё V
ссвР
ё V
= / / N .
(7)
Представим функции В(к) в виде разложения по параметрам а1 до второго порядка включительно
В/к) = В(к) аЦк) а^ к Ц
+....
(8)
] ] к
у ' ^¡¡к зависят только от пространственных координат и не
где функции ВО), В(к), В(к) зависят от времени.
Используя разложения вектора нормали и производных по нормали на возмущенной поверхности раздела жидкостей в ряд Тейлора и выражая все функции, входящие в кинематическое условия через их значения на невозмущенной поверхности раздела, исходную нелинейную задачу сведем к последовательному решению линейных краевых задач. С этой целью подставим разложения (6), (8) в (7) и, приравнивая выражения при одинаковых степенях параметров аг , после некоторых преобразований получим следующие краевые задачи для введенных выше функций:
V2 В^) = 0, в т,:
В?
ё V
■ = 0, на Б.
ёВ(0) ёВ^
г 0 = г 0 - = /г на Г
ё V
ё V
2 (к) В)
V2В к) = 0 , в тк, —И— = 0, на ; у к сv
ёВт ёВ(1)
-С- = = V' •(/ V 'В/0)) = V' • (/ V 'В®) на Г0
V2 Вк = 0, в Тк
ёвЦ к>
г к
ё V
= 0, на Бк ;
ёВ$_ СВ/к
ё V
1 = IV'-(/ V 'в, ) + к V ' в(к)) +1V'-(//к ¿V 2 1 гк ¡ 2 1 дх
(9) (10) (11) (12)
) на Г . (13)
г=1
Б
к
Далее воспользуемся динамическим граничным условием (4). Умножим уравнение (4) на функцию ^ и проинтегрируем по невозмущенной поверхности Г0. В результате получим бесконечную систему обыкновенных дифференциальных уравнений для обобщенных координат аг- по индексу /:
d 2 a
(ц(2) + g(N<2)2 -NPK + ££( j - j)а,
dt
+IK N(2) - N (k) j k
j k
dа j dak
dt dt
(2) - N(D
"UK j - j)akai
d 2a j
dt2 d2a,-
j k г
+ - j )ai jk = 0 С = 1'2'3'-)
dt2
(14)
где
И(k) = Pk J Bk)#)dГо, N}k)2 = pk J (#))2dГо;
(15)
j =Pk J (Bk) + Л(2) ^) j) d Го;
(16)
dB(k) dB( k) 1 d2 B(k) j =PkJ [ B« dB^+(dxr+2 f(2) ^d—^) fi) fj( j )]d Го
Г0 dxk dxk
(17)
Njk = Pk J (j) + 2VB(0)-VB®)f®dГо;
(18)
N(k = Pk i
J [2B(1 + f
(k)
(k) dj dx2
+vb (k) - VBkk) + VBk - VB(k) +
jk kl
(k ).VR( k).
k о
jl
+fk) A(VB(k) -VBkо))fk)]dГо 2 ax,
(19)
2. Определение гидродинамических коэффициентов уравнений движения для случая цилиндрической полости в неподвижном твёрдом теле
Рассмотрим сначала неподвижное твердое тело (т.е. U (t) = о) с цилиндрической полостью, целиком заполненной жидкостями. Пусть полость имеет форму кругового цилиндра радиуса r = го (см. рис. 2). Введем цилиндрическую систему координат x,r,п , связанную с декартовой x,y, z следующими формулами: x = x , y = r cos п, z = r sin п.
Для исследования поставленной задачи выделим две основные несимметричные гармоники, возбуждаемые в двух взаимно перпендикулярных плоскостях и определяемые обобщенными координатами и формами a( = a( (i = 1,2): a1 = a,
a2 = P , f = fa , f2 = f .
о
о
о
а
Рис. 2. Обозначения и системы коор- Рис. 3. Формы поверхностей раздела при возбу-
динат для случая прямого кругового ждении основных и вспомогательных гармоник
неподвижного цилиндра с двухслой- а, р
ной жидкостью Fig. 3. The shapes of interfaces in the course of
Fig. 2. Designations and coordinate generation of the fundamental and auxiliary
systems for the case of a straight circular harmonics а, p fixed cylinder with a two-layer fluid
В цилиндрической системе координат форму основного тона колебаний поверхности раздела жидкостей представим в виде, аналогичном представлению свободной поверхности жидкости при нелинейных колебаниях [23] (рис. 3).
/a(k)(y,z) = ф^япn, /ß(k)(y,z) = ф^Г^n , (k = l,2).
Представим искомые функции в виде
Bak0) = y( k )( X, r )sin n; Bß 0) = y( k )( X, r ) cos n;
\(k )<
r(k )
,(k ),
(20)
(2l)
= Y) - Y2) cos2n; = B^a) = Y(2k) sin 2n; B^^) = Y0) + Y2) cos 2n (22)
Краевая задача (9) совпадает с аналогичной задачей для малых колебаний жидкостей, поэтому воспользуемся результатами работ [17, 23] и запишем функцию в виде y(k) (х, г) = X%> (x)Ynm (r), k = 1,2,
X (1)( x) = Ch(knm (X - hl)) . X (2) (x) = ch(knm (X + h2)) . Y (r) = Jm (knmr) (23)
Xnm (x) =--k-h(k h ) ; Xnm (x) =~T-h(k , ) ; Ynm (r) = J~7-) , (23)
knm Sh(knmh1) knm Sh(knm"2) Jm (knmr0)
knm = 4nm / r0, 4nm является n-м корнем производной от функции Бесселя m-го порядка. Здесь индексы n, m принимают значения 1.1; 1.0; 1.2.
Имея решения краевых задач (9), нетрудно вычислить гидродинамические коэффициенты линейных уравнений движения тела с жидкостями (15), соответст-
Г
0
вующие случаю малых перемещений частиц жидкостей:
= Р:5„ , <2 = РЛ,, Цо(1) = Р1 ^, До(2) = Р28„ = п,
"п "ц 2^11
г® = -Л,, ЩкцН,), ^(2) = £„ Ш^А). (24)
Решения краевых задач (10) - (13), учитывающих нелинейные эффекты, будем искать в виде соответствующих разложений в ряд Фурье по функциям Бесселя нулевого и второго порядков, т.е.
^(х,г) = Xс(пк2хП2(х)¥пт(г), т = 0,2 . (25)
п
Коэффициенты этих разложений определяются следующими интегралами:
1 Го 1 г 2
^ = IК* + (7-кп)^12^ог^ ; = ; (26)
^ = ТТ^3°С^ц'2 -(-2 + к2)^2]!^ ; жП2)2 = . (27)
""11 Ж п0 0 Г 2Ь п2
Собирая значения функций (26) - (27), вычислим значения параметров ц1( к) и
(к)
ц 2 , определяющих нелинейность волновых движений жидкости, учет которых позволяет описать явление вращения поверхности раздела жидкостей:
Ц(к) = Рк4Ш/|[6^11'2 +[4-Зк2 '1 7п2]7п +ХГц' -к20Ш) +
4П1 0 ^ ( г ' п " п0
2с(к) С 2 1 1
+ Е^^п' +^-2 - к22 ^ УпГп]уигёГ; (28)
Ц2к) = Рк ^ {К - (7+к" ) ^ +Е к ^ +
+ С4 - кп221 ^^к^г. (29)
11 0
+ чг2
Перепишем уравнения для обобщенных координат волновых движений жидкостей на поверхности раздела в виде
а+ст2а + ё1 (а2 а+ а2 а+ ар р+ ар2) + ё2(Р2 а + 2рар-ара-2ар2) = 0; (30) р+ст2р + й?1(р2 р+а2 р + ара+рр2) + ^2(р2 р+2аар-ара-2ра2) = 0, (31)
где
о2 = 8(<2 -<2) ^ = (ц(2) -ц(1)) ^ = (ц22) -Ц21)) ц*=(ц(2) Ц(1)) * > "1 _ * , "2 _ * 5 И-0 —^И-0 М-0 > •
Ц0 Ц0 Ц0
3. Исследование устойчивости вынужденных и параметрических колебаний жидкостей в простейшем случае
Пусть полость совершает заданные движения в направлении оси Oz по закону U(t) = S cos rat, где U(t) - смещение заданного движения сосуда (см. рис. 4), S и га - амплитуда и частота возмущения.
*■ x
u(t) = S cos rot 4 *
O
Pi
"■4
hi
/ / 777 77777 TTTTTTTTTTTTTTTT
Рис. 4. Цилиндрический сосуд, совершающий поступательное движение Fig. 4. Cylindrical tank in translational motion
Соответствующая этому случаю система нелинейных уравнений приобретает вид
Z1(a,ß) = а+ с2а+^(а2 а+ а2 а + aßß+ aß2) + + d2 (ß2 а+ 2ßaß- aßa- 2aß2) - га2P cos rat = 0;
(32)
L2 (a,ß) = ß+ c2ß + di(ß2 ß+ a2 ß + aßa+ ßß2) + + d2 (ß2 ß+ 2aaß- aßa- 2ßa2) = 0,
(33)
где P = -
Sk
* ■
До
X = X(2) -X(i), k(k) = pknj r27i1dr:
Pk nro
§2i
k = 1,2.
Эта система приближенно описывает вынужденные и параметрические возбуждаемые колебания поверхности раздела жидкостей. В том случае, когда параметрические колебания в системе не возникают (в = 0), вынужденные колебания описываются нелинейным дифференциальным уравнением
L1(a) = a+ с2 a+d1(a2 a+ a2 a) = ra2P cos rat.
(34)
Приближенное решение этого уравнения найдем методом Бубнова - Галеркина, представив решения в виде
a(t) = a0 + ^ (ak cos krat + ak sin krat).
k=1
где ak и ak - неизвестные постоянные.
k
Г
г,
о
z
h
2
Удержав в (35) только основные гармоники
a(t) = A cos rat+A sin rat, (36)
получим уравнение
(с2 -1)A - m A3 = P, A = 0, (37)
где m1 = d1 /2, с2 =c2/ra2.
Уравнение (37) используется для определения амплитуд вынужденных колебаний двухслойной жидкости в зависимости от параметра P и ra. Положив в (37) P = 0 , для определения зависимости амплитуд свободных колебаний жидкостей от частоты получим формулу
(с2 -1) - mA2 = 0. (38)
В дальнейшем при определении установившихся режимов движения жидкостей и границ областей их устойчивости используются основные слагаемые в разложениях типа (35).
Остановимся на вопросе об устойчивости периодического решения (36). Ответ на него состоит в том, чтобы найти те значения параметров S и ra , при которых установившийся режим
a(t) = A cos rat, p = 0, (39)
описываемый системой нелинейных уравнений (32) - (33), физически реализуем. С этой целью наряду с движением (39), которое было принято за невозмущенное, рассмотрим также близкие к нему движения
a(t) = a(t) + 4(t), P(t) = P(t) + n(t), (40)
которые будем называть также возмущёнными, а величины 4(t) и n(t) - возмущениями.
В соответствии с общей теорией устойчивости составим уравнения в вариациях, соответствующие заданной системе нелинейных уравнений (32), (33). Подставляя выражение (39) в систему уравнений (32), (33) и учитывая, что a(t) является частным решением (34), получаем уравнения возмущенного движения в виде
(d a2 +1) 4 + 2d1aa^+ (с2 + 2d1 aa+ d1 a2 )4 + F1 (4,p, 4, p, 4, P) = 0; (41)
(1 + d2) n+ 2d 2 aan+ (с2 + ca a+ k3a2 )n + F2 (4, p, 4, p, 4, P) = 0, (42)
где F1 и F2 - функции, содержащие возмущения и их производные в степенях выше первой. Оставляя в (41) и (42) только линейные члены, с учетом (39) приходим к уравнениям в вариациях:
L3 (4) = (p + q cos 0t) 4- £j 4 sin 0t + (у - 5 cos 0t)4 = 0 ; (43)
L4 (n) = (p + q cos 0t) n- en sin 0t + (y - 5 cos 0t )n = 0, (44)
носящих названия уравнений первого приближения. Здесь приняты следующие обозначения:
q = (d1A2)/2, p = 1 + q, e1 = d1oA2, у = (с2 -(d1A2)/2)ra2, 5 = (3d1ra2A2)/2, p = 1 + q , q = (d2A2)/2, у = (с2 -(d2A2)/2)o>2,
5 = ю2A2(2d1 - 3d2), e = d2oA2, 0 = 2ю. (45)
Перейдем к построению областей неустойчивости решений уравнений (43), (44). Рассмотрим сначала уравнение (43) относительно функции ^(t), характеризующее возмущение периодического решения a(t) = A cos rat. Амплитуда A этого решения определяется уравнением (37). Необходимо установить, какая пара значений A и ю, удовлетворяющих уравнению (37), приводит к устойчивым решениям и какая пара к неустоичивым. Имея в виду построение основной области неустойчивости решения уравнения (43), представим его в виде
0t 0t
^(t) = a1 cos— + b1 sin—. (46)
Из уравнения Бубнова - Галеркина
4п 4п
Т 0t Т 0t
J L3(a)cos—dt = о, j L3(a)sin—dt = о, (47)
о 2 о 2
для определения границ области неустойчивости получим следующие отношения:
(с2 -1) - 3mA = о; (48)
(с2 -1) - mA2 = о. (49)
Уравнение (49) совпадает с уравнением скелетной линии (38). На рис. 5 эта линия соответствует кривой АВС. Уравнению (48) соответствует кривая AMO. Сопоставляя уравнение (48) с уравнением амплитудно-частотных характеристик (37), нетрудно заметить, что устойчивая ветвь резонансной кривой RNM отделяется от неустойчивой точкой M, в которой амплитудные кривые имеют вертикальную касательную. На устойчивой части ветви амплитудно-частотной характеристики, находящейся слева от линии АВС, производная должна быть обязательно положительной. В точке M она стремится к бесконечности, а на неустойчивой левой ветви она отрицательна. В области I, ограниченной кривыми AMO и АВС, решение a(t) = A cos rat неустойчиво. С физической точки зрения условие
устойчивости периодческого решения с периодом возмущающей силы означает, что амплитуда вынужденных колебаний с возрастанием внешней силы P увеличивается.
Рассмотрим далее уравнение (44) в вариациях относительно возмущения n(t) тривиального решения P(t) = о . Исследование решений уравнения (44) должно дать ответ на вопрос об устойчивости этого тривиального решения. В соотвествии с изложенным выше области неустойчивости уравнения (44) отвечают областям параметрически возбуждаемых колебаний (в Ф о), т.е. областям динамической неустойчивости режима движения (39). Для построения основной области неустойчивости положим
0t , • 0t
n(t) = a^os — + bjsin—. (5о)
Поставив выражения (50) в уравнения Бубнова - Галеркина
4п 4п
Т et Т et
Г Z4(n)sin—dt = 0, Г Z4(n)cos—-dt = 0, (51)
0 0 для определения границ областей динамической устойчивости получим следующие уравнения:
с2 = m2 Л2 +1; (52)
с2 = m Л2 +1, (53)
где m2 = (4d2 - d1). На рис. 5 уравнению (52) соответствует кривая линия ADE,
а уравнению (53) - кривая линия ABC. Следовательно, области неустойчивости решений уравнений (43), (44) непрерывно переходят одна в другую.
Рис. 5. Амплитудно-частотные характеристики и области неустойчивости вынужденных колебаний жидкостей в цилиндрическом баке при возбуждении основных гармоник а для случая рх=0 при S = 0.005, h = h2 = 2
Fig. 5. Frequency-response characteristics and instability regions of forced oscillations of fluids in a cylindrical tank in the course of generation of the fundamental harmonics а for
the case of p = 0 at S = 0.005, hj = h2 = 2
В области II рис. 5 и 6, ограниченной кривыми ABC и ADE, решение Р(/) = 0 неустойчиво. Вторая устойчивая ветвь амплитудно-частотной характеристики KLD примыкает к области II справа и отделяется от нее точкой D. В области динамической неустойчивой II установившийся режим, если он сушествует, описывается нелинейной системой уравнений (32), (33).
Рис. 6. Амплитудно-частотные характеристики и области неустойчивости вынужденных колебаний жидкостей в цилиндрическом баке при возбуждении основных гармоник а (а) для pj = 0.5, (b) для pj = 0.8 при S = 0.005,
h = h2 = 2
Fig. 6. Frequency-response characteristics and instability regions of forced oscillations of fluids in a cylindrical tank in the course of generation of the fundamental harmonics а for
the case of pj = (a) 0.5 and (b) 0.8 at S = 0.005, hj = h2 = 2
Рассмотрим теперь вопрос о построении установившихся режимов движения жидкостей, происходящих в основной области динамической неустойчивости II.
Предположим, что в области основного резонанса приближенное решение системы нелинейных уравнений (32) - (33) можно представить в виде
a(t) = A cos rat+A sin rat, P(t) = B cos rat + B sin rat. (54)
Применив к решению этой системы метод Бубнова - Галеркина, для постоянных A, A, B и B получаем следующие алгебраические соотношения:
(с2 -1) A - m A3 - m2 AB2 = P; (55)
(с2 -1)-m1 B2 -m2A2 = о; (A = о, B = о). (56)
Исключая В из (56) и подставляя результат в (55), находим уравнение для определения амплитудно-частотных характеристик в области II
(ст2 -1) A - m4A3 = m5P; где m4 = 2d2 , m5 = m1 /(m1 - m2). (57)
Решение (57) отвечает наблюдаемому в эксперименте режиму вращения поверхности раздела жидкостей. Соответствующие ему резонансные кривые представлены на рис. 5 и 6 линиями FGQ.
Из всех полученных соотношений на рис. 5 построены амплитудно-частотные характеристики (АЧХ) и области неустойчивости вынужденных колебаний жидкостей при отсутствии верхней жидкости р1 =о, которая полностью совпадала с результатом задачи для одной жидкости [22].
На рис. 6 представлены амплитудно-частотные характеристики (АЧХ) и области неустойчивости вынужденных колебаний двух жидкостей при различных соотношениях плотностей.
Области I и II являются областями неустойчивости вынужденных колебаний жидкостей, происходящих в плоскости действия возмущающей силы. В области II эта неустойчивость обусловлена неустойчивостью тривиального решения P(t) = о, т.е. возможно параметрическое возбуждение обобщенной координаты в . Линия FGQ соответствует вращательному движению поверхности раздела жидкостей, наблюдаемому в эксперименте.
Заключение
Теоретически исследованы нелинейные эффекты, возникающие в результате взаимодействия жидкостей с жестким сосудом, совершающим гармонические колебания. Наиболее интересным с практической стороны является случай колебаний жидкостей в окрестности самой низкой частоты собственных колебаний поверхности раздела. Здесь наблюдается ряд характерных существенно нелинейных особенностей движения жидкостей, среди которых можно указать на зависимость частоты колебаний от амплитуды, ограниченность амплитуд колебаний, подвижность узловых линий поверхности раздела и возникновение своеобразного вращения слоев жидкостей в некотором диапазоне частот возмущающей силы.
ЛИТЕРАТУРА
1. Сретенский Л.Н. Теория волновых движений жидкости. М.: Наука, 1977. 815 с.
2. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Т. VI: Гидродинамика. М.: Наука, 1986. 735 с.
3. Thorpe S.A. On standing internal gravity waves of finite amplitude // J. Fluid Mech. 1968. V. 32. No. 3. P. 489-528.
4. Camassa R., Hurley M.W., McLaughlin R.M., Passaggia P.-Y., Thomson C.F.C. Experimental investigation of nonlinear internal waves in deep water with miscible fluids // Journal of Ocean Engineering and Marine Energy. 2018. V. 4. P. 243-257.
5. La Rocca M., Sciortino G., Adduce C., Boniforti M.A. Experimental and theoretical investigation on the sloshing of a two-liquid system with free surface // Physics of Fluids. 2005. No. 17. P. 062101.
6. Калиниченко В.А., Секерж-Зенькович С.Я. Экспериментальное исследование вторичных стационарных течений в поверхностных волнах Фарадея // Изв. РАН. Механика жидкости и газа. 2008. № 1. С. 141-148.
7. Калиниченко В.А., Секерж-Зенькович С.Я. Экспериментальное исследование волн Фарадея максимальной высоты // Изв. РАН. Механика жидкости и газа. 2007. № 6. С. 103-110.
8. Порубов А.В. О локализации двумерных нелинейных внутренних волн в двухслойной жидкости // Журнал технической физики. 2005. Т. 75(7). С. 48-51.
9. Басинский К.Ю. Нелинейные волны на поверхности слоя вязкой жидкости // Изв. Са-рат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2015. Т. 15. № 3. С. 322-329.
10. Григорьев А.И., ФедоровМ.С., Ширяева С.О. Волновое движение в поле силы тяжести на свободной поверхности и на границе стратификации слоисто-неоднородной жидкости. Нелинейный анализ // Изв. РАН. Механика жидкости и газа. 2010. № 5. С. 130-140.
11. Kalinichenko V.A. Regularization of barotropic gravity waves in a two-layer fluid // Fluid Dynamics. 2019. V. 54. No. 6. P. 761-773.
12. Kalinichenko V.A. Effect of an upper layer of viscous liquid on breaking surface gravity waves // Journal of Physics: Conference Series. 2019. V. 1301. P. 012017.
13. Мерзляков А.В., Крюкова Е.А. Свободные колебания идеальной жидкости в прямоугольном сосуде с горизонтальной проницаемой перегородкой // Вестник Томского государственного университета. Математика и механика. 2019. № 60. C. 107-118. DOI: 10.17223/19988621/60/8.
14. Пожалостин А.А., Гончаров Д.А. Свободные осесимметричные колебания двухслойной жидкости с упругим разделителем между слоями при наличии сил поверхностного натяжения // Инженерный журнал: наука и инновация. МГТУ им. Н. Э. Баумана. Электрон. журн. 2013. № 12(24).
15. Гончаров Д.А. Динамика двухслойной жидкости, разделенной упругой перегородкой с учетом сил поверхностного натяжения // Наука и образование, МГТУ им. Н.Э. Баумана. Электрон. журн. 2013. № 11.
16. Науменко В.В., Стрельникова Е.А. Потенциал двойного слоя в задаче о свободных колебаниях пологой оболочки в идеальной несжимаемой жидкости // Исслед. по теор. пластин и оболочек. 1990. № 22. С. 122-133.
17. Вин Ко Ко, Темнов А.Н. Колебания дискретно-стратифицированных жидкостей в цилиндрическом сосуде и их механические аналоги // Вестник МГТУ им. Баумана. Сер. «Естественные науки». 2015. № 3. С. 57-69.
18. Вин Ко Ко, Темнов А.Н. Экспериментальное и теоретическое исследование колебаний твердого тела со слоистой жидкостью // Инженерный журнал: наука и инновации, МГТУ им. Н.Э. Баумана. Электрон. журн. 2018. № 04. С. 1-13.
19. Вин Ко Ко, Темнов А.Н. Теоретическое и экспериментальное исследования колебаний твёрдого полуцилиндра, имеющего полость, заполненную слоистой жидкостью // Инженерный журнал: наука и инновации, МГТУ им. Н. Э. Баумана. Электрон. журн. 2019. № 05.
20. Вин Ко Ко, Темное А.Н. Колебания вязкой трехслойной жидкости в неподвижном баке // Инженерный журнал: наука и инновации, МГТУ им. Н.Э. Баумана. Электрон. журн. 2019. № 07.
21. Win Ko Ko, Temnov A.N. Experimental and theoretical studies of oscillations of stratified fluid // IOP Conference Series: Materials Science and Engineering 468 (2018) 012031.
22. Лукоеский И.А. Введение в нелинейную динамику твердого тела с полостями, содержащими жидкость; отв. ред. В.А. Троценко; Ин-т математики АН УССР. Киев: Наук. Думка, 1990. 296 с. ISBN 5-12-001308-2.
23. Нариманов Г.С., Докучаев Л.В., Лукоеский И.А. Нелинейная динамика летательного аппарата с жидкостью. М.: Машиностроение, 1977. 208 с.
24. Микишев. Г.Н Экспериментальные методы в динамике космических аппаратов. М.: Машиностроение, 1987. 248 с.
25. Микишев. Г.Н., Рабинович Б.И. Динамика твердого тела с полостями, частично заполненными жидкостью. М.: Машиностроение, 1968. 532 с.
Статья поступила 25.12.2019
Win Ko Ko, Temnov A.N. (2021) A THEORETICAL STUDY OF OSCILLATIONS OF TWO IMMISCIBLE FLUIDS IN A LIMITED TANK. Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mekhanika [Tomsk State University Journal of Mathematics and Mechanics]. pp. 97-113
DOI 10.17223/19988621/69/8
Keywords: cylindrical tank, two-layer fluid, nonlinear oscillations, instability region, hydrodynamic coefficients.
In the paper, the nonlinear oscillations of a two-layer fluid that completely fills a limited tank are theoretically studied. To determine any smooth function on the deflected interface, the Taylor series expansions are considered using the values of the function and its normal derivatives on the undisturbed interface of the fluids. Using two fundamental asymmetric harmonics, which are generated in two mutually perpendicular planes, the differential equations of nonlinear oscillations of the two-layer fluid interface are investigated. As a result, the frequency-response characteristics are presented and the instability regions of the forced oscillations of the two-layer fluid in the cylindrical tank are plotted, as well as the parametric resonance regions for different densities of the upper and lower fluids. The Bubnov-Galerkin method is used to plot instability regions for the approximate solution to nonlinear differential equations.
At the final stage of the work, the nonlinear effects resulting from the interaction of fluids with a rigid tank that executes harmonic oscillations at the interface of the fluids are theoretically studied.
Ko Ko WIN (Candidate of Physics and Mathematics, Bauman Moscow State Technical University, Moscow, Russian Federation). E-mail: win.c.latt@gmail.com
Aleksandr N. TEMNOV (Candidate of Physics and Mathematics, Bauman Moscow State Technical University, Moscow, Russian Federation). E-mail: antt45@mail.ru
REFERENCES
1. Sretenskiy L.N. (1977) Teoriya volnovykh dvizheniy zhidkosti [Theory of wave motions in a fluid]. Moscow: Nauka.
2. Landau L.D., Lifshitz E.M. (1987) Course of Theoretical Physics. Volume 6: Fluid Mechanics. Oxford: Pergamon Press.
3. Thorpe S.A. (1968) On standing internal gravity waves of finite amplitude. Journal of Fluid Mechanics. 32(3). pp. 489-528.
4. Camassa R., Hurley M.W., McLaughlin R.M., Passaggia P.Y., Thomson C.F.C. (2018) Experimental investigation of nonlinear internal waves in deep water with miscible fluids. Journal of Ocean Engineering and Marine Energy. 4. pp. 243-257. DOI: 10.1007/s40722-018-0119-9.
5. La Rocca M., Sciortino G., Adduce C., Boniforti M.A. (2005) Experimental and theoretical investigation on the sloshing of a two-liquid system with free surface. Physics of Fluids. 17(062101). DOI: 10.1063/1.1922887.
6. Kalinichenko V.A., Sekerzh-Zen'kovich S.Ya. (2008) Experimental investigation of secondary steady flows in Faraday surface waves. Fluid Dynamics. 43(1). pp. 125-131. DOI: 10.1134/s10697-008-1014-1.
7. Kalinichenko V.A., Sekerzh-Zen'kovich S.Ya. (2007) Experimental investigation of Faraday waves of maximum height. Fluid Dynamics. 42. pp. 959-965. DOI: 10.1134/ S0015462807060117.
8. Porubov A.V. (2005) O lokalizatsii dvumernykh nelineynykh vnutrennikh voln v dvukhsloynoy zhidkosti [On localization of two-dimensional nonlinear internal waves in a two-layer fluid]. Zhurnal tekhnicheskoy fiziki. 75(7). pp. 48-51.
9. Basinskiy K.Yu. (2015) Nelineynye volny na poverkhnosti sloya vyazkoy zhidkosti [Nonlinear waves on the surface layer of a viscous liquid]. Izvestiya Saratovskogo universiteta. Novaya seriya. Seriya Matematika. Mekhanika. Informatika - Izvestiya of Saratov University. New series. Series: Mathematics. Mechanics. Informatics. 15(3). pp. 322-329.
10. Grigoryev A.I., Fedorov M.S., Shiryaeva S.O. (2010) Wave motion in a gravity field on the free surface and stratification interface of a fluid with stratified inhomogeneity. Nonlinear analysis. Fluid Dynamics. 45. pp. 793-802. DOI: 10.1134/S0015462810050121.
11. Kalinichenko V.A. (2019) Regularization of barotropic gravity waves in a two-layer fluid. Fluid Dynamics. 54(6). pp. 761-773. DOI: 10.1134/S0015462819060065.
12. Kalinichenko V.A. (2019) Effect of an upper layer of viscous liquid on breaking surface gravity waves. Journal of Physics: Conference Series. 1301(012017). DOI: 10.1088/17426596/1301/1/012017.
13. Merzlyakov A.V., Kryukova E.A. (2019) Svobodnye kolebaniya ideal'noy zhidkosti v pryamougol'nom sosude s gorizontal'noy pronitsaemoy peregorodkoy [Free oscillations of an ideal fluid in a rectangular vessel with a horizontal permeable membrane]. Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mekhanika - Tomsk State University Journal of Mathematics and Mechanics. 60. pp. 107-118. DOI: 10.17223/19988621/60/8.
14. Pozhalostin A.A., Goncharov D.A. (2013) Svobodnye osesimmetrichnye kolebaniya dvukhsloynoy zhidkosti s uprugim razdelitelem mezhdu sloyami pri nalichii sil poverkhnostnogo natyazheniya [Free axisymmetric oscillations of two-layered liquid with the elastic separator between layers in the presence of surface tension forces]. Inzhenernyy zhurnal: nauka i innovatsiya - Engineering Journal: Science and Innovation. 12(24).
15. Goncharov D.A. (2013) Dinamika dvukhsloynoy zhidkosti, razdelennoy uprugoy peregorodkoy s uchetom sil poverkhnostnogo natyazheniya [Dynamics of two-layer liquid, divided by an elastic baffle considering surface tension forces]. Nauka i obrazovanie: nauchnoe izdanie MGTU im. N.E. Baumana - Science and Education: Scientific Periodical of the Bauman MSTU. 11.
16. Naumenko V.V., Strel'nikova E.A. (1990) Potentsial dvoynogo sloya v zadache o svobodnykh kolebaniyakh pologoy obolochki v ideal'noy neszhimaemoy zhidkosti [The potential of a double layer in the problem of free vibrations of a flat shell in an ideal incompressible fluid]. Issledovaniyapo teoriiplastin i obolochek. 22. pp. 122-133.
17. Win Ko Ko, Temnov A.N. (2016) Kolebaniya diskretno-stratifitsirovannykh zhidkostey v tsilindricheskom sosude i ikh mekhanicheskie analogi [Oscillations of immiscible liquids in a stationary cylindrical vessel and their mechanical analogs]. Vestnik MGTU im. N.E. Baumana. Seriya «Estestvennye nauki» - Herald of the Bauman Moscow State Technical University. Series: Natural Sciences. 3(66). pp. 57-69. DOI: 10.18698/1812-33682016-3-57-69.
18. Win Ko Ko, Temnov A.N. (2018) Eksperimental'noe i teoreticheskoe issledovanie kolebaniy tverdogo tela so sloistoy zhidkost'yu [Experimental and theoretical study of vibrations of a rigid body with a layered liquid]. Inzhenernyy zhurnal: nauka i innovatsii - Engineering Journal: Science and Innovation. 4. pp. 1-13.
19. Win Ko Ko, Temnov A.N. (2019) Teoreticheskoe i eksperimental'noe issledovaniya kolebaniy tverdogo polutsilindra, imeyushchego polost', zapolnennuyu sloistoy zhidkost'yu [Theoretical and experimental studies of vibrations of the solid half-cylinder having a cavity filled with a layered fluid]. Inzhenernyy zhurnal: nauka i innovatsii - Engineering Journal: Science and Innovation. 5.
20. Win Ko Ko, Temnov A.N. (2019) Kolebaniya vyazkoy trekhsloynoy zhidkosti v nepodvizhnom bake [Oscillations of a three-layer viscous fluid in a stationary tank]. Inzhenernyy zhurnal: nauka i innovatsii - Engineering Journal: Science and Innovation. 7.
21. Win Ko Ko, Temnov A.N. (2018) Experimental and theoretical studies of oscillations of stratified fluid. IOP Conference Series: Materials Science and Engineering. 468(012031). DOI: 10.1088/1757-899X/468/1/012031/meta.
22. Lukovskiy I.A. (1990) Vvedenie v nelineynuyu dinamiku tverdogo tela s polostyami, soderzhashchimi zhidkost' [Introduction to nonlinear dynamics of a solid with cavities filled with a fluid]. Kyiv: Naukova Dumka.
23. Narimanov G.S., Dokuchaev L.V., Lukovskiy I.A. (1977) Nelineynaya dinamika letatel'nogo apparata s zhidkost'yu [Nonlinear dynamics of an aircraft with a fluid]. Moscow: Mashinostroenie.
24. Mikishev G.N. (1987) Eksperimental'nye metody v dinamike kosmicheskikh apparatov [Experimental methods in spacecraft dynamics]. Moscow: Mashinostroenie.
25. Mikishev G.N., Rabinovich B.I. (1968) Dinamika tverdogo tela s polostyami, chastichno zapolnennymi zhidkost'yu [Dynamics of a solid body with cavities partially filled with a fluid]. Moscow: Mashinostroenie.
Received: December 25, 2019
