Highest fusion performance without harmful edge energy bursts in tokamak.

The path of tokamak fusion and International thermonuclear experimental reactor (ITER) is maintaining high-performance plasma to produce sufficient fusion power. This effort is hindered by the transient energy burst arising from the instabilities at the boundary of plasmas. Conventional 3D magnetic perturbations used to suppress these instabilities often degrade fusion performance and increase the risk of other instabilities. This study presents an innovative 3D field optimization approach that leverages machine learning and real-time adaptability to overcome these challenges. Implemented in the DIII-D and KSTAR tokamaks, this method has consistently achieved reactor-relevant core confinement and the highest fusion performance without triggering damaging bursts. This is enabled by advances in the physics understanding of self-organized transport in the plasma edge and machine learning techniques to optimize the 3D field spectrum. The success of automated, real-time adaptive control of such complex systems paves the way for maximizing fusion efficiency in ITER and beyond while minimizing damage to device components.

Improved Particle Confinement with Resonant Magnetic Perturbations in DIII-D Tokamak H-Mode Plasmas.

Experiments on the DIII-D tokamak have identified a novel regime in which applied resonant magnetic perturbations (RMPs) increase the particle confinement and overall performance. This Letter details a robust range of counter-current rotation over which RMPs cause this density pump-in effect for high confinement (H mode) plasmas. The pump in is shown to be caused by a reduction of the turbulent transport and to be correlated with a change in the sign of the induced neoclassical transport. This novel reversal of the RMP induced transport has the potential to significantly improve reactor relevant, three-dimensional magnetic confinement scenarios.

AC compensation of 3D magnetic diagnostic signals in DIII-D and National Spherical Torus Experiment-Upgrade (NSTX-U) for real-time application.

A time domain algorithm has been developed to remove the vacuum pickup generated by both coil current (DC) and induced vessel current (AC) in real time from three dimensional (3D) magnetic diagnostic signals in the National Spherical Torus Experiment-Upgrade (NSTX-U) and DIII-D tokamaks. The possibility of detecting 3D plasma perturbations in real time is essential in modern and future tokamaks to avoid and control MHD instabilities. The presence of vacuum field pickup, due to toroidally asymmetric (3D) coils or to misalignment between sensors and axisymmetric (2D) coils, pollutes the measured plasma 3D field, making the detection of the magnetic field produced by the plasma challenging. Although the DC coupling between coils and sensors can be easily calculated and removed, the AC part is more difficult. An algorithm based on a layered low-pass filter approach for the AC compensation and its application for DIII-D and NSTX-U data is presented, showing that this method reduces the vacuum pickup to the noise level. Comparison of plasma response measurements with and without vacuum compensation shows that accurate mode locking detection and plasma response identification require precise AC and DC compensations.

Optimizing the differential connection schemes for detecting 3D magnetic perturbations in DIII-D.

Optimizing the differential pair connection scheme (i.e., the set of pairs) of a toroidal array of magnetic sensors dedicated to measuring slowly rotating asymmetric fields can enhance the mode number detection capability and failure-resilience. In this work, the condition number obtained from singular value decomposition of the design matrix is used as a metric to evaluate the quality of a connection scheme. A large number of possible pair connections are usually available, so evaluating all of them may require extensive use of computational resources and can be very time-consuming. Alternative methods to reduce the number of pairs evaluated without losing the capabilities of toroidal mode detection are presented in this paper. Three examples of the applications of such analysis for the 3D magnetic diagnostic system of DIII-D are also presented: the addition of two new toroidal arrays with n > 3 detection capabilities, the modification of an existing toroidal array in the low field side of the machine to accommodate the addition of a helicon antenna, and the design of changes in several toroidal arrays in the high field side to accommodate the addition of a lower hybrid current drive antenna on the center post.

Neoclassical Tearing Mode Seeding by Nonlinear Three-Wave Interactions in Tokamaks.

We report the experimental observation of seed magnetic island formation by nonlinear three-wave coupling of magnetic island triplets. In this experiment, disruptive 2,1 islands are seeded by the coupling of 4,3 and 3,2 tearing modes to a central 1,1 sawtooth precursor. Three-wave interactions between these modes are conclusively identified by bispectral analysis, indicating fixed phase relationships in agreement with theory. This new observation of this seeding mechanism has important implications for future reactors that must operate in stable plasma equilibria, free of disruptive 2,1 islands.

Electromagnetic torque measurements in DIII-D using internal/external magnetic field decomposition.

Given spatially resolved measurements of normal and tangential components of the magnetic field just outside the surface of a magnetically confined plasma, the field at the measurement location can be uniquely decomposed into contributions from the plasma and from external sources. This principle allows direct measurement of the electromagnetic torque on the plasma without knowledge of the distribution of the internal and external currents, similar to the more well-known formalism using the Maxwell stress tensor. The internal/external field decomposition also enables a mixed approach that incorporates any explicitly known current distributions (e.g., from non-axisymmetric coils). We discuss the requirements and limitations of such an approach to torque measurements. Experimental measurements of the torque evolution as a rotating tearing mode locks to the wall in the DIII-D tokamak are consistent with a simple model.

Quasisymmetric Optimization of Nonaxisymmetry in Tokamaks.

