Асимптотический анализ электромагнитного излучения незаряженной сфероидальной капли, осциллирующей в электростатическом поле Текст научной статьи по специальности «Физика»

Асимптотический анализ электромагнитного излучения незаряженной сфероидальной капли, осциллирующей в

электростатическом поле

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

Ярославский государственный университет им. П.Г. Демидова, ул. Советская, 14, г. Ярославль, 150000, Россия, e-mail: [email protected]

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

Ключевые слова: незаряженная капля, электростатическое поле, осцилляции, излучение электромагнитных волн.

УДК 551.594

ВВЕДЕНИЕ

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

В общем случае электромагнитное излучение объекта состоит из трех компонент: дипольной (самая интенсивная компонента), квадрупольной (вторая по интенсивности компонента, на много порядков величины слабее дипольной) и магнит-но-дипольной (самая слабая, на много порядков слабее квадрупольной компоненты) [2, с. 230]. Когда имеется дипольное излучение, квадру-польное на его фоне незаметно [2]. Когда нет дипольного излучения, самое интенсивное излучение - квадрупольное. Именно оно и предсказано в [3] в асимптотических расчетах первого порядка малости по амплитуде осцилляций для заряженной капли (дипольное излучение в расчетах первого порядка малости в случае заряженной капли запрещено условием неподвижности центра масс). С учетом поправок [4] квадрупольное излучение и принимается за электромагнитное излучение от заряженной капли (см. также [5]). Но, как показано в [6-7], электромагнитное излучение от заряженной капли можно рассчитать и в нелинейных асимптотических приближениях, когда запрет на дипольное излучение снимается и интенсивность

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

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

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

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

Рассмотрим каждую из половинок капли. Проинтегрировав поверхностную плотность заряда V по поверхности каждой половинки капли (51 и S2), получим величины разноименных поляризационных зарядов в виде:

© Ширяева С.О., Григорьев А.И., Колбнева Н.Ю., Электронная обработка материалов, 2017, 53(3), 33-40.

