Движение пузырька газа в жидкости с учетом искажения его сферической формы Текст научной статьи по специальности «Физика»

УДК 532.51:534.12

ДВИЖЕНИЕ ПУЗЫРЬКА ГАЗА В ЖИДКОСТИ С УЧЕТОМ ИСКАЖЕНИЯ ЕГО СФЕРИЧЕСКОЙ ФОРМЫ

© М.А.Ильгамов, Л.А.Косолапова, В.Г.Малахов

В статье построена математическая модель динамики пузырька газа в жидкости с учетом искажения его сферической формы и перемещения в пространстве. Предполагается, что искажения формы пузырька и скорость его пространственного движения малы. Соотношения получены с точностью до членов второго порядка относительно малых величин. Приближенно учтены вязкость и сжимаемость жидкости. Приведен ряд примеров, иллюстрирующих влияние свойств жидкости на динамику воздушного пузырька в воде.

Ключевые слова: газовый пузырек, нелинейные колебания, потенциальное течение, численное моделирование.

1. Введение

Рассматривается нелинейная динамика газового пузырька в жидкости с учетом искажений формы его поверхности и пространственного перемещения. Как следует из обзоров, посвященных этой проблеме [1; 2], в большинстве работ жидкость считается идеальной и несжимаемой, а искажения формы и скорость движения пузырька малыми. Так, в работе А.А.Бошкоу [1] получены уравнения квадратичной модели, а в исследовании SJ.Show [2] - уравнения модели третьего порядка относительно малых величин. При этом в последней работе жидкость считается идеальной и несжимаемой, а в построениях А.А.Бо-т1коу вязкость учтена только в уравнениях для сферических колебаний и поступательного движения пузырька. Вязкость при поступательном движении в исследовании А.А.Бо1шкоу учтена согласно вычислениям Б.У.ЬеуюЬ.

В настоящей работе, являющейся развитием выводов, изложенных в статье А.А.Аганина, М.А.Ильгамова, Л. А.Косолаповой, В.Г.Малахова [4], построена квадратичная относительно искажений сферической формы пузырька и скорости его пространственного перемещения математическая модель динамики пузырька газа в жидкости с учетом движения жидкости и приближенном учете ее вязкости и сжимаемости.

2. Постановка задачи

Изучается осесимметричная динамика газового пузырька в неограниченном объеме жидкости. Движение жидкости предполагается потенциальным. Вязкость учитывается через динамическое граничное условие на поверхности пузырька. Жидкость в ближней к пузырьку области г < ге (ближнее поле) считается несжимаемой. Тогда движение жидкости описывается уравнениями

АФ = o, Фt +1 ^Ф)2 +P = G (t), 2 Po

(1)

где t - время, p - давление в жидкости, po -плотность жидкости, Ф - потенциал скорости, G (t) - произвольная функция времени. Наряду с неподвижной цилиндрической системой координат х, z, ф вводится подвижная сферическая система координат r, Є, ф с началом в центре пузырька, связанная в осесимметричном случае с неподвижной системой уравнениями х = r sin Є, z = zc (t) + r cos Є, где zc - координата центра пузырька в неподвижной системе отсчета в момент времени t .

В подвижной системе координат второе уравнение (1) запишется в виде

, . - sin Є, .

Ф - ^^c (Фr cos е-------Фе ) +

r

(2)

+-(ф 2+

J_ ф2) + JL = G (t), r Po

где - скорость движения начала координат подвижной системы.

Считается, что газ в пузырьке идеальный, с однородным давлением, изменяющимся по закону

,Г V - ьу0 У

- o P = Pb

_o_

V - bVo

o

(3)

р0 - начальное давление газа в пузырьке, , У

- начальный и текущий объемы пузырька, Ь -постоянная, у - показатель адиабаты.

Предполагается, что на бесконечном удалении от пузырька жидкость движется со скоростью , определяемой выражением

ию = U (t) = Ui (t)ez

(4)

где ez - направляющий вектор оси z, а давление жидкости рда на уровне z = o изменяется по закону

Рда= Po - Ар sin rat, (5)

где Po - статическое давление в жидкости; Ар, га - амплитуда и частота колебаний.

На поверхности пузырька r = rs (е, t) в подвижной сферической системе координат ставится кинематическое условие в виде

Ft +(Фr - ^^c cosе)) +

(б)

Г Фе

r

zc sin е IГ Fе = o

рмоникам

да

1 + XSi(t) Pi (n)

i=2

