中国物理B ›› 2017, Vol. 26 ›› Issue (8): 85204-085204.doi: 10.1088/1674-1056/26/8/085204

• PHYSICS OF GASES, PLASMAS, AND ELECTRIC DISCHARGES • 上一篇    下一篇

Fast parallel Grad-Shafranov solver for real-time equilibrium reconstruction in EAST tokamak using graphic processing unit

Yao Huang(黄耀), Bing-Jia Xiao(肖炳甲), Zheng-Ping Luo(罗正平)   

  1. 1 Institute of Plasma Physics, Chinese Academy of Sciences, Hefei 230031, China;
    2 School of Nuclear Science & Technology, University of Science & Technology of China, Hefei 230027, China
  • 收稿日期:2017-02-22 修回日期:2017-04-01 出版日期:2017-08-05 发布日期:2017-08-05
  • 通讯作者: Zheng-Ping Luo E-mail:zhpluo@ipp.ac.cn
  • 基金资助:

    Project supported by the National Magnetic Confinement Fusion Research Program of China (Grant No. 2014GB103000), the National Natural Science Foundation of China (Grant No. 11575245), and the National Natural Science Foundation of China for Youth (Grant No. 11205191).

Fast parallel Grad-Shafranov solver for real-time equilibrium reconstruction in EAST tokamak using graphic processing unit

Yao Huang(黄耀)1, Bing-Jia Xiao(肖炳甲)1,2, Zheng-Ping Luo(罗正平)1   

  1. 1 Institute of Plasma Physics, Chinese Academy of Sciences, Hefei 230031, China;
    2 School of Nuclear Science & Technology, University of Science & Technology of China, Hefei 230027, China
  • Received:2017-02-22 Revised:2017-04-01 Online:2017-08-05 Published:2017-08-05
  • Contact: Zheng-Ping Luo E-mail:zhpluo@ipp.ac.cn
  • About author:0.1088/1674-1056/26/8/
  • Supported by:

    Project supported by the National Magnetic Confinement Fusion Research Program of China (Grant No. 2014GB103000), the National Natural Science Foundation of China (Grant No. 11575245), and the National Natural Science Foundation of China for Youth (Grant No. 11205191).

摘要:

To achieve real-time control of tokamak plasmas, the equilibrium reconstruction has to be completed sufficiently quickly. For the case of an EAST tokamak experiment, real-time equilibrium reconstruction is generally required to provide results within 1ms. A graphic processing unit (GPU) parallel Grad-Shafranov (G-S) solver is developed in P-EFIT code, which is built with the CUDATM architecture to take advantage of massively parallel GPU cores and significantly accelerate the computation. Optimization and implementation of numerical algorithms for a block tri-diagonal linear system are presented. The solver can complete a calculation within 16 μs with 65×65 grid size and 27 μs with 129×129 grid size, and this solver supports that P-EFIT can fulfill the time feasibility for real-time plasma control with both grid sizes.

关键词: tokamak, Grad-Shafranov equation, equilibrium reconstruction, GPU parallel computation

Abstract:

To achieve real-time control of tokamak plasmas, the equilibrium reconstruction has to be completed sufficiently quickly. For the case of an EAST tokamak experiment, real-time equilibrium reconstruction is generally required to provide results within 1ms. A graphic processing unit (GPU) parallel Grad-Shafranov (G-S) solver is developed in P-EFIT code, which is built with the CUDATM architecture to take advantage of massively parallel GPU cores and significantly accelerate the computation. Optimization and implementation of numerical algorithms for a block tri-diagonal linear system are presented. The solver can complete a calculation within 16 μs with 65×65 grid size and 27 μs with 129×129 grid size, and this solver supports that P-EFIT can fulfill the time feasibility for real-time plasma control with both grid sizes.

Key words: tokamak, Grad-Shafranov equation, equilibrium reconstruction, GPU parallel computation

中图分类号:  (Tokamaks, spherical tokamaks)

  • 52.55.Fa
52.55.-s (Magnetic confinement and equilibrium) 52.65.Kj (Magnetohydrodynamic and fluid equation)