О влиянии сил взаимодействия капель на динамику эмульсии

ВЕСТНИК ЮГОРСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

_2019 г. Выпуск 1 (52). С. 59-65_

DOI: 10.17816/byusu20190159-65 УДК 621.892.2

Т. В. Пронькина

О ВЛИЯНИИ СИЛ ВЗАИМОДЕЙСТВИЯ КАПЕЛЬ НА ДИНАМИКУ ЭМУЛЬСИИ

Исследована динамика осаждения составных капель эмульсии под действием силы тяжести. Учитывается взаимодействие капли с ее включением. Рассмотрены осесимметрич-ная и асимметричная задачи осаждения составных капель. Найдены выражения для относительных и абсолютных скоростей капель составной эмульсии. На основе численного моделирования получены траектории относительного и абсолютного движения капель.

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

T. V. Pronkina

ABOUT THE INFLUENCE OF THE FORCES OF INTERACTION BETWEEN THE DROPLETS ON THE DYNAMICS OF EMULSION

The dynamics of deposition of compound droplets of the emulsion under the action of gravity is investigated. The interaction of the droplet with its inclusion is taken into account. Axisymmetric and asymmetric problems of deposition of compound droplets are considered. Expressions for the relative and absolute velocities of the compound emulsion droplets are found. Based on numerical modeling, the trajectories of the relative and absolute motion of the droplets are obtained.

Keywords: emulsion, composite droplet, hydrodynamic interaction, relative velocity, absolute velocity.

Введение

В последние годы возрос интерес к направлениям исследования, связанным с созданием материалов, обладающих заданными характеристиками: наноэмульсии и составные эмульсии. Составная эмульсия представляет собой дисперсную среду, в которой капли вязкой жидкости содержат в себе включения в виде других капель. Количество таких включений может быть различным - от нескольких капель до сотен включений в одной капле [1, 2]. Наличие структуры в каплях эмульсии делает более сложным моделирование и прогнозирование свойств таких систем.

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

В эмульсии, как в любой многофазной системе (суспензии, аэрозоли), существует два различных механизма взаимодействия капель. Один из них связан с силами взаимодействия между самими каплями. Другой механизм связан с взаимодействием капель через движение несущей фазы (жидкость, газ): распределение скорости и давления жидкости вблизи какой-либо капли зависит от расположения других капель. Движение одной капли влияет на движение всех остальных и наоборот. Такое взаимодействие называется гидродинамическим и

оказывает влияние на все процессы, происходящие в жидких дисперсных системах, в том числе и в эмульсии [3-5].

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

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

В неограниченной несжимаемой жидкости, вязкость которой равна , находится жидкая частица А вязкости 7]а. Частица А, в свою очередь, содержит жидкое включение В вязкости 7ь. Обозначим через а размер капли А и через Ъ размер включения В. Размеры капель А и В достаточно малы, чтобы уравнения движения жидкости в каждой из фаз были линейными. Для определения скорости и и давления р в несущей жидкости и внутри каждой из капель записывались соответствующие уравнения в приближении Стокса с граничными условиями: 1) равенства скоростей жидкостей на границе капли А ; 2) условие равенства касательных напряжений на поверхности капли А; 3) условие непротекания несущей жидкости внутрь частицы А. Аналогичные граничные условия были записаны на поверхности включения В. Кроме того, учитывался тот факт, что скорость возмущения обращается в нуль при бесконечном удалении от капель.

В результате были получены распределения скорости и давления вокруг частиц [6-8], которые далее использовались для вычисления сил Стокса.

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

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

где и - скорость капли , иии - скорость капли , , , , - функции безразмерных параметров ] ,/]ь, т]а/]ь, а/Ъ, г /Ъ.

Также на капли действуют сила тяжести и сила Архимеда:

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

Осесимметричная задача об осаждении составных капель

4

Д,А = -7г(а3(ра - рг) + Ь3 (рь - ра))д, 4

