CN106647330B - 一种应用于托卡马克的等离子体位形实时重建系统及方法 - Google Patents

一种应用于托卡马克的等离子体位形实时重建系统及方法 Download PDF

Info

Publication number
CN106647330B
CN106647330B CN201611206935.0A CN201611206935A CN106647330B CN 106647330 B CN106647330 B CN 106647330B CN 201611206935 A CN201611206935 A CN 201611206935A CN 106647330 B CN106647330 B CN 106647330B
Authority
CN
China
Prior art keywords
plasma
real
time
shape
position shape
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201611206935.0A
Other languages
English (en)
Other versions
CN106647330A (zh
Inventor
夏凡
许光俊
毛瑞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwestern Institute of Physics
Original Assignee
Southwestern Institute of Physics
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southwestern Institute of Physics filed Critical Southwestern Institute of Physics
Priority to CN201611206935.0A priority Critical patent/CN106647330B/zh
Publication of CN106647330A publication Critical patent/CN106647330A/zh
Application granted granted Critical
Publication of CN106647330B publication Critical patent/CN106647330B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B17/00Systems involving the use of models or simulators of said systems
    • G05B17/02Systems involving the use of models or simulators of said systems electric

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Plasma Technology (AREA)

Abstract

一种应用于托卡马克的等离子体位形实时重建系统及方法,包括等离子体诊断系统,等离子体诊断系统连接实时位形采集工控机,实时位形采集工控机上拥有对应的数据采集卡和一张反射内存卡,等离子体诊断系统的信号进入到实时位形采集工控机的数据采集卡上,统一时序系统提供统一的时序信号,同样接入到实时采集工控机上,实时采集工控机通过采集卡采集到的诊等离子体诊断信号的数值信息通过其上的反射内存卡,并发送到反射内存交换机,反射内存交换机将这些诊断信号的数值同步到实时位形重建高级工作站上的反射内存卡中,将计算得到的等离子体位形信息再次通过连接在其上的反射内存卡发送回反射内存交换机,位形控制器通过这些信息进行计算。

Description

一种应用于托卡马克的等离子体位形实时重建系统及方法
技术领域
本发明属于实时重建系统领域,具体涉及一种应用于托卡马克的等离子体位形实时重建系统及方法。
背景技术
在受控核聚变实验研究中,用于约束等离子体的磁约束装置主要是托卡马克(或仿星器)装置。由于等离子体的飘移,单纯依赖纵场是不可能实现带电粒子的约束的;因此必须引入一个极向磁场,使总场形成一种螺旋状的结构。这样粒子在等离子中的漂移方向将随着时间的变化而变化,从而总的漂移相互抵消,达到粒子被磁场较好地约束的目的。
因此,托卡马克这种特定的聚变磁约束装置对的等离子体的实时控制提出了比较高的要求,即要实时控制将等离子体约束在真空室内部,使之整个放电过程中不与真空室相接触,同时需要让等离子体保持一定的形状,具备一定的拉长比和三角形形变,以获得更高的等离子体物理参数,因此对于等离子体实时控制而言,必须首先能够实时得到等离子体的位置和形状(以下简称位形)。
但是目前并没有直接测量等离子体形状的诊断手段,即等离子体的形状并不能直接通过测量得到,需要利用其它的诊断数据进行位形重建。在离线情况下,重建的方式可以通过Grad-Shafranove方程,通过计算得到等离子体的磁通剖面分布,从而确定等离子体的最外层磁面,等离子体最外层磁面即为等离子体的形状,确定了等离子体形状,其位置可以直接推导出来,Grad-Shafranov方程本质上是一个泊松方程,只有数值解,因此需要经过大量的迭代才能够得到,每次计算完一个时间片需要耗时大概几百毫秒到1秒的时间,而等离子体实时控制的控制周期是1ms到几毫秒的量级,并且每次计算只能计算一个时间片,无法连续计算。
要想将这个离线计算的方式用于实时控制,就必须搭建起一个实时的等离子体位形重建系统,包括实时操作系统、实时数据采集、实时数据传输和实时位形计算,在整个等离子体放电实验过程中连续稳定地计算,并且每次计算的耗时限制在一个控制周期的几个毫秒量级范围内。
发明内容
本发明的目的在于提供了一种应用于托卡马克的等离子体位形实时重建系统及方法,能够在实时状态下,利用快速解Grad-Shafranov方程重建位形,根据划分的网格大小不同,在1ms之内到几个毫秒的时间范围内得到当前时刻的等离子体的位置和形状,并且连续地计算,将结果不断地传给位形控制器进行反馈控制,直到放电实验结束。
本发明的技术方案如下:一种应用于托卡马克的等离子体位形实时重建系统,包括等离子体诊断系统,等离子体诊断系统连接实时位形采集工控机,实时位形采集工控机上拥有对应的数据采集卡和反射内存卡A,等离子体诊断系统的信号进入到实时位形采集工控机的数据采集卡上,统一时序系统提供统一的时序信号,同样接入到实时采集工控机上,实时采集工控机通过数据采集卡采集到的诊等离子体诊断信号的数值信息通过其上的反射内存卡A,并发送到反射内存交换机,反射内存交换机将这些诊断信号的数值同步到实时位形重建工作站上的反射内存卡中B,实时位形重建工作站的经过实时计算后,将计算得到的等离子体位形信息再次通过连接在其上的反射内存卡B发送回反射内存交换机,反射内存交换机将这些位形信息同步到位形控制器上的反射内存卡C中,位形控制器通过这些信息进行计算。
所述等离子体诊断系统用于获取等离子体的电磁信息。
一种应用于托卡马克的等离子体位形实时的重建方法,包括以下步骤:
1)实时位形重建工作站进行系统初始化;
2)实时位形采集工控机接到统一时序系统发送的表征实验开始的TTL电平,开始实时位形采集;
3)实时位形采集工控机通过对统一时序系统发送的基准时序方波的处理,得到控制周期信息;
4)实时位形采集工控机在每个控制周期时刻开始阶段,通过其上的采集卡实时采集从等离子体诊断系统过来的电压信号,并将之转化为对应的物理量,通过判断等离子体电流数值大小,确定等离子体放电是否结束,如果没有结束,就转向5),如果结束,则转向2);将数据写入反射内存卡A,并发送网络中断信号,通知实时位形重建工作站进行计算;
5)实时位形重建工作站根据上一个控制周期的磁通分布,重新计算等离子体的格林函数;
6)实时位形重建工作站通过读取其上的反射内存卡B收到的诊断数据,结合5)计算的等离子体格林函数,得到新的能够表征等离子体电流分布的参数;
7)实时位形重建工作站通过6)得到的等离子体电流分布参数,得到等离子体磁通剖面;
8)实时位形重建工作站在7)计算得到的磁通剖面,构建出网格上的双三次样条曲面插值的系数;
9)实时位形重建工作站通过8)计算得到的双三次样条曲面插值系数,得到等离子体的最外层磁面对应的磁通值和等离子体磁轴中心的磁通值;
10)实时位形重建工作站通过8)得到的双三次样条曲面插值系数的基础,结合9)计算得到的等离子形状信息,计算得到等离子体位形控制所需要的控制点处的磁通以及等离子体其他的位置信息。
所述步骤1)中,初始化包括配置文件的读取、格林函数文件的读取、内存的分配。
所述步骤3)中,得到频率为1kHz的控制周期信息。
所述步骤6)中,通过奇异值分解SVD的方法,得到新的能够表征等离子体电流分布的参数。
所述步骤7)中,利用并行技术,构建边界条件,然后利用快速bunuman方法解出Grad-shafranov方程,得到等离子体磁通剖面。
所述步骤8)中,实时位形重建工作站通过并行算法,在7)计算得到的磁通剖面,构建出网格上的双三次样条曲面插值的系数。
所述步骤9)中,通过双三次样条曲面插值算法,得到等离子体的X点位置,并结合计算等离子体与limiter交点的方法,得到等离子体的最外层磁面对应的磁通值和等离子体磁轴中心的磁通值。
还包括步骤11),实时位形重建工作站通过8)得到的双三次样条曲面插值系数的基础,结合9)计算得到的等离子形状信息,计算得到等离子体位形控制所需要的控制点处的磁通以及等离子体其他的位置信息。
本发明的显著效果在于:系统的可靠性好,测量精度高,与离线的位形重建结果高度吻合,适用于等离子体实时位形控制。
附图说明
图1是本发明所述的应用于托卡马克的等离子体位形实时重建系统示意图。
图中:1等离子体诊断系统、2实时位形采集工控机、3数据采集卡、4反射内存卡A、5统一时序系统、6反射内存交换机、7反射内存卡B、8实时位形重建工作站、9反射内存卡C、10位形控制器。
具体实施方式
一种应用于托卡马克的等离子体位形实时重建系统,包括等离子体诊断系统1,用于获取等离子体的电磁信息,等离子体诊断系统1连接实时位形采集工控机2,实时位形采集工控机2上拥有对应的数据采集卡3和反射内存卡A4,等离子体诊断系统1的信号进入到实时位形采集工控机2的数据采集卡3上,统一时序系统5提供统一的时序信号,同样接入到实时采集工控机2上,实时采集工控机2通过数据采集卡3采集到的诊等离子体诊断信号的数值信息通过其上的反射内存卡A4,并发送到反射内存交换机6,反射内存交换机6将这些诊断信号的数值同步到实时位形重建工作站8上的反射内存卡中B7,实时位形重建工作站8的经过实时计算后,将计算得到的等离子体位形信息再次通过连接在其上的反射内存卡B7发送回反射内存交换机6,反射内存交换机6将这些位形信息同步到位形控制器上10的反射内存卡C9中,位形控制器10通过这些信息进行计算。
一种应用于托卡马克的等离子体位形实时重建方法,包括以下步骤:
1)实时位形重建工作站8进行系统初始化,初始化包括配置文件的读取、格林函数文件的读取、内存的分配;
2)实时位形采集工控机2接到统一时序系统5发送的表征实验开始的TTL电平,开始实时位形采集;
3)实时位形采集工控机2通过对统一时序系统5发送的基准时序方波的处理,得到频率为1kHz的控制周期信息;
4)实时位形采集工控机2在每个控制周期时刻开始阶段,通过其上的采集卡3实时采集从等离子体诊断系统1过来的电压信号,并将之转化为对应的物理量,通过判断等离子体电流数值大小,确定等离子体放电是否结束,如果没有结束,就转向5),如果结束,则转向2);将数据写入反射内存卡A4,并发送网络中断信号,通知实时位形重建工作站8进行计算;
5)实时位形重建工作站8根据上一个控制周期的磁通分布,重新计算等离子体的格林函数;
6)实时位形重建工作站8通过读取其上的反射内存卡B7收到的诊断数据,结合5)计算的等离子体格林函数,通过奇异值分解SVD的方法,得到新的能够表征等离子体电流分布的参数;
7)实时位形重建工作站8通过6)得到的等离子体电流分布参数,利用并行技术,构建边界条件,然后利用快速bunuman方法解出Grad-shafranov方程,得到等离子体磁通剖面;
8)实时位形重建工作站8通过并行算法,在7)计算得到的磁通剖面,构建出网格上的双三次样条曲面插值的系数;
9)实时位形重建工作站8通过8)计算得到的双三次样条曲面插值系数,通过双三次样条曲面插值算法,得到等离子体的X点位置,并结合计算等离子体与limiter交点的方法,得到等离子体的最外层磁面对应的磁通值和等离子体磁轴中心的磁通值;
10)实时位形重建工作站8通过8)得到的双三次样条曲面插值系数的基础,结合9)计算得到的等离子形状信息,计算得到等离子体位形控制所需要的控制点处的磁通以及等离子体其他的位置信息;
11)实时位形重建工作站8将本次计算结果写入到其上的反射内存卡B7指定位置,并发送网络中断,通知位形控制器10进行控制处理。

