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

НОВЫЕ ТЕОРИИ

электростатический механизм образования радиоактивных облаков

Ю.Л. Ратис, И.А. Селезнева Самарский государственный аэрокосмический университет имени академика С.П. Королева Институт систем обработки изображений РАН

Аннотация

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

Введение

В цикле работ, суммированных в монографии [1], была развита радиоизотопная теория шаровой молнии (ШМ). Логика происхождение гипотезы о радиоизотопной природе ШМ такова:

1. Шаровая молния всегда оставляет после себя запах серы, озона и окислов азота.

2. Сера может образоваться только в результате в-распада фосфора.

Постоянная скорости в-распада сильно зависит от степени ионизации распадающегося радио-

дп

3.

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

4. Радиофосфор является распространенным в природе элементом. Он обнаружен в дождевой воде в макроскопических количествах.

Согласно [1] феноменологическая модель шаровой молнии, описывается системой кинетических уравнений (1):

—1 + У(пі^) = -Х1п1 +Уп1 +і3П2 -Р!2ПіП2 + Р22П2 +V ді

X ВііУпі-Пі X ^УП;

V і=1

дп „ „

—2 + У(;П2:У2) = -ХоП2 +Р 21П2П1 -в 22П2 +У ді

X °2іУПі-П2 X ^2;УП;

і=і

і=і

і=і

+ 42

+ 4і

(1)

дпз

д;

+ V(n3V3) = Х1п1 +(Х1 + Х2)п2 +У

X Вз;уП;-Пз X 8з;уП;

V і=1

где п1 - плотность невозбужденных атомов

32 р І?" 1

(“Р), п2 - число возбужденных атомов 15Р (ЦР), с Ь - оболочки которых удален один (или более)

2.

V, -

электрон, а п3 - плотность паров атомов серы скорости макроскопических течений соответствующих компонент смеси. В нашем случае

V я у2 я у3 я уа1г, где уа1г - скорость ветра.

Кроме того, в системе уравнений (1) использованы следующие обозначения: q1 - интенсивность

внешнего источника атомов радиофосфора (шаровая молния находится во внешней среде, содержащей

1 Всюду далее, где говорится о парах радиофосфора, подразумевается, что мы имеем дело либо с пятиокисью фосфора Р205, либо с фосфорной кислотой Н3Р04. Кроме того, термин

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

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

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

і=і

эти атомы, и захватывает их во время полета); 42 -интенсивность ионизации К - и Ь -оболочек атомов радиофосфора за счет механизмов, отличных от резонансной фотоионизации (например, столкновений с энергичными электронами, эффекта Пеннинга и т.п.); Х1 - постоянная р-распада в непрерывный спектр; X2 - постоянная индуцированного р-распада в связанное состояние (в дискретный спектр, на вакансию в Ь -оболочке); X3 - вероятность спонтанного перехода (в единицу времени) электрона на Ь -оболочку возбужденного (либо ионизованного) атома I2" с М,^... - оболочек или из непрерывного спектра; Р12 - константа скорости фотовозбуждения и фотоионизации атомов 1352 Р («обдирания» Ь -

оболочки атома 3? Р ) у -квантами, образующимися

при индуцированном р-распаде радиофосфора; Р21 - перекрестный коэффициент (очевидно, что Р12 =Р21); коэффициент Р22 - константа скорости выбытия возбужденных атомов фосфора из-за взаимодействия между собой.

В системе (1) учтены потоки паров радиофосфора, «ионов» радиофосфора и серы ^. В общем случае эти потоки можно представить в виде суммы линейной и нелинейной компонент:

3Ь =-2 »|]Уп],

\

3 Г =2х цА ^

(2)

(3)

].к

где Б. - коэффициенты линейной диффузии компонент, а Хцк - тензор третьего ранга, причем линейная комбинация БГЬ =-^%,.кпк

имеет смысл

нелинейного коэффициента диффузии [2].

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

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

Основной формализм

Облака паров радиофосфора образуются в верхних слоях атмосферы ядерных реакций

2481(а,р)12Р

1268(а,2р)12р

11р(й,р)12р

328(р,2р)125Р:

15Р(п, у )12р

12* ^ 17^*4“*^,/ 12*

