Scientific journal
Advances in current natural sciences
ISSN 1681-7494
"Перечень" ВАК
ИФ РИНЦ = 0,736

SURFACE GRAVITATIONAL ELECTROCAPILLARY WAVES

Taktarov N.G.
Исследовано распространение нелинейных поверхностных гравитационных электрокапиллярных волн на поверхности жидкого проводника. Библиогр. 6 назв.
Investigation of non-linear surface gravitational electrocapillary waves on surface of liquid conductor had been described.

Библиография работ, посвященных поверхностным волнам в жидких средах, взаимодействующих с электрическим полем, весьма обширна. В связи с этим здесь не представляется возможным привести сколько-нибудь подробный обзор этих работ. Впервые задача о распространении гравитационных электрокапиллярных волн на поверхности жидкого проводника (в линейном по амплитуде приближении) была решена, по-видимому, Я.И. Френкелем [3]. Авторы последующих работ по этой теме ограничивались линейным приближением. Лишь сравнительно недавно появились публикации, в которых при помощи методов возмущений [6] учитываются более высокие приближения. Так, например, в [2] в задаче об электрокапиллярных волнах на поверхности идеальной жидкости учитываются первые два приближения по амплитуде волны. Однако рассматриваемый нами анализ задачи с точностью до первых трех приближений позволяет более точно описать явление и объяснить новые эффекты.

Рассматривается распространение бегущих волн по заряженной поверхности бесконечно глубокого слоя жидкого проводника в поле тяжести. Несжимаемая и однородная жидкость граничит со средой пренебрежимо малой плотности (газ, вакуум).

Систему уравнений движения запишем в виде [3,4]

f; (1)

f, f (f ),

f, f,

где f- плотность, f - скорость,p* - давление, f - ускорение свободного падения, f- напряженность электрического поля, φ*- электрический потенциал, ε - диэлектрическая проницаемость. Звездочкой здесь и далее обозначены (в необходимых случаях) размерные величины, чтобы отличать их от безразмерных, обозначенных теми же буквами без звездочек. Отметим, что электрическое поле имеется лишь вне жидкого проводника.

Невозмущенная поверхность жидкости совпадает с плоскостью z*=0, жидкость находится в области f, ось z* направлена против f, ось x*- по направлению распространения волны.

Граничные условия на поверхности жидкости

f f       ,

f, f   (2)

f

где u - нормальная скорость поверхности,f - нормаль к поверхности, внешняя к области, занятой жидкостью, К - средняя кривизна поверхности, f- поверхностная плотность заряда, f- давление в атмосфере, γ- коэффициент поверхностного натяжения.

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

Уравнение поверхности жидкости запишем в виде f. На поверхности: f, в атмосфере f, где f - возмущение,f - невозмущенное поле. Аналогично: f, f.

Для установившихся бегущих волн предполагаем, что возмущения зависят от f, где с - фазовая скорость. В качестве малого параметра примем f, где f - максимальная амплитуда поверхности, λ- длина волны, предполагаемая заданной величиной.

Введем безразмерные величины

f  f, f,

f, f, f,(3)

f, f, f.

Уравнения (1) принимают вид

f, f,

f, f.

Граничные условия (2) должны быть записаны в безразмерных обозначениях. Необходимо добавить также условия периодичности и симметрии волны относительно вертикали, проходящей через вершину волны, а также условие расположения оси x* на плоской поверхности жидкости.

В результате получим нелинейную краевую задачу для нахождения величин (3): f, p, φ, Ē, σ, ξ, с. Для решения этой задачи применим метод малого параметра, использованный в [1] для исследования гидродинамических волн. Основная идея этого метода заключается в том, что многие дифференциальные уравнения, соответствующие конкретным физическим задачам, допускают введение безразмерного малого параметра δ, имеющего различный смысл в разных физических задачах, так что решение при f может быть найдено легко. Тогда решение при f ищется в виде ряда (не обязательно сходящегося) по степеням δ, такого, что нулевой член этого ряда соответствует решению краевой задачи при f.

Затем находим выражения для f, f, f, f, f, p*, f, фазовой скорости с.

Форма поверхности определяется соотношением

f      (4)

где A1, A2, A3 - некоторые коэффициенты, зависящие от двух безразмерных параметров взаимодействия, характеризующих относительную величину капиллярных и электрических сил по сравнению с гравитационными:

f, f     (5)

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

f, f,

где L - глубина частицы жидкости, c0 - фазовая скорость в линейном приближении (δ=0 ). Величина f известна как переносная скорость Стокса.

