автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Реконструкция эволюции равновесия тороидальной плазмы

На правах рукописи

ООЬи^--

Сучков Егор Петрович

Реконструкция эволюции равновесия тороидальной плазмы

Специальность 05.13.18 - математическое моделирование, численные методы и

комплексы программ

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

1 7 2012

Москва - 2012

005043087

Работа выполнена на кафедре автоматизации научных исследований факультета Вычислительной математики и кибернетики Московского государственного университета имени М.В. Ломоносова.

Научный руководитель: доктор физико-математических наук,

профессор Зайцев Федор Сергеевич.

Официальные оппоненты: доктор физико-математических наук,

профессор Мухин Сергей Иванович кандидат технических наук, Евстигнеев Николай Михайлович

Ведущая организация:

Национальный исследовательский центр "Курчатовский институт"

Защита состоится 30 мая в 15:30 на заседании диссертационного совета Д.501.001. при Московском государственном университете имени М.В. Ломоносова по адресу: 119991, г. Москва, ГСП-1, Ленинские горы, МГУ имени М.В. Ломоносова, 2ой учебный корпус, факультет ВМК, ауд. 685.

С диссертацией можно ознакомиться в научной библиотеке факультета Вычислительной математики и кибернетики МГУ имени М.В. Ломоносова.

Автореферат разослан 26 апреля 2012 года.

Ученый секретарь диссертационного совета, доктор физико-математических наук, профессор

Захаров Евгений Владимирович

Общая характеристика работы

Актуальность работы

В последние десятилетия проводятся интенсивные исследования по создаг аию термоядерного реактора. Одно из наиболее перспективных направлений в эбласти управляемого термоядерного синтеза (УТС) - замкнутые магнитные системы типа токамак. Успехи, достигнутые на установках токамак, позволяют рассматривать их в качестве способного обеспечить потребности человечества источника энергии. В конце 80-х годов прошлого века началось проектирование международного экспериментального термоядерного реактора ITER. (International Thermonuclear Experimental Reactor), целью которого является демонстрация экономической целесообразности использования УТС для производства энергии. Сейчас проект ITER, находится в фазе активного строительства.

Существенный прогресс в исследованиях на установках токамак поставил задачу детального сопоставления эксперимента с теорией. Такое сопоставление требует применения методов математического моделирования, поскольку многие величины, фигурирующие в теории, не измеряются непосредственно в эксперименте. Одной из базовых моделей в проблеме УТС является магнитногид-родинамическая (МГД) модель, позволяющая в соответствующем приближении решать задачу расчёта равновесной конфигурации плазмы.

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

Цели и задачи исследования

В диссертационной работе поставлены и решены следующие задачи:

• Создание методики, позволяющей исследовать с заданной точностью струк-

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

• Разработка метода реконструкции границы плазмы по измерениям магнитного поля.

• Программная реализация вышеописанных методик. Разработка графического интерфейса пользователя, позволяющего облегчить процесс поиска решений обратной задачи.

• Исследование роли данных MSE (motional Stark effect; диагностики, основанной на изменении контуров спектральных линий под воздействием электрического поля, возникающего при взаимодействии заряженных частиц), поляриметрии (диагностики, основанной на измерении степени поляризации света и угла поворота плоскости поляризации при прохождении через плазму) и интерферометрии (диагностики, основанной на измерениях фазы волны, прошедшей через плазму) в задаче реконструкции плотности тока и коэффициента запаса устойчивости.

• Исследование структуры множества решений обратной задачи о реконструкции плотности тока и коэффициента запаса устойчивости с неточно заданными входными данными для параметров установок MAST (Mega Ampere Spherical Tokamak, Великобритания), JET (Joint European Torus, Евросоюз, Великобритания), ITER. (Франция).

• Применение полученных теоретических результатов и созданного програм\ ного обеспечения для обработки экспериментальных данных и исследоваг ния эффективности измерительной аппаратуры на установке JET. Оценка с помощью вычислительного эксперимента необходимой точности измерений и интервала доверия реконструкций.

Научная новизна

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

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

