Search NASA⌕ Search

SEARCH · Search NASA

Results for “fusion impurity control”

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.

At least 73 records · Page 4

Stellarator Theory and Non-Axisymmetric Shaping (Final Technical Report)

Divertors are a critical concept for all magnetically confined fusion systems. The W7-X stellarator has an island divertor as do many of the privately funded stellarator reactor designs, but non-resonant divertors have many advantages, which are defined in the paper OSTI ID: 3021689, “Stellarators with enhanced tritium confinement and edge radiation control” and discussed in OSTI ID: 3021687, 3021688, 3021692, and 3021695. The three basic advantages are: (1) Resilience, the location at which the diverted plasma reaches the walls is insensitive to plasma conditions, unlike island divertors. (2) The width of the region at the plasma edge that flows into the divertor has a controllable width. It can be made sufficiently broad to avoid neutrals impending on the plasma from charge exchanging with high energy plasma particles. The resulting high energy neutrals erode the walls. (3) The controllable width and confinement of non-resonant divertors allow high-Z impurities to be added. The shortness of their confinement time limits their diffusion into the plasma interior. High-Z impurity radiation at the plasma edge seems the only way to sufficiently spread the outcoming power over the walls. The lack of experiments is a major impediment to application of non-resonant divertors, although the small experiment STAR_LITE is being built at Hampton University to study them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Time dependent model for non-inductive ECH X-I current ramp-up for SHPD tokamak facility

Non-inductive (NI) plasma current start-up and ramp-up is an important research topic for spherical tokamak (ST) based reactors and fusion pilot plant (FPP). For a compact FPP, the OH flux availability is highly restricted due to its compact geometry. Efficient fundamental extraordinary mode (X-I) electron cyclotron heating (ECH) current start-up and ramp-up regime was identified for a reactor-like high toroidal magnetic field range which has more than a hundred times higher current drive efficiency compared to more conventional ECH methods for the relevant start-up temperature range. High current drive efficiency is possible due to the strong X-I fundamental ECH interaction only with unidirectional passing electrons constrained by the wave accessibility conditions. Here, we extend the X-I electron cyclotron current drive (ECCD) investigation to a time dependent model to simulate the non-inductive current ramp-up to 10 MA for the Sustained High-Power Density (SHPD) facility. As the X-I ECCD driven current I EC rises, due to the back EMF driven negative current, the net plasma current Ip rises more slowly with the current resistive time scale. For tokamak confinement time (both L-mode and H-Mode) which tends to rise with I p , a positive feedback results and even with constant applied ECH power, T e0 , I EC , and I p can continue to rise to very high values. However, in a realistic situation, the T e0 rise should saturate due to a number of factors such as enhanced core radiation and increased power loss at high temperature. To simulate this effect, we adopt a maximum T e0 model which would limit the temperature rise to certain T e0 . With this model, we investigated the current ramp-up for various ECH power levels and the maximum T e0 of 15, 20, and 25 keV. We find that while the power required is reduced with increasing T e0 limit due to increased current drive efficiency, the time to reach 10 MA tends to go up due to the reduced plasma resistivity for the higher T e0 limit. We also find for a given Te0 limit, the time to reach 10 MA tends to be reduced by increasing the applied ECH power by over driving the current ramp-up where I EC is driven at significantly higher level than 10 MA. While the current ramp-up time may not be an issue for the steady-state reactor systems, if it is desirable to minimize the current ramp-up time, it is prudent to have a sufficient ECH power for current over-drive and have some T e0 limiting tools such as impurity seeding for enhanced radiation. In conclusion, a well-controlled NI ECH start-up also has a potential of improving the tokamak start-up reliability and avoid run-away electrons while the NI off-axis current drive could enhance MHD stability and plasma performance improvements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Scenario optimization for the tokamak ramp-down phase in RAPTOR: Part B. safe termination of DEMO plasmas

