Баллонные моды и энергобаланс в токамаке с адиабатическим сжатием тема автореферата и диссертации по физике, 01.04.08 ВАК РФ

ВВЕДЕНИЕ.

ГЛАВА I. АНАЛИТИЧЕСКИМ КРИТЕРИИ УСТОЙЧИВОСТИ БАЛЛОННЫХ МОД В ТОКАМАКЕ С НЕКРУГЛЫМ СЕЧЕНИЕМ

§ I.I. Исходные уравнения для идеальных баллонных мод в токамаке

§ 1.2. Метод эквивалентных гармоник для исследования высоких мод баллонной неустойчивости

§ 1.3. Локальный критерий устойчивости для круговых магнитных поверхностей

§ 1.4. Локальный критерий устойчивости для некруговых магнитных поверхностей

§ 1.5. Анализ критерия устойчивости

§ 1.6. В ы в о д ы.

ГЛАВА П. РАДИАЛЬНАЯ СТРУКТУРА ВЫСОКИХ МОД . БАЛЛОННОЙ НЕУСТОЙЧИВОСТИ

§2.1. Постановка задачи

§ 2.2. Вариационная формулировка задачи о радиальной структуре

§ 2.3. Радиальная структура беспороговой диссипативной неустойчивости баллонных мод

ГЛАВА Ш. ОБОБЩЕННЫЙ КРИТЕРИИ УСТОЙЧИВОСТИ ИДЕАЛЬНЫХ БАЛЛОННЫХ МОД В ТОКАМАКЕ

§3.1. Постановка задачи

§ 3.2. Метод получения критерия и его анализ.

§3.3. Моделирование устойчивости плазмы при большом давлении.

§ 3.4. В ы в о д ы.

ГЛАВА 1У. ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ЭНЕРГОБАЛАНСА ДЛЯ ТОКАМАКА С КОМБИНИРОВАННЫМ АДИАБАТИЧЕСКИМ СЖАТИЕМ

§4.1. Физическая модель и исходные уравнения

§4.2. Эволюция плазменного шнура на начальной стадии разряда до сжатия.

§4.3. Модель дополнительного нагрева

§ 4.4. Результаты моделирования энергобаланса плазмы в режимах с комбинированным сжатием

§ 4.4.1. Влияние степени рециркуляции частиц на параметры плазмы в разряде

§ 4.4.2. Зависимость конечных параметров плазмы при сжатии от величины дополнительного нагрева.

§ 4.4.3. Определение режима с оптимальной средней плотностью плазмы перед сжатием

§ 4.4.4. Зависимость параметров разряда от величины начального тока перед сжатием

§ 4.5. Применение обобщенного критерия устойчивости для учета влияния баллонных мод на энергобаланс в токамаке с адиабатическим сжатием . П

Актуальность проблемы. Успешное развитие работ по магнитному удержанию плазмы в установках типа токамак /I/ в последнее время создало условия для перехода к новой стадии исследований по управляемому термоядерному синтезу. Эта стадия характеризуется проектированием и сооружением целого ряда больших установок, которые входят в строй в настоящее время и будут запущены в ближайшие годы. Основной задачей исследований на этих экспериментальных установках является изучение физических закономерностей поведения плазмы с параметрами, близкими к термоядерным.

Одной из важнейших и актуальных проблем, возникающих при исследовании плазмы с термоядерными параметрами, является проблема, связанная с ограничением на величину давления плазмы, которое могут накладывать баллонные моды в токамаке. Важность изучения вопроса о предельном достижимом значении параметра = 8эт<р>/&о » характеризуемого отношением давления плазмы к давлению тороидального магнитного поля, определяется тем, что от величины этого параметра существенно зависит возможность создания экономически выгодного термоядерного реактора.

Получить информацию о свойствах термоядерной плазмы можно и на установке не очень значительного масштаба, используя способ организации эксперимента на токамаке с комбинированным адиабатическим сжатием плазмы по малому и большому радиусу /2, 3/. Адиабатическое сжатие является эффективным способом дополнительного нагрева плазмы в токамаке. Однако оно может приводить к появлению в разрядах крутых профилей давления, локально неустойчивых относительно баллонных мод, а развитие неустойчивости баллонных мод сопровождается увеличением потерь энергии из плазмы и уменьшением ее энергетического времени жизни. Это обстоятельство необходимо учитывать при анализе энергобаланса плазмы в токамаке с адиабатическим сжатием. Поэтому актуальной задачей является теоретическое исследование устойчивости плазмы по отношению к баллонным модам, а также проведение численного моделирования энергобаланса для токамака с комбинированным адиабатическим сжатием с учетом эффектов, связанных с взаимным влиянием профилей параметров плазмы и зависимостью от них энергетического времени жизни, которое может существенно снижаться при получении в разрядах профилей давления, неустойчивых относительно баллонных мод, то есть не являющихся предельными.

Цель настоящей работы состоит в теоретическом исследовании устойчивости плазмы относительно баллонных мод и моделировании их влияния на процессы переноса в токамаке с комбинированным адиабатическим сжатием плазменного шнура по малому и большому радиусу.

Диссертация состоит из четырех глав. В первых трех главах диссертации представлены результаты теоретического анализа устойчивости плазмы в токамаке относительно баллонных мод желобковой неустойчивости, которые могут приводить к усилению выноса энергии из плазмы при некотором критическом значении давления в ней. Сам факт возникновения ограничения на величину давления был впервые обнаружен в численных расчетах /4, 5/. Однако интерпретация этих результатов была затруднена из-за несовершенства аналитической теории, имевшейся в то время.

