ФИЗИКА ПЛАЗМЫ, 2014, том 40, № 7, с. 575-590
МАГНИТНЫЕ ЛОВУШКИ
УДК 533.9
НЕКОТОРЫЕ ЗАДАЧИ ОСЕСИММЕТРИЧНОГО ТОРОИДАЛЬНОГО РАВНОВЕСИЯ В ОРТОГОНАЛЬНЫХ ПОТОКОВЫХ КООРДИНАТАХ
© 2014 г. В. В. Арсенин, А. А. Сковорода
Национальный исследовательский центр "Курчатовский институт", Москва, Россия
e-mail: Arsenin_VV@nrcki.ru Поступила в редакцию 05.11.2013 г. Окончательный вариант получен 19.02.2014 г.
Даны постановки задач осесимметричного равновесия в ортогональных потоковых координатах с граничными условиями, не предполагающими фиксацию формы границ, как для ловушек типа ле-витрона с омываемым плазмой проводником, так и для конфигураций без проводника, включая то-камак. Приводятся иллюстративные примеры численного решения этих задач. В случае с проводником демонстрируется, как, используя потоковые координаты, можно найти геометрию равновесия с изодинамической магнитной поверхностью.
DOI: 10.7868/S0367292114070014
ВВЕДЕНИЕ
В ряде задач об устойчивости плазмы и о переносах, где важно поведение магнитного поля в окрестности магнитной поверхности, естественно при описании равновесия пользоваться системой координат, одной из которых служит именно метка магнитной поверхности (магнитные, иначе потоковые, см. [1, 2], или "обращенные" [3, 4] переменные). Хотя методам теории равновесия, использующим потоковые координаты, посвящена обширная литература, их развитие применительно к различным системам удержания и геометриям продолжается. Так, в недавней работе [5] (к ней и отсылаем за указаниями на многочисленные публикации, включая обзоры, по этим методам) улучшено описание области вблизи сепара-трисной поверхности в токамаке. В настоящей статье мы развиваем один из предложенных ранее подходов — подход с использованием ортогональных потоковых координат [3, 4], имея в виду, в частности, отыскание в рамках этого подхода конфигураций с изодинамическими магнитными поверхностями, вдоль которых модуль магнитного поля не изменяется. Такие поверхности привлекают внимание тем, что через них нет магнитного дрейфа частиц. Существуют осесимметрич-ные тороидальные равновесия, в которых изодинамичны все магнитные поверхности в плазме [6]. Однако в них на магнитной оси не только полоидальное, но и тороидальное поле, а также величина запаса устойчивости q обращаются в нуль; равновесие неустойчиво [6, 7]. Могут представлять интерес конфигурации с конечным q, в которых изодинамичность достигается на одной магнитной поверхности [8]. Как мы убедимся, использование обращенных ортогональных
переменных вполне адекватно задаче нахождения равновесий с изодинамической поверхностью в случае, когда имеется токонесущий внутренний тороидальный проводник.
Как видно из [5], разнообразие постановок задач равновесия в потоковых координатах со временем расширяется. В нашей работе изложены несколько нетипичных постановок задач осесимметричного равновесия с использованием ортогональных координат и даны иллюстративные примеры их решения, в том числе относящиеся к ловушкам с изодинамической магнитной поверхностью. Перед этим мы приводим в разд. 1 вывод самой системы уравнений равновесия в обращенных переменных, основанный на простых геометрических соотношениях и, как нам представляется, более наглядный, чем в [3, 4]. Следующий раздел посвящен обсуждению граничных условий. Они определяют вид конфигурации. В разд. 3 даны примеры расчетов при различных граничных условиях. Рассмотрены геометрии с внутренним проводником и без него, когда магнитная ось находится в плазме (как в токамаке). В случае с проводником строятся конфигурации с изодина-мической магнитной поверхностью. Заключительный раздел содержит обсуждение ряда важных особенностей наших постановок задач. В Приложении 1 приведен для полноты еще один, отличающийся от [3] и разд. 2 нашей работы (и более громоздкий, чем в разд. 2), вывод уравнений равновесия в потоковых координатах, исходя из употребительных потокового и токового представлений магнитного поля. Приложение 2 раскрывает детали отыскания приближенных аналитических решений вблизи магнитной оси, используемых в качестве стартовых в расчетах
(a)
х = 0
V = V0
: c°nst_
(в)
0
z/Ro 0.05
0.04
0.03
0.02
0.01
0
0.02
0.04
z
z
n
r/R0-1
Рис. 1. Ортогональные координаты у, х в меридиональном сечении тора (а); сетка линий у = const, х = const около магнитной оси в случае эллиптической формы сечений (б); углы, фигурирующие в формулах (1.22)-(1.24) (в); сетка неортогональных координат у, О в случае эллиптических сечений (г).
разд. 3. В Приложении 3 рассматривается влияние неизодинамичности магнитной поверхности на движение частиц около нее.
1. УРАВНЕНИЯ ОСЕСИММЕТРИЧНОГО РАВНОВЕСИЯ В ПОТОКОВЫХ КООРДИНАТАХ
Набор координат составляют потоковая переменная у — деленный на 2 п поток полоидального поля B (d у = Bprdn, где r — расстояние от главной оси, dn — элемент длины по нормали к магнитным поверхностям у = const), тороидальный угол q и полоидальная координата х, аналог обычного полоидального угла; в случае ортогональных координат, когда линии х = const в меридиональном сечении тора искривлены, см. рис. 1, х не совпадает с полоидальным углом. Величина у растет от значения у = 0 на поверхности тора до у = у 0 на магнитной оси, которая считается единственной, или до у = уc на поверхности омываемого плазмой токонесущего кольца, если таковое присутствует; в последнем случае плазма
трубчатая. Будем отсчитывать переменную х по часовой стрелке от внутреннего обвода тора и в большей части выкладок полагать ее изменяющейся от 0 до 2 п.
Система уравнений равновесия в потоковых координатах включает три уравнения. Два из них описывают связь этих координат с цилиндрическими r, z и имеют геометрическое происхождение. Третье есть собственно уравнение силового баланса. Уравнения равновесия в ортогональных координатах у, с, х получены, исходя из уравнения Грэда—Шафранова [9, 10], в работе [3], другой вывод — с непосредственным использованием геометрических соотношений — приведен в [11] (а для ловушек с чисто полоидальным полем был дан, включая ситуацию анизотропного давления плазмы, ранее в [12]). В общем случае, когда поверхности у = const и х = const могут быть неортогональными, уравнения равновесия выписаны в работе [13] и обзорной статье [4]. Мы проследим (дополняя тем [13, 4]) геометрический смысл и этих уравнений.
Напомним сначала, как нужные уравнения получаются в более простом случае ортогональ-
ных координат. Согласно рис. 1а, А г / A l = = 8r/8n = sin a, Ar/ Al = -8г/8n = cos a, где Al — длина отрезка (с проекциями Ar, Аг) силовой линии полоидального поля между близкими координатными поверхностями х = const, 8n — расстояние между близкими магнитными поверхностями у = const по нормали, a — угол между вектором Bp и радиальным направлением. Отсюда следуют уравнения для r(y, х) и г(у, х)
Gr dr = дг
ду
дг
ду
Gr ^ = -
дХ dr
где
G = B,
81 ;5х.
(1.1) (1.2)
(1.3)
Через О выражается якобиан преобразования кол — 1 2 ординат / = Вх д1 /дх = О/Вх. Расписанное в координатах у, с, х условие силового равновесия
—(Vx B) х B = Vp 4п
(1.4)
(коэффициенты Ламе равны Hv = (Bxr) 1, Hq = Hx = B7,lG) дает уравнение для G(y, х)
r,
G dG + 4п ду
или
F dF
dp
d y 4nr2dу
dr
dX.
+
= 0, (1.5)
dG
dp
+ -
F dF I 2 r X
■ + 4n
Gdy [ d y 4nr¿ d у
^ dr Л + ( дг^
vdyJ [dy
(1.6)
= 0.
Здесь р(у) — давление плазмы, а функция В (у) связана с тороидальным полем: Вд = В/г. При этом плотность тороидального тока в плазме
i_ =±-(Vx B)c = crl +
da F dF
(1.7)
4n \dy 4nr dyy
(c — скорость света), плотность полоидального
тока
ix = -Вг —. z 4я d у
Для полоидального поля имеем
(1.8)
j)2 _ т>2 _
Bp - Вх -
G2
(dr/дх)2 + (дг/дх)2
G2
(dr/dx)2 + G 2r 2(дr/ду)2
r 2[(дг/ду)2 + (dr / ду)2]
(1.9)
(1.10)
а для его г, ^-компонент
в =-в 2г -дг =_О_¿Т
г х ду (дг/дх)2 + (дг/дх)2 дх'
Вг = вХг =--О-2 .
г ду (дг/дх)2 + (дг/дх)2 дх Функции р(у) и В (у) считаются заданными.
Тройка (1.1), (1.2), (1.5) — это система уравнений равновесия в ортогональных потоковых координатах. К ней можно также прийти более сложным путем, используя формально технику магнитных координат [2], см. Приложение 1. Отметим, что силовые линии В в рассмотренных координатах не выпрямлены.
Вместо (1.2) можно использовать эллиптическое уравнение для г(у, х)
д(Ог ^) + А.(± &| = 0,
. . . , , (1.11) ду^ дуу dxv Gr дху вытекающее из (1.1), (1.2).
Хотя, в отличие от уравнения Грэда-Шафра-нова в цилиндрических координатах, уравнение силового равновесия (1.5) в обращенных переменных — это уравнение первого порядка, вся система уравнений равновесия имеет в этих переменных более высокий порядок.
Через величину G на поверхности у = const выражается значение тороидального тока I внутри этой поверхности. Пусть плазма располагается в слое 0 < у < у с. Тогда
- Ic _ ^i^dS _ jjiqdndl _
Vc
=ÍÍ
; dy Gdx _
X _
(1.12)
«ÍÍ
F dF \ G
dp +_
d y 4nr2 dy) B2
d yd x,
где 1С — ток в торе у = у С. Если границей у = у С служит сверхпроводник, то 1С — это ток, текущий по его поверхности. Принимая во внимание (1.9) и (1.6), получаем
В -1- = -
— Í [—dydх 4п ^ J dy
(1.13)
— 4п
cf [G(y, X) - G(yc, X)]dх.
С другой стороны, I можно записать через циркуляцию полоидального поля
I s(¥) = — Í Bdl = — ff Bx ^d х =
4п 1
4n¿ дх
— 4п
cf G(y, x)dx
(1.14)
V
V
¥
1
Из (1.1З), (1.14) следует
G I 0 не зависит от х, когда координатой х слу-
Ic =
4п
(J G(v, х)
V^Vc
dX
(l.l5)
В случае, когда проводника с током Ic нет, у магнитной оси у = у 0 должно быть
fx)!v^VodX^ 0,
и значит, ввиду положительности функции G (1.3), она должна удовлетворять требованию
GUVo^ 0 (1.16)
на всем интервале х.
В этом случае имеется связь между величинами у 0, Iq(0) и большим радиусом шнура R. Пусть, например, шнур тонкий (аспектное отношение R/a велико), имеет сечение, близкое к круглому с радиусом a , и плотность тороидального тока почти постоянна по сечению. Тогда у0 « Iq(0)R/c, а G(0, х) « const = 2у 0/R. При иных распределениях тока и форме сечения отношение cy0/(Iq(0)R) не будет, вообще говоря, близким к едини
Для дальнейшего прочтения статьи необходимо приобрести полный текст. Статьи высылаются в формате PDF на указанную при оплате почту. Время доставки составляет менее 10 минут. Стоимость одной статьи — 150 рублей.