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

ВЕСТНИК ТГГПУ. 2011. №3(25)

ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ

УДК 532.5:534-14

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

© А.А.Аганин, А.И.Давлетшин

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

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

1. Введение

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

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

По мере сближения пузырьков взаимодействие между ними усиливается. В качестве меры близости двух пузырьков принято считать безразмерный параметр 5 = (Я1 + Я2) / й < 1, где Я1, Я2 - радиусы пузырьков, й - расстояние между их центрами. К настоящему времени наиболее исследованным является дальнее взаимодействие, когда расстояние между пузырьками относительно велико (когда 5 < 0,1). Изучение дальнего взаимодействия началось еще на рубеже 19 и 20 веков [2]. Для этого применялись в основном математические модели второго порядка точности относительно 5. Однако для адекватного описания наблюдаемых экспериментально кластеров и стримеров такой точности недостаточно, поскольку расстояния между пузырьками в этих структурах могут быть меньше тех, что характерны для дальнего взаимодействия. Для изучения более близкого взаимодействия были

разработаны модели, имеющие относительно 5 третий и четвертый порядки точности [1; 3; 4; 5].

В настоящей работе предлагаются

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

пузырьков и их пространственных положений и первого порядка для температур газа в каждом из пузырьков. В уравнения для радиусов пузырьков и их пространственных положений входят коэффициенты представления потенциала скорости жидкости в виде суммы частных

решений уравнения Лапласа в локальных сферических координатах пузырьков, которые определяются из системы линейных

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

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

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

— +1 VO .VO+ Р Р” = 0.

д t 2 рп

Здесь р0 - плотность жидкости, р - давление в произвольной точке жидкости, рда - давление жидкости на бесконечном удалении от пузырьков.

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

dt

+ VO.VFk = 0 .

отсчета.

Рассматривается взаимодействие двух газовых пузырьков при гармоническом изменении давления жидкости. Сила тяжести не учитывается. Пузырьки считаются сферическими, а газ в них - невязким гомобарическим. Такие допущения можно использовать при достаточно малых размерах пузырьков. В частности, можно ожидать, что форма взаимодействующих пузырьков будет близкой к сферической, если их равновесные радиусы меньше критического Я* относительно устойчивости сферической формы отдельного пузырька при его радиальных колебаниях в рассматриваемых условиях. Например, для воздушного пузырька в воде при комнатных условиях, когда давление жидкости изменяется гармонически с амплитудой 0,8 бар и частотой 20,6 кГц, имеем Я* « 16 мкм [7].

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

В неподвижной декартовой системе координат х, у, z, осью z которой является ось симметрии задачи (рис.1), уравнения поверхностей пузырьков можно записать в виде ^к(х, у, z, ¿) = 0, где t - время. Здесь и далее параметры, относящиеся к первому пузырьку (к = 1), будем помечать нижним индексом 1, а параметры, относящиеся ко второму (к = 2) - индексом 2.

Без учета эффектов вязкости и сжимаемости жидкости, теплообмена между пузырьками и жидкостью уравнения динамики жидкости в терминах потенциала скорости Ф выглядят следующим образом

V2 Ф = 0, (1)

а динамическое - следующим 2а

Р = Рк - Т •

(3)

(4)

Здесь Як - радиус к-го пузырька, а - коэффициент поверхностного натяжения, рк - давление газа в к-м пузырьке, изменяющееся по закону

Рк =

p0 +

2а

R

0, к

R0J V Rk J

(5)

где р0 - статическое давление жидкости, R0, к -начальный радиус пузырька. Для изотермического пузырька Гк = 1, а для адиабатического Гк = к, к - показатель адиабаты.

На бесконечном удалении от пузырьков давление жидкости рда совершает гармонические колебания с амплитудой Ар и частотой ю по закону р»= Ро -АРsin at • (6)

3. Уравнения для радиусов пузырьков и координат их центров При выводе уравнений взаимодействия пузырьков наряду с неподвижной декартовой системой координат х, y, z используются также и подвижные сферические системы отсчета r1, 01, ф1 и r2, 02, ф2 с радиальными координатами, отсчитываемыми от центров пузырьков (рис.1). Расстояние между центрами пузырьков d определяется выражением d = z1 - z2, где z1, z2 - координаты центров (z1, > z2). Неподвижная и подвижные системы отсчета связаны соотношениями х = х(гк, 0к), У = y(rk, 0к), z = z(rk, 0к, t), так что для произвольной функции f=f (х, y, z, t) имеем f = Df д t Dt

