ФИЗИКА ПЛАЗМЫ, 2014, том 40, № 8, с. 687-696
МАГНИТНЫЕ ЛОВУШКИ
УДК 533.9.072
ПЛАЗМОСТАТИЧЕСКИЕ МОДЕЛИ ЛОВУШЕК-ГАЛАТЕИ С МАГНИТОПРОНИЦАЕМЫМИ ГРАНИЦАМИ
© 2014 г. К. В. Брушлинский*, **, А. С. Гольдич**
* Институт прикладной математики им. М.В. Келдыша РАН, Москва, Россия ** Национальный исследовательский ядерный университет МИФИ, Москва, Россия e-mail: brush@keldysh.ru, alexdeph@gmail.com Поступила в редакцию 28.01.2014 г.
Представлены математические модели и результаты расчетов равновесия плазмы в круглом цилиндре с погруженными в нее тремя винтовыми или прямыми проводниками с током — распрямленном аналоге тороидальных ловушек-галатей. Равновесие описывается двумерными краевыми задачами с плоским и винтовым аналогами уравнения Грэда—Шафранова для скалярной функции магнитного потока. Рассмотрены задачи с граничными условиями первого рода, которые соответствуют прозрачной для магнитного поля границе цилиндра, и второго рода — в задачах с заданной величиной электрического тока, протекающего по плазме помимо токов в проводниках. В расчетах исследованы деформации магнитоплазменных конфигураций в цилиндре, обязанные указанным разновидностям постановки задач.
DOI: 10.7868/S0367292114080022
1. ВВЕДЕНИЕ
Одно из направлений в проблеме управляемого термоядерного синтеза — разработка и исследование ловушек-галатей, в которых проводники с током, создающие магнитное поле, погружены в плазменный объем. Перспективы таких ловушек и ожидаемая эффективность удержания связаны, во-первых, с большим разнообразием геометрии поля при различных взаимных расположениях проводников и, во-вторых, — с возможностью удерживать проводники в равновесии без дополнительных устройств, компенсируя силу тяжести магнитным давлением. Идеология галатей (и их название) предложены А.И. Морозовым [1], первому этапу исследований посвящен обзор [2]. Примером современных разработок галатей может служить "Тримикс-3М" [3]. Галатеи предполагаются, как правило, традиционной для ловушек тороидальной формы, однако многие теоретические положения и даже некоторые экспериментальные исследования, связанные с ними, возможны и удобны в более простом распрямленном варианте — цилиндре, однородном или периодическом в осевом направлении. В качестве простого примера приведем ловушку "Пояс" [4] — плазменный тор с двумя кольцевыми проводниками внутри и его прямой аналог — цилиндр с двумя проводниками, параллельными оси. Исследования "Пояса" первоначально опирались на опыт работ С.И. Сыроватского [5] с нейтральным токовым слоем. Теоретические [4] и экспериментальные [6] исследования "Пояса" дополнены его математическими моделями:
плазмостатическими [7—9], описывающими возможные равновесные конфигурации в нем, и плазмодинамическими [10—12], связанными с их формированием. К обобщениям и развитию "Пояса" следует отнести "Трилистник" — цилиндр с тремя параллельными токами [8] и его более интересную разновидность — "стелларатор-гала-тею" (СГ) [13], в которой эти три тока являются винтовыми, обвивающими ось цилиндра. Изучение их проведено пока только численно на языке плазмостатических моделей [14, 15, 8], результаты которых могут стать полезными в становлении будущей теории. МГД-модели обсуждаемых равновесных магнитоплазменных конфигураций основаны на численном решении двумерных краевых задач с известным уравнением Грэда—Шафранова [16, 17] в круговом сечении цилиндра. Во всех рассмотренных случаях [8, 14, 15] граница цилиндра предполагается непроницаемой для магнитного поля, т.е. в граничных условиях первого рода искомая функция магнитного потока ¥ считается постоянной, а ее значение определяется конкретными обстоятельствами задачи. Модели при этом имеют относительно простую и замкнутую логику, не зависящую от того, что происходит за границами цилиндра.
Однако в интересах полноты теории и перспектив приложений к экспериментам в более сложных условиях имеет смысл рассмотреть краевые задачи с более сложными граничными условиями: первого рода — с непостоянной по азимуту функцией ¥ (ф) на границе, и второго рода — с заданием нормальной производной д^/ди, которая
(с обратным знаком) равна тангенциальной компоненте магнитного поля на границе. В обоих случаях граница становится прозрачной для магнитного поля, что может быть связано, например, с влиянием дополнительных проводников с током, расположенных вне цилиндра. Кроме того, задание поля и тем самым его циркуляции вдоль границы, означает задание полного электрического тока (в проводниках и в плазме) в ловушке, величина которого раньше определялась в расчетах и зависела от других параметров.
Постановка обоих типов задач и их реализация в расчетах осуществлена в работе [9] применительно к "Поясу". Результаты расчетов продемонстрировали разнообразие равновесных плазменных конфигураций и изменение их параметров в зависимости от магнитного потока сквозь границу и от величины полного электрического тока.
В настоящей работе представлены результаты расчетов краевых задач с упомянутыми граничными условиями в моделях указанных выше ловушек с тремя проводниками в плазменном цилиндре. В рассматриваемых конфигурациях СГ плазма сосредоточена у сепаратрисы магнитного поля, отделяющей окрестности проводников от области в центре и от периферии цилиндра. Они предпочтительнее для удержания плазмы по сравнению с конфигурациями в центральной области [15].
2. ПОСТАНОВКА ЗАДАЧИ И МЕТОД ЧИСЛЕННОГО РЕШЕНИЯ
Математические модели равновесных плазменных конфигураций в ловушках строятся в терминах уравнений магнитной газодинамики, от которых в состоянии покоя (д¡дt = 0, V = 0) остается только одно
Vp = 1 jpl х H, c
(1.1)
¡pi
где ] = у + ] — сумма индуцированного в плазме ("собственного") и заданного в проводнике ("внешнего") токов. Вместе с уравнениями Максвелла
v - H = 0, ^ j = v х H c
(1.2)
они составляют систему трех уравнений с неизвестными давлением р, напряженностью магнитного поля Н и плотностью электрического тока у Распределения этих величин в пространстве ловушки, полученные в результате решения краевой задачи с уравнениями (1.1), (1.2), полностью характеризуют искомую равновесную конфигурацию. Распрямленным аналогом СГ является круглый плазменный цилиндр, в который погружены три винтовых проводника с током. Есте-
ственно ввести дополнительное предположение о винтовой симметрии всей конфигурации: пусть ее характеристики зависят только от двух независимых переменных в цилиндрических координатах (r, ф, z):
r, 0 = ф - az, где а = 2я/h, h — шаг винта проводников.
В результате задача в круговом сечении цилиндра становится двумерной: уравнения (1.1), (1.2) сводятся к одному скалярному уравнению второго порядка для функции магнитного потока ¥ — винтовой разновидности уравнения Грэда—Ша-франова [18, 19, 14]
д **¥ + g(¥) = 0,
(1.3)
где
д**Р = 1 Ml +1 д2p
r dr vn dr J r 30
2 2
П = 1 + a r ,
g (p) = dp - 8na i + ^ = о
c dP cn2 c
Функция ¥ связана с магнитным полем уравнениями
rH = ^, H0 = H - arHz = .
r de 0 ф г dr
(1.4)
Давление p — функция от ¥, т.е. постоянно на каждой магнитной поверхности Т = const. Оно задается в соответствии с какими-либо требованиями к искомой конфигурации и в рассматриваемых задачах принимается вида
Р = Рое 4 ' ' , (1.5)
так что плазма сосредоточена около поверхности
^ = ^ о.
Величина I в уравнении (1.3) — функция электрического тока — положена постоянной I = 3Jc|h, равной отношению суммы токов в трех проводниках к шагу винта в соответствии со сделанными предположением о винтовой симметрии задачи [14].
Чтобы избежать постановки и решения задачи в многосвязной области — круге с вырезанными из него сечениями проводников конечного диаметра, в уравнение (1.3) введено слагаемое
2 (г-го )2+п(8-2ип/3)2
= Г (г, 0) = 7оXе " , (1.6)
п=0
соответствующее току в проводниках. Здесь (г0, 2пп/3) — координаты трех проводников, гс < г0 — условный радиус проводника, а коэффициент ]0 подобран из условия
JJ jex (r, Q)rdrdQ = 3Jc,
(1.7)
l <R
(а)
(б)
Рис. 1. Силовые линии магнитного поля - ^0 = const) и линии уровня давления p/p0 в базовых вариантах конфигурации: а) — в трилистнике ; б) — в СГ.
где Jc — величина тока в каждом из проводников [14, 20, 8].
Краевая задача с уравнением (1.3) в круге r < R включает в себя граничные условия при r = R. В предыдущих работах [14, 15, 8] они соответствовали изолированному кожуху Hn = Hr = 0, т.е. согласно (1.4), ¥ = ¥ г = const. Значение константы может быть произвольным, в расчетах она положена равной нулю, а параметр ¥ 0 в формуле (1.5) подбирается так, чтобы он соответствовал сепаратрисе, отделяющей окрестности проводников от центральной области и периферии ловушки (рис. 1). В интересах дальнейшего развития теории ловушек-галатей в настоящей работе рассматриваются два типа более сложных граничных условий:
1) Условие первого рода с заданной при г = Я функцией
Ч = Ч Г(0), (1.8)
оно позволяет магнитному полю проникать сквозь границу цилиндра, а именно,
d Ч.
Hn = Hr = ф 0. n r rdQ
(1.9)
= -Hr(6),
(1.10)
2) Условие второго рода с заданной нормальной производной искомой функции ¥
дЧ =
ди дг
где НГ = Н0(Я, 0) — "азимутальная" в винтовых координатах компонента магнитного поля. Как известно, вторая краевая задача с уравнением Пуассона разрешима только при выполнении необходимого согласования потока через границу с
интенсивностью источника. В рассматриваемой задаче аналогичное условие возникает, если проинтегрировать уравнение (1.3) по кругу г < r и использовать граничное условие (1.10):
П(К)
| HГ дое= Л g(x¥)гdгdе.
(1.11)
0 г <к
Равенство (1.11) при заданной функции ИГ(0) вместе с уравнением (1.3) однозначно определяют искомую функцию ¥ и максимальное значение давления p0 в формуле (1.5).
Поскольку электрический ток Jc во всех проводниках один и тот же, а их расположение симметрично относительно осей 0 = п п/3 , поставленная задача по физическому смыслу также симметрична относительно этих осей. Поэтому ее достаточно ставить и решать только в секторе
п
0 < г < ^ 0 < 0 <
3
(1.12)
дополнив очевидными условиями симметрии
(1.13)
— = 0 при 0 = 0, п.
5
Для дальнейшего прочтения статьи необходимо приобрести полный текст. Статьи высылаются в формате PDF на указанную при оплате почту. Время доставки составляет менее 10 минут. Стоимость одной статьи — 150 рублей.