f(x,vx) = fos(x,vx) + fes(x,vx), (1.1.6)
где fos и fes - соответственно функции распределения захваченных (осциллирующих) и уходящих частиц. При сделанном допущении интеграл столкновений зависит только от вида функции fos и уравнение (1.1.4) распадается на два уравнения для функций fos и fes
(1.1.7)
(1.1.8)
где h(y) = 1, при y > 0 и h(y) = 0, при y < 0.
Уравнение (1.1.7) можно подвергнуть дальнейшему упрощению при использовании следующего подхода. Из выполнения условия (3.5) следует, что характерная частота столкновений n много меньше частоты колебаний в потенциальной яме nos, поскольку n ~ <v> / l, а nos ~ <v> / R, т.е. большую часть времени частицы совершают свободные колебания. Столкновения, при которых меняется энергия частицы и, соответственно, амплитуда колебаний a (переход к амплитуде e = ej(a) возможен при соблюдении симметрии, j(x) = j(-x), что обычно выполняется), редки, и поэтому нахождение распределения fos(x,vx) можно разбить на два этапа, сначала рассмотреть сравнительно медленный процесс релаксации частиц по энергии или амплитуде, определить соответствующую функцию распределения f(a) , а затем выразить fos(x,vx) через f(a), считая, что частицы совершают свободные колебания. Для функции, зависящей от амплитуды, левая часть (1.1.7) oбращается в 0 и (1.1.7) принимает вид
S[a,j(a)] h(ej(a))=0. (1.1.9)
В работе [83] проведена подробная процедура преобразования кинетического уравнения в фазовом пространстве к уравнению в пространстве интегралов движения. В настоящей работе связь между fos(x,vx) и f(a) и интегральное уравнение для f(a) получим непосредственно из физических соображений.
Частица, находясь в точке х, будет иметь скорость в промежутке vx ¸ vx +dvx в том случае, если ее амплитуда находится в промежутке a¸ a + da, где a, da связаны с vx, dvx следующими соотношениями
, (1.1.10)
. (1.1.11)
Относительное время нахождения такой частицы в промежутке dx
, при ½x½< a, (1.1.12)
где dt(a,x)=dx/V(a,x), а T(a)=4 - период колебаний. Таким образом, для числа частиц в элементе dxdvx фазового пространства получаем соотношение
. (1.1.13)
Подставляя в (1.1.13) соотношения (1.1.10), (1.1.11), получим
, (1.1.14)
где a(x,vx) = j-1(j(x) + Mvx2/2e), j-1 - функция, обратная к j(x). Интегрируя fos(x,vx) по всем возможным скоростям, получим распределение по поперечной координате
, (1.1.15)
где Vm = V(R,x) = V(a,x)a=R .
Уменьшение числа частиц в интервале a ¸ a + da происходит за счет столкновений и равно n(a)f(a)da, где n(a) - усредненная по периоду колебаний частота столкновений. Новые частицы появляются как за счет столкновений, в результате которых частица, имевшая амплитуду a', попадает в промежуток a ¸ a + da, так и за счет ионизации. В общем виде интегральное уравнение можно записать следующим образом
-n(a)f(a) + G(a)+òW(a,a¢)f(a¢)n(a¢)da¢ = 0, (1.1.16)
где G(a) - распределение частиц, образующихся в единицу времени в результате ионизации, W(a,a') - плотность вероятности для частицы, имевшей амплитуду a', приобрести в результате столкновения амплитуду a. Отметим, что , т.к. есть вероятность получeния частицей в результате столкновения энергии, достаточной для преодоления потенциального барьера. Образование новых частиц происходит, как правило, при ионизации электронным ударом и определяется параметрами распределения быстрых электронов. Принимая, что в момент появления частицы имеют максвелловское распределение по скоростям, для Gi(a) получаем следующее выражение
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.