= -ттЬ3(рь — Ра)9>

где р I - плотность несущей жидкости, ра - плотность капли А, р й - плотность капли В.

Скорости капель К и 1/К определяются из равенства нулю сил, действующих на них (квазистатическое приближение).

(рА +рА = 0

и +: О

Для определения закона относительного движения капель г ( необходимо решить дифференциальное уравнение

Нг _,

Ц-= / 1 I-/? I I

Введем систему координат с началом отсчета в центре капли А. Пусть направление оси Ох совпадает с направлением вектора д. Находим проекцию уравнений системы (1) на ось О х

-6щ га (ДАУ 1 1 + /2А(Ш 1 1 - V11)) + ±п(а3 (ра - рг) + Ъ3(рь - ра) )д = О,

-6щаЬ{Г*у\\ + /2В(ЖИ - 7И)) + \пЬ\ръ - ра)д = 0.

Решая систему (2) относительно величин

и

, получаем

2д

7И =

ЖИ = -

9а(/2лАв-/1л/2>а^

~ ЙШрг ~ РаХа3 ~ Ь3) + аЬ2(Ал - /2л)?7г(Рь ~ Ра) 9а(//Ав - ААГ2в)ПаШ '

Таким образом, имеем обыкновенное дифференциальное уравнение

^ = 2д

М " 9а(///в - Г^Лат

Уравнение решалось численным способом при различных начальных условиях.

Рисунок 1 - Относительное движение включения В

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

Асимметричная задача осаждения составных капель

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

пА _ р\\А р±А

1 I 11 ~ 1 I '

рв _ р\\В , р1В

1 I 11 ~ 1 I '

и\\А и\\в

где , ,

, вычисляются по формулам

Здесь /Д Дв, /2Л, /25, /3Л, /35, /4Л, /4Л - функции безразмерных параметров ] г/]ь, ]а/]ь, а/Ъ,г/Ъ.

Как и в осесимметричной задаче, на капли действуют сила тяжести и сила Архимеда:

4

^ = -тг( а 3(ра - Р г) + Ъ3(рь - Ра) )#,

Я5 =-лЬ3(рь - ра)д.

Для определения скорости капель и и иии будем решать систему уравнений

+ £в = 0 .

Введем систему координат с началом отсчета в центре капли . Пусть направление оси совпадает с направлением вектора , а ось направлена перпендикулярно оси . Находим проекции уравнений системы (3) на оси О X! и О х2, учитывая, что

Ы I = II , 71 = £171, I = ^7 | | , 7± = £271,

^ 5 ^ 5 А2 у. 5 Л-2 '

ж£ 1 = ^ж 1 ^ ж/ = ^ж\ ж£ 1 = -ж 1 1 , ж/ 1

Х-1 у, ^ Л1 у ^ Х-2 у ^ л2 у

Система принимает следующий вид