(S)

Ф(е, r, t) =

Bo

B,

ulr-

B1

2

P

(9)

^-7+Г P (cos е).

i=2'

и динамическое граничное условие в виде

Р + = Р"- 2сН + тпп , (7)

где с - поверхностное натяжение, Н - средняя кривизна поверхности, тпп - вязкие напряжения,

2 2

тпп = 2цд Ф / дп , ц - динамический коэффициент вязкости. Условие (7) представляет собой равенства нулю суммы проекций сил на нормаль п к деформированной поверхности пузырька. Предположение о потенциальности течения жидкости исключает возможность использования динамического граничного условия по касательной, поскольку система становится переопределенной. Поэтому в настоящей работе требуется выполнение только условия (7). Сравнение моделей с различными способами учета вязкости показывает возможность такого подхода [5].

3. Метод решения Уравнение поверхности пузырька в сферической системе координат запишем в виде ряда по сферическим га

^ (0, {) = я Ц)

Дальнейшие соотношения получены в предположении малости искажений сферической формы 8; , скорости движения пузырька 2С и несферических составляющих потенциала скорости жидкости. При этом в уравнениях сохранялись члены второго порядка относительно малых величин.

Строгий учет сжимаемости является довольно сложной задачей. Поэтому в настоящей работе используется методика, которая предполагает, что несферические возмущения пузырька оказывают заметное влияние на движение жидкости лишь в некоторой области г < ге (ближнее поле), в которой жидкость считается несжимаемой [6]. При г ^ ге (дальнее поле) сжимаемость учитывается в рамках одномерной линейной акустики. Из условия сопряжения решений на границе

ближнего и дальнего полей r = re ражение для функции G (t) в виде

G(t) _-P^ + — ul - zcui +---------------------Bo

Po 2 cf

1

получаем вы-

(io)

Здесь P - полином Лежандра степени i, П = cos е . Параметр R будем называть радиусом пузырька, S, - искажением сферической формы по гармонике с номером i .

Потенциал скорости Ф берется в виде

При учете сжимаемости принималось дополнительное предположение

|Я/С7| ~ 8 . (11)

Подставим разложения (8), (9) в условия (6), (7) с учетом (2), (10). Применяя при интегрировании полученной системы по координате 0 процедуру метода Бубнова-Галеркина, приходим к системе дифференциальных уравнений первого порядка относительно неизвестных

в ________

я, 8;, Во, 5Ь 5;, 2с, ; = 2,3,... где 5; = —1-,; = 1,да .

я1

Исключив из нее неизвестные Во, 51,5;, запишем систему обыкновенных дифференциальных уравнений второго порядка относительно радиуса Я(?), искажений 8; (/) и координаты центра пузырька хС (/) в виде

1+^ S2 -±

_ 2 cf

R+а iiRSiS i- (1+ 3 a ii s2 )

PoR 2

RR+^i-

2

Poc/R

R2SiS'i ——ww =

2c

f

3

2

1 -■

R

\

3c

f

R 2 - a a (R2 s2 + 4 (Rs, S, + R 2S

(2i +1) 4(i +1)'

a.

■(9 R 2Si

(12)

2 1 б RiRsiS, + R 2s ,2 )

О | Г'!

Ро

и Я '

С{

Л ^ У

( -Рю)+)(р -рда)-) 2+(1-;2)<а;;82

Г л I

Ро я

[4Я - ай (б/б,-

-(/ + 1)( - 2)Я8; ){

С/

1+

я

С/

4

ар1-3Ур1 IяЯ8*■8j + (^у-1-Уу1 )я2'8;8]- 31 1 -582 Iя™ = --яI 1 + 582 Iи1+ ^Ч2] (2 р1

-а -

3Р р1

2(; +1)( ]+1)

■уг

Р*/1

■2^/1

+-

%1

О +1)( р +1)

( 2!в//1 -( + 1)р' + 3)ар1 )(3Я8г- + Я8г- )8)

+ -

Ро я

(; +1)(/ +1) ■/ * / [(58г8(+ 38;8/)ур1-] (+

V У

( - / - / 2 ) ау!8;8/- 3ц [(3(<Ху1- (/1)Я8г- 2Уу1Я8г- )8)- ( р1- (Зу1) )

(13)

к - 2_ Г 2к + 5 _ .

к+1СХкк8к +[“к+Г<а/- 3у/к ]8;8/ а кк

ЯЯ - Я

сх кк

к:^8:к- (а / -у/к)8 л

= -(Х ;к11Л1Я81 + —^ Я (3^^8к +5Я8 к ) + 3^^2 к +1

к + 3 _

УУк +па/к -

Я

3Ррк

сх 1кГ

35

1'к

Л

2(к +1)

8- =

2(к +1)

2(/ +1)( ] +1)

Я2

8i8 .■------------X

1 ] 2

сх

рк

в рк

О + 1)( 7 + 1)

■2^)к

8г■8 / + ЯЯ?

(388 р + 58г-8 / ))

5к +11

к +1

'а рк+

3Ррк

О' +1)( р +1)

8;8 р

1

Ро Я

(к - 1)(к + 2)акк8к + 2а/к Р1 -;-;2)8г-8р -2ц[-<акк (3кЯ8к + (к + 2)Я8к) +

(14)

+3(Iерк- 4«рк )Я88р - (2Ррк-Р + 1) Р + 3)<%к ) (Р + )87