Важнейшей задачей проблемы удержания плазмы в тороидальных системах является задача определения формы и положения плазмы. Информация о границе плазмы используется в качестве входной многими кодами, рассчитывающими или восстанавливающими характеристики плазмы, а также системами управления плазмой. Задача восстановления границы плазмы по измерениям магнитного поля также является некорректной, тем не менее разработаны различные методы ее решения, например, реализованные в кодах XLOC и EFIT. Однако код EFIT решает слишком общую задачу, восстанавливая не только границу, но и внутренние параметры плазмы. При этом для достижения достаточной точности в дополнение к магнитным измерениям алгоритм кода EFIT требует данные других диагностик. Недостатком кода XLOC является существенное использование конструктивных особенностей установки JET и типичной формы плазмы. В диссертационной работе представлена формулировка общей математической задачи о восстановлении границы плазмы по магнитным

измерениям и изложен алгоритм ее решения.

Теоретическая и практическая значимость

• Предложены новые постановки задач о реконструкции эволюции равнове-j сия тороидальной плазмы, разработаны и обоснованы методы их решения.

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

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

• Полученные теоретические результаты и созданное программное обеспечение используется на установке JET для анализа экспериментов.

Апробация работы

Основные результаты диссертации докладывались на следующих конференциях и научно-исследовательских семинарах:

1. Международная конференция "35th EPS Plasma Physics Conference" (о. Крит, Греция, 2008 год).

2. VI Курчатовская молодежная школа (г. Москва, 2008 год).

3. Международная летняя школа "ITER" (г. Экс-ан-Прованс, Франция, 2009 год).

4. XVII Международная конференция "Ломоносов-2010" (г. Москва, 2010 год).

5. XI Международный семинар "Супервычисления и математическое моделирование" (г. Саров, 2009 год).

6. XII Международный семинар "Супервычисления и математическое моделирование" (г. Саров, 2010 год).

7. Семинар ИКИ РАН "Вычислительные технологии в естественных науках. Перспективные компьютерные системы: устройства, методы и концепции" (г. Таруса, 2011 год).

8. Международная конференция "7-th Workshop on Data Processing Validation and Analyses" (Фраскати, Италия, 2012 год).

Публикации

По теме диссертации опубликовано 11 научных работ, в том числе 3 - в статьях журналов списка ВАК. Список публикаций по теме диссертации приведен в конце автореферата.

Структура и объем работы

Диссертация состоит из введения, четырех глав, заключения и списка цитируемой литературы. Объем работы - 113 страниц.

Содержание работы

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

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

В первом разделе рассматривается самосогласованная модель эволюции равновесия плазмы. Дается математическая постановка задачи, включающая уравнение равновесия Грэда-Шафранова с эллиптическим оператором и нелинейной правой частью,

{_9p(i, ф) 1 dF2(t, ф) „ м

дф 2 шй дф

А ч

Ji(t)i(A - (Z - Я) вне Tp(t)

усредненную продольную (по отношению к магнитному полю) компоненту обобщённого закона Ома и дополнительные ограничения. Использованы следующие обозначения: (R,r},Z) - цилиндрические координаты с осью Z, направленной вдоль оси аксиальной симметрии; ф{Ь, R, Z) - функция полоидального потока, равная ковариантной компоненте Av векторного потенциала магнитного поля В — V х А; Р(ф) - функция полоидального тока, связанная с тороидальным магнитным полем Bt0T = F/R; р&ф) - кинетическое давление; Гр(£) - граница плазмы, определяемая как замкнутая магнитная поверхность ф ~ const максимальной ширины; Цг) - токи в соленоиде, катушках полоидального поля и стенках камеры; jv(t, R, Z) - плотность тороидального тока. Искомыми функциями являются ф,р,Р и Гр.

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

Третий раздел главы посвящен описанию задачи реконструкции плотности тороидального тока. Задача восстановления плотности тока в плазме по измерениям магнитного поля является сильно некорректной. Помимо отсутствия или

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

В последнем разделе первой главы рассматривается задача реконструкции равновесия плазмы с учетом закона Ома. В разделе приводится математическая формулировка задачи, включающая уравнение ГрэдагШафранова, усредненный продольный закон Ома и ряд дополнительных ограничений.

