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

УДК 532.59, 537.29

С. О. Ю р ч е н к о

ВЛИЯНИЕ СЛАБОГО ЭЛЕКТРИЧЕСКОГО ПОЛЯ НА НЕЛИНЕЙНЫЕ ВОЗМУЩЕНИЯ ПОВЕРХНОСТИ РАЗДЕЛА КОНДЕНСИРОВАННЫХ СРЕД

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

E-mail: st.yurchenko@mail.ru

Ключевые слова: нелинейная теория, электрогидродинамика, уравнение Кор-тевега - де Фриза - Бенджамина - Оно, неустойчивость, конденсированные среды.

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

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

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

Рис. 1. К постановке задачи

(рис. 1); полупространство над жидкостью - вакуум. Жидкость предполагается несжимаемой, невязкой, идеально проводящей электрический ток, а движение жидкости - потенциально и в безразмерном виде описывается следующими уравнениями и граничными условиями:

z > eÇ:

eÇ > z > -1:

z —> œ:

z = -1:

z = eÇ: ,-2

Vzz + ^Ф = 0; Wzz + = 0;

Ф = -1, V^ = 0; V = 0; ф = 0;

eÇt - ¡r2vz + eV1ÇV1v = 0;

^ t+2 ((2+ivi42

Y

eju2 Ai£

E1

8npgh

( 2 + S

pgh (i\v4)

ViH2 )- A

pgh

где введены малые параметры задачи е = alh, /л = hlX. Здесь Ç - возмущение поверхности по вертикали от плоского уровня; h - толщина невозмущенного слоя; р - плотность жидкости; E - модуль напряженности невозмущенного электрического поля;

g - ускорение свободного падения; у, ф - потенциалы скорости и электрического поля; р0 - статическое давление на поверхности жидкости в отсутствие волнового возмущения; у - коэффициент поверхностного натяжения; У1, Д1 - двумерные операторы Гамильтона и Лапласа, действующие в плоскости ху; I, г - индексы, обозначающие частное дифференцирование по соответствующим переменным. Для приведения к безразмерному виду выбраны следующие единицы измерения:

[ ] = = [м] = ЕК

[£] = a, [г ] = Л, [x ] = [у ] = Л,

где а - амплитуда возмущения X - характерный горизонтальный

размер возмущения; с0 = ^[gh - скорость длинных гравитационных волн малой амплитуды.

Вывод квадратично-нелинейных по амплитуде уравнений, описывающих распространение поверхностных волн в исследуемой системе проводится при помощи разложения величин ф, у в ряд по степеням параметра е [9]. В приближении Буссинеска учитываются влияние квадратичной нелинейности и кубической дисперсии. Тогда электрическое поле влияет на дисперсионные компоненты, но не на слагаемые, не линейные по амплитуде. Эволюционное уравнение для волны, распространяющейся в положительном направлении [9], имеет вид

и+(3 и+Со) их+виххх— аЬ [их ]=0;

«=ß=c-0-pg 6

{ In,

3y 12ne

h2 +

2 2 Pg P g

Ь [и (х, I)] = | |к| | и (х', I)ехр (х - х'^^^х'ёк. (1)

—ад —ад

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

Вырождение эволюционного уравнения. Нелинейное эволюционное уравнение (1) может быть преобразовано к более удобному для анализа виду. Для этого стоит отметить связь между линейным оператором Ь[и] и оператором Гильберта Н[и], представленную в работе [10]:

^], Н [и ] = 1УР Г ^^, дх п J х — х'

L [u ] = -2nH

где интеграл понимается в смысле главного значения. Тогда уравнение эволюции примет вид уравнения Кортевега - де Фриза -Бенджамина - Оно:

+(с-+!"Ъ+ан[и~]=0 (2)

В приближении квадратичной дисперсии в = 0 уравнение (2) вырождается в уравнение Бенджамина - Оно [11]:

и + ( ео + 3 их + аН [ихх ] = (3)

Уравнение (3) впервые было выведено в приложении к эволюции волн на границе раздела двух глубоких жидкостей. Таким образом, прослеживается аналогия между влиянием поверхностного заряда и верхнего слоя жидкости в аналогичной системе (две несмешиваю-щиеся жидкости в поле сил тяжести). Физически такая аналогия связана, прежде всего, со спектрами линейных волн. В случае заряженной поверхности одной жидкости и двух несмешивающихся тяжелых жидкостей конечной толщины в дисперсионной зависимости ш(к) появляется квадратичное слагаемое по волновому числу к. В силу пространственной изотропии задачи, квадратичное слагаемое приходится писать в виде к|к|, откуда при обратном преобразовании Фурье появляется оператор Гильберта.

Уравнение Бенджамина - Оно (3) имеет точное решение в виде стационарной уединенной волны [11]. После перехода к бегущей переменной г = х - е^ + VI и однократного интегрирования уравнение (3) принимает вид

3

аН[и ] +—и2 + Vu + А = 0.

1 2 J 4

При условии А = 0 (для уединенной волны) и с учетом особенностей преобразования Гильберта находим решение в виде «рационального» солитона:

4а т, 2а

и =.+ 2 2, V = а (4)

4 + аа г 3а

где а - амплитуда нелинейной волны.

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

и

1/z2,

т. е. намного медленнее, чем в случае солитонного решения уравнения Кортевега - де Фриза [11], для которого справедлива асимптотика

u ~ exp(-z).

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

Существование солитона при отсутствии верхней жидкости обусловлено взаимодействием электрических и капиллярных эффектов. Солитонные решения, соответствующие полному кубическому спектру, являются электрокапиллярными, поэтому солитон вида (4) можно считать вырожденным электрокапиллярным. Значение параметра электрического поля а определяется, исходя из необходимого условия вырождения кубического спектра в = 0, а потому зависит от параметров, характеризующих систему:

Асимптотическое разложение Пуанкаре для периодических волн. Вернемся к нелинейному эволюционному уравнению (2). Обратимся к бегущей переменной г = х - (с0 + У)^, однократное интегрирование приводит к уравнению стационарных волн

При А Ф 0, а = 0 решения уравнения (6) - кноидальные волны, выражающиеся в эллиптических функциях Якоби. Вопрос о решениях уравнения (6) при произвольных а в настоящее время остается открытым, хотя уравнение типа Кортевега - де Фриза - Бенджамина - Оно уже было получено в работах [6, 9].

При малых а обнаруживается интересная внутренняя структура исходного уравнения (5). Так, в случае слабого электрического поля а << 1 и решение (6) ищут в виде асимпмтотического ряда:

ßu" + 3u2 - Vu + A + aH[u'] = 0.

(5)

В результате перехода к безразмерным переменным г ^ г^ЩУ, и ^ 2Уи, а^рУ ^ а, уравнение (5) приобретает вид

(6)

и (г) = Уо (г) + £а>п (г). (7)

п = 1

Подстановка решения вида (7) в уравнение (6) приводит к бесконечной системе неоднородных уравнений Шрёдингера:

< +2 ^2 -V о + А = 0;

Ж +(3у° -1) =-°п (г), п > 1;

01 (г) = Н[ц/0], Оп (г) = Н[уП-1 ] + ! "„-к, п > 2. (8)

к=1

Каждое уравнение Шрёдингера можно проинтегрировать, откуда для искомых функций уп(г) находим

Уп (г) = /((Ж0Оп (г)-2 йг. (9)

Периоды функций и 0п(г) являются кратными, откуда при повторном интегрировании возникают секулярные слагаемые. Все решения (9) становятся неустойчивыми при г ^ ю.

Солитонное решение получается, как видно, путем подстановки в (8) и (9) невозмущенного профиля ~ 8ееЬ2(г/2). При этом решение (9) сходится: при больших г ведет себя как ехр[-г], а интеграл расходится ~г. Таким образом, при г ^ ю возмущение солитонного решения ведет себя ~гехр[-г], а множители слагаемых ~ап в случае солитонного решения - сходятся как гп ехр[-г]. В случае периодического решения при г ^ ю интегралы становятся осциллирующими, а по модулю расходятся ~г.

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

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

ад ад

г = г'Х аЧя, и (г) = у ° (г ) + £а>я (г). (10)

п = 0 п = 1

Подстановка решения (10) в (6) приводит к уравнениям

ш0 + -ш2 -Шо + А = 0, ш -Щ-ш„ = -Оп (*), (11)

2 Шо

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

О (г) = Н[ш0]-2® Щ,

Шо

ш'" 3

О- (г) = - ( +2®-) ш0 - 2®1 Щ-ш1 +- Ш? + Н [ш0 + ш1 ].

Шо 2

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

Невозмущенное уравнение Кортевега - де Фриза можно представить [12] как

u0 = — dn s

f ^ Л , ь - ь.

-Г1П

2

2a , 2

-z, s

V 2 s у

+ К a = s =

- Ь2

ь - Ы'

где bj, Ь2, Ь3 - действительные корни кубического уравнения

u03 - u02 + Au0 + B = 0

в порядке убывания, 5 - модуль эллиптической функции Якоби dn(z, 5), а - амплитуда нелинейной волны. Для расчетов оказывается удобным использовать известное тригонометрическое разложение эллиптической функции Якоби [13]:

п 2п

^ nnz Л

, / \ /L ¿./L -г—1 I/ /LII*.

dn(z,s) =-— + —— > -— cos ——

V ' ' 2K (s) K (s) ^=11 + q2n УK (s)

q = exp(-nK(s')/K(s)), s' = VT-

s2.

(12)

где К(5) - полный эллиптический интеграл первого рода; 5' - сопряженный модуль. Модуль 5, как известно, может служить мерой нелинейности волны: предел 5 << 1 соответствует тригонометрическим функциям, а 5 ^ 1 - гиперболическим функциям. Использование разложения (12) позволяет находить преобразование Гильберта Н[м0'] с любой наперед заданной точностью в виде тригонометрического ряда [12].

Для кноидальной волны справедливо приближенное разложение

dn [z, ^] « 1 - 8q sin2

nz

, q <<1;

2К _

(г) = аёп2 [г,s]. (13)

Решение представляется в следующем виде:

У1 (г) = у0 (г )<Ь + А ) ) dz, (14)

где константа А определяется из условия = 0:

I 2 — К

А = -Иш ,

г +1

(г) = < (2 )dz ))0 )-2 w2 (г) = <1« )-2 Вычисления с учетом (13) приводят к следующим результатам:

= 11 =

= -B

пт

K

-2q(1 -4q)l 3-cos

(1 -8д-16д2 )со8

г = + а®1 + О (а2 ))■ Найденные константы равны:

2пт K

(15)

K п(1 - 16q + 96q -256q + 352q4)

0 =—, A = —--.

1 8п 3K

B =

Ka (1 - 4q) 24nq '

Функция у1(т) для кноидальных волн с различными модулями ^ представлена на рис. 2.

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

Пространственный эффект влияния слабого электрического поля состоит в «расплывании» возмущений по горизонтали и связан с тем, что > 0. Длина кноидальной волны в слабом электрическом поле

, ч [20 f aK(s

Х = 2sK (s). -^ 1 +-.

\aV , 8nJßV ,

8n4ßV,

íiM 1 ^^0,05 Oy 0,5 \Vo,7 \Г\\/0,8 \ 0,85

\ -i ц // ^ЛУЛ I h Т

' -1

Рис. 2. Вид функций у1(г) при различных значениях модуля ж (цифры у кривых)

Возмущенное решение определяется выражением u (г) = |0 (т) + а\| (т) + 0 (а2) « « a(dn2 [т, ^] + sv/1 (т)) = au(т), aK (1 - 4q)

s=

24nqJßV '

(16)

Выражение (16) для й(т) справедливо при в > 0, если в < 0, из исходного уравнения (5) получаем

?(т) = -dn2 [т,^] + s\j/1 (т).

(17)

На рис. 3-6 приведены профили нелинейных волн с различными модулями и возмущениями, рассчитанными согласно (16) и (17). Видно, что электрическое поле приводит к заострению профиля нелинейной волны в случае в > 0. Нелинейная волна с модулем ^ = 0,8 (см. рис. 6) подвергается такому влиянию несколько меньше, чем волна с модулем ^ = 0,5 (см. рис. 3).

В случае в < 0 слабое электрическое поле, наоборот, оказывает сглаживающее воздействие на нелинейные волны. Так, на рис. 5 отчетливо видно, что при ^ = 0,5 с увеличением параметра е меняется амплитуда нелинейной волны, после чего на месте минимумов

fi(T)

Рис. 3. Профили волн с модулем ж = 0,5 представлены кривыми (1-3) с соответствующими значениями параметра £ = 0; 0,05; 0,1 (в > 0)

"Ml

-1,5 -1 -0,5 0,5 1 1,5 Т

Рис. 4. Профили волн с модулем ж = 0,85 представлены кривыми (1-3) с соответствующими значениями параметра £ = 0; 0,05; 0,1 (в > 0)

й(т) \ -0,8 1 / \ / 2

vNT ' -^У! 4

-1,5 -1 -0,5 0,5 1 1,5 т

Рис. 5. Профили волн с модулем ж = 0,5 представлены кривыми (1-4) с соответствующими значениями параметра г = 0; 0,07; 0,1; 0,15 (р < 0)

в(т)

-1,5 -1 -0,5 0,5 1 1,5 Т

Рис. 6. Профили волн с модулем ж = 0,8 представлены кривыми (1-3) с соответствующими значениями параметра г = 0; 0,07; 0,15 (р < 0)

исходного профиля может появиться максимум возмущенного профиля. Нелинейные волны с модулем s = 0,8 менее подвержены влиянию электрического поля, но характер влияния поля остается прежним - профиль волны сглаживается.

Здесь еще следует отметить, что кривая на рис. 2, отвечающая возмущению волны с модулем s = 0,85, напоминает форму возмущения, найденную для солитонного решения [9], которому отвечает предел s ^ 1. Решение (9) совпадает с возмущенным решением для солитона Кортевега - де Фриза в указанном предельном случае, так как является более общим.

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

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

Работа выполнена при поддержке РФФИ (проекты № 12-0831104 мол а, 12-08-33112 молавед).

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

1. Ж а к и н А. И. Нелинейные волны на поверхности заряженной жидкости. Неустойчивость, ветвление и нелинейные равновесные формы заряженной поверхности II Изв. АН СССР. МЖГ. - 1984. - № 3. - C. 94-102.

2. E a s w a r a n C. Solitary waves on a conducting fluid layer II Phys. Fluids. - 1988. -Vol. 31. - № 11. - P. 3442, 3443.

3. G o n z a l e s A., C a s t e l l a n o s A. Kortweg-de Vries-Burgers equation for surface waves in nonideal conducting liquids II Phys. Rev. E. - 1994. -Vol. 49. - № 4. - P. 2935-2940.

4. G o n z a l e s A., C a s t e l l a n o s A. Nonlinear electrohydrodynamic waves on films falling down an inclined plane II Phys. Rev. E. - 1996. - Vol. 53. - № 4. -P. 3573-3578.

5. T i l l e y B., P e t r o p o u l o s P., P a p a g e o r g i o u D. Dynamics and rupture of planar electrified liquid sheets II Phys. Fluids. - 2001. - Vol. 13. - № 12. -P. 3547-3563.

6. A new aplication of the Korteweg-de Vries Benjamin-Ono equation in interfacial electrohydrodynamics I H. Gleeson, P. Hammerton, D. Papageorgiou et al. II Phys. Fluids. - 2007. - Vol. 19. - № 3. - 031703. - P. 3.

7. Effect of an electric field on film flow down a corrugated wall at zero Reynolds number / D. Tseluiko, M. Blyth, D. Papageorgiou et al. // Phys. of Fluids. - 2008. -Vol. 20. - № 4. - 042103, 19 p.

8. E l - S a y e d M. Nonlinear analysis and solitary waves for two superposed streaming electrified fluids of uniform depths with rigid boundaries // Arch. Appl. Mech. -2008. - № 78. - Р. 663-685.

9. А л и е в И. Н., Ю р ч е н к о С. О. О нелинейных волнах, распространяющихся на поверхности идеальной проводящей жидкости в электрическом поле // Известия РАН. Механика жидкости и газа. - 2009. - № 5. - С. 139-150.

10. Б е й т м е н Г., Э р д е й и А. Таблицы интегральных преобразований: В 3-х т. Преобразования Бесселя. Интегралы от специальных функций. - М.: Наука. -1970. - Т. 2. - 328 с.

11. А б л о в и ц М., С и г у р Х. Солитоны и метод обратной задачи: пер. с англ. -М.: Мир. - 1987. - 479 с.

12. К а р п м а н В. И. Нелинейные волны в диспергирующих средах. - М.: Наука. -1973. - 175 с.

13. Я н к е Е., Э м д е Ф., Л ё ш Э. Специальные функции (формулы, графики, таблицы). - М.: Наука. - 1964. - 344 с.

Статья поступила в редакцию 30.05.2012