Кк+ 2)3Я((к - (к +1)Урк )8г8/ - (к + 1)ЯijkЯ8i2

(Р 8+ Я8г)

Рьк^ (3^Я8*■ + Я8*■)- [3 (Ряк - 4«гк1)

2 р +1)

Ро Я

1 3 I 3 (к + 2) 5

-1 + 582 I^8г------^5Ш^’81

к = 2,3,. .

В (12)-(14) обозначено

<а;к=| Р (п)Рк (п)^

-1

/ = к

1

1

арк = | рРрРкЛя Ррк = |рРРк Р1 -п2)dn, -1 -1

dPi (п) ~ ~

Р =--------, 5рк = 2арк - Ррк,

Урк

d п

(; + 1) Р + 2) <%к - Iе рк

р +1) (к + 1)

Выражение для текущего объема пузырька с принятой точностью имеет вид

К = 3пЯ3 ^1 + |82а,.|. (15)

В (12)-(15) для сокращения записи принято соглашение о суммировании: если слагаемое содержит два одинаковых индекса, то по ним предполагается суммирование.

Отметим, что соотношения работы [1] получаются как частный случай (12)-(14) при исключении квадратичных слагаемых вида с8г■8р и

слагаемых, учитывающих сжимаемость и вязкость жидкости, кроме членов 4цЯ / Ро Я и 4ц№ / Ро Я в уравнениях (12), (13).

Система (12)-(14) с начальными условиями

2

Я(0) = % Я(0) = 0, 8г (0) = 8г°, єг- (0) = 0, г(0) = zc°, w (0) = ^0

решается методом Дормана-Принса седьмого порядка точности [7]. При вычислениях в представлении (8), (9) бесконечные верхние пределы сумм заменяются на число М, тогда система уравнений содержит 2М+2 неизвестные функции.

4. Результаты решения задач

4.1. Для проверки полученных соотношений была решена задача о движении несферического пузырька воздуха в воде при постоянном "на бесконечности" давлении в жидкости. В качестве возмущения в момент времени і = 0 задается ма-

лое искажение формы по третьей гармонике. Давление внутри пузырька изменяется по закону

(3). Результаты получены при следующих исходных данных:

Яо = 5омкм, у = 1,4, Ь = о, Су = 15оом/с,

Ру = 998кг/м3, с = о,о725кг/с2, ц = 1о-3кг/мс,

рда =о,1МПа, ^о=о, 8° = о,1, М = 4.

Решение данной задачи приведено в [1]. На рис.1 тонкими сплошными линиями изображено решение, полученное с использованием соотношений [1].

Я-Я0 , мкм 0

“I '----1----'--1---1 Г

0 2 4 а 6 8 іш/2л

103 8

2

4Н

0

-4Н

г , мкм с’ _____

-0.03

1---------'--------1---------1-------Т"

0 10 20 30 40 і ю/2л

Ь

Рис.1. Зависимости от времени величин Я - Яо (а), 2С (Ь), 82 (с), 83 (й), 84 (е). Пунктирные линии -1/су = о, ц = о, толстые линии - 1/су ф о, ц ф о (настоящая работа), тонкие линии - по уравнениям [1]

0