Claims (7)

1.一种应用于托卡马克的等离子体位形实时重建方法,其特征在于:包括以下步骤:
1)实时位形重建工作站(8)进行系统初始化;
2)实时位形采集工控机(2)接到统一时序系统(5)发送的表征实验开始的TTL电平,开始实时位形采集;
3)实时位形采集工控机(2)通过对统一时序系统(5)发送的基准时序方波的处理,得到控制周期信息;
4)实时位形采集工控机(2)在每个控制周期时刻开始阶段,通过其上的采集卡(3)实时采集从等离子体诊断系统(1)过来的电压信号,并将之转化为对应的物理量,通过判断等离子体电流数值大小,确定等离子体放电是否结束,如果没有结束,就转向5),如果结束,则转向2);将数据写入反射内存卡A(4),并发送网络中断信号,通知实时位形重建工作站(8)进行计算;
5)实时位形重建工作站(8)根据上一个控制周期的磁通分布,重新计算等离子体的格林函数;
6)实时位形重建工作站(8)通过读取其上的反射内存卡B(7)收到的诊断数据,结合5)计算的等离子体格林函数,得到新的能够表征等离子体电流分布的参数;
7)实时位形重建工作站(8)通过6)得到的等离子体电流分布参数,得到等离子体磁通剖面;
8)实时位形重建工作站(8)在7)计算得到的磁通剖面,构建出网格上的双三次样条曲面插值的系数;
9)实时位形重建工作站(8)通过8)计算得到的双三次样条曲面插值系数,得到等离子体的最外层磁面对应的磁通值和等离子体磁轴中心的磁通值;
10)实时位形重建工作站(8)通过8)得到的双三次样条曲面插值系数的基础,结合9)计算得到的等离子形状信息,计算得到等离子体位形控制所需要的控制点处的磁通以及等离子体位置信息。
2.根据权利要求1所述的一种应用于托卡马克的等离子体位形实时重建方法,其特征在于:所述步骤1)中,初始化包括配置文件的读取、格林函数文件的读取、内存的分配。
3.根据权利要求1所述的一种应用于托卡马克的等离子体位形实时重建方法,其特征在于:所述步骤3)中,得到频率为1kHz的控制周期信息。
4.根据权利要求1所述的一种应用于托卡马克的等离子体位形实时重建方法,其特征在于:所述步骤6)中,通过奇异值分解SVD的方法,得到新的能够表征等离子体电流分布的参数。
5.根据权利要求1所述的一种应用于托卡马克的等离子体位形实时重建方法,其特征在于:所述步骤7)中,利用并行技术,构建边界条件,然后利用快速bunuman方法解出Grad-shafranov方程,得到等离子体磁通剖面。
6.根据权利要求1所述的一种应用于托卡马克的等离子体位形实时重建方法,其特征在于:所述步骤8)中,实时位形重建工作站(8)通过并行算法,在7)计算得到的磁通剖面,构建出网格上的双三次样条曲面插值的系数。
7.根据权利要求1所述的一种应用于托卡马克的等离子体位形实时重建方法,其特征在于:所述步骤9)中,通过双三次样条曲面插值算法,得到等离子体的X点位置,并结合计算等离子体与limiter交点的方法,得到等离子体的最外层磁面对应的磁通值和等离子体磁轴中心的磁通值。
CN201611206935.0A 2016-12-23 2016-12-23 一种应用于托卡马克的等离子体位形实时重建系统及方法 Active CN106647330B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611206935.0A CN106647330B (zh) 2016-12-23 2016-12-23 一种应用于托卡马克的等离子体位形实时重建系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611206935.0A CN106647330B (zh) 2016-12-23 2016-12-23 一种应用于托卡马克的等离子体位形实时重建系统及方法