”8(п,р)”Р , 10Аг(р,4Ве)12Р и 2239Си(р,12Р)12Р , а также 12Р(п,у^Р, 118(п,р)11Р, 17С1(у,2р)ЦР,

”С1(г,а)”Р , 18Аг(р,2а)12Р и ^(р,^)^ .

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

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

3

нимает вид :

дг+^у 1)=

(4)

= -^1п1 +Уп2 + ВИУ 2п1 - 8цУ(п1Уп1) Рассчитаем коэффициенты у и 811 в электростатическом приближении исходя из того, что из-за интенсивного выброса электронов при в-распаде в непрерывный спектр радиоактивное облако является сильно заряженным. За 1 секунду оно выбрасывает

из себя отрицательный заряд примерно Дq~10-7 -10-2 Кл (оценки приводятся для облака суммарной активностью 100-200 кюри, то есть, имеющего параметры, превышающие критические значения). Вследствие этого электрический потенциал поверхности облака за 1 секунду возрастает примерно на

1 Д(| 9 • 109 ■ 10-7 101 „

Дф = --------------------------В.

4ле0 г г г

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

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

(Н3Р04)- и пятиокиси фосфора (Р202) , а также

нейтральные полярные молекулы этих соединений.

Предположим, что центральная часть облака паров радиофосфора имеет заряд О = еГ0, где е - заряд электрона. Тогда уравнение движения ионов (Н3Р04)- и (Р202) имеет вид:

н2Г е2! Г —Г

М—г = -кехр(-цг)—г°— — (5)

— гг

где Г - радиус-вектор заряженной частицы (считается, что центр облака находится в начале координат), М - эффективная масса ионов4, п - коэффициент вязкости воздуха, р - его плотность, к— -эффективный радиус заряженной капельки, ц - постоянная экранировки, по порядку величины равная Ц -I-1 (1е - длина свободного пробега в воздухе электронов с энергией Ее ~ 1 МэВ). Учитывая, что силы жидкого трения пропорциональны скорости, мы легко приходим к выражению для установившейся скорости заряженных капель (у—), содержащих радиофосфор:

к е2Г Г

-ехр(-цг) 0

г2 Г

(6)

2ппК—

В облаке, содержащем Г1 атомов радиофосфора, происходит Х1Г1 распадов в секунду, то есть формально заряд нарастает по закону О = еА 1Г1. Однако рост заряда имеет место только на начальной стадии формирования облака. На самом деле положительный заряд облака в нулевом приближении

к

3 В настоящей работе мы пренебрегаем гравитационными поправками в уравнении (5.1) и далее, поскольку пары соединений радиофосфора падают на Землю под действием силы тяжести очень медленно.

4 Ионы являются центрами конденсации водяных паров, по-

этому на высотах менее 10 км имеет место соотношение

М >> т , где т - масса отдельной молекулы соединений радио-

фосфора

можно считать пропорциональным величине А1Г1, поскольку его прирост постоянно компенсируется за счет притока свободных электронов и отрицательных ионов из окружающего пространства. По этой причине логично считать, что О = тА1еГ1, где Т -время релаксации. Учтем тот факт, что заряд центральной части облака распределен в пространстве, и представим уравнение (6) в следующем виде:

кА1т е21

У— =-----;— е2 Г —г'п, X

— 2ппЯ— V 1

ехр(-ц |г - г'|)? - ?' , <(?')—;—,2—!- ,

г - г' Г - Г

(7)

причем интегрирование ведется по объему облака V.

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

1 = п1 • V — , (8)

поэтому число ядер радиофосфора в объеме облака V нарастает по закону:

дГ1

д1

= -[аз'пДГ')? —(Г').

(9)

Подставляя (7) в (9), получаем:

2 Г —гп1 (Г) X

дГ. ктА. 2

—— =---------1—е

д1 2лпК—

( ,,, (Г') ехр(-Ц|?'- Г) ?'- г <1 аз'пДг')-------------'-'----

(10)

Г - Г

Введем обозначение Ь : следующие преобразования: д]Г,

г - г

ктА

2япК

1 е2 и проделаем

с л е

1 = - ь ] аГп1 (?) Г аг'п1 (Г ' )—

Г ■ - ? I? - ?|

-ц?

= - ь Г агп (? )| а? [1 (?+г ')]—- 5? ' =

V 8 (г)

= -ь Г ат 1 (? )| а? V'| [1 (? + ? ' )]—) ? ' | =

V V I (г') ]

= - ь Г а?п1 (? )| а? 'п1 (Г + ? ' ^ ' \ ^? ' I - (11)

V V 1(г') ]

-цг'

-ЬГ —гп1 (г)Г а?' (—( • V 'п1 (г + ?')) =

V V (г)

= цЬГ а?п1 (Г)Г а? п (? + ? ') ехр(~>ЦГ') -

V V (?')

-ьГ а?п1 (?)Г а?' ехр(-ЦГ ) (• V 'п1 (? + ? ')).

V V (г')

С другой стороны, согласно (4) прирост количества радиофосфора в объеме облака равен:

дГ г с

я а?п1 (?) -811Г аrv(l(Г)Vnl(Г)). (12)

Формально приравнивая эти выражения (имеющие различную степень точности и разные границы области применимости), получаем:

Ьц| аг^ (?) | а?'п1 (Г+?' )ех|(^|)|ЦГ )-

V V

-ь|а?п1(г)|а?' ехр^>)|ЦГ') ( •v'nl(f+?')) я (13)

V V

я а?п2 (Г) -811Г а^(п1 (Г )Vnl (?)).

V V

Разложим подынтегральные выражения в (13) в ряд по Г' и используем тот факт, что экранированный кулоновский потенциал имеет конечный радиус действия. Тогда для первого слагаемого в левой части соотношения (13) имеем:

Г а?^ (?) Г а?'^ (?+?')-^-г я (г')

V V

■\ агпд? )| аг'пд?) (-

V V

(г')

(14)

+Г а?п1 (?) Г а?' ( vnl (?))-^-2.

(г')

V V

В силу сферической симметрии задачи второе слагаемое в (14) зануляется:

Г а?'(г' vnl(Г ))ехр^,|ЦГ') - 0, V (г)

(15)

а интеграл в первом слагаемом вычисляется элементарно:

Г а?

ехр(-цг') = 4л

(г') Ц

еме

у я Ьц| а?'

(16)

Отсюда немедленно следует, что ,' ехр(-ЦГ') =

: 4пЬ :

(••'Г

2ктА1е2

3ПКа

(17)

Совершенно аналогично находится величина 811 : 811Г —^(пдг )Vnl(Г ))я

V

я ь Г а?п1 (? )| а?' ехр(-ЦГ ) ( • V' п1 (?+?')) я

V V

(?')

’> ь Г а?п1 (? )Г а? 'ехр( 11? );

■* ^ (г ' Г1

V V

с(г' • [((? ) + г ' V 2п1( ?)]) =

= ь I „г^ (г )| „? 'ехр( ,цг ) V 2п1 (г ) =

V V г

= ь| „гпд? )у2п1(? )| а?' ехр(-цг ) = (18)

4пЬ

Ц V

I „ГпД? )У 2Пі(? ),

и в рамках точности сделанных приближений мы приходим к выражению:

8 =^

а11 2 .

Ц

(19)

Таким образом, нелинейные члены уп1 и 8llV(nlVnl) в приближении малых градиентов

концентрации описывают процесс самоорганизации облака, происходящей за счет притяжения отрицательных ионов молекул (Р202) и (Ы3Р04) к центральной части облака химических соединений радиофосфора, имеющей большой положительный электрический заряд из-за потери электронов, вылетающих из него при в -распаде в непрерывный

спектр. Именно эти электроны с энергией Ее ~ 1 МэВ ионизируют окружающий воздух. В результате в воздухе появляется большое количество свободных электронов, которые захватываются нейтральными молекулами Р202, и Н3Р04, которые с огромной силой притягиваются к центральной положительно заряженной части облака. Таков микроскопический механизм самоорганизации паров радиофосфора в облака и затем в шаровые молнии.

Проанализируем полученные выражения.

Из соотношения (17) следует, что коэффициент у зависит от вязкости воздуха и времени электрической релаксации положительного заряда центральной части облака паров радиофосфора и не зависит от параметра ц -I-1, обратно пропорционального длине свободного пробега электронов.

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

над уровнем моря как т / п .

Коэффициент вязкости п согласно [3] равен:

п = 1nmXv. 3

(20)

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

Произведение птХ в формуле (20) от высоты над уровнем моря практически не зависит, поскольку определяется, в основном, сечением рассеяния ионов серы на молекулах воздуха. Средняя скорость

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

Проанализируем свойства коэффициента 811 . Для этого представим выражение (19) в следующем виде:

•у 12.

(21)

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

Наконец, определенный интерес представляет анализ влияния нелокальности реального кулонов-ского взаимодействия.

Согласно (4) и (11) в рамках сделанных предположений плотность паров химических соединений радиофосфора подчиняется следующему интегро-дифференциальному уравнению:

5п1

■ + V(n1 V1) = -Х1п1 + В11V п1 +

+Ьп,(г )| „Г еХр*=цГ> х

V (г )

х [цПі (г + г ') - (• V 'Пі (г + г '))] ,

(22)

где п' = ?' / г' - единичный радиус-вектор.

На очень больших высотах уравнение (22) должно быть дополнено внешним источником:

5п1

"дГ

+ V(n1V1) = -Х1п1 +D11V2n1 +

+ьп,(г )| „г х

V (г )

х [ЦПі (г + г ') - (• V 'Пі (г + г ')) + ^ ,

(23)

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

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

5п1

"дТ

Чп1 + Ьп1 (Г)п1 + Чі

(24)

где п1 - средняя концентрация паров радиофосфора. Решение этого уравнения тривиально:

n1 (t) = exp ((bn1 — X1 )t) x

t

Пі (0) + J d£ exp (—(bnl — Xl )^) (Jl ©

(25)

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

Поскольку Ь = ^/А1 / Б11 , то из (32) с очевидностью следует, что ускорение в -процессов за счет реакции индуцированного в -распада в связанное состояние и вызванное этим ускорение процессов ионизации приводят к быстрому самосжатию облака радиофосфора. Именно таков механизм формирования тела шаровой молнии на стадии, предшествующей ее «самовозгоранию».

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

Полученное в [1] соотношение уп1"“ =А1 в совокупности с соотношением (17) позволяет оценить критическую концентрацию паров радиофосфора, начиная с которой зародыш радиоактивного облака превращается в самоорганизующийся объект.

3nRd

2kTe2

(26)

Дадим численную оценку этой величины, опираясь на предположение о том, что процесс формирования облака радиофосфора происходит в тропосфере, и с определенной натяжкой можно использовать данные о вязкости воздуха при нормальных условиях (п = 1,8 10-2 кг-с-1-м-1). Полагая

Ка~10-9 м (мельчайшая капелька воды, сконденсировавшаяся на ионе (Н3Р04) ), т ~1 секунды, и

отбрасывая все несущественные факторы типа 3/2, приходим к оценке величины:

10—2 • 10—

1010 • 10—

(27)

что не противоречит другим оценкам и здравому смыслу.

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

Заключение

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

1. Дана оценка величины коэффициента у и коэффициента нелинейной диффузии Sjj .

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

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

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

Благодарности Работа выполнена при поддержке Министерства образования и науки РФ, правительства Самарской области и Американского фонда гражданских исследований и развития (CRDF Project SA-014-02) в рамках российско-американской программы «Фундаментальные исследования и высшее образование» (BRHE).

Авторы выражают особую благодарность ректорату Самарского государственного аэрокосмического университета, а также Фонду содействия экономическому развитию СГАУ (Фонду Лукачева) за финансовую поддержку.

Авторы благодарят В.А. Сойфера, Н.Л. Казанского, В.И. Фурмана, И.П. Завершинского, В.С. Казакевича, В.В. Ленивкина, Л.В. Туманова за многочисленные стимулирующие дискуссии, в результате которых на свет появилась эта работа.

Литература

1. Ратис Ю.Л. Шаровая молния как макроскопическое квантовое явление // Монография, СГАУ, ИСОИ РАН, Изд. СНЦ РАН. Самара, 2004. 132 с.

2. Таланов В.И. Стимулированная диффузия и кооперативные эффекты в распределенных кинетических системах // Нелинейные волны. Самоорганизация. Сборник трудов ИПФ АН СССР, М.: Наука, 1983. 264 с.

3. Сивухин Д.В. Общий курс физики, Т.2, Термодинамика и молекулярная физика // М.: Наука, 1975. 551 с.

n