Толстыми сплошными и пунктирными линиями приведены результаты настоящей работы: пунктирные - решение для идеальной несжимаемой жидкости, толстые линии - с учетом свойств жидкости. Видно, что характер решений качественно меняется в зависимости от используемых соотношений. Поскольку в данной задаче возмущение радиуса Я - Яо и искажения сферической формы пузырька 8; с / Ф 3 малы (они воз-

никают в результате проявления нелинейного взаимодействия), разница в решении обусловлена, в основном, способом приближенного учета вязкости жидкости. При использовании соотношений настоящей работы с течением времени скорость пузырька обращается в нуль, а амплитуды всех гармоник затухают. Отметим, что решения [1] и настоящей работы для идеальной несжимаемой жидкости (тонкие сплошные и пунктирные линии) также качественно различны. Та-

ким образом, учет в силах поверхностного натяжения слагаемых второго порядка малости оказывает существенное влияние на решение данной задачи.

4.2. Движение в воде пустой несферической полости с заданной начальной скоростью. Давление в жидкости "на бесконечности" постоянное. Брались следующие исходные данные:

Яо = 5омкм, Су = 15оом/с, ру = 998кг/м3, с = о,о725кг/с2, ц = 1о-3кг/мс,рда =о,1МПа, рЬ = о, ^о=о,1м/с, 82 = о .

Искажение сферической формы пузырька допускалось в виде второй сферической гармоники. В начальный момент времени в сферической полости при заданной начальной скорости в процессе движения возникает искажение 82, обусловленное нелинейным взаимодействием между скоростью пространственного перемещения пузырька и искажением по второй сферической гармонике. На рис.2. приведены зависимости от времени отношения Я/Яо, искажения 82, координаты центра пузырька хС и скорости ее изменения для 4 вариантов (номер кривых соответствует номеру варианта): 1 - идеальная несжимае-

мая жидкость без учета поверхностного натяжения, 2 - с учетом поверхностного натяжения, 3 -с учетом вязкости и сжимаемости и 4 - с учетом поверхностного натяжения, вязкости и сжимаемости жидкости. Во всех вариантах на достаточно коротком интервале времени наблюдается резкий рост искажения 82. Следует отметить, что приведенные результаты можно получить, используя упрощенные уравнения, когда в соотношениях пренебрегается слагаемыми порядка