-6Щ1а(Г1АУ" ^ + //(Ж11 - К11) ^ + /з^11 + /Д И'1 - |) + + 4тт( а 3(Ра -Р г) + Ъ3(Рь -Ра) )# = 0,

бтт^аСД^П ^ + /Д^П - К") ^ + /3ЛКХ21 + /Ди^ - = 0,

-вщаЪ{-^у\\^ + - 7") ^ - /3ВКХ ^ + /ДИ^

^7) +

+ -7гЬ3(рь - Ра)д = 0,

-еЩаъ{-^у\\| + /Д^п - к") | _ ^ + /Ди^ - = о.

Решая систему (4) относительно величин , , , , учитывая при этом, что

получаем

7" = -

2дгх1

,3 ¿-в

9а(//Лв - ААГ2в)г211а111 ж 1 1 = ( 2 ДТХ!(а 3(ЛВ - /2В) 7 а (Р г - Ра) + Ъ3(ДВ - /2В) 7 а (Ра - Рь) ) + +аЬ2(/1л - Г2аШРь - ра))/(9а(/2л/1в - ^Й)г2ЦаГ}1),

2дгх2

,3 ¿в

2 сА^

9а(///3в - ЬАивУ2ЛаЛ1 V 1 = (2дгх2(а- //)77 а(рг - ра) + Ъ3(П - //)77 а(ра - ръ)) +

+аЬ2(ЬА - иАШрь - ра))/(9а(///3в - ^ПУ2^). Чтобы определить относительную скорость, запишем систему обыкновенных дифференциальных уравнений

— = УК11 — + IV1 — — уЧ — — V1 — М г г г г '

— = ук" — — V]/1- — — уН — + V1 —

М г г г г '

где г

= 4*2 +

Хп

Из системы (5) легко получить уравнение траектории относительного движения включения В

^ х т/1/11 т/1/1 т/11 т/1 £1

У\/\\ ±2. _ у]/1 *1 _ у II *2 +

Г Г Г Г

г г г г

Данное уравнение также решалось численным способом при различных начальных условиях. Ниже приведены графики траектории относительного движения включения при различных значениях параметров а/Ъ, 71/ца, Р1/ра, Цъ/Ца, Рь/Ра-

В начальный момент времени центр частицы В находится в точке с координатами х± = О , х2 = 0 . 1. Начало координат находится в центре капли А.

Пусть а = 5 , — = 0 . 5 , — = 0 . 7 , — = 1. 5. Значение параметра — будем изменять.

Ь Т]а ра Т]а Ра

Сплошной линией изображена траектория движения при — = 5, пунктирной - при

Ра

— = 2. 5 , штрихпунктирной - при — = 0 . 5.

Ра Ра

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

0.1002

0.10015

0.1001

0.10005

0.1 0.2 0.3 0.4 0.5

Рисунок 2 - Траектории относительного движения капли В

Х2

Рисунок 3 - Траектории абсолютного движения капли В

Х2

Рисунок 4 - Траектории абсолютного движения капли А

Выводы

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

Литература:

1. Muschiolika, G. Multiple emulsions for food use [Text] / G. Muschiolika // Current Opinion in Colloid and Interface Science. - 2007. - V. 12, Is. 4-5. - P. 213-220.

2. Pal, R. Viscosity models for multiple emulsions [Text] / R. Pal // Food Hydrocolloidse. -2008. - V. 22, Is. 3. - P. 428-438.

3. Волны в жидкостях с пузырьками [Текст] / А. А. Губайдулин, А. И. Ивандаев, Р. И. Нигматулин [и др.] // Итоги науки и техники. Механика жидкости и газа. ВИНИТИ. - 1982. -Т. 17. - С. 160-249.

4. Нигматулин, Р. И. Основы механики гетерогенных сред [Текст] / Р. И. Нигматулин. -Москва : Наука, 1978. - 336 с.

5. Нигматулин, Р. И. Динамика многофазных сред [Текст] : в 2 ч. / Р. И. Нигматулин. -Москва : Наука, 1987. - Ч. 1. - 464 с.

6. Нигматулин, Р. И. Динамика многофазных сред [Текст] : в 2 ч. / Р. И. Нигматулин. -Москва : Наука, 1987. - Ч. 2. - 360 с.

7. Мартынов, С. И. Моделирование взаимодействия капель в линейном потоке и вязкость эмульсии [Текст] / С. И. Мартынов, Т. В. Пронькина // Информационные системы и технологии. - 2010. - № 3. - С. 86-90.

8. Мартынов, С. И. Составная капля эмульсии в однородном потоке вязкой жидкости [Текст] / С. И. Мартынов, Т. В. Пронькина // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. Математика. - 2010. - № 2. - С. 85-93.

9. Пронькина, Т. В. Обтекание составной капли эмульсии линейным потоком вязкой жидкости [Текст] / Т. В. Пронькина // Вестник Югорского государственного университета. -2009. - Вып. 2 (13). - С. 77-79.