Научная статья на тему 'ВЛИЯНИЕ ПРОСТРАНСТВЕННОЙ НЕОДНОРОДНОСТИ ПОДЛОЖЕК И ЭЛЕКТРИЧЕСКОГО ПОЛЯ НА ДИНАМИКУ ЗАЖАТОЙ КАПЛИ'

ВЛИЯНИЕ ПРОСТРАНСТВЕННОЙ НЕОДНОРОДНОСТИ ПОДЛОЖЕК И ЭЛЕКТРИЧЕСКОГО ПОЛЯ НА ДИНАМИКУ ЗАЖАТОЙ КАПЛИ Текст научной статьи по специальности «Физика»

CC BY
27
6
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ВЫНУЖДЕННЫЕ КОЛЕБАНИЯ / ЭЛЕКТРОСМАЧИВАНИЕ / КОНТАКТНАЯ ЛИНИЯ / КРАЕВОЙ УГОЛ

Аннотация научной статьи по физике, автор научной работы — Алабужев А. А., Пьянкова М. А.

В данной статье исследуются вынужденные колебания капли несжимаемой жидкости под действием неоднородного переменного электрического поля. Капля окружена несжимаемой жидкостью другой плотности и зажата между двумя неоднородными параллельными пластинами. В равновесии капля имеет форму круглого цилиндра, ограниченного в осевом направлении этими пластинами. Внешнее электрическое поле действует как внешняя сила, вызывающая движение линии контакта. Для описания движения контактной линии используется модифицированное граничное условие Хокинга: скорость движения линии контакта пропорциональна отклонению краевого угла от его равновесного значения и скорости быстрых релаксационных процессов, частота которых пропорциональна удвоенной частоте электрического поля. Использование этого уравнения позволяет качественно описать экспериментальную зависимость краевого угла в зависимости от напряжения в отличие от уравнения Юнга-Липмана. Решение задачи представляется в виде ряда Фурье по собственным функциям оператора Лапласа. Полученная система неоднородных уравнений для неизвестных амплитуд решается численно. Построены графики амплитудно-частотных характеристик и динамической формы капли при различных значениях параметров задачи. Параметр смачивания зависит от полярного угла a , т.е. коэффициент взаимодействия между пластиной и жидкостью (линия контакта) является функцией координат.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

INFLUENCE OF SURFACE PROPERTIES ON THE DYNAMICS OF A CLAMPED DROP IN AN INHOMOGENEOUS ELECTRIC FIELD

This article investigates the forced oscillations of an incompressible liquid drop under the action of an inhomogeneous alternating electric field. The drop is surrounded by an incompressible fluid of a different density and is clamped between two inhomogeneous parallel plates. In equilibrium, the drop has the shape of a round cylinder bounded in the axial direction by these plates. The external electric field acts as an external force causing the contact line to move. To describe the motion of the contact line, a modified Hocking boundary condition is used: the velocity of the contact line is proportional to the deviation of the contact angle from its equilibrium value and the rate of fast relaxation processes, the frequency of which is proportional to the doubled frequency of the electric field. The use of this equation makes it possible to qualitatively describe the experimental dependence of the contact angle as a function of stress, in contrast to the Young-Lippmann equation. The solution of the problem is represented as a Fourier series in terms of eigenfunctions of the Laplace operator. The resulting system of inhomogeneous equations for unknown amplitudes is solved numerically. Graphs of the amplitude-frequency characteristics and the dynamic shape of the drop are plotted for various values of the problem parameters. The wetting parameter depends on the polar angle a , i.e. the coefficient of interaction between the plate and the fluid (contact line) is a function of the coordinates.

Текст научной работы на тему «ВЛИЯНИЕ ПРОСТРАНСТВЕННОЙ НЕОДНОРОДНОСТИ ПОДЛОЖЕК И ЭЛЕКТРИЧЕСКОГО ПОЛЯ НА ДИНАМИКУ ЗАЖАТОЙ КАПЛИ»

ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА

2022

• ФИЗИКА •

Вып. 2

УДК 532.6, 532.5 PACS 47.55.D-, 47.55.dr

Влияние пространственной неоднородности подложек и электрического поля

А. А. Алабужев1'2, М. А. Пьянкова1*

1 Институт механики сплошных сред УрО РАН, Пермь, Россия

2 Пермский государственный национальный исследовательский университет, Пермь, Россия

* [email protected]

В данной статье исследуются вынужденные колебания капли несжимаемой жидкости под действием неоднородного переменного электрического поля. Капля окружена несжимаемой жидкостью другой плотности и зажата между двумя неоднородными параллельными пластинами. В равновесии капля имеет форму круглого цилиндра, ограниченного в осевом направлении этими пластинами. Внешнее электрическое поле действует как внешняя сила, вызывающая движение линии контакта. Для описания движения контактной линии используется модифицированное граничное условие Хокинга: скорость движения линии контакта пропорциональна отклонению краевого угла от его равновесного значения и скорости быстрых релаксационных процессов, частота которых пропорциональна удвоенной частоте электрического поля. Использование этого уравнения позволяет качественно описать экспериментальную зависимость краевого угла в зависимости от напряжения в отличие от уравнения Юнга-Липмана. Решение задачи представляется в виде ряда Фурье по собственным функциям оператора Лапласа. Полученная система неоднородных уравнений для неизвестных амплитуд решается численно. Построены графики амплитудно-частотных характеристик и динамической формы капли при различных значениях параметров задачи. Параметр смачивания зависит от полярного угла а , т.е. коэффициент взаимодействия между пластиной и жидкостью (линия контакта) является функцией координат.

