A numerical analysis of nonlinear behavior of m=1 resistive internal kink instability in a Tokamak is presented in a cylindrical approximation. A two-dimensional nonlinear code which solves incompressible resistive MHD equations with helical symmetry has been developed. For the case of small longitudinal wave number k, the magnetic axis moves outward until it reaches the critical surface where the value of helical magnetic flux has the same value as the one at the magnetic axis. On the other hand, as k is increased, the shift of the magnetic axis saturates before touching this critical surface. This result gives a possible explanation of suppression of internal disruptions observed in high power neutral beam injection experiment to Tokamaks.
トカマクにおけるm=1抵抗性内部キンク不安定性の非線形挙動の数値解析を、円柱近似で提示する。ヘリカル対称性を有する非圧縮性抵抗性MHD方程式を解く二次元非線形コードを開発した。小さな縦波数kの場合、磁気軸は、ヘリカル磁束の値が磁気軸上の値と等しくなる臨界表面に到達するまで外側へ移動する。一方、kを増加させると、磁気軸のシフトはこの臨界表面に接触する前に飽和する。この結果は、高電力中性粒子ビーム入射実験において観測される内部ディスラプションの抑制の可能な説明を与える。