Engineering topics

Park, J. -K. (ORCID:0000000324198667)

Publications and source records attributed to Park, J. -K. (ORCID:0000000324198667).

10 records

DOE OSTI2024

Quantification of locked mode instability triggered by a change in confinement

Abstract This work presents the first analysis of the disruptive locked mode (LM) triggered by the dynamics of a confinement change. It shows that, under certain conditions, the LM threshold during the transient is significantly lower than expected from steady states. We investigate the sensitivity to a controlled n = 1 error field (EF) activated prior to the L-H transition in the COMPASS tokamak, at q 95 ∼ 3, β N ∼ 1, and using EF coils on the high-field side of the vessel. A threshold for EF penetration subsequent to the L-H transition is identified, which shows no significant trend with density or applied torque, and is an apparent consequence of the reduced intrinsic rotation of the 2/1 mode during this transient phase. This finding challenges the assumption made in theoretical and empirical works that natural mode rotation can be predicted by global plasma parameters and urges against using any parametric EF penetration scaling derived from steady-state experiments to define the EF correction strategy in the entire discharge. Furthermore, even at EFs below the identified penetration threshold, disruptive locking of sawtooth-seeded 2/1 tearing modes is observed after about 30% of L-H transitions without external torque.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Plasma performance and operational space with an RMP-ELM suppressed edge

Abstract The operational space and global performance of plasmas with edge-localized modes (ELMs) suppressed by resonant magnetic perturbations (RMPs) are surveyed by comparing AUG, DIII-D, EAST, and KSTAR stationary operating points. RMP-ELM suppression is achieved over a range of plasma currents, toroidal fields, and RMP toroidal mode numbers. Consistent operational windows in edge safety factor are found across devices, while windows in plasma shaping parameters are distinct. Accessed pedestal parameters reveal a quantitatively similar pedestal-top density limit for RMP-ELM suppression in all devices of just over 3 × 10 19 m −3 . This is surprising given the wide variance of many engineering parameters and edge collisionalities, and poses a challenge to extrapolation of the regime. Wide ranges in input power, confinement time, and stored energy are observed, with the achieved triple product found to scale like the product of current, field, and radius. Observed energy confinement scaling with engineering parameters for RMP-ELM suppressed plasmas are presented and compared with expectations from established H and L-mode scalings, including treatment of uncertainty analysis. Different scaling exponents for individual engineering parameters are found as compared to the established scalings. However, extrapolation to next-step tokamaks ITER and SPARC find overall consistency within uncertainties with the established scalings, finding no obvious performance penalty when extrapolating from the assembled multi-device RMP-ELM suppressed database. Overall this work identifies common physics for RMP-ELM suppression and highlights the need to pursue this no-ELM regime at higher magnetic field and different plasma physical size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Effects of electron viscosity on resonant layer responses to non-axisymmetric magnetic perturbations

The resonant field penetration to magnetic islands is the central MHD mechanism of non-axisymmetric plasma responses in a tokamak such as disruptive locking or favorable ELM stabilization. The resonant field penetration can be induced by any non-ideal processes as manifested in the delicate balance under the generalized Ohm's law. Here, we show that the viscous effects by electrons are not ignorable in the field penetration unlike previous presumption, even if the electron viscosity is as small as the square root of its mass compared to the ions. It is clear that its effects become only bigger if the electron viscosity becomes anomalously large. The work strictly follows the three-field model in the linear regime targeting the prediction of the onset of the field penetration and successfully extending it with electron viscosity and identifying new regimes. The results also indicate that the error field thresholds become more strongly dependent on plasma density than ones predicted in the linear regimes without the electron viscosity, which is consistent with experimental observations and thus a significant implication.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Real time detection of multiple stable MHD eigenmode growth rates towards kink/tearing modes avoidance in DIII-D tokamak plasmas

Abstract Real time detection of time evolving growth rates of multiple stable magnetohydrodynamic (MHD) eigenmodes has been achieved in DIII-D tokamak experiments via multi-mode three-dimensional (3D) active MHD spectroscopy. The measured evolution of the multi-modes’ growth rates is in good accordance with the variation of the plasma β N . Using experimental equilibria, resistive MARS-F simulations found the two least stable modes to have comparable growth rates to those experimentally measured. Real time and offline calculations of the modes’ growth rates show comparable results and indicate that cleaner system input and output signals will improve the accuracy of the real time stability detection. Moreover, the shortest real time updating time window of multi-mode eigenvalues can be about 2 ms in DIII-D experiments. This real time monitoring of stable, macroscopic kink and tearing modes thus provides an effective tool for avoidance of the most common causes of tokamak disruption.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Transition in particle transport under resonant magnetic perturbations in a tokamak