Во второй главе представлены разработанные численные методы решения обратных задач эволюции равновесия. В диссертации задача реконструкции эволюции равновесия разбивается на две: восстановление границы плазмы и реконструкция параметров внутри плазмы.

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

В данной части второй главы изложены алгоритмы построения е-сетей полиномов и сплайнов. Для полиномов предложены два способа построения сети: первый основан на рассмотрении специальных наборов коэффициентов полиномов; второй - на построении сетки таким образом, чтобы проходящие через ее

1 А.Н. Колмогоров, В.М. Тихомиров. "Е-энтропия и е-емкость множеств в функциональных пространствах". Успехи мат. наук. 1956, т. XIV, вып. 2(86), с. 3-86.

узлы полиномы образовывали е-сеть. Также описан способ построения комбинированной е-сети полиномов. Данный способ заключается в использовании двух вышеописанных методов, каждый из которых применяется к определённой чат сти исследуемого множества, а не ко всему множеству. Для сплайнов предложен способ построения е-сети на сетке.

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

В заключении раздела приводятся стратегии применения е-сетей. Описываются способы выбора полосы, внутри которой происходит построение сети, степени полинома или числа узлов сетки для сплайнов, а также точности е. Приводятся различные способы конструирования е-сетей полиномов или сплайнов с некоторыми заранее заданными свойствами.

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

Функция тр ищется в виде многочлена

N М

= P = R2-B%, z = Z-ZO,

i=0 о

а граница плазмы - в параметрическом виде

{ R = Щ£) . г , 1 v, £ е [-к,ж].

Уравнения для определения границы плазмы включают следующие условия:

значения функции полоидального потока ф в точках измерения, уравнение Грэда-Шафранова вне плазмы, условия для ф на границе плазмы, вблизи стенок камеры и пассивных элементов. Для каждой кривой Z{£)) из рассматриваемого класса данная система уравнений дает возможность определить функцию ip(R,Z) вне плазмы. Для её решения применяется метод регуляризации Тихонова. В качестве искомой границы принимается та, на которой достигается минимум невязки системы.

В общем случае перебор границ осуществляется с помощью специальной е-сети сплайнов, построенной отдельно для и для Z(£). В разделе проанализирована точность общего метода и простейшей трехпараметрической аппроксимации границы

Г = i?o + 7 cos £ — 6 sin2 £ \ Z(£) — Zq + A7sin£,

где 7 - полуширина границы плазмы, А - её эллиптичность, д - треугольность.

Третий раздел главы посвящен изложению алгоритма восстановления равновесной плотности тороидального тока. Алгоритм разбивается на три этапа. На первом этапе определяется правая часть уравнения ГрэдагШафралова на основе закона Ампера. На втором этапе восстанавливается давление р и функция полоидального тока F с использованием теории е-сетей. На третьем этапе по найденным производным р и F вычисляем значения равновесной плотности тороидального тока.

Принципиально новым в алгоритме является использование информации о положении границы плазмы Гр, найденной с помощью алгоритма, предложенного в предыдущем разделе данной главы.

Далее в диссертационной работе приводится пример расчета плотности тороидального тока для параметров типичного разряда в установке MAST с использованием вышеизложенного алгоритма. Сравнение с решением прямой задачи, показало, что алгоритм решения обратной задачи обладает приемлемой точностью. Изложенный алгоритм использует данные о положении границы,

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

В следующих двух разделах описываются алгоритмы восстановления напряжения тороидального электрического поля и равновесия с учетом уравнения эволюции. В задаче восстановления напряжения рассматривается уравнение, полученное из уравнение ГрэдагШафранова взятием от обеих частей производной по времени £ при фиксированной границе плазмы. Метод решения данной задачи аналогичен методу, применявшемуся на первом этапе алгоритма восстановления равновесной плотности тока. В задаче восстановления равновесия с учетом уравнения эволюции, помимо уравнения Грэда-Шафранова, учитывается также усредненный продольный закон Ома, на основе которого восстанавливается функция полоидального тока Р. В шестом разделе второй главы приводится описание алгоритма решения данной задачи.

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

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

Предложенные численные алгоритмы реализованы в программном комплек-

ce SDSS, написанном на языках Fortran 2003 и Java. Описанию данного комплекса посвящена третья глава диссертации.

