Моделирование электрических атмосферных разрядов высокой мощности с помощью ускорителей протонов средних энергий

Язык труда и переводы:
УДК:
537.52
Дата публикации:
27 ноября 2022, 03:56
Категория:
Перспективные направления исследования необратимых физических процессов
Авторы
Смыслов Денис Юрьевич
МГТУ им. Н.Э. Баумана
Горев Владимир Васильевич
НИЦ «Курчатовский институт» — ИРЕА
Аннотация:
Молния остается до сих пор до конца не изученным явлением не только из-за сложности ее регистрации в природе, но и в силу невозможности создания в лабораторных условиях полноценного аналога происходящих в ней процессов. Такая ситуация требует принципиально новых подходов к изучению электрического атмосферного разряда. Впервые предложен метод, объединяющий естественные и лабораторные условия проведения эксперимента. Показано, что такой подход может устранить существующие недостатки в современных методах исследования. Проведены количественные оценки параметров, используемых в предложенной схеме генерации электрических атмосферных разрядов с помощью ускорителей протонов средних энергий. Результаты расчета показали практическую применимость рассмотренного метода.
Ключевые слова:
электрический атмосферный разряд, ускорители протонов средних энергий, двухфазная система, водность облаков, интенсивность осадков, многократное рассеяние заряженных частиц, время зарядки
Основной текст труда

Введение

На сегодняшний день не существует теории, полностью описывающей явление молнии. Сложность ее изучения в природе обусловлена хаотическим разбросом областей поражения. На каждый квадратный километр поверхности земли молния ударяет в среднем 2–4 раза в год на территории Европы. Высотные же объекты поражаются преимущественно восходящими молниями, оставляя затруднительным исследование нисходящего разряда. К примеру, Останкинская башня (560 м) поражается молниями 25–30 раз за год, среди которых только 2–3 оказываются нисходящими.

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

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

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

На рис. 1 показана примерная схема предлагаемого эксперимента. Лабораторной установкой является линейный ускоритель протонов средних энергий с отклоняющими магнитами на выходе, позволяющими направлять пучок заряженных частиц в нужную точку пространства. В момент приближения дождевого облака, созданного естественным или искусственным путем (распылением с самолета), ускоритель генерирует пучок высокоэнергетических протонов, которые, проходя через атмосферу, попадают в облако под некоторым рассеянным углом \theta и остаются в нем, распределяясь на поверхности капель воды. Облучение происходит в течении времени, необходимого для достижения электрического поля с напряженностью E\geq 30{\text{ кВ/см}} , при которой возникает пробой в сухом воздухе в промежутке облако – земля или облако – облако. Проведем оценку параметров дождевого облака, используемого в описанной схеме.

 

Оценка размеров облака

1.  Рассмотрим для простоты расчетов сферическое облако радиусом r, расположенное на высоте h = 500 м над землей. Пусть протоны с энергией E= 400 МэВ вылетают из ускорителя под углом 60º. Пробег протона в воздухе при такой энергии составляет tвоз = 93,73 г/см2 [2]. Переходя к единицам длины в метрах:

R_{\text{воз}}={\frac {t_{\text{воз}}}{\rho _{\text{воз}}}}=726{\text{ м}},\qquad                                                                 (1)

где \rho _{\text{воз}}=1,29\cdot 10^{-3}{\text{г/см}}^{3} — плотность воздуха.

Так как протоны должны остаться в облаке, R_{\text{воз}}=L+2r,

где L — часть пути, пробегаемой протоном в атмосфере; D = 2r — часть пути, приходящаяся на диаметр рассматриваемого облака.

Параметры L иr можно выразить из геометрических характеристик:

\sin 60^{\circ }={\frac {500+r}{L+r}}.

Решив систему из этих двух условий, получим оценку параметров L и r для пробега протона в воздухе:

{\begin{cases}L=588{\text{ м;}}\\r=69{\text{ м.}}\end{cases}}                                                                                                    (2)

2.  Для определения длины пробега протона в облаке рассмотрим двухфазную систему, состоящую из воздуха и капелек воды. Плотность такой системы выражается формулой [3]:

\rho _{\text{обл}}=\rho _{\text{вод}}e+\rho _{\text{воз}}(1-e),                                                                             (3)

где e=n\cdot V_{\text{к}} — объемная доля капель воды в облаке:

e={\frac {n\pi d^{3}}{6}};                                                                                                           (4)

\rho _{\text{вод}}=1{\text{ г/см}}^{3};  Vк — объем сферической капли; n — концентрация капель в облаке; d — диаметр капли.