Abstract Nonlinear 3D MHD simulations and validations reveal that the hybrid particle-MHD transport is a key process for driving the pump-out in the presence of Resonant Magnetic Perturbations (RMPs) in the KSTAR tokamak. Particle transport and the resulting density pump-out by RMPs are shown to be composed of not only the classical flow convection near magnetic islands due to polarization but also the neoclassical ion diffusion across perturbed magnetic surfaces. The latter is known as the Neoclassical Toroidal Viscosity (NTV) and is integrated into nonlinear MHD simulations here for the first time, revealing that the two-stage pump-outs observed in KSTAR experiments are reproduced only with such integrated nonlinear MHD and transport evolution. Near-resonant responses, which have received less attention than the resonant response, play distinct roles in the pump-out along with the island formation. In addition, this modeling is used to investigate the pump-outs in double-null-like plasmas and numerically capture the effect of the double-null shape on the pump-outs, which may explain the difficulty of Edge Localized Mode (ELM) suppression access in double-like plasmas. This reveals new aspects of the impact toroidal geometry and mode coupling have on 3D physics and reveals the importance of near-resonant components in suppressing ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Validation of density pump-out by pedestal-foot magnetic island formation prior to ELM suppression in KSTAR and DIII-D tokamaks

Abstract TM1 nonlinear two-fluid simulations reveal key characteristics of density pump-out caused by pedestal-foot island formation when applying resonant magnetic perturbation (RMP). These characteristics include a bifurcation in pump-out with a low RMP coil current threshold, a sensitivity of pump-out magnitude to q 95 , and the magnitude of pump-out scaling as the square root of the RMP coil current. Dedicated experiments are carried out in the KSTAR and DIII-D tokamaks to validate these features and show that: (1) a staircase bifurcation in density pump-out is observed when slowly ramping up the n = 1 RMP current, and the density fluctuations are found to be slightly decreased from the pedestal-foot to the pedestal-top; (2) the magnitude of density pump-out becomes weaker when decreasing q 95 from 5.5 to 4.9, and a partial recovery of density pump-out is observed when q 95 is ramped down to lower than 4.9; and (3) analysis of a DIII-D database of n = 3 RMP edge-localized mode control experiments finds that the magnitude of density pump-out is proportional to the square root of RMP coil current Δ n e / n e ∝ I RMP 0.5 . These experimental observations are consistent with TM1 simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Overview of the COMPASS results *

COMPASS addressed several physical processes that may explain the behaviour of important phenomena. This paper presents results related to the main fields of COMPASS research obtained in the recent two years, including studies of turbulence, L–H transition, plasma material interaction, runaway electron, and disruption physics: Tomographic reconstruction of the edge/SOL turbulence observed by a fast visible camera allowed to visualize turbulent structures without perturbing the plasma. Dependence of the power threshold on the X-point height was studied and related role of radial electric field in the edge/SOL plasma was identified. The effect of high-field-side error fields on the L–H transition was investigated in order to assess the influence of the central solenoid misalignment and the possibility to compensate these error fields by low-field-side coils. Results of fast measurements of electron temperature during ELMs show the ELM peak values at the divertor are around 80% of the initial temperature at the pedestal. Liquid metals were used for the first time as plasma facing material in ELMy H-mode in the tokamak divertor. Good power handling capability was observed for heat fluxes up to 12 MW m –2 and no direct droplet ejection was observed. Partial detachment regime was achieved by impurity seeding in the divertor. The evolution of the heat flux footprint at the outer target was studied. Runaway electrons were studied using new unique systems—impact calorimetry, carbon pellet injection technique, wide variety of magnetic perturbations. Radial feedback control was imposed on the beam. Forces during plasma disruptions were monitored by a number of new diagnostics for vacuum vessel (VV) motion in order to contribute to the scaling laws of sideways disruption forces for ITER. Current flows towards the divertor tiles, incl. possible short-circuiting through PFCs, were investigated during the VDE experiments. The results support ATEC model and improve understanding of disruption loads.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Optimization of 3D controlled ELM-free state with recovered global confinement for KSTAR with n = 1 resonant magnetic field perturbation

