В.С.Сухомлинов, А.С.Мустафаев DOI 10.18454/PMI.2016.4.611
Пространственное распределение энерговыделения...
Геонаноматериалы
УДК 02.2, 04.1
ПРОСТРАНСТВЕННОЕ РАСПРЕДЕЛЕНИЕ ЭНЕРГОВЫДЕЛЕНИЯ ПРИ РАСПРОСТРАНЕНИИ ПУЧКА БЫСТРЫХ ЭЛЕКТРОНОВ В ВОЗДУХЕ
В. С. СУХОМЛИНОВ1, А.С.МУСТАФАЕВ2
1 Санкт-Петербургский государственный университет, Россия
2 Санкт-Петербургский горный университет, Россия
Работа посвящена разработке аналитической теории для оценки пространственного распределения энерговыделения при распространении пучка быстрых электронов в газе и, в частности, в воздухе при энергиях электронов 1-100 кэВ. Развивается подход, который был применен авторами [2, 3] при рассмотрении неупругого торможения электронов в воздухе. Основываясь на том, что неупругое взаимодействие, в основном, приводит к релаксации энергии, а упругое - к изотропизации распределения по направлениям, в работе вначале решается задача о нахождении функции распределения электронов с учетом только упругих столкновений. В заключительной части находится аналитическое решение поставленной задачи с учетом обоих видов торможения электронов в воздухе.
Проведенные расчеты показывают, что учет упругих столкновений приводит к росту пространственной плотности энерговыделения и сужению области, где выделяется основная энергия быстрых электронов, по сравнению с расчетами, в которых учитывается только неупругое торможение.
Ключевые слова: торможение пучка электронов с энергией 1-100 кэВ в газе, упругое и неупругое взаимодействие электронов, релаксация энергии и импульса, функция распределения электронов по скоростям, кинетическое уравнение Больцмана, численное моделирование методом Монте-Карло.
Как цитировать эту статью: Сухомлинов В.С. Пространственное распределение энерговыделения при распространении пучка быстрых электронов в воздухе / В.С.Сухомлинов, А.С.Мустафаев // Записки Горного института. 2016. Т.220. С.611-621. DOI 10.18454/PMI.2016.4.611
Введение. При создании гиперзвуковых летательных аппаратов актуальна задача релаксации быстрых электронов с энергиями до 100 кэВ в газе. Для оптимизации внешнего аэродинамического обтекания, а также улучшения характеристик двигателей таких летательных аппаратов предлагается использовать МГД-управление газовыми потоками. При этом для создания МГД-эффекта можно применить высокоэнергетичный электронный пучок. Расчет степени ионизации возникающей плазмы сводится к вычислению пространственного распределения выделяемой в газе энергии. В свою очередь, это распределение рассчитывается с использованием метода Монте-Карло. Поскольку начальные энергии электронов велики, а дифференциальное сечение столкновения электронов с атомами и молекулами газа быстро падает с увеличением переданной атому энергии, потери энергии электронов при каждом столкновении малы. Это обстоятельство приводит к необходимости учета многих тысяч столкновений, что даже при современном быстродействии вычислительной техники требует больших затрат машинного времени.
Настоящая работа посвящена разработке аналитической теории для оценки пространственного распределения энерговыделения пучка быстрых электронов с энергией 1-100 кэВ при распространении в воздухе. Развивается подход, который был применен авторами [2, 3] при рассмотрении неупругого торможения электронов в воздухе.
Вначале решается задача о нахождении функции распределения электронов с учетом только упругих столкновений. В заключительной части находится аналитическое решение поставленной задачи с учетом обоих видов торможения электронов.
1. Решение уравнения Больцмана с учетом только упругих столкновений. Воспользуемся физической моделью, развитой в работе [2]. Кроме того, сохраним все обозначения, использованные в этой работе. Напомним, что мы рассматриваем приближение, в котором при неупругих столкновениях электронов с молекулами газа происходит изменение только энергии, а при упругих - направления движения. Как показано в [2], в рассматриваемых условиях это оправдано. Дифференциальное сечение упругого рассеяния электрона на молекулах, приведенное в [8], описывается формулой
da^ = r2 Z 2(1 -р2) dQ р4(1 - cos 0 + 2^)2 '
где r0 - классический радиус электрона; а = Z /137;р = v /c; Z, v и c - заряд ядра атома - мишени, скорость электрона и скорость света соответственно; © - полярный угол рассеяния. Параметр экранирования в приближении Мольер определяется соотношением [8]
ё В.С.Сухомлинов, АС.Мустафаев
Пространственное распределение энерговыделения.
Л = 1,7-10 5ZtilzP!!
Р2
( а2 ^
1,33 + 3,76—
Р
(1а)
Уравнение Больцмана с учетом только упругих столкновений имеет вид
Я/' 2 п 1
ц^ + Л в (Е)/ = 2^ 11 /(Е, ц', г)а в (Е, ц0^ф, (2)
дг 0 -1
где Ле (Е) - обратная длина свободного пробега относительно упругих столкновений; ц0 = cos0 =
= цц ' + tJ 1 — ц'2 д/1 — ц2 cosф - косинус угла рассеяния.
Поскольку функция распределения не зависит от азимутального угла ф, проводя интегрирование по этой переменной в правой части (2), получим
Ц§ + f = |f(E, ц', Z)ge (E, Ц ', (3)
~ , , , ч 2А(л)(1 + 2л —цц')
где ~ = zA е (E); g (^ц'^ц) = --V ' ^
{(ц — ц ' )2 + 4г|(1 + л — цц ' )}1,5
А(л) = л(1 + л); A e (E) = яг02 Z
При этом выполняется условие нормировки
1 1
J ge (E, ц ' ^ цМц ' =j ge (E, ц ' ^ ц)^ц = 1 .
—1 —1
Рассмотрим модельную индикатрису рассеяния gme (E, ц ' ^ ц):
gme (E, ц ' ^ц) = a (E)5(ц' — ц) + a2 (E). (4)
Функции a1( E), a2( E) выберем из условия равенства двух первых моментов модельной и точной индикатрис. Таким образом, получим
1 — a1 ц
(5) ц
где ц = J ge (E, ц ' ^ ц)ц'Уц ', а коэффициенты a1 (E), a2 (E) оказываются не зависящими от ц.
e
—1
Подставляя индикатрису gme (Е,ц ' ^ ц) в уравнение (3) и используя первое из равенств (5), получим
д& 1 1 ~
+1 = -{¡(г', ц'¥ц' при г' = 2а2~. (6)
дг 2 -1
Таким образом, мы свели задачу к кинетическому уравнению Больцмана со сферически симметричной индикатрисой рассеяния. Уравнение (6) следует дополнить граничным условием
/ (г', ц) г=0 =5(ц-1). (6а)
ц>0
Для решения задачи (6), (6а) можно воспользоваться результатами, приведенными в работах [3, 4]. Согласно этим результатам, решение уравнения (3), ограниченное на бесконечности, можно представить в следующем виде:
/(г', ц) = Л + } А(и)ехр^(ц)dм . (7)
Используя граничное условие (6а), можно получить
a2 =
2
ё В.С.Сухомлинов, АС.Мустафаев
Пространственное распределение энерговыделения.
Л0 = 5,036; Л(и) = 1,679-
[ и + р(и)5(и —1)1
и — ц \
C (и) Р(и)
где функции, входящие в (7), определяются формулами [4,5]:
1
Ло = 3|цР(ц) f (0, ц)ф;
о
1
Л(и) = [С(и)Р(и)]— /цР(ц) f (0, ^ (ц)ф;
(8)
(9)
С (и) = и {р 2(и) + л2 и2} при р(и) = 2 + 1п| 1—и |; 2 V 1 + и )
( , -
Р(ц) = ц ехр
— | агС£
( лt ^
V Р(0 )
Л
t + ц
а агС§(Л) изменяется от —л до 0, когда t изменяется от 0 до 1. Функция Gu (ц) задается формулой
Gu (ц) = Рк
и — ц
+ р(и)5(и — ц)
Выражение Рк функционала:
и — ц
понимается как обобщенная функция, определяемая с помощью интегрального
I Рк
и
Я (и)йи = Нт
и — ц 6^0
7 иж ёи +1 иж и
—1 и — ц ц+е и — ц
при произвольной функции Р(и).
Опуская громоздкие вычисления, приведем конечные формулы для величин N (z'), ц+ (z'),
ц— (z'), ц+ 2( z'):
N(z') = 2Л0 + ЯДz') + С1(z'); ц+ (z') =
^ , _ 0,5Л0 + Л1(z').
_ (— 0,5Лр + Д(z') , — 2 ( ,
ц—(z) = —. . ^ , ,ч-; ц (z) =
Л
Л0 + ЯД z') 2Л
Л0 + С1( z')
ц+ 2( z') = 3
3[2 Л0 + ЯД z') + СД z')]' 0 + Gl( z')
(10)
Л0 + ЯД z')
При этом
ЛД 2') = | у )и( у )ехр
( >
V и(У))
1>; Я1(2') = | у)ехр
( >
С1( 2') = | ЯД у )ехр
( ^
V и(У))
ад
Gl( 2') = I Р4( у )ехр
^ А(2') = | ЯДу)ехр
V и(У))
( ^
( 2' ^
V и(У))
V и(У)) ;; и(у) = 1 — ехр(—у),
(11)
и
и
ад
ад
ад
ад
ё В.С.Сухомлинов, АС.Мустафаев
Пространственное распределение энерговыделения.
а функции F0( г'), F1( г'), F2( г '), F3( г') рассчитываются по формулам:
р (г , = -1 679 к (г ') -1 + и(г ')1п(и(г ') ехр(г '))]и(г '), 0( ) , Р(и(г'))С(и(г')) ;
^(г ') = -1 679 [Р:(г') + и(г')1п(и(г')ехр(г'))]и(г'),
Р(и( г '))С (и( г')) г ') = 1,679 и 2( г' )1п(и( г' )ехР( г'));
Р(и(г'))С(и(г')) (12)
^(г ^ = 1 679 [1 + и(г')1п(и(г')ехр(г'))]и2(г'),
Р(и( г '))С(и( г'))
[и 2( г ')1п(и( г ')ех
F4( г ') = -1,679
1 679 [и2(г ') 1п(и(г ') ехр(г '))- и(г ') - 0,5 + и(г ')рх (г ')]и2(г ');
Р(и( г '))С (и( г')) р1 (г' ) = 2 - и( г' )[г' + 1п(1 + и( г'))]. Кроме того, выполняется соотношение
А( г') - АД г') - Ао. (13)
Формулы (7), (8), (9)-(13) и дают решение поставленной задачи. Для конкретных вычислений моментов функции /(г ', ц) оказывается полезной следующая приближенная формула, имеющая погрешность несколько процентов:
Р(и) « Р1(и) = 1,676 -1,095(1 - и)1,1. (14)
Очевидно, что переход от переменной г ' (расстояние, выраженное в длинах пробега, рассчитанных по изотропной части сечения упругого рассеяния), к переменой ~ (расстояние, выраженное в длинах пробега, рассчитанных по полному сечению упругого рассеяния) осуществляется по формуле
г ' = ~ • 2а2( Е) = . (15)
I ц )
Кроме того, для численных расчетов оказывается полезной аппроксимация, точность которой составляет несколько процентов:
п-1
°2(Е) « а2(Е) =
Е2
1 + 0,057-
1 + 0,057Е
(16)
2. Полное решение задачи о распространении высокоэнергетического пучка электронов в воздухе. Как показывают расчеты [2], в интересующем нас диапазоне энергий электронов сечение неупругого рассеяния значительно превосходит сечение упругого рассеяния. Так, при энергии 2 кэВ, их отношение равно приблизительно 20 и растет с увеличением энергии, поскольку, как отмечалось ранее, сечение неупругого рассеяния, в отличие от упругого, слабо зависит от энергии. Таким образом, между сравнительно редкими упругими столкновениями, в результате которых происходит изменение направления движения электронов, электроны движутся по прямой, испытывая неупругие столкновения, при которых происходит потеря энергии.
Рассмотрим сначала, как влияет на полученные результаты энергетическая релаксация между последовательными упругими столкновениями. Для оценки этого эффекта будем считать, что электроны испытывают только упругие столкновения, но перед каждым следующим столкновением их энергия уменьшается и равняется средней энергии потока Е (г), рассчитанной с учетом только неупругих столкновений в работе [2]. При этом все расстояния будем выражать в длинах пробега относительно неупругих столкновений, т.е. перейдем к переменной г0 (Е0 )> г . Очевидно, что теперь во всех формулах следует в качестве пространственной переменной писать величину X(г):
X(г,Е0) = |«2[Е(гЖ . (17)
«(Е0) 0
ё В.С.Сухомлинов, АС.Мустафаев
Пространственное распределение энерговыделения.
Вычисляя интеграл в этой формуле, получим
X (X, Ео) = - 1
Т'( г, Е о) + 4,098Т (X, Е 0) К(Ео) - ' ' ' 0'
Т1( X, Ео) = arctg[0,12(2 Ео +1)] - аге^[0,12(2 EоX( X, Ео))];
Т (X, Е о) = 1п
1 + о,о57 Ео(1 + Ео)
1 + о,о57 Е 0t (X, Ео )(1 + Е^ (X, Ео))
(17а)
- Ко(Ео)
t( х, Ео) = 1 + х
о Е о
Рассчитаем, при каких значениях X происходит полная изотропизация углового распределения потока электронов, т.е. выполняется соотношение
ц+[Х(X)]«Ц-[Х(X)]« 1.
С учетом полученных выше результатов имеем, что соответствующая величина Е1 (X, Ео) находится как решение трансцендентного уравнения
X к(Ео)]= 4. (18)
Величина 20 (Ео) оказывается много больше единицы и, например, для Ео = 2 кэВ выполняется го = 9о, а для Ео = 100кэВ параметр 20 = 3,8 • 1о4 .
В работе [2] было показано, что зависимость Е(X) хорошо описывается линейной функцией. При этом практически полной релаксации по энергиям соответствует расстояние от плоскости вылета:
^о^КЕт <19)
где функция К (Ео) определена в [2].
Расчеты величины г(Е0) = z0 (Ео) / Ze (Ео) показывают, что при начальных энергиях электронов 2 кэВ < Ео < 100 кэВ ее значение находится в диапазоне от 18 до 35 %. При этом оценки средней
энергии электронов в области X < z0 (Ео) с учетом неупругих столкновений показывают, что ее отличие от начальной энергии не превосходит 0,2 %. Это свидетельствует о том, что изотропизация углового распределения пучка электронов происходит в указанном диапазоне энергий практически без релаксации энергии электронов и на расстояниях от плоскости вылета, заметно меньших размеров области существенного энерговыделения. Но тогда можно предложить следующее решение задачи с учетом неупругих и упругих столкновений электронов с молекулами газа.
Рассмотрим сначала область X > z0 (Ео) . Как мы видели, при X = z0 (Ео) распределение электронов - практически изотропное. Далее электроны испытывают неупругие столкновения с потерей энергии и упругие с изменением направления. Но тогда легко доказать, что при любом X > z0 (Ео) функция распределения электронов изотропна. В этом случае, очевидно, упругие столкновения с изменением только направления не приводят к изменению функции распределения и их можно не учитывать. Кроме того, в рамках модели, предложенной в [2], будем считать, что полное сечение неупругих столкновений слабо зависит от энергии электрона. Кинетическое уравнение для функции распределения электронов в этом случае записывается в виде
Ц^(У, + /(X, X, ц) = а} /(X, X', ц)спй (X - X')йХ', (20)
C'2 о
где ц - косинус угла между внутренней нормалью к плоскости и направлением движения электрона;
х
X = 1 -х; х = Е/Ео; а = 1/(1,Ео); (х,Ео) = Ео{с„й(х,х-х')<яХ'.
о
Как было показано в работе [1], в данной ситуации для решения кинетического уравнения можно применить так называемое ^-приближение. Полагая, как и в работе [2], что сечение с па (X - X') имеет экспоненциальную форму:
-й- В.С.Сухомлинов, АС.Мустафаев 00110.18454/РМ1.2016.4.611
; Пространственное распределение энерговыделения...
а пй ((- О = 2 „ (1, Е0 )Ь( Е0) ехр[- Ь( - О], (20а)
(формулы для величин, входящих в (20а), приведены в [2]), представим / (г, t, ц) в виде
/ (г, ^ ц) = 2 [/(0) (г, t) + 3ц/(1) (г, t)]. (21)
Подставим (21) в (20) и проинтегрируем по всем направлениям движения электронов, затем подставим (21) в (20), умножим на ц и также проинтегрируем. В результате получим следующую систему линейных интегродифференциальных уравнений для функций /(0) (г, t), /(1)(г, t):
+ /(0) (г, t) = а | /(0) (г, t ' ^ (t -1 ' )Л';
дг
0
t
1 ^^^О + / (1) (_, t) = « ^ / (1) (_, t ' -1 ' )Л' . 3 дг 0
(22)
Граничное условие для системы (22) записывается в виде
/(0)(г0,0 = Фо(t); /(1)(г0,t) = ФДО , (23)
где Ф),Ф^) - некоторые, пока неизвестные функции.
Исключая функцию / (0)( г, t) из системы (22) и применяя к полученному уравнению преобразование Лапласа по переменной t, получим уравнение для образа Лапласа рДг, р) функции /(1)(г, 0 :
д 2 р
-41 -3[1 -а/(р)]2Р1 = 0; рх(г0,р) = Ох(р), (24)
дг
где G1(р) и I(р) - преобразование Лапласа соответственно от функции Ф^) и апд^).
Полагая, что сечение а пд ^) определяется формулой (20а) и применяя обратное преобразование Лапласа к решению уравнения (24), с учетом (23) имеем
/(1) (г, t)г>г0 = Ф1(t) ехр{-л/3(г - г0)} + ехр{-л/3(г - г0»л/а(г - г0) х
х{ Ф1(t - у (2л/а<г-г0хУ) ехр(-ау) ду. (25)
0 л/У
Рассмотрим теперь область г0 (Е0) > г . Выполняя интегрирование уравнения Больцмана для функции распределения электронов с учетом неупругого и упругого торможений [2], для области г0 (Е0 )> г получим систему уравнений:
/д(М+/(0) (г, t) = Ь( Е0) | /(0) (г, t ' Кд (t -1 ' д'; дг 0
д[ц2(г/0)(г,0] + /(1)(_,t) = ь(Е0)|/(г, t ')апда -1 ')Л',
дг 0
1 1(1 \ . 2
(26)
где ц 2( г) = |ц2 / (г,0, ц)дц / | / (г,0, ц)дц
Действуя аналогично предыдущему, получим
д [ця_(2г)р0]-[1 -Ь(Е0)/(р)]2Р0 = 0; р0(г0,р) = Gо(p), (27)
дг
где G0( р) - преобразование Лапласа от функции Ф ).
Как показывают оценки, уравнение (27) удовлетворяет условиям применимости ВКБ-приближе-ния [6]. В этом приближении решение уравнения (27) имеет вид
ё В.С.Сухомлинов, АС.Мустафаев
Пространственное распределение энерговыделения.
gо(г,р) = ^(р)^Ё21ехр{-г}ехр\ь(Е0)1(р)\^== [ . (28)
л/ц 2( г) I Ч ц 2( У) ]
Из (26) аналогично случаю г0 (Е0) < г получим
2 И
gl (г, р) = Gl (р) ехр{- (г - 2о )}ехр<| Ь(Ео )1 (р) {
Г
Применяя обратное преобразование Лапласа, для области г0 (Е0) > 2 окончательно имеем /(1)(2,0,-<20 = ФДОехр\- М=! + ехр{- } Иг
--- г I 1 j ,---
Ц 2(г) ] [ 20 4 ц 2(г) 4Ь(Ей)(? - 20)}Ф1 (X - У)11 (27Ь(Е0)(г - ¿0 )у )ехр(- ^У) Иу. (29)
л/У
ИЕ (г, Е0)
Напомним, что цель настоящей работы - вычисление величин Е(г, Е0); ---. При этом
Иг
кроме интегрирования функции распределения по переменной Е (с весом Е) следует провести интегрирование по всем направлениям движения с весом, равным ц . Очевидно, что когда функция распределения определяется выражением (21), слагаемое, в которое входит /(0), обратится в ноль, поэтому нет необходимости в вычислении симметричной части функции распределения электронов.
Теперь следует найти функцию, которая в соответствии с определением имеет вид
0 С 20 Иг 1
ФДО = {ц-/(0,X,ц)Иц + 5(0ехр{-. (29а)
Покажем, что в диапазоне начальных энергий электронов Е0 > 2 кэВ первым слагаемым в выражении для Ф 1(Х) (29а) можно пренебречь. Из изложенного выше ясно, что поток электронов, движущихся из полупространства к плоскости вылета, слабо анизотропен. Тогда в кинетическом уравнении при ц < 0 можно пренебречь упругими столкновениями. Электроны, движущиеся к плоскости вылета, появляются из-за изменения направления движения электронов в результате упругих столкновений. Ясно, что чем меньше величина ц > 0 перед упругим столкновением, тем больше у этого электрона вероятность иметь после столкновения направление движения к плоскости вылета, т.е. иметь после столкновения ц < 0 . Поскольку мы хотим провести оценку сверху, будем считать, что все электроны с ц > 0 движутся с минимально возможным средним косинусом между нормалью к плоскости вылета и направлением скорости, т.е. с ц = 0,5, соответствующим изотропному распределению. Интегрируя кинетическое уравнение по ц в пределах от 0 до -1 с весом 1 и ц, получим систему уравнений:
(1) (^, 0 + /-(0) (2, X) = Ь( Е0)} /-(0) (2, X ')с „и (X - X ' + «1 (2, X);
л— ^ — \ у , \ К) ' Л — V ' Пи ^ ' 1 / г-1 \
& 0 „(Е0)
СТе (Е0 )
0 „ ( Е0
1 + /- (!) (г, X) = Ь( Е0) | /-(1) (2, X ')с „И (X - X уЛ' + «0 ^^ /+ (2, X),
3 & 0 „(Е0)
(30)
где а1 = |ge( Е0,1 —> ц,1 |цИц; /+ (2,X) рассчитывается по формуле (23а) работы [2];
-1
«0 = | ge ^Е0^ — ц |Иц .
X
ё В.С.Сухомлинов, АС.Мустафаев
Пространственное распределение энерговыделения.
Применяя к системе (30) преобразование Лапласа и исключая функцию р -(0), получим дифференциальное уравнение второго порядка для функции р - (1). Решение этого уравнения имеет вид
Р-(1) (а0 - 3«1 )ехр[Р(р)г], (31)
4Р(р) ап(Е0)
где р-(1) - преобразование Лапласа функции /-(1)(г,t); Р(р) = 1 -Ь(Е0)1(р).
Когда 1(р) определена формулой 1(р) = Ь(Е0)/[р + Ь(Е0)], применение обратного преобразования Лапласа дает
/-(1) (0, t) =1 ^^ (а0 - 3«1 )[§(0 + Ь(Е0)]. (32)
4 а п(Е0)
Соответственно, для потока энергии из полупространства г > 0 через плоскость вылета имеем
Ь( Е0)
Е (0) = Е-^Е) (а0 - 3а!)
4 а п (Е0)
1 +-
(33)
2
Расчеты показывают, что величина Е- (0) / Е0 при Е0 = 2 кэВ составляет единицы процентов от Е0 и быстро падает с увеличением начальной энергии электронов. Физически это объясняется тем, что, во-первых, угол рассеяния при упругих столкновениях сравнительно мал, во-вторых, сечение упругого рассеяния мало по сравнению с сечением неупругого рассеяния. Таким образом, при расчете энерговыделения в газе при г = 0 можно не учитывать функцию распределения электронов, движущихся к плоскости вылета. Тогда с указанной точностью выполняется соотношение
Ф^) = 5(0ехргММ . (34)
I Ч ц 2(г) ]
Используя формулы (29) и (34), для средней энергии и энерговыделения в области значений г0 (Е0) > г можно получить следующие соотношения:
Епе (г) = Еп (Y (г, Е0)); (35)
дЕпе (г) = 1 дЕп {у (г, Е0)).
7ц2 (X (г, Е0)) ду '
у (г, Е0) = 1 д
0
2 [X(г, Е0)] '
где величина X(г,Е0) рассчитывается по формулам (17а), (18) настоящей работы, а величины
- _ дЕ (г)
Еп (г) и —п--по формулам соответственно (24) и (25) работы [2].
дг
При численных расчетах для сокращения времени счета удобно пользоваться следующей аппроксимацией для функции ц 2(г):
Г Г И2
I П 1П I
ц 2( г) = <¡0,578
0,313
1 +
(1 + 0,55г )4,2
Отметим, что поскольку при г0 (Е0) < г с точностью до десятых долей процента выполняется соот-
00
ношение
л/ц2 (X (г, Е0))
то формулы (35) справедливы во всем полупространстве. Следует иметь в виду, что формулы (35) настоящей работы дают возможность учитывать упругие столкновения вне зависимости от метода,
ё В.С.Сухомлинов, АС.Мустафаев
Пространственное распределение энерговыделения.
которым рассчитывается релаксация энергии электронов за счет неупругих столкновений (в том числе, и численно).
Основываясь на вышеизложенных результатах, для величины зоны энерговыделения можно получить формулу
Я
hne (E0) = Z0 +"
К
E0
1 -
Zo( Ео)К( Eo)
E0
Y
(36)
3. Обсуждение полученных результатов. На рис.1 приведены расчеты средней энергии пучка электронов с учетом обоих видов торможения Епе (г, Е0) и с учетом только неупругих столкновений Еп (г, Е0) для начальных энергий вылета электронов соответственно 10 и 5 кэВ. Видно, что учет упругих столкновений приводит к более быстрой энергетической релаксации, что, естественно, должно приводить к более высоким пиковым значениям энерговыделения в газе. На рис.2 для начальной энергии 2 кэВ приведены величины dEn (г, Е0) / ^; dEne (г, Е0) / dz , рассчитанные нами и определенные численно методом Монте-Карло (dEne (г, Е0)/ dZ) [9], а также методом разложения функции по числам столкновений (dEn (г, Е0) / dZ) [9]. Из этих данных видно, что аналитические расчеты удовлетворительно совпадают с результатами численных расчетов [9]. Кроме того, зависимости, приведенные на рис.2, свидетельствуют о том, что, во-первых, действительно, величина dEne (г, E0)/dZ превосходит dEn (г, E0) / dZ. Во-вторых, эта величина сначала заметно растет, а затем, пройдя максимум, падает до нуля.
Рассмотрим физические причины указанных закономерностей. При учете только неупругих столкновений энергетическая релаксация пучка электронов происходит при движении по нормали к плоскости вылета г = 0 , т.е. в элементе объема dV = dxdydz, находящемся в точке А1(х1,у1,г1), выделяется только часть энергии электронов, которые вылетели с элемента плоскости dS = dydx из точки А2(х1,у1,0). Учет упругих столкновений, вызывающих изменение направления движения электронов, приводит к тому, что в элемент объема dV попадают электроны, которые вылетают из других точек плоскости г = 0 с координатами х ф х1; у ф у1 (рис.3). Это, естественно, приводит к увеличению энерговыделения в элементе объема dV . При этом ясно, что чем меньше средний косинус угла между нормалью к плоскости вылета и направлением движения электронов (г!), тем (при фиксированной координате г1) больше часть плоскости вылета, «видимая» из точки А1. Таким образом, уменьшение (г1) приводит к более быстрой энергетической релаксации электронов. Это и является причиной того, что сначала dEne (г, E0) / dZ растет. После того, как величина (г1) достигает своего предельного значения, равного 0,5, рост удельного энерговыделения прекращается.
« Ene (Z), En (Z ),кэВ 10 -
5 -
б En(Z),Ene(Z),K3B
1-1—I I I I I 11-1-1—I I I I I 11-1-1—I I I I I 11-1-1—гт
1 10 100 1000 z
10
1-1-1—I—I I I I I
100
~I-1-1—I—I I I I I
1000 z
Ene (Z)
En (Z)
Рис. 1. Зависимость средней энергии электронов от числа столкновений, рассчитанной с учетом упругих и неупругих столкновений £ (г) и только неупругих столкновений ЁП (г) при начальной энергии электронов 10 кэВ (а) и 5 кэВ (б)
0
ё В.С.Сухомлинов, АС.Мустафаев
Пространственное распределение энерговыделения.
ИЕ / кэВ/ст. 0,01 -
1Е-3 -
1Е-4
~1-1-1 I I I I |-
10
100
ИЕ„ /И!, расчет по формулам работы [2]
ИЕ„ /и!, по данным [9]
ИЕ„е / и!, по данным [9]
ИЕ„е /и!, расчет по формулам (35)
Рис.2. Зависимость удельного энерговыделения от числа столкновений (ст.), рассчитанная с учетом упругих и неупругих столкновений ИЕ„е / Иг и только неупругих столкновений ИЕ„ / Иг при начальной энергии электронов 2 кэВ
(хЬ У1, 21)
Рис.3. Пояснение к соотношению величин ИЕ„е (г )/ Иг и ИЕ„ (г )/ Иг
х
1
„е / ¿г )тах,кэВ/ст.
0,015
0,010
0,005 -
0,000 -,—,-1-,-1-,-1-,-1-,-1-,
10 20 30 40 50
Е0, кэВ
Рис.4. Зависимость максимального удельного энерговыделения от начальной энергии электронов
КЕ0), Н(Е0)
10
10
о ГЕ0)
20
........................
30 40 50
Е0, кэВ
Н(Е0)
Рис.5. Сравнение относительной ширины зоны энерговыделения г(Е0) и Н(Е0), соответственно измеренной автором работы [7] и рассчитанной нами
В работе [7] было экспериментально показано, что зависимость размеров зоны энерговыделения R(Е0) и начальный участок зависимости Е1(г, Е0) в области начальных энергий 5 кэВ < Е0 < 54 кэВ хорошо описываются следующими соотношениями:
Я(Е0) = «1Е01,66; ЕДг,Е0) = Е,
(
0/ _ -^0
1 —
г >
R( Е0)
0' У
Из второго соотношения (37) легко получить
ИЕ,(0,Е0)
= Е1( г, Е0) г=0 =«2 Е0
-0,66
(37)
(38)
где «1, «2 - некоторые численные коэффициенты. 620
1
0,6
В.С.Сухомлинов, А.С.Мустафаев
Пространственное распределение энерговыделения...
На рис.4 приведены сравнительные данные о е(Е0), ЩЕ0) максимальном удельном энерговыделении в газе 1 при различных начальных энергиях электронов, а на рис.5, 6 - сравнительные данные расчета величин, не зависящих от констант а1, а2, соответственно:
0,5 -
r(Eo) = R(Eo)/R(5keV); H(E0) = hne(E0)/hne(5keV);
_ f _ f
e( E0) = Ene (Z E0) z=0 / Ene (Z 5keV W
. t
D(Eo) = E (z,E0)г=о/Ei (Z,5keV)?=o.
0,2 -
10
T—I—I II I | I I I I I I I I | I I I 11 I 1111111IIII |llll
20 30 40 50
Е0, кэВ
e(Ec)
D(E„)
Рис.6. Сравнение удельного энерговыделения электронов е(Е0) и D(Е0), соответственно измеренного автором работы [7] и рассчитанного нами
Видно, что аналитические расчеты данной работы хорошо описывают экспериментальные результаты [7]. Кроме того, с ростом E0 величина
dEne (г, E0)/ dZ уменьшается, а размер зоны (рассчитанный для энергии электронов E0 и выраженный в числе неупругих столкновений), где в основном выделяется энергия, - растет. Это связано с тем, что при увеличении энергии вылета падает средняя
энергия, переданная при одном столкновении электрона с атомом газа. При этом данное уменьшение происходит по нелинейному закону. Кроме того, при увеличении энергии вылета электронов растет размер области, в которой происходит изотропизация функции распределения электронов. Указанные закономерности и приводят к наблюдаемой зависимости величины dEne (г, E0) / dZ от E0.
Заключение. Перечислим основные результаты, полученные в работе.
• Предложена физическая модель энергетической релаксации быстрых электронов в газе с учетом упругих и неупругих столкновений. Рассчитанная массовая тормозная способность воздуха хорошо описывает известные экспериментальные данные.
• Обнаружено, что учет упругих столкновений важен при энергетической релаксации пучка электронов. В частности, это приводит к значительному увеличению энерговыделения в газе и уменьшению области, в которой, в основном, происходит выделение энергии.
• Получены аналитические выражения для средней энергии по потоку и значения знерговыделения в газе, что существенно (в десятки раз) сокращает время расчета параметров зоны энерговыделения по сравнению с расчетами методом Монте-Карло.
• Проведены расчеты знерговыделения при различных начальных энергиях вылета электронов. Результаты находятся в удовлетворительном согласии с известными данными экспериментов и численных расчетов по методу разложения функции распределения по числам столкновений.
ЛИТЕРАТУРА
1. Марчук Г.И. Методы расчета ядерных реакторов. М.: Госатомиздат, 1961. 667 с.
2. Сухомлинов В.С. Аналитическая теория релаксации энергии при распространении пучка быстрых электронов в газе / В.С.Сухомлинов, А.С.Мустафаев // Теплофизика высоких температур. 2016, в печати.
3. Сухомлинов В.С. Влияние неупругих столкновений на релаксацию энергии пучка быстрых электронов в воздухе / В.С.Сухомлинов, А.С.Мустафаев // Записки Горного института. 2016. Т.219. С.392-402. DOI 10.18454/PMI.2016.3.392
4. Черчиньяни К. Математические методы в кинетической теории газов. М.: Мир, 1973. 245 с.
5. Черчиньяни К. Теория и приложения уравнения Больцмана. М.: Мир, 1978. 495 с.
6. Froman N. JWKB Approximation / N.Froman, P.O.Froman; Inst. of the Theoretical Phys., University of Uppsala. Sweden. 1965. 168 p.
7. Grun A.E. Lumineszenz-photometrische Messungen der Energieabsorption im Strahlungsfeld von Elektronenquellen Eindimensionaler Fall in Luft // Z. Naturforsch A. 1957. Vol. 12a. P.89-95.
8. Mayol R.R. Total and transport cross section for elastic scatterig of electrons by atoms / R.R.Mayol, F.Salvat // Atomuc Data and Nuclear Data Tables. 1997. Vol.65. N 21. P.55-154.
9. Sheikin E.G. Calculation of Space Disribution of Energy Deposited by E-Beam for Flow Control Applications / E.G.Sheikin, V.S. Sukhomlinov //AIAA Paper. 2006. P.1369-1374.
Авторы: В.С.Сухомлинов, д-р физ.-мат. наук, доцент, prima-ivs@mail.ru (Санкт-Петербургский государственный университет, Россия), А.С.Мустафаев, д-р физ.-мат. наук, профессор, alexmustafaev@yandex.ru (Санкт-Петербургский горный университет, Россия).
Статья принята к публикации 18.01.2016.