Journal article

KSTAR equilibrium operating space and projected stabilization at high normalized beta

Along with an expanded evaluation of the equilibrium operating space of the Korea Superconducting Tokamak Advanced Research, KSTAR, experimental equilibria of the most recent plasma discharges were reconstructed using the EFIT code. In near-circular plasmas created in 2009, equilibria reached a stored energy of 54 kJ with a maximum plasma current of 0.34 MA. Highly shaped plasmas with near double-null configuration in 2010 achieved H-mode with clear edge localized mode (ELM) activity, and transiently reached a stored energy of up to 257 kJ, elongation of 1.96 and normalized beta of 1.3. The plasma current reached 0.7 MA. Projecting active and passive stabilization of global MHD instabilities for operation above the ideal no-wall beta limit using the designed control hardware was also considered. Kinetic modification of the ideal MHD n = 1 stability criterion was computed by the MISK code on KSTAR theoretical equilibria with a plasma current of 2 MA, internal inductance of 0.7 and normalized beta of 4.0 with simple density, temperature and rotation profiles. The steep edge pressure gradient of this equilibrium resulted in the need for significant plasma toroidal rotation to allow thermal particle kinetic resonances to stabilize the resistive wall mode (RWM). The impact of various materials and electrical connections of the passive stabilizing plates on RWM growth rates was analysed, and copper plates reduced the RWM passive growth rate by a factor of 15 compared with stainless steel plates at a normalized beta of 4.4. Computations of active RWM control using the VALEN code showed that the n = 1 mode can be stabilized at normalized beta near the ideal wall limit via control fields produced by the midplane in-vessel control coils (IVCCs) with as low as 0.83kW control power using ideal control system assumptions. The ELM mitigation potential of the IVCC, examined by evaluating the vacuum island overlap created by resonant magnetic perturbations, was analysed using the TRIP3D code. Using a combination of all IVCCs with dominant n = 2 field and upper/lower coils in an even parity configuration, a Chirikov parameter near unity at normalized poloidal flux 0.83, an empirically determined condition for ELM mitigation in DIII-D, was generated in theoretical high-beta equilibria. Chirikov profile optimization was addressed in terms of coil parity and safety factor profile.


Related material