An optimized plasma current ramp-down strategy is critical for safe and fast termination of plasma discharges in a tokamak demonstration fusion reactor (DEMO), both in planned and emergency scenarios, avoiding plasma disruptions and excessive heat loads to the first wall. Plasma stability limits and machine-specific technical requirements constrain the stable envelope through which the plasma must be navigated. Large amounts of auxiliary heating are required throughout the ramp-down phase, to avoid a radiative collapse in the presence of intrinsic tungsten and seeded xenon impurities, as quantitatively estimated in this work. As the plasma current is reduced, the current density becomes increasingly peaked, reflected by a growing value of the internal inductance $\ell$ $i3$ , resulting in reduced controllability of the vertical position of the plasma. The feasibility of different plasma current ramp-down rates is tested by applying an automated optimization framework embedding the RAPTOR core transport solver. Optimal time traces for plasma current $I$ $p$ $(t)$ and plasma elongation $κ(t)$ are proposed, to satisfy an $I$ $p$ -dependent upper limit on the plasma internal inductance, as obtained from vertical stability studies using the CREATE-NL code, as well as a constraint on the time evolution of $q$ 95 , to avoid an ideal MHD mode. A negative current density near the plasma edge is observed in our simulations, even for the most conservative $I$ $p$ ramp-down rate, indicating significant transient dynamics due to a large resistive time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Time-dependent outgassing and impurities in the NASA Lewis Bumpy Torus

To investigate the importance to fusion devices of desorption of gas from walls or cryosurfaces under long time operation, preliminary mass spectrometric investigations were carried out on the NASA Lewis Bumpy Torus with a 0.1-2A discharge in D2 at a pressure of 4.0 to 6.6 x 0.00001 torr. During the initial tens of minutes of discharge operation large quantities of a component with mass number 28, believed to be nitrogen, were released. The decrease with time of the nitrogen density after its initial maximum, indicated a diffusion controlled outgassing from cryodeposits or wall material. When cold surfaces were allowed to warm up, large quantities of gas were released. The integrated amount of N2 was typically the same as that released during a discharge, while the amount of H2O was orders of magnitude larger.

Persson, H.↗

Advances in the long-pulse steady-state high beta H-mode scenario with active controls of divertor heat and particle fluxes in EAST

Since the last IAEA-Fusion Energy Conference, the Experimental Advanced Superconducting Tokamak (EAST) research program has been, in support of ITER and CFETR, focused on development in terms of the long-pulse steady-state (fully noninductive) high beta H-mode scenario with active controls of the stationary and transient divertor heat and particle fluxes. The operational domain of the steady-state H-mode plasma scenario has been significantly extended with ITER-like tungsten mono-block divertor, plasma control and heating schemes. EAST has achieved several important milestones in the development of high β p H-mode scenario and its key physics and technologies. A 60 s-scale long-pulse steady-state high β p H-mode discharge with the major normalized plasma parameters similar to the designed performance of the CFETR 1GW fusion power operation scenario has been successfully established and sustained by pure RF heating and current drive. Several feedback control schemes have been developed for a sustained detachment with good core confinement. This includes control of the total radiation power, target electron temperature, and particle flux measured using divertor Langmuir probes or a combination of the control of target electron temperature and AXUV radiation near the X point. The detachment feedback control schemes have been integrated with small-ELM regimes and high β p scenario via neon seeding, enabling a core and edge compatible integrated high-beta scenario applicable to long-pulse operations. ELM suppression has been achieved using various methods, including resonant magnetic perturbations and impurity seeding. Full suppression of ELMs by using n=4 RMPs has been demonstrated for ITER for the first time in low input torque plasmas in EAST. Finally, EAST has been operated with helium to support the ITER research requirements for the first time. For a long-pulse, high bootstrap current fraction operation, a new lower tungsten divertor with active water-cooling has been installed, along with improvements in the heating and current drive capability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Implementation of AI/DEEP learning disruption predictor into a plasma control system