Vf,

(7)

где д // д ґ, Б // Б ґ - производные по времени при постоянных х, у, z и гк, 0к, фк, соответственно, wk = ^к к - скорость подвижной системы координат, к - единичный вектор оси z. Имеем wk = Zk, где Zk - скорость пространственного перемещения центра к-го пузырька. Здесь и далее верхней точкой обозначается производная по времени.

С учетом (7) интеграл Коши-Лагранжа (2) и кинематическое условие (3) примут следующий вид

- »УФ +1 УФ»УФ +

Р - Рда

=0, (8)

Щ.

Бґ

- Уф.У^к = 0 .

(9)

Уравнение поверхности пузырька в подвижной сферической системе координат записывается в виде

рк (гк, т) = гк - Як (т) = 0. (10)

Потенциал скорости жидкости при наличии в ней двух пузырьков принимается следующим [5]

, В..,, Р,

Ф =

да В Р

іу1^ У1

¿—І г У+1

7=0 А1

у=0

у 2 у 2 гУ+1

(11)

Здесь Ву к = Ву к (ґ), Ру к = Ру(оо8 0к) - полином Лежандра степени у от ео8 0к. Первое слагаемое в правой части (11), обусловленное наличием 1-го пузырька, записано в системе координат 1-го пузырька, а второе, обусловленное наличием 2-го пузырька, - в системе координат 2-го.

Для представления потенциала (11) в системе координат к-го пузырька воспользуемся соотношением [8]

Р

у 3-к

у+1

= 1

С гд Р

к1 д к у+д ^5+^+1

(12)

где Су ? = (-1)? (у + д)! / (у!д!), 5: 2 = 1, ^2 1 = -1,

Гк < й.

С учетом (12) потенциал (11) в окрестности к-го пузырька (при гк < й) запишется следующим образом

да В Р

ф = ^ _ ,к Рк

„У+1

■I

у,д=0

С д ВУ 3-к гк Рд к *1+-^’

(13)

В0к = ЯкЯк ,

В = -Як В1к = 2

да СВ ^

I у1 У3-к

В = пЯк

Впк =

V

2 п+1

г I

А V=0

1 с^+1 иу+2

У=0 ^ к 3-ки

СВ 3 к

уп у3-к

(14)

сУ+п ЛУ+п+1 7=0 дк3-кМ

, где

к = 1,2, п = 2, 3, ...

Система уравнений для радиусов пузырьков и координат их центров находится с использованием динамических условий (4), которые с учетом (8) становятся следующими

Б Ф Бґ

- -щ,, .УФ+1 УФ.УФ

Л

= Р да

2с

Рк - Я

(15)

к у

Подставив в (15) выражение потенциала (13) и проведя ряд преобразований с учетом ортогональности полиномов Лежандра, получим следующие обыкновенные дифференциальные уравнения

-I

у=0

В0к = Як "

(у + 1)В.

Рк - Рда

2а р0 Як

2Я

2 у+4

Су 0 Ву 3-к

*кУ3-к И ^

СВ Т

у3-И 3-к *!+-к^ ^2

(16)

-I

дСС в^в^Я?-2

2*У+? ИУ+^+2д+2 у,с,д=0 ^к3-к а

І=я z -'т?£!Мі

Г)2 к к л пУ+д+4

як у, д=0 2Як

-I

у=0

(

С1 В! 3-к Як

* У+1 И ^+2

к 3-к и

2СУ2 В3-кЯкТ3-к

+ 2;

У, С, 5=0

©Г ССд ВУкВС3-к^г"3

(17)

- I

у, д, С, 9=0

©д® С СВ В я д+9-2

д9 уд ^9 у3-к С 3-к ІХ-к

2 *у+д+С+9 лу+С+д+9+2 *

^к3-к а

Здесь числа 0^ ются выражениями

®пт = пт ауд! +

<н^пт

уд 5

«у с д, Ру с д определя-

ем = пт а„д1 -р

у=0 'к

Подставляя выражения поверхности пузырьков ^к (гк, т) из (10) и выражения потенциала Ф из (13) в кинематические условия (9) и учитывая ортогональность полиномов Лежандра, получим следующую систему линейных алгебраических уравнений для определения коэффициентов в представлении потенциала в виде (11)

а » =

у^д

+ Руд1 ..

^ 1РР Р-Рп

'уд1

уд1 :

Р,«= ^}(1 -П2)Р,РР*1 ,

£ -1

где п - аргумент полиномов Лежандра.

Дифференцируя первые два уравнения (14), подставляя в полученные выражения В0к , В1к из (16), (17) и проделав некоторые преобразования, получим следующие обыкновенные дифференциальные уравнения второго порядка относительно Як и zk

ЯкЯк + -Я2 - ^ - Рк Рда 2 к 4 рп

2а р0 Як

да СВ т:

= Ъ-к^к

2* У+1 И у+2 у=0 *к3-к а

да СВ

Е^у^ у 3-к

2* У+1 И у+2 у=0 *к3-к а

♦I

У=2

(т+1) В2

2ЯІ7+4

р

0

0

у=0

С 0 В 3 . С , В 3 к т3 к

у0 у3-к у1 у3-к 3-к

7У+1

*К+-кИ -2

I

дСУдССд£У3_к£С3-кЯГ2

9„у+С ,/у+С+2д+2 у,С, д=0 А!>к3-к а

■ V 0к + Лк

як тк + 3ТкЯк =

36В1кВ2 к 5Я7

да 07+1д+1 В В да

+ I уд ук дк +1

/_і я ї+д+4 <

у, д=2

у=0

С.

3Ву3-кЯк + Ву3-кЯк

*у+1 Иу+2 Дк 3-к“

2СУ2Ву3-кЯк (^^к - 3Т3-к )

I у+3

2©дд+1ССдВУквс3-кЯгу^3

(19)

С+д лС+д+1

12 , Т2, Я 1, Т1, Я 2 , Т2 ,

2) ■