Ключевые слова: вынужденные колебания; электросмачивание; контактная линия; краевой угол

Поступила в редакцию 23.05.2022; после рецензии 09.06.2022; принята к опубликованию 10.06.2022

Influence of surface properties on the dynamics of a clamped drop in an inhomogeneous electric field

A. A. Alabuzhev1'2, M. A. Pyankova1*

1 Institute of Continuous Media Mechanics UB RAS, Perm, Russia

2 Perm State University, Perm, Russia * [email protected]

This article investigates the forced oscillations of an incompressible liquid drop under the action of an inhomogeneous alternating electric field. The drop is surrounded by an incompressible fluid of a different density and is clamped between two inhomogeneous parallel plates. In equilibrium, the drop has the shape of a round cylinder bounded in the axial direction by these plates. The external electric field acts as an external force causing the contact line to move. To describe the motion of the

© Алабужев А. А., Пьянкова М. А., 2022

на динамику зажатой капли

распространяется на условиях лицензии

Creative Commons Attribution 4.0 International (CC BY 4.0).

contact line, a modified Hocking boundary condition is used: the velocity of the contact line is proportional to the deviation of the contact angle from its equilibrium value and the rate of fast relaxation processes, the frequency of which is proportional to the doubled frequency of the electric field. The use of this equation makes it possible to qualitatively describe the experimental dependence of the contact angle as a function of stress, in contrast to the Young-Lippmann equation. The solution of the problem is represented as a Fourier series in terms of eigenfunctions of the Laplace operator. The resulting system of inhomogeneous equations for unknown amplitudes is solved numerically. Graphs of the amplitude-frequency characteristics and the dynamic shape of the drop are plotted for various values of the problem parameters. The wetting parameter depends on the polar angle a , i.e. the coefficient of interaction between the plate and the fluid (contact line) is a function of the coordinates.

Keywords: forced oscillations; electrowetting; contact line; contact angle

Received 23 May 2022; revised 09 June 2022; accepted 10 June 2022

doi: 10.17072/1994-3598-2022-2-56-65

1. Введение

Движение капли жидкости на подложке встречается во многих промышленных процессах, в связи с чем привлекает большое внимание исследователей [1-10]. Особую сложность в этой области вызывает описание движения линии контакта трех сред и изменения краевого угла [11-19]. Изменение угла контакта под воздействием внешней силы является одним из способов изменения смачивания твердой поверхности жидкостью. Роль такой силы может выполнить электрическое поле - электросмачивание (electrowetting, EW). Если проводящая пластина (подложка) покрыта слоем диэлектрика, то это - электросмачивание на диэлектрике (electrowetting-on-dielectric, EWOD). Существует множество моделей, описывающих изменение краевого угла при электросмачивании, которые, так или иначе, основываются на моделях движения линии контакта начиная с первой модели - уравнения Юнга-Липмана [20-22]. Общее обсуждение эффективных граничных условий EWOD можно найти в статьях [23-30].

В данной работе мы продолжаем исследования поведения зажатой капли в переменном электрическом поле, начатые в работах [31-34]. Собственные колебания капли исследовались ранее в нескольких работах: как одинаковые [35] и разные [32] однородные пластины, так и одинаковые [36] и разные [34] неоднородные пластины. Капля находится в плоском конденсаторе, т.е. зажата между двумя разными неоднородными проводящими подложками, покрытыми диэлектриком. Внешнее электрическое поле играет роль движущей силы, действующей на линию контакта, и приводит к динамическому изменению краевого угла. Кроме того, мы предполагаем, что электрическое поле изменяется гармонически со временем по закону cos (a*t*). Следовательно, краевой угол

меняется по закону cos(2а*t*) согласно [20-22].

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

I1 = ±Л* (f^ A*cos(2aY)] , (1)

г*

где Q - отклонение поверхности капли от поло* к *

жения равновесия, z - осевая координата, Л -феноменологическая константа (так называемый параметр смачивания или параметр Хокинга [3-5, 13, 14, 35-38]), имеющий размерность скорости, A* - эффективная амплитуда, а* - частота электрического поля. Вторым слагаемым в граничном условии (1) является внешнее действие, которое записывается так же, как и в уравнении Юнга-Липмана [20-22]: cos 3 = Ew = EwV2, cos 30 = 0 в нашей задаче, таким образом,

с ~ ctg 3 = Ewlfi-E =

= EwV2 Д/1 - EWV4 *

* EwV2 + O(V6)* EwV2~ V2,

cos30 = а- (°ip -°cp ) ,

Ew = 0.5CVV-, C = ss0d- ,

30 - контактный угол без приложенного напряжения - равновесный краевой угол, который определяется уравнением Юнга, а - межфазное поверхностное натяжение между каплей проводящей жидкости (c), окружающей жидкостью (i) и диэлектрической поверхностью (p), Ew - электрокапиллярное число (оно представляет собой отноше-

ние электростатической энергии жидкости к межфазной энергии, окружающей жидкость), С - емкость на единицу площади, V - значение приложенного напряжения постоянного электрического тока, d - толщина диэлектрического слоя, е0 и е - коэффициенты диэлектрической проницаемости вакуума и подложки, соответственно.