Наибольшего значения скорость f достигает на поверхности жидкости (L=0 ) . Если L→∞, то f. Наличие поверхностного натяжения приводит к увеличению переносной скорости. Переносная скорость приводит к разомкнутости траекторий частиц жидкости, которые наряду с колебательным движением обладают также равномерным движением в направлении распространения волны. С учетом выражения линейной фазовой скорости  запишем

f,

f, f.

Отсюда следует, что при определенных значениях параметров f и f величина f обращается в нуль. При f поверхностная волна переходит в гравитационную [1].

Для каждого значения F существует такая длина волны λ, для которой f= 0. И наоборот, для каждого значения λ существует такое значение F, при котором для этой волны f= 0.

При Kc>>1 (короткие волны), выражение f примет вид: f. Вводя периоды колебаний волны f и частицы f, можно записать: f Аналогично для частот: f. Видно, что период колебаний частицы превышает период волны.

Выражение для высоты волны (расстояние по вертикали от уровня впадины при f до уровня вершины при x=0) имеет вид

f(6)

При f выражение (6) совпадает с полученным в [1], при этом A2+A3= 0,5. В случае капиллярных волн (Kc >>1) можно пренебречь гравитационной и электрической силами, тогда

f.

Следовательно, при увеличении поверхностного напряжения γ, амплитуда волны уменьшается. При f величина A2+A3 ≈ -2,2.

Если  ff , а f, то

f.

Выражение для фазовой скорости (при f) имеет смысл только при f. При увеличении fвеличина A2+A3 возрастает, начиная от A2+A3 = 0,5 при Ke=0. Следовательно, высота волны возрастает при увеличении напряженности f, что хорошо согласуется с известными экспериментами [5]. При достаточно большом значении параметров f величина A2+A3 становится отрицательной, что указывает на уменьшение высоты волны по сравнению со случаем Kc=Ke=0. Если Kc>Ke, то при f амплитуда волны будет уменьшаться, при этом A2+A3 ≥ 6,5.

Выражения для амплитуды вершины

f

и впадины

f

находятся на основании (4). Очевидно, что f. Разность между этими амплитудами равна

f                 (7)

Из (7) следует, что f не зависит от электрического поля. При f разность f положительна, т.е. амплитуда вершины больше, чем впадины; при этом вершина уже (т.е. заострена), а впадина шире. При f рассматриваемая теория неприменима. Если f, то амплитуда у вершины меньше, чем у впадины.

Рассмотрим влияние электрического поля на фазовую скорость волны

f      (8)

где f - параметр, зависящий от Kc и Ke. Из (8) следует, что с зависит от квадрата амплитуды волны. При f выражение с принимает вид [1]: f. При f, f выражение f имеет вид

f.

Отсюда при Kc=0 следует f. Если f, то f.

Для Kc=0, Ke≠0 ( Ke<1) имеем:

f.

Отсюда видно, что при увеличении Ke величина f возрастает по модулю, оставаясь отрицательной. При Ke= 0 имеем f= -0,5.

Из (8) следует, что при возрастании f величина с при фиксированном k уменьшается и обращается в нуль при f. При заданных γ и f величина c0 достигает минимума при f:

f.

При f величина f обращается в нуль. Значению f соответствует величина поверхностной плотности заряда

f

При f поверхность становится неустойчивой - происходит ее разрушение. Различие между плотностью поверхностного заряда на вершине (x=0 ) и впадине (f ) равно

f, f,

где параметр В зависит только от Kc. Видно, что разность f не зависит от электрического поля. Для Kc=0 имеем

f

т.е. плотность заряда на вершине больше, чем во впадине. Это приводит к тому, что электрическая сила на вершине волны больше, чем во впадине; т.е. электрические силы создают неустойчивость поверхности, стремясь вытянуть вверх острия на вершинах волн. Поскольку f, разность f будет положительной при всех реальных значениях параметров в выражении Kc.

СПИСОК ЛИТЕРАТУРЫ

  1. Алешков Ю.З. Теория волн на поверхности тяжелой жидкости. - Л.: Изд-во ЛГУ, 1981. - 196 с.
  2. Белоножко Д.Ф., Климов А.В., Григорьев А.И. // Сб. докладов VII Международной научной конференции «Современные проблемы электрофизики и электрогидродинамики жидкостей». - СПб: Из-во СпбГУ, 2003. - 316 с.
  3. Ландау Л.Д., Лифшиц Е.М. Электродинамика сплошных сред. - М.: Наука, 1982. - 624 с.
  4. Ландау Л.Д., Лифшиц Е.М. Гидродинамика. - М.: Наука, 1986. - 736 с.
  5. Мелчер Дж. Р. // Магнитная гидродинамика. 1974. №2. С.3.
  6. Найфэ А.Х. Методы возмущений. - М.: Мир, 1976. - 456 с.