В.Н.Докука

advertisement
На правах рукописи
ДОКУКА Владимир Николаевич
ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ
ФИЗИЧЕСКИХ ПРОЦЕССОВ В ПЛАЗМЕ
ТОКАМАКОВ TCV, KTM, JUST-T
Специальность: 01.04.08 – физика плазмы
АВТОРЕФЕРАТ
Диссертации на соискание ученой степени
кандидата физико-математических наук
Троицк - 2008
Работа выполнена в ГНЦ РФ Троицкий институт инновационных и термоядерных исследований (ТРИНИТИ), г. Троицк, Московской обл.
Научный руководитель:
кандидат физико-математических наук
Р.Р. Хайрутдинов
Официальные оппоненты:
доктор физико-математических наук, профессор
В.А. Курнаев
кандидат физико-математических наук
Ю.Ю. Пошехонов
Ведущая организация:
ФГУП НИИ электрофизической аппаратуры
им. Д.В. Ефремова (г. Санкт-Петербург)
Защита состоится « 05 »
декабря 2008 г. в 15-30
часов
на заседании диссертационного совета ДС 201.004.01 в ГНЦ РФ ТРИНИТИ
по адресу: 142190, Московская обл., г. Троицк, ул. Пушковых, владение 12
С диссертацией можно ознакомиться в библиотеки ГНЦ РФ ТРИНИТИ
Автореферат разослан « 31 » октября 2008 г.
Ученый секретарь
диссертационного совета,
кандидат физико-математических наук
2
С.А. Казаков
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность работы. Наиболее значимые результаты в решении проблемы
управляемого термоядерного синтеза были достигнуты на установках типа
токамак. В настоящее время ведутся работы по реализации во Франции (Кадараш) проекта Интернационального термоядерного экспериментального реактора-токамака (ИТЭР) для демонстрации научной и технологической осуществимости и использования реакции синтеза D-T для мирных целей.
Достижение проектных параметров плазмы в установках токамак невозможно без создания моделей и комплексов программ, направленных как на
моделирование процессов в плазме, так и управление этими процессами.
Неотъемлемой частью исследований на установках токамак становится выполнение предварительных численных экспериментов – изучение плазменных
процессов методами математического моделирования – предполагаемого
сценария работы установки до проведения реального физического эксперимента. Имеющиеся расчетно-теоретические модели плазмы и база экспериментальных данных позволяют проводить «численные эксперименты» еще на
стадии проектирования токамаков следующего поколения.
В этой связи является актуальным разработка и верификация на действующих токамаках плазмофизических кодов для проведения численного моделирования управляемых разрядов, что требует объединения численных моделей плазмы, диагностических, исполнительных (рис. 1,а) и управляющих
устройств в единый комплекс «компьютерный токамак» (рис. 1,б). При этом
важно исследование различных режимов токамака как в разомкнутой (рис.
1,а), так и в замкнутой системе управления (рис. 1,б). Выполненные исследования по численному моделированию экспериментальных разрядов на токамаке TCV (Швейцария, г. Лозанна) и сценариев разрядов установок КТМ (Казахстан, г. Курчатов), JUST-T (Россия) [1-21] с учетом систем управления
плазмой обозначили в диссертационной работе новый подход к численному
эксперименту – созданию «компьютерного токамака» - Numerical Tokamak
Project.
Токамак TCV (Tokamak Configuration Variable) проектировался специально для исследования функциональных преимуществ вытянутой плазмы в широком диапазоне плазменных конфигураций. Имеющаяся система магнитного
управления TCV обеспечивает большое разнообразие возможных равновесных конфигураций плазмы. На данном токамаке верификация разработанных
в диссертации численных моделей плазмы по экспериментальным данным
оказывается эффективной.
Создание токамаков-реакторов невозможно без разработки и испытаний
новых материалов для защиты первой стенки и приемных диверторных пластин. Предполагается, что сооружаемый в настоящее время токамак КТМ
(Казахстанский материаловедческий токамак) станет базовой установкой для
проведения системных исследований образцов материалов первой стенки и
3
дивертора при воздействии потоков энергии от 0,1 до 20 МВт/м 2 в широком
диапазоне временных экспозиций.
Поддержание устойчивой диверМодель
Входы
Модельплазмы
плазмы Выходы
торной конфигурации плазмы и
токамака
токамака
обеспечение мощных корпускулярных потоков на диверторные
а
пластины в процессе работы
установки приводит к необхоЧисленная модель Выходы
Входы
димости разработки соответЧисленная модель
плазмы
плазмытокамака
токамака
ствующих сценариев разряда в
КТМ. Требуется также проведеКомпьютерный токамак
ние численного моделирования
Управляющее
работы системы магнитного
Управляющее
устройство
устройство
управления током, положением
и формой плазмы, способной
обеспечить такие режимы, при
б
использовании разработанных
Рис. 1. Модель плазмы в токамаке а) в разомкнумоделей плазмы. Проверка и
той и б) замкнутой системе управления - комплекс
обоснование работоспособности
«компьютерный токамак»
систем магнитного управления
плазмой является задачами, обуславливающих актуальность темы диссертации.
Переход от демонстрационного токамака-реактора ИТЭР к термоядерной
энергетической установке возможен через этап создания на основе токамака
объемного источника нейтронов (ОИН) для решения задач ядерной энергетики. Российская концепции ОИН на базе сферического токамака JUST-T (Joint
Upgraded Spherical Tokamak for Transmutation), требует обоснования возможности комбинированного индуктивно - неиндуктивного подъема тока и достижения 100% неиндуктивного поддержания стационарного тока плазмы,
являющегося основным рабочим режимом будущего термоядерного реактора.
Указанные выше задачи предопределили тематику диссертации, в основу
которой положены работы по численному моделированию управляемых процессов в плазме токамаков, проведенных в период 2000-2007 гг. лично автором диссертации, либо при его непосредственном участии [1-21]. В области
численного моделирования управляемых плазменных процессов в токамаках
сделан большой вклад как отечественными специалистами (Л. Захаров, Г.
Переверзев, Р. Хайрутдинов, В. Лукаш, А. Кавин и др.), так и зарубежными
учеными (A. Portone, R. Albanese, D. Humphreys, M. Walker, J. Lister, L. Pearlstein, L. LoDestro и др.). В диссертации проведен краткий обзор их работ, это
позволило выбрать направление диссертационной работы, которое определило логику проведенных исследований, их новизну и практическую значимость.
4
Целью работы является разработка комплекса «компьютерного токамака»,
настройка и проведение его верификации в условиях разрядов в плазме токамака TCV, создание методики разработки сценариев разрядов в токамаке
КТМ с использованием комплекса «компьютерного токамака», а также обоснование концепции компактного сферического токамака JUST-T в качестве
объемного источника нейтронов для трансмутации минорных актинидов отработавшего ядерного топлива.
Методы исследования. Для анализа работы комплекса «компьютерного токамака» использовались методы современной технологии математического
моделирования совместно с методом последовательных приближений в соответствии с выбором набора критериев качества в настройке и верификации
моделей по накопленным данным физического эксперимента.
Научная новизна работы. Новизна диссертационной работы определяется
характером разработкой математических моделей, программ и методов построения комплекса «компьютерного токамака».
1. Разработаны и впервые включены в программно-вычислительный комплекс ДИНА модули для расчета коэффициентов переноса на основе бомовской/гиробомовской, МММ (Multi-Mode Model) 5.10 и МММ95 транспортных моделей.
2. Разработан новый комплекс «компьютерный токамак» и проведено моделирование экспериментальных разрядов токамака TCV с реально действующей системой магнитного управления плазмой.
3. Разработана оригинальная методика построения сценариев разряда в установках токамак. Предложены базовые сценарии разряда в проекте токамака КТМ с использованием разработанной методики.
4. Впервые разработана плазмофизическая часть концепции ОИН на базе
компактного токамака JUST-T.
Практическая значимость:
1. Разработанный на основе кода ДИНА комплекс «компьютерного токамака» используется для моделирования экспериментальных разрядов установки TCV, применяется при отработке систем магнитного управления
плазмой и при планировании будущих экспериментов в установках токамак. Комплекс используется в учебном процессе при подготовке специалистов по управлению плазмой в токамаке.
2. Созданные сценарии разряда в токамаке КТМ внесли определяющий
вклад в ряд проектных решений и были использованы при проектировании и сооружении установки. Результаты вошли в раздел «Физическое
обоснование токамака КТМ» проектной документации КТМ. Комплексное
исследование влияния моделей переноса на параметры плазмы и длительность разряда в токамаке КТМ необходимо при планировании будущих
экспериментов.
3. Разработанная плазмофизическая часть концепции применения токамака
JUST-T в качестве объемного источника нейтронов используется для тех5
нико-экономического обоснования сооружения такой установки. Данные
результаты, направленные на повышение экологической чистоты атомных
электростанций, имеют особую значимость для ядерной энергетики. Методика, разработанная для физического обоснования токамака как объемного источника нейтронов, может быть применена также при разработке
проектов термоядерных реакторов.
Личный вклад автора. Основным личным вкладом автора являются полученные научные результаты по численному моделированию на программном
комплексе ДИНА физических процессов в плазме токамаков TCV, КТМ и
JUST-T. Результаты получены для моделей плазмы в разомкнутом и замкнутом контурах систем магнитного управления плазмой «компьютерных токамаков», которые созданы автором в международном стандарте среды математического моделирования MATLAB/SIMULINK. Исследована динамика
плазмы для различных сценариев разрядов установки TCV в обратной связи
системы управления, в которой использовалась модель управляющего
устройства, разработанная J.B. Lister (EPFL-CRPP, Швейцария). Это позволило верифицировать модели плазмы по экспериментальным данным установки
TCV с высокой степенью достоверности. Автором разработана оригинальная
методика разработки сценариев разряда в установках токамак и по ней созданы базовые сценарии разрядов для проекта токамака КТМ. Для базовых сценариев разрядов в КТМ исследована и обоснована работоспособность системы магнитного управления (А.А. Кавин - НИИЭФА им. Д.В. Ефремова) в
замкнутом контуре с моделью плазмы. Методами численного моделирования
автором получены основные физические и технические характеристики плазмы с целью физического обоснования концепции сферического токамака
JUST-T в качестве объемного источника нейтронов.
Положения, выносимые на защиту:
1. Создание и включение в код ДИНА программных модулей для расчета
коэффициентов переноса на основе бомовской/гиробомовской, МММ 5.10
и МММ95 транспортных моделей.
2. Разработка S-блоков кода ДИНА в системе MATLAB-Simulink для решения задач моделирования плазмы «компьютерного токамака» в замкнутом
контуре управления. Результаты моделирования разрядов плазмы на установке TCV и их сравнение с экспериментальными данными, в том числе и
полученными по восстановительному коду LIUQE.
3. Создание методики планирования сценариев разряда в токамаке и ее применение для разработки базовых сценариев разрядов в КТМ. Проведение
тестирования системы магнитного управления плазмой КТМ с ее полномасштабной моделью.
4. Разработка эффективных сценариев работы токамака JUST-T, обеспечивающих необходимые потоки 14 МэВ нейтронов на внутреннюю стенку
бланкета с минорными актинидами.
6
Апробация работы. Основные результаты диссертации докладывались и
обсуждались на научных семинарах в Институте ядерного синтеза РНЦ
«Курчатовский институт», ТРИНИТИ, на Международных научных конференциях, симпозиумах и совещаниях:

Международные симпозиумы по технологии термоядерного синтеза
(SOFT-22, г. Хельсинки, Финляндия, 2002, SOFT-23, г. Венеция, Италия,
2004);

Международные конференции европейского физического общества по
физике плазмы и управляемому термоядерному синтезу (EPS-29, г. Монтрё, Швейцария, 2002; EPS-30, г. Санкт-Петербург, Россия, 2003; EPS-31,
г. Лондон, Англия, 2004; EPS-32, г. Тарагона, Испания, 2005; EPS-33, г.
Рим, Италия, 2006; EPS-34, г. Варшава, Польша, 2007);

Международные конференции МАГАТЭ по термоядерной энергии
(FEC-20, г. Виламура, Португалия, 2004; FEC-21, г. Ченду, Китай, 2006).
Публикации. По результатам диссертационной работы опубликовано 21 печатная работа, из которых 11 представлены в ведущих отечественных и зарубежных журналах: Вопросы атомной науки и техники - серия “Термоядерный синтез” [6, 14, 17], серия "Электрофизическая аппаратура" [13]; Plasma
Devices and Operations [12, 18, 20]; Plasma Physics and Controlled Fusion [5];
Fusion Engineering аnd Design [2, 8, 9]; остальные работы опубликованы в
трудах международных конференций [1, 3, 4, 7, 10, 11, 15, 16, 19, 21].
Структура и объем диссертации. Диссертация состоит из введения, четырех
глав, заключения и списка цитируемой литературы из 156 наименований. Работа содержит 153 страниц, включает 80 рисунков и 10 таблиц.
СОДЕРЖАНИЕ ДИССЕРТАЦИИ
Во введении приведена краткая история исследований по тематике диссертационной работы, обосновывается актуальность темы диссертации, формулируются цель работы, дана краткая аннотация диссертации по главам. Отражены новизна примененных научно-технических решений, также практическая
ценность работы и положения, вынесенные на защиту.
В первой главе представлено описание физической модели плазмы токамака, численно реализованной в коде ДИНА. Код позволяет вести расчеты
динамики плазмы токамака со свободной границей во внешних магнитных
полях в масштабах резистивного времени пассивной структуры токамака:
110-3-110-1 сек. Представлено описание полоидальной магнитной системы
токамака, включающей активные полоидальные катушки и структуру пассивной стабилизации. В основе физической модели плазмы лежит двумерное
равновесие плазмы со свободной границей во внешних магнитных полях,
7
одномерный (усредненный по магнитным поверхностям) перенос частиц,
тепла и полоидального магнитного потока.
Геометрия магнитных поверхностей определяется из решения двумерного
аксиально-симметричного уравнения Грэда-Шафранова:
20 r j , (r , z )  S p ,
 2   2  1  
L
 2 

r r 20  r I  (r  r ) ( z  z ),
z 2
r
i i
i
i
i 1

(r , z )  S ,
 dp 0 dF 2  .

j  2r 

 d 4r d 
Уравнение связывает распределения давления плазмы р и полоидального
электрического тока F с распределением полоидального магнитного потока 
таким образом, чтобы выполнялось условие для баланса сил в плазме, занимающей область Sp. Здесь L, (ri, zi), Ii – соответственно, число внешних стационарных проводников, их положение и сила тока в них. Токи подчиняются
 dI

 RI  U , где
