Файл: Вычислительные методы в физике плазмы..pdf

следуя Кауфману [11], мы определим критический угол в про­ странстве скоростей для потерь ионов и электронов. Обозначим величину магнитного поля в средней плоскости через В 0, а у про­ бок — через Вмакс. Пусть потенциал плазмы в средней плоскости равен ер и уменьшается до нуля у пробок. Пусть W — кинетиче­ ская энергия частицы; тогда полная энергия Н = W + е ф — постоянная величина в случае движения заряженной частицы. Мы предположим, что магнитный момент К = W±]В также является константой движения, где W ^ = Ѵ2 тѵ\ и гр — компонента ско­ рости, перпендикулярная магнитному полю. Критический угол а в средней плоскости определяется из соотношения

О,

Это выражение использовалось также в работе [7] при исследова­ нии стационарных состояний.

Использовалась следующая процедура определения vk. На каж­ дом шаге по времени из уравнений (20) и (24) рассчитывались величины пе (т) и тіі (т) и их разность пі (т) — пе (т). В про­ цессе образования плазмы электроны теряются из установки быстрее, чем ионы. Поскольку мы требуем выполнения равенства пе (т) = nt (т), то разность щ (т) — пе (т) сравнивалась с задан­ ным малым числом. Если она превышала это число, то vk увели­ чивалась на величину Avk, чтобы уменьшить интенсивность потерь электронов, и все вычисления для данного шага по времени повто­ рялись. Этот процесс повторялся до тех пор, пока разность

434 Гл. 10. Решение уравнения Фоккера — Планка

Заменим уравнения (16) и (17) следующими неявными разност­

где V 2 ^ P ^ 1. Нужно решить эти уравнения относительно не­ известных f f +1 и gl+1, j = 0 , 1 , 2 ,

Перепишем эти разностные уравнения как систему алгебраи­ ческих уравнений

для у = 1, 2, . . J — 1. Коэффициенты определяются форму­ лами

— рДт£»"+ 1 - (1 - р) Дт£>”,

- рДтВ"+1 - (1 - р) ДтL".

В уравнениях (39) и (40) неизвестные /"+1, g"+1, / = 1, . . J — 1

стоят в левой части уравнений, а известные величины — в пра­ вой части. Нам нужно решить эти уравнения для внутренних точек / = 1, - - -, J — 1, так как граничные условия при х = 0 и X = X] определяют решения для / = 0 и j = J. Следовательно, нас не должны беспокоить сингулярности в коэффициентах при

436Гл. 10. Решение уравнения Фоккера — Планка

(24)и (25). Если уравнения (16) и (17) решаются для произволь­ ного начального распределения без членов, соответствующих ис­ точникам и потерям, то все упомянутые выше интегралы должны быть константами, что обеспечивает удобный способ проверки программы. Мы нашли, что наилучшие результаты получаются, когда параметр р, входящий в разностные уравнения, полагается равным 1 .

§ 4. Д в у м е р н ы е ( н е гізо т р о п н ы е ) з а д а ч и

1.Граничные условия и независимые переменные

впространстве скоростей

Рассмотрим теперь задачи, в которых не используется предпо­ ложение об изотропности или квазиизотропности распределения ионов. Мы наложим только два упоминавшихся в § 2 ограниче­ ния на это распределение: распределение должно обладать ази­ мутальной симметрией и зеркальной симметрией относительно средней плоскости. В этом случае наиболее подходящи две пары независимых переменных: и и Ѳ или ѵ и ц, где и = cos Ѳ. В обла­ стях пространства скоростей, для которых предполагается суще­ ствование угла потерь, различие между этими парами незначи­ тельно. В этом случае можно потребовать, чтобы функция распре­ деления обращалась в нуль при Ѳ= Ѳп или р. = cos Ѳп [5]. В са­ мом деле, если пробочное отношение мало, то наиболее существен­ ны значения Ѳ около я/2 и cos Ѳ да я/2 — Ѳ. Тогда одну систему координат от другой отличает почти линейное преобразо­ вание.

Однако задачи, цель которых проследить за релаксацией рас­ пределения в отсутствие конуса потерь или других ограничений, должны решаться во всем пространстве скоростей. Более того, при наличии амбиполяриого потенциала распределение электро­ нов простирается до углов Ѳ = 0 и Ѳ = я при низких скоростях. В этих случаях нет четкого граничного условия для функции

Первое условие очевидно, поскольку случаю ѵ = 0 соответ­ ствует одна и та же точка (начало координат) для всех Ѳ. Послед­ ние два условия — результат требования азимутальной симме­ трии распределения. Используя преобразование от декартовой