According to recent DIII-D experiments (Logan et al 2024 Nucl. Fusion64 014003), injecting edge localized electron cyclotron current drive (ECCD) in the counter-plasma-current (counter-Ip) direction reduces the n = 3 resonant magnetic perturbation (RMP) current threshold for edge-localized mode (ELM) suppression, while co-Ip ECCD during the suppressed ELM phase causes a back transition to ELMing. This paper presents nonlinear two-fluid simulations on the ECCD manipulation of edge magnetic islands induced by RMP using the TM1 code. In the presence of a magnetic island chain at the pedestal-top, co-Ip ECCD is found to decrease the island width and restore the initially degraded pedestal pressure when its radial deposition location is close to the rational surface of the island. With a sufficiently strong co-Ip ECCD current, the RMP-driven magnetic island can be healed, and the pedestal pressure fully recovers to its initial ELMing state. On the contrary, counter-Ip ECCD is found to increase the island width and further reduce the pedestal pressure to levels significantly below the peeling-ballooning-mode limited height, leading to even stationary ELM suppression. These simulations align with the results from DIII-D experiments. However, when multiple magnetic island chains are present at the pedestal-top, the ECCD current experiences substantial broadening, and its effects on the island width and pedestal pressure become negligible. Further simulations reveal that counter-Ip ECCD enhances RMP penetration by lowering the penetration threshold, with the degree of reduction proportional to the amplitude of ECCD current. For the ∼1 MW ECCD in DIII-D, the predicted decrease in the RMP penetration threshold for ELM suppression is approximately 20%, consistent with experimental observations. These simulations indicate that edge-localized ECCD can be used to either facilitate RMP-driven ELM suppression or optimize the confinement degradation.
ISSN: 1741-4326
Nuclear Fusion is the acknowledged world-leading journal specializing in fusion. The journal covers all aspects of research, theoretical and practical, relevant to controlled thermonuclear fusion.
Open all abstracts, in this tab
Q.M. Hu et al 2024 Nucl. Fusion 64 046027
Sehila M. Gonzalez de Vicente et al 2022 Nucl. Fusion 62 085001
In the absence of official standards and guidelines for nuclear fusion plants, fusion designers adopted, as far as possible, well-established standards for fission-based nuclear power plants (NPPs). This often implies interpretation and/or extrapolation, due to differences in structures, systems and components, materials, safety mitigation systems, risks, etc. This approach could result in the consideration of overconservative measures that might lead to an increase in cost and complexity with limited or negligible improvements. One important topic is the generation of radioactive waste in fusion power plants. Fusion waste is significantly different to fission NPP waste, i.e. the quantity of fusion waste is much larger. However, it mostly comprises low-level waste (LLW) and intermediate level waste (ILW). Notably, the waste does not contain many long-lived isotopes, mainly tritium and other activation isotopes but no-transuranic elements. An important benefit of fusion employing reduced-activation materials is the lower decay heat removal and rapid radioactivity decay overall. The dominant fusion wastes are primarily composed of structural materials, such as different types of steel, including reduced activation ferritic martensitic steels, such as EUROFER97 and F82H, AISI 316L, bainitic, and JK2LB. The relevant long-lived radioisotopes come from alloying elements, such as niobium, molybdenum, nickel, carbon, nitrogen, copper and aluminum and also from uncontrolled impurities (of the same elements, but also, e.g. of potassium and cobalt). After irradiation, these isotopes might preclude disposal in LLW repositories. Fusion power should be able to avoid creating high-level waste, while the volume of fusion ILW and LLW will be significant, both in terms of pure volume and volume per unit of electricity produced. Thus, efforts to recycle and clear are essential to support fusion deployment, reclaim resources (through less ore mining) and minimize the radwaste burden for future generations.
I.A.M. Datta et al 2024 Nucl. Fusion 64 066016
The FuZE sheared-flow-stabilized Z pinch at Zap Energy is simulated using whole-device modeling employing an axisymmetric resistive magnetohydrodynamic formulation implemented within the discontinuous Galerkin WARPXM framework. Simulations show formation of Z pinches with densities of approximately 1022 m−3 and total DD fusion neutron rate of 107 per µs for approximately 2 µs. Simulation-derived synthetic diagnostics show peak currents and voltages within 10% and total yield within approximately 30% of experiment for similar plasma mass. The simulations provide insight into the plasma dynamics in the experiment and enable a predictive capability for exploring design changes on devices built at Zap Energy.
G. Federici et al 2024 Nucl. Fusion 64 036025
High temperature superconductors (HTSs) offer the promise of operating at higher magnetic field and temperature. Recently, the use of high field magnets (by adopting HTS) has been promoted by several groups around the world, including new start-up entries, both to substantially reduce the size of a fusion power reactor system and as a breakthrough innovation that could dramatically accelerate fusion power deployment. This paper describes the results of an assessment to understand the impact of using high field magnets in the design of DEMO in Europe, considering a comprehensive list of physics and engineering limitations together with the interdependencies with other important parameters. Based on the results, it is concluded that increasing the magnetic field does not lead to a reduction in device size with relevant nuclear performance requirements, because (i) large structures are needed to withstand the enormous electromagnetic forces, (ii) thick blanket and n-shield structures are needed to protect the coils from radiation damage effects, and (iii) new divertor solutions with performances well beyond today's concepts are needed. Stronger structural materials allow for more compact tokamaks, but do not change the conclusion that scalability is not favourable when increasing the magnetic field, beyond a certain point, the machine size cannot be further reduced. More advanced structural support concepts for high-field coils have been explored and concluded that these solutions are either unfeasible or provide only marginal size reduction, by far not sufficient to account for the potential of operating at very high field provided by HTS. Additionally, the cost of high field coils is significant at today's price levels and shows to scale roughly with the square of the field. Nevertheless, it is believed that even when not operated at high field and starting within conventional insulated coils, HTS can still offer certain benefits. These include the simplification of the magnet cooling scheme thanks to increased temperature margin (indirect conduction cooling). This in turn can greatly simplify coil construction and minimize high-voltage risks at the terminals.
J. Mailloux et al 2022 Nucl. Fusion 62 042026
The JET 2019–2020 scientific and technological programme exploited the results of years of concerted scientific and engineering work, including the ITER-like wall (ILW: Be wall and W divertor) installed in 2010, improved diagnostic capabilities now fully available, a major neutral beam injection upgrade providing record power in 2019–2020, and tested the technical and procedural preparation for safe operation with tritium. Research along three complementary axes yielded a wealth of new results. Firstly, the JET plasma programme delivered scenarios suitable for high fusion power and alpha particle (α) physics in the coming D–T campaign (DTE2), with record sustained neutron rates, as well as plasmas for clarifying the impact of isotope mass on plasma core, edge and plasma-wall interactions, and for ITER pre-fusion power operation. The efficacy of the newly installed shattered pellet injector for mitigating disruption forces and runaway electrons was demonstrated. Secondly, research on the consequences of long-term exposure to JET-ILW plasma was completed, with emphasis on wall damage and fuel retention, and with analyses of wall materials and dust particles that will help validate assumptions and codes for design and operation of ITER and DEMO. Thirdly, the nuclear technology programme aiming to deliver maximum technological return from operations in D, T and D–T benefited from the highest D–D neutron yield in years, securing results for validating radiation transport and activation codes, and nuclear data for ITER.
K.C. Shaing et al 2024 Nucl. Fusion 64 066014
Transport consequences of the wave–particle interactions in the quasilinear plateau (QP) regime are presented. Eulerian approach is adopted to solve the drift kinetic equation that includes the physics of the nonlinear trapping (NT) and QP regimes. The localization of the perturbed distribution simplifies the test particle collision operator. It is shown that a mirror force like term responsible for the flattening of the distribution in the NT regime is subdominant in the QP regime, and controls the transition between these two regimes. Transport fluxes, flux-power relation, and nonlinear damping or growth rate are all calculated. There is no explicit collision frequency dependence in these quantities; however, the width of the resonance does. Formulas that join the asymptotic results of these two regimes to facilitate thermal and energetic particle transport, and nonlinear wave evolution of a single mode are presented.
Vignesh Gopakumar et al 2024 Nucl. Fusion 64 056025
Predicting plasma evolution within a Tokamak reactor is crucial to realizing the goal of sustainable fusion. Capabilities in forecasting the spatio-temporal evolution of plasma rapidly and accurately allow us to quickly iterate over design and control strategies on current Tokamak devices and future reactors. Modelling plasma evolution using numerical solvers is often expensive, consuming many hours on supercomputers, and hence, we need alternative inexpensive surrogate models. We demonstrate accurate predictions of plasma evolution both in simulation and experimental domains using deep learning-based surrogate modelling tools, viz., Fourier neural operators (FNO). We show that FNO has a speedup of six orders of magnitude over traditional solvers in predicting the plasma dynamics simulated from magnetohydrodynamic models, while maintaining a high accuracy (Mean Squared Error in the normalised domain ). Our modified version of the FNO is capable of solving multi-variable Partial Differential Equations, and can capture the dependence among the different variables in a single model. FNOs can also predict plasma evolution on real-world experimental data observed by the cameras positioned within the MAST Tokamak, i.e. cameras looking across the central solenoid and the divertor in the Tokamak. We show that FNOs are able to accurately forecast the evolution of plasma and have the potential to be deployed for real-time monitoring. We also illustrate their capability in forecasting the plasma shape, the locations of interactions of the plasma with the central solenoid and the divertor for the full (available) duration of the plasma shot within MAST. The FNO offers a viable alternative for surrogate modelling as it is quick to train and infer, and requires fewer data points, while being able to do zero-shot super-resolution and getting high-fidelity solutions.
P. Rodriguez-Fernandez et al 2022 Nucl. Fusion 62 042003
The SPARC tokamak project, currently in engineering design, aims to achieve breakeven and burning plasma conditions in a compact device, thanks to new developments in high-temperature superconductor technology. With a magnetic field of 12.2 T on axis and 8.7 MA of plasma current, SPARC is predicted to produce 140 MW of fusion power with a plasma gain of Q ≈ 11, providing ample margin with respect to its mission of Q > 2. All tokamak systems are being designed to produce this landmark plasma discharge, thus enabling the study of burning plasma physics and tokamak operations in reactor relevant conditions to pave the way for the design and construction of a compact, high-field fusion power plant. Construction of SPARC is planned to begin by mid-2021.
Semin Joung et al 2024 Nucl. Fusion 64 066038
A neural network, BES-ELMnet, predicting a quasi-periodic disruptive eruption of the plasma energy and particles known as edge localized mode (ELM) onset is developed with observed pedestal turbulence from the beam emission spectroscopy system in DIII-D. BES-ELMnet has convolutional and fully-connected layers, taking two-dimensional plasma fluctuations with a temporal window of size 128 µs and generating a scalar output which can be interpreted as a probability of the upcoming ELM onset. As approximately labeled inter-ELM broadband () fluctuations are given to the network, BES-ELMnet learns by itself ELM-related precursors arising before the onsets through supervised learning scheme. BES-ELMnet achieves the gradually increasing ELM onset probabilities between two consecutive ELMs during the inter-ELM phases and can forecast the first ELM onsets which occur after the high confinement mode transition. We further investigate the network generality in terms of the selected frequency band to ensure the use of BES-ELMnet for various operation regimes without changing the trained architecture. Therefore, our novel prediction method will enhance a proactive high confinement mode control of fusion-grade plasmas.
M. Hoelzl et al 2021 Nucl. Fusion 61 065001
JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.
Open all abstracts, in this tab
S. Varoutis et al 2024 Nucl. Fusion 64 076011
The present work presents a 2D and 3D modeling of the neutral gas flow in the sub-divertor region of the W7-X. The investigations have been done using the DIVGAS code. The complex 2D and 3D geometries of the divertor components in the sub-divertor region have been considered and the Standard and High-Iota magnetic configurations have been numerically simulated. The main objective of this study is to investigate the dynamics of neutral particles in the sub-divertor region including the effects due to geometry and toroidal and poloidal leakages located at the divertor targets and baffles on the achieved pumping efficiency. A sensitivity analysis has been performed for the effect of various geometrical and flow parameters on the pumping performance, under different plasma scenarios. The considered incoming fluxes in the sub-divertor range between 1020 to 1022 (H2 s−1). The main conclusions, which can be extracted from the present numerical analysis could be summarized as follows; a large fraction of incoming neutral particle flux i.e. >70% on the low iota side and >40% for the high iota side is leaked back to the main divertor region, while higher incoming neutral fluxes facilitate the increase of the pumped flux as well as the decrease of the outflux. It has been estimated that a small fraction ∼3%–4% of the incoming neutral flux is being pumped via the turbo-molecular pumps. The closure of the toroidal leakages as well as the inclination of the pumping gap panel by 9o facilitate the increase of the pumped flux, but considering the all the engineering constraints, the latter option seems to be more easy to be implemented. For low incoming neutral fluxes (∼1020 H2 s−1) and for the case of AEH section, free molecular flow conditions are estimated and therefore neutral-neutral collisions could be neglected. For higher incoming neutral fluxes and for both AEH and AEP sections neutral-neutral collisions play a significant role in the flow establishment. A comparison with available experimental measurements of the neutral pressure in the sub-divertor has been performed for Standard and High-Iota plasma discharges. The 3D DIVGAS simulations predict qualitatively the experimental data with relative deviation between 25 and 63%. All the above numerical findings will actively support the optimization of the W7-X particle exhaust, in view of the experimental campaign OP2.
P. Liu et al 2024 Nucl. Fusion 64 076007
This paper reports global nonlinear gyrokinetic simulations that couple meso-scale reversed shear Alfvén eigenmodes (RSAEs) driven by energetic particles (EPs) and ion temperature gradient (ITG) microturbulence driven by thermal plasma, using equilibrium and profiles from DIII-D discharge #159243. In simulations focusing only on the ITG, electrostatic ITG drives a huge thermal ion heat transport, which is reduced by a factor of to a level close to the experimental value in electromagnetic simulation due to finite effect. In the simulations coupling the RSAE and ITG, ITG can scatter the resonant EP nonlinearly trapped by the RSAE and damp the zonal flows generated by the RSAE. The regulation of the RSAE by the ITG greatly reduces the initial saturation amplitude of the RSAE but increases the RSAE amplitude and associated EP transport to experimental levels in the quasi-steady state. The RSAE effects on the ITG, specifically the stronger zonal flows generated by the RSAE and the RSAE frequency modulation of the ITG-induced thermal ion heat transport, in turn, leads to a reduction of the thermal ion heat transport by more than a factor of . For a stronger background ITG, the regulation of the RSAE by the ITG is stronger, while the RSAE effects on the ITG are weaker. This work highlights the importance of cross-scale coupling in the dynamics of the AE turbulence and EP transport.
Xingquan Wu (伍兴权) et al 2024 Nucl. Fusion 64 076008
The magnetic coherent modes (MCM) with toroidal mode number n about 1 (Chen R. et al 2018 Nucl. Fusion58 112004) frequently appear in the edge pedestal of high-confinement tokamak plasmas on EAST in the absence of energetic particles. Although these modes are experimentally compatible with the steady-state operation of the pedestal, the driving mechanism without energetic particles of MCM is a long-standing mystery. To reveal the excitation mechanism, a fluid-drift kinetic hybrid local linear model has been developed. It is found that MCM is a new Alfvén eigenmode with a gap frequency much lower than the ideal Toroidal Alfvén Eigenmodes (TAEs) with two significant properties: (1) due to the unique steep pressure gradient in the pedestal region, the diamagnetic frequency becomes comparable to the ideal TAE frequency, which makes the Alfvén continuum in this region move significantly in the ion diamagnetic direction and form a gap of lower frequency; (2) due to the bounce frequencies of thermal electrons becoming also comparable to the ideal TAE frequency in the pedestal region, the free energy of the pressure gradient can be fed into the MCM through the thermal electron bounce resonance excitation, which is essentially the coupling between the shifted TAEs and low- trapped electron modes. The low- MCM is proved to be a shifted Alfvén gap mode in the pedestal region, which is anticipated to exist in low collisional plasmas of future fusion reactors. It is of great significance to carry out relevant physical model research to enhance the understanding of pedestal physics.
Bo S. Schmidt et al 2024 Nucl. Fusion 64 076009
Fast-ion loss detectors (FILDs) are crucial for analyzing fast-ion dynamics in magnetically confined fusion plasmas. A core challenge is to derive an accurate ion velocity distribution, requiring treatment of thousands of remapped camera frames for a full discharge. The ill-posed nature of this task necessitates regularization with a well-chosen regularization parameter and computationally efficient methods. In this work, we introduce the 'resolution principle,' a novel criterion for selecting the optimal regularization parameter, providing a distinction between genuine features and artefacts smaller than the diagnostic resolution in the reconstruction, thereby preventing misinterpretations. This principle, coupled with three iterative reconstruction techniques—Kaczmarz's method, coordinate descent, and Cimmino's method—demonstrates enhanced reconstruction capabilities compared to conventional methods like Tikhonov regularization. Utilizing these techniques allows rapid processing of measurements from full discharges, removing the computational bottleneck and facilitating between-discharge reconstructions. By reconstructing 6000 camera frames from an ELMy H-mode discharge at ASDEX Upgrade, we capture the temporal evolution of gyroradii and pitch angles, unveiling a direct correlation between pitch-angle behavior and changes in the toroidal magnetic field for a specific subset of lost ions accelerated by edge-localized modes (ELMs) to energies approximately twice that of the injection energy.
K. Ogawa et al 2024 Nucl. Fusion 64 076010
Energetic ion anisotropy was observed by tangential sightline compact neutron energy spectrometers (CNESs) in tangential neutral beam heated deuterium plasmas in Large Helical Device. Significant upper and lower energy shifts in D–D neutron energy from 2.45 MeV were measured according to the beam ion injection directions and CNES sightline using a conventional liquid scintillation detector with the unfolding technique and a novel Cs2LiYCl6:Ce with a 7Li-enrichment (CLYC7) scintillation detector without unfolding. The observed neutron energy spectrum was compared with that predicted by a numerical simulation based on orbit following models. Numerical simulation revealed that the Doppler shift in D–D neutron energy results from energetic ion anisotropy.
Open all abstracts, in this tab
G.D. Conway et al 2022 Nucl. Fusion 62 013001
Geodesic acoustic modes (GAMs) are ubiquitous oscillatory flow phenomena observed in toroidal magnetic confinement fusion plasmas, such as tokamaks and stellarators. They are recognized as the non-stationary branch of the turbulence driven zonal flows which play a critical regulatory role in cross-field turbulent transport. GAMs are supported by the plasma compressibility due to magnetic geodesic curvature—an intrinsic feature of any toroidal confinement device. GAMs impact the plasma confinement via velocity shearing of turbulent eddies, modulation of transport, and by providing additional routes for energy dissipation. GAMs can also be driven by energetic particles (so-called EGAMs) or even pumped by a variety of other mechanisms, both internal and external to the plasma, opening-up possibilities for plasma diagnosis and turbulence control. In recent years there have been major advances in all areas of GAM research: measurements, theory, and numerical simulations. This review assesses the status of these developments and the progress made towards a unified understanding of the GAM behaviour and its role in plasma confinement. The review begins with tutorial-like reviews of the basic concepts and theory, followed by a series of topic orientated sections covering different aspects of the GAM. The approach adopted here is to present and contrast experimental observations alongside the predictions from theory and numerical simulations. The review concludes with a comprehensive summary of the field, highlighting outstanding issues and prospects for future developments.
L. Marrelli et al 2021 Nucl. Fusion 61 023001
This paper reviews the research on the reversed field pinch (RFP) in the last three decades. Substantial experimental and theoretical progress and transformational changes have been achieved since the last review (Bodin 1990 Nucl. Fusion 30 1717–37). The experiments have been performed in devices with different sizes and capabilities. The largest are RFX-mod in Padova (Italy) and MST in Madison (USA). The experimental community includes also EXTRAP-T2R in Sweden, RELAX in Japan and KTX in China. Impressive improvements in the performance are the result of exploration of two lines: the high current operation (up to 2 MA) with the spontaneous occurrence of helical equilibria with good magnetic flux surfaces and the active control of the current profile. A crucial ingredient for the advancements obtained in the experiments has been the development of state-of-art active feedback control systems allowing the control of MHD instabilities in presence of a thin shell. The balance between achievements and still open issues leads us to the conclusion that the RFP can be a valuable and diverse contributor in the quest for fusion electricity.
Mohamed Abdou et al 2021 Nucl. Fusion 61 013001
The tritium aspects of the DT fuel cycle embody some of the most challenging feasibility and attractiveness issues in the development of fusion systems. The review and analyses in this paper provide important information to understand and quantify these challenges and to define the phase space of plasma physics and fusion technology parameters and features that must guide a serious R&D in the world fusion program. We focus in particular on components, issues and R&D necessary to satisfy three 'principal requirements': (1) achieving tritium self-sufficiency within the fusion system, (2) providing a tritium inventory for the initial start-up of a fusion facility, and (3) managing the safety and biological hazards of tritium. A primary conclusion is that the physics and technology state-of-the-art will not enable DEMO and future power plants to satisfy these principal requirements. We quantify goals and define specific areas and ideas for physics and technology R&D to meet these requirements. A powerful fuel cycle dynamics model was developed to calculate time-dependent tritium inventories and flow rates in all parts and components of the fuel cycle for different ranges of parameters and physics and technology conditions. Dynamics modeling analyses show that the key parameters affecting tritium inventories, tritium start-up inventory, and tritium self-sufficiency are the tritium burn fraction in the plasma (fb), fueling efficiency (ηf), processing time of plasma exhaust in the inner fuel cycle (tp), reactor availability factor (AF), reserve time (tr) which determines the reserve tritium inventory needed in the storage system in order to keep the plant operational for time tr in case of any malfunction of any part of the tritium processing system, and the doubling time (td). Results show that ηffb > 2% and processing time of 1–4 h are required to achieve tritium self-sufficiency with reasonable confidence. For ηffb = 2% and processing time of 4 h, the tritium start-up inventory required for a 3 GW fusion reactor is ∼11 kg, while it is <5 kg if ηffb = 5% and the processing time is 1 h. To achieve these stringent requirements, a serious R&D program in physics and technology is necessary. The EU-DEMO direct internal recycling concept that carries fuel directly from the plasma exhaust gas to the fueling systems without going through the isotope separation system reduces the overall processing time and tritium inventories and has positive effects on the required tritium breeding ratio (TBRR). A significant finding is the strong dependence of tritium self-sufficiency on the reactor availability factor. Simulations show that tritium self-sufficiency is: impossible if AF < 10% for any ηffb, possible if AF > 30% and 1% ⩽ ηffb ⩽ 2%, and achievable with reasonable confidence if AF > 50% and ηffb > 2%. These results are of particular concern in light of the low availability factor predicted for the near-term plasma-based experimental facilities (e.g. FNSF, VNS, CTF), and can have repercussions on tritium economy in DEMO reactors as well, unless significant advancements in RAMI are made. There is a linear dependency between the tritium start-up inventory and the fusion power. The required tritium start-up inventory for a fusion facility of 100 MW fusion power is as small as 1 kg. Since fusion power plants will have large powers for better economics, it is important to maintain a 'reserve' tritium inventory in the tritium storage system to continue to fuel the plasma and avoid plant shutdown in case of malfunctions of some parts of the tritium processing lines. But our results show that a reserve time as short as 24 h leads to unacceptable reserve and start-up inventory requirements. Therefore, high reliability and fast maintainability of all components in the fuel cycle are necessary in order to avoid the need for storing reserve tritium inventory sufficient for continued fusion facility operation for more than a few hours. The physics aspects of plasma fueling, tritium burn fraction, and particle and power exhaust are highly interrelated and complex, and predictions for DEMO and power reactors are highly uncertain because of lack of experiments with burning plasma. Fueling by pellet injection on the high field side of tokamak has evolved to be the preferred method to fuel a burning plasma. Extrapolation from the DIII-D penetration scaling shows fueling efficiency expected in DEMO to be <25%, but such extrapolations are highly uncertain. The fueling efficiency of gas in a reactor relevant regime is expected to be extremely poor and not very useful for getting tritium into the core plasma efficiently. Gas fueling will nonetheless be useful for feedback control of the divertor operating parameters. Extensive modeling has been carried out to predict burn fraction, fueling requirements, and fueling efficiency for ITER, DEMO, and beyond. The fueling rate required to operate Q = 10 ITER plasmas in order to provide the required core fueling, helium exhaust and radiative divertor plasma conditions for acceptable divertor power loads was calculated. If this fueling is performed with a 50–50 DT mix, the tritium burn fraction in ITER would be ∼0.36%, which is too low to satisfy the self-sufficiency conditions derived from the dynamics modeling for fusion reactors. Extrapolation to DEMO using this approach would also yield similarly low burn fraction. Extensive analysis presented shows that specific features of edge neutral dynamics in ITER and fusion reactors, which are different from present experiments, open possibilities for optimization of tritium fueling and thus to improve the burn fraction. Using only tritium in pellet fueling of the plasma core, and only deuterium for edge density, divertor power load and ELM control results in significant increase of the burn fraction to 1.8–3.6%. These estimates are performed with physics models whose results cannot be fully validated for ITER and DEMO plasma conditions since these cannot be achieved in present tokamak experiments. Thus, several uncertainties remain regarding particle transport and scenario requirements in ITER and DEMO. The safety standard requirements for protection of the public and release guidelines for tritium have been reviewed. General safety approaches including minimizing tritium inventories, reducing tritium permeation through materials, and decontaminating material for waste disposal have been suggested.
Boris N. Breizman et al 2019 Nucl. Fusion 59 083001
Of all electrons, runaway electrons have long been recognized in the fusion community as a distinctive population. They now attract special attention as a part of ITER mission considerations. This review covers basic physics ingredients of the runaway phenomenon and the ongoing efforts (experimental and theoretical) aimed at runaway electron (RE) taming in the next generation tokamaks. We emphasize the prevailing physics themes of the last 20 years: the hot-tail mechanism of runaway production, RE interaction with impurity ions, the role of synchrotron radiation in runaway kinetics, RE transport in presence of magnetic fluctuations, micro-instabilities driven by REs in magnetized plasmas, and vertical stability of the plasma with REs. The review also discusses implications of the runaway phenomenon for ITER and the current strategy of RE mitigation.
M.K.A. Thumm et al 2019 Nucl. Fusion 59 073001
In many tokamak and stellarator experiments around the globe that are investigating energy production via controlled thermonuclear fusion, electron cyclotron heating and current drive (ECH&CD) are used for plasma start-up, heating, non-inductive current drive and magnetohydrodynamic stability control. ECH will be the first auxiliary heating method used on ITER. Megawatt-class, continuous wave gyrotrons are employed as high-power millimeter (mm)-wave sources. The present review reports on the worldwide state-of-the-art of high-power gyrotrons. Their successful development during recent years changed ECH from a minor to a major heating method. After a general introduction of the various functions of ECH&CD in fusion physics, especially for ITER, section 2 will explain the fast-wave gyrotron interaction principle. Section 3 discusses innovations on the components of modern long-pulse fusion gyrotrons (magnetron injection electron gun, beam tunnel, cavity, quasi-optical output coupler, synthetic diamond output window, single-stage depressed collector) and auxiliary components (superconducting magnets, gyrotron diagnostics, high-power calorimetric dummy loads). Section 4 deals with present megawatt-class gyrotrons for ITER, W7-X, LHD, EAST, KSTAR and JT-60SA, and also includes tubes for moderate pulse length machines such as ASDEX-U, DIII-D, HL-2A, TCV, QUEST and GAMMA-10. In section 5 the development of future advanced fusion gyrotrons is discussed. These are tubes with higher frequencies for DEMO, multi-frequency (multi-purpose) gyrotrons, stepwise frequency tunable tubes for plasma stabilization, injection-locked and coaxial-cavity multi-megawatt gyrotrons, as well as sub-THz gyrotrons for collective Thomson scattering. Efficiency enhancement via multi-stage depressed collectors, fast oscillation recovery methods and reliability, availability, maintainability and inspectability will be discussed at the end of this section.
Open all abstracts, in this tab
den Harder et al
A low beamlet divergence is crucial for the efficiency of the ITER-NBI systems, since it affects the transmission of the beam through the duct. There is a requirement of 7 mrad e-folding divergence for the ITER Heating Neutral Beam. Significantly higher divergences (10-15 mrad) have been observed in RF-source based experiments albeit at low beam energy. This could be the consequence of a broad perpendicular velocity distribution of the H-/D- particles before extraction. This paper explores this hypothesis and its implications for ITER. To estimate H-/D- perpendicular temperatures in the RF-driven BATMAN Upgrade test facility, spatially resolved measurements of the beam power density are compared with IBSimu calculations. The estimated perpendicular temperatures show a strong dependence on the source filling pressure, decreasing from approximately 4 eV at 0.3 Pa to 2 eV at 0.4 Pa. Ion-optics calculations of the ITER-HNB grid system are performed to evaluate whether the temperatures estimated in the BATMAN Upgrade test facility are tolerable in view of beam-grid interaction and beamline transmission. The beamline transmission is fairly insensitive to the perpendicular temperature, but the heat loads at the downstream grids increase with the perpendicular temperature.
Frerichs et al
A subspace of resonant magnetic perturbation (RMP) configurations for edge localized mode (ELM) suppression is predicted for H-mode burning plasmas at 15 MA current and 5.3 T magnetic field in ITER. Perturbation of the core plasma can be reduced by a factor of 2 for equivalent edge stability proxies, while the perturbed plasma boundary geometry remains mostly resilient. The striation width of perturbed field lines connecting from the main plasma (normalized poloidal flux < 1) to the divertor targets is found to be significantly larger than the expected heat load width in the absence of RMPs. This facilitates heat load spreading with peak values at an acceptable level below 10 MW m−2 on the outer target already at moderate gas fueling and low Ne seeding for additional radiative dissipation of the 100 MW of power into the scrape-off layer (SOL). On the inner target, however, re-attachment is predicted away from the equilibrium strike point due to increased upstream heat flux, higher downstream temperature and less efficient impurity radiation.
Žohar et al
To computationally support hydrogen and helium plasma discharges in the early stages of tokamak operation and to support the commissioning of the neutron detectors during these operational phases, creation of a realistic neutron and gamma ray particle source for Monte Carlo simulations will be needed. One of the most important parts of creating the particle source is calculating the reaction rates of the particle-emitting reactions to determine the emission profile in the plasma and the energy spectra of the emitted particles. In this paper the analysis and evaluation of cross sections for important neutron-emitting reactions, namely, 9Be(p,nγ)9B, 9Be(3He,nγ) 11C and charged-particle emission reactions 9Be(p,d)2α and 9Be(p,α)6Li that cause neutron emission in the next step of interactions are presented. The reaction cross sections were evaluated based on experimental measurements and empirical models describing the interaction of two charged particles. Evaluation of the associated uncertainties was also performed. The main goal of the work is to propose the newly evaluated cross sections for inclusion in the FENDL nuclear data library, thus making the cross section available to other researchers studying the above listed reactions.
Marini et al
Current profile reconstructions are obtained for high current (Ip = 550 kA) post-disruption runaway electron (RE) plateau plasmas in DIII-D. Two novel methods of measuring the RE current profile in high-current RE plateaus are introduced and compared: localization of the q = 2 rational surface using visible synchrotron emission (SE) imaging and the measurement of the polarization angle of line-integrated Ar-II line emission. The two methods are found to be consistent with each other within the data uncertainties. Different simulations of the RE current profile are compared with the measurements: the toroidal fluid RE model is found to best fit the data, within the measurement uncertainties. In addition to introducing two novel methods to measure the RE current profile and validating present simulation capabilities, this work demonstrates that instabilities can grow at q = 2 and q = 1 surfaces without necessarily causing a RE final loss instability. Numerical simulations are also presented to elucidate the role of these instabilities on synchrotron emission.
Carloni et al
One of the main objectives of ITER is to produce 500 MW of power from a D-T plasma for several seconds. This goal presents two inherent challenges: firstly, in-vessel components will require active cooling to remove the heat coming from the fusion reaction (i.e., mainly fast neutrons and alpha particles). Secondly, the materials exposed to the neutron flux will yield activated corrosion products (ACPs) in all primary cooling circuits of ITER. From a safety point of view, ACPs are one of the contributors to the Occupational Radiological Exposure (ORE), they represent a source of radiological waste and also contribute to the source term for accidental scenarios involving the loss of primary confinement. 
Therefore, ACPs assessment is key to estimate radiological impact for nuclear workers and the public. ITER nuclear safety engineers adopted OSCAR-Fusion v1.4.a code to assess the ACPs inventory in the Integrated Blanket ELMs and Divertor (IBED) cooling loop. This paper describes the selection of input data, the modelling of the circuits and the operational scenarios used in OSCAR-Fusion calculations. This study also examines the outcomes of such calculations, notably in terms of ACPs inventory, emphasizing the impact on the ORE and highlighting its driving parameters. Furthermore, this paper provides recommendations for better ACPs management in the context of the ITER project and in accordance with the ALARA principle
Open all abstracts, in this tab
S. Varoutis et al 2024 Nucl. Fusion 64 076011
The present work presents a 2D and 3D modeling of the neutral gas flow in the sub-divertor region of the W7-X. The investigations have been done using the DIVGAS code. The complex 2D and 3D geometries of the divertor components in the sub-divertor region have been considered and the Standard and High-Iota magnetic configurations have been numerically simulated. The main objective of this study is to investigate the dynamics of neutral particles in the sub-divertor region including the effects due to geometry and toroidal and poloidal leakages located at the divertor targets and baffles on the achieved pumping efficiency. A sensitivity analysis has been performed for the effect of various geometrical and flow parameters on the pumping performance, under different plasma scenarios. The considered incoming fluxes in the sub-divertor range between 1020 to 1022 (H2 s−1). The main conclusions, which can be extracted from the present numerical analysis could be summarized as follows; a large fraction of incoming neutral particle flux i.e. >70% on the low iota side and >40% for the high iota side is leaked back to the main divertor region, while higher incoming neutral fluxes facilitate the increase of the pumped flux as well as the decrease of the outflux. It has been estimated that a small fraction ∼3%–4% of the incoming neutral flux is being pumped via the turbo-molecular pumps. The closure of the toroidal leakages as well as the inclination of the pumping gap panel by 9o facilitate the increase of the pumped flux, but considering the all the engineering constraints, the latter option seems to be more easy to be implemented. For low incoming neutral fluxes (∼1020 H2 s−1) and for the case of AEH section, free molecular flow conditions are estimated and therefore neutral-neutral collisions could be neglected. For higher incoming neutral fluxes and for both AEH and AEP sections neutral-neutral collisions play a significant role in the flow establishment. A comparison with available experimental measurements of the neutral pressure in the sub-divertor has been performed for Standard and High-Iota plasma discharges. The 3D DIVGAS simulations predict qualitatively the experimental data with relative deviation between 25 and 63%. All the above numerical findings will actively support the optimization of the W7-X particle exhaust, in view of the experimental campaign OP2.
P. Liu et al 2024 Nucl. Fusion 64 076007
This paper reports global nonlinear gyrokinetic simulations that couple meso-scale reversed shear Alfvén eigenmodes (RSAEs) driven by energetic particles (EPs) and ion temperature gradient (ITG) microturbulence driven by thermal plasma, using equilibrium and profiles from DIII-D discharge #159243. In simulations focusing only on the ITG, electrostatic ITG drives a huge thermal ion heat transport, which is reduced by a factor of to a level close to the experimental value in electromagnetic simulation due to finite effect. In the simulations coupling the RSAE and ITG, ITG can scatter the resonant EP nonlinearly trapped by the RSAE and damp the zonal flows generated by the RSAE. The regulation of the RSAE by the ITG greatly reduces the initial saturation amplitude of the RSAE but increases the RSAE amplitude and associated EP transport to experimental levels in the quasi-steady state. The RSAE effects on the ITG, specifically the stronger zonal flows generated by the RSAE and the RSAE frequency modulation of the ITG-induced thermal ion heat transport, in turn, leads to a reduction of the thermal ion heat transport by more than a factor of . For a stronger background ITG, the regulation of the RSAE by the ITG is stronger, while the RSAE effects on the ITG are weaker. This work highlights the importance of cross-scale coupling in the dynamics of the AE turbulence and EP transport.
Xingquan Wu (伍兴权) et al 2024 Nucl. Fusion 64 076008
The magnetic coherent modes (MCM) with toroidal mode number n about 1 (Chen R. et al 2018 Nucl. Fusion58 112004) frequently appear in the edge pedestal of high-confinement tokamak plasmas on EAST in the absence of energetic particles. Although these modes are experimentally compatible with the steady-state operation of the pedestal, the driving mechanism without energetic particles of MCM is a long-standing mystery. To reveal the excitation mechanism, a fluid-drift kinetic hybrid local linear model has been developed. It is found that MCM is a new Alfvén eigenmode with a gap frequency much lower than the ideal Toroidal Alfvén Eigenmodes (TAEs) with two significant properties: (1) due to the unique steep pressure gradient in the pedestal region, the diamagnetic frequency becomes comparable to the ideal TAE frequency, which makes the Alfvén continuum in this region move significantly in the ion diamagnetic direction and form a gap of lower frequency; (2) due to the bounce frequencies of thermal electrons becoming also comparable to the ideal TAE frequency in the pedestal region, the free energy of the pressure gradient can be fed into the MCM through the thermal electron bounce resonance excitation, which is essentially the coupling between the shifted TAEs and low- trapped electron modes. The low- MCM is proved to be a shifted Alfvén gap mode in the pedestal region, which is anticipated to exist in low collisional plasmas of future fusion reactors. It is of great significance to carry out relevant physical model research to enhance the understanding of pedestal physics.
Bo S. Schmidt et al 2024 Nucl. Fusion 64 076009
Fast-ion loss detectors (FILDs) are crucial for analyzing fast-ion dynamics in magnetically confined fusion plasmas. A core challenge is to derive an accurate ion velocity distribution, requiring treatment of thousands of remapped camera frames for a full discharge. The ill-posed nature of this task necessitates regularization with a well-chosen regularization parameter and computationally efficient methods. In this work, we introduce the 'resolution principle,' a novel criterion for selecting the optimal regularization parameter, providing a distinction between genuine features and artefacts smaller than the diagnostic resolution in the reconstruction, thereby preventing misinterpretations. This principle, coupled with three iterative reconstruction techniques—Kaczmarz's method, coordinate descent, and Cimmino's method—demonstrates enhanced reconstruction capabilities compared to conventional methods like Tikhonov regularization. Utilizing these techniques allows rapid processing of measurements from full discharges, removing the computational bottleneck and facilitating between-discharge reconstructions. By reconstructing 6000 camera frames from an ELMy H-mode discharge at ASDEX Upgrade, we capture the temporal evolution of gyroradii and pitch angles, unveiling a direct correlation between pitch-angle behavior and changes in the toroidal magnetic field for a specific subset of lost ions accelerated by edge-localized modes (ELMs) to energies approximately twice that of the injection energy.
K. Ogawa et al 2024 Nucl. Fusion 64 076010
Energetic ion anisotropy was observed by tangential sightline compact neutron energy spectrometers (CNESs) in tangential neutral beam heated deuterium plasmas in Large Helical Device. Significant upper and lower energy shifts in D–D neutron energy from 2.45 MeV were measured according to the beam ion injection directions and CNES sightline using a conventional liquid scintillation detector with the unfolding technique and a novel Cs2LiYCl6:Ce with a 7Li-enrichment (CLYC7) scintillation detector without unfolding. The observed neutron energy spectrum was compared with that predicted by a numerical simulation based on orbit following models. Numerical simulation revealed that the Doppler shift in D–D neutron energy results from energetic ion anisotropy.
Niek den Harder et al 2024 Nucl. Fusion
A low beamlet divergence is crucial for the efficiency of the ITER-NBI systems, since it affects the transmission of the beam through the duct. There is a requirement of 7 mrad e-folding divergence for the ITER Heating Neutral Beam. Significantly higher divergences (10-15 mrad) have been observed in RF-source based experiments albeit at low beam energy. This could be the consequence of a broad perpendicular velocity distribution of the H-/D- particles before extraction. This paper explores this hypothesis and its implications for ITER. To estimate H-/D- perpendicular temperatures in the RF-driven BATMAN Upgrade test facility, spatially resolved measurements of the beam power density are compared with IBSimu calculations. The estimated perpendicular temperatures show a strong dependence on the source filling pressure, decreasing from approximately 4 eV at 0.3 Pa to 2 eV at 0.4 Pa. Ion-optics calculations of the ITER-HNB grid system are performed to evaluate whether the temperatures estimated in the BATMAN Upgrade test facility are tolerable in view of beam-grid interaction and beamline transmission. The beamline transmission is fairly insensitive to the perpendicular temperature, but the heat loads at the downstream grids increase with the perpendicular temperature.
Heinke Frerichs et al 2024 Nucl. Fusion
A subspace of resonant magnetic perturbation (RMP) configurations for edge localized mode (ELM) suppression is predicted for H-mode burning plasmas at 15 MA current and 5.3 T magnetic field in ITER. Perturbation of the core plasma can be reduced by a factor of 2 for equivalent edge stability proxies, while the perturbed plasma boundary geometry remains mostly resilient. The striation width of perturbed field lines connecting from the main plasma (normalized poloidal flux < 1) to the divertor targets is found to be significantly larger than the expected heat load width in the absence of RMPs. This facilitates heat load spreading with peak values at an acceptable level below 10 MW m−2 on the outer target already at moderate gas fueling and low Ne seeding for additional radiative dissipation of the 100 MW of power into the scrape-off layer (SOL). On the inner target, however, re-attachment is predicted away from the equilibrium strike point due to increased upstream heat flux, higher downstream temperature and less efficient impurity radiation.
D. Terranova et al 2024 Nucl. Fusion 64 076003
The RFX-mod2 installation is planned to be completed by 2024 and the start of operations is expected in 2025. The high flexibility of the machine (already tested in the previous RFX-mod experiment) allows operation in Reversed Field Pinch and tokamak configuration as well as ultra-low q pulses. In this work we present predictive analysis on transport, performances and plasma control in RFX-mod2 in view of the first experimental campaigns.
D.P. Schissel et al 2024 Nucl. Fusion 64 076004
Full remote scientific operation of the DIII-D National Fusion Facility is now possible through significant advances in the computer science hardware and software infrastructure made over the last decade. Capabilities around information visualization, data movement, and communication have all been enhanced. The level of capability deployed to remotely operate DIII-D required an infrastructure advancement over what had previously been achieved in the fusion community. The large quantity of real-time data that is automatically displayed on DIII-D's control room screens can now be visualized by remote participants via web-based applications. New audio/video solutions using the VoIP and instant messaging application Discord have been implemented to mimic the dynamic and ad-hoc scientific conversations that are critical in successfully operating an experimental campaign. Discord's ability for a user to rapidly move between audio channels, text with images, and share screens is a significant enhancement over traditional videoconferencing tools. In addition, multiple combinations of broadcast audio are made available via a web-based application to allow remote participants to simultaneously listen to general announcements/sounds while conducting their own specific conversations. Secure methodologies have been put into place to allow remote control of hardware including DIII-D's plasma control system application. Secure methods also included the ability of the on-site team to closely coordinate their work with remote team members which has been enhanced through extensions to the wireless network and the use of tablet computers for audio/video/screen sharing. However, no amount of software can fully replace the need for 'hands on hardware.' This infrastructure was severely stress tested during the COVID-19 pandemic where occupancy of the DIII-D control room was restricted. Operational efficiency during the pandemic, measured in discharges per hour, remained high (3.8 ± 0.8) compared to values obtained pre-pandemic (3.7 ± 0.8).
A. Lasa et al 2024 Nucl. Fusion 64 076006
Integrated modeling of plasma-surface interactions provides a comprehensive and self-consistent description of the system, moving the field closer to developing predictive and design capabilities for plasma facing components. One such workflow, including descriptions for the scrape-off-layer plasma, ion-surface interactions and the sub-surface evolution, was previously used to address steady-state scenarios and has recently been extended to incorporate time-dependence and two-way information flow. The new model can address dynamic recycling in transient scenarios, such as the application presented in this paper: the evolution of W samples pre-damaged by helium and exposed to ELMy H-mode plasmas in the DIII-D DiMES. A first set of simulations explored the effect of ELM frequency. This study was discussed in detail in this conference's proceedings and is summarized here. The 2nd set of simulations, which is the focus of this paper, explores the effect of code-coupling frequency. These simulations include initial SOLPS solutions converged to the inter-ELM state, ion impact energy (Ein) and angles (Ain) calculated by hPIC2, and an improved heat transfer description in Xolotl. The model predicts increases in particle fluxes and decreases in heat fluxes by 10%–20% with the coupling time-step. Compared with the first set of simulations, the less shallow impact angle leads to smaller reflection rates and significant D implantation. The higher fraction of implanted flux (and deeper), in particular during ELMs, increases the accumulated D content in the W near-surface region. Future expansion of the workflow includes coupling to hPIC2 and GITR to ensure accurate descriptions of Ein and Ain, and W impurity transport.