RU2787571C1 - Способ быстродействующего определения формы плазмы в камере токамака в течение диверторной фазы плазменных разрядов - Google Patents

Способ быстродействующего определения формы плазмы в камере токамака в течение диверторной фазы плазменных разрядов Download PDF

Info

Publication number
RU2787571C1
RU2787571C1 RU2022117882A RU2022117882A RU2787571C1 RU 2787571 C1 RU2787571 C1 RU 2787571C1 RU 2022117882 A RU2022117882 A RU 2022117882A RU 2022117882 A RU2022117882 A RU 2022117882A RU 2787571 C1 RU2787571 C1 RU 2787571C1
Authority
RU
Russia
Prior art keywords
plasma
tokamak
signals
magnetic
separatrix
Prior art date
Application number
RU2022117882A
Other languages
English (en)
Inventor
Юрий Владимирович Митришкин
Артем Евгеньевич Коньков
Павел Сергеевич Коренев
Валерий Игоревич Кружков
Original Assignee
Федеральное государственное бюджетное образовательное учреждение высшего образования "Московский государственный университет имени М.В.Ломоносова" (МГУ)
Filing date
Publication date
Application filed by Федеральное государственное бюджетное образовательное учреждение высшего образования "Московский государственный университет имени М.В.Ломоносова" (МГУ) filed Critical Федеральное государственное бюджетное образовательное учреждение высшего образования "Московский государственный университет имени М.В.Ломоносова" (МГУ)
Application granted granted Critical
Publication of RU2787571C1 publication Critical patent/RU2787571C1/ru

Links

Images

Abstract

Изобретение относится к способу определения формы плазмы в камере токамака в течение диверторной фазы плазменных разрядов. На первом этапе способа по ряду уже проведенных плазменных разрядов и входо-выходным сигналам вне плазмы восстанавливают равновесие плазмы, т.е. распределение полоидального потока и тороидального тока. На основе этих распределений определяют положение плазменной сепаратрисы и тем самым определяют соответствие между входными сигналами (на магнитных петлях и зондах Мирнова, на поясах Роговского, измеряющих ток плазмы и токи в обмотках полоидального поля) и выходными сигналами, которыми являются расстояния между сепаратрисой и первой стенкой токамака (зазоры). На втором этапе настраивают средства оценки зазоров: статическую матрицу, наблюдатель состояния, нейронную сеть по входным и выходным магнитным сигналам на ряде разрядов в течение каждого разряда. После этого применяют данные средства для оценки зазоров и ударных точек в темпе наблюдений (в реальном времени) в течение разрядов плазмы, что позволяет оцененные сигналы применять для управления формой плазмы в обратной связи системы магнитного управления плазмой. Техническим результатом является возможность провести быстрое определение формы плазмы на диверторной фазе плазменных разрядов. 2 з.п. ф-лы, 9 ил.

Description

