基于紧致差分的Grad-Shafranov方程快速解法

doi:10.16568/j.0254-6086.201801002

核聚变与等离子体物理 ›› 2018, Vol. 2 ›› Issue (1): 8-14.DOI: 10.16568/j.0254-6086.201801002

基于紧致差分的Grad-Shafranov方程快速解法

胡金迪1, 2，肖炳甲1，罗正平1，黄耀1

(1. 中国科学院等离子体物理研究所，合肥230031； 2. 中国科学技术大学，合肥230031)

收稿日期:2016-12-21 修回日期:2017-10-22 出版日期:2018-03-15 发布日期:2018-03-15
通讯作者: 罗正平(1983-)，男，安徽池州人，博士，主要从事等离子体控制研究。
作者简介:胡金迪(1991-)，男，安徽蚌埠人，硕士，主要从事CUDA并行计算研究。
基金资助:
国家磁约束核聚变发展研究专项(2014GB103000)；国家自然科学基金(11575245)

Fast direct solver for Grad-Shafranov equation based on compact scheme

HU Jin-di1, 2, XIAO Bing-jia1, LUO Zheng-ping1, HUANG Yao1

(1. Institute of Plasma Physics, Chinese Academy of Science, Hefei 230031; 2. University of Science and Technology of China, Hefei 230026)

Received:2016-12-21 Revised:2017-10-22 Online:2018-03-15 Published:2018-03-15

摘要/Abstract

摘要：

在等离子体平衡重建迭代计算过程中，需要快速求解Grad-Shafranov方程(G-S方程)。构造了具有四阶精度紧致差分格式的离散方程，采用离散正弦变换技术对其进行快速求解并采用CUDATM实现GPU并行加速，将其应用到EAST等离子体平衡重建PEFIT代码中，实现基于紧致差分格式的快速G-S方程求解。结果表明，在65×65的网格下，给定方程右端项电流分布的前提下，使用GPU求解G-S方程所需时间为大约34μs。

关键词: 泊松方程；Grad-Shafranov方程；Solov&rsquo, ev方程；紧致差分格式；DST技术；CUDA

Abstract:

In the process of iterative computation of plasma equilibrium reconstruction, it is needed to fast solve Grad-Shafranov (G-S) equation. A fourth-order compact scheme was constructed with discrete sine transform (DST) technique to solve G-S equation, and CUDATM was used to realize parallel acceleration. This method will be used in EAST plasma equilibrium reconstruction PEFIT code to fast solve G-S equation based on compact scheme. It turns out that in 65×65 mesh, if the right side current distribution is known, the time needed to solve G-S equation based on GPU is about 34μs.

Key words: Poisson equation, Grad-Shafranov equation, Solov’ev equation, Compact scheme, DST technique, CUDA

中图分类号:

TL61+2.3

胡金迪1, 2，肖炳甲1，罗正平1，黄耀1. 基于紧致差分的Grad-Shafranov方程快速解法[J]. 核聚变与等离子体物理, 2018, 2(1): 8-14.

HU Jin-di, XIAO Bing-jia, LUO Zheng-ping, HUANG Yao. Fast direct solver for Grad-Shafranov equation based on compact scheme[J]. NUCLEAR FUSION AND PLASMA PHYSICS, 2018, 2(1): 8-14.

参考文献

[1] Lao L L, John H S, Stambaugh R D, et al. Reconstruction of current profile parameters and plasma shapes in tokamaks[J]. Nucl. Fusion, 1985, 25(11): 1611-1622.

[2] 田振夫. 求解泊松方程的高精度紧致差分方法[J]. 黄淮学刊(自然科学版), 1998, 14(4): 25-28.

[3] 田振夫. 泊松方程的高精度三次样条差分方法[J]. 西北师范大学学报(自然科学版), 1996, 32(2): 13-17.

[4] Wang Y, Zhang J. Sixth order compact scheme combined with multigrid method and extrapolation technique for 2D Poisson equation[J]. Journal of Computational Physics, 2009, 228(1): 137-146.

[5] Zhang J. Multigrid method and fourth-order compact scheme for 2D Poisson equation with unequal mesh-size discretization[J]. Journal of Computational Physics, 2002, 179(1): 170-179.

[6] 葛永斌, 吴文权, 卢曦. 求解泊松方程的六阶离散多重网格方法[C]. 中国工程热物理学会热机气动热力学学术会议, 2001.

[7] Hockney R W. A fast direct solution of Poisson's equation using Fourier analysis[J]. Journal of the ACM (JACM), 1965, 12(1): 95-113.

[8] 岳小宁, 肖炳甲, 罗正平. 基于CUDA的二维泊松方程快速直接求解[J]. 计算机科学, 2013, 40(10): 21-23.

[9] Haverkort J W. Axisymmetric ideal MHD tokamak equilibria[Z]. homepage.tudelft.nl, 2009.

[10] 吕涛. 区域分解算法[M]. 北京: 科学出版社, 1992.

[11] 罗正平, 肖炳甲, 朱应飞, et al. Online plasma shape reconstruction for EAST tokamak[J]. Plasma Science and Technology, 2010, 12(4): 412-415.

[12] Yue X N, Xiao B J, Luo Z P, et al. Fast equilibrium reconstruction for tokamak discharge control based on GPU[J]. Plasma Phys. and Contr. Fusion, 2013, 55(8): 85016-85024.

[13] Pataki A, Cerfon A J, Freidberg J P, et al. A fast, high-order solver for the Grad-Shafranov equation[J]. J. Comp. Phys., 2013, 243: 28-45.

[14] 张舒. GPU高性能运算之CUDA[M]. 中国水利水电出版社, 2009.

[15] Press W H, Teukolsky S A, Vetterling W T, et al. Numerical recipes in C[M]. Cambridge: Cambridge University Press, 1996.

[16] 岳小宁. 基于GPU的等离子体快速平衡重建研究[D]. 合肥: 中国科学技术大学, 2013.

基于紧致差分的Grad-Shafranov方程快速解法

Fast direct solver for Grad-Shafranov equation based on compact scheme

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 3

编辑推荐

Metrics

[1]	朱宗涛1，巩春志1，汪志健1，田修波1，杨士勤1，Ricky Fu2. 工件靶台复合加速离子注入系统中的靶台偏压效应[J]. 核聚变与等离子体物理, 2012, 32(3): 260-264.
[2]	朱廷儒，彭晓东，王刚，赵小明，唐昌建. 二维高阶CIP数值方法及其应用[J]. 核聚变与等离子体物理, 2012, 32(2): 109-116.
[3]	覃攀1，李祥1, 2，陶旭梅1，冉祎1，余徽1，戴晓雁1，印永祥1. 热平衡等离子体热力学函数理论计算[J]. 核聚变与等离子体物理, 2007, 27(4): 343-347.