Parallel computing of magnetohydrodynamic simulations based on conservative perturbed model

doi:10.16568/j.0254-6086.202101002

Nuclear Fusion and Plasma Physics ›› 2021, Vol. 41 ›› Issue (1): 8-15.DOI: 10.16568/j.0254-6086.202101002

• Plasma Physics • Previous Articles Next Articles

Parallel computing of magnetohydrodynamic simulations based on conservative perturbed model

TAN Yao1, 2, MA Jun1, XIANG Nong 1

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

Received:2019-10-14 Revised:2020-01-14 Online:2021-03-15 Published:2022-01-18

基于微扰守恒模型的磁流体模拟并行计算

谭耀 1, 2，马骏 1，项农 1

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

作者简介:谭耀(1992−)，男，湖南怀化人，硕士研究生，从事等离子体物理模拟研究。
基金资助:
国家自然科学基金(1197050421)

Abstract

Abstract: Based on the conservative perturbed model, the nonlinear magnetohydrodymamics code in
two-dimensional slab geometry is developed. It applies the ﬂux vector splitting method and the mapped weighted
essentially non-oscillatory reconstruction scheme. Numerical tests of the ideal Kelvin-Helmholtz instability and
the resistive single tearing mode show its capability and high accuracy in solving linear and nonlinear MHD
instabilities. To run the code on the computer cluster, a parallel computing architecture based on the MPI and
OpenMPlibraries is developed. In the process of parallel computing, load balancing is used to dynamically
allocate the computing amount of each node, which improves the computing efficiency.

Key words: Parallel computation, Weighted essentially non-oscillatory, Conservative perturbed MHD model;
Numerical simulation

摘要： 基于微扰形式的守恒磁流体方程，采用基于通量向量分裂(FVS)以及加权本质无振荡格式(WENO)的
有限差分方法，发展了平板位形下的非线性磁流体模拟程序。程序准确地模拟了开尔文-亥姆霍兹不稳定性和单
撕裂模不稳定性。采用 MPI 和 OpenMP 的并行计算框架，使得程序可以在多节点上运行。在并行计算过程中采
用负载均衡化，对各个节点计算量进行动态分配，提高了计算效率。

关键词: 并行计算；WENO；微扰守恒磁流体模型；数值模拟

CLC Number:

O534+.2

TAN Yao, , MA Jun, XIANG Nong . Parallel computing of magnetohydrodynamic simulations based on conservative perturbed model[J]. Nuclear Fusion and Plasma Physics, 2021, 41(1): 8-15.

谭耀, 马骏, 项农 . 基于微扰守恒模型的磁流体模拟并行计算 [J]. 核聚变与等离子体物理, 2021, 41(1): 8-15.

References

