Empirical scaling of the n = 2 error field penetration threshold in tokamaks
Abstract not provided.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Abstract not provided.
Abstract not provided.
Stellarators offer a promising path towards fusion reactors, but their design and construction are complicated by stringent tolerance requirements on highly complex 3D coils. A potential way to simplify the engineering requirements for stellarators is to use simple planar toroidal field coils along with permanent magnet arrays to generate shaping fields. In order to ensure sufficient field accuracy while minimizing engineering complexity and system cost, new techniques are required to correct the field produced by the permanent magnet arrays to within requirements set by plasma physics. This work describes a novel correction method developed for this purpose. This analysis is applied to the design of a quasi-axisymmetric stellarator that employs a combination of permanent magnets and planar toroidal field coils to generate its magnetic field. Analysis techniques and initial results using the method for error correction on a proposed permanent magnet stellarator are shown, and it is demonstrated that the method successfully meets the design requirements of the project.
Abstract not provided.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Commonwealth Fusion Systems (CFS) is deploying SPARC, a net-energy tokamak, by 2025. A notable challenge facing all tokamak approaches to fusion energy production is maintaining a stable plasma and thereby steady energy production. “Disruptions” to plasma operation are often encountered in devices when they operate at high normalized pressure, high normalized density, or high normalized current. SPARC is unique in that it is expected to demonstrate Q>=2 in a plasma with low normalized pressure and density, but with a more modest current limit buffer sufficient to avoid disruptions. Due to both the high normalized current and high magnetic field of the SPARC design, it is expected to decrease its relative resilience to instabilities driven by non-axisymmetry in the externally applied magnetic field, termed ‘error fields’. To prevent error field driven instabilities, a primary line of defense is strict engineering tolerances in magnet fabrication and installation. A secondary line of defense is a purpose-built set of magnetic coils to correct these errors, termed error field correction coils (EFCCs). The overarching goals of this INFUSE project were to evaluate the effectiveness of the EFCCs in the SPARC design and provide guidance as to how to improve this effectiveness.
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.
Bi-2212 superconducting cables are being considered for use in the high field magnets needed for the next generation of particle accelerators. Magnetization in these cables and the decay of that magnetization lead to field error and field-error drift, respectively, which need to be compensated. To study this, a segment of the winding pack was extracted from a racetrack coil made from Bi-2212 Rutherford cable. Using a Hall probe measurement technique, we measured the response of the cable's magnetization and its magnetization decay to changes in the applied magnetic field. The effect of adjustments to the cycling of the magnetic field was studied, intended to simulate the preinjection cycles of an accelerator magnet. Three M vs. H loops were constructed by sweeping the magnetic field applied to the sample from 0 to 2.5 T, then to a preinjection field 'x' (where x = 0, 0.25, 0.75 T), and finally up to 1 T. The applied field was then held at 1 T for 1500 s, and the magnetization decay was measured. Finally, the decay was found to vary from 8% to 14% after 1500 s, depending on the preinjection field cycle.
Following up on earlier work which demonstrated an improved numerical stellarator coil design optimization performance by the use of stochastic optimization (Lobsien et al., Nucl. Fusion, vol. 58 (10), 2018, 106013), it is demonstrated here that significant further improvements can be made – lower field errors and improved robustness – for a Wendelstein 7-X test case. This is done by increasing the sample size and applying fully three-dimensional perturbations, but most importantly, by changing the design sequence in which the optimization targets are applied: optimization for field error is conducted first, with coil shape penalties only added to the objective function at a later step in the design process. A robust, feasible coil configuration with a local maximum field error of 3.66 % and an average field error of 0.95 % is achieved here, as compared to a maximum local field error of 6.08 % and average field error of 1.56 % found in our earlier work. These new results are compared to those found without stochastic optimization using the FOCUS and ONSET suites. The relationship between local minima in the optimization space and coil shape penalties is also discussed.
Resonant error fields with toroidal harmonics n = 1 and n = 2 in otherwise axisymmetric tokamak plasmas can drive disruptive locked modes when their strength exceeds a threshold value. This paper presents the sensitivity of error field threshold in global parameters such as line-averaged density and toroidal magnetic field as well as toroidal harmonics dependence in KSTAR L-mode plasmas. The error field threshold for mode-locking increases along with density up to a roll-over density, which coincides well with the linear ohmic confinement (LOC) to saturated ohmic confinement (SOC) density for the studied KSTAR plasmas. As a result, this non-monotonic density correlation of the error field threshold in LOC–SOC can be explained by a modification of existing theoretical scaling in SOC, while plasma rotation can also play an important role in setting the resonant error field threshold.
Work supported by the grant, "Stability and Control of Burning Tokamak Plasmas with Resistive Walls" in the period Aug 2017 to 2023 can be cast in three general categories: (1) linear tearing regimes and the effects on the response to error fields and resistive walls: (2) nonlinear tearing modes in the presence of error fields, with locking; and (3) energetic particle effects. In the first category, we have completed a study of tearing mode regimes with extended MHD effects, including parallel dynamics and toroidal effects, notably favorable average curvature. This work has been presented at the 2018 Sherwood Conference. This work has been published in the ArXiv and is being prepared for publication in Physics of Plasmas. A systematic study of shaping in a tokamak plasma with resistive tearing layers and a resistive wall has been presented in the PhD dissertation by Dov Rhodes, a student of Andrew Cole at Columbia University. John Finn was on Rhodes's PhD committee, and was a co-author of a paper published in Physics of Plasmas on that work. A detailed study of nonlinear tearing modes in the presence of finite beta, favorable average curvature, and an error field has been nearly completed. In this study it was found that the response to an error field with real frequency a tearing layer agreed with analytic results, and that tearing mode saturation by flattening of the current profile as well as destabilization due to flattening of the pressure profile were important factors in the locking process. This work was presented at the 2017 APS/DPP Conference and at the 2018 Sherwood Conference. The effects of energetic ions on tearing modes was presented at the 2017 APS/DPP Conference and in the PhD dissertation of Michael Halfmoon. John Finn served as a committee member on Halfmoon's PhD committee. Extensions to this work including two fluid effects were performed and presented at the APS/DPP Conference in 2017.
ARC is designed to produce 400 MW of net electricity and prove the commercial feasibility of a fusion power plant. In order to achieve this goal ARC has to operate with optimal core performance in a stationary scenario that minimises wear on the first wall and divertor. This requires avoiding or mitigating magnetohydrodynamic (MHD) instabilities which have the potential to not only degrade the plasma core but also lead to deleterious transient heat loads on plasma facing components. Therefore, this work aims at characterising the MHD stability of the high performance ARC scenario and inform the design of error field correction coils. Firstly, simulations of vertical displacement events show that an in-vessel coil is not needed and instead the poloidal shaping coils can be used to control vertical stability. These simulations also inform the demands on the corresponding coil power supplies. Stability analysis of the ideal kink mode with or without a conducting wall and kinetic effects suggests that the ARC baseline scenario operates deeply in the stable region. Using RDCON, tearing modes at the 𝑚/𝑛 =2/1 and 3/2 surfaces (with poloidal mode number 𝑚, and toroidal mode number 𝑛) are shown to be linearly stable, and including thermal transport effects in the rational surfaces lead to further stabilisation. However, other transient plasma instabilities can seed neoclassical tearing modes (NTMs). The marginally stable width of NTMs in ARC strongly depends on the internal inductance and can fall below 0.1% of the normalised poloidal flux. Furthermore, an empirical cross-machine model of the 𝑛 =1 error field leading to a disruption predicts a critical error field larger than SPARC but smaller than ITER. Three-dimensional coils can be designed with the Generalised Purturbed Equilbium Code based on a simple model that calculates the maximum correctable error field that is limited by the neoclassical toroidal viscosity torque. Broad scans of different coil geometries identify a set of 2 rows of off-midplane coils to be a suitable solution. It is also determined that such a set of three-dimensional coils is capable of correcting 𝑛 =2 error fields to some degree and creating strong enough 𝑛 =2 or 𝑛 =3 edge resonant perturbation fields for the suppression of edge-localised modes at reasonable coil currents. The final design of the first ARC will be further informed by results from SPARC.
We report that the U.S. High-Luminosity LHC Accelerator Upgrade Project (HL-LHC AUP) is developing MQXFA magnets, a series of 4.5 m long 150 mm aperture high-field Nb 3 Sn quadrupole magnets for the HL-LHC upgrade at CERN. Five pre-series magnets, MQXFA03 through MQXFA07, have been developed. During the magnet assembly stage, we perform magnetic measurements on the coil-pack sub-assembly and magnets after loading to track the field quality for two purposes. First, it serves as a quality assurance tool to check if the magnetic field quality is on track to meet the acceptance criteria. Magnetic measurements are used to understand if magnetic shims are needed to compensate low-order field errors and to meet the field quality targets. Second, the measurements during the assembly stage can also help understand the field quality, especially the geometric field errors, for Nb 3 Sn accelerator magnets. Here we summarize the measurement results of the pre-series MQXFA magnets, including the magnetic axis and twist angle. The results will provide useful feedback for the series production of Nb 3 Sn magnets and on the optimization of field quality of accelerator magnets based on the wind-and-react Nb 3 Sn technology.
Nonaxisymmetric magnetic fields such as the intrinsic error field and the externally applied resonant magnetic perturbation (RMP) in a tokamak are known to influence the plasma momentum transport and flow evolution through plasma response, which itself strongly depends on the plasma flow as well. The nonlinear interaction between plasma response and flow has been previously modeled in the conventional error field theory with the “no-slip” condition, which has been recently extended to allow the “free-slip” condition. In this work, we further target this specific process and numerically simulate the nonlinear plasma response and flow evolution in the presence of a single-helicity RMP in a circular-shaped model tokamak configuration, based on the full resistive MHD model in the initial-value code NIMROD. Time evolution of the parallel (to k) flow or “slip frequency” profile and its asymptotic steady state obtained from the NIMROD simulations are compared with both conventional and extended nonlinear response theories. Here, k is the wave vector of the propagating island. Good agreement with the extended theory with free-slip condition has been achieved for the parallel flow profile evolution in response to RMP in all resistive regimes, whereas the difference from the conventional theory with the no-slip condition tends to diminish as the plasma resistivity approaches zero.
Errors in the Epperlein and Haines (1986 Phys. Fluids ) transport coefficients were recently found at low electron magnetizations, with new magnetic transport coefficients proposed simultaneously by two teams (Sadler et al 2021 Phys. Rev. Lett. and Davies et al 2021 Phys. Plasmas ); these two separate sets of updated coefficients are shown in this paper to be in agreement. Here, the importance of these new coefficients in laser-plasmas with either self-generated or applied magnetic fields is demonstrated. When an external magnetic field is applied, the cross-gradient-Nernst term twists the field structure; this twisting is reduced by the new coefficients in the low magnetization regime. For plasmas where only self-generated magnetic fields are present, the new coefficients are found to result in the magnetic field moving with the Righi–Leduc heat-flow, enhancing the impact of MHD. Simulations of Biermann battery magnetic fields around ICF hot-spot perturbations are presented, with cross-gradient-Nernst transport increasing spike penetration.
We present a framework for estimating the probability of locking to an error field in a rotating tokamak plasma. This leverages machine learning methods trained on data from a mode-locking model, including an error field, resistive magnetohydrodynamics modeling of the plasma, a resistive wall, and an external vacuum region, leading to a fifth-order ordinary differential equation (ODE) system. It is an extension of the model without a resistive wall introduced by Akçay et al. [Phys. Plasmas 28, 082106 (2021)]. Tearing mode saturation by a finite island width is also modeled. We vary three pairs of control parameters in our studies: the momentum source plus either the error field, the tearing stability index, or the island saturation term. The order parameters are the time-asymptotic values of the five ODE variables. Normalization of them reduces the system to 2D and facilitates the classification into locked (L) or unlocked (U) states, as illustrated by Akçay et al., [Phys. Plasmas 28, 082106 (2021)]. This classification splits the control space into three regions: L̂, with only L states; Û, with only U states; and a hysteresis (hysteretic) region Ĥ, with both L and U states. In regions L̂ and Û, the cubic equation of torque balance yields one real root. Region Ĥ has three roots, allowing bifurcations between the L and U states. The classification of the ODE solutions into L/U is used to estimate the locking probability, conditional on the pair of the control parameters, using a neural network. We also explore estimating the locking probability for a sparse dataset, using a transfer learning method based on a dense model dataset.
A rotating tokamak plasma can interact resonantly with the external helical magnetic perturbations, also known as error fields. This can lead to locking and then to disruptions. We leverage machine learning (ML) methods to predict the locking events. We use a coupled third-order nonlinear ordinary differential equation model to represent the interaction of the magnetic perturbation and the plasma rotation with the error field. This model is sufficient to describe qualitatively the locking and unlocking bifurcations. Here, we explore using ML algorithms with the simulation data and experimental data, focusing on the methods that can be used with sparse datasets. These methods lead to the possibility of the avoidance of locking in real-time operations. We describe the operational space in terms of two control parameters: the magnitude of the error field and the rotation frequency associated with the momentum source that maintains the plasma rotation. The outcomes are quan- tified by order parameters that completely characterize the state, whether locked or unlocked. We use unsupervised ML methods to classify locked/unlocked states and note the usefulness of a certain normalization of the order parameters. Three supervised ML classifiers are used in suite to estimate the probability of locking in the region of control parameter space with hysteresis, i.e., the set of control parameters for which both locked and unlocked states can exist. The results show that a neural network gives the best estimate of the locking probability. An analogy of the present locking model with the van der Waals equation of state is also provided.