уравнению Кирхгоффа для магнитных контуров:
–
I dt
I
квадратная матрица Якоби, а вектор токов I является вектором состояния
динамической модели объекта.
Давление плазмы p, входящее в правую часть уравнения равновесия,
определяется из уравнения баланса частиц и внутренней энергии. Для описания диффузии энергии, частиц и магнитного поля в программном комплексе
ДИНА используется система одномерных уравнений, полученная усреднением по тороидальным магнитным поверхностям.
Для многих установок экспериментально подтверждается неоклассический механизм ионной теплопроводности по теории Галеева-Сагдеева. Экспериментальные исследования показывают, что потоки энергии электронов и
потоки частиц в токамаках являются аномальными, поэтому пользоваться для
их описания неоклассическими выражениями нельзя. Теоретическое описание процессов переноса в плазме, во многом вызванное работой над проектом
ИТЭР, совершенствуется, однако прогнозирующая способность моделей еще
явно недостаточна. Для глобального энергетического времени жизни E плазмы накоплена обширная база данных. Поэтому в транспортных моделях часто
используются эмпирическими скейлингами – зависимостями  E от основных
параметров плазмы (Ip, a, R, k, ne, Bt, Paux и т.д.), полученных при анализе
экспериментальных данных. При моделировании нашли широкое применение
различные полуэмпирические бомовские/гиробомовские модели переноса.
Теоретические методы расчетов переноса базируются на анализе динамики плазмы с учетом возникновения в ней тех или иных неустойчивостей. К
числу теоретически обоснованных моделей можно отнести Multi-Mode Models (МММ 5.10 и МММ95), широко используемых в таких кодах, как
BALDUR, ASTRA, JETTO и CRONOS. В этих моделях коэффициенты пере8
носа рассчитываются из условия неустойчивости различных дрейфовых мод:
ионной градиентной моды (ITG – Ion Temperature Gradient), электронной градиентной моды (ETG – Electron Temperature Gradient), моды на запертых частицах (TEM – Trapped Electron Mode) и т.д.
В MMM 5.10 модели перенос энергии рассчитывается с учетом неустойчивости на запертых электронах (TEM) и i - моды, инициированной ионной
температурно-градиентной неустойчивостью (ITG), которые доминируют в
центре плазмы, а также из-за резистивных баллонных мод (RBM – Resistive
Ballooning Mode), которые определяют перенос на краю плазмы. Вклад от
запертой электронной моды рассчитывается по модели Dominguez и Waltz, от
ионной температурно-градиентной неустойчивости ( i -моды), соответственно, по модели Hamaguchi и Horton, и резистивные баллонные моды учитываются по модели Carreras и Diamond. Резистивная баллонная мода обеспечивает почти Bohm-диффузию, тогда как другие соответствуют gyroBohmдиффузии.
Дальнейшее развитие моделей типа МММ привело к созданию версии
МММ95 на основе модели Weiland-Nordman, самосогласованно учитывающей вклад ITG и TEM мод. Тороидальная версия ITG моды основывается на
комбинации ионного градиентного и центробежного дрейфов с одной стороны и E  B дрейфа с другой. В квазилинейном приближении, пренебрегая
возможным эффектом зацепления фурье-компонентов возмущения плотности
из-за неоднородностей плазмы, получено дисперсионное уравнение – полином четвертой степени, что объясняет возможность существования до двух
неустойчивых мод. Квазилинейные потоки частиц и тепла, возникающие изза турбулентности, усредняются по периоду колебаний и по пространственному изменению флуктуации. Вклад от всех неустойчивостей должен быть
просуммирован. В квазилинейном приближении полагается, что перенос
определяется наиболее быстрорастущей модой. Полученные коэффициенты
диффузии содержат значительную температурную зависимость для ионов и
электронов соответственно:
 R
R 

,

L
L
Ti
Ti
,th 

 R
R

 LTe LTe,th
i  TeTi1 2 

i  TeTi1 2 


,


где R- большой радиус, LT – обратный относительный градиент температуры.
При этом локальные пороги для ITG и TE-моды задаются как:
R
4 20



,

 (1  f ) 


LTi
3 n
9
t
(1  f t ) n
2(1  f t )
2 n2 (1  f t )
2
2 Ln ,
T ,
f ,
R
4
20  
1 ,


 1     e  n  R   t
