Search NASA⌕ Search

Engineering topics

Park, Jong-Kyu

Publications and source records attributed to Park, Jong-Kyu.

Analysis of neoclassical tearing mode stabilization experiment by electron cyclotron injection in KSTAR

We report the neoclassical tearing mode (NTM) stabilization experimental results by injection of the electron cyclotron (EC) beam in KSTAR, and its analysis with integrated modelling. In the KSTAR experiment, NTM was intentionally induced and the EC beam was injected by the plasma control system (PCS) to stabilize it. Here, the EC angle was controlled based on the minimum seeking algorithm for the island growth rate. Although the EC angle was not changed significantly, the island width significantly decreased towards the end of the experiment followed by the recovery of plasma performance. To assess the NTM stabilization experiment from a stability perspective, we performed integrated modelling with TRIASSIC utilizing the reconstructed equilibrium with magnetic diagnostics and motional stark effect diagnostics. As a result, we observed that the NTM is stabilized by EC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tailoring tokamak error fields to control plasma instabilities and transport

A tokamak relies on the axisymmetric magnetic fields to confine fusion plasmas and aims to deliver sustainable and clean energy. However, misalignments arise inevitably in the tokamak construction, leading to small asymmetries in the magnetic field known as error fields (EFs). The EFs have been a major concern in the tokamak approaches because small EFs, even less than 0.1%, can drive a plasma disruption. Meanwhile, the EFs in the tokamak can be favorably used for controlling plasma instabilities, such as edge-localized modes (ELMs). Here we show an optimization that tailors the EFs to maintain an edge 3D response for ELM control with a minimized core 3D response to avoid plasma disruption and unnecessary confinement degradation. We design and demonstrate such an edge-localized 3D response in the KSTAR facility, benefiting from its unique flexibility to change many degrees of freedom in the 3D coil space for the various fusion plasma regimes. This favorable control of the tokamak EF represents a notable advance for designing intrinsically 3D tokamaks to optimize stability and confinement for next-step fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Theoretical investigation of the triggering of neoclassical tearing modes by transient resonant magnetic perturbations in NSTX

The extended perturbed equilibrium code asymptotic matching code is used to simulate the triggering of n = 1 neoclassical tearing modes (NTMs) by a pulsed, rotating, n = 1, resonant magnetic perturbation (RMP) in two example NSTX discharges. Although the two discharges are significantly different from one another, the results of the two sets of simulations are quite similar. The critical n = 1 RMP pulse amplitude required to trigger an n = 1 NTM is minimized when the RMP pulse rotation frequency matches the linear natural frequency of an n = 1 tearing mode, resonant within the plasma, that is metastable to an NTM. However, if there is a frequency mismatch, then the seed magnetic island chain driven at the relevant resonant surface is forced to rotate with respect to the RMP, because the RMP pulse amplitude is nowhere near sufficient to lock the island chain to the RMP. This rotation causes the critical RMP pulse amplitude required to trigger an NTM to oscillate as the RMP pulse duration is varied. Furthermore, the critical amplitude is minimized when the RMP pulse duration is such that seed island chain executes a half-integer number of rotations with respect to the pulse. All of the minima have the same value.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stochastic fluctuation and transport of tokamak edge plasmas with the resonant magnetic perturbation field

Here we present that a statistical method known as the complexity–entropy analysis is useful to characterize a state of plasma turbulence and flux in the resonant magnetic perturbation (RMP) edge localized mode (ELM) control experiment. The stochastic pedestal top temperature fluctuation in the RMP ELM suppression phase is distinguished from the chaotic fluctuation in the natural ELM-free phase. It is discussed that the stochastic temperature fluctuation can be originated from the narrow layer of the field penetration on the pedestal top. The forced magnetic island can emit the resonant drift wave of comparable sizes (relatively low-k) in the RMP ELM suppression phase, and it can result in the generation of stochastic higher wavenumber fluctuations coupled to tangled fields around the island. The analysis of the ion saturation current measurement around the major outer striking point on the divertor shows that it also becomes more stochastic as the stronger plasma response to the RMP field is expected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Verification of neoclassical toroidal viscosity induced by energetic particles

The thermal particles contributed neoclassical toroidal viscosity (NTV) have been successfully developed and explored by many impressive works such as the study by Shaing et al. [Phys. Plasmas 10, 1443 (2003)] and Zhu et al. [Phys. Rev. Lett. 96, 225002 (2006)]. In this work, the scope of the NTV study is extended to explore the contribution of energetic particles (EPs) through both theory and experiments. In theory, the existence of the NTV torque due to the precessional drift resonance of trapped EPs is identified based on the equivalence between the NTV torque and the perturbed drift kinetic energy [J. Park, Phys. Plasmas 18, 110702 (2011)]. Toroidal modeling with the Magneto Resistive Spectrum - drift Kinetic code [Y. Liu, Phys. Plasmas 15, 112503 (2008)], based on this equivalence, indicates that trapped EPs can contribute a significant amount of the NTV torque. Meanwhile, this work also focuses on developing the dedicated DIII-D experiments in the presence of the n = 2 external magnetic perturbation to verify the EP induced NTV (EP-NTV) by measuring the change of the NTV torque while varying the angle and the voltage of the neutral beam injection. However, the developed experiments have been unable to create conditions necessary to clearly demonstrate the presence of EP-NTV. The main challenge is separating the resonant and non-resonant momentum transport responses in the plasma. Finally, the experience, gained from this study, can help the further exploration of EP-NTV in the future experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Non-disruptive error field identification based on magnetic island healing

