Since 18 of December 2019 conferences.iaea.org uses Nucleus credentials. Visit our help pages for information on how to Register and Sign-in using Nucleus.

17–22 Oct 2016
Kyoto International Conference Center
Japan timezone

Extension of numerical matching method to weakly nonlinear regime -- beyond the Rutherford theory of magnetic island evolution

18 Oct 2016, 08:30
4h
Kyoto International Conference Center

Kyoto International Conference Center

Takaragaike, Sakyo-ku, Kyoto 606-0001 Japan
Poster THS - Magnetic Confinement Theory and Modelling: Stability Poster 1

Speaker

Prof. Masaru FURUKAWA (Graduate School of Engineering, Tottori University)

Description

The theory of matched asymptotic expansion for resistive MHD is well established for linear modes [1] and for weakly nonlinear evolution [2]. Since then many applications of the Rutherford equation [2] have made much progress in fusion research, especially in the neoclassical tearing mode (NTM) studies [3]. However, the theoretical framework is still based on the Rutherford equation essentially. We have recently developed a new framework for linear stability analysis of resistive MHD, the numerical matching method [4]. This method utilizes a finite-width inner region around a resonant surface, instead of an infinitely thin inner layer. We devised the boundary condition for the direct matching at interfaces between the inner and outer regions. Then we succeeded to remove difficulties that remain in the numerical applications of the traditional matched asymptotic expansion even though sophisticated theories were developed [5]. We developed both an eigenvalue and an initial-value approaches. In this paper, we extend the initial-value approach of the numerical matching method to weakly nonlinear cases. In the presentation, we will explain the theory, and will show numerical results using reduced MHD in cylindrical plasmas that successfully reproduced Rutherford regime of magnetic island evolution. The computational cost is reduced, making inclusion of detailed physical effects easier because such a model requires more efficient solution method. Our new method will certainly aid understanding physics and will substantially contribute to the prediction of MHD activities such as NTMs in fusion plasmas. This work was supported by KAKENHI Grant No. 23760805 and No. 15K06647. [1] H. P. Furth, J. Killeen and M. N. Rosenbluth, Phys. Fluids 6, 459 (1963). [2] P. H. Rutherford, Phys. Fluids 16, 1903 (1973). [3] A. I. Smolyakov, Plasma Phys. Control. Fusion 35, 657 (1993); O. Sauter et al., Phys. Plasmas 4, 1654 (1997); R. J. La Haye, Phys. Plasmas 13, 055501 (2006). [4] M. Furukawa, S.Tokuda and L.-J. Zheng, Phys. Plasmas 17, 052502 (2010); M. Furukawa and S. Tokuda, Phys. Plasmas 18, 062502 (2011); Id., Phys. Plasmas 19, 102511 (2012). [5] A. Pletzer and R. L. Dewar, J. Plasma Phys. 45, 427 (1991); A. Pletzer, A. Bondeson and R. L. Dewar, J. Comput. Phys. 115, 530 (1994); S. Tokuda and T. Watanabe, Phys. Plasmas 6, 3012 (1999); S. Tokuda, Nucl. Fusion 41, 1037 (2001).
Country or International Organization Japan
Paper Number TH/P1-14

Primary author

Prof. Masaru FURUKAWA (Graduate School of Engineering, Tottori University)

Co-author

Dr Shinji TOKUDA (Research Organization for Information Science and Technology)

Presentation materials