Научная статья на тему 'К расчету напряженности электростатического поля осциллирующей незаряженной капли'

К расчету напряженности электростатического поля осциллирующей незаряженной капли Текст научной статьи по специальности «Физика»

CC BY
149
22
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
сфероидальная капля / однородное электростатическое поле / осцилляции / волновое и квазистационарное приближения / spheroidal drop / uniform electrostatic field / oscillation / wave and quasistationary approximations

Аннотация научной статьи по физике, автор научной работы — С. О. Ширяева, А. И. Григорьев, Н. Ю. Колбнева

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

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

The carried-out calculations afforded an analytical expression, in the linear approach on indignation amplitude, for the intensity of an electric field of an oscillating uncharged drop of a nonviscous conductive liquid, placed in an electrostatic field. Estimates of the field intensity values, both in the wave approach (at long distances from a drop) and in the quasistationary approach (at a drop surface) are carried out.

Текст научной работы на тему «К расчету напряженности электростатического поля осциллирующей незаряженной капли»

К расчету напряженности электростатического поля осциллирующей незаряженной капли

С. О. Ширяева, А. И. Григорьев, Н. Ю. Колбнева

Ярославский государственный университет им. П.Г. Демидова, ул. Советская, 14, г. Ярославль, 150000, Россия, e-mail: shir@uniyar.ac.ru,

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

Ключевые слова: сфероидальная капля, однородное электростатическое поле, осцилляции, волновое и квазистационарное приближения.

УДК 532.62: 541.24

ВВЕДЕНИЕ

Незаряженная сферическая капля радиуса Я с коэффициентом поверхностного натяжения с, помещенная в электростатическое поле напря-

женностью

поляризуется в нем и с точки

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

3E0 R

e =

W^ V О

В жидкости всегда существует капиллярное волновое движение, поэтому имеет место искажение равновесной поверхности раздела сред весьма малой (тепловой) амплитуды с характерной высотой гребней £ ~ VкТ / о , где к - постоянная Больцмана; Т - абсолютная температура [4]. Эти волны порождаются тепловым движением молекул жидкости. Волновое искажение поверхности капли сказывается и на создаваемом электрическом поле капли, воздействуя на ди-польную структуру поля и приводя к появлению зависимости напряженности поля от времени.

Расчет напряженности электрического поля такой капли представляет интерес в связи с исследованиями ее устойчивости во внешнем электростатическом поле [5-7] и ее нелинейных ос-цилляций [8]. Эта проблема вызывает интерес и при исследовании возможности зажигания коронного разряда в окрестности поляризованной капли [9] при напряженностях полного поля,

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

ПОСТАНОВКА ЗАДАЧИ

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

Будем решать задачу в сферической системе координат с началом в центре масс капли. Уравнение возмущенной волновым движением поверхности капли можно записать в виде:

г ( 0, t ) = г (0 ) + 4 (0, t),

где г(9) - равновесная форма капли, которую согласно сказанному выше будем считать сфероидальной; ^(9, t) - капиллярное волновое возмущение поверхности капли: max << min г (0).

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

Математическая формулировка задачи имеет вид:

dV. (Г, t) 1 - , ч - , ч

А ' = -—■VP. (Г,t); divV. (г,t) = 0; j = 1;2; dt р j

, 1 82Ё (Г, t) .

АЕ (Г, t ) = —-z^2; divE (Г, t ) = 0;

c et