Начало исследованиям баллонных мод было положено работой Фюрта, Киллена, Розенблюта и Коппи /6/. В этой работе переменная кривизна силовых линий магнитного поля токамака моделировалась в плоской геометрии введением эффективной силы тяжести. Было показано, что несмотря на выполнение критерия Сайдема /7/ плазма оказывается неустойчивой при давлении большем критического. Величина критического давления определялась соотношением длин участков силовых линий с "плохой" и "хорошей" кривизной. В то-камаке участки с "плохой" кривизной находятся на внешнем обводе тора, а участки с "хорошей" кривизной - на внутреннем обводе. Соотношение длин участков фиксировано. Поэтому возник вопрос о том, как получить критерий устойчивости, учитывающий баллонный эффект в случае реальной геометрии токамака.

Кадомцев и Погуце предложили метод расчета, использующий малость отношения полей в токамаке В3//В0 ( - магнитное поле тока, В0 - тороидальное магнитное поле), и получили следующий критерий устойчивости для токамака с круглым кожухом /8/: & -4 .о. (I)

А dr 2

Здесь S - шир, cj-rB0/RBj - коэффициент запаса устойчивости, ск --p'R^/Bo » R - большой радиус тора, U -магнитная яма в токамаке. В этом критерии первые два члена - сай-демовские, то есть они имеются и в цилиндрической геометрии. Видно, что шир стабилизирует рассматриваемые возмущения, а градиент давления - причина неустойчивости. В тороидальной геометрии появились два новых эффекта: дестабилизирующий баллонный, цропорци-ональный ( ^з=8зТ<р>/Вз , £ =г/R ), и средняя магнитная яма, играющая стабилизирующую роль.

Более аккуратные вычисления, проведенные Шафрановым и Юрчен-ко, показали, что происходит взаимное сокращение членов, связанных с тороидальной 1фивизной, и критерий устойчивости принимает вид /9/ i(H2) (2>

Независимо критерий (2) был получен из энергетического принципа Узйром и Хаазом /10/, 1фоме того Шафранов и Юрченко показали /9/, что этот критерий следует также и из общегеометрического критерия Мерсье /II/.

В критерии (2) баллонный эффект отсутствует, так как он целиком сократился с углублением магнитной ямы в токамаке за счет конечного давления плазмы, и поэтому при с^ > \ не существовало ограничения на допустимое давление. Шафранов и Юрченко сделали попытку получить ограничение на давление, учитывая более высокие члены разложения по параметру << \ /12/, но аккуратные вычисления Михайловского показали, что члены порядка являются стабилизирующими /13/. Таким образом, казалось, что баллонные моды не представляют опасности.

Ситуация резко изменилась после опубликования результатов расчетов /4, 5/, которые показали, что неустойчивость начинает развиваться при давлении плазмы выше критического значения. Необходимо было получить аналитический критерий устойчивости баллонных мод, чтобы понять физическую природу ограничения на давление плазмы.

Прогресс в аналитических исследованиях баллонных мод в этот период в значительной мере связан с методом расчета, развитым Коннером, Хасти и Тейлором /14/ и независимо Ли и Ван Дамом /15/, Этот.метод позволил отказаться от предположения о сильной локализации рассматриваемых возмущений, которое было существенным моментом всех предшествующих вычислений, и свел исследование системы уравнений в частных производных, к которой приводит задача об устойчивости плазмы, к решению одного обыкновенного дифференциального уравнения в фурье-пространстве. Вначале это уравнение, имеющее как периодические, так и непериодические коэффициенты, было решено численно /14/, а затем Погуце и Юрченко /16/, цред-ложив вариационный асимптотический способ для исследования этого уравнения, получили аналитический критерий устойчивости баллонных мод для круговых магнитных поверхностей:

С* / 1 \ Qi2 Я f- + -io(e >0. (3)

Этот критерий показывает, что при спадающем профиле тока, когда шир S > 0 , существует ограничение на давление плазмы, накладываемое неустойчивостью из-за баллонного эффекта, связанр ,2 / -vis1 ного с широм Ь<л . Дестабилизирующим является и член оСе , отражающий эффект перекрытия возмущений на близких рациональных магнитных поверхностях. Критерий (3) дает достаточно жесткое ограничение на величину давления плазмы. Однако этот критерий был получен для круговых магнитных поверхностей. В то же время хорошо известно, что за счет изменения формы магнитных поверхностей можно существенно улучшить устойчивость плазмы. В связи с этим возникает вопрос о получении аналитического критерия для некруговых магнитных поверхностей. Такой 1фитерий может существенно помочь в процессе оптимизации формы поверхностей в токамаке наряду с использованием численных кодов /17 - 19/.

Первая глава диссертации посвящена выводу и анализу аналитического критерия устойчивости плазмы относительно баллонных мод желобковой неустойчивости в токамаке с некруглым сечением /20/. В § I изложена процедура получения двумерного уравнения в частных производных, являющегося исходным при нахождении 1фитериев, дающих границу устойчивости плазмы относительно идеальных гидромагнитных неустойчивостей. В основу этой процедуры положен метод разложения исходных уравнений по малому параметру Bj / В0 > определяемому отношением магнитных полей в токамаке, предложенный впервые Кадомцевым и Погуце /8/. Полученное уравнение записано в симметричной форме, что позволяет легко доказать самосопряженность его дифференциального оператора. Поэтому можно написать для этого двумерного уравнения в частных производных функционал, который представляет собой потенциальную энергию возмущений. Показано, что этот функционал сводится к виду, полученному Погуце и Юрченко /21/ для потенциальной энергии прямой процедурой упрощения общего энергетического принципа методом разложения по малому параметру отношения полей. Во втором параграфе описывается, как, используя приближение эквивалентных гармоник /17/ при исследовании высоких мод баллонной неустойчивости ( П<^>> 1 , п -тороидальный номер гармоники возмущения, С^ - коэффициент запаса устойчивости), можно существенно упростить математическую задачу, сводя ее к анализу обыкновенного дифференциального уравнения в фурье-пространстве /14/, имеющего как периодические, так и непериодические коэффициенты, и соответствующего этому уравнению функционала. В § 3 приводятся результаты расчета критического давления плазмы, выполненного автором для установки токамак T-II /22/ с помощью аналитического критерия устойчивости относительно баллонных мод, полученного ранее Погуце и Юрченко /16/ для круговых магнитных поверхностей, и с помощью численного кода Захарова /17/. Отмечено соответствие результатов этих двух способов исследования устойчивости, показывающих, что баллонные моды же-лобковой неустойчивости дают достаточно жесткое ограничение на величину допустимого давления плазмы в случае токамака с круглым сечением,

В § 4 проводится вывод аналитического критерия устойчивости плазмы относительно баллонных мод желобковой неустойчивости в токамаке с некруглым сечением плазменного шнура, а в пятом параграфе этот критерий анализируется. Наибольший практический интерес, по-видимому, представляет случай совместного действия эллиптичности и треугольности ( D - образная форма сечения). Расчеты величины допустимого давления плазмы для токамака D «образной формы сечения, проведенные с помощью выведенного локального аналитического критерия устойчивости /20/, показывают, что хорошо подобранные профили тока и давления позволяют достигнуть относительно больших значений давления плазмы ( J5 ~ 4 6 % ), которые в два-три раза выше значений, получаемых для шнура 1фуглого сечения.

Во второй главе рассмотрен вопрос о радиальной структуре баллонных мод с высокими номерами /23 - 25/. Анализ радиальной структуры является вторым важным этапом при исследовании устойчивости плазмы после получения на первом этапе локальных критериев, о которых шла речь в первой главе. Для решения задачи о радиальной структуре идеальных баллонных мод Коннером, Хасти и Тейлором /26/ была разработана операторная теория возмущений, использующая малость отношения ширины огибающей гармоник возмущения к характерному размеру неоднородности. С использованием этой техники было показано, что локальный инкремент, получаемый на первом этапе исследования устойчивости, фактически является потенциалом радиального уравнения для огибающей гармоник, а наименьшее собственное значение для этого уравнения определяет истинный инкремент. В § 2 второй главы описывается предложенная вариационная формулировка задачи о радиальной структуре. Идея метода демонстрируется на уже известной задаче о радиальной структуре идеальных баллонных мод. В § 3 предложенный вариационный способ применяется для анализа радиальной структуры беспороговой дисси-пативной баллонной моды /27/, которая важна потому, что согласно локальной подходу она может существовать и при давлении плазмы ниже 1фитического значения, определяемого идеальными баллонными модами. Проведенное рассмотрение радиальной структуры этой неустойчивости показало, что, как и в случае идеальных баллонных мод, зависимость локального инкремента неустойчивости от радиуса определяет вид потенциала для радиального уравнения. Из конкретных расчетов следует, что локальный инкремент беспороговой диссипа-тивной баллонной моды нарастает к границе плазменного шнура, а потенциал радиального уравнения не имеет минимумов и монотонно спадает от центра плазмы к ее границе. Поэтому сделан вывод о том, что беспороговые дисеипативные баллонные моды не могут локализоваться в центральной области плазменного шнура, а "сваливаются" к границе плазмы. Для идеальных баллонных мод, возникающих при превышении порога по давлению, возможность локализации во внутренней области плазмы доказана.

В работах Мерсье /28/, Фильдинга и Хааза /29/, Захарова /17/ по идеальным баллонным модам было обнаружено, что при дальнейшем повышении давления выше критического плазма вновь становится устойчивой. Появление второй зоны устойчивости естественно связать с эффектом самостабилизации плазмы высокого давления, обнаруженным Михайловским и Шафрановым /30/. Для учета этого эффекта было необходимо получить аналитический критерий устойчивости баллонных мод с точностью до членов порядка cL^ /31/. Проведя последовательные расчеты, учитывающие члены такого порядка, Михайловский и Юрченко получили аналитический критерий устойчивости баллонных мод, связанных с широм, который описывает обе зоны устойчивости /32/. В этом же приближении был записан и критерий устойчивости типа Мерсье. Оба эти критерия являются необходимыми. Поэтому представлялось важным получить аналитический критерий устойчивости, который включал бы в себя ранее полученные критерии как частные случаи.

В третьей главе диссертации получен обобщенный аналитический критерий устойчивости плазмы в токамаке относительно идеальных баллонных мод желобковой неустойчивости /33, 34/. Этот критерий устойчивости является более жестким по сравнению с критерием устойчивости баллонных мод, связанных с широм /32/, и критерием типа Мерсье. В первых двух параграфах третьей главы дан вывод и анализ обобщенного 1фитерия устойчивости. В § 3 проведено модельное рассмотрение устойчивости плазмы при большом давлении, учитывающее эффект насыщения магнитной ямы при увеличении давления плазмы. Такое рассмотрение показало высокую чувствительность второй зоны устойчивости к изменению параметров равновесной конфигурации. Полученный в третьей главе обобщенный аналитический критерий устойчивости используется в следующей главе для учета влияния баллонных мод на энергобаланс в токамаке с комбинированным адиабатическим сжатием.

Четвертая глава посвящена численному моделированию энергобаланса для токамака с комбинированным адиабатическим сжатием плазмы на основе разработанной расчетной модели /35 - 38/. Эта работа была выполнена в связи с проводившимися в ИАЭ исследованиями по возможности применения адиабатического сжатия в токама-ках в качестве метода получения плазмы с термоядерными параметрами /2, 3/. Выводы работ /2, 3/ основаны на результатах расчетов, полученных с помощью "нульмерной" модели для описания энергобаланса плазмы, с адиабатическим пересчетом /39/ ее параметров при сжатии. Подход, основанный на "нульмерной" модели, обладает несомненным достоинством - простотой, дающей возможность проводить оценки вклада многих физических эффектов, а также учитывать инженерные проблемы. Поэтому такой подход полезен при первоначальном анализе энергобаланса в установках и часто испольуется /40, 41/. Тем не менее он требует уточнения и дополнения на основе использования более совершенных моделей для описания процессов переноса и в первую очередь одномерной модели.

Следует отметить, что впервые расчеты процессов переноса для токамака на основе одномерной модели были представлены в докладах Днестровского и Костомарова, а также Мерсье и Суббарамай-ера в 1969 году на конференции в Дубне /42/ по удержанию плазмы в замкнутых системах. В этих работах описаны результаты численного моделирования энергобаланса для разрядов в токамаке с омическим режимом работы. Предложенные модели энергобаланса оказались очень полезными для описания медленных эволюционных процессов в токамаке и поэтому получили широкое распространение и дальнейшее развитие /43 - 47/ с целью улучшения количественного согласия результатов расчета с экспериментом. Первые расчеты энергобаланса с включением эффектов поджатия по малому радиусу были выполнены в 1971 году Днестровским, Костомаровым и Павловой /48/. В этой работе было отмечено улучшение теплоизоляции плазмы при отрыве ее от стенок.

Использование сжатия по большому радиусу для нагрева плазмы было предложено Арцимовичем /49/. Такой способ нагрева плазмы исследовался на установке АТС в Принстоне с 1972 года /50/. В это же время проводились ранние эксперименты по поджатию плазмы по малому радиусу на установках типа "Туман" в ЛФТИ /51/.

Эксперименты по изучению возможности применения адиабатического сжатия в токамаках успешно ведутся и в настоящее время /52 - 55/. Параллельно с экспериментальными исследованиями проводятся и расчеты энергобаланса, причем в последнее время начинают использоваться транспортные модели, включающие расчет равновесия плазмы и эволюции магнитных поверхностей /56/. Модели для описания явлений переноса, лежащие в основе численных транспортных кодов, непрерывно совершенствуясь в соответствии с задачами, выдвигаемыми экспериментом, используются не только для анализа экспериментальных данных, но и применяются для прогнозирования результатов проектируемых установок. Примером тому могут являться работы по обоснованию международного проекта реактора-токамака ИНТОР /57, 58/.

Результаты, представленные > в четвертой главе диссертации, явились следствием численного моделирования энергобаланса, также выполненного с целью прогнозирования достижимых параметров плазмы в токамаке с комбинированным адиабатическим сжатием при различных режимах работы.

В § I четвертой главы описывается физическая модель, положенная в основу численных расчетов энергобаланса, приведены основные уравнения, получены граничные условия, учитывающие влияние геометрии системы при сжатии. Во втором параграфе излагается используемый способ описания эволюции плазменного шнура на начальной стадии разряда до сжатия, когда происходит пробой и подъем тока в плазме. Так как в токамаке с комбинированным сжатием для повышения эффективности сжатия предусмотрен дополнительный нагрев, в § 3 описана модель дополнительного нагрева инжекцией пучка нейтралов. Эта модель учитывает, что в системах относительно малых размеров и со слабым током ширина банана запертого горячего Иона оказывается порядка размеров системы, и его энергия существенным образом распределяется по радиусу, а не передается плазме лишь на той магнитной поверхности, где произошел акт ионизации нейтрала /59/. В § 4 представлены результаты расчетов динамики плазмы в режимах с комбинированным сжатием. Для построения общей картины энергобаланса прослеживается эволюция параметров плазмы на всех стадиях разряда. Рассмотрено влияние степени рециркуляции частиц на параметры плазмы в разряде. Показано, что во время сжатия по большому радиусу при полном рециклинге происходит ски-нирование плотности плазмы, обусловленное сгребанием диффузионных хвостов плазмы, так как плазменный шнур переходит в область с меньшим по размеру радиусом камеры. Анализируется зависимость конечных параметров плазмы при сжатии от величины мощности дополнительного нагрева, а также зависимость параметров разряда от величины начального тока перед сжатием. Отмечается существенное увеличение энергетического времени жизни при сжатии по малому радиусу. В результате расчетов найден режим с оптимальной средней плотностью плазмы перед сжатием, при которой можно получить наибольшую мощность термоядерного выхода. Проведенный анализ энергобаланса позволяет сделать вывод о возможности получения в токамаке с комбинированным адиабатическим сжатием плазмы с параметрами, близкими к термоядерным. При таких параметрах в плазме возможно возникновение баллонных мод желобковой неустойчивости, для которой энергетическим резервуаром является тепловая энергия плазмы. Приведенные в § 5 результаты численных расчетов энергобаланса, в которых моделировался эффект увеличения потерь энергии при возникновении неустойчивости баллонных мод путем резкого увеличения локального коэффициента электронной теплопроводности /60/ в момент нарушения обобщенного аналитического критерия устойчивости, дополняют представленный в предшествующих параграфах анализ энергобаланса для токамака с комбинированным адиабатическим сжатием.

Научная новизна. В диссертации рассмотрено влияние повышения давления плазмы на ее устойчивость. Получен аналитический критерий устойчивости плазмы относительно баллонных мод желобковой неустойчивости в токамаке с некруглым сечением.

Исследован вопрос о радиальной структуре беспороговой дис-сипативной баллонной моды. Показано, что эта неустойчивость стремится локализоваться на границе плазмы.

Выведен обобщенный локальный критерий устойчивости плазмы с закрепленной границей в токамаке относительно вдеальных баллонных мод желобковой неустойчивости. Полученный критерий уточнил границу области устойчивости и разъяснил вопрос о связи критерия Мерсье с критерием баллонных мод.

Разработана численная модель для анализа энергобаланса в токамаке с комбинированным адиабатическим сжатием, учитывающая влияние баллонных мод с помощью обобщенного аналитического критерия устойчивости. В результате численного моделирования обнаружено существование режима с оптимальной средней плотностью плазмы перед сжатием, позволяющей получить наибольший выход термоядерной мощности в момент максимального сжатия. Показано, что баллонные моды начинают приводить к снижению максимально достижимых параметров в токамаке с комбинированным сжатием при достаточно высоких значениях давления плазмы ( J5 ~2 % ).

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

Результаты, полученные в данной диссертации, были использованы при анализе возможности получения в токамаке с комбинированным адиабатическим сжатием плазмы с параметрами, близкими к термоядерным.

Автор выносит на защиту.

1. Аналитический критерий устойчивости плазмы относительно баллонных мод желобковой неустойчивости в токамаке с некруглым сечением. Анализ этого критерия показал, что в плазменном шнуре с D -образной формой сечения критические градиенты давления существенно больше, чем в шнуре круглого сечения, а допустимое давление плазмы может быть увеличено в два-три раза.

2. Вариационный подход для изучения радиальной структуры высоких мод баллонной неустойчивости, с помощью которого показано, что беспороговые диссипативные баллонные моды не локализуются в центральной области плазмы, а "сваливаются" к ее границе.

3. Обобщенный локальный критерий устойчивости плазмы с закрепленной границей в токамаке относительно идеальных баллонных мод, который является более жестким, чем критерий Мерсье и критерий баллонных мод, связанных с широм.

4. Результаты численного моделирования энергобаланса для токамака с комбинированным адиабатическим сжатием при учете влияния баллонных мод с помощью обобщенного критерия устойчивости.

Апробация работы. Основные результаты исследований, представленные в диссертации, докладывались на УШ конференции МАГАТЭ по физике плазмы и УТС (Брюссель, 1980), на X Европейской конференции по физике плазмы и УТС (Москва, 1981), на УП объединенном семинаре по вычислительной физике (Сухуми, 1981), на XI Европейской конференции по физике плазмы и УТС (Аахен, 1983) и опубликованы в работах /3, 20, 23-25, 33-38, 60/.

ЗАКЛЮЧЕНИЕ

Перечислим основные результаты диссертации:

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

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

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

4. Проведено моделирование энергобаланса для токамака с комбинированным адиабатическим сжатием на основе разработанной численной модели, учитывающей влияние баллонных мод с помощью обобщенного локального критерия устойчивости. Отмечено существенное увеличение полного энергетического времени жизни при сжатии по малому радиусу, что объясняется отжатием горячих областей плазмы от стенки и запиранием теплового потока по электронному каналу. Обнаружено появление скинирования профиля плотности при сжатии по большому радиусу. Этот эффект существенно зависит от степени рециркуляции частиц. Найден режим с оптимальной средней плотностью плазмы перед сжатием, позволяющей получить наибольший выход термоядерной мощности в момент максимального сжатия. Показано, что баллонные моды желобковой неустойчивости начинают приводить к снижению максимально достижимых параметров в токамаке с комбинированным адиабатическим сжатием при достаточно высоких значениях давления плазмы (J3 ~ 2%).

В заключение выражаю глубокую благодарность д.ф.-м.н. 0.П.Погуце за руководство работой и поддержку. Автор признателен д.ф.-м.н. В.А.Чуянову за постоянный интерес к работе по моделированию энергобаланса для токамака с комбинированным адиабатическим сжатием, а также коллективу сотрудников Отдела теории плазмы за обсуждение работы.

1. Кадомцев Б.В., Муховатов B.C., Шафранов В.Д. Магнитное удержание плазмы. - Шизика плазмы, 1983, т. 9, вып. 1. с. 5.

2. Демичев В.Ф., Чуянов В.А. Об адиабатическом сжатии плазмы с термоядерными параметрами в токамаках. Препринт ИАЭ-3222, М., 1979.

3. Todd A.M.M., Chance M.S., Greene J.M. et al. Stability limitations on high-beta tokamaks. Phys. Rev. Letters, 1977, vol. 38, II 15, p. 826.

4. Bateman G., Peng Y.-К.Ы. Hagnetohydrodynamic stability of flux-conserving tokamak equilibria. Phys. Rev. Letters, 1977, vol. 38, Ж 15, p. 829.

5. Purth H.P., Kill en J., Rosenbluth M.1T., Coppi B. Stabilization by shear and negative V". In: Plasma Physics and Controlled Nuclear Fusion Research, Culham, 1965, Vienna, IAEA, 1966, vol. 1, p. 103.

6. Сайдем Б.P. Устойчивость самосжатого линейного разряда.

7. В кн.: Физика горячей плазмы и термоядерные реакции. Труды второй международной конференции по мирному использованию атомной энергии, Женева, 1958, М., 1959, с. 89.

8. Кадомцев Б.Б., Погуце О.П. О желобковой неустойчивости плазмыв тороидальной геометрии. ДАН СССР, 1966, т. 170, № 4, с. 811.

9. Шафранов В.Д., Юрченко Э.И. Критерий желобковой неустойчивости плазмы в тороидальной геометрии. ЖЭТФ, 1967, т. 53, вып. 3, с. 1157.

10. Ю. Ware A.A., Haas P.A. Stability of a circular toroidal plasma under average magnetic v/ell conditions. Physics of Fluids, 1966, vol. 9, Ж 5, P. 956.

11. Mercier C. Un critere necessaire de stabilite hydromagnetique pour un plasma en symetrie de revolution. Huclear Fusion, 1960, vol. 1, I 1, p. 49.

12. Шафранов В.Д., Юрченко Э.И. Удержание плазмы в токамаке с J^»! при произвольном распределении тока. In : Plasma Physics and Controlled Nuclear Fusion Research, Vienna, IAEA, 1971, vol. 2, p. 519.

13. Mikhailovskio A.B. Tokamak stability at high plasma pressures. Uuclear Pusion, 1974, vol. 14, N 4, p. 483.

14. Connor J.W., Hastie R.J., Taylor J.B. Shear, periodicity, and plasma ballooning modes. Phys. Rev. Letters, 1978, vol. 40, N 6, p. 396.

15. Lee Y.C., Van Dam J.W. Kinetic theory of ballooning instabilities. In: Proceedings of the finite Beta Theory Workshop, Varenna, 1977, p. 93.

16. Погуце О.П., Юрченко Э.И. Дестабилизирующее влияние шира и предельное давление плазмы в токамаке. Письма в ЖЭТФ, 1978,т. 28, вып. 6, с. 344.

17. Zakharov L.E. Hi gh-vvave-number MHD-mode stability in high-pressure tokamaks. In: Plasma Physics and Controlled

18. Nuclear Fusion Research, Innsbruck, 1978, Vienna, IAEA, 1979, vol. 1, p. 689.

19. Вабщев П.Н., Дегтярев Л.М., Медведев С.Ю. Предельные значения J5 относительно неустойчивости баллонных мод в токамаке.- Физика плазмы, 1983, т. 9, вып. I, с. 163.

20. Bernard L.C., Helton F.J., Moore R.V/., Todd T.N. MHD beta limits: scaling laws and comparision with doublet III data. Nuclear Fusion, 1983, vol. 23, IT 11, p. 1475.

21. Погуце О.П., Чудин H.B., Юрченко Э.И. Устойчивость баллонных мод в токамаке с некруглым сечением. Физика плазмы, 1980, т. 6, вып. 3, с. 621.

22. Pogutse О.Р., Yurchenko E.I. Energy principle and kink instability in toroidal plasma with strong magnetic field.- Nuclear Fusion, 1978, vol. 18, N 12, p. 1629.

23. Mukhovatov V.S. Beam-heatihg research at the I.V.Kurchatov Institute. In: 2nd Joint Varenna Grenoble International Symposium on Heating in Tokamak Plasma, Como, 1982, vol. 2, p. 761.

24. Pogutse O.P., Chudin N.V. Radial structure of resistivethballooning modes. In: 11 European Conference on Controlled Fusion and Plasma Physics, Aachen, 1983, vol. 7D, B-49, P. 163.

25. Погуце О.П., Чудин H.B. Радиальная структура диссипативной неустойчивости баллонных мод. Препринт ИАЭ-3828/6, М., 1983.

26. Погуце О.П., Чудин Н.В. Радиальная структура диссипативных баллонных мод. Физика плазмы, 1984, т. 10, вып. 2, с. 339.

27. Connor J.W., Hastie R.J., Taylor J.В. High mode numberstability of an axisymmetric toroidal plasma. Proceedings of the Royal Society, A., 1979, vol. 365, p. 1.

28. Погуце О.П., Юрченко Э.И. Беспороговые диссипативные баллонные моды. Письма в ЖЭТФ, 1980, т. 31, вып. 8, с. 479.

29. Mercier С. MHD stability criteria for localized displacements. In: Plasma Physics and Controlled Nuclear Fusion Research, Innsbruck, 1978, Vienna, IAEA, 1979, vol. 1, p. 701.

30. Sykes A., Turner M.F., Fielding P.J., Haas P.A.

31. MHD theory of ballooning modes applied to JET and the effect of anisotropic pressure on ballooning modes. Ibidem, 1979, vol. 1, p. 625.

32. Михайловский А.Б., Шафранов В.Д. Эффект самостабилизации плазмы высокого давления в тороидальных ловушках, КЭТФ, 1974, т. 66, вып. I, с. 190.

33. Mikhailovskij А.В., Demchenko V.V. The shear effect onstability of ballooning modes in closed magnetic traps, th- In: 10 European Conference on Controlled Fusion and Plasma Physics, Moscow, 1981, vol. 1, B-11.

34. Михайловский А.Б., Юрченко Э.И. Аналитическая теория идеальной индуцируемой широм баллонной моды токамака. Препринт ИАЭ-3505/6, М., 1981.

35. Погуце О.П., Чудин Н.В., Юрченко Э.И. Обобщенный критерий устойчивости идеальных баллонных мод в токамаке. Физика плазмы, 1983, т. 9, вып. I, с. 156.

36. Погуце О.П., Чудин Н.В., Юрченко Э.И. Обобщенный критерий устойчивости идеальных баллонных мод в токамаке при конечном давлении плазмы. Препринт ИАЭ-3607/6, М., 1982.

37. Вихрев В.В., Погуце О.П., Чудин Н.В., Юрченко Э.И. Процессыпереноса в токамаке с адиабатическим сжатием. Препринт ИАЭ-3321/8, М., 1980.

38. Pogutse О.P., Chudin N.V., Chuynov V.A. Transportcalculation in a tokamak with adiabatic compression. In: th

39. European Conference on Controlled Fusion and Plasma Physics, Moscow, 1981, vol. 1, B-17.

40. Gribkov V.M., Morozov D.X., Pogutse O.P., Chudin If.V. numerical model for plasma heating by neutral injection in a high field tokamak with adiabatic compression. In: Ibidem, 1981, vol. 1, H-21.

41. Коvan I.A., Chudin H.V., Shurygin R.V. Numerical simulation of ion-cyclotron heating in D2+He++ mixure in tokamak. In: Ibidem, 1981, vol. 1, H-9.

42. Furth H.P., Yoshikawa S. Adiabatic compression of tokamak discharges. The Physics of Fluids, 1970, vol. 13, H 10, p. 2593.

43. Абрамов B.A., Кукушкин А.С., Пистунович В.И. Нульмерная модель расчета плазменных параметров гибридного термоядерного реактора-токамака. Препринт ИАЭ-3064, М., 1978.

44. Васильев Н.Н., Лукаш В.Э., Недоспасов А.В. и др. Термоядерная электростанция с высокотемпературным реактором типа токамак. в сб.: Доклады Всесоюзной конференции по инженерным проблемам термоядерных реакторов. - Л., НИИЭФА, 1977, т. I, с. 114.

45. Днестровский Ю.Н., Костомаров Д.П. Международная конференция по замкнутым системам. Дубна, 1969; I.iercier е., Soubbaramayer, там же (цит. по /44/).

46. Mercier С., Soubbaramayer. numerical models for plasmathevolution in tokamak devices. In: 5 European Conferenceon Controlled Fusion and Plasma Physics, 1972, Grenoble, vol. 2, p. 157.

47. Днестровский Ю.Н., Костомаров Д.П. Модели баланса энергии и частиц в установках токамак. В сб.: Вычислительные методы в физике плазмы. Дополнение. - М., Мир, 1974, с. 483.

48. Duchs D.F., Post D.E., Rutherford Р.Н. A computer model of radial transport in tokamaks. ITuclear Fusion, 1977, vol. 17, IT 3, P. 565.

49. Хоган Дж.Т. Многокомпонентные модели переноса в токамаке. В сб.: Управляемый термоядерный синтез. - М., Мир, 1980, с. 142.

50. Днестровский Ю.Н., Костомаров Д.П. В кн.: Математическое моделирование плазмы. -М.: Наука, 1982, с. 222.

51. Днестровский Ю.Н., Костомаров Д.П., Павлова Н.Л. Поджатие плазмы нарастающим продольным магнитным полем в токамаке. -Письма в ЖЭТФ, 197I, т. 13, вып. 12, с. 697.

52. Artsimovich L.A. Tokamak devices. ITuclear Fusion, 1972, vol. 12, IT 2, p. 215.

53. Bol K., Ellis R.A., Eubank H. et al. Adiabatic compression of the tokamak discharge. Phys. Rev. Letters, 1972, vol. 29, IT 22, p. 1495.

54. Каганский М.Г. Адиабатическое сжатие плазмы в токамаке. Л., Наука, 1979, с. 170.

55. Berezovsij E.L., Vil'dzhynas M.I., Gladushchak V.I. et al. Experiment on ohmic heating and compression of a plasmain the Tuman-2A tokamak. In: Plasma Physics and Controlled ITuclear Fusion Research, 1978, Vienna, IAEA, 1979, vol. 1, p. 349.

56. Rui-yan Lin, Kei-ichi Hirano et al. Simulation of adiabaticcompression plasma. Japan Institute of Plasma Physics, ifegoya University Report: IPPJ-494, 1980.

57. Cima G., Gowers C.W., Guire K.M.IvIc et al. Compression experiments on the Tosca tokamak. UKAEA, Abingdon, England, Culham Laboratory Report: CLM-R213, 1981.

58. Golant V.E., Kalmykov S.G. Experiments on adiabaticxi cLcompression of a tokamak plasma. In: 2 Joint Varenna Grenoble International Symposium on Heating in Tokamak Plasma, Como, 1982, vol. 2, p. 975.

59. Hogan !J.T. The accessibility of high-beta tokamak states. Nuclear Fusion, 1979, vol. 19, IT 6, p. 753.

60. International Tokamak Reactor. Zero Phase. Vienna, IAEA, 1980.

61. Днестровский Ю.Н., Кукушкин А.С., Лысенко C.E. и др. Баланс энергии и частиц в установке ИНТОР. В сб.: ИНТОР. Международный токамак-реактор (материалы СССР). Отв. редактор-составитель В.А.Абрамов. -М.: ИАЭ, 1980, ч. I, с. 29.

62. Днестровский Ю.Н., Костомаров Д.П., Лысенко С.Е. Моделирование нагрева плазмы в установке T-II при инжекции пучка нейтралов. Препринт ИАЭ-2978, М., 1978.

63. Погуце О.П., Чудин Н.В., Юрченко Э.И. Гидромагнитные неустойчивости и критическое J2> в токамаке. Препринт ИАЭ-3239/6, М., 1980.

64. Погуце О.П., Юрченко Э.И. Баллонные эффекты и устойчивость плазмы в токамаке. В кн.: Вопросы теории плазмы. - М.: Энергоатомиздат, 1982, вып. II, с. 56.

65. Соловьев Л.С., Шафранов В.Д. Замкнутые магнитные конфигурации для удержания плазмы. В кн.: Вопросы теории плазмы. - М., Атомиздат, 1967, вып. 5, с. 3.

66. Кадомцев Б.Б. Гидромагнитная устойчивость плазмы. В кн.: Вопросы теории плазмы. - М.: Госатомиздат, 1963, вып. 2, с. 132.

67. Shafranov v.d., Yurchenko E.I. Influence of ballooning effects on plasma stability in closed systems. Nuclear Fusion, 1968, vol. 8, N 4, p. 329.

68. Погуце О.П., Юрченко Э.И. Влияние шира на устойчивость плазмы в токамаке. Физика плазмы, 1979, т. 5, вып. 4, с. 786.

69. Todd A.M.LI., Manickam J., Okabayashi M. et al. Dependence of ideal MHD kink and ballooning modes on plasma shape and profiles in tokamaks. Nuclear Fusion, 1979, vol. 19, П 6, P. 743.

70. Юрченко Э.И. О равновесии симметричных плазменных конфигураций. ЖТФ, 1967, т. 37, № 8, с. 1458.

71. Михайловский А.Б. К теории баллонных мод в токамаке с большим

72. Jb . Физика плазмы, 1978, т. 4, вып. 6, с. 1226.

73. Брагинский С.И. Явления переноса в плазме. В кн.: Вопросы теории плазмы. -М.: Госатомиздат, 1963, вып. I, с. 183.

74. Hinton F.L., Hazeltine R.D. Theory of plasma transport in toroidal confinement systems. Reviews of Modem Physics, 1976, vol. 48, N 2, p. 239.

75. Кукушкин А.С. Программа dkl76 для расчета плазмофизических параметров реактора-токамака в одномерной транспортной модели. Препринт ИАЭ-3140, М., 1979.

76. Кадомцев Б.Б., Погуце О.П. Турбулентные процессы в тороидальных системах. В кн.: Вопросы теории плазмы. - М., Атомиздат,1967, вып. 5, с. 209.

77. Галеев А.А., Сагдеев Р.З. Явления переноса в разреженной плазме в тороидальных магнитных ловушках. ЖЭТФ, 1967, т.53, вып. I, с. 348.

78. Коврижных Л.М. Процессы переноса в тороидальных магнитных ловушках. ЖЭТФ, 1969, т. 56, вып. 3, с. 877.

79. Шафранов В.Д. О классической теплопроводности в тороидальном плазменном шнуре. Атомная энергия, 1965, т. 19, вып. 2,с. 120.

80. Duchs D.F., Furth Н.Р., Rutherford Р.Н. Radial transport in tokamak discharges. In: Plasma Physics and Controlled Nuclear Fusion Research, 1971, Madison, Vienna, IAEA, 1971, vol. 1, p. 369.

81. Post D.E., Goldstein R.J., Grimm R.C. et al. Computational studies of impurity effects, impurity control and neutral beam injection in large tokamaks, In: Plasma Physics and Controlled Nuclear Fusion Research, 1978, Vienna, IAEA, 1979, vol. 1, p. 471.

82. Preiffer W., Waltz R.E. Empirical scaling laws for energy confinement in ohmically-heated tokamaks. Uuclear Fusion, 1979, vol. 19, N 1, p. 51.

83. Waltz R.E., Guest G.E. Maximum energy confinement time in joule-heated tokamaks. Phys. Rev. Letters, 1979, vol. 42, H 10, p. 651.

84. Абрамов B.A., Вихрев B.B., Погуце О.П. Начальная стадия развития разряда в токамаке. Физика плазмы, 1977, т. 3, вып. 3, с. 512.

85. Козлов Б.Н. Скорости термоядерных реакций. Атомная энергия, 1962, т. 12, $ 3, с. 238.

86. Пистунович В.И. Токамак с инжекцией быстрых атомов. Препринт ИАЭ-2209, М., 1972.

87. Коган В.И., Мигдал А.Б. Зависимость спектра тормозного излучения от электронной температуры плазмы. В кн.: Физика плазмы и проблемы управляемых термоядерных реакций. - М., Изд-во АН СССР, 1958, т. I, с. 174.

88. Трубников Б.А. Универсальный коэффициент выхода циклотронного излучения из плазменных конфигураций. В кн.: Вопросы теории плазмы. -М., Атомиздат, 1973, вып. 7, с. 274.

89. Абрамов В.А., Погуце О.П., Юрченко Э.И. Пробой газа в токамаке. Физика плазмы, 1975, т. I, вып. 4, с. 536.

90. Днестровский Ю.Н., Костомаров Д.П., Переверзев Г.В., Тарасян К.Н. Ионизация и подъем тока в токамаке.- Препринт ИАЭ-2980, М., 1978.

91. Днестровский Ю.Н., Костомаров Д.П., Переверзев Г.В., Тарасян К.Н. Перезамыкание магнитных поверхностей и перемешивание тока в начальной стадии разряда в токамаке. Физика плазмы, 1978, т. 4, вып. 5, с. 1001.

92. Dnestrovskij Yu.N., Kostomarov D.P., Karetkina U.V. et al. Equilibrium and energy balance in tokamaks. In: VI European Conference on Controlled Fusion and Plasma Physics, Moscow, 1973, vol. 1, p. 41.

93. Riviere A.C. Penetration of fast hydrogen atoms into a fusion reactor plasma. ITuclear Fusion, 1971, vol. 11, N 4, p. 363.

94. Гришанов H.B., Днестровский Ю.Н., Кареткина H.B., Костомаров Д.П. Передача энергии быстрых ионов плазме. Физикаплазмы, IS76, т. 2, вып. 2, с. 260.

95. Stodick W., Goldston R., Sauthoff П. et al. Transport studies in the Princeton Large Torus. In: Plasma Physics and Controlled Nuclear Fusion Research, Brussels, 1980, Vienna, IAEA, 1981, vol. 1, p. 9.

96. Stix Т.Н. Fast wave heating of a two-component plasma. Nuclear Fusion, 1975, vol. 15, N 5, p. 737.

97. Азизов Э.А., Алексеев Ю.А., Бревнов Н.Н., Велихов Е.П. и др. Основные физические и инженерные проблемы создания токамака с сильным магнитным полем и адиабатическим сжатием плазмы. -Атомная энергия, 1982, т. 52, вып. 2, с. 108.

98. Connor J.W., Taylor J.В., Turner M.F. Ideal MHD ballooning instability and scaling law for confinement. Nuclear Fusion, 1984, vol. 24, N 5, p. 642.