В первом разделе третьей главы сформулированы функциональные возможности кода SDSS.

Далее в диссертации приводится описание графического интерфейса пользователя, реализованного на языке Java. Данный интерфейс позволяет автоматизировать построение и визуализацию е-сетей, задание объемных входных данных и запуск кода SDSS.

В третьем и четвертом разделах приводится описание входных и выходных данных кода SDSS, представлена методика его использования, изложенная в виде поэтапого руководства.

Последний раздел третьей главы посвящен описанию организации параллельных вычислений. В данном разделе продемонстрированы преимущества использования языка Фортран 2003 для организации параллельных вычислений при решении задач проблемы У ТС. Также в разделе приводится описание технологии MPI и сформулированы некоторые рекомендации по её использованию.

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

В первом разделе продемонстрированы результаты восстановления плотности тока и коэффициента запаса устойчивости с помощью алгоритма, описанного в седьмом разделе второй главы, для параметров установок MAST, JET и ITER. Построение решения проводилось с помощью кода SDSS. Для каждой из рассмотренных установок показано, что использования данных только магнитной диагностики недостаточно для однозначной идентификации физического процесса внутри плазмы. В рамках заданной погрешности входных данных, условиям задачи могут удовлетворять существенно различные (как качественно, так и количественно) решения, соответствующие различным режимам удержания плазмы.

В диссертации исследована роль диагностики MSE в выделении решения, адекватного физическому процессу. Во втором разделе показано, что использование MSE-ограничения, заданного с погрешностью, сравнимой с экспериментальной, позволяет выделить решение, соответствующее реальному физическому процессу. Ширина интервала доверия для плотности тороидального тока с использованием диагностики MSE оказалась сравнима с погрешностями измерения магнитных диагностик и MSE.

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

MSE.

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

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

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

Основные результаты

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

2. Новые подходы реализованы в программном комплексе SDSS, написанном на языках Fortran 2003 и Java. Комплекс имеет развитый графический интерфейс пользователя, позволяющий автоматизировать процесс поиска решений обратной задачи и визуализировать большое количество связанных с ней массивов данных.

3. Исследована структура множества решений задачи о реконструкции плотности тока и коэффициента запаса устойчивости с различными дополнительными ограничениями для параметров установок MAST, JET и ITER. Показано, что обратная задача, основанная только на магнитных измерениях, является сильно неустойчивой по входным данным и имеет решения, соответствующие различным режимам удержания плазмы. Изучена роль измерений MSE (motional Stark effect), поляриметрии и интерферометрии в выделении решения.

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

Публикации по теме диссертации

(Научные статьи в журналах и изданиях, включенных в перечень ВАК, выделены курсивом):

1. Д.П. Костомаров, Ф.С. Зайцев, Е.П. Сучков "Построение сильно различающихся решений некоторого класса некорректных задач с неточно заданными входными данными", ДАН, 2011, т-437, №, с. 316-S20.

2. Д.П. Костомаров, Ф.С. Зайцев, А.Г. Шишкин, Д.Ю. Сычугов, С.В. Степанов. Е.П. Сучков. Программное обеспечение библиотеки "Виртуальный токомак". // Вестник Моск. ун-та. Сер. 15, Вычислительная ма-тем. и киберн. 2011, с. 48-543. F.S. Zaitsev, D.P. Kostomarov, Е.Р. Suchkov, V.V. Drozdov, E.R. Solano, A.

Murari, S. Matejcik, N.C. Hawkes and JET-EFDA Contributors. "Analyses of substantially different plasma current densities and safety factors reconstructed from magnetic diagnostics data", Nuclear Fusion, 2011, v. 51, 103044, 10 p.

4. F.S. Zaitsev, S. Matejcik, A. Murari, E.P. Suchkov and JET EFDA Contributors A new method to identify the equilibria compatible with the measurements using the technique of the e-nets. // 7-th Workshop on Data Processing Validatic and Analyses. Frascati, Italy. March 2012, 1 p.

