Modelling sawteeth in tokamak plasma as a sequence of Multi-region Relaxed MHD equilibria
ORAL
Abstract
The Multiregion Relaxed MHD [1] was successful in the construction of equilibria in 3D configurations, bridging the gap between Taylor relaxation, which allows relaxation but only globally, and ideal MHD, which includes no relaxation at all but infinite constraints. In MRxMHD, the plasma is sliced into sub-volumes separated by ideal interfaces, each undergoes relaxation. Stepped Pressure Equilibrium Code (SPEC) [2] was developed to solve MRxMHD equilibria numerically.
A sawtooth [3] is a spontaneous energy relaxation process in tokamak plasma. In the early phase of a sawtooth crash with a relatively small resistivity, the reconnection velocities are not too fast compared with the Alfven speed. One can model the plasma globally as a sequence of equilibria, treating the reconnection layer as an infinitely thin current sheet [4]. Starting with an unstable axisymmetic MRxMHD equilibrium, we found that a lower energy state exists with a helical core. The two interfaces bounding the volume containing the q=1 surface clash together to form a “cusp-type” current sheet [5]. We then construct a sequence of equilibrium by removing the interfaces one by one. The result is compared to a M3D-C1 [6] simulation.
[1] M. Hole, S. Hudson, and R. Dewar, Nucl. Fusion 47, 746 (2007).
[2] S.R. Hudson, R.L. Dewar et al., Phys. Plasmas 19, 112502-1–18, (2012).
[3] S. von Goeler, W. Stodiek, and N. Sauthoff, Phys. Rev. Lett. 33, 1201 (1974).
[4] R.M. Kulsrud, Phys. Plasmas 18, (2011).
[5] D.A. Uzdensky and R.M. Kulsrud, Phys. Plasmas 4, 3960 (1997).
[6] S.C. Jardin, N. Ferraro, J. Breslau, and J. Chen, Comput. Sci. Discov. 5, 014002 (2012).
A sawtooth [3] is a spontaneous energy relaxation process in tokamak plasma. In the early phase of a sawtooth crash with a relatively small resistivity, the reconnection velocities are not too fast compared with the Alfven speed. One can model the plasma globally as a sequence of equilibria, treating the reconnection layer as an infinitely thin current sheet [4]. Starting with an unstable axisymmetic MRxMHD equilibrium, we found that a lower energy state exists with a helical core. The two interfaces bounding the volume containing the q=1 surface clash together to form a “cusp-type” current sheet [5]. We then construct a sequence of equilibrium by removing the interfaces one by one. The result is compared to a M3D-C1 [6] simulation.
[1] M. Hole, S. Hudson, and R. Dewar, Nucl. Fusion 47, 746 (2007).
[2] S.R. Hudson, R.L. Dewar et al., Phys. Plasmas 19, 112502-1–18, (2012).
[3] S. von Goeler, W. Stodiek, and N. Sauthoff, Phys. Rev. Lett. 33, 1201 (1974).
[4] R.M. Kulsrud, Phys. Plasmas 18, (2011).
[5] D.A. Uzdensky and R.M. Kulsrud, Phys. Plasmas 4, 3960 (1997).
[6] S.C. Jardin, N. Ferraro, J. Breslau, and J. Chen, Comput. Sci. Discov. 5, 014002 (2012).
*This research was undertaken with the assistance of resources and services from the National Computational Infrastructure (NCI), which is supported by the Australian Government.This work was supported by a grant from the Simons Foundation/SFARI (560651, AB). This work is partly funded by Australian Research Council project DP170102606.
–
Presenters
-
Zhisong Qu
- Mathematical Sciences Institute, Australian National University
- Mathematical Sciences Institute, The Australian National University
- Australian Natl Univ
- Australian National University