Mitigation of deleterious heat flux from edge-localized modes (ELMs) on fusion reactors is often attempted with 3D perturbations of the confining magnetic fields. However, the established technique of resonant magnetic perturbations (RMPs) also degrades plasma performance, complicating implementation on future fusion reactors. Here, we introduce an adaptive real-time control scheme on the KSTAR tokamak as a viable approach to achieve an ELM-free state and simultaneously recover high-confinement (βN ~ 1.91, β p ~ 1.53, and H 98 ~ 0.9), demonstrating successful handling of a volatile complex system through adaptive measures. We show that, by exploiting a salient hysteresis process to adaptively minimize the RMP strength, stable ELM suppression can be achieved while actively encouraging confinement recovery. This is made possible by a self-organized transport response in the plasma edge which reinforces the confinement improvement through a widening of the ion temperature pedestal and promotes control stability, in contrast to the deteriorating effect on performance observed in standard RMP experiments. These results establish the real-time approach as an up-and-coming solution toward an optimized ELM-free state, which is an important step for the operation of ITER and reactor-grade tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Nonlinear two-fluid modeling of plasma response to RMPs for the ELM control in the ITER baseline

Numerical modeling, combining the toroidal ideal MHD code GPEC and the nonlinear two-fluid MHD code TM1, was used for comprehensive studies of the plasma response to resonant magnetic perturbations (RMPs) with toroidal mode number n = 1–5 for controlling edge-localized modes (ELMs) in ITER for the standard operation scenario (15 MA Q = 10). Several issues related to RMP ELM control are investigated, including the optimization of the RMP coils configuration, the evaluation of the magnitude of density pump-out and the q95 windows of ELM suppression. GPEC calculates the magnetic response, which consistently includes the very important edge kink/peeling response to static magnetic perturbations. Furthermore, GPEC two-dimensional scans of the relative coil current phasing among the three rows of internal coils, at fixed coil current amplitude, reveal the optimal phasing for the RMP coil configuration with n = 1–5, respectively. The poloidal half wavelength of resonant mode at the edge of plasma calculated by GPEC indicates that the midplane row coils have the best resonant coupling with the plasma for n = 2, while the upper and lower row coils have the best resonant coupling with the plasma for n = 3. Based on the plasma kinetic equilibrium and the GPEC calculations of the magnetic response, TM1 was used to simulate the conditions for RMP field penetration in the ITER pedestal. TM1 shows magnetic island formation at the foot of ITER pedestal with RMP coil current threshold ranging from 4 kAt to 8 kAt with n = 2 to 4. These magnetic islands at the pedestal-foot lead to density pump-out, the magnitude of which scales as ${I}_{\text{RMP}}^{0.5}$ and ranges from 5% to 20% at the pedestal-top when scanning the coil current from 4 to 60 kAt. The density pump-out is found to be weaker for higher n RMP. The nonlinear TM1 simulations also show field penetration at the pedestal-top, where the threshold of RMP coil current depends on the q 95 . The alignment of the magnetic island and the location of the pedestal-top decreases the height and width of the pedestal to suppress ELMs. Simulations by two-dimensional scans of RMP coil current and q 95 reveal the accessible q 95 windows of ELM suppression for both n = 3 and 4 RMPs. The predicted q 95 windows of ELM suppression are very similar to the ones in currently operating tokamaks and the required RMP coil current for ELM suppression is less than 40–50 kAt, which is well within the designed capability for ITER. In addition, the simulations indicate that wide q 95 windows of ELM suppression may be accessible in ITER by operating with dominant n = 4 (or n = 5) RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Identification of multiple eigenmode growth rates towards real time detection in DIII-D and KSTAR tokamak plasmas

Abstract The successful application of three-dimensional (3D) magnetohydrodynamic (MHD) spectroscopy enables us to identify the multi-mode eigenvalues in DIII-D and KSTAR tokamak experiments with stable plasmas. The temporal evolution of the multi-modes’ stabilities have been detected. The new method is numerically efficient allowing the real time detection of MHD modes’ stabilities during the discharge. The method performs active detection of the plasma stability by utilizing the upper and lower rows of internal non-axisymmetric coils to apply a wide variety of 3D fields. Multi-mode eigenvalues are extracted using subspace system identification of the plasma response measured by 3D-field magnetic sensors distributed at different poloidal locations. The equivalence of this new method with the one introduced by Wang (2019 Nucl. Fusion 59 024001) has been numerically corroborated. The more robust and efficient calculation developed here will enable real time monitoring of the plasma stability based on the extracted eigenvalues of stable modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