5. Ф.С. Зайцев, Д.Ю. Сычугов, А.Г. Шишкин, В.Э. Лукаш, Ю.В. Митриш-кин, P.P. Хайрутдинов, В.Н. Докука, И.Б. Семенов, А.А. Лукьяница, И.В. Зотов, В.В. Нефёдов, С.В. Степанов, Е.П. Сучков, С.А. Унучек. Концепция комплекса имитационного моделирования «Виртуальный токамак» с системами управления плазмой. // XII Международный семинар "Супервычисления и математическое моделирование 11-15 октября 2010. - Саров: РФЯЦ-ВНИИЭФ, 2011, с. 194-195.

6. Ф.С. Зайцев, А.Г. Шишкин, Д.Ю. Сычугов, В.Э. Лукаш, Ю.В. Митриш-кин, P.P. Хайрутдинов, С.В. Степанов, Е.П. Сучков. Структура и функци-

ональные возможности комплекса имитационного моделирования "Виртуальный токамак". // Вычислительные технологии в естественных науках. Перспективные компьютерные системы: устройства, методы и концепции. М.: ИКИ РАН, 2011, с. .

7. Ф.С. Зайцев, Д.П. Костомаров, Е.П. Сучков "Восстановление границы тороидальной плазмы по данным магнитной диагностики". Труды XII Международного семинара "Супервычисления и математическое моделирование", г. Саров, 2010, с.193-194.

8. Е.П. Сучков "Алгоритм построения сильно различающихся решений некоторого класса обратных задач". Сборник тезисов XVII Международной научной конференции "Ломоносов-2010", секция "Вычислительная математика и кибернетика", Москва, 2010, с. 149-150.

9. Ф.С. Зайцев, Д.П.Костомаров, Е.П. Сучков. Анализ существенно различных решений обратной задачи магнитной диагностики тороидальной плазмы. // XI Международный семинар "Супервычисления и математическое моделирование". - Саров: РФЯЦ. Октябрь 2009. С. 67.

10. F.S. Zaitsev, D.P. Kostomarov, Е.Р. Suchkov "Existence of substantially different solutions in an inverse problem of plasma equilibrium reconstruction". 35th EPS Plasma Physics Conference. Crete, Greece, 2008. P-1.091, 4p.

11. Ф.С. Зайцев, Е.П. Сучков. Существование решений обратной задачи для уравнения Грэда-Шафранова с сильно различающимися q. // 6-я Курчатовская молодёжная школа. - Москва: РНЦ "Курчатовский институт". Ноябрь 2008. С. 93-94.

Напечатано с готового оригинал-макета

Издательство ООО "МАКС Пресс" Лицензия ИД N 00510 от 01.12.99 г. Подписано в печать 26.04.2012 г. Формат 60x90 1/16. Усл.печл. 1,0 Тираж 100 экз. Заказ 203. Тел. 939-3890. Тел./факс 939-3891. 119992, ГСП-2, Москва, Ленинские горы, МГУ им. М.В. Ломоносова, 2-й учебный корпус, 527 к.

61 12-1/793

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ

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

ФАКУЛЬТЕТ ВЫЧИСЛИТЕЛЬНОЙ МАТЕМАТИКИ

И КИБЕРНЕТИКИ

Реконструкция эволюции равновесия тороидальной плазмы

Специальность 05.13.18 - математическое моделирование, численные методы и

комплексы программ

На правах рукописи

Сучков Егор Петрович

Диссертация на соискание ученой степени кандидата физико-математических наук

Научный руководитель: д.ф-м.н. Зайцев Федор Сергеевич

Москва - 2012

Оглавление

Введение ...............................................................3

1 Формулировки математических задач эволюции равновесия .... 15

1.1 Самосогласованные модели эволюции равновесия плазмы . 15

1.2 Задача о восстановлении границы плазмы....................22

1.3 Задача о реконструкции плотности тороидального тока . . 23

1.4 Реконструкция равновесия с учетом закона Ома............26

2 Численные методы решения обратных задач эволюции равновесия 28

2.1 Построение сильно различающихся решений некорректной задачи методом е-сетей ........................................28

2.2 Реконструкция границы плазмы по данным магнитной диагностики ...............................................49

2.3 Реконструкция равновесной плотности тока..................53

2.4 Пример расчета и анализ точности ..............' 58

2.5 Реконструкция напряжения тороидального электрического поля ..........................................................60

2.6 Реконструкция эволюции равновесия ........................61

