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

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 11.04.2024

Просмотров: 259

Скачиваний: 6

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Дополнение

4 89

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

§ 6*. Д и ф ф у з и я ч а с т и ц

Диффузионный поток частиц DS в уравнении (4) можно пред­ ставить в виде

DS = D I (1

 

 

 

 

(18>

где

 

 

 

 

Фг

dj, еЧ~3 V»г, е

di = 1,33,

de=

l , l l .

1+У ^Гв ’

 

 

 

 

Член с квадратными скобками соответствует

термодиффузии.

В «банановой»

области (а;, е <С 1) фг. е ~

dif е і

фг. е — 3/г < О

и термодиффузионный член дает поток частиц по направлению градиента температуры. Коэффициент диффузии D равен

D = ^ l e .

(19)

Дрейфовый поток частиц геГ в том случае, когда электроны находятся в «банановом» режиме, определяется пинчеванием запертых частиц в скрещенных полях Ez и II0,

f =

Г = 1,5 ^ еѴ2 _ І _ J L ln (я2ц),

(20)

а

аг

Tq ^

 

 

(Г = 2, W

Er

 

 

« i ) .

 

В режиме плато запертых частиц нет и Г = 0. Ряд авторов отме­ чает, что пинчевание частиц приводит к дополнительным членам в уравнении (2):

<?е = <?др = 0 , 2 ^ Г е^ (г ш Г ) + 1,2Г дТ

(21)

 

 

дх

 

На периферии плазмы

 

 

 

^др

 

ПрТе (эрг)

(22)

<?Дж

e1/2ß/<l,

ßi=: (1/8я) Щ (Э)

 

Однако в центральной части шнура первый член в формуле- (21) возрастает и 0 ДР может оказаться сравнимым с (?дж.


490 Дополнение

§ 7. Н е й т р а л ы о ст а т о ч н о го га за

Ионизация атомов остаточного газа играет существенную роль в балансе частиц. В балансе энергии роль перезарядки и иони­ зации до сих пор была невелика, однако с ростом температуры ионов и переходом в «банановый» режим уход энергии через ней­ тралы может стать основным источником потерь. Если через N (ж, t) обозначить плотность нейтралов в единицах 1013 см-3, через Тп (ж, t) — температуру нейтралов в электронвольтах, то уход энергии с нейтралами и увеличение плотности за счет иониза­ ции могут быть описаны дополнительными членами в уравнениях

(1) — (4):

_

 

& = - 5 6 N ( T , - T n) ( j / ^ оех+ і , 07 У ¥ еае) ,

Q e = - 60N {Тё- Т п) У Т ео*,

(23)

Р ^ Р п = 60пЫѴтеое,

 

где аехЛЛ0~15 см2 — сечение

перезарядки, ае-10~16

см2 — сече­

ние ионизации атомов электронами. При выводе (23) предполо­ жено, что Тп < Те,і.

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

распределения нейтралов / (ж, ѵ) в

слое — а < х <

а будем

иметь следующую задачу:

 

 

v j L + s f - F N

1

 

 

2у0 [ М ^ ) + 6( ^ ) ]

(24)

 

N о

при И Зс 0.

(25)

f ( ± a , v)

г0'6(Ч^)

 

 

 

Здесь

 

 

 

s = s (ж) = 0,04оехпѵі +

10-Заегане,

 

F = F (х) = 0,0івехпѵі, Ѵі — Ѵі (ж) = 1,4- ІО6 j /" у

,

Ѵе = Ѵе(Ж )=0,6.1081/Ге,

 


Д ополчение

491

и Vq — плотность и тепловая скорость падающих из вакуума нейтралов. Из выражений (24) —(25) нетрудно получить интегралъ-

СО

ное уравнение для плотности нейтралов N (х) = ^ / (х , ѵ) dv:

а

00

N (х) = N 0 (х) + j К (X, I) N Ц) dl,

(26)

о

 

где

 

К{х’

^ =4--£§-{Ф(^

X

sign(*-£))+

 

+ Ф(0, I, Ѵі (і ))Ф(0,

X, Vi(x))},

N 0(x) = № { Ф(ж, а , ѵ„)+ Ф(0,

а,

у0)Ф(0, x , v q)},

 

Ф (l, x,i>) = exp

j s(x')dx'J .

 

 

i

 

Учитывая наличие ионизации, уравнение (26) естественно решать методом последовательных приближений:

оо

а

 

'N (X) = 2 N k (х),N h+1 (х) = (

К (X, I) N h (I) dl.

(27)

fe=0

о

 

Нулевой член этого ряда N 0 (х) дает плотность нейтралов без учета вторичных нейтралов перезарядки, а каждый следующий член описывает нейтралы, появляющиеся при ^-кратной переза­ рядке. При плотности плазмы п ~ 2—5 член Nu+i отличается от N k (при /с > 1) почти на полпорядка. Таким образом 4—5 итера­ ций обеспечивают достаточную точность вычислений. По плот­ ности N (х) можно установить пространственное распределение

температуры

нейтралов

 

 

а

 

N (X)

Тп (X) = N 0 (X) Тп0-\- j К (X, I) N (I) Ть Ц) dl,

(28)

 

О

 

где Тпо = тѵ\І1 — температура падающих нейтралов.

зави­

Плотность

и температура нейтралов параметрически

сят от времени через функции F , уг, s и №. В плазме с малой плот­

ностью (п <

0,5) нейтралы

хорошо

проникают

в центральную

часть N (0) ~

0,5іѴ (а). При

п > 2

основная

часть нейтралов

ионизуется и испытывает перезарядку в периферийном слое: при

п ~ 5, N(0) - О,ОШ (а).

Нетрудно оценить долю нейтралов р, отражаемых плазмой (альбедо плазмы). Если и Тп0 известны, а У (0) и Тп (0) полу-


492 Д ополчение

чены из решения уравнения (26) и формулы (28), то

 

 

 

Ѵ2

ІѴ(0) ГП(0)

(29>

 

р = ф

NO - ‘П

\ № Т п0

где /+ и

— падающий и выходящий* /

потоки атомов

 

§ 8.

М етоды

д о п о л н и т е л ь н о г о н а грева п л а з м ы

 

1. Поджатие магнитным полем

Адиабатическое увеличение основного магнитного поля Н при­ водит к поджатию плазмы и дополнительному нагреву частиц. Если при этом большой радиус тора R не изменяется, то в гидро­ динамическом приближении плотность плазмы п должна возра­ стать с полем линейно, а Тг и Те — как Н 2/з. Л. А. Арцимович высказал идею о возможности нагрева плазмы при уменьшении R и перемещении шнура в область с более сильным магнитным полем (Н ~ R -1). Эта идея недавно была реализована на установке «Токопрессор» в Принстоне. В этом случае из гидродинамики сле­

дует, что п ~ R~2, Т і и Те ~ і?_4/з.

Для расчета баланса энергии и частиц при поджатии систему

(1) —(4) следует дополнить членами

*

Ь ’‘ = { 1 } т ^ ' + Т х * Т Г ’

 

р = ^ 2 } а п + т х ^ '

(30>

^ = { _ ? } а!Х +Т * - | г ’

 

где а = (1/Н) (dH/dt). В фигурных скобках верхние коэффициенты соответствуют случаю R const, нижние— случаю R ~ Н~х.

2.Использование тора с некруглым сечением

Вобычном токамаке с круглым сечением плазменного шнура величина продольного тока ограничивается условием Крускала — Шафранова

q (a ) ~ ' n ( a ) ~~ 0 ,2 В / > 1 -

(31^

В современном эксперименте пользуются еще более жестким условием q > 3. В противном случае приходится принимать спе­ циальные меры по подавлению второй гармоники гидродинамиче­ ских винтовых колебаний. Условие (31) означает, что при одном


Д ополнение

493

обходе тора прокручивание силовой линии вокруг магнитной оси не должно превышать 360°. Для некруглого сечения прокру­ чивание меньше и условие (31) становится более слабым. В част­ ности, для эллиптического сечения с большой полуосью, парал­ лельной главной оси тора, условие (31) принимает вид

q ~ 0,2R I 2 ГДѲ а ~ а

а и Ъ — длины полуосей эллипса. А. А. Галеев показал, что коэф­ фициенты теплопроводности ионов и электронов при этом умень­ шаются пропорционально множителю

кі - « - 2)

Е(ф /і —а -2)

где К и Е — полные эллиптические интегралы.

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

1)увеличение допустимой плотности тока,

2)уменьшение коэффициентов теплопроводности. Количественное описание этих факторов для каждого типа сече­ ния требует детального анализа.

3. Инжекция пучка быстрых нейтралов

Инжекция пучка горячих нейтралов — один из перспектив­ ных методов нагрева плазмы в больших установках до темпера­ туры зажигания термоядерной реакции. Для того чтобы пучок играл заметную роль в общем энергетическом балансе плазмы, его мощность должна быть сравнимой с мощностью джоулева нагрева Wgw « 2n2Ra2Qj^m. Заметим, что в современном эксперименте ГГдж ~ 100—200 кВт. Для нагрева большого токамака-реактора

с

параметрами R = 500, а = 150 до температуры T it е ~ 10 кэВ

с

запасом энергии ІО7 Дж при тв ~ 10 с необходим пучок мощ­

ностью W b = Е0і ~ 103 — ІО4*кВт (в дальнейшем Е 0 — энергия частиц пучка в кэВ, і — эквивалентный ток пучка в амперах).

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

1) Прежде всего надо определить количество быстрых атомов, испытавших ионизацию или перезарядку на частицах плазмы вдоль траектории пучка. Если Т іше ~ 1—10 кэВ, Е 0 ~ 20— 100 кэВ, то Ѵі ѵъ Ѵе (ѵь — скорость частиц пучка) и с разум­ ной точностью можно считать, что

{ о ех I Ѵі — ѵь I ) = o exv b, ( o t \vt — Vb \) =

= 7iVb,c (O e IVe Vb I>= O eVe ,