CC BY
25
УДК 532.591
DOI 10.21685/2072-3040-2016-4-5
К. Ю. Басинский
ОБРАЗОВАНИЕ ВОЛН НА ПОВЕРХНОСТИ ЖИДКОСТИ ПРИ ОБТЕКАНИИ КРУГОВОГО ЦИЛИНДРА
Аннотация.
Актуальность и цели. Цель работы - исследование процесса образования волн при обтекании потоком идеальной жидкости бесконечной глубины кругового цилиндра, находящегося под поверхностью жидкости.
Материалы и методы. Для решения задачи действие цилиндра заменено действием диполя, который заменял бы собой цилиндр при его обтекании неограниченным потоком. Применен метод теории возмущений - метод малого параметра. Разложением неизвестных функций в ряды по степеням малого параметра исходная задача разбита на задачи по порядку малости.
Результаты. С точностью второго приближения получены выражения для потенциала скорости жидкости и ординаты свободной поверхности.
Выводы. Исследовано влияние различных параметров задачи на форму свободной поверхности жидкости.
Ключевые слова: поток идеальной жидкости, обтекание цилиндра, метод последовательных приближений, поверхностные волны.
K. Yu. Basinskiy
OCCURRENCE OF WAVES ON THE SURFACE OF FLUID FLOWING IN A CIRCULAR CYLINDER
Abstract.
Background. The aim of the work is to study the process of formation of waves of ideal fluid of infinite depth flowing around a circular cylinder located below a liquid surface.
Materials and methods. To solve the problem the action of a cylinder was replaced by the action of a dipole, which would replace a cylinder with its unlimited flow stream. The author applied a method of the perturbation theory - the small parameter method. By the decomposition of unknown functions in powers of the small parameter the initial task was broken down into tasks in the order of smallness.
Results. With the accurace of the second approximation the author obtained expressions for the potential velocity of liquid and the ordinates of free surface.
Conclusions. The researcher has studied the effect of various parameters of the problem on the form of free surface of the liquid.
Key words: ideal fluid flow, flowing around a cylinder, method of successive approximations, surface waves.
Введение
Решение задачи об образовании волн при обтекании кругового цилиндра, находящегося под поверхностью жидкости, обычно находится либо для волн малой амплитуды [1-3], когда условия на поверхности жидкости сводятся к линейным в предположении, что скорость потока значительно превосходит скорость возмущения жидкости, либо с применением численных методов [4-6]. При аналитическом решении нелинейных задач о поверхност-
ных волнах как правило используется метод малого параметра, подробное описание которого можно найти, например, в [2, 7]. Однако рассматриваемые при этом задачи обычно отличаются параметрами среды, а влияние на волновое движение погруженных в жидкость тел не исследуется. В данной работе метод малого параметра применяется для решения нелинейной задачи об образовании волн, вызванном обтеканием кругового цилиндра потоком идеальной жидкости.
1. Задача об обтекании цилиндра
Пусть под поверхностью идеальной жидкости бесконечной глубины расположен цилиндр радиуса а, на который набегает поток жидкости со скоростью с (рис. 1). Введем систему координат так, что ось Ох совпадает с невозмущенной свободной поверхностью, а ось Оу направлена вертикально вверх. Центр цилиндра находится в точке (0, -И), а его образующие перпендикулярны плоскости Оху. Предположим, что движение жидкости потенциально. Кроме того, в силу постановки задачи движение является установившимся. Тогда компоненты вектора скорости возмущенного движения жидкости можно представить в виде
Эф Эф
и = с + —, V = —, Эх Эу
где функция ф(х, у) - потенциал скорости возмущения основного потока, который удовлетворяет уравнению Дф = 0 .
Рис. 1. Геометрия задачи
На свободной поверхности должно выполняться кинематическое условие, которое для случая плоских установившихся волн запишется в виде
=(сЗЖ у=*, (1)
Эу ^ Эх) ах
здесь £(х) - ордината свободной поверхности.
Также вдоль свободной поверхности давление р сохраняет постоянное значение, равное значению атмосферного давления:
р = Ра, у = £.
Тогда для точек свободной поверхности интеграл Бернулли запишется в виде
c +
Эф Эх
2 ( Эф ^
Эу
+ g % = const, y =
(2)
Исключим из условий (1) и (2) функцию , продифференцировав второе условие по переменной x и подставив в него вместо производной функции ее выражение из первого условия. Получим
Эф Э2ф ( Эф ^ Эф Э2ф
c + — I + lc + — I—-— + .
Эх) Эх2
Эф
Эх ) Эу ЭхЭу ° Эу
= 0, y = %.
(3)
На поверхности цилиндра должна равняться нулю нормальная производная потенциала скорости:
^ = 0. Эп
При заглублении обе составляющих возмущения скорости должны стремиться к нулю:
Эф Эх
= 0,
Эф ¥
= 0,
y
2. Линейная задача
Применение метода последовательных приближений требует введения малого параметра, значение которого определяется характеристиками волнового движения. В данной задаче его величину можно задать, исходя из условия малости скорости возмущения относительно скорости потока c. Для этого потребуется сначала решить задачу в линейном приближении, когда условия (1) и (3) принимают вид
^^ ЭШ
у = 0. (4)
Эф ¥
2 Эф
= c—, c2-±- +.
а% ёх'
Эх2
= 0, Эу ,
Действие цилиндра заменим действием диполя, который заменял бы собой цилиндр при его обтекании неограниченным потоком [2]. Из этого допущения следует, что вблизи центра цилиндра потенциал Ш должен принимать следующий вид [8]:
ф
a2cх
+ (у + h )
2
Поэтому функцию Ш представим в следующем виде:
ф(х, у ) = Ф(х, у ) + -
a 2cх
a 2cх
2
(5)
х2 + (у + И) х2 +(у - И)
Здесь функция Ф(х, у) является потенциальной и удовлетворяет согласно второму условию (4) следующему граничному условию:
2
х
2 Э2Ф ЭФ 4са2ghx
с ~гГ+gir = 1-, у = 0. (6)
Эх2 Эу (x2 + h2 )2
Функцию Ф(х, y) будем искать в следующем виде:
Ф = J Aeky sin kxdk, (7)
0
соответствующем установившимся волновым движениям. При этом условие затухания волнового движения при заглублении в выражении (7) уже учтено. Подставляя этот вид функции Ф в условие (6), получим равенство для определения функции A:
| AI k-g I к sin kxdk = --
2
4а ghx
0 4 с (x2 + h2 )2
Применим к этому равенству обратное преобразование Фурье:
AÍk -JL1k = J 4aghx„ sinkxdx = -kh.
C ' 0 с
(x2 + h 2 )2
Отсюда получается выражение для функции А:
= 2са2g е-¿И А =( - ск ) ,
подстановка которого в интеграл (7) позволяет определить выражение для потенциала Ф:
^ о 2 Г к (y -h) sin kx 4na^g g(y h) gx Ф = 2са g I e w '---dk =--— eC cos-s (8)
0 (g- Ck) с с2
4^2g (у -h)
Тогда выражение для искомой функции ф согласно (5) принимает вид
4па2g ~(у-И) gx а2сх а2сх
ф =--- ес + .
с2 х2 + (у + И)2 х2 + (у - И)2
Подставив его в первое условие (4), получим выражение для формы ординаты свободной поверхности далеко за обтекаемым цилиндром:
2 - мИ 2
„ 4па м 2 . мх 2а И
£ = ---е с + -2, (9)
с с2 х2 + И2
которое совпадает с полученным в [2].
3. Метод малого параметра
Для приближенного решения исходной нелинейной задачи воспользуемся методом малого параметра. Малый параметр выберем, исходя из предположения малости скорости возмущения жидкости относительно скорости потока, т.е. малости отношения
Эф Эх
2 ^ Эу
Далеко за обтекаемым цилиндром выполняется приближенное соотно-
шение
шах,
Эф
Эх
2
Эу
+
4ка
2 2
с с4 "
Поэтому в качестве малого параметра выберем величину
2 2
£ = ■
а
4
(10)
Тогда выражение (5) для функции ф запишется в виде
5
£С
ф(х, у ) = Ф(х, у ) + —
х2 +(у + h )2 х2 +(у - h )2
(11)
Неизвестные функции Ф и Е, будем искать в виде степенных рядов по малому параметру е:
Ф = £Ф1 + £2Ф2 + £3Ф3 + ... , ^ = £^1 +£2^2 +£3^3 +...
(12)
Сведя условия (1) и (3) разложением в ряд Маклорена входящих в них функций к условиям на фиксированной поверхности у = 0 [2, 7] и подставив выражение (11) и ряды (12) в уравнение и граничные условия, можно выписать задачи первого и второго порядка малости по параметру е .
Задача первого порядка малости имеет вид
ЭФ
Эу
1 d L 1 = с—— + -
4c5hx
dx
:(х2 + h 2)
ЭФ1
АФ1 = 0,
2
2 Э Ф] ЭФ]
с2-^ +" 1
ЭФ1
4с5кх
Эх2
Эу
г( х2 + h 2 )2
у = 0,
Эх
= 0,
ЭФ
Эу
1 = 0, у .
с
Для второго порядка малости задача примет следующий вид:
ДФ 2 = 0,
ЭФ2 _ d$2 э2Ф1 $ ЭФj d$1
■ _ c
Эу dx Эу2
$1 +■
Эх dx
у _ 0,
(13)
2 Э2Ф2 ЭФ2 „ -2 + g—- _-$r
Эх2
Эу Эу
2 Э2Фт ЭФт c2-г1 + "-1
Эх 2
Эу
- 2c
ЭФ1 Э2Ф1 Эх Эх2
ЭФ1 Э2Ф, -с—1-1 + -
4с6Их Э2Ф1 4c6h (h2 - 3х2) эф1
Эу ЭхЭу „2 (х2 + h2 )2 ЭхЭу g2 (х2 + h2 )3 Эу
48Лх (х2 - h2) 16c11h2х (3х2 - h^
+-*-Л1 +--$1, у _ 0,
:(х2+h2)
ЭФ2
(х2 + h2 )5
ЭФ2
_ 0, —2 _ 0, у
(14)
Эх ' ду
Задача первого порядка малости решена выше, и функции Ф^ и ^ согласно выражениям (8), (9) и (10) запишутся в виде
^ 4яг3 4(у-h) gx „ 4^c -jgh . gx 2c Ф1 _--ec cos^, $1 _--e c sin^- +
(х2 + h 2)
. (15)
Подставив выражения (15) в условия (13) и (14), получим задачу для определения неизвестных функций Ф2 и ^2 , которая решается аналогично линейной задаче, вследствие чего ее решение здесь опускается. Функции Ф2 и ^2 определяются в виде
. 4кс3 Ч(у_к) . 2&х 0 4(у_к) &х
Ф 2 =--ес 8ш—+ Вес О08^т,
g с2 с2
_ 2ф
4 КС 2 2 &Х =-(с _ 2к)е с 008—+
& с2
B--
8nc4
^ gh
:(х2 + h 2)
__ 2 2 gh
2 . gx 8л c 2
e c sin—2 +-e c ,
c2 g
где
B_
2c9 л
3g
4,3
3 + % + 3I + 21 2
c I c J I c
+ 32c3I 2лg c2 | g 1 c2h gh
+11 - ^ л
2 gh
2 gh
Здесь Ег(т) - интегральная показательная функция [9]. Собирая вместе решения задач первого и второго порядка малости, получим приближенное решение исходной задачи:
, 4лс3 "Г(у -h) gx с5
ф = -£-ес cos^- + £—
g с2 g2
+ (у + h )2 х2 +(у - h )
+ £
( 4лс3 ¥(у-h) . 2gx D 4(-h) gx ^ --е с sin^^ + Вес cos^r-
g с2 с2
-gh 4 4лс ^ с2 g^ , 2с
-е с sin^r- +
2 „2
(х2 +h2)
+ £
-2gh
4лс 2 2 рх -(с-2л)е с cos^2 +
В--
8лс4
(х 2 + h 2)
2 2 gh
2 . рх 8л с 2
е с sin-2— Н--е с
с2 g
Для иллюстрации на рис. 2-4 приведены графики формы свободной поверхности при различных значениях радиуса цилиндра (рис. 2), скорости набегающего потока (рис. 3) и глубине расположения цилиндра (рис. 4).
Рис. 2. Форма свободной поверхности при разном радиусе цилиндра (с = 5 м/с, И = 1 м)
'С
^с = 5 м/с /
/уУ 65 \ 70 / 5 \
с = 3 м/с X X
Рис. 3. Форма свободной поверхности при разной скорости потока (а = 0,5 м, И = 1 м)
Рис. 4. Форма свободной поверхности при разной глубине расположения цилиндра (с = 5 м/с, a = 0,5 м)
Из рис. 2-4 видно, что увеличение радиуса цилиндра приводит к возрастанию амплитуды волны. Большей скорости потока соответствуют волны большей длины. Возрастание глубины расположения цилиндра приводит в волновым возмущениям меньшей амплитуды.
Заключение
Таким образом, методом малого параметра получено решение задачи об образовании волн при обтекании кругового цилиндра потоком идеальной жидкости. Проиллюстрировано влияние различных параметров задачи на форму свободной поверхности жидкости.
Список литературы
1. Елизаров, А. М. Капиллярно-гравитационные волны при циркуляционном обтекании подводного цилиндра в канале конечной глубины / А. М. Елизаров, К. В. Кириллин, С. И. Филиппов // Ученые записки Казанского университета. Сер.: Физико-математические науки. - 2011. - Т. 153, № 1. - С. 147-154.
2. Сретенский, Л. Н. Теория волновых движений жидкости / Л. Н. Сретенский. -М. : Наука, 1977. - 816 с.
3. Хаскинд, М. Д. О поступательном движении тел под свободной поверхностью тяжелой жидкости конечной глубины / М. Д. Хаскинд // Прикладная математика и механика. - 1945. - Т. 9, № 1. - С. 67-78.
4. Стурова, И. В. Численные расчеты в задачах генерации плоских поверхностных волн / И. В. Стурова. - Красноярк, 1990. - 48 с. - (Препринт № 5 / ВЦ СО АН СССР).
5. Терентьев, А. Г. Численные методы в гидродинамике / А. Г. Терентьев, К. Е. Афанасьев. - Чебоксары : Изд-во Чувашского гос. ун-та, 1987. - 94 с.
6. Romate, J. E. The numerical simulation of nonlinear gravity waves / J. E. Romate // Engineering Analysis with Boundary Elements. - 1990. - Vol. 7, № 4. - P. 152-166.
7. Алешков, Ю. З. Теория волн на поверхности тяжелой жидкости / Ю. З. Алешков. - Л. : Изд-во Ленингр. ун-та, 1981. - 196 с.
8. Кочин, Н. Е. Теоретическая гидромеханика / Н. Е. Кочин, И. А. Кибель, Н. В. Розе. - М. : Физматгиз, 1963. - Ч. 1. - 584 с.
9. Градштейн, И. С. Таблицы интегралов, сумм, рядов и произведений / И. С. Градштейн, И. М. Рыжик. - М. : Физматгиз, 1963. - 1100 с.
References
1. Elizarov A. M., Kirillin K. V., Filippov S. I. Uchenye zapiski Kazanskogo universiteta. Ser. «Fiziko-matematicheskie nauki» [Proceedings of Kazan University. Series: Physical and mathematical sciences]. 2011, vol. 153, no. 1, pp. 147-154.
2. Sretenskiy L. N. Teoriya volnovykh dvizheniy zhidkosti [The theory of wave motion of fluid]. Moscow: Nauka, 1977, 816 p.
3. Khaskind M. D. Prikladnaya matematika i mekhanika [Applied mathematics and mechanics]. 1945, vol. 9, no. 1, pp. 67-78.
4. Sturova I. V. Chislennye raschety v zadachakh generatsii ploskikh poverkhnostnykh voln [Numerical calculations in problems of flat surface waves generation]. Krasnoyark, 1990, 48 p. (Preprint № 5 / VTs SO AN SSSR).
5. Terent'ev A. G., Afanas'ev K. E. Chislennye metody v gidrodinamike [Numerical methods in hydrodynamics]. Cheboksary: Izd-vo Chuvashskogo gos. un-ta, 1987, 94 p.
6. Romate J. E. Engineering Analysis with Boundary Elements. 1990, vol. 7, no. 4, pp. 152-166.
7. Aleshkov Yu. Z. Teoriya voln na poverkhnosti tyazheloy zhidkosti [The theory of waves on the gravity solution's surface]. Leningrad: Izd-vo Leningr. un-ta, 1981, 196 p.
8. Kochin N. E., Kibel' I. A., Roze N. V. Teoreticheskaya gidromekhanika [Theoretical hydromechanics]. Moscow: Fizmatgiz, 1963, part 1, 584 p.
9. Gradshteyn I. S., Ryzhik I. M. Tablitsy integralov, summ, ryadov i proivedeniy [Tables of integrals, sums, orders and products]. Moscow: Fizmatgiz, 1963, 1100 p.
Басинский Константин Юрьевич
кандидат физико-математических наук, доцент, кафедра математического моделирования, Тюменский государственный университет (Россия, г. Тюмень, ул. Володарского, 6)
E-mail: [email protected]
Basinskiy Konstantin Yur'evich Candidate of physical and mathematical sciences, associate professor, sub-department of mathematical modeling, Tyumen State University (6 Volodarskogo street, Tyumen, Russia)
УДК 532.591 Басинский, К. Ю.
Образование волн на поверхности жидкости при обтекании кругового цилиндра / К. Ю. Басинский // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2016. - № 4 (40). -С. 51-59. БОТ 10.21685/2072-3040-2016-4-5