= { = | У(1Б1; Ч- = \ = | ;

(1)

где - положительным и - отрицательным поляризационные заряды, причем |'+| = |.

Хотя, как уже отмечалось, заряды непрерывно распределены по поверхности половинок капли, каждому из них можно поставить в соответствие «эквивалентный» заряд, расположенный внутри капли. Положение таких «эффективных» зарядов

Я рассчитывается по формуле [8]:

- 1 г (2)

где г - радиус-вектор; - (в нашем случае) площадь одной из половинок капли, по которой ведется интегрирование; dq - заряд элемента поверхности капли.

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

Учтем, что на поверхности капли существует капиллярное движение, возбуждаемое уже тепловым движением молекул [9]. Амплитуда таких осцилляций весьма мала: -V кТ / о, где к - постоянная Больцмана; Т - абсолютная температура. При температурах порядка комнатной такая амплитуда для всех жидкостей меньше ангстрема. Однако амплитуды некоторых из мод осцилляций под влиянием внешних условий могут увеличиваться. Поскольку индуцированные заряды распределяются по возмущенной осцилляциями поверхности капли, положения «эффективных» зарядов при указанных осцилля-циях тоже осциллируют, но с другими амплитудами и частотами. Будут осциллировать со временем и величины поляризационных зарядов, так как при осцилляциях меняется и величина поверхностной плотности зарядов на половинках капли. В результате при осцилляциях поверхности капли как величины, так и центры индуцированных зарядов осциллируют со временем возле положений равновесия и, следовательно, излучают электромагнитные волны.

Не ограничивая общности, будем рассматривать только осесимметричные осцилляции, что

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

Интенсивность дипольного излучения для ускоренно движущегося заряда ' определяется выражением [2]:

= _2_ (*)

3с3 дх2 '

где d - дипольный момент, определяемый формулой:

d (X ) = ' (х )• Я' (X); где Я (х) - вектор положения центра «эффективного» заряда (2).

Для того чтобы определить величины и положения центров «эффективных» зарядов, рассмотрим эволюцию во времени формы осциллирующей поверхности незаряженной капли. Все расчеты задачи будем проводить в сферической системе координат (г, 9, ф) с началом в центре масс капли в безразмерных переменных, в которых Я = р = с = 1.

Выражение для формы равновесной поверхности капли, совпадающее с точностью до квадрата эксцентриситета е2 с уравнением вытянутого сфероида, имеет вид [10]:

г (0)« 1 + А е2 Р2 (ц) + о(е4 );

е =,

9Е2

' 16п

Равновесный потенциал Ф(г, 0) в окрестности проводящего незаряженного жесткого сфероида, помещенного в коллинеарное внешнее однородное электростатическое поле Е0, может быть получен, например, переходом к сферическим координатам с последующим разложением по степеням эксцентриситета е известного выражения, полученного в [11, с. 11], в сфероидальных координатах. Точно такое же выражение для Фщ (г, 0) можно получить и прямым решением

соответствующей электростатической задачи, используя разложение по степеням е в окрестности исходной сферы, на основе метода возмущений [12].

Положения центров «эквивалентных» зарядов и их величины для жесткого сфероида легко вычисляются в виде (например, для положительного заряда):

Я(+)= 2яГ1 + Ае2\ '+= 1Я2Ео (1 + £е2].

Если рассмотреть не жесткий сфероид, а осциллирующую каплю, то положение Rz ' и

величина q+ будут изменяться, и можно сформулировать задачу об электромагнитном излучении при осцилляциях поверхности капли (при изменении R^ и q+).

МАТЕМАТИЧЕСКАЯ ФОРМУЛИРОВКА ЗАДАЧИ

Пусть в начальный момент времени t = 0 равновесная сфероидальная форма капли с эксцентриситетом e претерпевает виртуальное осесим-метричное возмущение 5(0, t) фиксированной амплитуды в, существенно меньшей радиуса капли. Примем, что уравнение линии, вращением которой получается поверхность капли, в сферической системе координат с началом в центре масс капли в безразмерных переменных имеет вид:

r(9, t) = r(9) + 5(0, t) = 1 + 3e2P2 (r) + 5(0, t) -

- 1 + e2 h(0) + 5(0, t);

|5(0, t)| << 1; r - cos 0, (3)

где P«(r) - осесимметричный полином Лежандра порядка n.

Движение жидкости в капле будем полагать потенциальным, и примем, что поле скоростей движения жидкости V (r, t) = Vy (r, t) полностью

определяется функцией потенциала скорости y(r, t), величина которого имеет тот же порядок

малости, что и амплитуда осцилляций поверхности капли y(r, t) ~ 5(0, t) ~ s.

Математическая формулировка задачи имеет вид [3, 5]:

Ay(r, t) = 0; ДФ(г, t) = 0; (4)

r ^ 0: y(r, t) ^ 0; r ф(г, t) ^-E0r cos 0;

r = r (0) + 5(0, t):

a5 (0, t) ay(r, t) -

dt dr

(5)

1 dy(r, t)

ae

dr (e) + a5 (e,t)'

ae

ae

AP-1 (W )2 + PE = Pc; dt 2V ^ E °

Ф(г, t) = Ф 5 (t).

(6)

(7)

Дополним выписанную систему естественными интегральными условиями: неизменностью полного объема капли (следствие несжимае-

(8)

(9)

мости жидкости), неподвижностью ее центра масс, а также незаряженностью капли:

f 2 4

I r dr sin 9йШф = — п;

к 3

| r • r 2dr sin 9d9dф = 0;

к

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

—ф(й, VO)d^ = 0;

4п S

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

В выражениях (6)-(7), (9) введены следующие обозначения: AP - перепад постоянных давлений внутри и вне капли в состоянии равновесия; ФД0 - постоянное вдоль поверхности капли значение ее электрического потенциала; PE = (VФ)2 /8п - давление электрического поля,

Pa = divn - капиллярное давление; n - единичный вектор нормали, определяемый выражением:

n =

v( r - r (e, t)- 5(e, t))

|v(r - r (e, t)- 5(e, t) )|

r=r (e,t)

В рассматриваемой задаче имеются два малых параметра: эксцентриситет равновесной поверхности капли е и амплитуда в ее осцилляций. Для определенности установим фиксированное соотношение между малыми параметрами, то есть примем, что е2 ~ в, сводя тем самым задачу к одному малому параметру в. При дальнейшем анализе задачи будем учитывать в разложениях слагаемые, имеющие порядок малости не выше, чем ~ в . Это означает, что слагаемые ~ Е0 • е2 • в ~ 85/2 будут отброшены. Учитываемый

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

Ввиду того что e2

s, а E0 ~ e ~ s1/2, для

удобства дальнейших разложений, чтобы иметь возможность в конечных выражениях легко вернуться к физическим величинам Е0 и е2, введем формальные параметры-маркеры рЕ, ре в соответствии со следующими соотношениями:

Ео - вЕ • 81/2,

е2 - ве • 8, где Ре = ве = 1.

Электростатический потенциал Ф (г, 0) в д2 , ч 2 / ч (15)

е"У / -^Мп (х ) + ®„м„ (х ) = ^

окрестности равновесного незаряженного сферо- дх

ида во внешнем электростатическом поле при 2 = , 2)

этом примет вид: п \ Л )>