Publications (2)

Publication Number Publication Date
CN106647330A CN106647330A (zh) 2017-05-10
CN106647330B true CN106647330B (zh) 2019-09-17

Family

ID=58827188

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611206935.0A Active CN106647330B (zh) 2016-12-23 2016-12-23 一种应用于托卡马克的等离子体位形实时重建系统及方法

Country Status (1)

Country Link
CN (1) CN106647330B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109996036B (zh) * 2019-03-18 2020-12-11 合肥工业大学 用于托卡马克的图像采集与高速处理传输方法
CN111935892B (zh) * 2019-05-13 2022-11-22 中科智云科技有限公司 用于测量等离子体状态的方法和设备
CN112366007B (zh) * 2020-11-11 2023-01-24 核工业西南物理研究院 一种用于托卡马克装置的人身安全联锁系统
CN112559075B (zh) * 2020-12-18 2022-02-01 核工业西南物理研究院 一种ioc远程控制系统
US20230268084A1 (en) * 2021-03-15 2023-08-24 Dalian University Of Technology Simulation method for electron temperature evolution caused by east tokamak radiofrequency wave
CN113161020B (zh) * 2021-04-20 2023-01-24 核工业西南物理研究院 一种用于托卡马克装置的多系统联合等离子体控制平台

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1028152A1 (ru) * 1982-01-25 1985-08-30 Ордена Ленина Физико-Технологический Институт Им.А.Ф.Иоффе Голографический способ исследовани фазовых объектов
CN1696940A (zh) * 2004-05-25 2005-11-16 中国科学院等离子体物理研究所 计算机辅助聚变堆概念设计系统及方法
CN1722194A (zh) * 2004-06-23 2006-01-18 三星Sdi株式会社 等离子体显示器件及其图像处理的方法
CN1794516A (zh) * 2005-12-29 2006-06-28 上海交通大学 自重构等离子体天线
CN102841964A (zh) * 2012-08-17 2012-12-26 大连理工大学 一种用于等离子体刻蚀剖面演化的三维计算方法
CN203859921U (zh) * 2014-05-04 2014-10-01 核工业西南物理研究院 一种测量等离子体电流分布的极化ece诊断系统
CN205450466U (zh) * 2015-12-30 2016-08-10 中国科学院西安光学精密机械研究所 基于表面等离子体双稳态的弱光信号重构装置

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1028152A1 (ru) * 1982-01-25 1985-08-30 Ордена Ленина Физико-Технологический Институт Им.А.Ф.Иоффе Голографический способ исследовани фазовых объектов
CN1696940A (zh) * 2004-05-25 2005-11-16 中国科学院等离子体物理研究所 计算机辅助聚变堆概念设计系统及方法
CN1722194A (zh) * 2004-06-23 2006-01-18 三星Sdi株式会社 等离子体显示器件及其图像处理的方法
CN1794516A (zh) * 2005-12-29 2006-06-28 上海交通大学 自重构等离子体天线
CN102841964A (zh) * 2012-08-17 2012-12-26 大连理工大学 一种用于等离子体刻蚀剖面演化的三维计算方法
CN203859921U (zh) * 2014-05-04 2014-10-01 核工业西南物理研究院 一种测量等离子体电流分布的极化ece诊断系统
CN205450466U (zh) * 2015-12-30 2016-08-10 中国科学院西安光学精密机械研究所 基于表面等离子体双稳态的弱光信号重构装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
EAST装置等离子体放电位形快速识别研究;黄勤超等;《物理学报》;20060131;第55卷(第1期);第281-286页
HT-7托卡马克上实时EFIT程序的数据采集实现;刘有高等;《计算机测量与控制》;20090331;第17卷(第3期);第603-604、616页
J-TEXT装置等离子体平衡重建与集成运算;张帆等;《核聚变与等离子体物理》;20160630;第36卷(第2期);第131-135页