2.7 Построение существенно различных решений ..............62

2.8 Распараллеливание алгоритмов ..............................66

3 Программное обеспечение ........................................69

3.1 Концепция и функциональные возможности кода SDSS . 69

3.2 Графический интерфейс пользователя........................70

3.3 Входные и выходные данные..................................71

3.4 Технология использования ....................................73

3.5 Организация параллельных вычислений......................74

4 Результаты вычислительного эксперимента..........................81

4.1 Построение существенно различных решений обратной задачи для параметров установок MAST, JET и ITER ... 81

4.2 Роль измерений MSE в задаче реконструкции ..............83

4.3 Влияние данных поляриметрии и интерферометрии .... 86

4.4 Сравнение значимости данных MSE, поляриметрии и ин-треферометрии............................................87

4.5 Реконструкция эволюции равновесия в установке MAST . 88

Заключение ...............................................90

Литература.............................................................92

Приложение №1 "Иллюстрации"......................................99

Введение

Актуальность работы

В последние десятилетия проводятся интенсивные исследования по созданию термоядерного реактора. Одно из наиболее перспективных направлений в области управляемого термоядерного синтеза (УТС) - замкнутые магнитные системы типа токамак [1, 2]. Успехи, достигнутые на установках токамак, позволяют рассматривать их в качестве способного обеспечить потребности человечества источника энергии. В конце 80-х годов прошлого века началось проектирование международного экспериментального термоядерного реактора ITER (International Thermonuclear Experimental Reactor), целью которого является демонстрация экономической целесообразности использования УТС для производства энергии [3]. Сейчас проект ITER, находится в фазе активного строительства.

Существенный прогресс в исследованиях на установках токамак поставил задачу детального сопоставления эксперимента с теорией. Такое сопоставление требует применения методов математического моделирования, поскольку многие величины, фигурирующие в теории, не измеряются непосредственно в эксперименте. Одной из базовых моделей в проблеме УТС является магнитногид-родинамическая (МГД) модель, позволяющая в соответствующем приближении решать задачу расчёта равновесной конфигурации плазмы.

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

Цели и задачи исследования

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

• Разработка метода реконструкции границы плазмы по измерениям магнитного поля.

• Программная реализация вышеописанных методик. Разработка графического интерфейса пользователя, позволяющего облегчить процесс поиска решений обратной задачи.

• Исследование роли данных MSE (motional Stark effect; диагностики, основанной на изменении контуров спектральных линий под воздействием электрического поля, возникающего при взаимодействии заряженных частиц) [4, 5], поляриметрии (диагностики, основанной на измерении степени поляризации света и угла поворота плоскости поляризации при прохождении через плазму) и интерферометрии (диагностики, основанной на измерениях фазы волны, прошедшей через плазму) [4, 6] в задаче реконструкции плотности тока и коэффициента запаса устойчивости.

• Исследование структуры множества решений обратной задачи о реконструкции плотности тока и коэффициента запаса устойчивости с неточно заданными входными данными для параметров установок MAST (Mega Ampere Spherical Tokamak, Великобритания), JET (Joint European Torus, Евросоюз, Великобритания), ITER (Франция).

• Применение полученных теоретических результатов и созданного программного обеспечения для обработки экспериментальных данных и исследова-

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

Научная новизна

Одной из центральных задач в проблеме УТС является задача восстановления плотности тока в плазме по измерениям магнитного поля. Однако эта задача является сильно некорректной. Определение корректности [7] может нарушаться сразу по всем трем пунктам. Помимо отсутствия или неединственности решения может иметь место неустойчивость по входным данным, т.е. близким измерениям могут соответствовать существенно различные плотности тока.

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

Важнейшей задачей проблемы удержания плазмы в тороидальных системах является задача определения формы и положения плазмы. Информация о границе плазмы используется в качестве входной многими кодами, рассчитывающими или восстанавливающими характеристики плазмы, а также системами управления плазмой. Задача восстановления границы плазмы по измерениям магнитного поля также является некорректной, тем не менее разработаны различные методы ее решения, например, реализованные в кодах XLOC[8 - 10] и

