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

Journal of Siberian Federal University. Engineering & Technologies, 2016, 9(8), 1153-1161

УДК 534.551.5

Using the Effects of Acoustic wave Interaction with Charged by Atmospheric Aerosols to As-sess the Electrical State of the Atmosphere

Ilya E. Kuznetsov, Mark V. Biletov* and Dmitry V. Bulgin

Military Education and Research Centre of Military-Air Forces

«Military-Air Academy Named After Professor N.E. Zhukovsky and Yu.A. Gagarin» 54а Starykh Bolshevikov Str., Voronezh, 394064, Russia

Received 20.02.2016, received in revised form 15.07.2016, accepted 24.09.2016

The paper proposes a theoretical and experimental substantiation of the possibility of determining the electrical state of the atmosphere radio. Is developing a method of assessment of the charge of the investigated volume of the atmosphere, characterizing its electrical activity, based on the physical effect of radiation of electromagnetic energy moving with acceleration of charged aerosol particles under the influence of initiated by the sonar acoustic waves.

Keywords: electric state of the atmosphere, acoustic effects, charged particles.

Citation: Kuznetsov I.E., Biletov M.V., Bulgin D.V. Using the effects of acoustic wave interaction with charged by atmospheric aerosols to assess the electrical state of the atmosphere, J. Sib. Fed. Univ. Eng. technol., 2016, 9(8), 1153-1161. DOI: 10.17516/1999-494X-2016-9-8-1153-1161.

© Siberian Federal University. All rights reserved Corresponding author E-mail address: m.biletov@yandex.ru

*

Использование эффектов взаимодействия акустической волны

с заряженными атмосферными аэрозолями для оценки электрического состояния атмосферы

И.Е. Кузнецов, М.В. Билетов, Д.В. Булгин

Военный учебно-научный центр Военно-воздушных сил

«Военно-воздушная академия имени профессора Н.Е. Жуковского и Ю.А. Гагарина» Россия, 394064, Воронеж, ул. Старых Большевиков, 54а

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

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

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

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

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

Недостатком контактных способов является низкая пространственно-временная дискретность, невозможность получения информации в больших объемах пространства, искажение

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

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

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

Предлагаемый метод основан на зависимости интенсивности электромагнитного излучения аэрозольных частиц, частиц облаков и гидрометеоров (в дальнейшем просто заряженных частиц) от величины их заряда и ускорения движения. Известно [2-4], что при движении с ускорением электрических зарядов (заряженных частиц) они излучают электромагнитные волны. Если определенный объем атмосферы, содержащий заряженные частицы, подвергнуть воздействию акустических волн, то частицы будут вовлечены в движение атмосферного воздуха, вызванного акустическими волнами. Вследствие этого заряженные частицы, находящиеся в пределах данного объема, станут источниками излучения электромагнитных волн. Напряженность электрического поля электромагнитных волн, излучаемых совокупностью N частиц, находящихся в пределах объема, подвергшегося акустическому воздействию, является суперпозицией напряженностей E,(t):

R

N q,a, (t--'-)sin р

E(t) = 2--*-aR', (1)

i=1 4 ne0c R'

где qi - величина заряда частицы; a,(t) - ускорение движения частицы; е0 - электрическая постоянная (диэлектрическая проницаемость воздуха), е0 = 8,9 •lO-12 Ф/м; с - скорость распространения электромагнитной волны в вакууме; Pi - угол между направлением движения частицы и направлением на точку, в которой наблюдается электромагнитная волна; S - коэффициент ослабления напряженности электрического поля волны.

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

d x ,, dx .—. ,,

т— + f— = Fx(t), (2)

dt dt

где m - масса частицы; f - динамическая вязкость воздуха; x - координата частицы вдоль направления действия силы Fx(t).

Обе силы, входящие в левую часть уравнения (2), при постоянных в течение действия силы Рх(ф) плотности и форме частицы будут зависеть от её размеров. Размеры частиц атмосферного аэрозоля, облаков, туманов, твердых и жидких осадков являются случайными величинами, имеющими плотность распределения /(г), где г - размер частиц. Кроме того, величина зарядов частиц также случайная величина, описываемая законом распределения/^). Поэтому Е(ф) - случайный процесс. В силу этого оценка суммарного электрического заряда аэрозольных частиц, находящихся в объеме атмосферы, подвергнутому акустическому воздействию, может быть найдена с использованием характеристик случайного процесса Е(ф). Полное описание случайного процесса Е(ф) дается многомерными плотностями вероятности. Для оценки суммарного электрического заряда аэрозольных частиц объема атмосферы будем использовать математическое ожидание и корреляционную функцию Е(ф). Такой подход основывается на свойствах процесса Е(ф), зависящих от характеристик акустического воздействия. Зондирующие сигналы акустического локатора представляют собой последовательность импульсов, имеющих одинаковую длительность и период повторения. В каждом из импульсов мгновенные значения колебаний давления определяются соотношением