1/2 ( 1 Л где с - безразмерная частота собственных коле-

Фе'(г,0) = 8 РЕгР1 (ц)1~"4 + (10) баний поверхности незаряженной сферической

капли. Следует отметить, что при Е0 ^ 0 выпол-

2 1 „3/2, 31 Л

Г-вкве Р 00+ Т-2 Р3 М

V

5 r V 2 r

няется также е ^ 0, то есть все решение строится в окрестности сферы.

Решением уравнения (15) являются гармони-Выражения для возмущения равновесной ческие функции времени t: формы поверхности осциллирующей сфероидальной капли £(9, t) и решения уравнений (4) Mn (t) = an exp['(pnt + bn)] + к.с.; (n > 2), (l6)

для гидродинамического w (r, t) и электростати- ,

v ' где an и bn - вещественные константы, определя-

ческого Ф^ (r, t) потенциалов, происходящих емые из начальных условий; аббревиатура «к.с.» из-за капиллярного возмущения свободной означает слагаемые, комплексно сопряженные к

поверхности, представим в виде рядов по поли- выписанным.

номам Лежандра: Таким оф^ом используя (3), (11) для обра-

зующей формы поверхности осциллирующей £(9 t) = е VM (t)Р ( )• (11) незаряженной капли, находящейся во внешнем

n=0 ' '' однородном электростатическом поле, получим

аналитическое выражение:

W(r,0,t) = е V Dn (t)rnPn (ц); (l2) » (17)

П n=0 Л ) r(0,t) = l + веек(0) + еVMn(t)Pn(ц), (l7)

n=2

Ф (r 0 t) = ез/2 V F (t)r"(n+l)p (ц) (l3) в котором амплитудные коэффициенты Mn(t) ' ' > n=0 nW nW' определяются (14).

Представляя электростатический потенциал в РАСЧЕТ виде разложения

Для определения коэффициентов Dn, Mn в Ф(г , t) = Ф eq (r, 0) + Ф£ (r, 0, t), (l8)

решениях (11), (12) из уравнений (5), (6), (8) в 3/2

первом порядке малости по в получаем систему: выделим краевую задачу порядка в для определения электрического потенциала Ф^(г, 9, t), свя-

r = i: 5£(0 t) (r 0 t) занного с возмущением равновесной поверхно-

' ——— =-5-• сти капли, осциллирующей во внешнем электри-

_/r 0 t) ческом поле. Система уравнений для определе-

wv , , / = ( + Lq))(0 t)• ния коэффициентов Fn в (13) получается из (4),

5t (7), (9) и (14) группировкой слагаемых: ~ в3/2.

п п Используя разложения (11), (13), получим

Í£(0,t)sin0d0 = 0; í£(0,t)cos0sin0d0 = 0, следующие соотношения между коэффициен-

0 0 ' тами:

F (t) = 0; Ф^ = 0;

Lq Иsin0^1, F" (t) = 3вE (ц+_M_1 (t) + ц_+M+1 (t)); (n >l). (l9)

sin 0 50 V 50) v '

+ n +1 _ n

которая позволяет вывести следующие соотно- ц" =-, цп =-.

шения для коэффициентов: 2n +1 2" +1

Наконец, подставив (19) в (13) и перейдя от