В метеорологии используются понятия интенсивности осадков и водности облаков.

Водность облака — это масса сконденсированной воды mвод приходящаяся на единичный объем воздуха Vобл [4]:

Q={\frac {m_{\text{вод}}}{V_{\text{обл}}}}={\frac {n\cdot m_{\text{вод}}}{N_{\text{к}}}}=n\rho _{\text{вод}}V_{\text{к}}={\frac {n\rho _{\text{вод}}\pi d^{3}}{6}},                                        (5)

где N_{\text{к}} — количество капель воды в облаке.

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

I={\frac {\pi d^{3}}{6}}(v-u),                                                                            (6)

где v — скорость падающей капли, м/с; u — скорость восходящих потоков, м/с.

Из (5) и (6) получим связь водности с интенсивностью:

Q={\frac {I\rho _{\text{вод}}}{v-u}}.                                                                                       (7)

Подставив (4) в (3) с учетом (5), получим формулу, в которой плотность облака будет зависеть только от его водности:

\rho _{\text{обл}}=\rho _{\text{воз}}+Q\left(1-{\frac {\rho _{\text{воз}}}{\rho _{\text{вод}}}}\right).                                                             (8)

Заряженные частицы, попадая в облако, будут скапливаться на поверхности капель воды, поэтому для наиболее эффективной зарядки необходимо, чтобы оно обладало максимальной водностью. В работах [4, 5] показано, что I\sim u и d\sim u , а v\sim d . Возьмем для расчета плотности облака максимальные значения параметров, при скорости восходящего потокаu = 1 м/с, что довольно часто встречается в природе [5]:

  • I = 8 мм/мин = 1,3 \cdot 10–4 м/с;
  • d = 4 мм = 4 \cdot 10–3 м;
  • v = 8,83 м/с;
  • \rho _{\text{воз}}=1,29\cdot 10^{3}{\text{ г/м}}^{3};
  • \rho _{\text{вод}}=10^{6}{\text{г/м}}^{3}.

Тогда максимальная водность облака

Q={\frac {I\rho _{\text{вод}}}{v-u}}={\frac {1,3\cdot 10^{-4}\cdot 10^{6}}{7,83}}=16,6{\text{ г/м}}^{3}.                              (9)

Плотность облака

\rho _{\text{обл}}=\rho _{\text{воз}}+Q\left(1-{\frac {\rho _{\text{воз}}}{\rho _{\text{вод}}}}\right)=1,29\cdot 10^{3}+16,6\cdot 0,9987=1.31\cdot 10^{3}{\text{ г/м}}^{3}.                                       (10)

3. Поскольку на общую длину пробега протона в воздухе Rвоз приходится k диаметров облака D:

R_{\text{воз}}=kD\Rightarrow k={\frac {R_{\text{воз}}}{D}}={\frac {L+2r}{2r}}=5,26.                                         (11)

Величина пробега протона в воздухе t в г/см2 (1), приходящаяся на долю k, составляет

{\frac {\rho _{\text{воз}}R_{\text{воз}}}{k}}={\frac {93,73}{5,26}}=17,82{\text{ г/см}}^{2}.                                       (12)

Тогда длина пробега протона в облаке

\rho _{\text{обл}}R_{\text{обл}}=17,82\Rightarrow R_{\text{обл}}={\frac {17,82}{1,31\cdot 10^{-3}}}=136{\text{ м}}.                                                                (13)

Так как цель задачи состоит в зарядке облака, то его диаметр D не должен быть меньше длины пробега протона в нем:

D\geq R_{\text{обл}}.                                                              (14)

Многократное рассеяние протонов

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

1. В работе [6] средний квадрата угла многократного рассеяния заряженных частиц выражается формулой

 

\langle \theta ^{2}\rangle =0.157\left({\frac {Z(Z+1)}{A}}\right)\left({\frac {z^{2}t}{(pv)^{2}}}\right)\ln \left[1,13\cdot 10^{4}{\frac {Z^{\frac {4}{3}}z^{2}t}{A\beta ^{2}}}\right].                                      (15)