Predictive 3D optimization reveals a novel approach to modify a nonaxisymmetric magnetic perturbation to be entirely harmless for tokamaks, by essentially restoring quasisymmetry in perturbed particle orbits as much as possible. Such a quasisymmetric magnetic perturbation (QSMP) has been designed and successfully tested in the KSTAR and DIII-D tokamaks, demonstrating no performance degradation despite the large overall amplitudes of nonaxisymmetric fields and strong response otherwise expected in the tested plasmas. The results indicate that a quasisymmetric optimization is a robust path of error field correction across the resonant and nonresonant field spectrum in a tokamak, leveraging the prevailing concept of quasisymmetry for general 3D plasma confinement systems such as stellarators. The optimization becomes, in fact, a simple eigenvalue problem to the so-called torque response matrices if a perturbed equilibrium is calculated consistent with nonaxisymmetric neoclassical transport.

Wide Operational Windows of Edge-Localized Mode Suppression by Resonant Magnetic Perturbations in the DIII-D Tokamak.

Edge-localized mode (ELM) suppression by resonant magnetic perturbations (RMPs) generally occurs over very narrow ranges of the plasma current (or magnetic safety factor q_{95}) in the DIII-D tokamak. However, wide q_{95} ranges of ELM suppression are needed for the safety and operational flexibility of ITER and future reactors. In DIII-D ITER similar shape plasmas with n=3 RMPs, the range of q_{95} for ELM suppression is found to increase with decreasing electron density. Nonlinear two-fluid MHD simulations reproduce the observed q_{95} windows of ELM suppression and the dependence on plasma density, based on the conditions for resonant field penetration at the top of the pedestal. When the RMP amplitude is close to the threshold for resonant field penetration, only narrow isolated magnetic islands form near the top of the pedestal, leading to narrow q_{95} windows of ELM suppression. However, as the threshold for field penetration decreases with decreasing density, resonant field penetration can take place over a wider range of q_{95}. For sufficiently low density (penetration threshold) multiple magnetic islands form near the top of the pedestal giving rise to continuous q_{95} windows of ELM suppression. The model predicts that wide q_{95} windows of ELM suppression can be achieved at substantially higher pedestal pressure in DIII-D by shifting to higher toroidal mode number (n=4) RMPs.

Spatial and temporal analysis of DIII-D 3D magnetic diagnostic data.

An extensive set of magnetic diagnostics in DIII-D is aimed at measuring non-axisymmetric "3D" features of tokamak plasmas, with typical amplitudes ∼10-3 to 10-5 of the total magnetic field. We describe hardware and software techniques used at DIII-D to condition the individual signals and analysis to estimate the spatial structure from an ensemble of discrete measurements. Applications of the analysis include detection of non-rotating MHD instabilities, plasma control, and validation of MHD stability and 3D equilibrium models.

Observation of a multimode plasma response and its relationship to density pumpout and edge-localized mode suppression.

Density pumpout and edge-localized mode (ELM) suppression by applied n=2 magnetic fields in low-collisionality DIII-D plasmas are shown to be correlated with the magnitude of the plasma response driven on the high-field side (HFS) of the magnetic axis but not the low-field side (LFS) midplane. These distinct responses are a direct measurement of a multimodal magnetic plasma response, with each structure preferentially excited by a different n=2 applied spectrum and preferentially detected on the LFS or HFS. Ideal and resistive magneto-hydrodynamic (MHD) calculations find that the LFS measurement is primarily sensitive to the excitation of stable kink modes, while the HFS measurement is primarily sensitive to resonant currents (whether fully shielding or partially penetrated). The resonant currents are themselves strongly modified by kink excitation, with the optimal applied field pitch for pumpout and ELM suppression significantly differing from equilibrium field alignment.

Pedestal bifurcation and resonant field penetration at the threshold of edge-localized mode suppression in the DIII-D Tokamak.

Rapid bifurcations in the plasma response to slowly varying n=2 magnetic fields are observed as the plasma transitions into and out of edge-localized mode (ELM) suppression. The rapid transition to ELM suppression is characterized by an increase in the toroidal rotation and a reduction in the electron pressure gradient at the top of the pedestal that reduces the perpendicular electron flow there to near zero. These events occur simultaneously with an increase in the inner-wall magnetic response. These observations are consistent with strong resonant field penetration of n=2 fields at the onset of ELM suppression, based on extended MHD simulations using measured plasma profiles. Spontaneous transitions into (and out of) ELM suppression with a static applied n=2 field indicate competing mechanisms of screening and penetration of resonant fields near threshold conditions. Magnetic measurements reveal evidence for the unlocking and rotation of tearinglike structures as the plasma transitions out of ELM suppression.

An upgrade of the magnetic diagnostic system of the DIII-D tokamak for non-axisymmetric measurements.

The DIII-D tokamak magnetic diagnostic system [E. J. Strait, Rev. Sci. Instrum. 77, 023502 (2006)] has been upgraded to significantly expand the measurement of the plasma response to intrinsic and applied non-axisymmetric "3D" fields. The placement and design of 101 additional sensors allow resolution of toroidal mode numbers 1 ≤ n ≤ 3, and poloidal wavelengths smaller than MARS-F, IPEC, and VMEC magnetohydrodynamic model predictions. Small 3D perturbations, relative to the equilibrium field (10(-5) < δB/B0 < 10(-4)), require sub-millimeter fabrication and installation tolerances. This high precision is achieved using electrical discharge machined components, and alignment techniques employing rotary laser levels and a coordinate measurement machine. A 16-bit data acquisition system is used in conjunction with analog signal-processing to recover non-axisymmetric perturbations. Co-located radial and poloidal field measurements allow up to 14.2 cm spatial resolution of poloidal structures (plasma poloidal circumference is ~500 cm). The function of the new system is verified by comparing the rotating tearing mode structure, measured by 14 BP fluctuation sensors, with that measured by the upgraded B(R) saddle loop sensors after the mode locks to the vessel wall. The result is a nearly identical 2/1 helical eigenstructure in both cases.