Affiliation:
1. Sun Yet-sen University
2. Zhejiang University
Abstract
Abstract
Aiming at the long-term and high-precision simulation of the magnetohydrodynamic (MHD) instabilities in the tokamak model, we developed a parallelized solver based on a fully implicit difference scheme. A 4th-order precision difference scheme and the Newton-Krylov method are employed in the proposed solver for both the flow and the electromagnetic field. To achieve high parallel efficiency, we adopt a strategy based on the spatial domain decomposition to partition the large Jacobian matrices in the iteration, and a buffer area based on the grid density is utilized to minimize the memory and time consumption. The accuracy of the methodology is verified, and the numerical results are validated by comparison with recognized results. The numerical results of the tearing mode instability in the tokamak model have demonstrated the precision and reliability of the algorithm, and the high parallel efficiency has been proven by the scalability test on the platform with up to 1280 threads, showing significant potential in the large-scale simulation of MHD problems.
Publisher
Research Square Platform LLC
Reference35 articles.
1. Chapter 3: MHD stability, operational limits and disruptions,";Mhd IPEGoD;Nuclear Fusion,1999
2. Higher stable beta by use of pre-emptive electron cyclotron current drive on DIII-D,";Haye RJL;Nuclear Fusion,2005
3. Stabilization effect of early ECCD on a neoclassical tearing mode in the JT-60U tokamak,";Nagasaki K;Nuclear Fusion,2003
4. Numerical study of magnetic island suppression by RF waves in large tokamaks,";Yoshioka Y;Nuclear Fusion,1984
5. Neoclassical tearing modes and their control,";Haye RJ;Physics of Plasmas,2006