Область техники, к которой относится изобретение
Изобретение относится к способам магнитной диагностики и магнитного управления плазмой в D-образных токамаках, а именно к способам восстановления равновесия плазмы по магнитным измерениям вне плазмы в реальном времени в токамаках (тороидальных камер с магнитными катушками). Заявляемый способ предназначен для быстродействующего определения состояния равновесия плазмы в D-образных токамаках и позволяет оценить зазоры между сепаратрисой и первой стенкой токамака, а также смещения ударных точек на диверторной фазе разряда в реальном времени. Зазоры и смещения ударных точек определяют форму плазмы.
В рамках заявляемого изобретения использованы следующие специальные термины:
Равновесие плазмы - состояние плазмы в токамаке, задаваемое распределениями полоиального магнитного потока и плотности тороидального тока.
Первая стенка токамака - внутренняя стенка камеры или бланкета токамака, ограничивающая плазму.
Диверторная фаза разряда - фаза разряда плазмы токамака, в которой граница плазмы не касается первой стенки токамака. Характеризуется наличием точки нулевого магнитного поля (X-точки) на границе плазмы и ударных (диверторных) точек на камере токамака или на диверторных пластинах, в которые попадают вылетающие из X-точки частицы.
Сепаратриса - термин для границы плазмы в течение диверторной фазы разряда. Может быть найдена как наибольшая замкнутая линия равного уровня полоидального потока.
Уровень техники
Конструкция токамака (тороидальная камера с магнитными катушками) представляет собой трансформатор, в первичную обмотку которого поступает импульс тока от внешнего источника энергии, а плазменный виток является вторичным витком трансформатора, в котором протекает ток, нагревающий плазму.
Для регулирования формы плазмы в камере токамака необходимо знать распределение плотности тороидального тока, полоидального потока, форму и положение плазмы. Однако эти параметры не могут быть измерены непосредственно в ходе эксперимента и должны быть определены из измеряемых значений полоидального потока вне плазмы, полного тока плазмы и токов в обмотках полоидального поля токамака. Данная задача известна как задача восстановления равновесия плазмы. Имея значения указанных магнитных сигналов, с учетом известных уравнений физики плазмы (в т.ч. уравнение Града-Шафранова), можно восстановить значения параметров внутри плазмы (в частности, распределение тороидального тока). В общем случае, известны такие решения, обеспечивающие восстановление плазмы, однако все они характеризуются значительным (в рамках времени диверторной фазы плазменных разрядов) временем вычисления, что затрудняет применение таких решений при необходимости управления плазмой в режиме реального времени (в темпе наблюдения).
Так, из уровня техники известны следующие средства и методы определения распределения полоидального магнитного потока и тороидального тока по магнитным измерениям вне плазмы.
В источнике (Lao L., John H., Stambaugh R., Kellman A., Preiffer W. Reconstruction of current profile parameters and plasma shapes in tokamaks // Nuclear Fusion. 1985. Vol. 25, N. 11. P. 1611-1622) раскрыт способ решения обратной краевой задачи, плохо обусловленной по Адамару, позволяющий посредством итераций Пикара решать уравнение Града-Шафранова и найти распределение полоидального потока и тороидального тока. Эти распределения дают возможность определить расположение плазменной сепаратрисы и восстановить положение плазмы в камере.
В источнике (Y.V. Mitrishkin, A.A. Prokhorov, P.S. Korenev and M.I. Patrov. Hierarchical robust switching control method with the equilibrium reconstruction code based on improved Moving Filaments approach in the feedback for tokamak plasma shape. Fusion Engineering and Design, 2019, v. 138, pp. 138-150.
https://doi.org/10.1016/j.fusengdes.2018.10.031) раскрыт способ определения равновесия плазмы при помощи аппроксимации тока плазмы подвижными токовыми кольцами (филаментами).
Как уже отмечено, описанные аналоги не позволяют реализовать способ, обеспечивающий максимально возможное быстродействие при определении положения сепаратрисы плазмы: из-за итерационного характера определения равновесия посредством итераций Пикара или из-за процедуры вычисления положения филаментов и токов в них. Таким образом, использование известных способов в режиме реального времени затруднительно.
Источник (Braams B.J., Jilge W. and Lackner K., et al. Fast determination of plasma parameters through function parametrization // Nuclear Fusion. - 1986. - Vol. 26, N 6. - P. 699-708) раскрывает способ регулирования формы плазмы, основанный не на решении уравнения равновесия плазмы, а на регрессионном анализе большого количества экспериментальных данных. Это не дает физического понимания получаемых решений, что затрудняет применение результата для управления формой плазмы в реальном времени.
Наиболее близким к заявляемому изобретению по технической сущности является способ определения формы плазмы в D-образном токамаке, направленный на восстановление равновесия плазмы, т.е. распределения полоидального потока и тока плазмы, раскрытый в источнике (Коренев П.С., Митришкин Ю.В., Патров М.И. Реконструкция равновесного распределения параметров плазмы токамака по внешним магнитным измерениям и построение линейных плазменных моделей. Мехатроника, автоматизация и управление. Том 17, №4, 2016, с. 254-265. DOI: 10.17587/mau.17.254-266).
Однако данный способ имеет относительно низкое быстродействие (время вычислительного шага при использовании данного способа составляет около 200 микросекунд (мкс)), т.к. для восстановления равновесия необходимо применять итерации Пикара при решении уравнения Града-Шафранова.
Технической проблемой заявляемого изобретения является необходимость преодоления технических недостатков, присущих аналогам, за счет создания быстродействующего способа, позволяющего оценить положение и форму плазмы в реальном времени.
Раскрытие изобретения
Техническим результатом заявленного изобретения является сокращение времени определения величины зазоров между первой стенкой токамака и сепаратрисой плазмы, а также смещение ударных точек сепаратрисы в течение диверторной фазы в режиме реального времени.
При известных значениях вычислительного шага прототипа в 200 мкс в результате применения заявляемого изобретения удалось достичь сокращения времени указанного шага до значений 5-50 мкс (в зависимости от используемого типа токамака), что позволит более оперативно реагировать на изменения положения плазмы посредством подачи корректирующего воздействия, а также освобождает дополнительное время внутри каждого шага дискретного управления для применения более сложных и совершенных алгоритмов управления формой плазмы.
Заявляемый технический результат достигается тем, что способ быстродействующего определения формы плазмы в камере токамака в течение диверторной фазы плазменных разрядов включает
- предоставление данных входных экспериментальных сигналов плазменных разрядов, характеризующих ток плазмы и токи в обмотках полоидального поля, имеющих однозначную взаимосвязь с величиной зазоров между сепаратрисой и первой стенкой токамака с учетом ударных точек сепаратрисы,
- восстановление равновесия плазмы посредством преобразования указанных значений входных экспериментальных сигналов в выходные сигналы, характеризующие зазоры между сепаратрисой и первой стенкой токамака, а также смещение ударных точек сепаратрисы, с использованием метода итераций Пикара и последующим формированием зависимости между указанными экспериментальными входными и полученными выходными сигналами,
- настройку средства отображения входных магнитных сигналов в выходные сигналы, характеризующие положение сепаратрисы, с использованием полученной на предыдущем шаге зависимости, при этом в течение диверторной фазы на вход средства отображения подают входные магнитные сигналы, измеряемые в режиме реального времени вне плазмы, по результатам преобразования которых в значения, характеризующие положение сепаратрисы, судят о положении и форме плазмы в камере токамака в течение диверторной фазы плазменных разрядов токамака. В качестве экспериментальных входных сигналов плазменных разрядов, а также входных магнитных сигналов, измеряемых в режиме реального времени, используют значения магнитных потоков на магнитных петлях, сигналы на магнитных зондах Мирнова и сигналы на поясах Роговского. В качестве средства отображения входных магнитных сигналов в положение сепаратрисы используют статическую матрицу, искусственную нейронную сеть, наблюдатель состояния с применением их в обратной связи системы магнитного управления формой плазмы в реальном времени токамака.
Основная идея предлагаемого способа состоит в том, что сначала по данным достаточно большого количества зарегистрированных значений магнитных сигналов плазменных разрядов восстанавливают равновесие плазмы посредством известных методов восстановления равновесия плазмы типа итераций Пикара. При этом получают для этих разрядов связь потоков через магнитные петли, магнитные зонды Мирнова, тока плазмы и токов в обмотках полоидального поля с зазорами между сепаратрисой и первой стенкой и смещений ударных точек. Полученную взаимосвязь используют для получения реальных значений зазоров между сепаратрисой и первой стенкой и значений смещения ударных точек по измеряемым вне плазмы в реальном времени значениям входных магнитных сигналов.
После этого настраивают выбранное средство отображения, использующее полученную зависимость между входными сигналами и параметрами, характеризующими форму и положение плазмы. Например, наблюдатель состояния, статическая матрица или искусственная нейронная сеть. Такой способ после настройки может быть применен в реальном эксперименте для управления формой плазмы в обратной связи по зазорам и ударным точкам.
Краткое описание чертежей
Изобретение поясняется чертежами, где
на фиг. 1 - расположение катушек полоидального поля и магнитных петель в вертикальном сечении токамака ГЛОБУС-М2; плазменная граница с ударными точками g1, g2 и зазорами g3-g6;
на фиг. 2 - полоидальный поток на магнитной петле токамака Глобус-М2 в разряде №31648 (а), токи в обмотках PF1 и PF2 (б) полоидального поля и ток плазмы (в) в разряде №31648;
на фиг. 3 - определение полоидального магнитного потока в токамаке;
на фиг. 4 - линии равного уровня восстановленного распределения полоидального потока в токамаке ГЛОБУС-М2 и точки измерения: а - полоидального потока на сепаратрисе; б - зазоров между первой стенкой и сепаратрисой;
на фиг. 5 - восстановленное распределение полоидального потока для момента времени
Figure 00000001
мс разряда №31648 токамака Глобус-М2;
на фиг. 6 - восстановленное распределение плотности тороидального тока плазмы для момента времени
Figure 00000002
мс разряда №31648 токамака Глобус-М2;
на фиг. 7 - оценка величины зазоров в разряде №37320 сферического токамака Глобус-М2 с помощью статической матрицы;
на фиг. 8 - сравнение вариаций зазоров δg, полученных с помощью модели LPV (Linear Parameter Varying - линейное изменение параметров), полученной из кода FCDI (сплошная линия) и оценкой вариаций зазоров
Figure 00000003
, полученной с помощью робастного наблюдателя, синтезированного с помощью LMIs (Linear Matrix Inequalities - линейные матричные неравенства) (прерывистая линия). Разряд токамака Глобус-M2 №37263;
на фиг. 9 - оценка величины зазоров в разряде №37255 сферического токамака Глобус-М2 с помощью нейронной сети.
Осуществление изобретения
В настоящее время в мире действует ряд вертикально вытянутых D-образных токамаков, на которых изучается физика высокотемпературной плазмы и проводятся исследования систем управления плазмой1. К таким наиболее прогрессивным токамакам можно отнести DIII-D (США, год пуска 1986, R= 1,66; a=0,67), JET (Англия, 1992, 3,00; 1,25-2,1), TCV (Швейцария, 0.88; 0,25-0,7), ASDEX Upgrade (Германия, 1991 1,65; 0,5-0,8), EAST (Китай, 2006, 1.75; 0,43), KSTAR (Южная Корея, 2008, 1,8; 0,5), Глобус-М2 (Россия, 1999, 0,36, 0,24)
Для управления формой плазмы, т.е. расположением в пространстве камеры токамака границы плазмы, необходимо по магнитным измерениям вне плазмы (поскольку внутри высокотемпературной плазмы измерения невозможны) восстанавливать распределение полоидального магнитного потока и положение сепаратрисы плазмы.
Предлагаемый способ состоит из двух частей.
В первой части для определения формы плазмы способ основан на итерационном решении дифференциального эллиптического уравнения равновесия плазмы Града-Шафранова (Шафранов В.Д. Равновесие плазмы в магнитном поле // Вопросы теории плазмы. 1963. Вып. 2. С. 92-131.) в частных производных методом функций Грина. Проводят накопление данных изменений входных магнитных сигналов определенного токамака. На основе восстановленных распределений плотности тороидального тока плазмы и полоидального потока находят сепаратрису плазмы (зазоры и смещения ударных точек). Дальнейшее описание действий способа представлено на следующем примере, не ограничивающем возможность применения заявляемого способа, а лишь раскрывающем один из вариантов его реализации.
Способ определения формы плазмы рассматривается на примере магнитной диагностики токамака Глобус-М2 по следующим входным сигналам:
- 21 магнитная петля, обеспечивающие измерение магнитного полоидального потока в различных точках пространства вокруг плазмы,
- пояса Роговского, обеспечивающие измерение полного тока плазмы, суммарный ток плазмы и ток по вакуумной камере токамака,
- а также токи в 8-и полоидальных обмотках токамака Глобус-М.
Расположение катушек полоидального поля, магнитных петель и элементов вакуумной камеры токамака Глобус-М2 представлены на фиг. 1. Типичное показание магнитной диагностики - величина полоидального потока на магнитной петле - приведено на фиг. 2а; токи в основных обмотках управления токамака - на фиг. 2б, ток плазмы на фиг. 2в.
Равновесие плазмы в токамаке описывается распределением полоидального магнитного потока
Figure 00000004
в цилиндрической системе координат
Figure 00000005
с учетом аксиальной симметрии плазмы. Полоидальный поток
Figure 00000006
в точке P определяется как отношение величины магнитного потока на радиан горизонтальной окружности с центром на оси токамака, проходящей через точку P и ограничивающей площадь
Figure 00000007
(фиг. 3):
Figure 00000008
,
где
Figure 00000009
вектор индукции магнитного поля,
Figure 00000010
- скалярное произведение индукции и нормали элемента поверхности. Из определения полоидального потока следует, что он должен обращаться в ноль на аксиальной оси и на бесконечности (граничные условия).
Граница плазмы (сепаратриса) в полоидальной плоскости (вертикальной плоскости, проходящей через ось z) может быть найдена как наибольшая замкнутая линия уровня полоидального потока (Ariola M., Pironti A. Magnetic control of tokamak plasmas, Springer-Verlag, 2008). В случае, если сепаратриса касается камеры токамака, конфигурацию плазмы называют лимитерной. Если плазма не касается стенок токамака, то говорят, что плазма находится в диверторной фазе. Диверторная фаза также характеризуется наличием на сепаратрисе X-точки, в которой магнитное поле обращается в ноль. Сепаратриса имеет ветви, которые «входят» в первую стенку токамака, которой может являться первая стенка вакуумной камеры, или диверторные пластины. Вдоль этих ветвей из плазмы выходят примеси и поглощаются диверторными пластинами. Таким образом, для определения формы и положения магнитной оси плазмы достаточно найти распределение полоидального потока.
Из уравнений Максвелла в дифференциальной форме (Матвеев А.Н. Электричество и магнетизм. М.: Высшая школа, 1983. 463 с.):
Figure 00000011
где
Figure 00000012
- вектор плотности тока,
Figure 00000013
- магнитная постоянная, и аксиальной симметрии тороидальной плазмы
Figure 00000014
можно получить дифференциальное уравнение в частных производных второго порядка, связывающее распределение плотности тороидального тока плазмы
Figure 00000015
с распределением создаваемого им полоидального потока
Figure 00000016
(Матвеев А.Н. Электричество и магнетизм. М.: Высшая школа, 1983. 463 с.)
Figure 00000017
.(1)
Используя уравнение (1), находят распределение полоидального потока при заданном распределении плотности тока плазмы, однако правая часть уравнения (1) неизвестна и должна быть найдена из других условий. Краевая задача является обратной, т.к. в задаче известно следствие, т.е. значения потока вне плазмы, а неизвестной является причина: распределение потока и тороидального тока плазмы. Более того, решение задачи не является единственным.
Комбинируя уравнение (1) с условием баланса между силами Ампера магнитного поля и градиента давления в плазме
Figure 00000018
,
можно получить уравнение Града-Шафранова (Шафранов В.Д. Равновесие плазмы в магнитном поле // Вопросы теории плазмы. 1963. Вып. 2. С. 92-131), которому подчиняется равновесие плазмы в токамаке:
Figure 00000019
где плотность тороидального тока плазмы
Figure 00000020
выражена через газокинетическое давление плазмы p и функцию F, пропорциональную полоидальному току плазмы
Figure 00000021
Функции
Figure 00000022
и
Figure 00000023
также неизвестны, но зависят только от полоидального потока. Ввиду этого, составляющие плотности тороидального тока плазмы
Figure 00000024
могут быть аппроксимированы функциями полоидального потока:
Figure 00000025
Для решения задачи восстановления равновесия плазмы используется итерационный метод Пикара. На каждой итерации плотность тока плазмы аппроксимируется линейной комбинацией базисных функций
Figure 00000026
от полоидального потока, полученного на предыдущей итерации. Полоидальный поток в этом случае вычисляется как сумма, определяемая выбранными базисными функциями:
Figure 00000027
где
Figure 00000028
Здесь G - функция Грина для уравнения (1), имеющая физический смысл полоидального потока в точке
Figure 00000029
, создаваемого единичным бесконечно тонким кольцевым током с координатами
Figure 00000030
, распределение плотности тока которого имеет вид дельта функции
Figure 00000031
, и выражаемая через эллиптические интегралы первого K и второго рода E (Ariola M., Pironti A. Magnetic control of tokamak plasmas, Springer-Verlag, 2008):
Figure 00000032
Коэффициенты
Figure 00000033
,
Figure 00000034
и токи в катушках
Figure 00000035
подбираются на каждой итерации через минимизацию суммы квадратов разностей между рассчитанными и измеренными токами и потоками
Figure 00000036
Здесь
Figure 00000037
и
Figure 00000038
- рассчитанный и измеренный полные токи плазмы,
Figure 00000039
- измеренные токи в обмотках токамака,
Figure 00000040
,
Figure 00000041
,
Figure 00000042
- погрешности измерений полоидального потока, тока плазмы и тока в обмотках соответственно.
На фиг. 4 показаны примеры равновесий плазмы в виде линий равного уровня полоидального потока системы вложенных магнитных поверхностей, полученных применением рассмотренного способа к экспериментальным данным токамака Глобус-М2. В эксперименте на Глобусе-М2 наблюдались различные виды равновесий: с нижней Х-точкой (фиг. 4а) и верхней Х-точкой (фиг. 4б). На этих же рисунках приведены точки измерений как полоидальных потоков на сепаратрисе (фиг. 4а), так и зазоров между сепаратрисой и первой стенкой (фиг. 4б), которые необходимы для управления формой плазмы.
Таким образом, в результате применения рассмотренного итерационного способа восстановления равновесия плазмы к достаточно большому количеству параметров, характеризующих осуществленные плазменные разряды, т.е. в результате определения формы плазмы по экспериментальной базе данных, получают взаимосвязь входных магнитных сигналов и величин, характеризующих форму плазмы в камере токамака.
На следующем шаге заявляемого способа полученную взаимосвязь применяют для настройки одного из известных средств отображения, обеспечивающих преобразование реальных входных сигналов (сигналы на магнитных петлях, на магнитных зондах Мирнова, на поясах Роговского для токов в обмотках полоидального поля) в выходные сигналы (зазоры между сепаратрисой и первой стенкой, смещения ударных точек). В качестве одного из таких средств отображения может использоваться статическая матрица с ее соответствующей настройкой с учетом вышеуказанной полученной взаимосвязи. Принцип статической матричной оценки формы плазмы заключается в том, что в любой момент времени j
Figure 00000043
[1; s i] любого разряда i
Figure 00000043
[1, N], оценка формы плазмы
Figure 00000044
может быть получена путем умножения измеряемых сигналов в конкретный момент времени
Figure 00000045
и матрицы
Figure 00000046
, суммирующимися со значениями базовых зазоров
Figure 00000047
. При этом матрица K и вектор
Figure 00000048
являются постоянными для всех плазменных разрядов:
Figure 00000049
Коэффициенты матрицы K и вектора
Figure 00000048
настраиваются посредством метода наименьших квадратов по сумме квадратов разностей между выходными сигналами по зазорам и оцененными сигналами посредством итераций Пикара (Y.V. Mitrishkin, P.S. Korenev, A.E. Konkov, V.I. Kruzhkov, N.E. Ovsyannikov. New identification approach and methods for plasma equilibrium reconstruction in D-shaped tokamaks. Mathematics. 2022, 10, 40, p. 1-23. https://doi.org/10.3390/math10010040).
Правомерность использования в качестве средства отображения статической матрицы показана на стенде реального времени (Y.V. Mitrishkin. Plasma magnetic control systems in D-shaped tokamaks and imitation digital computer platform in real time for controlling plasma current and shape. Advances in Systems Science and Applications. 2022, Vol 22, No 1, pp. 1-15) на достаточно большом массиве матриц плазменных разрядов токамака Глобус-М2 (фиг. 7).
По аналогии со статической матрицей может применяться искусственная нейронная сеть (Фиг. 9) (Y.V. Mitrishkin, P.S. Korenev, A.E. Konkov, V.I. Kruzhkov, N.E. Ovsyannikov. New identification approach and methods for plasma equilibrium reconstruction in D-shaped tokamaks. Mathematics. 2022, 10, 40, p. 1-23. https://doi.org/10.3390/math10010040, Митришкин Ю.В., Кружков В.И. Быстродействующее восстановление формы плазмы в реальном времени в сферическом токамаке Глобус-М2 нейронной сетью по обработанным данным алгоритмом идентификации итерациями Пикара, Ломоносовские Чтения - 2022, Секция Физики, МГУ). Полносвязная нейронная сеть состоит из нескольких матриц и нелинейных функций активации. Элементы матриц настраиваются методом градиентного спуска, минимизируя сумму квадратов разностей между зазорами, восстановленными с помощью итераций Пикара и оцененными нейронной сетью. Сравнивая результаты работы настроенной на тестовых сигналах с различными функциями активации, размерами и числом скрытых слоев - выбирается наилучшая структура сети.
Предлагаемый способ определения формы плазмы был проверен посредством и другого подхода, а именно, с использованием наблюдателя состояния (Y.V. Mitrishkin, P.S. Korenev, A.E. Konkov, V.I. Kruzhkov, N.E. Ovsyannikov. New identification approach and methods for plasma equilibrium reconstruction in D-shaped tokamaks. Mathematics. 2022, 10, 40, p. 1-23. https://doi.org/10.3390/math10010040)]. Для этого получалась модель плазмы в токамаке в пространстве состояний в непрерывном времени:
Figure 00000050
Затем эта модель дискретизировалась с периодом T s
Figure 00000051
По полученной модели в пространстве состояний в дискретном времени строился наблюдатель состояния
(https://ru.wikipedia.org/wiki/Наблюдатель_(динамические_системы)), который по входным и выходным сигналам объекта определял в темпе наблюдения состояния объекта
Figure 00000052
в число которых входили и зазоры. Ошибка оценивания состояния
Figure 00000053
определяется по формуле
Figure 00000054
.
На фиг. 8 приведены графики восстановления равновесия плазмы посредством наблюдателя состояния в темпе наблюдения на стенде реального времени.
Использование средств отображения применительно к оценке зазоров по данным восстановленных равновесий приводит к существенному снижению времени выполнения задачи (Task Execution Time - TET): для статической матрицы это - 6 мкс, для наблюдателя состояния - 10 мкс, для способа оценки посредством неподвижных филаментов - 30 мкс. Выигранное время особенно важно при относительно коротких плазменных разрядах, длящихся, например, на токамаке Глобус-М2 (Россия), порядка 200-300 миллисекунд (мс) с диверторной фазой в 20-30 мс, на токамаке Damavand (Иран) 50-60 мс. При таких временных промежутках сокращение времени, в течение которого можно получить оценку положения и формы плазмы с целью дальнейшей корректировки и управления, особенно важно. Таким образом, заявляемое изобретение в результате сокращения времени получения информации о состоянии плазмы в токамаке способствует сокращению времени реакции на изменение состояния плазмы.

Claims (6)

1. Способ быстродействующего определения формы плазмы в камере токамака в течение диверторной фазы плазменных разрядов, включающий
- предоставление данных входных экспериментальных сигналов плазменных разрядов, характеризующих ток плазмы и токи в обмотках полоидального поля, имеющих однозначную взаимосвязь с величиной зазоров между сепаратрисой и первой стенкой токамака с учетом ударных точек сепаратрисы,
- восстановление равновесия плазмы посредством преобразования указанных значений входных экспериментальных сигналов в выходные сигналы, характеризующие зазоры между сепаратрисой и первой стенкой токамака, а также смещение ударных точек сепаратрисы, с использованием метода итераций Пикара и последующим формированием зависимости между указанными экспериментальными входными и полученными выходными сигналами,
- настройку средства отображения входных магнитных сигналов в выходные сигналы, характеризующие положение сепаратрисы, с использованием полученной на предыдущем шаге зависимости, при этом в течение диверторной фазы на вход средства отображения подают входные магнитные сигналы, измеряемые в режиме реального времени вне плазмы, по результатам преобразования которых в значения, характеризующие положение сепаратрисы, судят о положении и форме плазмы в камере токамака в течение диверторной фазы плазменных разрядов токамака .
2. Способ по п.1, отличающийся тем, что в качестве входных сигналов плазменных разрядов используют значения магнитных потоков на магнитных петлях, сигналы на магнитных зондах Мирнова, сигналы на поясах Роговского.
3. Способ по п.1, отличающийся тем, что в качестве средства отображения входных магнитных сигналов в положение сепаратрисы используют статическую матрицу, искусственную нейронную сеть, наблюдатель состояния с применением их в обратной связи системы магнитного управления формой плазмы в реальном времени токамака.
RU2022117882A 2022-06-30 Способ быстродействующего определения формы плазмы в камере токамака в течение диверторной фазы плазменных разрядов RU2787571C1 (ru)

Publications (1)

Publication Number Publication Date
RU2787571C1 true RU2787571C1 (ru) 2023-01-11

Family

ID=

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116167247A (zh) * 2023-04-21 2023-05-26 中国科学院合肥物质科学研究院 一种基于冯·哈格诺夫方法的gs方程数值计算方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012064767A1 (en) * 2010-11-08 2012-05-18 Lawrence Livermore National Security, Llc Inertial confinement fusion chamber
RU2633517C1 (ru) * 2016-07-25 2017-10-13 Российская Федерация, от имени которой выступает Государственная корпорация по атомной энергии "Росатом" Устройство для измерения характеристик спектральных линий плазмы в реакторе-токамаке
RU2702137C1 (ru) * 2018-09-17 2019-10-04 Федеральное государственное бюджетное образовательное учреждение высшего образования "Московский государственный университет имени М.В. Ломоносова" (МГУ) Способ формирования модели магнитного управления формой и током плазмы с обратной связью в токамаке
KR102070544B1 (ko) * 2019-04-17 2020-01-29 주식회사 기가레인 플라즈마 안테나 및 이를 포함하는 플라즈마 처리장치
RU2756460C1 (ru) * 2020-10-28 2021-09-30 Федеральное государственное бюджетное образовательное учреждение высшего образования «Московский государственный университет имени М.В.Ломоносова» (МГУ) Способ определения характеристик поверхностных электромагнитных волн в плазменных образованиях конечной длины

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012064767A1 (en) * 2010-11-08 2012-05-18 Lawrence Livermore National Security, Llc Inertial confinement fusion chamber
RU2633517C1 (ru) * 2016-07-25 2017-10-13 Российская Федерация, от имени которой выступает Государственная корпорация по атомной энергии "Росатом" Устройство для измерения характеристик спектральных линий плазмы в реакторе-токамаке
RU2702137C1 (ru) * 2018-09-17 2019-10-04 Федеральное государственное бюджетное образовательное учреждение высшего образования "Московский государственный университет имени М.В. Ломоносова" (МГУ) Способ формирования модели магнитного управления формой и током плазмы с обратной связью в токамаке
KR102070544B1 (ko) * 2019-04-17 2020-01-29 주식회사 기가레인 플라즈마 안테나 및 이를 포함하는 플라즈마 처리장치
RU2756460C1 (ru) * 2020-10-28 2021-09-30 Федеральное государственное бюджетное образовательное учреждение высшего образования «Московский государственный университет имени М.В.Ломоносова» (МГУ) Способ определения характеристик поверхностных электромагнитных волн в плазменных образованиях конечной длины

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Коренев П.С.и др. Реконструкция равновесного распределения параметров плазмы токамака по внешним магнитным измерениям и построение линейных плазменных моделей. Мехатроника, автоматизация и управление. Том 17, N4, 2016, с. 254-265. DOI: 10.17587/mau.17.254-266. *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116167247A (zh) * 2023-04-21 2023-05-26 中国科学院合肥物质科学研究院 一种基于冯·哈格诺夫方法的gs方程数值计算方法

Similar Documents

Publication Publication Date Title
CN104217074B (zh) 一种基于矩阵指数的电磁暂态隐式降阶仿真方法
Beddek et al. Adaptive method for non-intrusive spectral projection—application on a stochastic eddy current NDT problem
Huang et al. GPU-optimized fast plasma equilibrium reconstruction in fine grids for real-time control and data analysis
Trinh et al. Iterative matrix pencil method for power system modal analysis
Wai et al. Neural net modeling of equilibria in NSTX-U
RU2787571C1 (ru) Способ быстродействующего определения формы плазмы в камере токамака в течение диверторной фазы плазменных разрядов
CN108462181B (zh) 考虑稀疏性的智能配电网潮流雅可比矩阵鲁棒估计方法
Sanpei et al. A data-assimilation based method for equilibrium reconstruction of magnetic fusion plasma and its application to reversed field pinch
Chen et al. Benchmarking gyrokinetic simulations in a toroidal flux-tube
Bao et al. Tests of the real-time vertical growth rate calculation on EAST
Pruski et al. Assessment of Polish Power System angular stability based on analysis of different disturbance waveforms
Blum et al. An inverse problem in plasma physics: The identification of the current density profile in a tokamak
Horton et al. Self‐consistent plasma pressure tensors from the Tsyganenko magnetic field models
Xu et al. Polynomial-chaos-based decentralized dynamic parameter estimation using Langevin MCMC
CN103576120A (zh) 三次谐波分量准同步信息传输校验与自愈算法
Hommen et al. A fast, magnetics-free flux surface estimation and q-profile reconstruction algorithm for feedback control of plasma profiles
CN110034559B (zh) 基于切换系统模型的电力系统融合状态估计方法
Wróblewski Neural network evaluation of tokamak current profiles for real time control
Lazareva et al. Computational experiment for solving the Stefan problem with nonlinear coefficients
Farahani et al. Simulation of open-loop plasma vertical movement response in the Damavand tokamak using closed-loop subspace system identification
Zotov et al. Optimization of electromagnetic diagnostic system of the T-15 tokamak
Dokuka et al. Study of Globus-M tokamak poloidal system and plasma position control
Ambrosino et al. On-line plasma shape identification via magnetic measurements
Mulay et al. On the development of adaptive random differential quadrature method with an error recovery technique and its application in the locally high gradient problems
Morrow-Jones et al. Use of tokamak dynamics models for digital filtering and control