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

Автор: Анас Ильгизович Давлетшин

Организация: Институт механики и машиностроения ФИЦ Казанский научный центр РАН

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

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

Для исследования взаимодействия пузырьков наиболее эффективными являются модели, в которых разрешающие соотношения представляются в виде системы обыкновенных дифференциальных уравнений (ОДУ). При получении таких уравнений используется метод разложения по малому параметру R / D, где R – характерный радиус пузырьков, D – характерное расстояние между их центрами. При таком подходе по мере уменьшения расстояния между пузырьками для адекватного описания их совместной динамики нужно использовать модели, имеющие все более высокий порядок точности относительно указанного параметра. Для исследования особенностей взаимодействия подвижных пузырьков в основном применяются модели от третьего до пятого порядка точности относительно R / D. В большинстве из этих моделей пузырьки считаются сферическими.

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

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

 

Работа выполнена за счет гранта РНФ (проект № 21-11-00100, https://rscf.ru/project/21-11-00100/).

 

1. Aganin A.A., Davletshin A.I. Equations of interaction of weakly non-spherical gas bubbles in liquid // Lobachevskii Journal of Mathematics. 2018. V. 39. № 8. P. 1047-1052.

2. Luther S. Theoretische Beschreibung und experimentelle Untersuchung raum-zeitlicher Strukturbildung in akustischen Kavitationsblasenfeldern. – Niedersächsische Staats- und Universitätsbibliothek, 2000.

3. Xu Y., Butt A. Confirmatory experiments for nuclear emissions during acoustic cavitation // Nuclear Engineering and Design. 2005. V. 235. P. 1317-1324.