Abstract This paper reports on advances in the state‐of‐the‐art deep learning disruption prediction models based on the Fusion Recurrent Neural Network (FRNN) originally introduced in a 2019 NATURE publication [ https://doi.org/10.1038/s41586‐019‐1116‐4 ]. In particular, the predictor now features not only the “disruption score,” as an indicator of the probability of an imminent disruption, but also a “sensitivity score” in real time to indicate the underlying reasons for the imminent disruption. This adds valuable physics interpretability for the deep learning model and can provide helpful guidance for control actuators now implemented into a modern plasma control system (PCS). The advance is a significant step forward in moving from modern deep learning disruption prediction to real‐time control and brings novel AI‐enabled capabilities relevant for application to the future burning plasma ITER system. Our analyses use large amounts of data from JET and DIII‐D vetted in the earlier NATURE publication. In addition to “when” a shot is predicted to disrupt, this paper addresses reasons “why” by carrying out sensitivity studies. FRNN is accordingly extended to use more channels of information, including measured DIII‐D signals such as (i) the “n1rms” signal that is correlated with the n = 1 modes with finite frequency, including neoclassical tearing mode and sawtooth dynamics; (ii) the bolometer data indicative of plasma impurity control; and (iii) “q‐min”—the minimum value of the safety factor relevant to the key physics of kink modes. The additional channels and interpretability features expand the ability of the deep learning FRNN software to provide information about disruption subcategories as well as more precise and direct guidance for the actuators in a PCS.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A review of lithium application for the plasma-facing material in EAST Tokamak

As a low-Z material, lithium is widely tested and used in tokamak devices for particle recycling and impurity control to improve plasma performance. Lithium application as plasma-facing materials via different methods, such as coating, injection and flowing liquid limiter, has been investigated more than 10 years on EAST. This study presents the development of the lithium techniques applied on EAST, and reviews the progresses on the investigation of its interaction with plasma and its effects on the improvement on the plasma performances. In addition, this paper further discusses the potential applications of lithium technique in future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Divertor detachment and heat exhaust mitigation control in KSTAR with tungsten divertor

KSTAR has recently undergone an upgrade to use a new tungsten divertor to run experiments in ITER-relevant scenarios. Even with a high melting point of tungsten, it is important to control the heat flux impinging on tungsten divertor targets to minimize sputtering and contamination of the core plasma. Heat flux on the divertor is often controlled by increasing the degree of detachment of scrape-off layer plasma from the target plates. In this work, we have demonstrated successful divertor detachment and heat exhaust dissipation control experiments using two different methods. The first method uses attachment fraction as a control variable which is estimated using ion saturation current measurements from embedded Langmuir probes in the divertor. The second method uses a novel machine-learning-based surrogate model of 2D UEDGE simulation database, DivControlNN. We demonstrated running inference operation of DivControlNN in realtime to estimate heat flux at the divertor and use it as the control variable in a feedback loop with impurity gas flow. We present interesting insights from these experiments including a systematic approach to tuning controllers and discuss future improvements in the control infrastructure and control variables for future burning plasma experiments.

KSTAR tungsten divertor operations↗

Design Considerations to Ensure Robustness of the ITER Diagnostics Residual Gas Analyzer

Increasing robustness of the ITER diagnostic residual gas analyzer (DRGA) is critical for the potential control of plasma heating and fuel-cycle processing. Robustness is a requirement for a diagnostic to have a control function. The DRGA is a multisensor diagnostic system capable of resolving isotopic compositions of hydrogen and helium as well as other heavier elements and compounds. The divertor-specific DRGA system is intended to measure the composition of gases in the ITER subdivertor region and midplane. Its analysis station will be located in a port cell at the divertor level. From there, it will sample a slip stream of gas from the cryogenic pump duct. It will then exhaust into a shared roughing line where helium or other light gas impurities, some potentially from other diagnostic systems, are likely to be present. In order to provide reliable measurements, the DRGA must be robust in areas such as plasma optical emission source geometry for a compact design, mitigation of back-streaming from light gases and resilience to the ITER port cell environment, and radiation hardening of the DRGA electronics. By incorporating robustness into the design, areas of plasma heating and fuel-cycle control may be explored with the DRGA for ITER and next-generation fusion devices.

Quinlan, Brendan↗

A reaction mechanism for plasma electrolysis of AgNO 3 forming silver nanoclusters and nanoparticles

In plasma-driven solution electrolysis (PDSE), gas-phase plasma-produced species interact with an electrolytic solution to produce, for example, nanoparticles. An atmospheric pressure plasma jet (APPJ) directed onto a liquid solution containing a metallic salt will promote reduction of metallic ions in solution, generating metallic clusters that nucleate to form nanoparticles. In this article, results from a computational investigation are discussed of a PDSE process in which a radio-frequency APPJ sustained in helium impinges on a silver nitrate solution, resulting in growth of silver nanoparticles. A reaction mechanism was developed and implemented in a global plasma chemistry model to predict nanoparticle growth. To develop the reaction mechanism, density functional theory was used to generate probable silver growth pathways up to Ag 9 . Neutral clusters larger than Ag 9 were classified as nanoparticles. Kinetic reaction rate coefficients for thermodynamically favorable growth pathways were estimated based on an existing, empirically determined base reaction mechanism for smaller Ag particle interactions. These rates were used in conjunction with diffusion-controlled reaction rate coefficients that were calculated for other Ag species. The role of anions in reduction of Agn ions in forming nanoparticles is also discussed. Oxygen containing impurities or admixtures to the helium, air entrainment into the APPJ, and dissociation of saturated water vapor above the solution can produce additional reactive oxygen species in solution, resulting in the production of anions and O 2 – in particular. For a given molarity, delivering a sufficient fluence of reducing species will produce similar nanoparticle densities and sizes for all applied power levels. In conclusion, comparisons are made to alternate models for nanoparticle formation, including charged nanoparticles and use of direct current plasmas.

36 MATERIALS SCIENCE↗

Overview of results from the 2023 DIII-D negative triangularity campaign

Negative triangularity (NT) is a potentially transformative configuration for tokamak-based fusion energy with its high-performance core, edge localized mode (ELM)-free edge, and low-field-side divertors that could readily scale to an integrated reactor solution. Previous NT work on the TCV and DIII-D tokamaks motivated the installation of graphite-tile armor on the low-field-side lower outer wall of DIII-D. A dedicated multiple-week experimental campaign was conducted to qualify the NT scenario for future reactors. During the DIII-D NT campaign, high confinement ($H_\mathrm{98y,2}\gtrsim$ 1), high current ($q_\mathrm{95} \lt$ 3), and high normalized pressure plasmas ($\beta_\mathrm{N} \gt$ 2.5) were simultaneously attained in strongly NT-shaped discharges with average triangularity $\delta_\mathrm{avg}$ = −0.5 that were stably controlled. Experiments covered a wide range of DIII-D operational space (plasma current, toroidal field, electron density and pressure) and did not trigger an ELM in a single discharge as long as sufficiently strong NT was maintained; in contrast, to other high-performance ELM-suppression scenarios that have narrower operating windows. These strong NT plasmas had a lower outer divertor X-point shape and maintained a non-ELMing edge with an electron temperature pedestal, exceeding that of typical L-mode plasmas. Also, the following was achieved during the campaign: high normalized density ($n_\mathrm{e}$/$n_\mathrm{GW}$ of at least 1.7), particle confinement comparable to energy confinement with $Z_\mathrm{eff}\sim2$ , a detached divertor without impurity seeding, and a mantle radiation scenario using extrinsic impurities. These results are promising for a NT fusion pilot plant but further questions on confinement extrapolation and core-edge integration remain, which motivate future NT studies on DIII-D and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Observation of fully detached divertor integrated with improved core confinement for tokamak fusion plasmas

Integration of divertor detachment with a high-performance (β N ~ 3, β p > 2, H 98 ~ 1.5) core plasma has been demonstrated in DIII-D high-β p (poloidal beta) plasmas associated with a sustained core internal transport barrier (ITB) and an H-mode edge transport barrier (ETB). Such good core-edge integration has been achieved for both neon and nitrogen seeding, for both favorable and unfavorable B-field directions, independently from the impurity puffing locations, though these variations play important roles on divertor characteristics. Compared to the standard H-mode plasmas, the high-β p plasma exhibits a much wider window of detachment compatible with high confinement core. Fully detached divertor plasmas with low plasma temperature (Te < 5 eV), low particle flux, and low heat flux across the entire divertor target plate were obtained by using nitrogen seeding. This detached high-β p plasma is compatible with a newly developed detachment control system which can help optimize the nitrogen gas flow rate. Several features, i.e., the high edge safety factor in the high-β p scenario, impurity injection, closed divertor and reduced heating power requirement due to the high confinement, facilitate the achievement of full divertor detachment at lower density. Instead of degrading global performance, the divertor detachment facilitates the access to an even stronger ITB at large radius with a relatively weak ETB through self-organized synergy between ITB and ETB, leading to sustained high confinement. The strengthening of the large-radius ITB compensates for the ETB degradation associated with divertor detachment. In addition, a weak ETB naturally has smaller edge localized modes (ELMs). In particular, with neon injection, a long-period no-ELM H-mode phase has been achieved simultaneously with high-performance core and partially detached divertor plasmas. Furthermore, these results demonstrate the possibility of integrating excellent core plasma performance with an effective divertor solution, an essential step toward steady-state operation of reactor-grade plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Overview of physics results from the ADITYA-U tokamak and future experiments

The ADITYA upgrade (ADITYA-U), a medium-sized (R 0 = 75 cm, a = 25 cm) conventional tokamak facility in India, has been consistently producing experiments findings by usingcircular and shaped-plasmas. Recognizing the plasma parameters aligning closely with the design parameters of circular limited plasmas, ADITYA-U shifted its focus toward exploring the operational regime for experimentation on saw-tooth and MHD phenomena. Moreover, ADITYA-U has made consistent advancements toward conducting preliminary plasma shaping experiments through the activation of top and bottom divertor coils utilizing hydrogen as well as deuterium fuels. Confinement is improved by a factor of ~1.5 in D 2 plasmas when compared to H 2 plasmas of ADITYA-U. Further, ADITYA-U operations emphasize preventing disruptions and runaway electrons (REs) to ensure safe operations for future fusion devices. Significant suppression of REs has been achieved in ADITYA-U with the application of pulsed localized vertical magnetic field (LVF) perturbation, thereby establishing the technique's independence from the tokamak device. The successful RE mitigation requires a critical threshold of LVF pulse magnitude, which is approximately 1% of the toroidal magnetic field, and a minimum duration of ~5 ms. Apart from this, several novel findings have been achieved in the ADITYA-U experiments, including the modification of sawtooth duration through gas-puff, the emergence of MHD-induced geodesic acoustic mode-like oscillations, the propagation of fast heat pulses induced by MHD activity, the control of RE dynamics through Gas-puffs, the propagation of pinch-driven cold-pulses, the transport and core accumulations of argon impurities, the mass dependency of plasma toroidal rotation and the detection of 'RICE' scaling, as well as the characterization of edge plasma using wall conditioning methods, such as glow discharge cleaning using a combination of Ar-H 2 mixture, localized wall cleaning by electron cyclotron resonant plasma, and the development of machine learning-based disruption predictions, will be discussed in this paper.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The importance of properties in modeling

Casting and welding of superalloys, stainless steel and titanium alloys are processes which can be improved through modeling of heat flow, fluid flow, residual stress development, and microstructural evolution. These simulations require inputs of thermophysical data, some of which involves the partially or totally liquid state. In particular, these processes involve melting, flow in the liquid, and solidification. Modeling of such processes can lead to an improved understanding of defects such as shrinkage, inclusions, cracks, incomplete filling (or penetration), macrosegregation, improper grain structure, and deviations from dimensional specifications. Effective modeling can shorten process development time and improve quality. An approach to these problems is to develop efficient models; validate through correlations with thermal, distortion, and microstructural data; run parametric studies; extract knowledge based rules; and apply to adaptive closed loop control systems. With the appropriate pre- and post-processing, such analyses can be made 'user friendly'. This would include graphical user interfaces as well as realistic images and color maps. In such form, these models can be used for sensitivity analyses, which are useful in defining appropriate sensors and in the development of control strategies. Such modeling can be done at several levels, e.g., the MARO level, modeling large scale phenomena such as heat and fluid flow or material deformation; the MICRO level, modeling the development of dendrites, grains or precipitates; or at the NANO level, modeling point defects, dislocations, stacking faults, etc. There are many computational issues associated with these simulations, e.g. computational efficiency and accuracy. In addition, there are many materials issues, not the least of which is the availability of accurate high temperature thermophysical data for complex alloys. This would include latent heat of fusion, temperature dependent heat capacity and thermal conductivity (for liquid and solid), viscosity, surface tension, thermal expansion, mechanical properties, etc. Preliminary data is frequently gathered from the literature; however, this is often not available for modern alloys. If additional data are required, measurements can be used; however, these are costly, time consuming and can be erroneous due to a lack of testing standards or impure materials. Microstructural predictors can be extracted from thermal information, e.g. cooling rate and thermal gradient; the prediction of microstructure is dependent on solidus and liquidus temperature, mushy zone permeability, the solidification curve, volume changes, phase transformations, alloying effects (such as surface tension or viscosity), mold/metal reactions, metal/environment reactions, etc. Defect maps may be needed to predict the onset of shrinkage, hot cracking or 'freckling'. Constants may be needed for stress relaxation, dendrite coarsening, vaporization, etc. Visualization was used as a tool to better comprehend complex data sets associated with the analysis of directional solidification (including crystal growth) and welding. Examples include not only isotherms, but also cooling rate, growth rate and thermal gradient. The latter two are not single valued scalars, but rather time and space dependent vector fields. Efficient models were developed for both casting and welding to predict heat flow and the relationship to dendrite and grain growth. These codes include many of the non-linear effects, e.g. radiation, which dominate these processes. The home-built FDM code(s) were designed to be useful not only to the scientist, but also to the process engineer. Special output can be requested to compare directly to experimental data. Visualization procedures were developed to visualize critical results, e.g. fusion zone width at the surface opposite that where the arc is applied ('penetration'). Both elaborate and simplified distortion analyses were carried out. It is clear that extensive mechanical property data are critical in order to accurately predict residual stress patterns. A scheme is currently being developed to integrate these modeling tools into a set of control algorithms; however, the success of this approach is critically dependent on the availability of accurate high temperature thermophysical data.

Giamei, A. F.↗

Time-dependent SOLPS-ITER simulations of the tokamak plasma boundary for model predictive control using SINDy *

Abstract Time-dependent SOLPS-ITER simulations have been used to identify reduced models with the sparse identification of nonlinear dynamics (SINDy) method and develop model-predictive control of the boundary plasma state using main ion gas puff actuation. A series of gas actuation sequences are input into SOLPS-ITER to produce a dynamic response in upstream and divertor plasma quantities. The SINDy method is applied to identify reduced linear and nonlinear models for the electron density at the outboard midplane n e , s e p O M P and the electron temperature at the outer divertor T e , s e p d i v . Note that T e , s e p d i v is not necessarily the peak value of T e along the divertor. The identified reduced models are interpretable by construction (i.e. not black box), and have the form of coupled ordinary differential equations. Despite significant noise in T e , s e p d i v , the reduced models can be used to predict the response over a range of actuation levels to a maximum deviation of 0.5% in n e , s e p O M P and 5%–10% in T e , s e p d i v for the cases considered. Model retraining using time history data triggered by a preset error threshold is also demonstrated. A model predictive control strategy for nonlinear models is developed and used to perform feedback control of a SOLPS-ITER simulation to produce a setpoint trajectory in n e , s e p O M P using the integrated plasma simulator framework. The developed techniques are general and can be applied to time-dependent data from other boundary simulations or experimental data. Ongoing work is extending the approach to model identification and control for divertor detachment, which will present transient nonlinear behavior from impurity seeding, including realistic latency and synthetic diagnostic signals derived from the full SOLPS-ITER output.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Power handling in a highly-radiative negative triangularity pilot plant

Abstract This work explores detailed power handling solutions for a class of high-field, highly-radiative negative triangularity (NT) reactors based around the MANTA concept (Rutherford et al 2024 Plasma Phys. Control. Fusion ). The divertor design is kept as simple as possible, opting for a standard divertor with standard leg length. FreeGS is used to create an equilibrium for the boundary region, prioritizing a short outer leg length of only ∼50 cm (∼40% of the minor radius). The UEDGE code package is used for the boundary plasma solution to track plasma temperatures and fluxes to the divertor targets. It is found that for P SOL = 25 MW and n sep = 0.96 × 10 20 m −3 , conditions consistent with initial core transport modeling, little additional power mitigation is necessary. For a fixed impurity fraction of just 0.13% Ne in the plasma, the peak heat flux density at the more heavily loaded outer targets falls to 7.8 MW m −2 , while the electron temperature T e remains just under 5 eV. Scans around the parameter space reveal that even at densities lower than in the primary operating scenario, P SOL can be increased up to 50 MW, so long as a slightly higher fraction of extrinsic radiator is used. With less than 1% neon (Ne) impurity content, the divertor still experiences less than 10 MW m −2 at the outer target. Design of the plasma-facing components includes a close-fitting vacuum vessel with a tungsten inner surface as well as FLiBe-carrying cooling channels fashioned into the VV wall directly behind the divertor targets. For the seeded heat flux profile, Ansys Fluent heat transfer simulations estimate that the outer target temperature remains at just below 1550 ∘ C. Initial scoping of advanced divertor designs shows that for an X-divertor, detachment of the outer target becomes much simpler, and plasma fluxes to the targets drop considerably with only 0.01% Ne content.

Miller, M. A. (ORCID:0000000265406533)↗

Physics insights from a large-scale 2D UEDGE simulation database for detachment control in KSTAR

A large-scale database of two-dimensional UEDGE simulations has been developed to study detachment physics in KSTAR and to support surrogate models for control applications. Nearly 70 000 steady-state solutions were generated, systematically scanning upstream density, input power, plasma current, impurity fraction, and anomalous transport coefficients, with magnetic and electric drifts across the magnetic field included. The database identifies robust detachment indicators, with strike-point electron temperature at detachment onset consistently T e,target ~ 3-4 eV, largely insensitive to upstream conditions. Scaling relations reveal weaker impurity sensitivity than one-dimensional models and show that heat flux widths follow Eich’s scaling only for uniform, low D and χ. Distinctive in–out divertor asymmetries are observed in KSTAR, differing qualitatively from DIII-D. Complementary time-dependent simulations quantify plasma response to gas puffing, with delays of 5 - 15 ms at the outer strike point and ∼ 40 ms for the low-magnetic-field-side radiation front. These dynamics are well captured by first-order-plus-dead-time models and are consistent with experimentally observed detachment-control behavior in KSTAR (Gupta et al 2025 Plasma Phys. Control. Fusion (submitted)).

Physics - Plasma physics↗

Aligning thermal and current quenches with a high density low-Z injection

The conventional approach for thermal quench (TQ) mitigation in a tokamak disruption is through a high- Z impurity injection that radiates away the plasma’s thermal energy before it reaches the wall. The downside is a robust Ohmic-to-runaway current conversion due to the radiatively clamped low post-thermal-quench electron temperature. An alternative approach is to deploy a low- Z (either deuterium or hydrogen) injection that aims to slow down the TQ, and ideally aligns it with the current quench (CQ). This approach has been investigated here via 3D MHD simulations using the PIXIE3D code. By boosting the hydrogen density, a fusion-grade plasma is dilutionally cooled at approximately the original pressure. Energy loss to the wall is controlled by a Bohm outflow condition at the boundary where the magnetic field intercepts a thin plasma sheath at the wall, in addition to Bremsstrahlung bulk losses. Robust MHD instabilities proceed as usual, while the collisionality of the plasma has been greatly increased and parallel transport is now in the Braginskii regime. The main conclusion of this study is that the decreased transport loss along open field lines due to a sufficient low- Z injection slows down the TQ rate to the order of 20 ms, aligned with the CQ timescale for a 15 MA ITER plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