Запишем необходимые параметры для вычисления среднеквадратичного угла многократного рассеяния:

  • заряд протона z = 1;
  • заряд рассеивающей частицы — азота Z = 7 (составляет 80 % атмосферы)
  • атомная масса рассеивающей частиц A = 14 а.е.м.
  • расстояние, пробегаемое протоном в воздухе t=L\rho _{\text{воз}}=588\cdot 10^{2}\cdot 1,29\cdot 10^{-3}=75,85{\text{ г/см}}^{2};
  • кинетическая энергия протона E = 400 МэВ
  • масса протона mp = 1,67 \cdot 10–27кг;
  • скорость света c = 2,99 \cdot 108 м/с;
  • \beta ={\frac {v}{c}}={\sqrt {1-{\frac {1}{\left(1+{\frac {E}{m_{p}c^{2}}}\right)^{2}}}}}={\sqrt {1-{\frac {1}{1,43^{2}}}}}=0,72;
  • pv={\frac {mv^{2}}{\sqrt {1-\beta ^{2}}}}={\frac {2E}{\sqrt {1-\beta ^{2}}}}={\frac {800}{\sqrt {1-(0,72)^{2}}}}=1152,78{\text{ МэВ}};

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

{\sqrt {\langle \theta ^{2}\rangle }}={\sqrt {0.157\cdot {\frac {56}{14}}\cdot {\frac {75,85}{1328901,73}}ln\left[1,13\cdot 10^{4}{\frac {13,39\cdot 75,85}{14\cdot 0,52}}\right]}}=0,023{\text{ рад}}=1,29^{\circ }.                 (16)

2. Вычислим ширину S области на облаке, образуемой рассеиваемым пучком протонов (рис. 2):

\tan \left({\frac {\theta }{2}}\right)={\frac {S}{2L}}\Rightarrow S=2\cdot L\cdot \tan \left({\frac {\theta }{2}}\right)=2\cdot 588\cdot 0,0113=13,3{\text{ м}}.                                 (17)

Ширина полученной области рассеяния S на облаке значительно меньше его предполагаемого диаметра D, полученного в (14), что удовлетворяет условию попадания всех частиц в облако.

Рис. 2. Конус, образованный рассеянным пучком протонов

В результате проведенных оценок оптимальным диаметром сферического облака представляется значение D = 150 м.

Оценка времени зарядки облака

Для полученного размера облака рассчитаем время, необходимое на зарядку до значения электрического потенциала, при котором возникает пробой воздуха. Напряженность такого электрического поля должна быть E\geq 30{\text{ кВ/см}} [1]. Тогда потенциал сферического облака с радиусом r = 75 м должен быть равен

{\vec {E}}=-{\vec {\bigtriangledown }}\varphi \Rightarrow E={\frac {\varphi }{r}}\Rightarrow \varphi =Er=3\cdot 10^{4}\cdot 75\cdot 10^{2}=225{\text{ МэВ}}.                                  (18)

Вычислим, сколько заряда нужно сообщить облаку, чтобы произошел пробой:

\varphi ={\frac {q}{4\pi \varepsilon _{0}r}}\Rightarrow q=\varphi 4\pi \varepsilon _{0}r=225\cdot 10^{6}\cdot 10^{-10}\cdot 75=1,7{\text{ Кл}}.                                             (19)

Ток пучка в рассматриваемом типе ускорителей достигает величины I = 0,1 А. Тогда искомое время зарядки облака можно получить из определения силы тока:

i={\frac {q}{\tau }}\Rightarrow \tau ={\frac {q}{i}}={\frac {1,7}{0,1}}=17{\text{ с}}.                                                                         (20)

Заключение

В результате проведенных расчетов установлено, что при энергии протонов E = 0,4 ГэВ радиус сферического облака на высоте h = 500 м должен составлять r = 75 м. Для образования электрического разряда в атмосфере этому облаку необходимо сообщить заряд q = 1,7 Кл, на что при токе пучка I = 0,1 А потребуется 17 с.

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

Литература
  1. Базелян Э.М., Райзер Ю.П. Физика молнии и молниезащиты. Москва, ФИЗМАТЛИТ, 2001, 320 с.
  2. Кикоин И.К. (ред.). Таблицы физических величин: справочник. Москва, Атомиздат, 1976, 1008 c.
  3. Коныгин С.Б., Крючков Д.А. Макроскопические свойства дисперсных систем. Самара, Самар. гос. техн. ун-т, 2007, 22 с.
  4. Роджерс Р.Р. Краткий курс физики облаков. Ленинград, Гидрометеоиздат, 1979, 231 с.
  5. Шишкин Н.С. Облака, осадки и грозовое электричество. Ленинград, Гидрометеоиздат, 1964. 404 с.
  6. Ишханов Б.С., Капитонов И.М., Юдин Н.П. Частицы и атомные ядра. Москва, ЛЕНАНД, 2019, 672 c.
Ваш браузер устарел и не обеспечивает полноценную и безопасную работу с сайтом.
Установите актуальную версию вашего браузера или одну из современных альтернатив.