Here a technique to identify intrinsic error fields (EFs) in tokamaks with minimized risk of disruption is demonstrated on the DIII-D tokamak. The method extends the conventional driven magnetic island ‘compass scan’ approach by modifying asynchronous control waveforms to enable prompt healing of the island instability. Healing of the island is achieved by reducing the imposed non-axisymmetric coil current and raising the density (here via gas fueling). The method is also shown to support multiple island threshold measurements per pulse, thus reducing the number of dedicated pulses necessary to conduct an EF identification. Non-linear modeling with the TM1 code reproduces the experimental results and approximately recovers the critical density required for island healing. Island healing is explained in the non-linear modeling by an increase in the viscous coupling between the static island and the nearby flowing plasma, thus healing the island as it accelerates into the plasma frame. Due to both simplicity and risk minimization, this technique is suitable for plasma-based EF identification in the early commissioning stages of future disruption-averse tokamaks such as ITER and SPARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Parametric dependencies of resonant layer responses across linear, two-fluid, drift-MHD regimes

Non-axisymmetric magnetic fields arising in a tokamak either by external or internal perturbations can induce complex non-ideal MHD responses in their resonant surfaces while remaining ideally evolved elsewhere. This layer response can be characterized in a linear regime by a single parameter called the inner-layer Δ, which enables outer-layer matching and the prediction of torque balance to non-linear island regimes. In this work, we follow strictly one of the most comprehensive analytic treatments including two-fluid and drift MHD effects and keep the fidelity of the formulation by incorporating the numerical method based on the Riccati transformation when quantifying the inner-layer Δ. The proposed scheme reproduces not only the predicted responses in essentially all asymptotic regimes but also with continuous transitions as well as improved accuracies. In particular, the Δ variations across the inertial regimes with viscous or semi-collisional effects have been further resolved, in comparison with additional analytic solutions. The results imply greater shielding of the electromagnetic torque at the layer than what would be expected by earlier work when the viscous or semi-collisional effects can compete against the inertial effects, and also due to the intermediate regulation by kinetic Alfvén wave resonances as rotation slows down. These are important features that can alter the non-axisymmetric plasma responses including the field penetration by external fields or island seeding process in rotating tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Parametric dependencies of resonant layer responses across linear, two-fluid, drift-MHD regimes

Non-axisymmetric magnetic fields arising in a tokamak either by external or internal perturbations can induce complex non-ideal MHD responses in their resonant surfaces while remaining ideally evolved elsewhere. This layer response can be characterized in a linear regime by a single parameter called the inner-layer Delta, which enables outer-layer matching and the prediction of torque balance to non-linear island regimes. Here, we follow strictly one of the most comprehensive analytic treatments including two-fluid and drift MHD effects and keep the fidelity of the formulation by incorporating the numerical method based on the Riccati transformation when quantifying the inner-layer Delta. The proposed scheme reproduces not only the predicted responses in essentially all asymptotic regimes but also with continuous transitions as well as improved accuracies. In particular, the Delta variations across the inertial regimes with viscous or semi-collisional effects have been further resolved, in comparison with additional analytic solutions. The results imply greater shielding of the electromagnetic torque at the layer than what would be expected by earlier work when the viscous or semi-collisional effects can compete against the inertial effects, and also due to the intermediate regulation by kinetic Alfven wave resonances as rotation slows down. These are important features that can alter the nonaxisymmetric plasma responses including the field penetration by external fields or island seeding process in rotating tokamak plasmas.

MHD instability↗

NSTX-U theory, modeling and analysis results

Here, the mission of the low aspect ratio spherical tokamak NSTX-U is to advance the physics basis and technical solutions required for optimizing the configuration of next-step steady-state tokamak fusion devices. NSTX-U will ultimately operate at up to 2 MA of plasma current and 1 T toroidal field on axis for 5 s, and has available up to 15 MW of neutral beam injection power at different tangency radii and 6 MW of high harmonic fast wave heating. With these capabilities NSTX-U will develop the physics understanding and control tools to ramp-up and sustain high performance fully non-inductive plasmas with large bootstrap fraction and enhanced confinement enabled via the low aspect ratio, high beta configuration. With its unique capabilities, NSTX-U research also supports ITER and other critical fusion development needs. Super-Alfvénic ions in beam-heated NSTX-U plasmas access energetic particle (EP) parameter space that is relevant for both α-heated conventional and low aspect ratio burning plasmas. NSTX-U can also generate very large target heat fluxes to test conventional and innovative plasma exhaust and plasma facing component solutions. This paper summarizes recent analysis, theory and modelling progress to advance the tokamak physics basis in the areas of macrostability and 3D fields, EP stability and fast ion transport, thermal transport and pedestal structure, boundary and plasma material interaction, RF heating, scenario optimization and real-time control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design and experimental demonstration of feedback adaptive RMP ELM controller toward complete long pulse ELM suppression on KSTAR