Ti
1  ft
LTi 3 n 9 2   n 
где ft – доля запертых электронов и Ln- обратный относительный градиент
плотности. Такая зависимость коэффициентов диффузии от температуры ведет к увеличению переносов с ростом температур. Отметим, что для TEM и
ITG мод в МММ 5.10 также имеется зависимость  Te3 2 . Кроме модели
9
Weiland, МММ95 включает в себя модель Gusdar-Drake для учета вклада в
перенос дрейфово-резистивной (DR – Drift Resistive) и кинетической баллонных (KB – Kinetic Ballooning Mode) мод.
Система уравнений переноса дополняется граничными и начальными
условиями. Для уравнений равновесия на границе ставится краевое условие Iрода, которое определяется из уравнения диффузии магнитного поля. Токи в
обмотках управления и пассивной стабилизации, а также в проводящих элементах вакуумной камеры рассчитываются с учетом приложенных внешних и
индуцированных процессами в плазме напряжений.
Код ДИНА реализует численное решение уравнения равновесия во внешних магнитных полях совместно с транспортом кинетических параметров
внутри плазмы и уравнением Кирхгоффа для системы магнитных контуров.
Уравнение Грэда-Шафранова решается методом Бунемана на прямоугольной
сетке для определения граница плазмы {r, z}b и методом “обращения” переменных в полярных координатах {,} для нахождения координат магнитных
поверхностей для расчета метрических коэффициентов. Уравнения переноса
энергии для электронов и ионов решаются матричным вариантом метода потоковой прогонки, разработанного специально для сильно связанных между
собой уравнений. Нелинейности реализуются с помощью итерационных циклов.
Вторая глава посвящена тестированию Simulink-блоков, разработанных
на основе кода ДИНА для моделирования управления плазмой в токамаке
TCV, спроектированного для исследования функциональных преимуществ
вытянутой по вертикали плазмы в широком диапазоне плазменных конфигураций. Полоидальное сечение TCV показано на Рис. 2. Гибкость в управлении
формой плазмы обеспечивается независимым питанием токов в обмотках
управления формой и обмотках индуктора. Система магнитной диагностики
включает 38 магнитных зондов и 38 датчиков измерения полоидального магнитного потока. Установка TCV снабжена системой электронного циклотронного нагрева (ECH – Electron Cyclotron Heating) и возбуждения токов
увлечения (ECCD - Electron Cyclotron Current Drive). Гиротроны обеспечивают полную мощность ECH нагрева  4,2 МВт с длительностью импульса  2
сек. Для разработки систем управления широко применяется пакет прикладных программ Simulink-MATLAB, ориентированный на моделирование
сложных систем и основанный на модульном принципе организации моделирования. Система используется на TCV при подготовке и анализе сценариев
разрядов, а также при разработке и тестировании систем управления с использованием упрощенных моделей плазмы RZIP, CREATE-L и т.п. Модульная структура кода ДИНА позволяет выделить из него модуль управляющего
устройства системы магнитного управления (Рис. 1,б), а оставшуюся часть
кода, снабдив необходимыми входами-выходами, организовать в виде Sблока Simulink, который рассматривается как физическая модель плазмы в
10
токамаке. Объединением S-блока кода ДИНА с блоком управляющего
устройства в замкнутый контур создается вычислительный комплекс, названный «компьютерный токамак» и предназначенный для проведения расчетного сопровождения экспериментов.
Необходимые начальные данные определяются с помощью восстановительного кода LIUQE. При моделировании
разряда с ECH нагревом плазмы и генерацией токов увлечения ECCD на соответствующий вход S-блока ДИНА подаются текущие параметры каждого
пучка EC: центр (R, Z) и ширина зоны
поглощения (R, Z), поглощаемая
мощность PECCH, ток увлечения IECCD,
рассчитанные лучевым кодом TORAY
по данным восстановительного кода
LIUQE. Верификация комплекса ДИНА
достигается настройкой «внутренних»
параметров модели переноса с использованием глобального скейлинга RLW
(Rebut-Lallia-Watkins) и двухпараметрической модели внутреннего транспортного барьера. На динамический вход
Рис. 2. Вакуумная камера TCV, полоидальные магнитные катушки, датчики
модели плазмы поступают также средмагнитного поля внутри вакуумной каменяя плотность плазмы вместе с настраиры (обозначены знаком “”), датчики
ваемым параметром пикированности
полоидального магнитного потока (обоnb/n0 (отношением плотностей на гранизначены знаком “”)
це и в центре плазмы). Настройка параметров ведется методом последовательных приближений, ставя целью достижения разумного согласия результатов моделирования и эксперимента.
Новые подходы расчета переноса связаны с концепцией «открытой архитектуры» (open architecture version) путем включения в комплекс «компьютерный
токамак» дополнительного S-блока для моделирования переноса. Развитие
указанной концепции привело к созданию комплекса ДИНА-CRONOS, объединившего код ДИНА с транспортным кодом CRONOS (V. Basiuk et al).
Моделирование разряда № 12610 с чисто омическим нагревом плазмы
включает в себя стадии подъема плазменного тока, квазистационарной фазы
и вывода тока из плазмы. Для предотвращения расхождения между моделируемыми и экспериментальными первичными токами индуктора из-за различий проводимости плазмы в численном моделировании и эксперименте введен контур искусственной обратной связи, корректирующий возможное расхождение первичных токов обмоток OH1 и OH2 через регулирование проводимости плазмы. Продемонстрировано приемлемое согласие результатов мо11
делирования и эксперимента. Моделирование показывает несколько большую величину вытянутости плазмы по вертикали в сравнении с экспериментом, что ведет к ухудшению управлению по вертикали.
В разряде № 19692 с нецентральным ЭЦР нагревом плазмы управление
формой плазмы отсутствует, а плазменный шнур управляется только системой обратной связи по положению. Результаты моделирования вместе с данными LIUQE представлены на Рис. 3.
При нецентральном нагреве формируется более широкий профиль
плазменного тока и снижается величина внутренней индуктивности
шнура, что ведет к увеличению вытянутости плазмы из-за почти постоянного квадрупольного компонента полоидального вакуумного
поля. На стадии ЭЦР нагрева плазмы искусственная обратная связь,
корректирующая расхождение первичных токов в индукторе, не исРис. 3. Результаты моделирования разряда №
19692: Ip - ток плазмы, (Rmag, Zmag) - положение
пользуется. Динамика вытянутости
магнитной оси, 95 – вертикальная вытянушнура после включения мощности
тость. Данные ДИНА (-) и LIUQE (--)
ЭЦР – нагрева лучше согласуется с
данными эксперимента при улучшенном времени удержания. Изучение влияния нецентрального ЭЦР нагрева на вытянутость плазмы продолжились в
разряде № 22832 при отсутствии управления формой плазмы и в разряде №
24884, соответственно, с управлением вытянутостью плазмы. Приемлемое
согласие с экспериментом для разряда № 22832 получается путем соответствующей настройки пикированности плотности плазмы и фактора улучшения времени удержания. Предполагается формирование внутреннего транспортного барьера на радиусе  ITB  0.6 . Выявлена чувствительность вытянутости к профилю плотности. Вытянутость оказывается слишком высокой при
«пикированном» профиле плотности и низкой для «плоского» профиля.
Управление вытянутостью плазмы ведется в разряде № 24884, и тенденция
увеличения вытянутости 95 из-за уменьшения индуктивности li компенсируется уменьшением квадрупольного поля. Результаты эксперимента и моделирования показали, что на стадии ЭЦР нагрева изменения 95 оказываются минимальными.
Проведено моделирование разряда № 22895 с ЭЦР нагревом и высокой
долей бутстреп-тока fbs ~ 60 %. Динамика тока плазмы согласуется с экспериментальными данными, а величина бутстреп тока совпадает результатами
кода PRETOR. Однако при устойчивом уровне мощности ЭЦР нагрева отмечены колебания токов увлечения до  50 %, что отражает трудности расчета
токов увлечения по восстановленным профилям. В разряде № 24828 индук12
тивное поддержание тока плазмы принудительно снижается до нуля путем
управления изменением тока в индукторе через обратную связь, реализуя 100
% неиндукционное поддержание тока, после чего один из гиротронов отключается и через  50 мсек включается другой.
На TCV выполнен цикл экспериментов по изучению характеристик мод,
локализованных на краю плазмы (ELM - Edge Localized Modes). Краевые колебания обычно сопровождают режим улучшенного удержания плазмы (Hмоду) и рассматриваются как характерный признак режимов с улучшенным
удержанием. Возбуждение ELM связано с возмущениями плотности тока на
краю плазмы. Энергия, выделяющаяся за одно колебание ELM, может создать
недопустимую тепловую нагрузку на компоненты установки. Существование
линейной зависимости энергии ELM от интервала между ними telm позволяет предположить, что их разрушительное воздействие можно ослабить, искусственно повысив частоту их возбуждения. На TCV с этой целью использовались электромагнитные возмущения для модулирования плотности тока на
периферии плазмы и изменения интервала между колебаниями. Для возбуждения использовалась пара обмоток G контура стабилизации плазмы по вертикали (Рис. 2). Последовательность прямоугольных импульсов напряжения
продолжительностью 1 мс с варьируемой задержкой между импульсами
направлялась на выход контура обратной связи вертикальной стабилизации.
Возмущение плотности тока на краю плазмы в ответ на вынужденные вертикальные смещения плазмы можно оценить следующим образом. Внешний
полоидальный поток определяется как  ext      G , где  является полоидальным потоком от токов в обмотках E, F и OH, скорость изменения которого обеспечивает постоянное напряжение по обходу шнура - VL, и поток
G, определяемый токами G-обмоток и, наведенными ими токами в вакуумной камере. Допуская, что плазма движется как твердое тело с вертикальной
скоростью uz t  и что     G , уровень возмущения по отношению к
омическому компоненту плотности тока на краю можно записать как:
  
 J V
1   G

 
 uz