Also Published As

Publication number Publication date
CN106647330A (zh) 2017-05-10

Similar Documents

Publication Publication Date Title
CN106647330B (zh) 一种应用于托卡马克的等离子体位形实时重建系统及方法
Liu et al. Data-driven event detection of power systems based on unequal-interval reduction of PMU data and local outlier factor
Desjacques et al. Non-Gaussian Halo Bias Re-examined: Mass-dependent Amplitude<? format?> from the Peak-Background Split and Thresholding
Militello et al. Characterisation of the L-mode scrape off layer in MAST: decay lengths
CN107807860B (zh) 一种基于矩阵分解的电力故障分析方法及系统
CN110739091B (zh) 一种核电厂次临界度直接测量系统及硬件平台和测量方法
CN103198184A (zh) 一种电力系统低频振荡特征类噪声辨识方法
De Santis et al. On the issues related to compliance assessment of ICNIRP 2010 basic restrictions
Bierwage et al. Dynamics of low-n shear Alfvén modes driven by energetic N-NB ions in JT-60U
Tasakos et al. Deep learning-based anomaly detection in nuclear reactor cores
Böckenhoff et al. Reconstruction of magnetic configurations in W7-X using artificial neural networks
Wai et al. Neural net modeling of equilibria in NSTX-U
CN101566662A (zh) 电磁电器测试系统
CN111935892B (zh) 用于测量等离子体状态的方法和设备
Chang et al. Voltage fault diagnosis of a power battery based on wavelet time-frequency diagram
Beg et al. Evolution, structure, and topology of self-generated turbulent reconnection layers
Harrison et al. Characteristics of X-point lobe structures in single-null discharges on MAST
Wei et al. Reconstruction of tokamak plasma safety factor profile using deep learning
Wei et al. MHD mode tracking using high-speed cameras and deep learning
CN109030976A (zh) 一种电力电子设备监测诊断系统和方法
Verdoolaege et al. Modeling fusion data in probabilistic metric spaces: Applications to the identification of confinement regimes and plasma disruptions
CN103809020B (zh) 互联电网低频振荡频率和阻尼估计值联合置信区间的确定方法
Van Milligen et al. Filaments in the edge confinement region of TJ-II
Xingwei et al. Research on cognitive mechanism and brain-computer interface application in visual-auditory crossmodal stimuli
Meyer et al. Isotope effect on blob-statistics in gyrofluid simulations of scrape-off layer turbulence

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant