Амбиполярная ловушка является одной из возможных систем с магнитным удержанием плазмы, на основе которой можно реализовать реактор для проведения управляемого термоядерного синтеза. При этом, будучи открытой системой, пробкотрон обладает рядом преимуществ по сравнению с замкнутыми системами, в которых силовые линии магнитного поля не пересекают границ плазмы: в нем эффективно используется удерживающее плазму магнитное поле (достижимы значения? ~ 1 — отношение давления плазмы к давлению магнитного поля рм = В2/&п), он может работать в стационарном режиме (нет проблемы накопления продуктов реакции), он мало чувствителен к накоплению примесей, он проще замкнутых систем в топологическом отношении [1−3]. Центральные ячейки существующих амбиполярных ловушек представляют собой соленоиды с осесимметричным полем, с длиной, гораздо большей длины концевых ячеек или пробок, которые имеют квадрупольные магнитные катушки, обеспечивающие макроскопическую устойчивость. Эти концевые ячейки имеют структуру с минимумом магнитного поля, которая обладает благоприятной кривизной силовых линий и компенсирует дестабилизирующий вклад от остальных участков ловушки.
В результате многолетних экспериментальных исследований были достигнуты значительные успехи на крупных амбиполярных ловушках, таких как TMX-U (Livermore, USA) [4,5], GAMMA-10 (University of Tsukuba, Japan) [6,7], изучены различные способы улучшения удержания плазмы и увеличения ее параметров. Экспериментальные результаты, полученные на крупнейших современных амбиполярных ловушках, и перспективы создания МГД устойчивой, осесимметричной ловушки с высоким ?3 подробно описаны в обзоре [8]. Достигнута МГД устойчивость плазмы в тандеме ловушек с неустойчивой геометрией центрального пробкотрона (соленоида). В концевых проб-котронах, пристыкованных к соленоиду, были сформированы амбиполярные барьеры для продольного удержания ионов в соленоиде (ионные барьеры) до 1,5 кВ в ТМХ-11 (1984) до 2 кВ в САММА-Ю (1992). Продольное время удержания частиц амбиполяр-ными барьерами было существенно увеличено по сравнению с удержанием магнитными пробками в обычном пробкотроне Будкера-Поста. На ТМХ-11 было достигнуто 50-кратное, а на САММА-Ю — 1000-кратное увеличение продольного времени удержания. В то же время, плотность амбиполярно удерживаемой плазмы в соленоиде была ограничена из-за относительно небольшой плотности горячих ионов в ионных барьерах и из-за низкого предела величины (3 по МГД устойчивости в центральном соленоиде САММА-Ю на уровне 0,5%. Исследовались микронеустойчивости в концевых пробко-тронах в диапазоне ионно-циклотронных частот. Были подробно изучены популяций горячих электронов, формируемые в концевых пробкотронах ЭЦР-нагревом. Изучаются различные способы контроля потоков электронов из ловушки. В центральном пробкотроне исследовались поперечные неоклассические потери плазмы с контролем через радиальное поле в плазме. Исследуется стабилизация низкочастотных колебаний путем формирования заданного радиального профиля потенциала подачи с помощью радиально секционированных торцевых проводящих пластин.
Тем не менее достигнутые результаты пока далеки от максимальных возможностей открытых ловушек. Одним из существенных препятствий на пути увеличения параметров плазмы было то, что для обеспечения крупномасштабной магнитогидродина-мической (МГД) устойчивости плазмы применялись неосесимметричные магнитные поля. Как выяснилось [9,10], они могут легко становиться источником усиленных поперечных потерь плазмы, которые становятся особенно ощутимыми в условиях, когда продольные потери подавлены. Со временем стало ясно, что в амбиполярных ловушках с мультипольными ячейками, обеспечивающими «средний минимум В», необходимый для МГД устойчивости, очень трудно ограничить эти неоклассические потери термоядерной плазмы. Кроме того, при сооружении осесимметричных ловушек существенно упрощаются инженерные проблемы и появляется возможность значительно повысить магнитное поле в пробках. Это стимулировало интенсивный поиск МГД устойчивых осесимметричных конфигураций. Хорошо известно, что в параксиальной геометрии плазма низкого давления ((3 1) с естественной функцией распределения по энергии (т.е. / является монотонно убывающей функцией энергии д//де < 0) неустойчива по отношению к желобковым возмущениям. Возможности стабилизации связаны с отказом от параксиальности, переходом к высоким давлениям, (3 ~ 1, и использованием анизотропных распределений частиц по скоростям.
В 80-е годы был предложен ряд способов МГД стабилизации и выполнено значительное количество теоретических работ и ряд экспериментов, которые описаны, например, в обзорах [11,12]. Однако большинство предложенных способов имеют ограничения, в частности, по величине /3, или пригодны только для полой структуры плазмы, с относительно небольшой толщиной плазмы. Выяснилось, что, по-видимому, нет простых магнитных конфигураций, обеспечивающих устойчивость основного объема плазмы с естественным распределением по энергии и при любом давлении, от /3 = 0 до ?3 ~ 1. В качестве основного метода МГД стабилизации высокотемпературной плазмы в осесимметричной амбиполярной ловушке предлагается стабилизация проводящим кожухом в комбинации с эффектом конечного ларморовского радиуса (КЛР) при /3 ~ 1. При этом допускается существование значительного зазора между плазмой и проводящей стенкой. Этот метод базируется на стабилизации плазмы горячими частицами, в первую очередь, электронными кольцами и дисками [12,13].
Для изучения возможности получения плотной МГД устойчивой плазмы в осесимметричной геометрии была спроектирована установка АМБАЛ-М. Установка АМБАЛ-М изображена на рис. 0.1 а ее подробное описание представлено в [14,15].
Рис. 0.1: Установка АМВАЛ-М (вид сверху). 1 — центральный пробкотрон, 2 — один концевой пробкотрон, 3 — полукасп, 4 — кольцевой приемник плазмы, 5 — горизонтальные атомарные инжекторы, 6 — поглотители пучков.
На ней планируется осуществить следующие экспериментальные задачи:
• МГД стабилизация высокотемпературной плазмы в осесимметричной геометрии.
• приближение потерь плазмы в концевых пробкотронах к классическим.
• реализация длительного поддержания термобарьеров для электронов соленоида.
• исследование поперечных потерь плазмы в длинном осесимметричном соленоиде.
• повышение параметров высокотемпературной плазмы в соленоиде: пс > 1013 см-3, Тс ~ 1 кэВ, ионный барьер Д<£>" ~ 5 кВ, пт ~ 1013 см3сек.
Сооружение установки было запланировано таким образом, чтобы осуществить ее запуск в две очереди. В этой работе мы не будем упоминать все известные схемы стабилизации плазмы в осесимметричной геометрии, а рассмотрим лишь те, которые уже используются на концевой системе АМБАЛ-М и те, которые будут применяться на полномасштабной установке АМБАЛ-М [14,16]. В первую очередь запущена концевая система, которая представляет собой концевой пробкотрон и полукасп в качестве МГД якоря. На этой системе запланировано провести основные эксперименты по МГД стабилизации плазмы. Затем планируется установить 2-ю концевую систему и соленоид, осуществить запуск установки в целом и выполнить всю намеченную программу экспериментов.
На 1-ой очереди установки предусмотрена МГД стабилизация плазмы в концевом пробкотроне двумя способами: с помощью полукаспа [16−18] и популяцией горячих электронов с высоким /?, создаваемой в средней плоскости пробкотрона внутри проводящего кожуха [13,19]. МГД устойчивость плазмы к поперечному смещению в целом (по глобальной моде) обеспечивается, в соответствии с энергетическим принципом, вкладом непараксиального члена, который в полукаспе является существенным, в потенциальную энергию возмущения. Давление плазмы в полукаспе поддерживается с помощью атомарной инжекции, производимой перпендикулярно к оси системы, в промежутке между осевой пробкой полукаспа и поверхностью минимумов магнитного поля в полукаспе. МГД устойчивость плазмы по высоким перестановочным модам обеспечивается эффектом конечного ларморовского радиуса. Одновременно со стабилизацией МГД якорем планируется осуществить электронно-циклотронный резонансный нагрев в пробкотроне. Предусматривается создание цилиндра горячих электронов с энергией до 400 кэВ, радиусом около 10 см и длиной ~ 40 см. Эта анизотропная популяция с высоким (3]н (отношение давления горячих электронов к давлению магнитного поля), должна обеспечивать стеночную стабилизацию для низших перестановочных мод т > 1. Стартовая плазма в концевой системе создается с помощью газоразрядного источника плазмы с кольцевым каналом [20]. В 1-ой очереди предусмотрено использовать кольцевые газоразрядные ячейки на приемниках плазмы для дополнительной МГД стабилизации плазмы во время старта за счет «вмороженности» в торцы [16].
В настоящее время в концевом осе симметричном пробкотроне установки АМБАЛ-М с магнитным полем 7,5 кГс в центральной плоскости (рис. 0.2) получена горячая стартовая плазма, которая имеет электронную температуру Те ~ 50 эВ, ионную температуру ~ 200 эВ и плотность около 1013 см-3 в центральной плоскости ловушки. Горячая стартовая плазма создается с помощью струи плазмы от газоразрядного источника, расположенного в области перед входной пробкой.
Механизм создания горячей стартовой плазмы и общие экспериментальные результаты описаны в [22,23]. В данной работе рассматриваются экспериментальные результаты по МГД устойчивости горячей стартовой плазмы в концевой системе АМБАЛ-М, производится их сравнение с численными оценками и предлагается модель, которая может объяснять наблюдаемую устойчивость плазмы при отключенном полукаспе. Для МГД устойчивости в осесимметричных открытых ловушках надо соблюдать довольно жесткие ограничения на параметры плазмы и их распределения в различных секциях установки. Поэтому каждая установка с новой конфигурацией магнитного поля дает полезную и необходимую экспериментальную информацию о реальной ра.
1 5 4.
Рис. 0.2: Схема концевой системы АМБАЛ-М и силовая линия, выходящая из источника плазмы. Ниже — магнитное поле на оси. 1 — катушки пробкотрона, 2 — соленоид источника плазмы, 3 — источник плазмы, 4 — плазмоприемник, 5 — катушка полукаспа, 6 и 7 — входная и выходная пробки, соответственно. Стрелками на рисунке отмечены координаты сечений, в которых производились измерения магнитным зондом. ботоспособности и эффективности заложенных в ее конструкцию принципов стабилизации плазмы. АМБАЛ-М позволяет изучать поведение бесстолкновительной плазмы в пробкотроне при стабилизации полукаспом и с помощью описанных выше популяций горячих частиц. Ввиду важности МГД устойчивости для успешной реализации намеченной программы исследований мы уделили этому вопросу особое внимание, провели ряд экспериментов по изучению поведения плазмы в разных условиях и сделали анализ полученных данных. Хотя проделанное исследование нельзя назвать полным, мы получили удовлетворительное описание поведения стартовой плазмы в рамках выбранных моделей для анализа устойчивости.
МГД устойчивость плазмы в системе осесимметричных пробкотронов была рассмотрена в [18], где было показано, что в такой составной системе все азимутальные моды кроме т = 1 (смещение плазменного шнура как целого) могут быть стабилизированы эффектами КЛР ионов, а мода то = 1 стабилизируется с помощью непараксиальной осесимметричной ячейки. Непараксиальная ячейка в виде каспа была подробно изучена в работе [17] в связи со стабилизацией моды т = 1, где изучалось влияние геометрии магнитного поля и распределений давления плазмы на устойчивость. В настоящей работе мы используем этот подход для вычисления энергии возмущений для первой моды, используя реальную геометрию магнитного поля и радиальное распределение давления плазмы, измеренное в полукаспе.
Для изучения стабилизирующей роли полукаспа были проведены эксперименты с его отключением. При этом оказалось, что температура, плотность, диамагнетизм плазмы практически не изменились. Время распада плазмы после выключения источника тоже было примерно такое, как при включенном полукаспе. Эти наблюдения указывают, что-либо характерное время развития желобковой неустойчивости одного порядка или велико по сравнению с измеренным временем распада плазмы, либо есть какой-то дополнительный механизм стабилизации плазмы пробкотрона даже при отключенном полукаспе. Для объяснения устойчивости были сделаны некоторые оценки, которые не противоречат данным наблюдений. Построение точной модели поведения плазмы и выяснение причин ее устойчивости в данном случае затруднительно из-за наличия нескольких возможных стабилизирующих факторов, учет каждого из которых требует тщательного и длительного экспериментального исследования. Поскольку МГД неустойчивость не проявилась и не стала препятствием для проведения экспериментов по изучению накопления и нагрева стартовой плазмы в пробкотроне, то несмотря на то, что сам факт устойчивости при отключенном полукаспе был безусловно любопытен, мы не стали подвергать его детальному изучению, а ограничились информацией от небольшого числа диагностик и разумными численными оценками явления. При отключении полукаспа система становится параксиальной, и подобные конфигурации ранее были изучены во многих теоретических и численных работах. Мы используем результаты работы [21], в которой получено уравнение устойчивости с учетом эффектов КЛР, проводящей боковой стенки и «вмороженности» силовых линий через холодное плазменное гало. Подставляя измеренные профили плотности плазмы, температуры и электрического потенциала и учитывая реальную геометрию магнитного поля в установке, мы решали это дифференциальное уравнение второго порядка, усредненное вдоль силовых линий, чтобы найти собственные частоты и собственные функции радиального смещения, и обнаружили, что в условиях нашего эксперимента неустойчивость может быть сильно ослаблена, если плазма окружена идеально проводящим плазменным гало.
Одним из важных вопросов для удержания плазмы в открытой ловушке является ограничение поперечных потерь. Эта проблема становится особенно актуальной при использовании мощного ИЦР и ЭЦР нагрева, а также при создании плазмы в длинном центральном соленоиде. В настоящее время эти пункты программы АМБАЛ-М еще не реализованы полностью. Проводятся первые серии экспериментов по ИЦР нагреву на частоте 10 МГц с помощью антенны, установленной в переходной области между выходной пробкой концевого пробкотрона и полукаспом. Исследуется структура возникающих при нагреве высокочастотных полей в плазме, изучаются различные режимы нагрева: при напуске газа в различные участки установки, при включении нагрева до и после работы газоразрядного источника плазмы, при различной длительности ИЦР нагрева вплоть до квазистационарного режима с длительностью до 80 мс [24].
Тем не менее стартовая плазма, создаваемая газоразрядным источником в концевом пробкотроне, также представляет собой интересный объект в связи с изучением поперечного переноса частиц, в частности, потому, что в струе плазмы, заполняющей пробкотрон из источника, наблюдаются значительные колебания плотности и потенциала (до 100%), что может приводить к существенному турбулентному потоку плазмы поперек силовых линий магнитного поля. Классический поперечный поток частиц, рассчитанный по теории парных столкновений, сильно ослабляется магнитным полем, Da± ~ Da\/(ujcaTa)2 ос 1/.В2, где Da±, Daц — коэффициенты диффузии частиц сорта «а» поперек и вдоль магнитного поля, соответственно, шса — циклотронная частота и та — время между столкновениями частиц одного сорта друг с другом [25]. В реальных экспериментальных условиях шса та 1, и наблюдаемые поперечные потоки гораздо больше потоков, оцениваемых по классическим формулам. Это связано с тем, что даже небольшой уровень флуктуаций, всегда присутствующий в плазме, приводит к существенному увеличению поперечного переноса, так что классическим поперечным переносом, который оказывается на несколько порядков меньше турбулентного, можно полностью пренебречь.
Механизм возникновения аномальных переносов можно понять из следующего рассмотрения [26]. Плотность и температура плазмы, находящейся в магнитном поле, как правило, неоднородны в пространстве. Наличие пространственных градиентов поперек магнитного поля, либо присутствие электрического поля приводят к возникновению в плазме дрейфовых токов. Частицы, поддерживающие эти токи, движутся в направлении, перпендикулярном направлению среднего магнитного поля и градиента плазмы (в направлении Уп х В, УТ х в или У^ х в). Если такое дрейфовое движение оказывается неустойчивым, то возникает движение частиц поперек магнитного поля. Смещения частиц поперек поля имеют компоненту вдоль направления градиента плазмы и могут во много раз превышать их гирорадиус. В результате участия частиц плазмы в коллективных колебаниях системы возникает эффект переноса плазмы вдоль направления макроскопического градиента. Можно сказать, что в макроскопическом масштабе проявление плазменных эффектов связано с возникновением коллективной диссипации в форме аномальной диффузии частиц (тепловой энергии) вдоль градиента плотности плазмы (градиента температуры). Коллективный перенос не зависит от парных соударений и по величине может превышать столкновительный на много порядков. Явления переноса возникают также в случае, когда какие-либо внешние источники возбуждают в плазме электрические или магнитные поля, по свойствам близкие к коллективным (например, при высокочастотном нагреве). При достаточно малых давлениях аномальный перенос в основном определяется механизмом, связанным с коллективными электрическими полями. Эти низкочастотные электрические поля можно считать электростатическими (т.е. потенциальными), что существенно упрощает проблему. Аномальный перенос плазмы на открытых ловушках АМБАЛ-М и МАЛ, работающих в ИЯФ СО РАН (г.Новосибирск), связанный с электростатической турбулентностью, развивающейся в струе плазмы, вытекающей из газоразрядного источника с кольцевым каналом, рассматривался в [27,29]. Было показано, что колебания плотности плазмы и азимутального электрического поля сильно коррелированы друг с другом, и радиальный поток плазмы, определенный по этим измеренным флук-туациям, соответствует пространственному распределению плотности в ловушках.
В экспериментах на концевой системе открытой ловушки АМБАЛ-М был обнаружен продольный электронный ток величиной до 1,6 кА [32,33], протекающий в стартовой плазме, и перенос плазмы поперек магнитного поля с коэффициентом диффузии, на два порядка превышающим классический. Эффект протекания продольного электронного тока по плазме при использовании источника плазмы с кольцевым разрядным каналом наблюдался и раньше на установках МАЛ и АМБАЛ-Ю [34]. В этих экспериментах проводились измерения тока поясами Роговского и было выяснено, что на внутренних радиусах плазменного столба электронный ток протекает от источника к противоположному торцу установки, а на периферийных радиусах имеет обратное направление. Такая картина распределения тока по сечению плазмы обусловлена тем, что часть тока разряда выносится из источника вдоль силовых линий и замыкается поперек силовых линий в плазме и на плазмоприемнике. Этот ток регистрируется на всей длине установки, поглощается торцом установки, затем по вакуумной камере и по периферии плазмы возвращается на анод источника плазмы.
Измерения радиального распределения плотности тока в нескольких сечениях магнитным зондом и плоским двусторонним ленгмюровским зондом во входной пробке пробкотрона показали, что вблизи источника плазмы профиль плотности тока имеет кольцевую структуру, а затем заполняет внутреннюю полость по мере удаления от источника.
Интерес к изучению магнитных флуктуаций не только снаружи плазмы, но и внутри нее возник после того, как теоретически было показано, что магнитные флуктуации могут приводить к существенным потокам частиц и тепла поперекмагнитного поля из-за дрейфа частиц вдоль возмущенных флуктуациями магнитных поверхностей. Появилось много экспериментальных работ, посвященных прямому измерению флуктуаций магнитного поля, изучению их спектрального и модового состава (см., например, [35−37]). В работе [35] описан один из первых экспериментов по контактному измерению флуктуаций магнитного поля в плазме токамака с помощью миниатюрных многовитковых катушек, помещенных в изолирующую оболочку. Применимость таких контактных диагностик ограниченна, поскольку вводимые в высокотемпературную плотную плазму зонды либо разрушаются под действием потоков горячих частиц, либо вызывают неприемлемые возмущения плазмы. Однако в относительно холодную (Те ,$ 100 эВ), неплотную (п ^ 1013 см-3), а тем более импульсную плазму магнитные зонды, помещенные в изолирующие корпуса, могут быть легко введены без заметного возмущения ее параметров и без угрозы разрушения зондов под воздействием плазмы.
На различных установках (на токамаках, стеллараторах, пинчах с обращенным полем) проводятся эксперименты по прямому измерению потоков частиц и тепла путем определения корреляций между определенными параметрами плазмы, такими как плотность, продольный ток, электрическое и магнитное поле и т. д. (см. обзоры [38−40], а также работы [41−43]). Насколько нам известно, на открытых ловушках такие подробные локальные измерения не проводились. Значительный ток, текущий в плазме АМБАЛ-М, вызывает возмущения магнитного поля, что может приводить к поперечному переносу частиц вследствие их движения по деформированным магнитным поверхностям. Поэтому такие измерения представляют для нашей установки особый интерес в связи с определением радиального потока и изучением пространственных и спектральных свойств магнитных флуктуаций, которые могут дать важную дополнительную информацию о процессах, протекающих в плазме АМБАЛ-М. Таким образом, возникло естественное желание произвести непосредственные локальные измерения радиального потока частиц, вызванного флуктуациями магнитного поля. Магнитные флуктуации станут особенно существенными при повышении /3 до величин, когда давление плазмы станет сравнимо с давлением вакуумного магнитного поля. Тогда равновесие плазмы будет устанавливаться самосогласованным образом и ситуация будет сильно отличаться от той, что реализована в настоящее время, когда плазма незначительно искажает вакуумное магнитное поле, и ее влияние может рассматриваться как возмущение. Поэтому при увеличении ?3 роль магнитных флуктуации должна возрастать, и потребуется их тщательное исследование.
В экспериментах по изучению свойств продольного тока, магнитных флуктуаций и определения вклада последних в поперечный перенос плазмы использовались подвижные магнитные зонды, измеряющие интересующие компоненты магнитного поля. Анализ спектров и радиальных профилей флуктуаций магнитного поля и плотности показал, что между этими величинами есть корреляция, которая объясняется наличием неустойчивости Кельвина-Гельмгольца (НКГ) в струе плазмы. Выполнены оценки параметров характерных неоднородностей тока в струе.
Диссертационная работа состоит из Введения, пяти глав и Заключения. В 1-й главе описана концевая система установки АМБАЛ-М и диагностики, которые используются для измерения параметров плазмы. 2-я глава посвящена оценке запаса МГД устойчивости, который обеспечивает полукасп для 1-й азимутальной моды колебаний в системе «пробкотрон + полукасп». В 3-й главе рассматривается МГД устойчивость стартовой плазмы в концевой системе при отключенном МГД стабилизаторе (полука-спе), предлагается возможный механизм, обеспечивающий устойчивость, приводятся результаты численных расчетов и обсуждаются недостатки данной модели. В 4-й главе описаны результаты измерений и численных оценок флуктуаций магнитного поля в установке, связанных с заполнением ловушки плазмой и с протеканием значительного продольного тока. Эти результаты используются в 5-й главе для определения радиаль.
ЗАКЛЮЧЕНИЕ
.
Сформулируем основные результаты диссертационной работы:
1. Измерены параметры горячей плазмы в полукаспе. С использованием экспериментальных данных определен запас устойчивости концевой системы, обеспечиваемый полукаспом.
2. На осонове экспериментально полученных данных показано, что плазма в проб-котроне при отключенном полукаспе устойчива. Проанализированы возможные причины устойчивости, и для ее объяснения предложена модель стабилизации плазмы в пробкотроне холодной периферийной плазмой.
3. Обнаружены и измерены флуктуации магнитного поля в стартовой плазме. Предложен механизм возникновения этих флуктуаций, и выполнены оценки, объясняющие их величину.
4. Выполнена оценка роли магнитных флуктуаций в поперечном переносе плазмы в концевой системе. Измерен радиальный профиль поперечного потока частиц, связанного с магнитными флуктуациями, и определен коэффициент диффузии.