.
J  VL
VL  t
z 
Моделирование разрядов велось с учетом системы управления, использованной в эксперименте. Возмущения следовали каждые 10 мс на протяжении
( 0,2 сек) всей стадии возбуждения колебаний. Для получения типичного
результата от накладываемых возмущений в присутствии внутренних срывов
и медленного действие управления формой, циклы возмущений были когерентно усреднены с целью увеличения отношения сигнал – шум. Моделировались разряды № 20333 и № 20334 с SNL (Single Null Lower) диверторными
конфигурациями в верхней части вакуумной камеры (z = + 0,23 м) и противоположными знаками возмущения напряжений и подобный им разряд № 22678
13
с перевернутым равновесием (Х-точка сверху) - SNU (Single Null Upper) диверторной конфигурацией в нижней части камеры (z = - 0,23 м).
Результаты
моделирования
представлены на Рис. 4. У всех
трех разрядов наиболее вероятное возбуждение ELM происходило при движении плазмы в том направлении, в котором увеличивалась плотность
тока на краю. На поверхности
плазмы индуцируется отрицательное напряжение и, следовательно, отрицательный ток,
при движении плазмы в
направлении обмотки, создающей X-точку. В разряде №
20333, максимум вероятности
наблюдается за фронтом возмущения напряжения. В разрядах № 20334 и № 22768
максимум вероятности имеет
место на передней фронте
сигнала возмущения. ПриемРис. 4. Результаты моделирования трех разрядов в
лемое согласие результатов
TVC с возмущением тока в G-обмотке. Динамика
моделирования и эксперименнапряжения VG, тока IG в G-обмотке и вертикального
та сделало возможным исположения z на протяжении одного импульса возпользовать комплекс для исмущения. Возмущение тока на краю Jedge нормализуется относительно среднего тока на краю. Конфигуследования механизмов возрации равновесия для разрядов (№ 20333, № 20334, №
буждения ELM в других уста22768) представлены слева направо
новках.
Результаты Главы 2 показывают, что при соответствующей настройке параметров транспортной модели удается достичь приемлемого согласия результатов моделирования и экспериментов, и дает основание к использованию S-блока кода ДИНА при разработке управляемых сценариев проектируемых установок.
Третья глава посвящена разработке сценариев разрядов в плазме токамака КТМ, который создается с целью исследования материалов, предназначенных для токамаков следующего поколения. В токамаке КТМ предполагается достичь удельных тепловых потоков на пластины дивертора, сравнимых
с потоками в ИТЭР. Минимизация времени экспозиции для получения запланированных интегральных потоков энергии возможна лишь при увеличении
длительности разрядов в условиях ограниченного запаса полоидального потока центрального соленоида (ЦС). Полоидальная система установки включа14
ет односекционный ЦС, шесть полоидальных магнитных катушек управления
и пару «быстрых» обмоток управления по вертикали. Все обмотки имеют
независимые источники питания.
Методика создания сценариев разряда в КТМ использует итерационную
процедуру нахождения программных значений (уставок) для управляемых
параметров путем проведения моделирования на каждом интервале
t0 , tk , k  1 : n с последующей коррекцией уставок при t  tk . В сценарии с чисто омическим нагревом плазмы планируется провести отладку работы систем управления токамака, систем поддержания равновесия плазменного
шнура обратными связями.
Базовый индуктивный сценарий с током плазмы Ip=0,75 МА разрабатывался с учетом стадии пробоя и формирования плазменного шнура, используя
в качестве начальных данных результаты, полученные при моделировании
начальной стадии ввода тока плазмы кодом TRANSMAK. На Рис. 5, 6 показаны уставки, расчетные данные управляемых параметров и соответствующих
управляющих токов в процессе моделирования разряда.
При моделировании использовалась модель переноса со скейлингом Т-11 для
электронов и неоклассикой для ионов (T11&nc). Разработанный сценарий позволяет обеспечить поддержание тока и
формы плазмы, отрабатывая программные токи в обмотках управления и
управляя током плазмы и положением
магнитной оси с обратной связью.
Управляющие токи, полученные при
Рис. 5. Программные (*) и расчетные (-)
координаты магнитной оси Rmag, Zmag и
моделировании разряда, близки к сооттоки управления положением магнитной
ветствующим программным значениям.
оси IPF3 и IHFC
Конфигурация плазмы на момент окончания стационарной стадии развития разряда представлена на Рис. 7. Резистивные потери полоидального потока связаны с электронной температурой
плазмы, поэтому было предпринято изучение влияние транспортных моделей на параметры индуктивного сценария разряда.
Моделирование проводилось для транспортных модулей с использованием неоалкаторного (NA) скейлинга для электронов и неоклассики для ионов
(NA&nc), NA-скейлинга по электронам и
Рис. 6. Программные (*) и расчетные (-)
ионам (NA&NA), полуэмпирической
токи плазмы Ip и индуктора ICS
бомовской/гиробомовской
модели
(Bohm&gyro-Bohm), транспортных моделей МММ 5.10 и МММ95. Длительность разрядов в индуктивном сценарии для транспортных моделей (T11&nc), (NA&nc) и (NA&NA) составила  1 сек, примерно такой она оказа15
лась в сценарии с моделью переноса (МММ95). Наименьшую продолжительность разряда  0,83 сек показали расчеты с моделями (Bohm&gyroBohm) и
(МММ 5.10). Оценка расхода полоидального потока    ind   res на
стадии подъема тока, определяется индуктивными  ind  L p I p и резистивными потерями  res  CEjima0 R0 I p .
Внешняя и внутренняя индуктивности
шнура определяются согласно формулам:
  8R0  li
 ,
L   R ln
 2
p
0
0

 
 a   2


li  ln 1.65  0.89q95  1 ,
где
q95 
Рис. 7. Конфигурация плазмы в момент
t=1сек. На левом верхнем графике
представлен профиль плотности
плазменного тока jp по радиусу шнура
5a 2 Bt 1   2 1  2 2  1.2 3  1.17  0.A65  .
2
R0 I p
2
1  A12


При значении коэффициента CEjima= 0,4
индуктивные и резистивные потери составляют  0,83 Вб и  0,34 Вб соответственно. Результаты моделирования дают ind=0,74 Вб res=0,39 Вб, при
этом обмотки управления формой обеспечивают расход 0,28 Вб. На стационарную стадию остается не более 0,58 Вб запаса полоидального потока ЦС, и
ее длительность  0,6 сек. Расход полоидального потока определяется как
1.3
резистивными потерями, так и долей бутстреп-тока f bs  Ccs   p . Значения


входящих в формулу Cbs и полоидальной беты p определяются согласно
формулам:
2
2
и  p  t Bt B p 2 ,
Cbs  1.32  0.235q95
 0.0185q95
где
t 
neTe  niTi и
I p MA .
Bp 
2
2
Bt 20
5a 12
В омическом разряде числитель в формуле определения t можно заменить
выражением Pohm E V p , где Vp  2 2a 2 R0 . Для сценария с моделью переноса
(T-11&nc) доля бутстреп-тока составляет  0,09. Основные параметры плазмы
для разных транспортных моделей в сценариях с омическим нагревом плазмы
представлены в Табл. 1.
Таблица 1. Параметры плазмы при индуктивном поддержании тока
Te, eV
E, ms
Ures, V
fbs
POhm,MW
T-11&nc
340
67
1,06
0,10
0,79
NA&nc
384
87
0,92
0,12
0,69
NA&NA
328
64
1,09
0,10
0,82
16
B&gyroB
278
46
1,32
0,08
0,98
МММ5.10
278
48
1,32
0,08
0,98
МММ95
310
58
1,17
0,09
0,87
Начальная фаза сценария базового разряда с ВЧ - нагревом плазмы в диапазоне частот ионно-циклотронного резонанса (ИЦР-нагрев), вплоть до
включения ВЧ – мощности, аналогична соответствующей стадии сценария с
индуктивным поддержанием тока. Моделирование базового сценария разряда
с ВЧ - нагревом плазмы в диапазоне частот ионно-циклотронного резонанса
проводилось для модели переноса с использованием скейлинга IPB98(2,y)
при различных значениях фактора улучшения удержания HH. Расчеты выполнялись для параболического профиля плотности плазмы с пьедесталом, равным 0,9. Предполагается, что мощность ИЦР - излучения PICRH , равная 5,0
МВт, поглощается в равных долях ионами и электронами вблизи магнитной
оси на ширине  а/2. В базовом сценарии с ИЦР-нагревом при HH=1,3 длительность стационарной стадии разряда превышает 4 сек. Допустимость использования скейлинга IPB98(2,y) в установках с аспектом A=2 требует экспериментального подтверждения. Изучение влияния моделей переноса на
параметры плазмы в сценарии с ИЦР- нагревом направлено на решение этой
проблемы. Перенос рассчитывался с использованием модели переноса со
скейлингом Т-11 для электронов и неоклассики для ионов, полуэмпирической
бомовской/гиробомовской модели, транспортных моделей МММ 5.10 и
МММ95. Основные параметры плазмы, полученные для разных транспортных моделей в сценарии с ИЦР – нагревом плазмы, представлены в Табл. 2.
Таблица 2. Параметры плазмы для разряда с ИЦР нагревом
Te, eV
E, ms
Ures, V
fBS
 tsc, s