ЕИТ[11 - 15]. Однако код ЕИТ решает слишком общую задачу, восстанавливая не только границу, но и внутренние параметры плазмы. При этом для достижения достаточной точности в дополнение к магнитным измерениям алгоритм кода ЕЕ1Т требует данные других диагностик. Недостатком кода ХЬОС является существенное использование конструктивных особенностей установки ЛЕТ и типичной формы плазмы. В диссертационной работе представлена формулировка общей математической задачи о восстановлении границы плазмы по магнитным измерениям и изложен алгоритм ее решения.

Предмет защиты составляют следующие положения:

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

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

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

Разработан программный комплекс, реализующий предложенные алгоритмы.

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

зано, что при достаточной точности дополнительных измерений их использование даёт возможность выделить одно решение обратной задачи с неточными входными данными, соответствующее реальному физическому процессу.

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

Апробация работы

Основные результаты диссертации были представлены на следующих конференциях и семинарах:

• Международная конференция "35th EPS Plasma Physics Conference", о. Крит, Греция, 9-13 июня 2008 года.

• VI Курчатовская молодежная школа, г. Москва, 17-19 ноября 2008 года.

• Международная летняя школа "ITER", г. Экс-ан-Прованс, Франция, 22-26 июня 2009 года.

• XVII Международная конференция "Ломоносов-2010", г. Москва, 12-15 апреля 2010 года.

• XI Международный семинар "Супервычисления и математическое моделирование", г. Саров, октябрь 2009 года.

• XII Международный семинар "Супервычисления и математическое моделирование", г. Саров, октябрь 2010 года.

• Семинар ИКИ РАН "Вычислительные технологии в естественных науках. Перспективные компьютерные системы: устройства, методы и концепции", Таруса, март 2011 г.

• 7-th Workshop on Data Processing Validation and Analyses. Frascati, Italy. March 2012

Публикации

Результаты диссертационной работы опубликованы в следующих статьях:

1. Д.П. Костомаров, Ф.С. Зайцев, Е.П. Сучков "Построение сильно различающихся решений некоторого класса некорректных задач с неточно заданными входными данными", ДАН, 2011, т.437, №3, с. 316-320.

2. Д.П. Костомаров, Ф.С. Зайцев, А.Г. Шишкин, Д.Ю. Сычугов, С.В. Степанов. Е.П. Сучков. Программное обеспечение библиотеки "Виртуальный токамак". // Вестник Моск. ун-та. Сер. 15, Вычислительная матем. и ки-берн. 2011, с. 48-54.

3. F.S. Zaitsev, D.P. Kostomarov, Е.Р. Suchkov, V.V. Drozdov, E.R. Solano, A. Murari, S. Matejcik, N.C. Hawkes and JET-EFDA Contributors. "Analyses of substantially different plasma current densities and safety factors reconstructed from magnetic diagnostics data", Nuclear Fusion, 2011, v. 51, 103044, 10 p.

4. F.S. Zaitsev, S. Matejcik, A. Murari, E.P. Suchkov and JET EFDA Contributors. A new method to identify the equilibria compatible with the measurements using the technique of the e-nets. // 7-th Workshop on Data Processing Validatio: and Analyses. Frascati, Italy. March 2012, 1 p.

5. Ф.С. Зайцев, Д.Ю. Сычугов, А.Г. Шишкин, В.Э. Лукаш, Ю.В. Митриш-кин, P.P. Хайрутдинов, В.Н. Докука, И.Б. Семенов, A.A. Лукьяница, И.В. Зотов, В.В. Нефёдов, C.B. Степанов, Е.П. Сучков, С.А. Унучек. Концепция комплекса имитационного моделирования «Виртуальный токамак» с системами управления плазмой. // XII Международный семинар "Супервычисления и математическое моделирование 11-15 октября 2010. - Саров: РФЯЦ-ВНИИЭФ, 2011, с. 194-195.

6. Ф.С. Зайцев, А.Г. Шишкин, Д.Ю. Сычугов, В.Э. Лукаш, Ю.В. Митриш-кин, P.P. Хайрутдинов, C.B. Степанов, Е.П. Сучков. Структура и функциональные возможности комплекса имитационного моделирования " Виртуальный токамак". // Вычислительные технологии в естественных науках. Перспективные компьютерные системы: устройства, методы и концепции. М.: ИКИ РАН, 2011, с. .