M0 (t) = 0; Ml (t) = 0; D0 (t) = 0; Dl (t) = 0; комбинации Р^в12 к величине E0, запишем выражение для компоненты электрического потенци-D (t) = i 5Mn(t); (" > 2) (14) ала Фц(г, 9, t) в виде:

" " 5T0 Ф£ (r, Q t) = (20)

и однородное дифференциальное уравнение для » , ( ) _ ( ) -("+1 ( )

нахождения коэффициентовMn(t) при n > 2: = 3Е0еV(°_iM"_i() + ц"+1М"+1(t))r Р" (И>

ВЕЛИЧИНЫ ПОЛЯРИЗАЦИОННЫХ ЗАРЯДОВ

Величины поляризационных зарядов каждой из половинок возмущенной поверхности капли г(9, г) определяются уравнениями (1), если принять, что V - v(9, г) - поверхностная плотность заряда на возмущенном сфероиде, а 51 и определены как:

Si -

S2 -

r = r (e, t );0 < e < —;0 < ф < 2n

r = r(e, t);— < e < п;0 < ф < 2n

q+ = Í

we,t) 2

-^f r sin ededф;

( er )

^t ) = (n^ V<

r=r (e,t)

(i-(sPeA (e)+5 (e, t)))

j sPe am), ee

ae

ae

(П, er )

r(e,t)

= 4П E0 [ P1 (^)+15 e2 (1 +15c°s (2e))P1 (n)+

q+ = -4 E0

i да

1+77 e2 + 2sjMn (t )G (

n=1

(25)

Gi

1 (n)=(+I)( n2Fn-1 +(n + 7)(n + 2)F

n/2

Fn = j Pn (ц)sin e

Рассмотрим положительно заряженную половинку капли, используя поверхностную плотность заряда V = v(9, г) на возмущенной поверхности капли г(9, г):

1

n-1

(-1) ^ (n -1)!

,n-1/ I n - 1

(n +1)

(n = 0); 2 (n = 2k + 1);

(21)

(22)

Подставив в (22) разложение (18), с учетом (10), (17) и (20), а также вектора нормали для возмущенной поверхности капли:

Г1 -1Г 8Ре 5Ш42'

21 50 50

0 (п = 2к).

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

ИЗЛУЧЕНИЕ ЭЛЕКТРОМАГНИТНЫХ ВОЛН

Интенсивность электромагнитного излучения от единичной капли I в соответствии с известным выражением [3] запишем в виде:

I

I = -

3c3

(( )■ q±(t))

2

(23)

получим выражение для поверхностной плотности индуцированного заряда. В (23) ег и

е0 - орты сферической системы координат. Затем, учитывая вид функции 5(9, г) (11), решение (18), (10), (20) для потенциала Ф(г, г) и

переходя от введенных формальных параметров рЕ • в12, рЕ • в к величинам Е0 и е2, запишем подынтегральное выражение в (21) на возмущенной поверхности вытянутой сфероидальной капли (17):

где индексом «max» обозначены максимальные значения функций.

Для определения положения центра индуцированного заряда половинки капли учтем, что радиальный орт er сферической системы координат связан с ортами декартовой системы координат соотношениями:

er = ex sin e cos ф + ey sin e sin ф + ez cos e.

Поскольку мы рассматриваем осесимметрич-ные осцилляции, то смещения центров зарядов капли в плоскости x, y не происходит:

Rqx = Rqy = 0.

(26)

Учитывая выражение (26), запишем проекцию вектора смещения центра положительного индуцированного заряда вдоль оси г согласно (2) в виде:

(24)

М 1 / \

+2вХМп (г)(ц) + (п + 1)(п + 2)Ри+1 (ц)

п=1 (2п +1) '

Подставив (24) в (21) и интегрируя по половинке сфероида 51, найдем величину положительного индуцированного заряда на возмущенной поверхности капли г(9, г):

f f,

(e,t) 3

-—-r c°s e sin ededф.

R =

Откуда, используя выражения для поверхностной плотности (22) индуцированного заряда и для вектора нормали (23), получим:

Rqz — ~

qz 3

i оо

1 + -е2 + еXMn (t)(3G2(n)-2G(n))

3 n—2

G2 (n) =

(n+1)

(2n+1)

n(n-1) (4n3 +10n2 + 2n-3) (и+2)(n+3)

-17 1-Fn 1-

(27)

((-1)

Fn-2-

(2n-1)(2n+3)

(2n+3)

Fn+2

X ^ — 1; е «1; — 0;

TE 3 8t

(28)

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

Удовлетворяя начальным условиям (28), для вещественных констант ап и Ьп в выражении (16) получим значения:

a —1 h.8 .;

n 2 3 n,3 '

— 0; (3 eE, n — 0,1,2,...); (29)

где 5п;- - символ Кронекера.

Подставив выражения (29) в (16), запишем амплитуды первого порядка малости в выражение для формы поверхности осциллирующей капли в виде:

M0 (t)= Mi (t) = 0;

Mn (t) = 5Я,jhj cos); ( gS; n > 2).

В итоге выражения для величины положительного индуцированного заряда (25) и для положения его центра (27) примут вид:

( 1 1 (30)

1 + е2 + 2е X hj cos)G1 (j)

15 jeS v

q+ — 3 Eo R2

Rqz ——R

qz 3

( 1

Í 1 +-е2 + е X h. cos(/)( (j) -2G (j))

(31)

j — Eo2R6e2

3c3

1 +-! е2 15

X hj®2 (3G2 (j )-2Gj (j))

V 3e

+ 2| 1 + 3 е2

X (3)

V 3e

Аналогичным образом для второй половинки сфероидальной капли определим выражение для смещения центра отрицательного заряда.

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

х = 0: £ (0) = 8X Нр ( ц);

Вычисляя вторую производную по времени от смещений центров индуцированных зарядов капли (31), пренебрегая в (30) слагаемыми ~ получим окончательное выражение для максимального значения интенсивности электромагнитного излучения незаряженной сфероидальной капли, осциллирующей во внешнем электростатическом поле:

ЧИСЛЕННЫЕ ОЦЕНКИ

Возможный источник электромагнитного излучения от облаков связан с осцилляциями конечной амплитуды мелких капель из диапазона наиболее часто встречающихся в облаке размеров от 3 до 30 мкм. Концентрация n таких капель в облаке ~ 10 см- [13-14]. Осцилляции большой амплитуды облачных капель могут быть вызваны различными причинами: коагуляцией; дроблением на более мелкие в результате процессов столкновения или реализации электростатической неустойчивости; гидродинамическим и электрическим взаимодействием близко пролетающих капель; аэродинамическим взаимодействием с развитой мелкомасштабной турбулентностью, характерной для грозовых облаков. Амплитуды колебаний, согласно данным натурных наблюдений [15-16], могут достигать десятков процентов от радиуса капли.

Проведем оценку интенсивности фонового дипольного электромагнитного излучения, когда смещение центров индуцированных зарядов связано с возбуждением моды j = 2. Для численных оценок примем: в = 0,1, h2 = 1, с = 73 дин/см, р = 1 г/см3, Е0 = 50 В/см (~ 2 • 10-4 Е0кр для R = 3 мкм и ~ 7 •Ю-4 Е0кр для R = 30 мкм; Е0кр - напряженность внешнего поля, критическая для реализации электростатической неустойчивости). Тогда из (32) для капли радиусом R = 30 мкм можно получить I ~ 2 • 10-31 ^.W на частоте « 100 кГц.

Интегральная интенсивность электромагнитного излучения облака, состоящего из N капель, будет в N раз больше, если все капли осциллируют синфазно. Это может иметь место при резком изменении напряженности внутриоблачного поля, что бывает при разряде молнии. Следует отметить, что с изменением размера капель мощность излучения от капли практически не меняется. Это можно увидеть и из (32), если расписать выражение для частоты осцилляций в размерном виде. Данным обстоятельством рассматриваемая задача отличается как от нелинейного (дипольного) излучения заряженной капли [6], так и от линейного (квадрупольного) излучения незаряженной капли во внешнем электростатическом поле [17]. Если капли осциллируют независимо друг от друга, интегральная интенсивность увеличится в VN раз.

На рис. 1 приведена зависимость интенсивности излучения от напряженности внешнего электрического поля. Видно, что с увеличением напряженности электростатического поля интенсивность излучения быстро увеличивается: при увеличении напряженности поля в 3 раза интенсивность излучения увеличивается на порядок.

60 90 120 150

П0> V/cm

Рис. 1. Зависимость интенсивности электромагнитного излучения единичной незаряженной капли, находящейся во внешнем слабом электрическом поле, от напряженности внешнего поля, рассчитанная при ] = 2, в = 0,1, к2 = 1, с = 73 дин/см, р = 1 г/см3, К = 30 мкм.

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

J_I_I_I_I_L

5 10 15 20 25 30

Я, цт

Рис. 2. Зависимость частоты электромагнитного излучения единичной незаряженной капли, находящейся во внешнем слабом электростатическом поле, от радиуса капли, рассчитанная при таких же значениях физических величин, что и на рис. 1, и Е0 = 50 В/см (~ 2 • 10-4 Е0кр для К = 3 мкм и ~ 7 • 10-4 Е0кр для К = 30 мкм).

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

наблюдений [13-14], достаточно высока ~ 103 м-3, а скорость их свободного падения имеет величину « 78 см/с. При такой скорости падения сквозь облако капель с радиусами от 3 до 30 мкм и с максимумом концентрации, приходящимся на диапазон от 3 до 7 мкм, гидрометеор будет испытывать ежесекундно около 22 столкновений, в результате чего в нем будут возбуждаться осцилляции мод с п е {2 * 30}.

Проведем оценку интенсивности излучения для второго источника излучения, полагая Е0 = 50 В/см (~ 1 • 10-3 Е0кр для К = 80 мкм и ~ 2 • 10-3 Е0кр для К = 150 мкм). Тогда интенсивность составит I ~ 1 • 10-31 на частоте около 8 кГц.

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

Е0, V/cm

Рис. 3. Зависимость интенсивности электромагнитного излучения единичной незаряженной капли, находящейся в слабом электростатическом поле, от напряженности внешнего поля, рассчитанная при ] = 20, в = 0,1, Н20 = 1,

90 110 130 150

R, urn

Рис. 4. Зависимость частоты электромагнитного излучения единичной незаряженной капли, находящейся во внешнем слабом электростатическом поле, от радиуса капли, рассчитанная при таких же значениях физических величин, что и на рис. 3, и E0 = 50 В/см (~ 1 • 10-3 Е0кр для R = 80 мкм и ~ 2 • 10-3 Е0кр для R = 150 мкм).

При оценке интегрального излучения из облака следует учесть, что концентрация гидрометеоров с R « 100 мкм ниже на шесть порядков, чем концентрация мелких капелек с R « 30 мкм. В итоге по порядку величины интегральная интенсивность излучения облака определится именно мелкими каплями.

ЗАКЛЮЧЕНИЕ

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

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

ЛИТЕРАТУРА

1. Белоцерковский А.В., Дивинский Л.И. и др.

Активно-пассивная радиолокация грозовых и грозоопасных очагов в облаках. СПб: Гидрометео-издат. 1992. 286 с.

2. Ландау Л.Д., Лифшиц Е.М. Теория поля. М.: Наука, 1973. 504 с.

3. Калечиц В.И., Нахутин И.Е., Полуэктов П.П. ДАН СССР. 1982, 262(6), 1344-1347.

4. Богатов Н.А. Сборник докладов VIМеждународной конференции «Солнечно-земные связи и физика предвестников землетрясений». Россия, Петропавловск-Камчатский, 9-13 сентября 2013. ДВО РАН. С. 22-26.

5. Grigoryev A.I., Kolbneva N.Yu., Shiryaeva S.O. Surf Eng Appl Electrochem. 2015, 51(6), 530-539.

6. Ширяева С.О. ЖТФ. 2002, 72(4), 15-19.

7. Григорьев А.И., Ширяева С.О., Жаров А.Н., Коро-мыслов В.А. ЭОМ. 2005, (4), 24-35.

8. Френкель Я.И. ЖЭТФ. 1936, 6(4), 348-350.

9. Фихтенгольц Г.М. Курс дифференциального и интегрального исчисления. Том 2. М: Наука, 1970. 800 с.

10. Cheng K.J. Phys Lett. 1985, A112(11), 392-396.

11. Ландау Л.Д., Лифшиц Е.М. Электродинамика сплошных сред. М.: Наука, 1982. 621 с.

12. Найфе А.Х. Методы возмущений. М.: Мир, 1976. 455 с.

13. Мазин И.П., Шметер С.М. Облака. Строение и физика образования. Л.: Гидрометеоиздат, 1983. 280 с.

14. Мазин И.П., Хргиан А.Х., Имянитов И.М. Облака и облачная атмосфера. Справочник. Л.: Гидро-метеоиздат, 1989. 647 с.

15. Beard K.V., Tokay А.А. Geophys Res Lett. 1991, El(12), 2257-2260.

16. Beard K.V. Rev Geophys. 1987, 25(3), 357-370.

17. Ширяева С.О., Колбнева Н.Ю., Григорьев А.И., Артёмова Т.К. ЖТФ. 2015, 85(4), 20-27.

Поступила 05.01.16 После доработки 04.02.16 Summary

The oscillations of an uncharged spheroidal drop of a conducting liquid in a uniform electrostatic field are considered. It is shown that the oscillations of the centers of inducting charges of different signs connected with the surface oscillation lead to the radiation of electromagnetic waves of the dipole type. A mathematical model of the electromagnetic radiation is constructed and an assessment of its intensity depending on the size of a drop and on the external field strength is carried out.

Keywords: uncharged drop, electrostatic field, oscillation, radiation of electromagnetic waves.