тк = ёк (Я1, (1, Я2, т2, Я1, Т1, Я2, Т2 ) .

Для этого нужно найти коэффициенты Вук из (14) и их производные Вук из (16), (17) и уравне-

ний

Вм =!

пС„

у=0

(п +1)

(Я п+1 ву 3-к

7У+п+1

Л'

(20)

В0а) = --?кЯк2, а > 0; В1(к0)=-

В(а) = Як

^^к-I

г- С.., В.

(а-1) Л

7=0 ‘*к 3-к

у1 у 3-к

у+1

И1+2

Я Тк 2 ,

, а > 1;

(21)

Впк = 0 , 0 < а < п,

(а-1)

„г>2 п+1 а-1 С В '

(а) пЯк уп^у3-к

в;г =

т I

А 0

у+піу+п+1 ’ 7=0 *к 3-кИ

а > п+1,

(22)

-5

У, С ,5=0

“ ©5®С СВ В Я5+9-2

+ у '-Ч9^У5^СЭДУ 3-кД С 3-к к + у

^ „у+5+С+Э йУ+С+5+9+2 +Ч'1к-

у, 5, С ,9=0 Лк 3-к Ы

В (18), (19) введены поправки у0 к, у к, Ак (к = 1, 2), учитывающие эффекты вязкости (у0 к, у к) и сжимаемости (Дк) жидкости. Они будут определены ниже.

4. Метод решения

Система (18), (19) решается численно высокоточным методом Дормана-Принса с переменным шагом интегрирования по времени [6]. Для его применения уравнения (18), (19) нужно представить в каноническом виде

Як = /к (К ^ Я2

где к = 1, 2, п > 2, а верхним индексом в круглых скобках обозначен номер приближения.

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

уравнений (16), (17), (20).

Таким образом, для получения решения ^-го порядка точности относительно 5 при определении коэффициентов Вук из (14) и их производных Вук из (16), (17), (20) методом последовательных

приближений делается N итераций.

5. Учет эффектов вязкости и сжимаемости жидкости, теплообмена между пузырьками и жидкостью

Эффект вязкости жидкости учитывается в предположении потенциальности движения жидкости [9]. В этом случае динамическое граничное условие (4) изменяется на р = рк -2(2цЯк +а)/Як , а в уравнениях (18),

(19) будем иметь (Лц = 1)

=

4у^

~Я~

Я

г. V СуА 2 л ?1-!

(23)

V 02 =■

к = 1, 2, п = 2, 3, ..., полученных дифференцированием последнего выражения (14) (штрих, как и точка, означает дифференцирование по ґ). Рассмотрим определение Вук и В,к более подробно.

Коэффициенты Вук находятся из (14) методом последовательных приближений. Он построен с учетом того, что, как следует из (14), Впк ~ 8”+' при п > 2, т.е. коэффициенты Впк при п > 2 имеют относительно 8 порядок п + 1, где 8 = (Яі + Я2) / И < 1. С учетом этого для определения решения с Ы-м порядком точности нужно выполнить N итераций по следующему алгоритму, реализующему метод последовательных приближений,

12 4л

Я

^ И1+2

у=0 и

4уЯ2

“ЯГ’

да О В Л

ч -I О1Ву1

у=0

(24)

И1+2

Влияние сжимаемости жидкости на динамику пузырьков учитывается без учета их взаимодействия друг с другом [5]. В результате в уравнениях для радиусов пузырьков будем иметь [10]

Л, = Я

(

ЯЛ

я„

4

Рк - Ро

Л

Р к - Р да

(25)

где с0 - скорость звука в жидкости.

Выражение для давления (5) соответствует либо чисто изотермическому (Гк = 1), либо чисто

адиабатическому (Гк = к, к - показатель адиабаты) пузырьку. В настоящей работе, следуя [11], Гк считается функцией числа Пекле, которое принимается в виде Рек = Якик/хк, где ик - характерная скорость радиального движения поверхности пузырька, хк = ^к / (рк Су), Хк - коэффициент теплопроводности газа в пузырьке, су -удельная теплоемкость газа при постоянном объеме, рк - плотность газа в пузырьке. Полагается, что Хк = ^0 (Тк / T0)0.5, рк = р0,к (Я0, к / Яk)3, р0,к = Р0,к / Я ТО, где 7 - температура жидкости, Тк - температура газа, ^0 - теплопроводность газа при температуре Т0, р0, к - начальная плотность газа, Р0,к - начальное давление газа, Яё - газовая постоянная. В качестве характерной скорости ик принимается характерная скорость свободных незатухающих колебаний адиабатического пузырька

и =

Я2

2 Р0

3р0

1-

1

к-1

ґ Я Л3

Я0,

V Як у

Ґ я Vк-3

Я0,к

V Як у

-к

Зависимость Гк от числа Рек аппроксимируется выражением [11]

Г к (Ре к ) = 1 + (к- 1)ехр [-5.8 (Рек )-0'6 ^.

В результате давление газа в пузырьках оп-

ределяется выражением

Рк

Р0 +

2а

Я0, к

3

Я,

V Як у

(26)

где температура газа Тк находится из уравнения

Тк =

3 (1 -Гк ) ЯкТ± + 3хк (7-тк)

Я

(27)

сжимаемость - согласно (25) в предположении Zk = 0, а пузырьки предполагались адиабатическими (Гк = к).

Ч “к

При использовании (26), (27) вместо (5) в правых частях уравнений для радиусов пузырьков (18) возникает зависимость от температуры Тк, а к системе уравнений (18), (19) добавляются уравнения (27). При этом описанный выше алгоритм численного решения остается неизменным.

Система уравнений (18), (19), (27), в которой коэффициенты Ву к определяются выражениями (14), а их производные Вук - выражениями (16),

(17), (20), поправки на вязкость и сжимаемость жидкости у0 к, у1 к, Дк - выражениями (23)-(25), а давление газа в пузырьках рк - равенством (26), и представляет собой предлагаемую модель взаимодействия пузырьков.

6. Верификация метода расчета

Для проверки правильности работы алгоритма расчета было проведено сравнение с результатами работы [5], где использовалась модель третьего порядка точности по 5, вязкость жидкости учитывалась согласно (23), (24) при Лц = 1,5,

£ со / 2п

Рис.2. Изменение положений центров двух пузырьков на линии их взаимодействия в условиях работы [5]. Штриховые линии - результаты настоящей работы, сплошные - результаты работы [5].

На рис.2 приведено сравнение изменений координаты центров пузырьков, рассчитанных с третьим порядком точности N = 3) по методике настоящей работы (штриховые кривые), с результатами работы [5] (сплошные кривые).

Входные данные задачи следующие:

с0 = 1500 м/с, р0 = 998 кг/м3, р0 = 1 бар,

а = 0,0725 Н/м, ц = 10-3 кг/(м с), к = 1,4, Др = 1,2 бар, ю/2л = 20 кГц. В начальный момент времени Я1 = Я2 = 0, Zl = ¿2 = 0 , Я1 = 5 мкм, Я2 = 7 мкм, z1 = 200 мкм, z2 = 0 (й = 200 мкм). Как видно, согласование решений настоящей работы и работы [5] весьма удовлетворительное.

7. Моделирование взаимодействия пузырьков Некоторые характерные особенности взаимодействия пузырьков в рассматриваемых ниже задачах демонстрируются на примере взаимодействия пузырьков одинаковых размеров. Результаты расчетов для одного периода колебаний давления жидкости представлены на рис.3. В верхней части рис.3 приведены временные зависимости радиусов пузырьков и параметра Гк, а в нижней - временные зависимости положений центров пузырьков на оси z и величины малого параметра 5. Решение получено с третьим ^ = 3) порядком точности относительно 5 при следующих входных данных: с0 = 1500 м/с,

р0 = 998 кг/м3, р0 = 1 бар, Т0 = 293,15 К,

а = 0,0725 Н/м, ц = 10-3 кг/(м с), к = 1,4,

су = 720 Дж/(кг К), ^0 = 0,0258 Вт/(м К),

Яё = 287 Дж/(кг К), Др = 1,2 бар, ю/2л = 20 кГц. При t = 0: Я1 = Я2 = 3 мкм, Я1 = Я2 = 0,

z1 = 300 мкм, z2 = 0 (й = 300 мкм), ^ = :^2 = 0 .

Рис.3. Изменение Як - радиуса пузырьков, zk - координат их центров, Гк - показателя политропы газа в пузырьках, зависящего от числа Пекле, и 5 - параметра, характеризующего близость пузырьков в задаче взаимодействия двух одинаковых пузырьков.

Большую часть времени пузырьки при взаимодействии ведут себя как изотермические (Гк « 1). Вместе с тем, кратковременно в окрестностях локальных экстремальных сжатий, возникающих в ходе радиальных колебаний пузырьков, следующих за их максимальным расширением, поведение пузырьков становятся близкими к адиабатическому (величина Гк заметно превышает 1, приближаясь к к = 1,4). Это означает, что пузырьки нельзя рассматривать как чисто изотермические, или адиабатические.

Для пространственных перемещений пузырьков характерно, что при расширении они смещаются значительно меньше и медленнее, чем при сжатии. Наибольшие смещения наблюдаются в ходе первого наиболее сильного сжатия, особенно в его заключительной высокоскоростной стадии, где смещение составляет около 6 мкм, что в 2 раза больше равновесного радиуса. Последующие пространственные перемещения в ходе затухающих радиальных колебаний являются менее выраженными, а к концу периода практически исчезают. Параметр 5, характеризующий близость взаимодействия, принимает максимальное значение (около 0,16) в момент максимального расширения пузырьков. Далее приводятся два примера при входных данных и начальных условиях, близких к тем, что использовались в предыдущем случае.

0 10 20 30 40 50

I ш / 2п

Рис.4. Изменение координат центров пузырьков и параметра 5 в задаче взаимодействия двух сферических пузырьков. Период пространственных колебаний пузырьков в паре совпадает с периодом внешнего возбуждения. Штриховые линии - решение, полученное относительно 5 с порядком точности N = 3, пунктирные - с порядком точности N = 4, сплошные - с N = 5.

На рис.4 представлены временные зависимости положения пузырьков на оси z и параметра 5 в ходе взаимодействия пузырьков на протяжении 50 периодов колебаний давления жидкости. В начале взаимодействия Я1 = 2,5 мкм,

Я2 = 4,5 мкм. Приведены решения, полученные относительно 5 с порядком точности N = 3 (штриховые кривые), 4 (пунктирные кривые) и 5 (сплошные кривые). Кривые с N > 5 совпадают. В ходе первых 15 периодов пузырьки находятся еще относительно далеко друг от друга, так что их взаимодействие можно описывать моделью третьего порядка. Затем примерно через 3 периода пузырьки формируют связанную пару, которая в последующем смещается как единое целое. При N = 3 расстояние между пузырьками в паре описывается удовлетворительно, а направление их смещения - неверно (пузырьки движутся противоположно тому, что должно быть). При N = 4 скорость смещения пары получается правильной по направлению, но неверной по величине. И лишь модели с N > 5 дают полностью правильное решение. Максимальное значение параметра 5 при сближении пузырьков и формировании связанной пары возрастает от 0,15 до

0,55, а затем остается постоянным.

На рис.5 представлены результаты для начальных значений Я1 = 3 мкм, Я2 = 4,2 мкм, z1 = 600 мкм, z2 = 0 (й = 600 мкм). Сходимость в этом примере достигается при N > 6. Пузырьки здесь также сначала сближаются, а затем образуют связанную пару. Однако теперь пузырьки в связанной паре колеблются не с периодом, равным периоду колебаний давления жидкости, а с периодом, в два раза превышающим его.

t со / 2тг

Рис.5. Изменение координат центров пузырьков и параметра 5 в задаче взаимодействия двух сферических пузырьков. После образования пары пузырьки двигаются с периодом в два раза превышающим период колебаний давления жидкости. Решение получено с шестым (Ы = 6) порядком точности относительно 5.

8. Заключение

Предложены математическая модель и метод расчета задач близкого взаимодействия двух сферических пузырьков газа в жидкости, когда точности имеющихся в литературе моделей третьего порядка относительно безразмерного параметра 5 = (Я1 + Я2) / й < 1 (Я1, Я2 - радиусы пузырьков, й - расстояние между их центрами) оказывается недостаточно. Порядок точности предлагаемой модели относительно 5 произволен. В ней учитываются эффекты вязкости и сжимаемости жидкости, влияние теплообмена пузырьков с окружающей жидкостью. Предлагаемая модель представляет собой систему обыкновенных дифференциальных уравнений второго порядка относительно радиусов пузырьков и пространственного положения их центров и первого порядка относительно температуры газа в пузырьках. В уравнения для радиусов пузырьков и координат их центров входят коэффициенты представления потенциала скорости жидкости в виде суммы частных решений уравнения Лапласа в сферических координатах пузырьков. Для иллюстрации применимости предлагаемой модели приведены два примера, различающиеся характером пространственных перемещений пузырьков после образования связанной пары. В первом примере период пространст-

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

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

1. Konovalova S.I., Akhatov I.Sh. Structure formation in acoustic cavitation // Multiphase Science and Technology. - 2005. - Vol.17. - №3. - Р.343-371.

2. Bjerknes V.F.K. Field of Force. - New York: Columbia University Press, 1906. - 173 p.

3. Кузнецов Г.Н., Щукин И.Е. Взаимодействие пульсирующих пузырьков в вязкой жидкости // Акустический журнал. - 1972. - T.18. - №4. -С.565-570.

4. Harkin A., Kaper T.J., Nadim A. Pulsation and translation of two gas bubbles // Journal of Fluid Mechanics. - 2001. - Vol.445. - Р.377-411.

5. Doinikov A. A. Translational motion of two interacting bubbles in a strong acoustic field // Physical Review E. - 2001. - Vol.64. - №2. -026301 (6 p.).

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

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

- С.521-533.

8. Hobson E.W. The theory of spherical and ellipsoidal harmonics. - Cambridge: Cambridge University Press, 1931. - 476 p.

9. Gavrilyuk S.L., Teshukov V.M. Drag force acting on a bubble in a cloud of compressible spherical bubbles at large Reynolds number // European Journal of Mechanics B/Fluids. - 2005. - Vol.24. - №4. -Р.468-477.

10. Doinikov A.A. Equations of coupled radial and translational motions of a bubble in a weakly compressible liquid // Physics of Fluids. - 2005. -Vol.17. - №12. - 128101 (4 p.).

11. Hilgenfeldt S., Grossmann S., Lohse D. Sonoluminescence light emission // Physics of Fluids.

- 1999. - Vol.11. - №6. - Р.1318-1330.

INTERACTION OF TWO SPHERICAL GAS BUBBLES IN A LIQUID IN AN ACOUSTIC FIELD

A.A.Aganin, A.I.Davletshin

A mathematical model and a technique for solving problems of close-range interaction of two spherical gas bubbles in a liquid with allowing for the effects of the liquid viscosity and compressibility, heat transfer between the bubbles and liquid are proposed. The model is a system of ordinary second order differential equations for the radii of the bubbles and their spatial positions and the first order ones for the gas

temperature in each bubble. The equations for the radii of the bubbles and their spatial positions contain coefficients of expansion of the fluid velocity potential in Legendre polynomials, which are determined from a system of linear algebraic equations. The order of accuracy of the proposed model relative to the ratio of sum of the radii of interacting bubbles to the distance between their centers is arbitrary. Some examples of application of the proposed model to describing close-range interaction between the bubbles are presented.

Key words: bubble interaction, acoustic field, velocity potential.

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

E-mail: aganin@kfti.knc.ru

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

E-mail: davanas@mail.ru

Поступила в редакцию 07.05.2011