r ^ 0: V1( r, t 0;

© Ширяева С.О., Григорьев А.И., Колбнева Н.Ю., Электронная обработка материалов, 2016, 52(2), 37-44.

г у2 (г, г)^ 0; Ё (г, г)^ Ё0 - Ё0е2;

г = г (0) + ^: ^ + (, = 0; ^ (г, 0, г)-г - г (0)-£ (0, г);

р - Р2 + Рё - Р. = 0; (г, Ё (г, г)) = 0.

Здесь ^ (г, г) - скорость движения жидкости в капле (/ = 1) и среде (/ = 2). Аналогично Р. (г,г) -гидродинамическое давление в капле (/ = 1) и среде (/ = 2); Ё (г, г) - напряженность электрического поля; т - орт касательной к возмущённой поверхности капли; Рё = еХЁ2 /8п - давление электрического поля; Р. = о п - капиллярное

давление. Орт нормали к поверхности капли п определяется выражением:

п =1

rotrot E +

_L aE

c2 at2 = 0;

div E = 0.

Представим напряженность электрического поля в виде:

Ё = Ё( 0)+ Ё(1), (1)

г (0)

где Ёу - напряженности поля нулевого порядка

малости по е2 и Ё(1) - добавка к напряженности первого порядка по е и

Для отыскания напряженности поля в первом

приближении Ё(1) следует решить систему уравнений неразрывности электрического поля и волнового:

1 а2 Е(1)

divE (1)= 0; aE (1)=-

c2 at2

(2)

с граничными условиями:

VF

Для замыкания выписанной системы уравнений введем условия неизменности полного объема (следствие несжимаемости жидкости), неподвижности центра масс, а также незаряженности капли соответственно: 4

r ^да: e(1) ^ 0

JdV = 3nR3; J rdV = 0;

V 3 V

V = [0 < r < r (0) + 4 (0, t); 0 < 0 < п;0 < ф < 2n];

cJ(n, E)s = 0,

S

S = [r = r (0) + 4(0, t), 0 < 0 < п,0 < ф < 2n].

Решения этой задачи будем искать, полагая эксцентриситет капли малым, в виде разложений по полиномам Лежандра:

да

4 (0, t ) = Х аи (t ();

n=2

an () = an exP (-iaj), где con - частота капиллярных колебаний поверхности, в общем случае комплексная; an - амплитуда, Рп(ц) - полином Лежандра, ц = cos 0.

Нижеследующие расчеты будем проводить с точностью до произведения двух малых параметров, включая задачи: квадрата эксцентриситета и малого возмущения (е2 •

НАПРЯЖЕННОСТЬ ПОЛЯ У ПОВЕРХНОСТИ ОСЦИЛЛИРУЮЩЕЙ ВО ВНЕШНЕМ ОДНОРОДНОМ ЭЛЕКТРОСТАТИЧЕСКОМ ПОЛЕ НЕЗАРЯЖЕННОЙ КАПЛИ

Поле вне капли удовлетворяет уравнениям:

г = г (0): (Т, Ё (1)(г, г )) = 0; ф( П, Ё(1)) ёБ = 0.

Представим вектор Ё(г, г) в виде разложения на три ортогональных вектора [10]:

Ё(1) (г, г) = #1<х>1 (г, г) + N Ф 2 (г, г) + ^Ф 3 (г, г), (3) где Ф/ (г, г) - неизвестные скалярные функции,

а векторные операторы-проекторы N■ имеют вид:

N -V, N - Ухг, N -Ух(Ухг),

^+--У, Щ- г XV, N3 -(г х у)х V и удовлетворяют условиям ортогональности:

N + • N = 0, при т Ф /,

/ т 5 1 и у

где символ "+" означает эрмитово сопряжение; а жирная точка между векторными величинами здесь и далее - скалярное умножение.

С учетом приведенного свойства ортогональности операторов первое из уравнений системы (2) преобразуется в уравнение Лапласа для функции ф1 (г, г):

VE(1) =-т E=

(1) _

= -//+•( + ^Ф2 + ^Ф3 ) =

= -N+ • Nf101 = V • V01 = АФ1 = 0.

(4)

Несложно убедиться, что операторы N■ коммутируют с оператором Лапласа, то есть

F =0

ЛЛА = АЛТ, благодаря чему второе уравнение

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

свойством коммутативности, приведем его к виду:

X N, 1=1

аф 1 (г, г )-

1 д2Ф 1 (Г, г)

с2 дг2

= 0.

. , 1 д 2ф . (г, г)

аф 1 (г, г)-—-= 0, 1 = 1,2,3.

дг2

ты

= о • ё. +

1 дФ 2 (г, г) дФ 2 (г, г)

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

--— • ё0---—- • ёф,

г 0 дф д9

1

Л3Ф3 (г, г ) = --Ь Ф3 (г, г )ёг +

1 д дФ3 (г, г)

+--г—^—'- е0 +

г дг д0

1 д дФ3 (г, г)

--г--——е

г 8т 0 дг дф ф'

где ь =

1

к2Ф 2 (г, г ) = --

Л3Ф3 (г, г ) = - - ¿Ф3 (г, г)

г

1 д дФ3 (г, г) +--г-—— е0.

г дг д0

е +

Функции Ф2 (г, г) и Ф3 (г, г) являются решениями уравнений Гельмгольца, имеющих вид:

АФ ? (г, г )-2 Ф, (, г ) = о, д = 2,3. (9)

Решение уравнений (9) для функций Ф (г,г)запишется в виде:

Фд =Х)ехр(-тяг)Ря (ц), д = 2,3, (10)

п

где кп2 (кг) - сферическая функция Бесселя третьего рода, выбираемая из тех соображений, чтобы на асимптотике (г ^ да) электромагнитная волна расходилась; к = (а>п / с)- волновое число.

Неизвестные константы в решениях (10)

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

Ф 2 (г, г), эти граничные условия приведут к следующим соотношениям:

г (0, г ) = г (0) + £ (0, г): Л2Ф2 (г, г )• т = о, (11) Л2Ф2 (г,г)• п • йХ = 0, (12)

где Т и п - единичные векторы касательной и нормали к поверхности капли.

Поскольку величина Е (1)(г, г), а, следовательно, согласно (3), и функции Ф. (г, г) имеют порядок малости ~ е£, (или ~ е2^), то выражения (11) и (12) для слабых полей можно использовать, исходя из сферической поверхности капли, то есть при г = Я. При этом вектором Т могут служить орты сферической системы координат е0 и е , а

вектор нормали п будет совпадать с ортом ег (в

используемом приближении слагаемыми ~ е2£, пренебрегаем).

Согласно (9) тороидальная компонента поля

(Л^Ф2 (г,г)) имеет лишь составляющую ~ еф, и,

значит, соотношения (11) и (12) при Т = е0 удовлетворяются тождественно при любых константах в решении (10). Исходя из условия (11) при Т = е получим

Умножая слева полученное выражение последовательно на и учитывая, что • ф 0, получим три скалярных уравнения:

(5)

Поскольку решается задача об излучении каплей электромагнитных волн, естественно принять Ф ~ ехр (-шпг), при этом уравнения (5)

сведутся к уравнениям Гельмгольца.

Из уравнения (5) при значении индекса 1 = 1 и уравнения (4) получим (2 / п2)Ф1 (г,г) = 0, и поскольку частота колебаний отлична от нуля, то, следовательно, скалярная функция Ф1 (г,г) = 0.

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

Е(1) (г, г ) = Л2Ф 2 (г, г) + Л3Ф3 (г, г).

В сферической системе координат компонен-Е(г, г) имеют вид:

л2ф 2 (г, г ) =

(6)

(7)

--1 0— I - угловая часть операто-

¡ип 0 5050 )

ра Лапласа в сферических координатах.

В силу того что рассматривается осесиммет-ричная задача, в выражениях (6) и (7) следует отбросить слагаемые, содержащие производные по ф, тогда получим:

дФ 2 (, г) е ---——е

д0 ф'

2

г = Я:

дФ2 (г, г)

д0

= 0.

Поскольку это соотношение должно быть справедливым при любом значении угла 0, то необходимо потребовать обращения в нуль всех констант в решении для ф2 (г, г). Таким образом, поле Ё(1)(г, г) полностью определяется скалярной функцией Ф3:

Ё(1) (г, г) = Ж,Ф 3 (г, г) - V х (V х г) Ф3 (г, г)

или с учетом (10) получим:

Ё (1)(г, г) =

= X ехр(-/ю„г) Ы3)г-1кП2)(кг )п(п + 1)Ри (ц)еГ +