7. Ф.С. Зайцев, Д.П. Костомаров, Е.П. Сучков "Восстановление границы тороидальной плазмы по данным магнитной диагностики". Труды XII Международного семинара "Супервычисления и математическое моделирование", г. Саров, 2010, с.193-194.

8. Е.П. Сучков "Алгоритм построения сильно различающихся решений некоторого класса обратных задач". Сборник тезисов XVII Международной научной конференций "Ломоносов-2010", секция "Вычислительная математика и кибернетика", Москва, 2010, с. 149-150.

9. Ф.С. Зайцев, Д.П.Костомаров, Е.П. Сучков. Анализ существенно различных решений обратной задачи магнитной диагностики тороидальной плазмы. //XI Международный семинар "Супервычисления и математическое моделирование". - Саров: РФЯЦ. Октябрь 2009. С. 67.

10. F.S. Zaitsev, D.P. Kostomarov, E.P. Suchkov "Existence of substantially different solutions in an inverse problem of plasma equilibrium reconstruction". 35th EPS Plasma Physics Conference. Crete, Greece, 2008. P-1.091, 4p.

11. Ф.С. Зайцев, Е.П. Сучков. Существование решений обратной задачи для уравнения Грэда-Шафранова с сильно различающимися q. // 6-я Курчатовская молодёжная школа. - Москва: РНЦ "Курчатовский институт". Ноябрь 2008. С. 93-94.

Структура и объем

Диссертация состоит из введения, четырех глав, заключения и списка цитируемой литературы.

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

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

Во второй главе описываются разработанные численные методы решения обратных задач эволюции равновесия. В первой части данной главы излагается метод поиска сильно различающихся решений некорректной задачи с использованием теории £-сетей, описываются способы построения £-сетей. В разделе дана общая формулировка задач, представленных в первой главе, и описан алгоритм их решения. Следующая часть главы посвящена изложению алгоритма решения общей задачи о восстановлении границы плазмы по магнитным измерениям. В третьей части второй главы описывается алгоритм восстановления равновесной плотности тороидального тока. Алгоритм излагается поэтапно -сначала описывается алгоритм реконструкции тороидальной компоненты тока на основе закона Ампера, затем описывается алгоритм нахождения давления плазмы р и функции полоидального тока F, и последним этапом, учитывающим предыдущие два, является описание алгоритма восстановления непосредственно равновесной плотности тока. Далее приводится пример моделирования типичного разряда в установке MAST. Проводится анализ точности расчета, подтверждающий пригодность предложенных алгоритмов для практического использования. В разделах 2.5 и 2.6 описываются алгоритмы реконструкции напряжения тороидального электрического поля и эволюции равновесия, соответственно. Седьмой раздел второй главы посвящен изложению способа построения существенно различных решений, в нём также приводятся критерии их отбора, излагается способ построения доверительных интервалов. В последнем разделе второй главы описываются способы распараллеливания предложенных алгоритмов построения s-сетей и реконструкции эволюции равновесия.

В третьей главе дается описание программной системы, реализующей предложенные в предыдущей главе алгоритмы. В первом разделе третьей главы

сформулированы функциональные возможности кода SDSS. Следующий раздел описывает графический интерфейс пользователя, реализованный на языке Java, который позволяет автоматизировать построение и визуализацию е-сетей, задание объемных входных данных и запуск кода SDSS. В разделах 3.3 и 3.4 приводится описание входных и выходных данных кода SDSS, а также представлена методика его использования. Последний раздел третьей главы посвящен описанию организации параллельных вычислений. В данном разделе продемонстрированы преимущества использования языка Фортран 90/95/2003 для организации параллельных вычислений при решении задач проблемы УТС. Также в разделе приводится описание технологии MPI и сформулированы некоторые рекомендации по ее использованию.

В четвертой главе продемонстрированы результаты восстановления плотности тока и коэффициента запаса устойчивости с помощью алгоритма, описанного в седьмом разделе второй главы, для параметров установок MAST, JET и ITER. Построение решения проводилось с помощью кода SDSS. Для каждой из рассмотренных установок показано, что использования данн