[1] Sovinec C R, Schnack D D, Pankin A Y, et al. Nonlinear extended magnetohydrodynamics simulation using high-order finite elements [J]. Journal of Physics Conference Series, 2005, 16(1): 355−386.
[2] Breslau J, Ferraro N, Jardin S. Some properties of the M3D-C1 form of the three-dimensional magneto-hydrodynamics equations [J]. Phys. Plasmas, 2009, 16(9):092503.
[3] Park W, Belova E V. Plasma simulation studies using multilevel physics models [J]. Phys. Plasmas, 1999, 6(5):1796−1803.
[4] Brio M, Wu C C. An upwind differencing scheme for the equations of ideal magnetohydrodynamics [J]. Journal of Computational Physics, 1988, 75(2): 400−422.
[5] Ma J, Guo W, Yu Z. Simulating magnetohydrodynamic instabilities with conservative perturbed MHD model using discontinuous Galerkin method [J]. Commun.Comput. Phys., 2017, 21(5): 1429−1448.
[6] Shu C W. Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws [M]. Springer, 1997.
[7] Jiang R L, Fang C, Chen P F. A new MHD code with adaptive mesh refinement and parallelization for astrophysics [J]. Computer Physics Communications,2012, 183(8): 1617−1633.
[8] Wong U H, Aoki T, Wong H C. Efficient magnetohydrodynamic simulations on distributed multi-GPU systems using a novel GPU Direct-MPI hybrid approach [J]. Computer Physics Communications, 2014, 185(7): 1901−1913.
[9] Brackbill J U, Barnes D C. The effect of nonzero ∇⋅B on the numerical solution of the magnetohydrodynamic equations [J]. Journal of Computational Physics, 1980, 35(3): 426−430.
[10] Evans C R, Hawley J F. Simulation of magneto-hydrodynamic flows: a constrained transport method [J].Astrophysical Journal, 2007, 332(2): 659−677.
[11] Dedner A, Kemm F, KröNer D, et al. Hyperbolic divergence cleaning for the MHD equations [J]. Journal of Computational Physics, 2002, 175(2): 645−673.
[12] Mignone A, Tzeferacos P, Bodo G. High-order conservative finite difference GLM-MHD schemes for cell-centered MHD [J]. Journal of Computational Physics, 2010, 229(17): 5896−5920.
[13] Fornberg, Bengt. Generation of finite difference formulas on arbitrarily spaced grids [J]. Mathematics of Computation, 1988, 51(184): 699−706.
[14] Lax P D. Systems of conservation laws [J]. Commun.Pure Appl. Math., 2010, 13(2): 217−237.
[15] Dahlburg R B, Karpen J T. Transition to turbulence in solar surges [J]. Astrophysical Journal, 1994, 434(2):766−772.
[16] Li X H, Zhang J, Yang S H, et al. Observing Kelvin-Helmholtz instability in solar blowout jet [J].Scientific Reports, 2018, 8(1): 8136.
[17] Keppens R, TóTH G, Westermann R H J, et al. Growth and saturation of the Kelvin-Helmholtz instability with parallel and antiparallel magnetic fields [J]. Journal of Plasma Physics, 1999, 61(1): 1−19.
[18] Miura A, Pritchett P L. Nonlocal stability analysis of the MHD Kelvin-Helmholtz instability in a compressible plasma [J]. Journal of Geophysical Research, 1982,87(A9): 7431−7444.
[19] Wesson J. Finite resistivity instabilities of a sheet pinch[J]. Phys. Fluids, 2004, 6(2): 130−134.
[20] Del Zanna L, Landi S, Papini E, et al. The “ideal” tearing mode: theory and resistive MHD simulations [J]. Journal of Physics: Conference Series, 2016, 719(1): 12−16.

Parallel computing of magnetohydrodynamic simulations based on conservative perturbed model

基于微扰守恒模型的磁流体模拟并行计算

PDF

Knowledge

Abstract

Cite this article

share this article

References

Related Articles 5

Recommended Articles

Metrics

[1]	TAN Sheng-jun, ZHANG Yang, YE Min-you, ZHANG Xiao-dong, LIU Lei, LUO Zheng-ping, GONG Xian-zu. Analysis and prediction of disruption induced by vertical instability on EAST [J]. NUCLEAR FUSION AND PLASMA PHYSICS, 2019, 39(2): 104-111.
[2]	CHEN Hui, ZHONG Qing-lin, LIU San-qiu. Dust-acoustic shock waves in strongly coupled dusty plasmas with charge fluctuating positive dust [J]. NUCLEAR FUSION AND PLASMA PHYSICS, 2016, 36(4): 315-321.
[3]	WANG Yi-ru, HU Shuang-hui, YAO Long-bao, CHEN Shu-ping. Discrete Alfvén eigenmodes in high performance discharges in the DⅢ-D tokamak [J]. NUCLEAR FUSION AND PLASMA PHYSICS, 2012, 32(2): 140-147.
[4]	YAO Long-bao, HU Shuang-hui, WANG Yi-ru, CHEN Shu-ping. Discrete Alfvén eigenmodes in present tokamaks [J]. NUCLEAR FUSION AND PLASMA PHYSICS, 2012, 32(1): 8-14.
[5]	DAI Yu-jie, LIU Jin-yuan, WANG Xue-hui. Effect of the plasma pressure on growth rate and eigenfunction of line-tied kink instability [J]. NUCLEAR FUSION AND PLASMA PHYSICS, 2010, 30(4): 306-311.