E/E,IPB98(y,2)
IPB(98,y)&HH=1.3
1321
34
0,15
0,33
4,63
1,33
T-11&nc
1173
33
0,18
0,32
4,47
1,29
B&gyroB
878
31
0,26
0,38
3,31
1,19
МММ5.10
876
27
0,29
0,31
2,90
1,05
МММ95
1137
31
0,20
0,33
3,10
1,21
Для разряда c ИЦР-нагревом плазмы числитель в формуле определения t
можно заменить выражением PICRH E V p . Оценка доли бутстреп-тока в разряде с моделью переноса (IPB(98,y)&HH=1.3) дает  0,33, что близко к расчетной. В сценариях разряда с ИЦР нагревом для транспортных моделей
(IPB98(2,y)&HH=1,3) и (T-11&nc) длительность разряда  4,5 сек, а моделей
(Bohm&gyroBohm), (МММ 5.10) и (МММ95) длительность  3 сек. У всех
использованных моделей время удержания выше его значения по скейлингу
IPB98(2,y).
Разработанные сценарии разряда использовались при тестировании системы магнитного управления плазмой в токамаке КТМ, которая в контуре с
обратной связью должна обеспечить выполнение программного сценария
разряда. Как видно на Рис. 8, структура системы управления током, положением и формой плазмы КТМ, разработанной в НИИЭФА (А.А. Кавин), строится по двухконтурной схеме, принятой для ИТЭР. Отдельно выделен «быстрый» контур управления для стабилизации скорости вертикального смещения
17
плазмы относительно нулевого значения и контур «медленного» управления
током и формой плазмы. Признано целесообразным использование разных
регуляторов для управления током и формой плазмы для лимитерной и диверторных конфигураций. Моделирование выполнялось в системе MATLABSimulink с использованием комплекса «компьютерного токамака».
g
Моделирование системы магнитZ
ного управления в омическом
МОДЕЛЬ
I
ТОКАМАКА
сценарии разряда начинается по
I
КTМ
I
окончанию формирования лимитерной плазменной конфигурации
на внешнем обводе вакуумной
Модель
Регулятор
камеры и током Ip ~ 50 кА. В соисточника
по вертикали
U
питания
Z
ответствии с программным сценарием на момент t  0,18 сек форДиагнос.
Модель
Регулятор по
U
задержка
мируется диверторная конфигураисточника
форме и току
питания
плазмы
ция с плазменным током 430 кА.
Проектного значения Ip = 750 кА
Рис. 8. Схема управления плазмой КТМ. Вектор
ток плазмы достигает на момент t
g содержит параметры формы плазмы
~ 0,3 сек. Моделирование омического сценария проводилось при использовании разных инженерных моделей
переноса энергии. Стабилизация вертикального положения плазмы обеспечивалась обмоткой HFC, источником питания которой служит инвертор напряжения, работающий в автоколебательном режиме с частотой до 1 кГц и амплитудой напряжения 360 В. Управление током и формой плазмы поддерживается токами остальных обмоток. На Рис. 9 приведены ошибки (отклонения
от программы) управляемых параметров g5 ,  , Z axis , I p на лимитерной стадии.
p
p
pf
HFC

HFC
0.04
0.04
0.02
0.02
err(K)
err(g5), m
pf17
0
-0.02
-0.04
0
0.05
0.1
t, s
0.15
-0.04
0.2
0
0.05
0.1
t, s
0.15
0.2
0.03
err(Ip), MA
0.04
err(Zax), m
0
-0.02
0.02
После формирования диверторной
конфигурации управление передается
«диверторному» регулятору, на вход
которого подаются ошибки управляемых параметров диверторной стадии g15 , I p , I CS , I PF1 PF 6 . Поведение оши-
бок пяти зазоров на диверторной фазе
показано на Рис. 10. Точность отработки регулятором программных задающих воздействий по зазорам составляет
Рис. 9. Ошибки управляемых параметров
примерно 1 см, причем в стационарном
err(g5), err(K), err(Zax) и err(Ip) на лимитерной стадии подъема тока
режиме точность достигает нескольких
миллиметров. Программа поддержания тока плазмы на этой фазе выполняется с точностью нескольких кА. Моделирование системы магнитного управления было проведено также в разряде с дополнительным нагревом плазмы. По
достижению стационарного значения тока плазмы включалась мощность ВЧ 0.02
0
-0.02
0
0.05
0.1
t, s
0.15
0.2
0.01
0
-0.01
0
0.05
0.1
t, s
0.15
0.2
18
нагрева. За времена порядка нескольких E температура плазмы возрастает и
наблюдается значительный рост бутстреп-тока до  200 кА. На стационарной
стадии разряда с ВЧ - нагревом регулятор отрабатывает программные значения зазоров g1-5 с точностью  1 см.
Программное значение тока плазмы на
этой фазе поддерживается с точностью
до нескольких кА. Хотя по сравнению
с омическим разрядом, параметры
плазмы в стационарной стадии разряда
значительно отличаются, один и тот же
«диверторный» регулятор обеспечивал
управление с высокой точностью в
обоих случаях. Это свидетельствует,
что используемые регуляторы обеспеРис. 10. Ошибки пяти зазоров err(g1-5) на
чивают достаточные запасы устойчидиверторной стадии в омическом режиме
вости замкнутой системы управления
при численном моделировании. Реализация системы в физическом эксперименте может потребовать расширения запасов устойчивости, т.к. в модели
токамака невозможно предусмотреть всех неопределенностей реального объекта.
Результаты Главы 3 показали, что разработанные сценарии разрядов в
КТМ обеспечивают развитие разрядов в соответствии с поставленными задачами, учитывая электромагнитные параметры полоидальных обмоток. Созданные сценарии разряда служили основой при разработке системы магнитного управления плазмой токамака КТМ и моделировании ее работы. Регуляторы обеспечили высокую точность стабилизации параметров плазмы, и токи
в обмотках не превышали наложенных ограничений.
0.04
1
2
3
4
5
0.03
err(g1-5), m
0.02
0.01
0
-0.01
-0.02
-0.03
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
t, s
Четвертая глава посвящена исследованию стационарных объемных источников нейтронов на базе токамака. Создание токамака JUST-T (R=2 м, A =
2 и  = 1,7), работающего в стационарном режиме и обеспечивающего плотность потока 14 МэВ нейтронов на внутреннюю стенку бланкета n0,3
МВт/м2, позволит применить его для трансмутации, т.е. превращения минорных актинидов из отработанного топлива тепловых реакторов в короткоживущие изотопы путем обработки потоком нейтронов. Излагается сценарий
двух - ступенчатого индукционно-неиндукционного выхода на стационарную
стадию развития разряда в токамаке. Показано, что применение соленоида стартёра, предполагающего использование запаса полоидального потока без
перемагничивания ЦС, обеспечит подъем тока плазмы до 40% от минимального уровня Ip,min, необходимого для удержания энергии -частиц в плазме.
На окончании индукционной стадии подъема тока плазмы формируется диверторная конфигурацию с Х-точкой в нижней части вакуумной камеры.
19
Дальнейшее наращивание тока плазмы и его квазистационарное поддержание предлагается осуществлять неиндуктивно бутстреп-током и токами
увлечения, генерируемых тангенциальной инжекцией нейтрального пучка
(NBI – Neutral Beam Injection) атомов дейтерия. Приводится математическое
описание физических процессов, сопровождающих инжекцию пучка быстрых
нейтралов в плазму, и численных моделей для расчета профиля поглощения
пучка нейтральных атомов в плазме, вычисления функции распределения
быстрых ионов, расчета токов увлечения. Исследована зависимость эффективности генерации токов увлечения от энергии горячих нейтралов и температуры фоновой плазмы. Эффективность генерации тока NB пучком
нейтральных атомов дейтерия Ab=2 c энергией Eb и прицельным параметром
b  Rtan g R0 в плазме с температурой Te и плотностью ne и эффективным зарядом ионов Zeff можно оценить по формулам:
 NBI  Abd