Условие (1) включает в себя два важных предельных случая: 1) линия контакта закреплена и неподвижна (краевой угол может изменяться) -Л* = 0, 2) краевой угол постоянный (линия контакта свободно скользит по поверхности пластины) - 5^7& = 0 .

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

Исследуется динамика капли жидкости плотно* ^ * ^

сти р и кинематической вязкости у1 , зажатой в плоском бесконечном конденсаторе толщиной к* (рис. 1). Его пластины являются проводящими и покрыты изнутри диэлектрической пленкой. Предполагается, что линейные размеры пластин конденсатора много больше размеров капли и его края не влияют на ее динамику. Капля окружена жидкостью другой плотности ре* и вязкости V*. Равновесная форма капли описывается круглым цилиндром радиуса Л* и высотой к*; краевой угол 50 между боковой поверхностью капли и твердой поверхностью равен 0.5^ . К пластинам прикладывается переменное напряжение, имеющее частоту ю , что позволяет возбуждать внутри него электрическое поле. Внешнее переменное электрическое поле действует как внешняя сила, имеющая

л * /-ч *

амплитуду А и частоту 2ю , которая вызывает движение контактной линии.

Рис. 1. Геометрия задачи: 1 - электрод;

2 - диэлектрический слой

Для описания движения контактной линии используется модифицированное условие (1). Параметр смачивания Л* будет функцией координат, т.к. поверхность пластин неоднородная [34], и различным для каждой поверхности, т.к. пластины разные: Л*и при г = 0.5к* и Л* при г = -0.5к*.

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

Будем рассматривать цилиндрическую систему

**

координат г , а , г , исходя из симметрии задачи. Следовательно, поверхность капли описывается

* т-1* *-* / *

следующим уравнением: г = Л0 (а,г ,/ I.

Далее жидкости рассматриваются как невязкие и несжимаемые, т.е. скорость течения будет описываться уравнением Эйлера, а изменение плотности - уравнением неразрывности. Предполагая потенциальное движение жидкости, введем потенциал скорости V* = Уф . Отметим, что в системе присутствует диссипация энергии благодаря использованию граничного условия Хокинга [3-5, 13, 14, 35-38], несмотря на использование приближения невязкой жидкости. Это позволяет нам рассматривать течение как потенциальное даже при наличии деформируемой поверхности раздела.

Принимая в качестве характерных величин длины - Л*, высоты - к*, плотности - ре* + р*, времени -

U{Р +Р)R*3,

потенциала скорости -

А*^(р+р)л*3) ,

давления - А*а( Л0*) и отклонения поверхности -

А*, переходим к безразмерным переменным и получаем следующую линейную задачу:

Pj =-PjVjt' AVj =0, j = е '

Л 1 д ( д Л 1 д2 ,2 д2

А=--1 r— I+ ——- + Ь2 —

r дг \ дг J r да dz

r =1: \pr ] = Q, =V , [р]=С+Саа+ bQ ,

z = + -: p = 0, 2

r = 1, z =+ -: 2

(2)

(3)

(4)

С, = +Л„,4 {a){Cz + af {a)cos{2®,)), (5)

где p - давление жидкости, p - потенциал скорости, f {a) - функция пространственной неоднородности электрического поля, а - азимутальный угол, a - эффективная амплитуда электрического поля, С - отклонение поверхности от равновесного положения, КиЬ - безразмерный параметр Хо-

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

Краевая задача (2)-(5) включает в себя 6 безразмерных параметров: геометрический параметр Ь = К0И~1, соотношение плотностей:

* / * *\_1 * / * *\_1 Р, = Р, [Ре + Р, ) , Ре = Ре [Ре + Р ) ,

капиллярный параметр:

А = Л*а^[Р + р* ) % , частоту электрического поля:

* -1/2 /7 * ю=юа ^(p* + p

*) R

и амплитуду электрического поля:

а = 0.5 A'Ca^J (p* + p ) R* .

3. Вынужденные колебания

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

1. неоднородные пластины, однородное электрическое поле (далее НПОП):

Л„,ъ (а) = К,ъ |sin {cos (а)) , f (а) = 1,

соответственно;

2. однородные пластины, неоднородное поле (ОПНП):

К,ъ (а) = К*, f (а) = |sin (cos (а))|;

3. неоднородные пластины, неоднородное электрическое поле (НПНП):

Л„,ъ (а) = А^ |sin (cos (а)) , f (а) = |sin (cos (а)) .

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

Решение краевой задачи (2)-(5) будем искать в виде рядов Фурье по собственным функциям оператора Лапласа (2). С учетом осевой симметрии решение уравнения Лапласа (2) с граничным условием (4) можно записать в виде:

(pt (г,а, z,t) = Re(i2ae'2mty¡ (г,а, z)) , (6)

Wi (r, а, z) = £ £ (amkR^ (r) cos(2knz) +

m=0k=0

+ aik Rk (r) sin ((2k +1) nz)) cos (та), cpe (г,а, z, t) = Re (i2ae'2m,ye (г,а, z)) , (7)

¥e (r, а, z) = £ £ (ьттЦRk* (r) cos(2knz) +

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

m=0k=0

+ ът) Rmk: (r) sin (( 2k +1) nz )) cos (та),

где Rf (r) = Im (2krir), Rm1) (r) = Im ((2k +1) rir),

Rmk: (r) = K (2П ) , R^ (r) = Kü ((2k +1) пъг ) , Im и Km - модифицированные функции Бесселя m -го порядка первого и второго рода, соответственно. Мы используем разложение только по чётным азимутальным функциям, т.к. функции неоднородности являются чётными.

Кинематическое условие на поверхности раздела (второе граничное условие в (3)) определяет выражение для отклонения поверхности:

с {а, z, t) = Re (#(а, z) e'2at)

Z (а z) = £ £ (42 cos (2knz) -

(8)

sin ((2k +1) nz)) cos (mа) + d00) cos j ъ 1 +

+ d01) sin jЫ + di(0)z2 cos(а) + dfzcos(а) +

m=2

z

jj dm0) ch

Vm2 -1

+ d{m] sh

Vm -1

2

cos

(mа).

/j

Последние слагаемые представляют собой общее решение однородного уравнения - скачок давления на границе раздела (3).

Подставив решения (6)-(8) в систему (2)-(5), получим следующую систему для нахождения неизвестных амплитуд атк, Ьтк, Стк и ёт:

2ia

£(-1)k (^ ±41k К0 jd00) cos(¿j±d01) sin[¿Jj + ^1 [-4di0) + 2d«) +

+¿m

d<0) ch

^ -1]+ d« sh

л/n2 -1

(±d<0)+ d«) + ^n

V72"-!

= ,ъ jjj

j Г2~: 1

Bm

St j ^ * [ £) + ^ cos j 2ъ

(9)

±d7°°> sh

#-1

+ d7(1) ch

4¡r-\ 2ъ

+ aC_

; jj

+

2

+

ъ

> (0) (0) ,(0) (0) Лткг (1) ,(1) (1) К^г (1) (0) (0)

ЬЮ = -<*>, Ы, = ат —, ЬЦ = а ' -ткг— у» ^ ' - ^ '

(1)

(1)

тат, = Ст0 + ят<>0) е

^ (1)=С2+ , мх (1)=42+, с(?=-/£Ч(0)- ю ^,

с(0)=-„(0)d(0) с(0)=4ю¿т00dmo) _(0)=4®2¿т,dmo)

'-00 1-т0 „(0)2 . 2 ' тк

0)2 - 4ю2

т0

0)2 - 4ю2 т,

с( 1) = т,

4°02 = т(т2-1) , ^ =■

т

-1 + (2якЪ )2

(1) Ре лт0кге (1)

(1)

^тк2 =-

т

4ю2 ¿а dml) ^тг - 4ю2'

-1 + (2к +1)2 тт2Ъ2

Лт) (1)

Лт) (1)

р,лтг (1) ~Релт£ (1)

лог (г)=^ (г), лткг (г)=(г), (г)=(г), (г)=^ (г).

dr

dr

0.9

и 0.6

0.3 Н

2.1

Си

15 СО 20

1 1 I 15 0) 20

1.2

сч 0.8

0.4 —I

1 1 I 15 СО 20

dr

1.5

и 1

0.5-

10

"Г-1—I

15 СО 20

90-,

15 СО 20

д

Рис. 2. Зависимость максимального отклонения: а - боковой поверхности Сл б - контактной линии Си; в - контактной линии г - поверхности капли д, е - отклонения краевого угла у от частоты т для трех разных значений ХЪ (Ъ = 1.0, а=5.0, р^ = 0.7, Х„ = 1, ХЪ = 0.1). НПНП -пунктирная, ОПНП - штриховая, НПОП - сплошная линия

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

1.2-, У -0.60-0.6-

п—I—1—I—1—I—1—I -1.2 -0.6 0 0.6 х 1 -2

1.2-1 У -0.60-0.6-

~1—I—1—I—1—I—1—I -1.2 -0.6 0 0.6 х 1 -2

б

0.25

1-1-1-1-1

0.5 0.75 Т 1

а б в

