CC BY
117
Новая дипольная динамическая модель шаровой молнии
В. Н. Сошников
Всероссийский институт научной и технической информации Российской Академии наук, ул. Усиевича, 20, г. Москва, 125315, Россия, e-mail: [email protected]
Впервые представлена дипольная динамическая модель (ДДМ) шаровой молнии (ШМ), которая просто, естественно и рационально объясняет большую часть наблюдений последней возможностью больших плотностей энергии, обусловленной дипольным разделением зарядов ШМ с притоком энергии атмосферного электрического поля за счёт поступательного движения ШМ. Теорема вириала, которая могла бы ограничивать энергию шаровой молнии, неприменима, так как ШМ не является замкнутой системой наподобие Солнечной или Галактической, и является сильно зависимой частью протяжённой во времени и пространстве большой системы.
УДК 53.03+537.523.9+537.84:533.98
ВВЕДЕНИЕ
Оценки плотности энергии в наблюдаемых шаровых молниях ведут в зависимости от их размеров - вплоть до (2-3) м и более - к значениям до ~ 1010 Дж/м3 и выше [1, 2]. Для объяснения таких огромных энергий и плотностей энергии предложено много изначально экзотических и даже фантастических моделей (см. также подробные обзоры [1] и [3]): может быть, ШМ является чёрной дырой; возможно наличие в ней ядерных реакций; в ШМ присутствует тёмное вещество; она содержит набор ридберговских атомов; является замкнутым круговым резонатором, заполненным световой энергией; обусловлена фазовым переходом к сверххолодной неидеальной плазме [4]; более реалистично - источником ее энергии являются коллективные электронно-ионные колебания [5] (что, вероятно, применимо к статичному плазмоиду, но не к движущейся типичной шаровой молнии); крайне усложнённые квантово-механические построения G.C. Dijkhius’a [1] и многочисленные другие гипотезы и модели этого рода. Ряд моделей связывают энергию ШМ с её заимствованием от линейной молнии [1, 2]: это модели динамических МГД-конденсаторов, электростатических микроконденсаторов, химические модели сгорания разреженных атомных кластеров Si, образовавшихся при ударе линейной молнии в почву. Главной трудностью этих моделей является неспособность объяснить иногда наблюдаемые огромные плотности энергии ШМ.
Газоплазменный диполь
Во всех случаях остаётся вне рассмотрения возможность больших энергий электростатического разделения зарядов с сильной поляризацией в облаке плазмы шаровой молнии атмосферным электрическим полем с образованием газоплазменного диполя (см. рисунок).
ПОСТАНОВКА ЗАДАЧИ
В разреженном газе ШМ (в огненном шаре fireball FB) c температурой молекул и ионов TFB несколько тысяч кельвинов и с электронными температурами до Te ~(8000-12000) K и более, типичными для термодинамически неравновесных разрядов низкого давления, возникающее отталки-
© Сошников В.Н., Электронная обработка материалов, 2012, 48(6), 54-64.
54
вающее действие положительного атмосферного электрического поля Eenv на ионы в задней части
ШМ может быть почти полностью скомпенсировано отдачей ионов от холодного окружающего плотного воздуха, что приведёт к движению ШМ в направлении атмосферного электрического поля, воздействующего и на ионы, и на электроны при наличии силы притяжения между ионами и электронами [6-8].
Простые оценки [7] показывают пренебрежимую малость потока убегающих ионов Eenvр д в
задней части ШМ, где рг- - подвижность ионов; n. - их концентрация, то есть наличие результирующей силы, двигающей вперёд молнию как единое целое. При поступательном движении ШМ во внешнем электрическом поле как электроны, так и ионы получают энергии, однако при умеренных скоростях движения ШМ только небольшая часть внутренней энергии идёт на поддержание ее движения. Основная же часть расходуется на нагрев частиц kTFB и кТе при хаотических столкновениях с
молекулами ШМ и друг с другом, принимая затем форму излучения. Благодаря дополнительной энергии внешнего электрического поля образуется энергетическая струя отдачи. В то же время для осуществления процессов ионизации и рекомбинации внутри ШМ необходимо значительное перекрытие ионного и электронного облаков, облегчаемое достаточно высокими температурами T. и Te.
Максимальный полный импульс, передаваемый окружающему воздуху позади ШМ за время t электрическим полем внешнего источника Eenv, равен ft, где f = |e|EenvN; e - заряд электрона,
N. - полное содержание ионов молнии. Внутри нее энергия должна распределяться путём электростатического дипольного взаимодействия между электронами и ионами соответственно энергиям их поступательного движения с одинаковой скоростью. Суммарная энергия, передаваемая струе отдачи и электронам ШМ, в одинаковых долях складывается из энергии, переносимой внешней струе отдачи и электронам молнии внешним источником.
Облако электронов, движущееся во внешнем электрическом поле, играет роль локомотива для движения ШМ с облаком ионов и нейтралов.
Вследствие больших сечений перезарядки при столкновениях ионов и нейтралов внутри ШМ она движется как единое целое с очень малым газодинамическим сопротивлением [8], получая энергию от поступательного движения электронов в атмосферном электрическом поле и теряя её главным образом при излучении, независимо от теоремы вириала. Последняя неприменима к ШМ, так как молния - всего лишь часть большой системы с протяжённым атмосферным электрическим полем, излучением и реактивной струёй отдачи холодного воздуха позади нее.
Приравнивая ускорение ионов вместе с нейтралами в дипольной ШМ под действием силы
притяжения ионов электронами e|NeN^(2rFB) на заднем крае молнии и ускорение ее электронной
части силой внешнего электрического поля |e|EenvNe, получаем очень простое соотношение, являющееся базовым в дипольной модели ШМ [6-8]:
Г
3
FB T
n m lelx2nL
e | | e L rrr
TFB
■ Е ,
env ’
(1)
где значения внешнего электрического поля Eenv , в котором находится ШМ, могут составлять от сотен В/см до нескольких кВ/см и более при электронных температурах Te ~(8000-12000) K и более в зависимости от радиуса молнии rFB; mav - средняя масса ионов и нейтралов ШМ при температуре TFB ; xe - молярная доля электронов при температуре электронов Te и давлении patm = 1 атм; Tenv = 293K - температура окружающего воздуха; nL - число Лошмидта при 293K; me, e, ne - соответственно масса, заряд и концентрация электронов, ne ^ nm, где nm - концентрация нейтралов в ШМ; ne ~ xenm . Это уравнение не исключает ускоренного движения молнии как целого.
Согласно уравнению (1) для удельной энергии ШМ получаем J ~ \e\nerFBEenv к, что придаёт правдоподобие возможности достижения очень больших плотностей энергии молнии (см., оценку предельных плотностей энергии в таблице).
55
Характерные ограничения на параметры ШМ
Тт = 1000 K, Т = 10 000 K, (yEe„v) = 1 кВ/см Тт =3000 K, Te =10 000K (YEenv ) =1 к^/см
г, см 10 100 1000 10 100
, см-3 х10-15 16,0 5,1 1,6 9,2 2,9
а 10,6 33,7 107 6,2 19,6
Espec, Д^/м3 5,15 х107 1,63х108 5,12 х108 2,9 х107 9,2 х107
е 6,8 х10-3 4,0 х10-2 ~ 2 х10-1 6,8 х10-3 4,0 х10-2
w, кВт 4,3 2,5 х102 1,2 х104 12,5 7,4 х102
w/SFB, Вт/см2 3,45 2,0 1,0 10 5,9
Максимальный тротиловый эквивалент ШМ 0,052 кг 0,164 Т 0,514 кТ 0,029 кг 0,092 Т
E = env = e [ПFB) ]/3, кВ/см 9,1 4,3 2,0 6,3 2,23
Y <0,11 <0,23 <0,51 <0,16 <0,45
Балансовые уравнения сохранения заряженных частиц (ионов и электронов), энергии ШМ, условие устойчивости термодинамически неравновесной молнии [7] (см. ниже), задание степени неравновесности а в замене ne ^ njа и уравнение (1) полностью определяют все параметры ШМ [7].
ИЗБЫТОК ПОЛОЖИТЕЛЬНОГО ЗАРЯДА
Вследствие амбиполярной диффузии шаровая молния обладает большим положительным зарядом [7], что эквивалентно поверхностному натяжению, удерживающему ионы и электроны от разлёта.
Избыток положительного заряда |e| AN. определяется коэффициентами диффузии электронов и ионов De, D. и их подвижностями рe, р. [9]:
D - D 1 lei AN.
AE ~ e . --------L!--2,
KI + Kl rFB (2грв)
где AE - электрическое поле амбиполярной диффузии, так что при произвольных разумных значениях Те ~ 10 000K, T- ~ 3000K, неравновесном пе ~2х1015см-3, rFB ~5см, ас°а ~ 2х10-16см2,
а“" ~ 2,5х10-15см2 получим De,|pe| » D.,|рг.|; |e|AN. ~4х10-2Кл при |e|N. ~ 0,17 Клс необходимостью
учесть AE в исходных уравнениях [7, 8] (см. уравнения (1) и (5) в [8]).
Отметим в связи с этим кажущееся парадоксальным положение: несмотря на положительный заряд ШМ, она движется в направлении положительного атмосферного электрического поля, притягивающего электронную часть молнии.
56
Оценки [7] показывают, что потери заряженных частиц при амбиполярной диффузии много меньше, чем их воспроизводство при ионизации в случае Te > (8000 f 10 000) K. Предполагается,
что большой положительный заряд, удерживающий электроны, должен приводить к устойчивости ШМ и при более низких температурах (Te < 8000 K).
УСТОЙЧИВОСТЬ ШМ
Показано [7], что шаровая молния устойчива только в случае сильного термодинамического неравновесия как по температуре Te » T = TFB , так и при много меньших, чем при равновесии, концентрациях электронов, соответственно существенно меньших светимостях.
Уравнение (1) может быть переписано в виде [7]
3 m E n 2 W л 3 m E n
____av env m ^2 = A__________ A —________av env m
n me lei n2e ne2 ’ n me lelw
rFB =
(2)
где nm - концентрация нейтральных атомов и молекул в ШМ; w - наблюдаемая термодинамически неравновесная мощность излучения молнии, принятая в опциях w равной (100-200)Вт [7, 8]; me, mav - соответственно масса электрона и средняя масса нейтрала в ШМ,
W = 4roFB ecTe4,
е ~ e(Te,10 см)
rFB
10 см’
(3)
а - постоянная Стефана-Больцмана, £ является здесь степенью черноты полусферического объёма [10], а2 — wjw = (njАеfbB ) , W- мощность равновесного излучения ШМ, ne — xe (Te)nm - концен-
трация электронов при температурах Tm, Te [11], ПввВ - действительная неравновесная концентрация электронов ШМ.
Предполагается, что баланс концентраций заряженных частиц brn2e = binenm , где br - коэффициент рекомбинации; bi - коэффициент ионизации; nm - концентрация нейтралей, то есть ne = nm b /br, ведёт к резкому уменьшению ne вследствие преобладающего по сравнению с br ~ const
резкого уменьшения bt с уменьшением температуры TFB < Te.
Отметим, что при выводе уравнения (2) пренебрегалось очень малым газодинамическим сопротивлением среды [8], которое воздействует на ШМ как целое.
В качестве примера рассмотрим опцию молнии для расчёта со следующими параметрами [7]:
T = 10 000 K; T = T = 1000 K; rFB = 10 см; е = 6,8х10-3;
e 7 i m 7 FB 77 7
vFB = 1,85 х10 2 м/с; nBFBB = nj а = 3,47 х1015 см
а = 49,3; Eenv = 46,6 В/см;
-3
; w = 200 Вт.
(4a)
(4б)
При фиксированном Eenv, заменяя обозначение n{FB) на ne, получаем обычное необходимое условие устойчивости ШМ на производную (при замене обозначения nBeFBB = xenm/а на xenm /а ^ ne с неравновесным ne):
0= dr^ = AW_ [Э( ln eTeA )ldTe - 2 d( ln Пе )/dTe ] = 0
dT n
1 - A
W d(ln е4в)
n
dr
'J ‘ T7\
и при полученной согласно (5) температуре T -12 500 K [7]:
1 - A
W d(ln erFB)
dr.
< 0,
2
n
e
(5)
(6)
d( ln er2B) д
dr
'Jl £7)
dr
'J‘ T7\
rfi)] = - 5
•[ln ( kFb )]
2r
2r 17
(7)
57
для случая больших оптических толщин с s (Te, rFB) <х т^в. Таким образом, уравнение (7) сводится к уравнению
d(lnеТ/)1дГе -2d(lnne)/dT = 0
(8)
с rFB = const. Условие устойчивости drFB/dTe ^ 0 эквивалентно условию минимума потенциальной энергии ШМ [7].
Этот результат можно интерпретировать как двухпараметрическое пространство (Te, а) или (Te, w) непрерывной последовательности устойчивых шаровых молний, в котором Te и w соотносятся параметрической траектории условия устойчивости (5).
Подобную процедуру можно применить при определении устойчивости ШМ в случае термодинамически равновесной плазмы, для которой
то есть
Т
n — n ~env, n « n ,
m L t ’ e m ’
± e
(9)
r
FB
BTl B = 1 mav Eenv 1
X2 ’ n me \e\ nLTenv
(10)
где xe = xe(Te) - молярная доля электронов [11], nL = 2,5х1019см 3. В этом случае, по крайней мере
при Te < 12000K,
drFB = B dT x
(
1 - 2T
e V
d ln xe dT
\
<< 0.
J
(11)
Таким образом, в ранее рассмотренном случае термодинамически равновесной плазмы [8], в дополнение к нереалистично высокой мощности излучения, ШМ абсолютно неустойчива.
В этих предположениях кажется возможным [7], что в зависимости от размеров шаровой молнии и концентрации электронов при её значениях в разумном интервале от njа~ (1015 —1016) cм-3,
удельная энергия ШМ может достигать значений до ~ 10-12 Дж/м3, и она может также сопровождаться присоединяющимся локальным разрядом линейной молнии.
В [7, см. табл. 4] кратко обсуждается также пример учёта реабсорбции излучения с переходом изменения черноты излучения £ от зависимости квадратного корня к линейной зависимости от оптической толщины.
СОПОСТАВЛЕНИЕ ОСНОВНЫХ ОСОБЕННОСТЕЙ ДДМ И В НАБЛЮДЕНИЯХ ШМ
Дипольная динамическая модель (ДДМ) согласуется с большей частью наблюдений [1, 2].
1. Движение шаровой молнии (ШМ) обычно прямолинейно, иногда с изломами, в широком интервале скоростей для различных ШМ и, по-видимому, определяется не скоростью ветра, а направлением и силой окружающего атмосферного электрического поля.
2. В ДДМ шаровая молния может двигаться только в положительном направлении напряжённости атмосферного электрического поля (электронным краем вперёд).
3. ШМ относительно “холодная” с температурой TFB ~ (1000 — 5000) K, но достаточной при радиусе rFB ~ (5-10) см для наблюдаемого прожигания дыр в оконных стёклах, с относительно слабым свечением мощностью ~(100 - 200) Вт, что свидетельствует о термодинамической неравновесности с электронной температурой Te » TFB .
4. Наблюдаемое прохождение ШМ сквозь неповреждённые оконные стёкла можно объяснить ее сильной дипольной поляризацией ШМ.
5. ШМ может появляться вне связи с линейной молнией.
6. Очень сильное атмосферное электрическое поле, необходимое для поддержания ШМ ДДМ, может приводить к смертельному поражению людей и животных электрическим током.
58
7. Позади ШМ ДДМ можно наблюдать поток холодного или слегка нагретого атмосферного воздуха реактивной струи отдачи, обусловленной силой f < eEenvNi.
8. Время жизни ШМ ДДМ определяется главным образом внешними факторами, изменением атмосферного электрического поля во времени и пространстве. Существуют и внутренние факторы в виде возможности возникновения неустойчивостей при ускоренном движении ШМ.
9. Разнообразие ШМ может объясняться наличием двух свободных параметров ДДМ, например радиуса и светимости.
10. Существует характерное расчётное расплывание порядка нескольких миллиметров как на переднем, так и на заднем крае ШМ.
11. Возникновение ШМ ДДМ обусловлено двумя простыми принципами самоорганизации: 1) наличие затравочного разряда с высокой степенью термодинамической неравновесности Te » TFB . Наиболее подходит для этого, возможно, микроволновой сфокусированный СВЧ-неравновесный разряд; 2) сильное протяжённое электрическое поле, выводящее затравочный разряд в состояние устойчивого поступательного движения. Не исключается возможность самопроизвольного возникновения затравочного разряда при локальном электрическом пробое в местах резко неоднородного атмосферного электрического поля.
12. Быстрое движение крупной ШМ (возможно, с ускорением?) может сопровождаться свистящим звуком, похожим на шум реактивного самолёта [2], вероятно, обусловленным реактивной струёй отдачи позади ШМ.
13. ШМ в области сильного атмосферного электрического поля может являться спусковым механизмом для сопутствующего локального атмосферного разряда, в большой степени усиливающего разрушительную мощь ее разряда.
14. Следует отметить, что, строго говоря, ШМ должна быть вытянута в направлении внешнего поля и иметь овальную форму, характерную для диполя, с неоднородной внутренней структурой, и рисунок является идеализацией первого приближения. Избыточный положительный заряд ШМ, удерживающий электроны, допускает значительную наблюдаемую изменчивость формы молнии при направленном движении с сохранением её устойчивости.
15. Из табл. 4 [7] следует, что при фиксированном rFB имеется тенденция к значительному уменьшению скорости ШМ с увеличением её температуры и энергосодержания.
16. Разрушительная сила взрыва в ДДМ определяется не только энергосодержанием ШМ, но и в не меньшей степени скоростью внезапного падения атмосферного электрического поля и ее энергоснабжения с результирующим схлопыванием диполя и последующим взрывным разлётом квазинейтральной плазмы.
Для экспериментального получения ШМ необходимо создать сильное направленное постоянное электрическое поле, воздействующее на сильно неравновесный затравочный разряд.
Импульс ШМ MBLvBL очень мал; разрушительная сила молнии обусловлена взрывом с высвобождением очень большой энергии электростатического разделения зарядов.
НЕКОТОРЫЕ ОБЩИЕ ЗАМЕЧАНИЯ
Оценки параметров дипольной динамической модели ШМ наряду с оценками [7, 8] показывают их качественное соответствие многочисленным и различным феноменологическим наблюдениям ШМ.
Шаровая молния рассматривается как особый вид движущегося электрического диполя с балансом ионизации в средней его части; ионы и электроны на концах диполя играют роль силовых оболочек огненного шара вместе с поступательным движением, удерживающим ШМ от рассеяния и схлопывания.
Тайна ШМ состоит в её уникальности, не имеющей аналогов среди плазменных объектов, с движением в атмосфере, сильной поляризацией, она является сильно неравновесным воздушным плазменным кластером и пополняет внутреннюю энергию за счёт очень протяжённого внешнего электрического поля. ШМ не имеет тесной аналогии со статическим неподвижным тлеющим разрядом в атмосферном электрическом поле [15, 16].
Несомненно, имеют право на существование также МГД- и кластерные (с присутствием реакционноспособных атомов Si, Fe и тому подобных) варианты моделей как не характерные короткоживущие низкоэнергетические плазменные объекты, тесно связанные с линейной молнией и заимствующие её энергию.
59
Имеется также играющий важную роль эффект притяжения избыточным положительным разрядом, эквивалентный поверхностному натяжению, с удержанием убегающих электронов.
Существование такого типа плазменных образований, не замкнутых в конечном пространстве, обусловлено их поступательным движением, обусловливающим приток энергии от внешнего электрического поля. Это делает неприменимой теорему вириала.
Потенциальная энергия шаровой молнии U < \e\EenvNes в атмосферном электрическом поле,
где Ne - содержание электронов в ШМ, s - произвольный пробег ШМ в электрическом поле, преобразуется в энергию излучения. Таким образом, в зависимости от s величина потенциальной энергии может становиться произвольно большой с отношением внутренней потенциальной энергии к внутренней кинетической энергии ШМ, никак не связанным с теоремой вириала. Отношение внутренних потенциальной и кинетической энергий ШМ определяется не теоремой вириала, а скоростью потерь энергии молнии, включая потери на излучение. Теорема вириала неприменима к ШМ, так как молния является только составной частью большой системы с протяжённым электрическим полем, излучением и реактивной струёй позади нее.
Ионизация ускоренными электронами происходит преимущественно в передней части шаровой молнии, рекомбинация преобладает в ее задней части, где скапливаются ионы. Могут возникать микропотоки, существенно искажающие Максвелловское распределение.
ШМ получает энергию вследствие ускорения электронов в атмосферном электрическом поле. Электроны очень быстро хаотизируются при столкновениях с нейтралами и друг с другом с увеличением температуры Те, что приводит к росту скорости рекомбинации и радиационных потерь. В то же
время ШМ ведёт себя как единое целое, а действующая на ионы сила внешнего электрического поля компенсируется реактивной струёй позади ШМ. Это указывает на то, что на формирование реактивной струи должна расходоваться энергия внешнего источника с мощностью |e| EenvNivF,B . Таким образом, атмосферное электрическое поле обеспечивает суммарное поступление энергии в ШМ с мощностью И EenvNeVFB , включая потери на излучение, и такую же мощность поступления энергии на создание реактивной струи. Нагрев электронного облака до некоторой температуры Te останавливает
скорость поступления энергии на некотором уровне, определяемом условием устойчивости ШМ.
Если m - масса всех электронов ШМ, а M - вся ее масса, то по закону сохранения полного импульса, то есть при равенстве импульса, получаемого электронами от внешнего поля (при их скорости v) до хаотизации, и импульса ШМ после хаотизации (со скоростью ШМ v2), mvx = Mv2 для энергии поступательного движения ШМ после хаотизации электронов получим
Таким образом, практически вся начальная энергия электронов e EenvNes, полученная на расстоянии
s передвижения шаровой молнии, переходит в тепловую энергию нагрева электронного облака с температурой Te.
Атмосферное электрическое поле в конечном счёте сообщает реактивной струе тот же импульс \e\EenvNet за то же время, что и импульс всей ШМ в прямом направлении, и, таким образом,
суммарный импульс, затрачиваемый внешним источником, равен нулю.
Механизм прекращения увеличения импульса ШМ можно объяснить передачей получаемого электронным облаком импульса образующимся при рекомбинации нейтралям, приобретающим импульс, направленный вперёд. Остаётся предположить, что эти нейтрали передают импульс молекулам холодного воздуха на переднем крае ШМ, образуя направленную вперёд струю отдачи, то есть создаётся дополнительное сопротивление движению молнии (помимо ранее отмечавшегося малого газодинамического тормозного сопротивления движению твёрдого тела) с парадоксальной направленной вперёд струёй отдачи холодного воздуха. При установившемся движении ШМ с постоянной скоростью атмосферное электрическое поле “тянет” электроны, связанные с ионами в одно целое в виде ШМ, с силой, равной силе сопротивления, зависящей от скорости молнии. Таким образом, энергия внешнего источника затрачивается приблизительно в равных долях на поддержание задней струи и на излучение с затратами энергии на поддержание передней струи. Однако при этом возможна пред-
ЭНЕРГИИ И ИМПУЛЬСЫ
60
положительная неполная (само)компенсация импульса, сообщаемого электронам ШМ внешним полем, и импульса торможения движению ШМ. Дисбаланс столкновительных процессов может приводить в той или иной степени к небольшому ускорению/замедлению ШМ с естественно конкурирующим укорочением времени жизни ШМ вследствие возникающей неустойчивости.
Отметим, что импульс, непосредственно передаваемый электронами молекулам окружающей среды при движении ШМ как целого с удержанием электронов положительным зарядом, очень мал вследствие малых сечений упругого рассеяния, малого отношения масс mjmmol и исчезающе малого газодинамического сопротивления. Таким образом, если перенос импульса тяжёлым нейтралям происходит только при их упругих столкновениях с электронами, то ШМ будет ускоряться силой f ~ вплоть до почти неограниченной скорости (вплоть до разрушения ШМ) без образования
передней струи, но с большим усилением задней струи. Разрушение шаровой молнии может происходить из-за возникающей неустойчивости при ускорении/замедлении движения ШМ.
ЗАДНЯЯ РЕАКТИВНАЯ СТРУЯ
Задняя реактивная струя отдачи может содержать незначительную примесь горячих молекул и убегающих ионов, несколько уменьшая скорость движения ШМ. Кроме того, можно предположить, что при больших значениях напряжённости атмосферного электрического поля (но ниже порога стримерного пробоя ~ (1v3) кВ/см) возможно возникновение нестримерной атипичной линейной молнии в виде ШМ со следом позади неё, образующим расширяющийся канал линейной молнии.
Позади ШМ должен возникать реактивный кильватерный след окружающего воздуха, обусловленный импульсом отдачи молнии. Однако в ряде случаев импульс отдачи может быть сильно переоценен, что требует дальнейшего объяснения и что может быть обусловлено неполным поляризационным воздействием атмосферного электрического поля на ионы с соответствующим уменьшением эффективной силы воздействия на электроны и ионы. Это приведёт к необходимости введения поляризационных поправок, увеличивающих электрическое поле, необходимое для поддержания ШМ, с соответствующим уточнением уравнения (2) и других соотношений, содержащих Eenv.
ПОЛЯРИЗАЦИОННАЯ И РЕКОМБИНАЦИОННАЯ ПОПРАВКИ
Можно предположить, что в вышеуказанных случаях следует ввести корректирующий фактор у в напряжённость электрического поля Eenv в уравнениях (1), (2) и всех других: Eenv ^ yEenv.
Ключевым параметром в величине у является отношение |e|Eenv = a1 к силе кулоновского взаимодействия заряженных частиц на среднем расстоянии между ними (или на длине свободного пробега) а2 (см. [8, табл. 7]. Если axja2 ~ Eenvj|e||^n(FB) J/3 < 1,то у ~ a1/a2, и если a2/a1 < 1, то y-1. Из табл. 7 [8] следует, что у-1 может достигать значений у-1 > 10. В случае больших диаметров, при возрастающих неравновесных пе в интервале ne ~ (1015 - 1016)см-3, возможно значительное увеличение Eenv до сотен и тысяч В/см с увеличением плотности энергии ШМ до j > 108 Дж/V.
Согласно уравнению (1) электрическое поле Eenv, в котором находится ШМ, пропорционально радиусу rFB при данной концентрации электронов пе. Из этого следует, что, если типичная концентрация (согласно примерам в [7, 8] при радиусе ШМ rFB~5 см) составляет пе ~ (1015 v 1016)см-3, удельная энергия в ШМ с объёмом VFM (при полной электростатической энергии Eelst )
Espec = EelstlVFB ~\e\neEenvrFB К nerFB ,
то есть при увеличении диаметра ШМ до 1,5 м и соответствующем значительном увеличении мощности излучения w предполагаемая удельная только электростатическая энергия разделения зарядов сама по себе может составить Espec > 1010 Дж/м3 (что можно сопоставить с примером в [8]
ESpeC > 108 Д^/м3). Это не противоречит наблюдениям Espec ~ 3,4х1010 Дж/м3 [12, 13] и
Espec ~ 3х1012 Д^/м3 [14], что представляется, однако, крайне завышенным*.
С учетом не столь огромного повышения реально наблюдаемой светимости w при больших размерах ШМ, следует предположить, что при сохранении электронных плотностей в интервале
*
Сравнение: например, удельная теплота сгорания бензина составляет ~ 3,1 х 1010 Дж/м3.
61
n
(FB
e
~ (1015 -1016)
см3 должна существенно уменьшиться температура (Te > 8000K и ниже) с сильным
уменьшением мощности излучения ~ zT4rlB и соответствующего ей параметра неравновесности а»1, n(FBB = ne/а. Следует также учесть очень сильную реабсорбцию излучения при размерах ШМ ~ 100 см и более. Если W0 - равновесная мощность излучения, соответствующая температуре Те в отсутствие реабсорбции, и Wr - мощность излучения при наличии реабсорбции, то при замене ne (Te) на n (Te )/ а можно получить соотношение, учитывающее реабсорбцию, в виде
( w А
w =-
W
V W0 у
_2
Значения W0 определяются линейной экстраполяцией имеющихся расчётных данных “кривой роста”
2
а
от очень малых толщин на большие оптические толщины.
Гипотетически нельзя исключать, однако, наличия двух ветвей ШМ: низкотемпературной с малыми а > 1 и высокотемпературной ветви с большими а >> 1. Открытие двух ветвей ШМ могло бы стать интригующим результатом, характеризующим многообразие форм молнии.
Отметим, что согласно (1) большим размерам ШМ соответствуют очень большие напряжённости поля E , обусловливающие воздействие на одиночные заряды молнии, преобладающее над
кулоновским взаимодействием соседних зарядов.
Эти соображения прокладывают путь к расчёту опций, характеризующих шаровую молнию больших размеров, включая трудоёмкие расчёты условий устойчивости.
ЗАКЛЮЧЕНИЕ
Предложенная ДДМ [6-8], по-видимому, является первой моделью этого рода и имеет эвристический характер. Отмечается принципиальная возможность расчёта основных параметров шаровой молнии, но полные расчёты в настоящее время вряд ли возможны из-за в высшей степени комплексного характера проблемы: необходимости знания параметров основных процессов, крайне кропотливых отсутствующих вычислений параметров неравновесной плазмы с необходимостью знания многих сечений соотносящихся друг с другом столкновительных процессов и коллективного установления степени черноты плазменного объёма. Эмпирически могут быть рассмотрены различные опции, не противоречащие субъективным оценкам наблюдаемых температур Tm, Te и светимостей молнии w. Самое главное: с этими версиями температур и светимостей ШМ может быть устойчивой только в случае термодинамически сильно неравновесной плазмы и, по крайней мере, при температурах Tm несколько тысяч кельвинов и Те > (6 -12) x103K. В дополнение в литературе приводятся наблюдения ШМ с очень большими энергиями и плотностями энергии до »1010 Дж/м3. В то же время
расчёты в рамках оценок предложенной приближённой теории по порядку величины по вышеуказанным формулам:
YEenv
п me | | n ;
о \e\rFB 2 ;
3 mav а nm
YEenv = 1 к^/см
а2 = W/ w ;
Espec = N YEenv ' 2F ' nelа
показывают, что, например, при ограничении Eenv напряжением стримерного пробоя
(1 < Eem < 10) кВ/cм и Tm ~ 1000 K, Te ~ 10 000 K, rFB = 100 cм получаем в ДДМ при jEern = 1 кВ/см
< 1,9x107ynjTFB = 5,1x1015 см-3;
Espec <6,0x10-2VvFB~ 1,6x109 Дж/м3 .
Величины nBeFBB = njа ~ xe[11]- nm/а, а2 = W/w и Espec ограничиваются как максимально допустимой
напряжённостью атмосферного электрического поля, так и приемлемыми из наблюдений значениями мощности излучения шаровой молнии w и условиями устойчивости. Результаты оценок, представленные в таблице, показывают ориентировочную максимальную плотность энергии ШМ при боль-
62
ших размерах порядка Espec < 109 Дж/м3 . При этом скорость движения молнии определяется отношением vFB = w/f , f = |e| yEenvNe и может быть крайне малой, менее (1 -10) см/с. Примеры быстро перемещающихся ШМ малых размеров с малым энергосодержанием, соответствующие малым у, температурам Te <10 000 K, приведены в табл. 5 и 6 работы [8] и в [7], а также в нижеприведённых резуль-
татах расчёта параметров типичной низкоэнергетической малой ШМ (n^fB') = njа):
е ~ 3,0 х10-
а=18,7
\е\ lE,nv = 1,09 х10 11 дн Vfb =1,62 м/с
Tm = 3000 K у = 8,3х10-3 nj а = 1,09 х1015 см-3 e2 (ne/а/3 = 2,43х10-9 дн wjSFB = 0,32 Вт/см2
Те = 9000 K nm = 2,44 х1018 см-3 yEenv = 6,8 В/см 1/у £ 2,2 х102 Epec =1,2х104 Дж/м3
w = 100 Вт ne = 2,0х1016 см-3 Ffb = 6,2 кгс Eenv > 1,5 кВ/см Efull = 5,6 Дж
Для скорости vFB существует простое соотношение (без учёта потерь энергии на создание передней струи отдачи):
3 есГ
m„,
2 2 2 п а е r„R
„( FB) У
Vfb =5 см
П
m
Кажущийся парадокс малых скоростей может разрешаться тем, что в случае небольших ШМ с наблюдаемой кажущейся постоянной скоростью движения можно было пренебречь силами инерции.
Крупные ШМ могут двигаться с ускорением в соответствии с выражением баланса энергии dE/dt
d_
dt
( Fs ) = FvFB + s
dF
---= w
dt
(t) + Wgd (t) + MFBVFB
a
acc
vL dM „
2 dt
+
(12)
= J aacc (t) dt; F = N YEenvNe '; Ne = ^Kb , S = J vFBdt ,
(13)
где MFB - масса ШМ; aacc - ускорение; s - проходимое ШМ расстояние; Wdd (t) - мощность, затрачи
v
0
0
ваемая на газодинамическое сопротивление ШМ, включая мощность передней струи. Так что можно получить некоторое ускорение при замене в уравнениях (12) или (13) w(t) + Wgd (t)-FvFB , например на
+£FvFB , 4 ^ 1. При этом трудно представить существование предельных режимов ШМ, аналогичных представленным в таблице, из-за кажущихся нереальными больших значений силы F электростатического разделения зарядов и существенно возрастающей роли сил инерции, что означает необходимость учёта ранее введённой поляризационной поправки 0,01 < у < 1, приводящей к уменьшению
силы F ^ F = у|e\EenvNe, а также необходимость дальнейшего снижения предельных плотностей энергии ШМ при крайней предполагаемой неустойчивости таких режимов.
Увеличение скорости шаровой молнии необходимо для компенсации дополнительных потерь энергии на ускорение ШМ. Это означает, что существует необходимость дополнительных неизвестных выражений для временных взаимосвязанных зависимостей n^eFB'} (t), rFB (t), w (t), F (t), которые
могут быть получены из базового уравнения (1-2) и условий устойчивости в параметрическом пространстве. Это означает также, что все ШМ могут двигаться с ускорением, по крайней мере на начальном участке, наличие которого определяет возникновение неустойчивостей. Ускорение может быть признаком неустойчивости и являться возможной внутренней причиной, помимо атмосферного окружения и изменений Eenv, определяющей время жизни ШМ.
Большие удельные электростатические энергии возможны при разумном увеличении яркости ШМ до одного и более Вт/см2. При наличии огромных молний (как с малыми, так и с большими удельными энергиями) время жизни очень больших ШМ с большой мощностью излучения может быть ограничено быстрым исчерпанием запасов окружающего атмосферного электричества.
Отметим, что энергия ШМ определяется балансом энергии, поступающей за счёт ее движения в атмосферном электрическом поле, и потерь энергии главным образом на излучение.
Типичная ШМ экспериментально не получена. Это привело к огромному количеству зачастую крайне экзотических моделей этой молнии с более или менее правдоподобными, но чаще весьма сомнительными допущениями. Поэтому единственным критерием истинности модели ШМ, который можно было бы до сих пор использовать, является естественность этих допущений (например, нерав-
63
новесные температуры T ^ Te обычного электрического разряда в разреженном газе и наблюдения “холодных” ШМ) и качественное соответствие модели многочисленным разрозненным наблюдениям.
Шаровая молния не может быть создана экспериментально - как неподвижное (статичное) плазменное образование.
Можно было бы сначала создать сильно неравновесный сфокусированный СВЧ-разряд и протяжённое очень сильное постоянное электрическое поле, “вытягивающее” нарождающуюся самоорганизующуюся ШМ, что являются само по себе сложной инженерной задачей. При экспериментальной реконструкции молнии основные усилия должны быть направлены на создание крупномасштабной системы направляющих электродов, воспроизводящих очень сильное атмосферное электрическое поле, способное привести в движение и обеспечить энергией объёмный затравочный разряд.
Признательность. Выражаю искреннюю признательность В.Л. Бычкову, руководителю Отделения физической электроники Московского государственного университета, за ряд критических замечаний и рекомендованную литературу. Особую признательность и благодарность выражаю профессору Ярославского государственного университета А.И. Григорьеву за действенную поддержку и ценные замечания.
ЛИТЕРАТУРА
1. Bychkov V.L., Nikitin A.I., Dijkhius G.C. Ball Lightning Investigations. Physics of Earth and Space Environment. Springer-Verlag, 2010. 373 pp.
2. Никитин А.И. Высокоэнергетичные шаровые молнии - путь к решению проблемы. Труды IXМеждународной конференции “Волновая электрогидродинамика проводящей жидкости. “Долгоживущие плазменные образования и малоизученные формы естественных электрических разрядов в атмосфере”, 30.06-04.07 2011, Ярославль, Россия. С. 138-154.
3. Никитин А.И. Удастся ли решить проблему шаровой молнии в 21-м веке? Химическая физика. 2006, 25(3), 18-37.
4. Норман Г.Э. Шаровая молния как переохлаждённая неидеальная плазма. Плазма ХХ век. Тезисы конференции ФНТП-98. Петрозаводский государственный университет, Россия, 1998.
5. Dvornikov M. Formation of Bound States of Electrons in Spherically Symmetric Oscillations of Plasma. Phys. Scr. 2010, 81, 055502.
6. Soshnikov V.N. Dipole Dynamical Model of Ball Ligtning. Proceed. IX International Conference ”Wave Electrohydrodynamics of Conducting Fluid. Long-lived Plasma Formations and Little-studied Natural Forms of Electrical Discharges in the Atmosphere”. 30.06-04.07 2011, Yaroslavl, Russia, pp. 203-205.
7. Soshnikov V.N. Comments to Support of the Dipole Model of Ball Lightning. Preprint ArXiv.org/physics.gen-ph/arXiv: 1007.4377 (2010).
8. Soshnikov V.N. Purely Electrical Nature of Ball Lightning (BL), its Elementary Equations, Calculated Parameters and Conditions of Possible BL Experimental Generation. Preprint ArXiv.org/physics/plasma physics/arXiv:0903.4977 (2009).
9. Франк-Каменецкий Д.А. Лекции по физике плазмы. М.: Атомиздат, 1964. С. 199-204.
10. Авилова И.В., Биберман Л.М., Воробьёв В.С. и др. Оптические свойства горячего воздуха. М.: Наука, 1970. 320 с.
11. Предводителев А.С. и др. Таблицы термодинамических функций воздуха (для температур от 6000 К до 12 000 К) и давлений от 0,001 до 1000 атм. М.: Изд-во АН СССР, 1957.
12. Egely G. Analysis of Hungarian Ball Lightning Observations. in. Keul A.G (ed.) Progress in Ball Lightning Research: Proc. of the VIZOTUM. Salzburg, Austria, 1993, p.22.
13. Григорьев А.И. Шаровая молния. Ярославль: ЯрГУ, 2006. 200 с.
14. Bychkov V.L. On Observational Properties of Ball Lightning. Proc. of the 9th ISBL 06. Eindhoven. Netherlands, 2006, pp. 18-25.
15. Finkelstein D., Rubinstein J. Ball Lightning. Phys. Rev. 1964, 135, A390-A396.
16. Powell J.R., Finkelstein D. Ball Lightning. American Scientist. 1970, 58(8), 262-280.
Поступила 23.03.12
Summary
We present the first dipole dynamic model (DDM) of ball lightning (BL), which explains in a is simple, natural and rational manner most of the observations of BL, with the possibility of very high energy densities provided by dipole charge separation in BL and with the inflow of energy of atmospheric electricity at BL translational movement. The virial theorem that can limit BL energy does not apply to BL, which is not a closed system like the Sun and Galaxy systems and is a part of a large system extended in time and space strongly dependent on it.
64