где
и
R
T
25 1016
1  f s  tan g e J  x, y  F  Zb , Z eff ,   a 2R0  ,
R0
R0 ne
F ( Z b , Z eff ,  ) 
1
1
1  G  Z eff ,    ,


Z b Z eff 
G  Z eff ,    1.55  0.85 Z eff
Множитель
G  Z eff ,  
тока. Величина

  0.20 1.55 Zeff   .
учитывает влияние запертых электронов на величину
x   Eb E c 
12
  Eb 10 AbTe 
12
относительная скорость нейтралов
пучка, и y  4Z eff 5 Ab - параметр рассеяния. Коэффициент Abd= 0,11А/Вт.
Величина fs определяет долю нейтрального пучка, прошедшего сквозь плазму
без поглощения. Для выражения J  x, y  справедлива приближенная формула: J  x, y   x2 4  3 y  x 2  x  1.39  0.61y 0.7  . При характерных значениях парамет

ров плазмы ne=11020м-3, Te=6 кэВ, Zeff=1,5 и энергии пучка Eb=140 кэВ имеем
эффективность генерации NBI  0,087 А/Вт. При фиксированном значении
параметра y эффективность генерации тока в D-T плазме максимальна при
некотором значении xmax(y)= 2,21, что соответствует энергии атомов пучка
2
 y AbTe , равной примерно 600 кэВ. Максимум функции J x, y  , в