0(8 ), за исключением нескольких членов вида ^82, ^82, которые обеспечивают взаимодействие между искажением формы и пространственным перемещением. Видно, что учет свойств жидкости влияет на время разрушения пузырька ^ ( |8^ ~ 1 ). В основном, время { сокращается за

счет влияния поверхностного натяжения (линия

2 по сравнению с линией 1), а учет вязкости и сжимаемости приводит к его увеличению (линия

3 по сравнению с линией 2). На рис.2Ь штриховой линией приведена форма полости в момент I = о и сплошной линией для варианта 4 в момент ? = 4,527 мкс , когда наблюдается заметное искажение сферической формы.

2с, МКМ

Рис.2. Зависимости от времени отношения Я /Яо , искажения 82 , координаты центра пузырька 2С и его скорости в случае: 1 - с = о, ц = о, 1/су = о; 2 - сФо, ц = о, 1/су = о; 3 - с = о, цФо, 1/су Фо ; 4 -

сФ о, цФ о, 1/су Ф о

4.3. Движение несферического пузырька воз- Яо= 2омкм, у = 1,4, Ь = 8,5, Су = 15оом/с,

духа в воде под действием периодически изме- 3 2 3

няющегося давления жидкости (5). Для расчета ру = 998кг/м , с = о,о725кг/с , ц =1о

взяты данные: ро=о,1МПа, Ар = о,обМПа, и /2п = 26,5кГц,

^о=о, 8о = о,3, 8го = о,; Ф 5; М = 1о.

ію/27

0-

-0^

ію/2п

Рис.3. Зависимости от времени отношения Я/Яо, координаты центра пузырька гС и искажений 83, 85 : а, Ь, с, й - идеальная несжимаемая жидкость; е, у g, к - с учетом свойств

жидкости

На рис.3 приведены зависимости от времени отношения Я/Яо, координаты центра пузырька гС и искажений 83, 85 для идеальной несжимаемой жидкости (линии а, Ь, с, й) и с учетом вязкости и сжимаемости (линии е,у, g, к). Видно, что в идеальной несжимаемой жидкости при (> 1, 5Т, где Т - период изменения давления жидкости, наблюдается резкий рост искажений и перемещения пузырька. Расчеты показывают, что учет сжимаемости мало влияет на решение задачи. И лишь учет вязкости приводит к качественным изменениям результатов. При этом радиус пузырька Я совершает периодические колебания с периодом Т. Заданное в начальный момент времени искажение 85 в течение двух периодов убывает до нуля (рис.3к), но за это время в результате нелинейного взаимодействия возникают другие гармоники и движение пузырька вдоль оси г. За период времени ~ 5Т устанавливаются периодические колебания всех гармоник. Наибольшая амплитуда наблюдается у третьей гармоники (рис^). Амплитуды других гармо-

ник значительно меньше. Одновременно с развитием искажений начинается и поступательное движение пузырька, которое затем переходит в его колебания относительно оси г.

5. Заключение Построена математическая модель нелинейной динамики несферического газового пузырька в жидкости с учетом искажения его сферической формы и пространственного перемещения. Движение жидкости считается потенциальным. При этом вязкость жидкости учитывается через динамическое граничное условие на поверхности пузырька. Сжимаемость жидкости учитывается лишь в дальней от пузырька области. Предполагается, что искажения сферической формы пузырька и скорость его движения малы. Соотношения получены с учетом слагаемых второго порядка относительно малых величин. Проведено сравнение полученных результатов с результатами работы [1]. Наблюдается их удовлетворительное совпадение. Решены задачи о движении пустой несферической полости в воде при постоянном "на бесконечности" давлении жидкости и

о колебаниях пузырька воздуха в воде под воздействием периодически изменяющегося давления "на бесконечности". Показано, что при пренебрежении сжимаемостью жидкости, ее вязкостью или поверхностным натяжением могут получиться качественно различные решения.

Работа выполнена в рамках программы РАН и при поддержке РФФИ (код проекта 08-01-00215).

1. Doinikov A.A. Translational motion of a bubble undergoing shape oscillations // J. Fluid. Mech. - 2004.

- V.501. - Р.1-24.

2. Show S.J. Translation and oscillation of a bubble under axisymmetric deformation // Physics of fluids. -2006. - 18. - 072104. - 15 p.

3. Levich B.V. Physicochemical Hydrodynamics. -Prentice-Hall, 1962.

4. Аганин А.А., Ильгамов М.А., Косолапова Л.А., Малахов В.Г. Нелинейные несферические колебания пузырька газа при периодическом изменении давления окружающей жидкости // Теплофизика и аэромеханика. - 2оо8. - Т.15. - №3. - С.521-533.

5. Нигматуллин Р.И., Ахатов И.Ш., Вахитова Н.К. О сжимаемости жидкости в динамике газового пузырька // Докл. РАН. - 1996. - Т.348. - №6. -С.768-771.

6. Ильгамов М.А., Аганин А.А., Косолапова Л.А. и др. Модели динамики несферического пузырька с учетом вязкости жидкости // Тр. мат. центра им. Н.И.Лобачевского. Модели механики сплошной среды. Материалы ХУ1-ой сессии Междун. школы по моделям механики сплошной среды. - Казань: Изд-во Казан. мат. об-ва, 2оо2. - Т. 16. -С.192-2Ы.

7. Хайрер Э., Нерсетт С., Ваннер Г. Решение обыкновенных дифференциальных уравнений. Нежесткие задачи. - М.: Мир, 199о. - 512 с.

TRANSLATIONAL MOTION OF A GAS BUBBLE IN A LIQUID WITH ALLOWING FOR DISTORTIONS OF ITS SPHERICAL SHAPE

M.A.Ilgamov, L.A.Kosolapova, V.G.Malakhov

A mathematical model of dynamics of a gas bubble in a liquid with allowing for distortion of its spherical shape and its translational movement has been developed. The distortions of the shape of the bubble and the velocity of its spatial movement are assumed to be small. The resolving equations have been obtained up to the terms of the second order of accuracy in the small values. Viscosity and compressibility of the liquid are taken into account approximately. A number of examples are presented, showing the influence of the liquid properties on the dynamics of a gas bubble in water.

Key words: gas bubble, nonlinear oscillation, potential flow, numerical simulation.

Ильгамов Марат Аксанович - доктор физико-математических наук, член-корреспондент РАН, заведующий лабораторией Института механики Уфимского научного центра РАН, главный научный специалист Института механики и машиностроения Казанского научного центра РАН.

E-mail: [email protected]

Косолапова Людмила Александровна - кандидат физико-математических наук, старший научный сотрудник Института механики и машиностроения Казанского научного центра РАН.

E-mail: [email protected]

Малахов Владимир Георгиевич - кандидат физико-математических наук, старший научный сотрудник Института механики и машиностроения Казанского научного центра РАН.

E-mail: [email protected]