р(0 = Р • Р(Г) • ъ1п(ю? + (3)

гдер(ф) - мгновенное значение давления; Р - амплитудное значение давления; А(ф) - огибающая импульса

Г1 при 0 < X <т„

Р(1) = ]0 ' . ^ " , (4)

[0 при X > ги

где ти - длительность импульса; ю=2 ж/- круговая частота изменения давления; /- частота; ф0 -начальная фаза, характеризующая давление в импульсе в момент времени ф = 0.

Распространение зондирующего сигнала вызовет колебательное движение воздуха, имеющее мгновенную скорость

К«) = Ж ,(5)

РСзе

где р - плотность воздуха; Сзв - скорость распространения акустической волны; рСзв - удельное акустическое сопротивление воздуха. В соответствии с выражениями (4) и (5)

Р • Р (/)

V(/) =-ът(®/ + ср0) = КаР(/) ът(®/ + ср0 ), (6)

с зв -Р

где V = - амплитуда скорости колебательного движения воздуха.

РСзв

Амплитудное значение скорости колебательного движения частиц воздуха можно связать с акустической мощностью локатора Ш:

V (t) = -1 R\

--R

We 2

■ „ • F (t) • sin(©t + (7)

2жрСзв

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

+ fV (t ) = Fv (t). (8)

dt

Если принять форму частиц сферической, а силу сопротивления - стоксовой, то уравнение (8) может быть приведено к виду

2 p-r^ Щ) + v^(t) = v(t), (9)

9 г) dt

где Vr(t) - мгновенная скорость колебательного движения частиц; V(t) - мгновенная скорость колебательного движения воздуха, определяемая формулой (7); р - плотность аэрозольных частиц; r - радиус частиц; ц - динамический коэффициент вязкости воздуха.

Уравнение (9) является линейным дифференциальным уравнением с постоянными коэффициентами. Его решение будет простым в том случае, когда V(t) будет гармонической функцией, т.е. будет меняться по синусоидальному или косинусоидальному закону на бесконечном интервале времени. Вид функции V(t) определяется видом зондирующего сигнала p(t) (формулы 3-4). Реальные зондирующие сигналы акустических локаторов имеют конечную длительность и форму, отличающуюся от гармонического колебания, но они могут быть полностью представлены своими спектрами. Линейность уравнения (9) и возможность представления зондирующих акустических сигналов своими спектрами позволяют свести его решение к решению алгебраического уравнения. Применив прямое преобразование Фурье, получим

'2 pr2 . ^

--J0 +1

9

Sv (©) = Sv(©), (10)

где Бу (ю ) и Бу(ю) - комплексные спектры У() и ¥(1) соответственно. Решив (10) относительно ¡&у (ю ), после преобразования получим

Буг (©) = К(©), (11)

где

1 л 2

ад = ——, г = 2 (12)

1 + ]у<о 9 ц

Из соотношения (11) видно, что объем атмосферы, подвергнутый акустическому воздействию, можно представить в виде четырехполюсника, характеризуемого передаточной функцией (12). Используя обратное преобразование Фурье, найдем скорость колебательного движения частиц

■] ^ 1 1

Vr (t) = — f K(a)SV (a)ejatda= — Г 1 e-Jarctgrca ■ Sv (a)ejmtda. (13) 171 - 1 + p2œ2

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

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

a, (t) =

dVri (t) dt '

(14)

где (г) - скорость колебательного движения частицы, определяемая выражением, например (14).

Импульсы имеют постоянную длительность, определяемую длительностью зондирующего сигнала. Амплитуда и фаза импульсов имеют случайные значения, так как коэффициент у, определяющий амплитуду Уг = у/1 + у2т2 и фазу атсг§ут, зависит от размера частиц - г. Кроме того, значения зарядов частиц q¡, связанные с амплитудой Е,() прямой зависимостью, являются случайными величинами, имеющими плотность распределения /^). Исследования, выполненные в [1], дают основание считать, что между величиной заряда q¡ и размерами части нет устойчивой связи. Величина заряда в первую очередь определяется физическими особенностями его возникновения. Случайный процесс Е(ф) флуктуационный. Его реализации имеют вид непрерывных функций времени. Корреляционная функция и энергетический спектр процесса Е(ф) будут определяться видом функции а,(ф), описывающей ускорение движения частицы.

