Theoretical Review of M N=1 1 Sawtooth
Total Page:16
File Type:pdf, Size:1020Kb
저작자표시-비영리-변경금지 2.0 대한민국 이용자는 아래의 조건을 따르는 경우에 한하여 자유롭게 l 이 저작물을 복제, 배포, 전송, 전시, 공연 및 방송할 수 있습니다. 다음과 같은 조건을 따라야 합니다: 저작자표시. 귀하는 원저작자를 표시하여야 합니다. 비영리. 귀하는 이 저작물을 영리 목적으로 이용할 수 없습니다. 변경금지. 귀하는 이 저작물을 개작, 변형 또는 가공할 수 없습니다. l 귀하는, 이 저작물의 재이용이나 배포의 경우, 이 저작물에 적용된 이용허락조건 을 명확하게 나타내어야 합니다. l 저작권자로부터 별도의 허가를 받으면 이러한 조건들은 적용되지 않습니다. 저작권법에 따른 이용자의 권리는 위의 내용에 의하여 영향을 받지 않습니다. 이것은 이용허락규약(Legal Code)을 이해하기 쉽게 요약한 것입니다. Disclaimer Master's Thesis Theoretical Review of m/n=1/1 Sawtooth Instability and Comparison with Experimental Observations in KSTAR Soomin Lee Department of Physics Graduate School of UNIST 2017 Theoretical Review of m/n=1/1 Sawtooth Instability and Comparison with Experimental Observations in KSTAR Soomin Lee Department of Physics Graduate School of UNIST Theoretical Review of m/n=1/1 Sawtooth Instability and Comparison with Experimental Observations in KSTAR A thesis/dissertation submitted to the Graduate School of UNIST in partial fulfillment of the requirements for the degree of Master of Science Soomin Lee 06. 09. 2017 Approved by _________________________ Advisor Hyeon-Keo Park Theoretical Review of m/n=1/1 Sawtooth Oscillation and Comparison with Experimental Result in KSTAR Soomin Lee This certifies that the thesis/dissertation of Soomin Lee is approved. 06. 09. 2017 signature ___________________________ Advisor: Hyeon-Keo Park signature ___________________________ Kyujin Kwak: Thesis Committee Member #1 signature ___________________________ Min Sup Hur: Thesis Committee Member #2 signature Abstract One of the main problems in nuclear fusion is plasma instability. This instability causes the loss of energy in tokamak, and results in shorter confinement time. By studying plasma instability, it can be predicted in more detail, and therefore, better controlled. In this thesis, sawtooth instability at the core of the tokamak is discussed. At the core of tokamak, the plasma density and temperature shows periodic temporal behavior. When the plasma is heated, the temperature increases linearly. As the core temperature increases and reaches a critical value, the core density and temperature rapidly drop and start to increase again. The periodic relaxation of the core plasma is observed as the sawtooth-like signals of the density and temperature. In this thesis, two famous models by Kadomtsev and Wesson to explain the “sawtooth instability” are introduced. Kadomtsev’s model explains that the fast crash is due to the magnetic reconnection. This model is called ‘the fast reconnection model’. Due to the internal kink instability, the plasma is concentrated near 푞 = 1 surface, and pressure increases in this region. Then the magnetic reconnection takes place near 푞 = 1 surface, and the core heat and temperature collapses through this region. Wesson’s model explains the sawtooth instability assuming a flat 푞-profile inside the 푞 = 1 surface. Due to the flat 푞-profile, magnetic shear near the core is very small, and that results in the formation of ‘hot crescent, cold bubble’ inside the 푞 = 1 surface. This model is called ‘the quasi-interchange’ model. In order to explain both models theoretically, the change of potential energy, 훿푊 for a finite plasma displacement 훏, is derived to determine stability. When 훿푊 > 0, the plasma equilibrium is determined to be stable whereas when 훿푊 < 0, the plasma equilibrium is determined to be unstable. In 푚/푛 = 1/1 internal kink mode, 훿푊 is composed of a second order term, 훿푊2, and a fourth order term, 훿푊4 . Kadomtsev’s model assumes d휉/훿푟 → 0 except 푞 = 1 surface, and Wesson’s model assumes 푞 ~ 1 inside the 푞 = 1 surface. These models have different approaches for explaining the fast disruption in the core. To validate the theoretical models, experimental observation is conducted by measuring several parameters such as temperature or density. In tokamak, it is impossible to measure the temperature, because there is no material to measure the plasma temperature at that high temperature. To measure the temperature, an indirect method by I measuring radiation intensity emitted by gyrating electrons is used in which the radiation is proportional to electron temperature. The ECEI system is a diagnostic system which can directly visualize the MHD instabilities by measuring electron temperature fluctuation. The KSTAR ECEI system measures 2- D/quasi 3-D electron temperature fluctuation which is proportional to the intensity of radiation emitted by electrons. According to Rayleigh-Jean’s law, the optical depth of the plasma is thick enough, the plasma temperature is proportional to the intensity of radiation. Sawtooth oscillations were observed in this study using ECEI system by measuring electron temperature fluctuations in 3D. The sawtooth period was observed as 휏sawtooth ~ 10 ms and collapse time was observed as 휏푐 ~ 100 휇s. The observed spatial structure at the core of the sawtoothing plasma resembles the 푚/푛 = 1/1 internal kink mode, which concludes that the Kadomtsev’s model gives a more relevant explanation of the sawtooth instability in the experiment than Wesson’s model. Keywords: Sawtooth Oscillation, MHD instability, KSTAR ECEI system II III Contents I. Introduction ---------------------------------------------------------------------------------------------- 1 1.1 Motivation ------------------------------------------------------------------------------------------- 1 1.2 Fusion energy and tokamak ------------------------------------------------------------------------ 2 1.3 Background of this study -------------------------------------------------------------------------- 6 1.4 Thesis Outline --------------------------------------------------------------------------------------- 9 II. Theoretical & Mathematical Development --------------------------------------------------------- 10 2.1 MHD Physics Review in Tokamak-------------------------------------------------------------- 10 2.1.1 MHD equations ------------------------------------------------------------------------- 10 2.1.2 Mathematical Description in Tokamak ----------------------------------------------- 13 2.1.3 The 푧-pinch configuration ------------------------------------------------------------ 15 2.1.4 The Screw Pinch Configuration------------------------------------------------------- 17 2.1.5 The Safety factor, 푞 ------------------------------------------------------------------- 19 2.1.6 The Beta, 훽 ----------------------------------------------------------------------------- 22 2.1.7 Properties of large aspect ratio in tokamak ------------------------------------------ 24 2.2 Reviews of MHD instability in tokamak ------------------------------------------------------- 27 2.2.1 Basic idea of MHD stability ----------------------------------------------------------- 27 2.2.2 Energy principle ------------------------------------------------------------------------ 28 2.3 Theories of sawtooth instability ----------------------------------------------------------------- 39 2.3.1 Internal kink, 푚/푛 = 1/1 mode ---------------------------------------------------- 39 2.3.2 Kadomtsev’s model -------------------------------------------------------------------- 41 2.3.3 Wesson’s model ------------------------------------------------------------------------- 46 2.3.4 Comparison of the two models -------------------------------------------------------- 49 IV III. Experimental Methods -------------------------------------------------------------------------------- 50 3.1 KSTAR ECEI System ---------------------------------------------------------------------------- 50 3.1.1 Physical principle of ECE ------------------------------------------------------------- 50 3.1.2 KSTAR ECEI System ------------------------------------------------------------------ 52 3.2 Experimental Set-up ------------------------------------------------------------------------------ 53 IV. Results & Discussions --------------------------------------------------------------------------------- 55 4.1 Experimental results ------------------------------------------------------------------------------- 55 4.1.1 Observation of periodic behavior ----------------------------------------------------- 55 4.1.2 Observation of precursor phase ------------------------------------------------------- 59 4.1.3 Observation of fast crash --------------------------------------------------------------- 62 4.2 Discussions ----------------------------------------------------------------------------------------- 63 V. Conclusion ---------------------------------------------------------------------------------------------- 65 Appendix A. Plasma Parameters --------------------------------------------------------------------------- 66 Appendix B. Vector Relations ------------------------------------------------------------------------------ 67 Appendix B1. Vector formulas ----------------------------------------------------------------------- 67 Appendix B2. Applications to the cylindrical coordinates (푟, 휃, 푧) ---------------------------- 68 Appendix C. Detail derivations ---------------------------------------------------------------------------- 69 Appendix C1. Integration of trigonometric function in section 2.1.4 ---------------------------- 69 Appendix C2. Kadomtsev’s collapse time calculation process ----------------------------------- 70 V References ---------------------------------------------------------------------------------------------------- 73 Acknowledgement ------------------------------------------------------------------------------------------ 77 VI List of figures