УЧЕНЫЕ ЗАПИСКИ КАЗАНСКОГО УНИВЕРСИТЕТА. СЕРИЯ ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ
2020, Т. 162, кн. 4 С. 426-440
ISSN 2541-7746 (Print) ISSN 2500-2198 (Online)
УДК 519.634
doi: 10.26907/2541-7746.2020.4.426-440
АСИМПТОТИКА МЕДЛЕННЫХ ДВИЖЕНИИ ПРЯМОУГОЛЬНОГО ЦИЛИНДРА В ЖИДКОСТИ ПОСЛЕ ОТРЫВНОГО УДАРА
М.В. Норкин
Южный федеральный университет, г. Ростов-на-Дону, 344 090, Россия
Аннотация
Рассмотрен процесс схлопывания каверны, образованной в результате отрывного удара прямоугольного цилиндра в идеальной, несжимаемой, тяжелой жидкости. В предположении, что скорость цилиндра мала, построены асимптотики основных характеристик удара. Возникающие при этом трудности связаны главным образом с тем, что динамика точек отрыва заранее неизвестна и зависит от малого параметра, которым является безразмерная скорость движения цилиндра. С помощью специальной замены переменных проблема сведена к исследованию задачи, в которой динамика точек отрыва соответствует главному приближению, не зависящему от указанного параметра. Это дает возможность определить второй член асимптотики, учитывающий нелинейные слагаемые в модели. В главном приближении сформулирована задача со свободной границей, которая в каждый фиксированный момент времени аналогична классической модели удара с отрывом. На основе первых двух членов асимптотики описан процесс схлопывания каверны с учетом подъема внутренней свободной границы жидкости и проведено сравнение с известными результатами.
Ключевые слова: идеальная несжимаемая жидкость, прямоугольный цилиндр, отрывной удар, динамика точек отрыва, число Фруда, число кавитации
Удар твердого тела, полностью погруженного в жидкость, приводит к отрыву частиц жидкости от его поверхности и к образованию присоединенной каверны за телом [1]. Форма каверны и конфигурация внешней свободной поверхности заранее неизвестны и подлежат определению в ходе решения задачи. При небольших скоростях движения цилиндра удается получить полное решение проблемы образования и схлопывания каверны, опираясь на задачи со свободной границей, близкие по своей структуре к классической модели удара [1], а также на асимптотический метод малого параметра. В работах [2, 3] изучалась динамика точек отрыва, вызванная медленными движениями кругового и прямоугольного цилиндров в жидкости после отрывного удара. Исследование этих задач проводилось с помощью специальной математической модели, основанной на полном пренебрежении квадратичными членами и возмущениями свободных границ жидкости. При этом зависимость точек отрыва от времени носила нелинейный характер. Предложенная модель была получена из физических соображений и подтверждалась хорошим согласованием с другими подходами. Было установлено, что при числах Фруда, меньших 0.5, линеаризованную модель можно использовать для изучения процесса схлопывания тонкой каверны и приближенного определения времени
Введение
схлопывания. Из физических соображений ясно, что данная модель должна соответствовать главному члену асимптотики при малых скоростях движения цилиндра. Однако соответствующий асимптотический анализ исходной нелинейной задачи не проводился. В предлагаемой статье строится асимптотика задачи по малому параметру, в качестве которого выступает безразмерная скорость движения цилиндра (фактически число Фруда). Подробно анализируются первые два члена асимптотики. На их основе изучается динамика точек отрыва и описывается процесс схлопывания тонкой каверны с учетом подъема внутренней свободной границы жидкости. Приводятся конкретные примеры решений.
В настоящей работе проводится некоторая аналогия со статьей [4], в которой рассматривается такая же физическая постановка задачи. При этом существенное отличие состоит в том, что в [4] асимптотика строится по другому параметру - малому времени, а скорость движения цилиндра не предполагается малой величиной. Отметим, что в работах, посвященных асимптотическому анализу на малых временах, главным образом рассматривались большие числа Фруда. При малых числах Фруда вследствие сильного увеличения скоростей движения точек отрыва асимптотики малых времен становятся не эффективными. Этим, в частности, объясняется интерес к задачам удара при малых числах Фруда.
Анализ полученных ранее результатов показывает, что при малых скоростях движения цилиндра процесс схлопывания каверны целесообразно рассматривать с учетом искусственной кавитации (в естественной ситуации, когда давление в каверне равно нулю, он происходит практически мгновенно). При давлении в каверне, близком к атмосферному, процесс схлопывания становится хорошо обозримым. Во многих случаях время схлопывания составляет сотые или десятые доли секунды. Именно такие времена в ударных задачах считаются малыми. При этом под мгновенным процессом удара следует понимать тысячные или десятитысячные доли секунды. Заметим, что исследование задачи на малых временах является актуальным в силу важности происходящих там физических явлений и процессов (образование волн, каверн, брызговых струй, существенное изменение реакции жидкости на тело и различные другие вопросы).
Среди близких работ отметим статью [5], в которой рассмотрен отрыв столба жидкости от днища цилиндрической емкости при ударе и изучено схлопыва-ние образующейся цилиндрической полости. В ней дан обзор современных работ, в которых теоретически и экспериментально подтверждается возникновение кавитации при ударе по цилиндрической емкости, частично заполненной жидкостью. Современные публикации по гидродинамическому удару свидетельствуют о том, что большой интерес проявляется не только к изучению динамических процессов, происходящих после удара, но и к самой классической модели удара. При этом наряду с кавитационным ударом также рассматривается и безотрывный удар [6-8]. В последних работах определяются присоединенные массы для некоторых новых конфигураций и приводятся данные экспериментов.
1. Постановка задачи
Рассматривается плоская задача об отрывном ударе прямоугольного цилиндра, полностью погруженного в идеальную, несжимаемую, тяжелую жидкость [1]. Предполагается, что после удара цилиндр движется вглубь жидкости с постоянной скоростью без вращения. Позади тела образуется присоединенная каверна, форма которой зависит от физических и геометрических параметров задачи. Требуется изучить процесс схлопывания каверны при малых скоростях движения цилиндра, соответствующих небольшим числам Фруда. Математическая постановка задачи,
S2(t)
Рис. 1. Постановка задачи
записанная в безразмерных переменных, в подвижной системе координат, связанной с цилиндром, имеет вид (рис. 1)
Др = 0, R e ü(t)
dn = -ny, R e Sn(t), lf + e2 lf + 0-5e2 (Vp)2 + b + С(x, T) - £2T - H - 0.5X = 0, R e Sl2 (t),
дт dy
edf + e = e^ + e^ Iе, R e S12(t), dy dx dx дт
lf + e2 d-f + 0.5e2 (Vp)2 + ^(x, т) = 0, R e Si{t), дт
(1) (2)
(3)
(4)
(5)
(6)
(7)
(8)
dy
e d¡~ = e + e-i ^ R e S2(t),
dy dx dx дт
d-¡T =0, y = -Hb + eV; ^ = 0, x = ±Hr, dy dx
p(x,y, 0) = fo(x,y), С(x, 0)=0, ^(x, 0)=0. Подробная запись условия непротекания (2) определяется равенствами dp
— = -1, y = b, с(т) < \x\ < 1; y = -b, -1 <x< 1; dy
dp =0, x = ±1, -b < y < b. dx
Начальное условие po(x,y), соответствующее классической модели удара, находится на основе решения смешанной краевой задачи теории потенциала с односторонними ограничениями на поверхности тела [1]:
Др0 = 0, R e П(0); p0 = 0, R e S2(0), dpn = -ny, po < 0, R e Sil(0), dfo > -ny, po = 0, R e Si2(0),
(9) (10) (11)
0 и 0 -LU ÎTO\
=0' y = -Hb; ~дх ' x = R' (12)
Безразмерные переменные, растянутое время т и функции Ф, п, С, подлежащие определению, вводятся по формулам
t' = — ■ t, X = ax, y' = ay, Ф' = амФ, p' = pv2p, v = л/да, v
t = ET, ф(X,y,T ) = £-1Ф(х,у,£Т ), Z (x,T ) = n(x,£T ), ф(х.Т ) = Ç(x,£T ).
где штрихами помечаются размерные величины.
Неподвижные координаты X, Y связаны с подвижными x, y соотношениями: X = x, Y = y + h(t). Ось y направлена вертикально вверх (против вектора силы тяжести), начало координат находится в центре прямоугольника. Картина течения жидкости предполагается симметричной относительно оси y.
В работе используются следующие обозначения: ф^,у,Ь) - потенциал скоростей абсолютного движения жидкости, записанный относительно подвижной системы координат; Q(t) - область течения жидкости; Sn(t) = {c(t) ^ \x\ ^ 1, y = 6}U{-1 ^ x ^ 1,y = — b} У {x = ±1, —b ^ y ^ b} - часть поверхности цилиндра, на которой не происходит отрыва частиц жидкости, Si2(t) = {y = = b + n(x,t), —c(t) < x < c(t)} - внутренняя свободная граница жидкости (граница каверны); x = ±c(t), y = b - координаты точек отрыва в момент времени t; y = H+C(x, t)+t - уравнение внешней свободной границы; V = (0, —vo) - скорость, приобретенная цилиндром в результате удара ( vo > 0 ); h(t) = —t - безразмерный закон движения цилиндра; 2а, 2b - длины горизонтальной и вертикальной сторон прямоугольника (после обезразмеривания для b сохраняется прежнее обозначение); R - радиус-вектор с координатами (x,y).
Задача (1)-(8) содержит следующие безразмерные параметры:
„ v0 2 Pa — Pc 2 Pa — Pc
E = Fro = , xo = 2-—, X = 2-,
лД— Pv2 Pv2
где Fro - число Фруда; xo - число кавитации (безразмерная разность давлений на внешней свободной поверхности и в каверне); x - число кавитации, связанное с фиксированной скоростью v ; pa - атмосферное давление; pc - давление в каверне; д - ускорение свободного падения; р - плотность жидкости. Скорость цилиндра vo предполагается малой величиной по сравнению с характерной скоростью v = л/д—.
В дальнейшем для сравнения с результатами работы [3] будут использованы следующие простые зависимости:
22 T = E to, x = E Xo,
где безразмерное время to определяется по аналогии с t (вместо v пишется vo ).
В точках пересечения внутренней свободной границы жидкости с поверхностью цилиндра (в точках отрыва) ставится условие Кутты- Жуковского, означающее, что скорость жидкости в этих точках должна быть конечной.
2. Асимптотическое решение задачи при малых е
Основные трудности, возникающие при решении поставленной задачи, заключаются в том, что функция е(Ь) = с(£т) > 0, определяющая в каждый фиксированный момент времени одну из двух симметричных точек отрыва, заранее неизвестна и содержит малый параметр £. Так как в главном приближении уравнение и краевые условия не зависят от £, то старшее приближение для точки отрыва,
обозначаемое через со(т), также не будет зависеть от этого параметра. Отметим, что с помощью прямого разложения искомых функций в асимптотические ряды по параметру е не удается найти второй член асимптотики. Для его определения необходимо сделать такую замену переменных, в результате которой функция с(ет) перейдет в функцию со(т), не зависящую от указанного параметра. Такую замену удобно сделать с помощью степенной функции:
где х = с(ет), у = Ь - координаты точки отрыва, соответствующие фиксированному (малому) числу е, а х = со(т), у = Ь - координаты точки отрыва, отвечающие главному приближению.
Отметим, что при замене (13) прямоугольник, расположенный в области х > > 0 (половина прямоугольного цилиндра), остается инвариантным. В самом общем виде (при замене х = х(х1,т), у = ух) преобразованная модель дана в [4, с. 110], где проводился асимптотический анализ исходной нелинейной задачи на малых временах. При этом нужно учесть, что в указанной работе по-другому вводились безразмерные переменные, поэтому малый параметр не стоит в граничных условиях при соответствующих членах. Кроме этого имеются некоторые отличия в обозначениях: ж ^ хх, Ь ^ т, е ^ Ь.
Отметим также, что при построении асимптотики по малому времени нужно определить только одну точку отрыва со , соответствующую начальному моменту времени. В отличие от этого, в настоящей статье возникает функция со(т), зависящая от времени. В силу симметрии задачи относительно оси у дальнейшие рассуждения проводятся при х > 0.
Решение преобразованной задачи будем искать в виде следующих асимптотических разложений (е ^ 0):
Подставляя (14)-(17) в преобразованную задачу, осуществляя с помощью формулы Тейлора перенос краевых условий с возмущенных участков границы на первоначально невозмущенный уровень и приравнивая коэффициенты при одинаковых степенях е, придем для определения функций и ^>2 к смешанным краевым задачам для уравнений Лапласа и Пуассона в области П(0). При дополнительном предположении, что точки отрыва монотонно сближаются, проинтегрируем динамические условия (в задачах для у>1 и ^2) по времени от 0 до т при любом фиксированном х € (0, со(т)). В результате для определения главного приближения придем к следующей задаче:
(13)
у>(х(хьт ),У1,т) = ¥>1(х1,у1 ,т) + е2^2 (хьуьт) +----,
С (х(хьт ),т) = е2С1(хьт) + е4С2(хьт) +----,
ф(х(х1, т),т) = е2ф1(х1 ,т) + еАф2(х1,т) + • • • , с(ет) = со(т)+ е2с1(т) + ••• .
(14)
(15)
(16) (17)
=0, Я € П(0)
у>1 = (н - Ь)т + 0.5хт, у1 = Ь, 0 < х1 < со(т),
(18) (19)
ду>1 ду1
1, у 1 = Ь, со(т) <х1 < 1; у1 = -Ь, 0 < х1 < 1, (20)
(22)
=0, ух = Н; —^ =0, у 1 = -Иъ.
(23)
При т = 0 приходим к задаче, полностью совпадающей с классической моделью удара (9)—(12), ее решение удовлетворяет следующему начальному условию: ^х(х1,у1, 0) = у>о(ж1,у1).
В силу того, что динамика точки отрыва заранее неизвестна, полученная в главном приближении задача является задачей со свободной границей. Неизвестная функция со (т) определяется из условия регулярности потенциала в точке отрыва, обеспечивающего выполнение условия Кутты- Жуковского. Кроме последнего условия должно также выполняться условие неотрицательности давления на смоченной поверхности тела, которое проверяется уже после решения задачи. В процессе решения задачи это условие проверить не удается, поскольку в краевое условие входит не давление, а интеграл от давления по времени.
В классической задаче гидродинамического удара формулируются граничные условия типа неравенств, которые равносильны условию Кутты-Жуковского. Эти неравенства оказываются очень полезными для численного решения задачи со свободной границей итерационными методами. В случае аналитического решения задачи, напротив, более удобным оказывается условие Кутты-Жуковского. Так как полученная в главном приближении задача исследуется численно, то естественно попытаться сформулировать односторонние ограничения на поверхности тела, которые хотя бы локально будут равносильны условию Кутты - Жуковского.
Сначала сформулируем неравенства в зонах контакта и отрыва, исходя из физических соображений. В реальных жидкостях давление в зоне контакта должно быть неотрицательным. Точнее, должно выполняться неравенство р ^ рс, где рс и 0 - давление насыщенных паров жидкости или газа. В случае искусственной кавитации и малых скоростей движения цилиндра это условие можно усилить, потребовав, чтобы давление в зоне контакта не опускалось ниже давления в каверне (р ^ рс). После интегрирования по времени придем к следующему неравенству:
В естественной ситуации (рс = 0) получаем совпадение с неравенством р ^ 0. Заметим, что при больших скоростях движения цилиндра и искусственной кавитации условие р ^ рс может нарушаться. Сделанные выводы подтверждаются в работах, посвященных асимптотическому анализу исходной нелинейной задачи на малых временах (соответствующий обзор дан в [4]).
В зоне отрыва формулируется приближенное условие
означающее, что жидкие частицы не могут входить внутрь твердого тела. При его выводе предполагается, что внутренняя свободная граница мало отличается от поверхности тела, и, следовательно, жидкие частицы, находящиеся на ней, фактически лежат на верхней стороне прямоугольника. Точность этого условия зависит от справедливости предположения о малости возмущения внутренней свободной границы жидкости.
Теперь рассмотрим вопрос о связи условий (24), (25) с условием Кутты-Жу-ковского. Перенесем начало координат в точку отрыва (Х1 ^ Х1 + со(т), у1 ^ у1 + + Ь) и введем в рассмотрение новую гармоническую функцию , связанную с
0.5^т - ¥>1(хь у1,т) + у>о(х1, у1) + (Н - у1)т > 0.
(24)
(25)
соотношением
= -У1 + 0.5\г + (Я - Ь) т + ф*
Поскольку точки отрыва монотонно сближаются, при т > 0 функция у>о(ж1,у!) будет равна нулю в некоторой окрестности точки отрыва при у1 = Ь. Учитывая это, получим для функции ф* справа и слева от точки отрыва однородные краевые условия Дирихле и Неймана. В окрестности точки отрыва решение задачи представляется в полярных координатах в виде разложения по функциям вида сптп+о'5 сов(п + 0.5)у>, п = 0,1, 2,... Из условия регулярности решения задачи в точке отрыва следует, что коэффициент со равен нулю. Таким образом, искомое разложение имеет вид
ф*
с1т1'5 сов1.5у>, с1 = с1(т).
Примеры решения классической задачи об ударе с отрывом показывают, что С1(0) = 0. На основании этого будем считать, что С1(т) = 0 на некотором промежутке (0,т*). Для таких т согласно последнему соотношению получим в достаточно маленькой окрестности рассматриваемой точки соответственно слева и справа от нее следующие неравенства:
ддм* С1 < 0 ^ -1— > 0, ^ < 0; С1 > 0 ^ -1— < 0, ^ > 0.
ду1 ду1
Заметим, что случай С1 < 0 при т = 0 полностью согласуется с классической моделью удара с отрывом, а случай с1 > 0, напротив, ей противоречит. В конкретных примерах справедливость неравенства с1 < 0 для всех т (вплоть до момента схлопывания) можно проверить численно. Таким образом, при С1 < 0 условие Кутты - Жуковского приводит к неравенствам, которые по крайней мере локально совпадают с (24), (25). Отметим также , что при численной реализации задачи (18)-(23) достаточно локальной формулировки краевых условий типа неравенств.
Теперь остановимся на определении второго приближения. Функцию , удовлетворяющую уравнению Пуассона, можно свести к гармонической при помощи подстановки:
< \ 1 дф1 с1(т)
ф2 = в(т)х11пХ1---+ и в(т) ^ ---,
дх1 со(т )1п со(т)
где функция и, не зависящая от з(т), определяется решением следующей задачи:
Аи = 0, (26)
и = С(ф1,Х1,у1,т), у1 = Ь, 0 < Х1 < со(т); у1 = Я, 0 < Х1 < Як, (27) ди
— = 0, у1 = Ь, со(т) < Х1 < 1; у 1 = -Ь, 0 < Х1 < 1, (28)
ду1
ди ди д2ф1 тг . .
75— = 0, Х1 = 1, -Ь<у1 <Ь; — =т, у1 = -Яь, (29)
дХ1 ду1 дХ\
ди
-— =0, Х1 = 0, -Яь <у1 < -Ь, Ь<у1 < Я; Х1 = Як, (30)
дХ1
т т т
С(ф1,Х1,у1, т) = 0.5 ! §1 ¿в - дт ! (да + 1) ¿в - ^(т - в)да ¿в,
( )
9т = (xi,yi,т). dyi
После решения краевых задач для fi и u коэффициент ci(r) находится из условия регулярности функции f 2 в точке отрыва. При этом учитывается поведение функций и u вблизи точки раздела краевых условий:
У>1 ~ const Xi — с0(т )]3/2 , u ^ const [x1 — с0 (т )]1/2 при yi = b, xi ^ со(т) + 0 .В итоге приходим к следующей формуле:
( ) Д2(т) (31)
с1(т) =--ГТ > (31)
ai(T)
а1(т)= lim \Jxi — со(т), о,2(т) = lim \Jxi — со(т).
x 1 —>c0 (т)+0 dx2 xi —(т) + 0 dxi
Отметим, что рассуждения, проведенные при выводе формулы (31), аналогичны тем, которые были сделаны в [4] при построении асимптотики по малому времени.
Асимптотики (14)—(17) позволяют последовательно рассматривать вопросы, связанные с динамикой точек отрыва и определением форм внутренней и внешней свободных границ жидкости. Ниже приводятся окончательные формулы для коэффициентов асимптотических разложений (15), (16):
(32)
С1(жьт )= VI, С2(Ж1,Т ) = У1 = ь; фг(х! ,т ) = VI, ^2(ж1 ,т ) = V2, У1 = Н,
т т т
V1 = / (д8 + 1) ¿в, V2 = з(т)х11п Х1 I дх! ¿8 + | ди ¿80 0 0
Отметим, что для корректности рассуждений, связанных с регулярностью решения задачи и выводом формулы (31), вообще говоря, нужна непрерывная диф-ференцируемость функции О = 0(^1,Х1,у1,т) при у1 = Ь, Х1 € [0, со(т)]. Краевое условие (27), содержащее функцию О, было получено на основании представления функции С1(х1 ,т) в виде (32). Вследствие этого О является непрерывной в точке отрыва, а ее первая производная по переменной Х1 имеет в ней корневую особенность. Чтобы сгладить эту особенность, нужно правильно определить функцию С1(х1, т) в маленькой окрестности точки отрыва. Другими словами, нужно построить погранслойное решение в точке отрыва и согласовать его с внешним разложением (15). Вкратце остановимся на решении этого вопроса. Преобразованное кинематическое уравнение внутренней свободной границы имеет вид
+ £ = ~2 ^ + £-1 dyi dxi VdxiJ dxi
d( d( f dx \ 1 dx
дт dxi V dxi I дт
Погранслойное решение, согласующееся с внешним разложением (15) и удовлетворяющее начальному условию при т = 0, можно представить в виде:
С(Х1,т)= £7/2 4Г)№) - Пхо)], X = , = . (33)
со(т) е е
Рис. 2. Согласование погранслойных решений и внешнего разложения при т = 0.625:
е = Fro = 0.3, хо = -5.(5)
Сначала оно ищется в форме (33) с Е(Хо) = 0. На основании регулярности решения задачи с односторонними ограничениями левая часть дифференциального уравнения заменяется асимптотикой (в главном приближении ф ^ ф1):
дф
---+1 - в(т)л/со(т) - Х1, Х1 ^ со(т) - 0,
ду1
Производная функции ф1 по Х1 , которая входит в правую часть дифференциального уравнения, стремится к нулю при у1 ^ Ь. Поэтому, чтобы получить главный член ее разложения, нужно сделать один шаг по формуле Тейлора (разложение ведется по у1 при у1 = Ь и фиксированном Х1 ), а затем учесть асимптотическое поведение функции ф1 в окрестности точки раздела краевых условий:
дф1 д2ф1 А ( ) + в(т) м )
1Г~ = я—^ (Х1,т) +-----—===== С (Х1,т).
дХ1 дХ1ду1 2у/со(т) - ж
После перехода в дифференциальном уравнении к погранслойной переменной в главном приближении по е (порядок е3/2) найдем функцию Е(ж):
F'(х) = Vx, F(ж) = - ж'
= f. Х3/2 3
Построенное таким образом решение необходимо подправить, так как в силу отрицательности коэффициента (с0(т) < 0, ß(r) > 0) оно уходит внутрь прямоугольника. Непосредственно можно убедиться, что для производной ZT(x1,т) (по-прежнему считаем, что F(Хо) = 0) будет выполнено условие сращивания с внешним разложением в промежуточном пределе (е ^ 0, xп фиксировано) [9]:
со(т) - xi 0 п(е)
xv = -^-, п(е) ^ 0,--
П(е) е
Используя установленный факт, получим искомую погранслойную функцию с помощью интегрирования по времени. Так как зона сращивания сама зависит от времени, оправданным будет только локальное интегрирование в окрестности фиксированного т. При этом возникает произвольная постоянная (зависящая от данного т), которую можно представить в виде e7/2ß(т)/c0(т) • y. С целью сохранения структуры погранслойного решения (временной коэффициент умножается на функцию, зависящую от времени только через погранслойную переменную) будем считать величину y не зависящей от времени. Определяя ее однозначно из начального условия, придем к представлению (33). Полученная формула хорошо подтверждается численными экспериментами для различных значений параметров задачи. На рис. 2 показано согласование погранслойного решения (сплошная линия) с внешним разложением (пунктирная линия) при т = 0.625 (е = 0.3, хо = = -5.(5)). В этом случае величины [3(т) и с0(т) приближенно равны 1.24 и -0.6. При дальнейшем уменьшении параметра е кривые практически сливаются.
Решение кинематического уравнения внутренней свободной границы, удовлетворяющее однородному условию при xi = со(т), будем искать в виде:
Z Ы,т )= eaF (Ж) + ••• , Ж = Со(т x,
где постоянные а и А подбираются таким образом, чтобы первые два члена в правой части дифференциального уравнения (нелинейный член и производная по времени) имели одинаковый порядок по е, а остальные члены имели более высокий порядок малости. Одинаковые порядки по е получаются при А = 2а + 4. В этом случае функция Я(Х) имеет вид
я (Х) = - ем
Возвращаясь к переменной Х1, получим для возмущения свободной границы вблизи точки отрыва асимптотическую формулу
С(Х1,т) ~ -е -2 2вт) л/со(т) - Х1.
Р(т)
Отметим, что при е ^ 0 возмущение свободной границы должно исчезать. Таким образом, множитель е- 2 означает, что последнюю формулу нужно рассматривать в очень маленькой окрестности точки отрыва (корень есть маленькая величина). Вследствие этого кривая выходит из точки отрыва под прямым углом и при малых е практически не отличается от вертикального отрезка прямой (рис. 2). В промежуточном пределе (п(е)/ех ^ ж) главные по е члены двух по-гранслойных решений совпадут, если п(е) = е8, А > 8 (например, А =10, а = 3), а фиксированная величина хп = [со(0) — со(т)]3в4(т)/[2со(т)]4.
Заметим, что особенность у производной в граничном условии (27) создается за счет функции — (дт +1)^1 (х1, т), которая выделяется из ее среднего члена (к функции дт нужно прибавить и вычесть единицу). Будем считать, что за исключением очень маленькой окрестности точки отрыва функция £1(х1,т) определена по первой формуле (32), а в указанной окрестности имеет вид погранслойного решения.
Обратим также внимание на то, что погранслойное решение, удовлетворяющее условию сращивания, после перехода к старой переменной Х1, будет сохранять вид внешнего разложения (е2 умножается на функцию, не зависящую от е), а второе погранслойное решение, напротив, не будет его сохранять. В связи с этим вторые члены разложений (14)—(17) могут получить дополнительную зависимость от е. Это не нарушает рассуждений, проведенных при выводе вторых членов асимптотики, так как отброшенные члены имеют более высокие порядки малости. Однако можно предложить другой способ сглаживания, который не приводит к дополнительной зависимости вторых членов асимптотических разложений от е . В нем вместо второго погранслойного решения нужно использовать сглаживающую функцию, сохраняющую вид внешнего разложения. Эта функция будет связывать первое погранслойное решение с точкой отрыва. Теоретически оба подхода приводят к эффективному приближенному решению задачи, удовлетворяющему условию Кутты-Жуковского. Различие в этих подходах (например, при определении с1(т)) может проявиться только при очень малых е, для которых вторые члены разложений теряют свою актуальность.
Следует отметить, что проведенное сглаживание практически не влияет на численное определение коэффициента с1(т) по формуле (31). Это объясняется тем, что для численного анализа задачи (26)-(30) применяется слабая вариационная постановка (в интегральной форме), которая сама по себе сглаживает особенности. Таким образом, при численной реализации можно использовать функцию О в первоначальном виде.
Табл. 1
¿0 0.05 0.1 0.15 0.2 0.23 0.25
авшр 0.77 0.64 0.48 0.24 0.02 0
Нпеаг 0.78 0.68 0.53 0.34 0.15 0
Табл. 2
¿0 0.05 0.1 0.15 0.16
автр 0.72 0.49 0.1 0
Ипеаг 0.73 0.53 0.17 0
Табл. 3
¿0 0.03 0.05 0.07 0.08 0.09
авшр 0.71 0.58 0.38 0.23 0.01
Ипеаг 0.72 0.59 0.4 0.28 0.06
3. Численная реализация и анализ результатов
Сначала рассмотрим задачу со свободной границей (18)-(25). При фиксированном т приходим к определению точки отрыва, отделяющей на границе тела область безотрывного обтекания от зоны отрыва. С помощью специального итерационного метода решение данной нелинейной задачи (с неизвестной зоной отрыва) сводится к последовательному решению линейных краевых задач с фиксированными точками раздела краевых условий первого и второго рода. Для начала итерационного процесса нужно выбрать положение точки отрыва так, чтобы в некоторой ее окрестности нарушалось условие (24). Этого всегда можно добиться, уменьшая начальную зону отрыва, то есть приближая со(т) к нулю. Точка, соответствующая локальному отрицательному минимуму функции (24), выбирается за следующее приближение к точке отрыва. Дальше процесс повторяется. На каждом шаге итерационного метода зона отрицательных значений функции (24) уменьшается. Процесс заканчивается, когда эта зона исчезает. Подробно этот алгоритм описан в работах автора, посвященных начальному этапу движения твердых тел в жидкости с учетом явления кавитации (библиографию можно посмотреть в [4]). Аналогично организуется итерационный процесс, основанный на уточнении кинематического условия (25). В этом случае главное приближение к зоне отрыва должно быть достаточно большим (величина со(т) приближается к единице). В каждом методе вторые условия в виде неравенств проверяются после окончания алгоритма. Отметим также, что линейные задачи, возникающие на каждом шаге итерационного процесса, решаются численно с помощью метода конечных элементов [10]. После нахождения главного приближения и определения функции со(т) решается задача (26)-(30). При этом интегралы по времени, входящие в динамические краевые условия, вычисляются с помощью формулы правых прямоугольников. Так как для определения функции дт при каждом фиксированном т нужно решить задачу со свободной границей (18)-(25), для вычисления этих интегралов, а следовательно, и для решения всей задачи нужно применить метод шагов по времени.
Таким образом, чтобы воспользоваться асимптотической формулой (17), нужно на каждом шаге по времени определить величину со(т), решая задачу со свободной границей (18)-(23), затем найти функцию и и приближенно определить коэффициент с1 (т) по формуле (31).
Численные результаты, найденные для точки отрыва с(Ь) = с(ет) сравниваются с данными, полученными на основе линеаризованной модели, предложенной в [2, 3]. Результаты сравнения при х = —0.5, Н = 1, Нь = —3, Нд = 4 приведены
Сп(1)
С0(0.9375)
С0(0.625)
С0(0.3125) С0(0)
Рис. 3. Подъем внутренней свободной границы при схлопывании каверны: е = Ег0 = 0.4, Хо = -3.125 , т = 6.25^о
в табл. 1-3. Во второй строке каждой из этих таблиц представлены результаты, найденные по асимптотической формуле (17). Третья строка содержит данные, полученные на основе специальной линеаризованной модели. Численные результаты, приведенные в табл. 1, получены при е = 0.5, хо = -2, т = 4Ьо. Табл. 2 соответствует случаю е = 0.4, хо = -3.125, т = 6.25Ьо. В табл. 3 собраны результаты, отвечающие е = 0.3, хо = -5.(5), т = 11.(1)Ьо. Заметим также, что для сравнения численных результатов используется связь между безразмерным параметром ¿о, введенным в работе [3] (там он обозначается через Ь) и параметром т: т = е-2Ьо . Подчеркнем, что разница между Ь и ¿о состоит в том, что при введении безразмерного времени вместо V пишется Уо.
Хорошо видно, что при уменьшении параметра е погрешность, полученная в результате сравнения, уменьшается. Таким образом, линеаризованная модель, предложенная в [2, 3], дает хорошее согласование с асимптотической теорией при малых скоростях движения цилиндра.
Теперь остановимся на определении формы внутренней свободной границы жидкости. Сравнение численных результатов, полученных на основе главного члена асимптотики и линеаризованной модели [3], показало их хорошее согласование во всем диапазоне изменения переменной Х1 . Проведенное сравнение подтверждает справедливость предложенной в [3] модели схлопывания тонкой каверны с учетом подъема внутренней свободной границы жидкости. Основной вывод состоит в том, что при малых скоростях движения цилиндра (и не очень большом давлении в каверне) точки отрыва монотонно сближаются, а свободная граница немного поднимается. На рис. 3 изображена динамика рассматриваемого процесса при е = = Его = 0.4, хо = -3.125, т = 6.25Ьо (остальные параметры не меняются). Анализ проведен без учета погранслойного решения, которое сглаживает углы, но при этом не влияет на форму каверны в целом. Отметим, что при дальнейшем уменьшении параметра е возмущение свободной границы становится менее заметным. Следовательно, процесс схлопывания каверны при малых скоростях движения цилиндра в основном сводится к изучению динамики точек отрыва. Заметим также, что использование двух членов асимптотики позволяет несколько расширить диапазон применимости рассматриваемой модели.
Заключение
Проведен асимптотический анализ задачи о медленных движениях прямоугольного цилиндра в жидкости после отрывного удара с учетом искусственной кавитации. На основании первых двух членов асимптотики изучена динамика точек отрыва внутренней свободной границы жидкости и проведено ее согласование со специальной линеаризованной моделью. Ввиду малости второго члена асимптотики в формуле (15) (порядок 0(е4), е ^ 0), анализ возмущения внутренней свободной границы жидкости выполнен в главном приближении. Проведенный асимптотический анализ подтверждает основной результат работы [3], состоящий в том, что монотонное сближение точек отрыва сопровождается небольшим подъемом внутренней свободной границы жидкости. При малых скоростях движения
цилиндра возмущением внутренней свободной границы можно пренебречь и свести основной вопрос к изучению динамики точек отрыва. Проведенное исследование показывает, что при малых скоростях движения цилиндра можно получить полное решение проблемы образования и схлопывания присоединенной каверны. При этом реальные числа Фруда, для которых применим предложенный метод, образуют довольно широкий диапазон: Fro < 0.5. При других числах Фруда динамика отрыва после удара изучена только на малых временах [4]. Таким образом, перспективы исследования таких задач связаны главным образом с изучением дальнейшей динамики каверны при больших и умеренных числах Фруда.
Литература
1. Седов Л.И. Плоские задачи гидродинамики и аэродинамики. - М.: Наука, 1966. -448 с.
2. Норкин М.В. Динамика точек отрыва при ударе плавающего кругового цилиндра // Прикл. механ. и техн. физика. - 2019. - Т. 60, № 5. -С. 19-27. - doi: 10.15372/PMTF20190503.
3. Норкин М.В. Динамика точек отрыва при вертикальном ударе плавающего прямоугольного цилиндра // Вестн. ЮУрГУ. Серия «Матем. моделирование и программирование». - 2020. - Т. 13, № 2. - С. 108-120. - doi: 10.14529/mmp200209.
4. Норкин М.В. Движение прямоугольного цилиндра в жидкости после удара на малых временах с образованием каверны // Сиб. журн. индустр. матем. - 2020. - Т. 23, № 2. - С. 106-118. - doi: 10.33048/SIBJIM.2020.23.208.
5. Аганин А.А., Ильгамов М.А., Мустафин И.Н. Ударная кавитация жидкости в цилиндрической емкости // Учен. зап. Казан. ун-та. Сер. Физ.-матем. науки. - 2020. -Т. 162, кн. 1. - C. 27-37. - doi: 10.26907/2541-7746.2020.1.27-37.
6. Hilmervik K.B., Tyvand P.A. Incompressible impulsive wall impact of liquid cylinders // J. Eng. Math. - 2017. - V. 103, No 1. - P. 159-171. - doi: 10.1007/s10665-016-9866-6.
7. Hilmervik K.B., Tyvand P.A. Impact of narrow plates on broader liquid bodies // Appl. Ocean Res. - 2019. - V. 87. - P. 247-255. - doi: 10.1016/j.apor.2019.04.002.
8. Savchenko Yu., Savchenko G., Semenov Yu.A. Impulsive motion inside a cylindrical cavity //Mathematics.-2020.-V. 8, No 2. - Art. 192, P. 1-11. - doi:10.3390/math8020192.
9. Коул Дж. Методы возмущений в прикладной математике. - М.: Мир, 1972. - 274 с.
10. Жуков М.Ю., Ширяева Е.В. Использование пакета конечных элеметов FreeFem+—+ для задач гидродинамики, электрофореза и биологии. - Ростов н/Д: Изд-во ЮФУ, 2008. - 256 c.
Поступила в редакцию 30.09.2020
Норкин Михаил Викторович, доктор физико-математических наук, профессор Института математики, механики и компьютерных наук Южный федеральный университет
ул. Мильчакова, д. 8а, г. Ростов-на-Дону, 344090, Россия E-mail: norkinmi@mail.ru
ISSN 2541-7746 (Print)
ISSN 2500-2198 (Online)
UCHENYE ZAPISKI KAZANSKOGO UNIVERSITETA. SERIYA FIZIKO-MATEMATICHESKIE NAUKI
(Proceedings of Kazan University. Physics and Mathematics Series)
2020, vol. 162, no. 4, pp. 426-440
doi: 10.26907/2541-7746.2020.4.426-440
Asymptotics of Slow Motions of a Rectangular Cylinder in a Liquid after a Separation Impact
M.V. Norkin
Southern Federal University, Rostov-on-Don, 344090, Russia E-mail: norkinmi@mail.ru
Received September 30, 2020 Abstract
The process of collapse of a cavern, which is formed following the separation impact produced by a rectangular cylinder in an ideal, incompressible, and heavy liquid, was studied. Assuming that the speed of the cylinder is low, asymptotics of the main characteristics of the impact were constructed. The difficulties arising in this case are mainly related to the fact that the dynamics of the separation points is not known in advance and depends on a small parameter, i.e., on the dimensionless velocity of the cylinder. With the help of a special change of variables, the matter was reduced to the study of a problem with the dynamics of separation points corresponding to the leading approximation independent of the indicated parameter. This enabled us to determine the second term of the asymptotics with account for the nonlinear terms in the model. In the leading approximation, a problem with a free boundary, which is similar to the classical model of an impact with separation at each fixed moment of time, was formulated. On the basis of the first two terms of the asymptotics, the process of collapse of the cavern was described with allowance for the rise of the inner free boundary of the liquid. Comparison with the previously obtained results was carried out.
Keywords: ideal incompressible fluid, rectangular cylinder, separation impact, dynamics of separation points, Froude number, cavitation number
Figure Captions
Fig. 1. Formulation of the problem.
Fig. 2. Matching the boundary layer solutions and the external expansion for e = Fro = 0.3, Xo = -5.(5).
Fig. 3. Lifting of the internal free boundary when the cavern collapses: e = 0.4, Xo = -3.125, t = 6.25to.
References
1. Sedov L.I. Ploskie zadachi gidrodinamiki i aerodinamiki [Two-Dimensional Problems of Hydrodynamics and Aerodynamics]. Moscow, Nauka, 1966. 448 p. (In Russian)
2. Norkin M.V. Dynamics of separation points upon impact of a floating circular cylinder. J. Appl. Mech. Tech. Phys., 2019, vol. 60, no. 5, pp. 798-804. doi: 10.1134/S0021894419050031.
т = 0.625: Fro = 0.4,
3. Norkin M.V. Dynamics of separation points during vertical impact of a floating rectangular cylinder. Vestn. YuUrGU. Ser. "Mat. Model. Program,.", 2020, vol. 13, no. 2, pp. 108-120. doi: 10.14529/mmp200209. (In Russian)
4. Norkin M.V. The movement of a rectangular cylinder in a liquid at short times after impact with formation of a cavity. J. Appl. Ind. Math.., 2020, vol. 14, no. 2, pp. 385-395. doi: 10.1134/S1990478920020155.
5. Aganin A.A., Ilgamov M.A., Mustafin I.N. Impact-induced cavitation in a cylindrical container with liquid. Uchenye Zapiski Kazanskogo Universiteta. Seriya Fiziko-Matema-ticheskie Nauki, 2020, vol. 162, no. 1, pp. 27-37. doi: 10.26907/2541-7746.2020.1.27-37. (In Russian)
6. Hilmervik K.B., Tyvand P.A. Incompressible impulsive wall impact of liquid cylinders. J. Eng. Math., 2017, vol. 103, no. 1, pp. 159-171. doi: 10.1007/s10665-016-9866-6.
7. Hilmervik K.B., Tyvand P.A. Impact of narrow plates on broader liquid bodies. Appl. Ocean Res., 2019, vol. 87, pp. 247-255. doi: 10.1016/j.apor.2019.04.002.
8. Savchenko Yu., Savchenko G., Semenov Yu.A. Impulsive motion inside a cylindrical cavity. Mathematics, 2020, vol. 8, no. 2, art. 192, pp. 1-11. doi: 10.3390/math8020192.
9. Cole D.J. Metody vozmushchenii v prikladnoi matemtike [Perturbation Methods in Applied Mathematics]. Moscow, Mir, 1972. 274 p. (In Russian)
10. Zhukov M.Yu., Shiryaeva E.V. Ispolzovanie paketa konechnykh elementov FreeFem+ + dlya zadach gidrodinamiki, elektroforeza i biologii [Using the FreeFem++ Finite Element Software for the Problems of Hydrodynamics, Electrophoresis, and Biology]. Rostov-on-Don, Izd. YuGU, 2008. 256 p. (In Russian)
Для цитирования: Норкин М.В. Асимптотика медленных движений прямоуголь-/ ного цилиндра в жидкости после отрывного удара // Учен. зап. Казан. ун-та. \ Сер. Физ.-матем. науки. - 2020. - Т. 162, кн. 4. - С. 426-440. - doi: 10.26907/25417746.2020.4.426-440.
For citation: Norkin M.V. Asymptotics of slow motions of a rectangular cylinder in / a liquid after a separation impact. Uchenye Zapiski Kazanskogo Universiteta. Seriya \ Fiziko-Matematicheskie Nauki, 2020, vol. 162, no. 4, pp. 426-440. doi: 10.26907/25417746.2020.4.426-440. (In Russian)