Б(3)г) + ^^^

дг

р(ц) , I (13)

д0

п =

1 -

1 дг (0) д4 (0, г)

г2 д0 д0

е., -

(15)

1 _д_

г д0

(г (0 ) + 4 (0, г )).

В силу осевой симметрии задачи орт касательной к параллелям совпадает с соответствующим ортом сферической системы координат тф = еф.

Орт касательной в меридиональном направлении найдем, воспользовавшись векторным соотношением

Т0 = п х Тф .

С учетом (15) выражение для орта касательной запишется в виде:

Г0 = 1 (г (0) + 4 (0, г ))еег +

г д0 (16) \ - 1. дтО д4 (0, г)4

г2 д0 д0

е0.

Подставляя (16) и (1) в условие эквипотенциаль-ности проводящей капли и раскладывая полученное выражение в окрестностях сферы, перепишем условие эквипотенциальности капли в виде:

ё<"

.дЕШ. 4 (0, г)+1 Ё«

дг ^ ' г д0 г

= 0, (17)

где Ё^г0) и Ё00) - соответствующие компоненты вектора напряженности поля нулевого порядка малости по е2 и

Ё(0) = Ё0 С08 (

(

1 + 2 Яг

г

V

3 \ (

е„ - Ё0

3

1 - Яг

V У

0. (18)

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

Подстановка выражений для Ё^ из (13)

Ёг0) и из (18) и возмущения поверхности

капли 4(0, г) в виде ряда по полиномам Лежандра в выражение (17) позволяют определить выра-

г>(3)

жение для константы Бп ' через амплитуды ап:

г(0)

Сферическую функцию Бесселя третьего рода при целых значениях индекса выразим через элементарные функции [11]:

^00 = ^ Ь-+1 (2п-т)! 1 . (14)

п \ г / 1 / \ | ■ / ч п-т

г т~0 (п -т)!т! (2г) Для определения неизвестной константы Б^3 в (13) учтем, что поверхность проводящей капли эквипотенциальна, то есть проекция вектора на орт касательной к поверхности т есть нуль:

(т, Ё (г, г )) = 0.

Чтобы найти орт т , сначала запишем орт нормали к возмущенной поверхности капли:

Б(3) = 3Ё0 У 0 У1У 2а п-Х-1 + У 0а„+1Пп+1 .

дг (2)(кг))

(19)

П+-

п +1

^(2п + 1)(2п + 3) ^(2п - 1)(2и +1) Подставляя (19) в (13), получим выражение для напряженности электрического поля Ё :

Ё «(?, г) = - ^ X Мп (г)

г п=0

К2 К. кг)