Eopt  10 xmax
которой присутствует зависимость от энергии пучка, довольно пологий. Значения этой функция лишь на 10% меньше максимального, начиная с энергии
пучка 260 кэВ. При мощности нейтральной инжекции 45 МВт обеспечивается
неиндукционный сценарий подъема и поддержания тока плазмы при заданном уровне нейтронной нагрузки на внутреннюю стенку бланкета n  0,3
МВт/м2. Показано, что при токе плазмы Ip =Ip, = 4,5 MA, необходимом для
удержания -частиц в плазме, и факторе улучшения удержания HIPB98(y,2) = 1,6
доля бутстреп-тока fbs  0,5. Эффективность генерации тока NB ~ 0,05 А/Вт.
Достигается неиндукционный режим поддержания тока fBS + fNB=1.
20
Изучался вопрос увеличения нейтронного выхода за счет D-T реакций на
пучке. Сечение ядерных взаимодействий быстрых дейтонов с ядрами трития
может быть значительно больше сечения DT-реакций максвелловской плазмы. Поскольку вероятность реакций на пучке пропорциональна плотности
ядер трития, было решено увеличить концентрацию трития до 70%, чтобы
максимально использовать вклад данного эффекта в суммарный нейтронный
выход. Показано, что эффективность генерации нейтронов на пучке Pnbeam/PNB
составляет  0,5.
Исследовалось влияние профиля плотности плазмы на параметры в стационарной стадии разряда. Расчеты показывают (Рис. 11), что для профиля
плотности без пьедестала обеспечивается плотность потока  0,42 МВт/м2,
тогда как при плоском профиле он снижается до величины  0,29 МВ/м2.
По мере уплощения профиля
плотности плазмы все представленные параметры: бутсреп-ток
Ibs, ток увлечения INB и ток плазмы
Ip снижаются на (20-25) %.
Влияние внутренних транспортных барьеров учитывалось в
рамках двухпараметрической модели путем задания: а) положения
барьера через нормализованный
радиус ITB, и б) увеличения коэффициента
теплопроводности
Рис. 11. Зависимость тока плазмы Ip, тока увле
в
области
сильной
турбулентITB
чения INB, бутстреп-тока Ibs и уровня удельной
нейтронной нагрузки n от величины пьедестала
ности. Исследования выполнены
профиля плотности плазмы
при заданном положении барьера
ITB и изменяющимся факторе ITB, а также при заданном ITB и изменяющимся положении барьера ITB. Энергия пучков нейтралов полагалась равной
140/300 кэВ, с мощностями – 25/20 МВт. Профиль плотности - параболический без пьедестала n(0):nb=1:0,1.
Влияние увеличения коэффициента теплопроводности ITB изучалось при
заданном положении барьера ITB=0,6. На Рис. 12 видно возрастание плотности бутстреп-тока jbs слева от барьера и снижение справа, а плотность токов
увлечения jNB монотонно снижается с ростом параметра ITB. Снижение тока
плазмы Ip при увеличении ITB вызвано уменьшением плотности бутстрептока на периферии плазмы. Это ведет к уменьшению энергетического времени E, средних температур и температуры на периферии шнура, и вызывает
уменьшение плотности нейтронного потока n на 20 %.
Изучение влияния положения транспортного барьера проводилось в диапазоне (0,51) значений параметра ITB при заданном коэффициенте увеличения теплопроводности ITB =3.
21
Смещение положения транспортного барьера в направление
магнитной оси ведет к уменьшению бутстреп-тока Ibs и токов
увлечения INB в периферийной
области из-за снижения температуры и ее градиента с внешней
стороны барьера. При дальнейшем смещении барьера к центру
шнура ITB (от 0,7 до 0,5) наблюдается рост бутстреп-тока Ibs и
токов увлечения. При ITB  0,7
Рис. 12. Профили плотностей тока плазмы jp, тодостигается минимум тока плазков увлечения jnb, бутстреп-тока jbs и плотности
плазмы ne при различных значения фактора ITB
мы, что влечет минимум энергетического времени E, и средних температур плазмы. Зависимость нейтронного выхода n от положения транспортного барьера также немонотонная, но
достаточно слабая (~10%).
Результаты Главы 4 показали, что при использовании инжекции быстрых
атомов дейтерия достигается стационарный режим работы токамака JUST-T с
неиндукционным поддержанием тока плазмы, выше уровня, необходимого
для удержанию -частиц в плазме. Удельная мощность нейтронного потока
n превышает величину 0,34 МВт/м2, требуемую для нейтронного источника.
В Заключении сформулированы основные результаты диссертационной
работы:
1. Созданы программные модули расчета коэффициентов переноса для
кода ДИНА, использованные при моделировании и анализе сценариев
разрядов в токамаке КТМ.
2. Разработаны и прошли тестирование на установке TCV «Simulinkмодели плазмы токамака» на основе кода ДИНА с включением в комплекс «компьютерный токамак» для задач управления плазмой токамака.
Результаты тестирования показали приемлемое согласие с данными экспериментов и восстановительного кода LIUQE.
3. Создана методика подготовки сценариев разрядов в токамаке. Методика
применена при разработке сценариев разрядов в токамаке КТМ. Проведены расчеты сценариев разрядов в токамаке КТМ с использованием
различных моделей расчета переноса. Тестирование системы магнитного
управления плазмой в токамаке КТМ показало ее работоспособность в
численном моделировании. Результаты, полученные автором по разработке сценариев разрядов, включены в материалы проектной документации токамака КТМ.
4. Расчеты сценария работы токамака JUST-T показали, что достигается
стационарная стадия развития разряда с неиндукционным поддержанием
22
тока плазмы и удельной мощностью нейтронного потока n выше 0,3
МВт/м2, что достаточно при использовании JUST-T в качестве объемного источника нейтронов в целях трансмутации. Результаты, полученные
в процессе разработки сценариев разряда токамака JUST-T – неиндуктивное поддержание тока, использование комбинации пучков горячих
нейтралов, исследование влияния профилей плотности плазмы и т.п.
могут использоваться при проектировании термоядерных установок токамаков следующего поколения.
Основные результаты работы изложены в следующих публикациях:
1. Raju D., Dokuka V.N., Favez J.-Y., Khayrutdinov R.R. et al. DINA simulations of TCV Electron Cyclotron Current Drive and Heating // Proc. 29 th EPS
Conference on Plasma Physics and Controlled Fusion, 17-21 June 2002. Montreux. ECA Vol. 26B, P-2.082.
2. Lukash V.E., Raju D., Dokuka V.N., Favez J-Y. et al. DINA simulations of
TCV electron cyclotron heating discharge // Fusion Eng. and Design. - 2003. Vol. 66-68 - P. 767-770.
3. Khayrutdinov R.R., Lister J.B., Dokuka V.N., Duval B.P. et al. An Open Architecture Version of the DINA 1.5D Simulation Code // Proc. 30 th EPS Conference on Plasma Physics and Controlled Fusion, 7-11 July 2003. StPetersburg. ECA Vol. 27A, P-3.163.
4. Lukash V.E., Lister J.B., Dokuka V.N., Khayrutdinov R.R. et al. Simulation
of TCV Equilibria Evolution using the DINA Code // Proceeding in 30 th EPS
Conference on Plasma Physics and Controlled Fusion, 7-11 July 2003. StPetersburg. ECA Vol. 27A, P-3.124.
5. Degeling A.W., Martin Y.R., Lister J.B., Dokuka V.N. et al. Magnetic triggering of ELMs in TCV. // Plasma Phys. Control. Fusion. - 2003. - Vol. 45 P. 1637 - 1655.
6. Лукаш В.Э., Докука В.Н., Хайрутдинов Р.Р. Программно-вычислительный комплекс ДИНА в системе MATLAB для решения задач управления
плазмой токамака // ВАНТ. Cер. «Термоядерный синтез». - 2004. - вып.
1. - С. 40 - 49.
7. S.H. Kim, Cavinato M., Dokuka V.N., Khayrutdinov R.R. et al. Comparing
magnetic triggering of ELMs in ASDEX Upgrade and TCV with DINA-CH
tokamak simulator // Proc. 33th EPS Conference on Plasma Physics and Controlled Fusion, 27 June – 1 July 2005. Tarragona, P2-072.
8. Lister J.B., Dokuka V.N., Khayrutdinov R.R., Lukash V.E. et al. Evolution of
the DINA-CH tokamak full discharge simulator // Fusion Eng. аnd Design. 2005. - Vol. 74. - P. 633 - 637.
9. Korotkov V.A., Azizov E.A., Cherepnin Yu.S., Dokuka V.N. et al. KTM Tokamak Conceptual Design and Basic Parameters // Fusion Engineering and
Design. 2001. – Vol. 56-57. P. 831 – 835
23
10. Azizov E.A., Barkalov A.D., Gladush G.G., Dokuka V.N. et al. Control of
the Plasma Fluxes into the Divertor Region of the Tokamak KTM // Proc. 19 th
IAEA Fusion Energy Conference, 14 - 19 Oct. 2002. Lyon. France. FT/P2-12.
11. Dokuka V.N., R.R. Khayrutdinov. Use of DINA code for testing controllers
for KTM device in the framework of Matlab-Simulink environment // Proc.
30th EPS Conference on Plasma Physics and Controlled Fusion. 7 - 11 July
2003. St.-Petersburg. ECA Vol. 27A. P-4.168pd.
12. Azizov E.A., Dokuka V.V., Dvorkin N.Ya., Khayrutdinov R.R., et al. Kazakhstan Tokamak for Material Testing // Plasma Devices and Operations. –
2003. - Vol. 11. - No. 3. P. 39-56.
13. Азизов Э.А., Беляков В.А., Бондарчук Э.Н., Гостев А.А., Докука В.Н. и
др. Казахстанский токамак материаловедческий (КТМ) // ВАНТ, серия
"Электрофизическая аппаратура". - 2005. – вып. 3. - С. 13-18.
14. Докука В.Н., Хайрутдинов Р.Р., Кавин А.А. Синтез и моделирование системы магнитного управления плазмой в токамаке КТМ // ВАНТ, серия
«Термоядерный синтез». - 2008. - вып. 1. - С. 12 - 20.
15. Azizov E.A., Dokuka V.N., Khayrutdinov R.R. Maintenance of stationary
stage of VNS JUST-T for transformation of minor actinides // Proc. 31 st EPS
Conference on Plasma Phys. 28 June – 2 July 2004. London. ECA Vol. 28G.
P-5.038.
16. Azizov E.A., Dokuka V.N., Khayrutdinov R.R. Achievement and maintenance of stationary stage of VNS JUST-T for transmutation of minor actinides
// Proc. 30th EPS Conference on Plasma Physics and Controlled Fusion. 7 - 11
July 2003. St.-Petersburg. ECA Vol. 27A. P-4.169pd.
17. Азизов Э.А., Гладуш Г.Г., Докука В.В., и др. Численное исследование
параметров плазмы токамака-реактора для трансмутации минорных актинидов отработавшего ядерного топлива // ВАНТ, серия «Термоядерный синтез», - 2007. - вып. 3. - С. 3 - 10.
18. Azizov E.A., Arefiev Yu.P., Buzhinskij O.I., Dokuka V.N. et al. Plasmaphysical and Electrophysical Aspects of the Compact Stationary Neutron
Source on Basis of a Tokamak // Plasma Devices and Operations. - 2005. Vol. 13. - No. 3. - P. 167 - 180.
19. Azizov E.A., Dokuka V.N., Khayrutdinov R.R. Stationary Compact VNS
Tokamak for Transmutation // 20th IAEA Fusion Energy Conference. 1-6 November 2004 / Vilamoura. FT/P7-2.
20. Azizov E.A., Arefiev Yu.P., Gladush G.G., Dokuka V.N. et al. The Concept
of the Volumetric Neutron Source on Basic of The JUST-T Tokamak for Minor Actinides Transmutation // Plasma Devices and Operations. - 2003. -Vol.
11 - No. 4. - P. 279 - 286.
21. Azizov E.A., Gladush G.G., Dokuka V.N., Filatov V.V. et al. Investigation of
Minor Actinides Transmutation on Based of Spherical Tokamaks with Aspect
Ratio A=2 // Proc. 19th IAEA Fusion Energy Conference, 14 - 19 Oct. 2002.
Lyon. FT/P1-23.
24
Download