Определим а;(/) для типового акустического воздействия, создаваемого акустическим локатором (формула (3)), используя формулу (14). Значение ¥гф найдем из выражения (13). Тогда, если ¡-я частица находится на удалении Я,, на основании формулы (7) имеем

R

V(t) = VF(t - -i-)sin(®t + Ро),

(15)

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

Подставляя выражение (15) в (13) и интегрируя, получим ускорение, приобретаемое заряженными атмосферными частицами:

ai (t) = —

м

W

2жрС.

"--Rt

-e 2 F (t

1

R,

С

(16)

cos(cot + (ро - arctg ym) + — sin((ро - arctg ym)e

7

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

a, (t) = A. (t )ej(№ -arCtgYa) ejmt,

(17)

где

A., (t) =

a

W

—R,

R,

Я,}! 2npC3

F (t - )

2

e

Тогда автокорреляционная функция будет иметь вид

K изл (О =

' ^2 -n-

Z q¿ái (t - 13i ) sin Pi

v1ж J

i=1

- , (18)

*

n

x E q . á . (t - . + z) sin fí . j = 1 J J 3 J J

~ 1

где Q =-2--постоянная величина.

4як0с 2 Л

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

Ризл(*) = Кшл(0) = ^^ Ё^А2(Г - г3,)аи2 рг. (19)

Величину суммарного заряда Q аэрозольных частиц можно получить по формуле

2л

2"

6 "i

P (г)

'излУЧ _. (20)

х ^i2(t -13¿) sin2 p. i=1

Практическая реализация предлагаемой методики была осуществлена по данным температурно-ветрового зондирования станции «Воронеж», характеризующим состояние атмосферы, и с учетом технических характеристик типового метеорологического акустического локатора: частоте акустической волны / = 2 кГц; длительности зондирующего импульса т = 250 мс; акустической мощности излучателя Ж = 1000 Вт. Полученные результаты представлены на рис. 1.

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

Ei, мкЕ/м

0,00004 0,00003 0,00002 0,00001 0

-0,00001 -0,00002 -0,00003 -0,00004

25 0

i ;--------

t, c

Рис. 1

0

0

3

энергетическим приемником равна Робн = 0,95 при вероятности ложной тревоги Рлт. = 102 и отношении сигнал/помеха по мощности Pc / Pn = 25.

Для экспериментальной проверки предлагаемого метода был создан измерительный комплекс, состоящий из генератора акустических колебаний, динамика, источника заряженных аэрозольных частиц, приемника электромагнитного сигнала, аналогово-цифрового преобразователя, электронно-вычислительной станции с программным обеспечением и звуковым редактором Gold Wave 5.55 (рис. 2).

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

Задачей эксперимента явилось подтверждение возможности определения заряда атмосферы радиоакустическим способом. Исследования выполнялись при различных условиях с наличием акустического сигнала и без него с различной частотой акустического зондирования и мощностью излучаемого сигнала. В качестве примера на рис. 3 приведен результат экспери-

Рис. 2

а б

Рис. 3. Без акустического воздействия (а); при акустическом воздействии с частотой 6 кГц на заряженный объем атмосферы (б)

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

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

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

[1] Современные исследования ГГО им. Воейкова. Т. 2 / под ред. М.Е. Берлянда, В.П. Ме-лешко. СПб.: Гидрометеоиздат, 2001. 345 c. [Modern research GGO im. Voeikov. V. 2. Ed. by M.E.Berliand, V.P.Meleshko. St. Petersburg, Gidrometeoizdat, 2001. 345 p. (in Russian)]

[2] Кузнецов И.Е., Билетов М.В. Сверхкраткосрочный прогноз опасных явлений погоды по радиолокационным данным с учетом энергетического состояния атмосферы. Наукоемкие технологии, 2012, 13, 3, 21-26 [Kuznetsov I.E., Biletov M.V. Very short-range forecasting of hazardous weather phenomena radar data with respect to the energy state of the atmosphere. High technologies, 2012, 13, 3, 21-26 (in Russian)]

[3] Пат. 2319981 (2008). РФ // Б. И. 2008. №8 [Patent 2319981 (2008) RF (in Russian)]

[4] Пат. 2491574 (2013). РФ // Б. И. 2013. №24 [Patent 2491574 (2013) RF (in Russian)]