Operation of a fusion power plant requires robust edge localized mode (ELM) suppression simultaneously with high plasma performance. In this study, we describe a novel feedback adaptive resonant magnetic perturbation (RMP) ELM controller designed to address this problem by achieving optimized ELM suppression through the advanced application of 3D RMPs. From real-time $D_α$ data, the controller is achievable robust ELM suppression while simultaneously minimizing the applied RMP in order to enhance plasma performance. In real time, the instantaneous ELM-frequency is analyzed with an adaptive feedback algorithm to determine amplitudes and phases of RMP coil currents that will maximize plasma performance while maintaining ELM suppression. When applied through the KSTAR plasma control system in several experiments using $textiti{n}$ = 1 RMPs, robust ELM suppression is achieved and sustained in feedback while reducing the RMP strength to ~65% of its initial value. Minimization of the RMP strength in this manner not only allows for operation of longer discharges due to a decrease in flux consumption but also allows for a strong recovery of up to ~60% of $β_N$ throughout the ELM-free period.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Overview of recent progress in 3D field physics in KSTAR

Various 3D field physics challenges of magnetically confined plasmas arise when the driving source comes from either externally applied non-axisymmetric 3D magnetic perturbations or plasma instabilities inside the plasma. Recently, several key outstanding topics of 3D field physics have been extensively studied in the Korean Superconducting Tokamak Advanced Research (KSTAR), such as edge-localized-mode (ELM) control by resonant magnetic perturbation (RMP), error field (EF) control, 3D field effects on rotation and transport, and RMP-induced alteration of divertor heat flux and detachment. KSTAR has a few physically unique features (i.e., high rotation and long-pulse plasmas with a low intrinsic EF) and machine/diagnostic capabilities (i.e., 3-row in-vessel control coil and state-of-the-art 2D/3D imaging diagnostics), which have been taken advantage of until now to address critical 3D field physics issues relevant to ITER and K-DEMO. Among many remarkable achievements are the robust access to and control of n = 1 RMP ELM suppression, along with a development of its physics basis tools, parameter expansion, optimization, and long-pulse control techniques. Nonetheless, a series of unresolved 3D physics themes, as well as limited coverage of 3D field operating regimes, have also been identified as future works for the 3D field research in KSTAR. In this paper, we provide an overview about the recent progress of KSTAR 3D field physics and present future plans of KSTAR 3D research toward a future fusion reactor.

3D field physics↗

Nonlinear modeling of the scaling law for the m/n = 3/2 error field penetration threshold

The scaling law for the m/n=3/2 error field (EF) penetration threshold is predicted numerically based on nonlinear single-fluid and two-fluid modeling using the TM1 code. The simulated penetration threshold of radial magnetic field b r at the plasma edge is scaled to the electron density n e , temperature T e , viscous time τ μ , toroidal field B t and the natural frequency ω in the form of b r /B t ∝n e αn T e αT τ μ αμ B t αB ω αω by scanning these parameters separately. Here, α n , α T , α μ , α B and α ω are the scaling coefficients on n e , T e , τ μ , B t and ω, respectively. Single-fluid modeling shows that the 3/2 EF threshold scales as b r /B t ∝n e 0.56 T e 0.6 τ μ -0.59 B t -1.15 ω, which is similar with the analytical scaling law in both the Rutherford and visco-resistive regimes. Yet, two-fluid modeling shows that the scaling law differs significantly in particular regarding the dependence on plasma rotation. In detail, the scaling coefficient α n on density decreases from 0.67 to 0.56 and α T on temperature decreases from 0.67 to 0.32, while α μ on viscous time is around -0.45 and α B on toroidal field decreases slightly from -1.15 to -1, when the ratio |ω E /ω *e | between plasma rotation frequency ω E and diamagnetic drift frequency ω *e varies from 0 to 10. Scans of the plasma rotation reveals that the penetration threshold linearly depends on the perpendicular rotation frequency (or natural frequency) ω ⊥e =ω E +ω *e , and there is a minimum in the required field amplitude when ω ⊥e 0. In addition, the enduring mystery of non-zero penetration threshold at zero plasma natural frequency in EF experiments is resolved by two-fluid simulations. We report that the very small island and smooth bifurcation in EF penetration near zero frequency is hard to detect in the experiment, leading to a finite penetration threshold within the capability of the experimental measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