Рис. 3. Форма контактной линии 1+0.1С(г,У (а, б) и значение краевого угла уи (в) на верхней поверхности в разные моменты периода колебаний капли (Ъ=1.0, а=5.0, р1 = 0.7, Х„ = 1, ХЬ = 0.1, т = 1). Линии на а, б: / = 0 - пунктирная, / = 0.125Т - штриховая, / = 0.25Т - штрих -пунктирная, / = 0.375Т - сплошная

0

0

0

а

в

0

г

е

0

Внешняя сила будет возбуждать азимутальные моды вынужденных колебаний капли, которые появляются из разложения в ряд функции неоднородности электрического поля /(а) (5). Следовательно, если поле и пластины однородные, то у нас будут возбуждаться только осесимметричные колебания. Если пластины одинаковые, т.е. Ли = Хь = Л, то будут только нечётные гармоники осесимметричной моды [31]; если разные, то и нечётные, и чётные [32].

Для дальнейшего описания динамики боковой поверхности введём следующие обозначения:

0.9

0.6

12-,

0.3-

15 СО 20

2.1

си

1.4 0.7 —I

15 СО 20

=1^(0,0.5,0) - амплитуда колебаний линии контакта на «верхней» пластине, С,ь = \С, (0, -0.5,0)| - на «нижней», = \С(0,0,0)| - в центре слоя и при г = 0.25 - =\С(0,0.25,0) ; значение краевого угла у на верхней пластине - у„, а на нижней -уь. Частоты собственных колебаний Отк, где т -азимутальное число, к - волновое число.

На рис. 2 показаны амплитуды колебаний боковой поверхности и краевых углов для рассматриваемых случаев неоднородности. Отметим, что мак-2.1

СИ

1.4

д

Рис. 4. Зависимость максимального отклонения: а - боковой поверхности Сл б - контактной линии Си,' в - контактной линии Сь; г - поверхности капли д, е - отклонения краевого угла у от частоты т для трех разных значений Хь (Ь=1.0, а=5.0, pi = 0.7, Хи = 10, Ль = 100). НПНП -пунктирная,ОПНП- штриховая, НПОП- сплошная линия

1.4-У

0.70-0.7-

~~I -0.7

I-1-1

0.7 х 1 -4

1.4-У

0.70-0.7-

~1—I—1—Г 1.4 -0.7 0

б

~~I-1-1

0.7 х 1 -4

1.4-У

0.70-0.7-

-1.4 -0.7

~I-1-1

0.7 х 1 -4

Рис. 5. Форма контактной линии 1+0ЛС,(гЛ) на верхней поверхности в разные моменты периода колебаний капли (Ъ = 1.0, а=5.0, р1 = 0.7, Хи = 10, Ль = 100, т = 1): а - НПОП; б - ОПНП; в -НПНП; t = 0 - пунктирная, t = 0.125Т- штриховая, t = 0.25Т - штрих - пунктирная, t = 0.375Т -сплошная линия

а

0

г

е

а

в

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

Параметр Хокинга описывает взаимодействие линии контакта с подложкой и, чем меньше его значение, тем сильнее взаимодействие. При Я = 0 - линия контакта неподвижна. Следовательно, при рассматриваемых значениях параметров Яи = 1 и Я = 0.1 взаимодействие линии контакта с подложкой велико и основное влияние оказывает именно неоднородность поверхности (рис. 2). В результате, в случаях НПОП (сплошная линия) и НПНП (пунктирная) появляются дополнительные резонансные пики по сравнению со случаем ОПНП (штриховая). Эти пики располагаются на частотах азимутальных мод, которые возбуждаются главным образом из-за неоднородности подложек. Неоднородность поля в данном случае не является основополагающей.

На рис. 3 показаны форма контактной линии капли на «верхней» пластине в разные моменты времени, а также значения краевого угла в течение периода для случаев НПОП и ОПНП. Формы линии контакта практически не отличаются (рис 3, а, б) в обоих случаях, и их указывают на наличие азимутальных мод в спектре вынужденных колебаний. Значения краевого угла отличаются незначительно (рис. 3, в).

С увеличением значения Я уменьшается сила взаимодействия линии контакта с подложкой, т.е. неоднородность электрического поля становится более значимой. С другой стороны, параметр Хо-кинга для рассматриваемой модели является эффективной амплитудой вибраций, что приводит к увеличению амплитуды колебаний. Таким образом, увеличение параметра Хокинга не вызывает в данном случае уменьшение влияния неоднородности подложки (рис. 4). В большинстве случаев видно наличие дополнительных пиков на частотах азимутальных мод. Это подтверждают и динамические формы линии контакта на рис. 5. Отметим, что в случае однородных пластин ОПНП (рис. 5, б) форма линии контакта наиболее близка к осесимметричной. Увеличение амплитуды в резонансе связано с уменьшением диссипации (см. рис. 2).

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

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

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

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

Исследование выполнено за счет гранта РФФИ № 20-31-90104.

Список литературы

1. Blake T. D., Clarke A., De Coninck J., De Ruijter M. J. Contact angle relaxation during droplet spreading: Comparison between molecular kinetic theory and molecular dynamics // Langmuir. 1997. V. 13 (7). P. 2164-2166.

2. Vafaei S., Podowski M.Z. Analysis of the relationship between liquid droplet size and contact angle // J. Colloid Interface Sci. 2005. V. 113 (2-3). P. 133-146.

3. Shklyaev S., Straube A. V. Linear oscillations of a hemispherical bubble on a solid substrate // Phys. Fluids. 2008. V. 20, 052102.

4. Fayzrakhmanova I. S., Straube A. V. Stick-slip dynamics of an oscillated sessile drop // Phys. Fluids. 2009. V. 21, 072104. DOI: 10.1063/1.3174446

5. Fayzrakhmanova I. S., Straube A. V., Shklyaev S. Bubble dynamics atop an oscillating substrate: Interplay of compressibility and contact angle hysteresis // Phys. Fluids. 2011. V. 23, 102105.

6. Erbil H. Y. The debate on the dependence of apparent contact angles on drop contact area or three-phase contact line: A review // Adv. Colloid Interface Sci. 2012. V. 170 (1-2). P. 67-86.

7. Иванцов А. О. Акустические колебания полусферической капли // Вестник Пермского университета. Серия: Физика. 2012. Вып. 3 (21). С. 16-23.

8. Snoeijer J. H., Andreotti B. Moving contact lines: scales, regimes, and dynamical transitions // Annu. Rev. Fluid Mech. 2013. V. 45. P. 269-292. DOI: 10.1146/annurev-fluid-011212-140734

9. Bormashenko E. Variational framework for defining contact angles: a general thermodynamic approach // J. Adhesion Sci. Tech. 2020. V. 34 (2). P. 219-230.

10. Ren J., Duan F. Recent progress in experiments for sessile droplet wetting on structured surfaces // Current Opinion in Colloid and Interface Science. 2021. V. 53, 101425.

11. Воинов О. В. Гидродинамика смачивания // Известия АН СССР. Механика жидкости и газа. 1976. № 5. С. 76-84.

12. Cox R. G. The dynamics of the spreading of liquids on a solid surface. Part 1. Viscous flow // J. Fluid Mech. 1986. V. 168. P. 169-194. DOI: 10.1017/S0022112086000332

13. Hocking L. M. The damping of capillary-gravity waves at a rigid boundary // J. Fluid Mech. 1987. V. 179. P. 253-266.

14. Hocking L. M. Waves produced by a vertically oscillating plate // J. Fluid Mech. 1987. V. 179. P. 267-281. DOI: 10.1017/S0022112087001526

15. Пухначев В. В., Семенова И. Б. Модельная задача о внезапном движении линии трехфазного контакта // Прикладная механика и техническая физика. 1999. Т. 40, № 4. С. 51-61

16. Картавых Н. Н., Шкляев С. В. О параметрическом резонансе полуцилиндрической капли на осциллирующей твердой подложке // Вестник Пермского университета. Серия: Физика. 2007. Вып. 1(6). С. 23-28.

17. Bonn D., Eggers J., Indekeu J., Meunier J. Wetting and spreading // Rev. Mod. Phys. 2009. V. 81. P. 739-805. DOI: 10.1103/RevModPhys.81.739

18.Xia Y., Steen P. H. Moving contact-line mobility measured // J. Fluid Mech. 2018. V. 841. P. 767783.

19. Amberg G. Detailed modelling of contact line motion in oscillatory wetting // npj Microgravity. 2022. V. 8, 1. DOI: 10.1038/s41526-021-00186-0

20. Lippmann G. Relations entre les phénomènes électriques et capillaries // Ann. Chim. Phys. 1875. V. 5. P. 494-549.

21. Berge B. Electrocapillarite et mouillage de films isolants par l'eau // Comptes Rendus Acad. Sci. II. 1993. V. 317. P. 157-163.

22. Quilliet C. and Berge B. Electrowetting: A recent outbreak // Curr. Opin. Colloid Interface Sci. 2001. V. 6. P. 34-39.

23.Mugele F., Baret J. C. Electrowetting: from basics to applications // J. Phys. Condens. Matter. 2005. V. 17. P. 705-774.

24. Mugele F. Fundamental challenges in elec-trowetting: from equilibrium shapes to contact angle saturation and drop dynamics // Soft Matter 2009. V. 5. P. 3377-3384.

25. Klarman D., Andelman D., Urbakh M. A Model of Electrowetting, Reversed Electrowetting, and Contact Angle Saturation // Langmuir 2011. V. 27. N. 10. P. 6031-6041.

26. Zhao Y. -P., Wang Y. Fundamentals and Applications of Electrowetting: A Critical Review // Rev. Adhesion Adhesives. 2013. V. 1. P. 114-174.

27. Royal M. W., Jokerst N. M., Fair R. B. Droplet-based sensing: optical microresonator sensors embedded in digital electrowetting microfluidics systems // IEEE Sensors Journal. 2013. V. 13. P. 4733-4742

28. Chen L., Bonaccurso E. Electrowetting - from statics to dynamics // Advances in Colloid and Interface Science 2014. V. 210. P. 2-12

29. Teng P., Tian D., Fu H., Wang S. Recent progress of electrowetting for droplet manipulation: from wetting to superwetting systems // Materials Chemistry Frontier. 2020. V. 4. no. 1. P. 140-154. DOI: 10.1039/c9qm00458k

30. Barman S. R., Khan I., Chatterjee S., etc. Electro-wetting-on-dielectric (EWOD): current perspectives and applications in ensuring food safety // Journal of Food and Drug Analysis. 2020. V. 28. N. 4. P. 595-621. DOI: 10.38212/2224-6614.1239

31. Alabuzhev A. A., Kashina M. A. The oscillations of cylindrical drop under the influence of a nonuniform alternating electric field // J. Phys.: Conf. Ser. 2016. V. 681, 012042.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

32. Alabuzhev А. А., Kashina М. А. Influence of surface properties on axisymmetric oscillations of an oblate drop in an AC electric field // Radiophysics and Quantum Electronics. 2019. V. 61. N. 8-9. P. 589-602

33. Алабужев А. А., Кашина М. А. Динамика зажатой капли в неоднородном электрическом поле // Вестник Пермского университета. Физика. 2019. № 4. С. 33-43. DOI: 10.17072/1994-35982019-4-33-43

34. Kashina М. А., Alabuzhev А. А. The forced oscillations of an oblate drop sandwiched between different inhomogeneous surfaces under AC vibrational force // Microgravity Science and Technology. 2021. V. 33. N. 3, 35. DOI: 10.1007/s12217-021-09886-4

35. Алабужев А. А., Любимов Д. В. Влияние динамики контактной линии на собственные колебания цилиндрической капли // Прикладная механика и техническая физика. 2007. Т. 48б. № 5. С. 78-86.

36. Alabuzhev A. A. Influence of heterogeneous plates on the axisymmetrical oscillations of a cylindrical drop // Microgravity Sci. Technol. 2018. V. 30, N. 1-2. P. 25-32.

37. Алабужев А. А. Осесимметричные колебания цилиндрической капли с подвижной контактной линией // Прикладная механика и техническая физика. 2016. Т. 57. № 6. С. 53-63

38. Алабужев А. А. О влиянии внутреннего давления на колебания цилиндрического газового пузырька // Вестник Пермского университета. Физика. 2020. № 4. С. 51-62. doi: 10.17072/1994-3598-2020-4-51-62

References

1. Blake T. D., Clarke A., De Coninck J., De Ruijter M. J. Contact angle relaxation during droplet spreading: Comparison between molecular kinetic theory and molecular dynamics. Langmuir, 1997, vol. 13 (7), pp. 2164-2166.

2. Vafaei S., Podowski M. Z. Analysis of the relationship between liquid droplet size and contact

64

A. A. Amôywee, M. A. nbHUKoea

angle. J. Colloid Interface Sci. 2005, vol. 113 (23), pp. 133-146.

3. Shklyaev S., Straube A. V. Linear oscillations of a hemispherical bubble on a solid substrate. Phys. Fluids. 2008, vol. 20, pp. 052102.

4. Fayzrakhmanova I. S., Straube A. V. Stick-slip dynamics of an oscillated sessile drop. Phys.Fluids. 2009, vol. 21, 072104. DOI: 10.1063/1.3174446

5. Fayzrakhmanova I. S., Straube A. V., Shklyaev S. Bubble dynamics atop an oscillating substrate: Interplay of compressibility and contact angle hysteresis. Phys. Fluids. 2011, vol. 23, 102105.

6. Erbil H. Y. The debate on the dependence of apparent contact angles on drop contact area or three-phase contact line: A review. Adv. Colloid Interface Sci. 2012, vol. 170 (1-2), pp. 67-86.

7. Ivantsov A. O. Acoustic oscillations of semispheri-cal drop. Bulletin of Perm University. Physics. 2012. no. 3 (21). pp. 16-23 (In Russian).

8. Snoeijer J.H., Andreotti B. Moving contact lines: scales, regimes, and dynamical transitions. Annu. Rev. Fluid Mech. 2013. Vol. 45, pp. 269-292. DOI: 10.1146/annurev-fluid-011212-140734

9. Bormashenko E. Variational framework for defining contact angles: a general thermodynamic approach. J. Adhesion Sci. Tech. 2020, vol. 34 (2), pp. 219-230.

10. Ren J., Duan F. Recent progress in experiments for sessile droplet wetting on structured surfaces. Current Opinion in Colloid and Interface Science. 2021, vol. 53, 101425.

11. Voinov O. V. Hydrodynamics of wetting. Fluid Dynamics. 1976, vol. 11, pp. 714-721.

12. Cox R. G. The dynamics of the spreading of liquids on a solid surface. Part 1. Viscous flow. J. Fluid Mech. 1986, vol. 168, pp. 169-194. DOI: 10.1017/S0022112086000332

13. Hocking L. M. The damping of capillary-gravity waves at a rigid boundary. J. Fluid Mech., 1987, vol. 179, pp. 253-266.

14. Hocking L. M. Waves produced by a vertically oscillating plate. J. Fluid Mech. 1987, vol. 179, pp. 267-281. DOI: 10.1017/S0022112087001526

15. Pukhnachev V. V., Semenova I. B. Model problem of instantaneous motion of a three-phase contact line. J. Appl. Mech. Tech. Phys. 1999, vol. 40. no. 4, pp. 594-603. DOI: 10.1007/BF02468433

16. Kartavih N. N., Shklyaev S, vol. On the parametric resonance of a semi-cylindrical drop on an oscillating solid substrate. Bulletin of Perm University. Series: Physics. 2007. no. 1(6), pp. 23-28. (In Russian).

17. Bonn D., Eggers J., Indekeu J., Meunier J. Wetting and spreading. Rev. Mod. Phys. 2009, vol. 81, pp. 739-805. DOI: 10.1103/RevModPhys.81.739

18. Xia Y., Steen P. H. Moving contact-line mobility measured. J. Fluid Mech., 2018, vol. 841, pp. 767783.

19. Amberg G. Detailed modelling of contact line motion in oscillatory wetting. npj Microgravity, 2022, vol. 8, 1. DOI: 10.1038/s41526-021-00186-0

20. Lippmann G. Relations entre les phénomènes électriques et capillaries. Ann. Chim. Phys., 1875, vol. 5, pp. 494-549.

21. Berge B. Electrocapillarite et mouillage de films isolants par l'eau. Comptes Rendus Acad. Sci. II. 1993, vol. 317, pp. 157-163.

22. Quilliet C. and Berge B. Electrowetting: A recent outbreak. Curr. Opin. Colloid Interface Sci., 2001, vol. 6, pp. 34-39.

23. Mugele F., Baret J. C. Electrowetting: from basics to applications. J. Phys.: Condens. Matter. 2005, vol. 17, pp. 705-774.

24. Mugele F. Fundamental challenges in electro-wetting: from equilibrium shapes to contact angle saturation and drop dynamics. Soft Matter 2009, vol. 5, pp. 3377-3384.

25. Klarman D., Andelman D., Urbakh M. A Model of electrowetting, reversed electrowetting, and contact angle saturation. Langmuir, 2011, vol. 27, no. 10, pp. 6031-6041.

26. Zhao Y. -P., Wang Y. Fundamentals and applications of electrowetting: a critical review. Rev. Adhesion Adhesives. 2013, vol. 1, pp. 114-174.

27. Royal M. W., Jokerst N. M., Fair R. B. Droplet-based sensing: optical microresonator sensors embedded in digital electrowetting microfluidics systems. IEEE Sensors Journal. 2013, vol. 13, pp. 4733-4742.

28. Chen L., Bonaccurso E. Electrowetting - from statics to dynamics. Advances in Colloid and Interface Science 2014, vol. 210, pp. 2-12.

29. Teng P., Tian D., Fu H., Wang S. Recent progress of electrowetting for droplet manipulation: From wetting to superwetting systems. Materials Chemistry Frontier. 2020, vol. 4. no. 1, pp. 140-154. DOI: 10.1039/c9qm00458k

30. Barman S. R., Khan I., Chatterjee S., et al. Elec-trowetting-on-dielectric (EWOD): Current perspectives and applications in ensuring food safety. Journal of Food and Drug Analysis. 2020, vol. 28. no. 4, pp. 595-621. DOI: 10.38212/22246614.1239

31. Alabuzhev A. A., Kashina M. A. The oscillations of cylindrical drop under the influence of a nonuniform alternating electric field. J. Phys. : Conf. Ser. 2016, vol. 681, 012042.

32. Alabuzhev A. A., Kashina M. A. Influence of surface properties on axisymmetric oscillations of an oblate drop in an AC electric field. Radiophysics and Quantum Electronics. 2019, vol. 61. no. 8-9, pp. 589-602.

33. Alabuzhev A. A., Kashina M. A. Dynamics of a drop under the influence of an inhomogeneous electric field. Bulletin of Perm University. Physics.

2019. no. 4. pp. 33-43. doi: 10.17072/1994-35982019-4-33-43.

34. Kashina M. A., Alabuzhev A. A. The forced oscillations of an oblate drop sandwiched between different inhomogeneous surfaces under AC vibrational force. Microgravity Science and Technology. 2021, vol. 33. no 3, 35. DOI: 10.1007/s12217-021-09886-4

35. Alabuzhev A. A., Lyubimov D. V. Effect of the contact-line dynamics on the natural oscillations of a cylindrical droplet. J. Appl. Mech. Tech. Phys. 2007, vol. 48, pp. 686-693. DOI: 10.1007/s10808-007-0088-6.

36. Alabuzhev A. A. Influence of heterogeneous plates on the axisymmetrical oscillations of a cylindrical drop. Microgravity Sci. Technol, 2018, vol. 30, no. 1-2, pp. 25-32.

37. Alabuzhev A. A. Axisymmetric oscillations of a cylindrical droplet with a moving contact line. Journal of Applied Mechanics and Technical Physics, 2016, vol. 57, no. 6, pp. 53-63.

38. Alabuzhev A. A. Effect of the internal pressure on oscillations of a cylindrical gas bubble. Bulletin of Perm University. Physics, 2020, no. 4, pp. 51-62. doi: 10.17072/1994-3598-2020-4-51-62 (In Russian).

Просьба ссылаться на эту статью в русскоязычных источниках следующим образом: Алабужев А. А., Пьянкова М. А. Влияние пространственной неоднородности подложек и электрического поля на динамику зажатой капли // Вестник Пермского университета. Физика. 2022. № 2. С. 56-65. doi: 10.17072/1994-3598-2022-2-56-65

Please cite this article in English as:

Alabuzhev A. A., Pyankova M. A. Influence of surface properties on the dynamics of a clamped drop in an in-homogeneous electric field. Bulletin of Perm University. Physics, 2022, no. 2, pp. 56-65. doi: 10.17072/19943598-2022-2-56-65

Сведения об авторах

1. Алексей Анатольевич Алабужев, канд. физ.-мат. наук, доц., с.н.с., Институт механики сплошных сред УрО РАН, ул. Академика Королева, 1, Пермь, 614013; профессор кафедры теоретической физики, Пермский государственный национальный исследовательский университет, ул. Букирева, 15, Пермь, 614990.

2. Марина Анатольевна Пьянкова, аспирант, инженер-исследователь, Институт механики сплошных сред УрО РАН, ул. Академика Королева, 1, Пермь, 614013.

Author information

1. Alexey A. Alabuzhev, Candidate of Physical and Mathematical Sciences, Senior Researcher, Institute of Continuous Media Mechanics UB RAS; 1, Akademika Koroleva st., Perm, 614013, Russia; Professor, Department of Theoretical Physics, Perm State University; 15, Bukireva st., Perm, 614990, Russia

2. Marina A. Pyankova, Postgraduate Student, Research Engineer, Institute of Continuous Media Mechanics UB RAS; 1, Akademika Koroleva st., Perm, 614013, Russia

i Надоели баннеры? Вы всегда можете отключить рекламу.