хп(п +1)рп (ц) +

д г ( гк^ \кг) )

дг (г42)(кг)) дРп (ц)

д г (гк(\кг)

д0

(20)

мп (г) = У0У1У2ап-1 ехр(/®п-1г)Ц+-1 + У0 а п+1 ехР(г'®п+1г)Ц -+1;

+ п +1 - п

Ц п

(2п +1)

Цп =

(2п +1)

Электрическое поле, создаваемое каплей, колеблющейся во внешнем электрическом поле Ё0, при выполнении условия Я/г << 1 принимает вид волнового поля при кг = га>п/с >> 1 и квазистатического поля при кг = га>п/с << 1, где г - расстояние от капли до точки наблюдения; к = — - вол-

с

новое число.

Сферические функции Бесселя третьего рода (сферические функции Ханкеля) при целых зна-

п

п

чениях индекса выражаются через элементарные функции [11]:

к?^ )=I ехр („ );т,»+1 (^т—)

- =0 (п -т)!:! (22)п

:=0 2 = кг.

£%,г) = ^Xмп (г)| ^ I ехр[-/(2 --0)]

3Е,

п=0

1 (2п - т)! (22): 2 :=0 (п - : )!:!

X (((и:) (2/20): (.+(п -:)

:=0

(п - т )!:!

п (п + 1)Рп (ц )

X (2п - :)! (Ъ2): (. + (п - :) 1

:=0 ( -:)!:! 1 - ) дРп (ц) . XX (2л-) (2.-0):+ _(п-:)] д0 '

:=0

(п - т )!:!

Е '1|(гг,' > - £ "(' ^)!

-ехр [-/ (2 - 20 )]>

XX М(2.2 (п + 1)Ря (ц)

' :=0

(п - :)!

:=0

(п - : )!

д0

3Е0

п=2

Е(1)(г, г) X ап ехр (Ы)

г п

ехР (2 - 20 )]

п!

(2 ( +1)

(Г^ <2Й«"+')(п+2>рп+1 (Ц )ег-

-( + ( +1 - :)) дР"+1 (Ц), V '' д0

(2( - 1))! !=0 ( - 1 - !)!т! {п(п - 1)Р_1 (ц) - ( 2 + (п -1 - т))дPn50<ЦТе

(21)

Подставляя (21) и выражение для производной дг (гк(] (кг))

д г ( к 2) (кг ) )=д 2 ( 2к^( 2) ) =

(2п - т)! 1 (. (п - тУ

= - ехр (-.2) X / т+1) Г —| / + '

т=0 (п -т)!т! (22)п-т I 2

в (20), получим напряженность электрического поля первого порядка малости:

Представляет интерес рассмотрение асимптотики этого выражения на больших и малых расстояниях от капли.

Рассмотрим случай волнового поля, то есть поля излучения на больших расстояниях от излучающей капли. При 2 >> 1 (кг = та>п/с >> 1), сохраняя в выражении (24) в суммах по т старшие слагаемые (с наибольшей степенью 2), для напряженности электрического поля запишем:

Е^аге (г, г) - — X ап ехр (¡Ы) I —

3Е,

п=2

ехр [-/ (2 - 20 )]

п!

(2 (п +1))!

<(2/2)п+1 ](п + 1)(п + 2)Рп+! (ц)г - ¡2^^^ей К

(22)

_ (п - 2)! , чп_1 | , ч , ч

+цп (2 (п - 1))! (2/2) |п (п - 1)Рп4 (ц)ег - '

д0

др- (ц).

де

где 20 = кЯ. Учитывая, что 20 << 1, и сохраняя в суммах, стоящих в знаменателях, лишь старшее слагаемое (т = 0), выражение (22) представим в виде:

(п -1)! ~п+1

" X мп ( \

г

Несложно заметить, что в рассматриваемом случае радиальной составляющей поля можно пренебречь по сравнению с его угловой и, учитывая, что 20 << 2, записать для напряженности поля в волновом приближении следующее выражение:

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

Ё^ае (г,г) - -^ X ап ехр[/ (Ы - 2)](2/ 20 )"

2г п=2

дРп+1 (ц )

-XX(т^-ГI/+«п-т^дР^1(23)

+ ц п

\2 + п__

( 0 )цп (2 (п +1))! д0

(п - 2)! дР -1 (ц)" (2 (п -1))! д0

(25)

Подставим в (23) явный вид коэффициентов Мп(г) и перегруппируем слагаемые в суммах по индексу п:

Рассмотрим случай квазистатического поля или поля в ближайшей окрестности капли. При 2 << 1 (кг = га>п/с << 1), сохраняя в выражении (24) в суммах по т слагаемые с наименьшей степенью 2 (т = 0) и принимая во внимание, что 20 << 2, для напряженности электрического поля запишем:

Еа(г,г) X а" ехр[/Ы - 2)]р

3Е,

цп (£0] + 2)( + 1)Рп+1 (ц)ёг-(¡2 + (п + l))ГPn±-Тц)ё0 ^ +

(п +1)

д0

(п -1)

'(п - 1)Pn_1 (ц)ёг - (¡'2 + (п - 1)

дРп-1 (ц) -

д0

х

Рис. 1. Зависимость амплитуды напряженности волнового поля (единицы измерения - 10-22 СОБЕ) от угла наблюдения и расстояния (единицы измерения - см), рассчитанная при п = 2, ст = 73 дин/см, р1 = 1 г/см3, р2 = 1,3-10"3 г/см3, ап = 0,1Я, Я = 30-10"4 см, Е0 = 0,1 СОБЕ.

Рис. 2. Зависимость амплитуды напряженности волнового поля от угла наблюдения и напряженности внешнего электрического поля Е0 (единицы измерения - СОБЕ): (а) - рассчитанная при тех же прочих значениях физических величин, что и на рис. 1, и г = 5000 см (единицы измерения - 10-22 СОБЕ); (б) - рассчитанная при тех же прочих значениях физических величин, что и на рис. 1 (единицы измерения - 10-23 СОБЕ).

Рис. 3. Зависимость амплитуды напряженности квазиста- Рис. 4. Зависимость амплитуды напряженности квазистати-

тического поля (единицы измерения - 10-11 СОБЕ) от угла ческого поля (единицы измерения - 10-11 СОБЕ) от угла

наблюдения и расстояния (единицы измерения - см), рас- наблюдения и напряженности внешнего электрического

считанная при тех же прочих значениях физических вели- поля Е0, рассчитанная при тех же прочих значениях физиче-

чин, что и на рис. 1. ских величин, что и на рис. 1, и г = 5 см.

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

^stat(r,t) -

X an exp (rnj) I — I exp [-/ (z - zo ) > ' n=2 V z )

[z0)2 Ц+ {(n + 2)Pn+l (ц) e, j +

+ -{ P ( )e P-t ( Ц) e Ц (26)

+ Цп jnPn-t ( Ц)---seee

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

Eüt (?, t) - ^ ехР [-'' (z - zo ))

ад

<E ((-tan-t exp (i®n-tt) + Ц-+1Нп+1 exp (n+tt))

n=2

n+t

(n+t)Pn (ц)

9Pn (ц) -

se

E[wie (r, t) - ^ a2 exp [i (®2t - ¿r ) (¿R)2 5r

sPt ( ц) (¿r)2 дрз ( ц)

se

30

se

ee.

?(t)

6E

(r,t) - —a2exp[i(ú)2t - kr))(() sin(e)ee.

1

E(t) =-144E

^wave c 0 n 2/o ,-! \

5 R rc2 (2P2 + 3Pt)

sin e.

Зависимости амплитуды напряженности волнового поля для основной моды от угла

наблюдения, расстояния от капли до точки наблюдения г и напряженности внешнего электрического поля Е0 (при разных величинах последнего) представлены на рис. 1 и 2. Видно, что периодичность по углу 0 с увеличением расстояния до точки наблюдения сглаживается, а с увеличением напряженности внешнего поля становится более очевидной.

Квазистатическое приближение (см. (26)):

t (r, t )>

3E0 ( R R

a2 exp [i (m2t - kr )] >

5 ^(ц),-SPMee j+( R {4P3 (ц)е

SP3 (P-) -

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

Учитывая явный вид коэффициентов ц* и заменяя г ^ кг и г0 ^ кЯ, запишем напряженность электрического поля в волновом и квазистатическом приближениях для случая возбуждения основной моды (п = 2). Волновое приближение (см. (25)):

se Vr) 5 [ ' ' se

Поскольку R << r, понятно, что основной вклад в напряженность электрического поля в квазистатическом приближении определяется первой фигурной скобкой:

^tat (r,t) - Що (f ) a2exp [' (®2t - kr )]

|2cos(e )er + sin (e )-e J. Амплитуда такого поля определена как

E« -6Eo íRT^(t + 3co^ 1

i0 V r

■) R (t+3cos2(e)))

Поскольку г0 = кЯ << 1, очевидно, что основной вклад в напряженность электрического поля в волновом приближении определяется первым слагаемым:

Подставив сюда явный вид частоты собственных колебаний капли

2 п(п- 1)(п + 1)(п + 2) О

П (р2п + (п + 1)Р1) Я3' амплитуду напряженности электрического поля в волновом приближении для основной моды запишем как

Зависимости амплитуды напряженности квазистатического поля е^ для основной моды от

расстояния от капли до точки наблюдения г, угла 0 и напряженности внешнего электрического поля Е0 представлены на рис. 3 и 4. Качественно зависимости такие же, как на рис. 1 и 2, но степень периодичности меняется.

Квазистатическое электрическое поле определяет электрическое давление на поверхность капли, которое ранее неоднократно рассчитывалось для различных геометрий заряженной поверхности жидкости исходя из других соображений (см., например, [4-6, 9]).

ЗАКЛЮЧЕНИЕ

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

Авторы выражают благодарность Н.А. Бога-тову.

Работа выполнена при поддержке грантов РФФИ № 14-01-00170-а и 14-08-00240-а.

ЛИТЕРАТУРА

1. О' Konski C.J., Thacher Н.С. The Distortion of Aerosol Droplets by an Electric Field. JPhys Chem. 1953, 57, 955-958.

2. Taylor G.I. Disintegration of Water Drops in an Electric Field. Proc R Soc Lond A. 1964, 280, 383-397.

3. Cheng K.J., Chaddock J.B. Deformation and Stability of Drops and Bubbles in an Electric Field. Phys Lett A. 1984, 106(1-2), 51-54.

4. Френкель Я.И. К теории Тонкса о разрыве поверхности жидкости постоянным электрическим полем в вакууме. ЖЭТФ. 1936, 6(4), 348-350.

5. Григорьев А.И., Синкевич О.А. К механизму развития неустойчивости капли жидкости в электростатическом поле. Изв. АН СССР. МЖГ. 1985, (6), 10-15.

6. Cheng K.J. Capillary Oscillations of Drop in an Electric Field. Phys Lett A. 1985, 112(8), 392-396.

7. Inculet I.I. Breakup of Large Water Droplets by Electric Fields. IEEE Transactions on Industry Applications. 1989, 25(5), 945-948.

8. Grimm R.L., Beauchamp J.L. Dynamics of Field-induced Droplet Ionization: Time-resolved Studies of Distortion, Jetting, and Progeny Formation from Charged and Neutral Methanol Droplet Exposed to

Strong Electric Fields. J Phys Chem B. 2005, 109, 8244-8250.

9. Ширяева С.О., Григорьев А.И., Волкова М.В. О возможности зажигания коронного разряда в окрестности нелинейно-осциллирующей во внешнем электростатическом поле электропроводной капли. ЖТФ. 2005, 75(7), 40-47.

10. Ширяева С.О., Григорьев А.И. Скаляризация векторных краевых задач гидродинамики. Ярославль: Изд. ЯрГУ им. П.Г. Демидова, 2010. 180 с.

11. Абрамовиц М., Стиган И. Справочник по специальным функциям. М.: Наука, 1979. 830 с.

Поступила 10.07.14

Summary

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

The carried-out calculations afforded an analytical expression, in the linear approach on indignation amplitude, for the intensity of an electric field of an oscillating uncharged drop of a nonviscous conductive liquid, placed in an electrostatic field. Estimates of the field intensity values, both in the wave approach (at long distances from a drop) and in the quasistationary approach (at a drop surface) are carried out.

Keywords: spheroidal drop, uniform electrostatic field, oscillation, wave and quasistationary approximations.

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