SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
This is the final technical report for the activities acrried out in support of the U.S. Burning Plasma Organization (USBPO). The report contains an overview of the administrative and organisational tasks that were conducted in the role of the Deputy Director of the USBPO.
The major goal of this project is to deliver a suite of validated energetic particle (EP) modeling tools for time-dependent integrated tokamak simulations (with the code TRANSP) that can be used for modeling plasma discharges featuring a substantial EP content, such as burning plasmas. During FY20, issues were resolved with an energetic particle transport model called TGLF-EP, with improvements made to the TGLF-EP workflow, generation of input files, and interpretation of runs on the NERSC supercomputing cluster. The recipient also gave input on simulations with the code TRANSP to reproduce the reference JET discharge, and assessed the accuracy and uncertainties of the simulated B field, current evolution, thermal plasma profiles, stored energy, etc. Input was also given on interpretation of fast-ion instability calculations and recommendations for modifications to JET afterglow discharge for the upcoming DT campaign.
In future tokamaks, the control of burning plasmas will require careful regulation of the plasma density and temperature. Along with the design of effective burn-control systems, understanding how the fusion power varies in the density-temperature space is vital for the operation of fusion power plants. Here in this work, the steady-state operational space of ITER is studied using a control-oriented core-plasma model coupled to a two-point model of the scrape-off-layer (SOL) and divertor regions. The two models are coupled through the exchange of input-output parameters. The deuterium and tritium recycling from the wall are output parameters of the SOL-divertor model that are used as input parameters in the core-plasma density balance. Furthermore, the separatrix temperature, which is an output parameter of the SOL-divertor model, is incorporated into the radial core-plasma temperature profiles. Therefore, the temperature-dependent power balance of the plasma core is intimately linked to the SOL-divertor model. Both the power entering the SOL from the core, as determined by the core-plasma power balance, and the separatrix density, as dictated by the core-plasma density balance, are input parameters to the SOL-divertor model. They are control knobs in the SOL-divertor model that can be regulated using the core-plasma actuators: auxiliary power and pellet injection. There are various operational limitations, such as the saturation of the aforementioned actuators, that will prevent ITER from accessing certain high-fusion plasma regimes. The achievable tritium concentration in the fueling lines and the maximum sustainable heat load on the divertor will impose further restrictions. By accounting for these limitations, the ITER operational space is computed based on the coupled core-SOL-divertor model and visualized using Plasma Operation Contour (POPCON) plots that map performance metrics, such as the fusion to auxiliary power ratio, over the density-temperature space. Comparisons are drawn between plasmas with different recycling, confinement, and SOL-divertor conditions.
Divertor design and choice of plasma-facing materials (PFM) will be essential to the success of next-generation fusion reactors as they operate under more powerful scenarios. Understanding and controlling interactions between the plasma and PFM is essential to making these choices. Within these plasma–material interactions and especially in tungsten (W), the interplay between the most abundant plasma species (hydrogen isotopes and helium, He) with the wall material alters fuel retention. However, this interplay is yet to be sufficiently understood to confidently project fuel retention levels to future fusion devices. The paper presents a series of integrated simulations of fusion plasmas and their interaction with tungsten. Specifically, this study assesses the impact of He plasma pre-exposure on hydrogenic species retention during 100 s of burning plasma operations (BPO) in ITER. Multiple pre-exposure scenarios are considered, including sub-surface damage resulting from exposures in the linear device PISCES and from early ITER He-operation. The predictions from these consecutive He-BPO exposures show that fuel content and spatial distribution in the material are largely determined by the He-induced damage, as manifest in: (i) changes in surface temperature expected during BPO have little effect on fuel retention in the presence of He-induced damage; (ii) gas content stabilizes quickly in substrates pre-exposed in PISCES, at levels set by the concentration of pre-existing vacancies, while it continues to increase in substrates initially pristine or pre-exposed to ITER He plasmas; (iii) the presence of He and He–V clusters in the near-surface region locally increases hydrogenic retention, but decreases its permeation; this results in hydrogenic species that remain closer to the surface in pre-damaged substrates, while the bulk content is higher for initially pristine cases. In summary, the interaction and binding of D and T with the pre-existing He–V clusters modifies retention and permeation of hydrogen species during ITER BPO.
In the push to higher performance fusion plasmas, two critical quantities to diagnose are α-heat deposition that can improve and impurities mixed into the plasma that can limit performance. In high-density, highly collisional inertial confinement fusion burning plasmas, there is a significant probability that deuterium–tritium (DT) fusion products, 14.1 MeV neutrons and 3.5 MeV α-particles, will collide with and deposit energy onto (“up-scatter”) surrounding deuterium and tritium fuel ions. These up-scattered D and T ions can then undergo fusion while in-flight and produce an up-scattered neutron (15–30 MeV). These reaction-in-flight (RIF) neutrons can then be uniquely identified in the measured neutron energy spectrum. Further, the magnitude, shape, and relative size of this spectral feature can inform models of stopping-power in the DT plasma and hence is directly proportional to α-heat deposition. In addition, the RIF spectrum can be related to mix into the burning fuel, particularly relevant for high-Z shell and other emerging National Ignition Facility platforms. The neutron time-of-flight diagnostic upgrades needed to obtain this small signal, ~10 –5 times the primary DT neutron peak, will be discussed. Results from several gain > 1 implosions will be shown and compared to previous RIF spectra. Finally, comparisons of experimental data to a simplified computational model will be made.
The EFIT-AI project is creating a modern advanced equilibrium reconstruction code suitable for tokamak experiments of burning plasmas. EFIT [1,2] was the first and is the most extensively used equilibrium reconstruction code in the world. This project builds on the production-level experience and adds key elements as follows. 1. A Model Order Reduction (MOR) version of the two-dimensional (2D) Grad-Shafranov equation solver (EFIT-MORNN) using physics-informed neural networks. 2. Improved optimization and data analysis capabilities using a Bayesian framework enhanced with machine learning. 3. A MOR version of the three-dimensional (3D) perturbed equilibrium reconstruction tool.
We review the physics of energetic particles (EPs) in magnetically confined burning fusion plasmas with focus on advances since the last update of the ITER Physics Basis (Fasoli et al 2007 Nucl. Fusion 47 S264). Topics include basic EP physics, EP generation, diagnostics of EPs and instabilities, the interaction of EPs and thermal plasma instabilities, EP-driven instabilities, energetic particle modes (EPMs), and turbulence, linear and nonlinear stability and simulation of EP-driven instabilities and EPMs, 3D effects, scenario optimization strategies based on EP phase-space control, EPs in reduced field scenarios in ITER before DT, and the physics of runaway electrons. We describe the simulation and modeling of EPs in fusion plasmas, including instability drive and damping as well as EP transport, with a range of approaches from first-principles to reduced models, including gyrokinetic simulations, kinetic-MHD models, gyrofluid models, reduced models, and semi-analytical approaches.
The initiation, termination, and control of internal transport barriers associated with E × B flow shear near local minima of magnetic shear are examined for burning plasmas to determine if the positive feedback loops between profiles, instability, transport, and flow shear operate in regimes with fusion self-heating. A five-field transport model for the evolution of profiles of density, ion and electron temperature, ion and electron fluctuations, and radial electric field is utilized to examine the efficacy of controls associated with external inputs of heat and particles, including neutral beam injection, RF, pellets, and gas puffing. The response of the plasma to these inputs is studied in the presence of self-heating. The latter is affected by the external inputs and their modification of profiles and is, therefore, not an external control. Provided sufficient external power is applied, internal transport barriers can be created and controlled, both in ion and electron channels. Barrier control is sensitive to the locations of power deposition and pellet ablation, as well as temporal sequencing of external inputs.
Abstract 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. Perturbations to 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 striated domain 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.
We are reporting on the modeling of fusion alpha particle transport in the planned ARC fusion device being designed by the CFS (Commonwealth Fusion Systems: https://cfs.energy). The ARC tokamak is designed to operate in a burning-plasma regime characterized by a substantial population of fusion-born alpha particles. Alfvén eigenmode (AE) stability is assessed both analytically and numerically, incorporating alpha-particle drive, ion Landau, and radiative damping from thermal species and collisional damping from trapped electrons. Regions of unstable and near-threshold AE activity are mapped across ARC’s operational parameter space. Linear stability analysis with NOVA indicates multiple, often marginally unstable AEs, extending to toroidal mode numbers up to n= 30. The present report focuses on the ARC flat-top operating point prior to the sawtooth event. Alpha-particle transport on timescales exceeding the neoclassical slowing-down time is assessed using the NUBEAM module [1][2] of the TRANSP code [3], employing transport coefficients derived from the RBQ quasilinear modeling (cf. Appendix B). These global simulations identify favorable and unfavorable operating regimes with respect to alpha confinement, pressure redistribution, and overall alpha-heating efficiency. We also evaluate additional transport mechanisms—including neoclassical tearing mode (TM)–induced stochasticity, sawtooth-driven redistribution, and toroidal-field ripple using the kick model (cf. Appendix C) which makes use of the guiding-center code ORBIT, see Section 5. The kick model is integrated into TRANSP to enable self-consistent predictions of alpha-driven current formation and sustainment within the ARC scenario. Sensitivity scans are performed over the mode frequency, rational-surface alignment, island width, mode amplitude, and proximity of the limiter to the plasma. Our study provides an initial, physics-based guidance for machine design, operational planning, and equilibrium control, ensuring adequate alpha confinement and robust self-heating performance in ARC. Our simulations mostly targeted worst case scenarios, e.g. for TMs and sawteeth. Overall, we expect benign effects for the ARC scenario investigated in this work on fusion alpha confinement and losses in the presence of AEs, tearing modes and sawteeth. This report addresses three thrusts identified at the outset. The first thrust focuses on analytic estimates of the parametric dependencies of EP relaxation based on local AE stability simulations (Section 3). The second thrust involves global evaluations of AE stability using the NOVA, RBQ, and NUBEAM codes (Section 4). Finally, we investigate alpha-particle transport driven by low-frequency instabilities associated with sawteeth and tearing modes (Section 5).
Abstract This work presents the PORTALS framework (Rodriguez-Fernandez et al 2022 Nucl. Fusion 62 076036), which leverages surrogate modeling and optimization techniques to enable the prediction of core plasma profiles and performance with nonlinear gyrokinetic simulations at significantly reduced cost, with no loss of accuracy. The efficiency of PORTALS is benchmarked against standard methods, and its full potential is demonstrated on a unique, simultaneous 5-channel (electron temperature, ion temperature, electron density, impurity density and angular rotation) prediction of steady-state profiles in a DIII-D ITER Similar Shape plasma with GPU-accelerated, nonlinear CGYRO (Candy et al 2016 J. Comput. Phys. 324 73–93). This paper also provides general guidelines for accurate performance predictions in burning plasmas and the impact of transport modeling in fusion pilot plants studies.
Here we report the progress of the development and optimization of operational scenarios for ITER and beyond, focusing upon baseline, hybrid, and steady-state scenarios since 2007. This includes advancements made by the integrated operation scenarios (IOS) topical group of the international tokamak physical activity as well as contributions from the broader tokamak community. The key area of research involves developing IOSs that encompass tokamak physics, operation, and technology by utilizing integrated modeling and control strategies. This requires leveraging available actuators to simultaneously control plasma position and shape, MHD activities that could lead to disruptions, transport, plasma-wall interaction and power exhaust, fuel cycle, fusion burn, and tritium breeding. The control extends from the plasma initiation phase, through the current ramp-up, flattop, start and end of the fusion burn, and current ramp-down, to the plasma termination phase. A review of the currently developed scenarios and modeling is provided in terms of (i) optimizing plasma initiation in ITER, (ii) preparing for the low activation phase to fully commission all tokamak systems and establish and validate physics and scenario conditions in preparation for deuterim-tritium (DT) operation, (iii) developing and preparing baseline and hybrid scenarios to demonstrate the feasibility of achieving these regimes within device constraints, (iv) exploring steady-state scenarios to meet ITER’s steady-state goals, (v) evaluating and preparing actuators for ITER, (vi) developing integrated control solutions using shared actuators. The most notable achievements include; (i) the development of ITER demonstration discharges by matching various dimensionless parameters, (ii) the development of scenarios in an ITER-like tungsten environment and DT operation, and (iii) the development of scenarios in superconducting tokamaks, enabling long-pulse operations with similar coil constraints to ITER. Along with these significant achievements, outstanding issues and recommendations for further research and development are provided. Importantly, this study goes beyond simply updating the ITER Physics Basis; it carries profound implications for the broader field of burning plasma research, offering valuable insights and guidance for the next generation of fusion experiments and devices.
At the National Ignition Facility, inertial confinement fusion experiments aim to burn and ignite a hydrogen plasma to generate a net source of energy through the fusion of deuterium and tritium ions. The energy deposited by α-particles released from the deuterium–tritium fusion reaction plays the central role in heating the fuel to achieve a sustained thermonuclear burn. In the hydrodynamic picture, α-heating increases the temperature of the plasma, leading to increased reactivity because the mean ion kinetic energy increases. Therefore, the ion temperature is related to the mean ion kinetic energy. In this report we use the moments of the neutron spectrum to study the relationship between the ion temperature (measured by the variance in the neutron kinetic energy spectrum) and the ion mean kinetic energy (measured by the shift in the mean neutron energy). We observe a departure from the relationship expected for plasmas where the ion relative kinetic energy distribution is Maxwell–Boltzmann, when the plasma begins to burn. Understanding the cause of this departure from hydrodynamic behaviour could be important for achieving robust and reproducible ignition.
A recent experiment conducted on the National Ignition Facility (NIF) described in the study by Abu-Shawareb et al. achieved a fusion yield output of 1.3 MJ from ~ 220 kJ of x-ray energy absorbed by the capsule, demonstrating remarkable progress in the field of laser driven inertial confinement fusion. In the study by A. R. Christopherson [“Effects of charged particle heating on the hydrodynamics of inertially confined plasmas,” Ph.D. thesis (2020)], the plasma conditions needed to claim the onset of ignition and burn propagation were outlined and multiple criterion were provided to assess progress in inertial fusion experiments. In this work, we modify the metrics from A. R. Christopherson to accurately calculate performance metrics for indirect-drive experiments on the NIF. We also show that performance metric trends observed in NIF data are consistent with theory and simulations. This analysis indicates that all the identified criterion for ignition and burn propagation have been exceeded by experiment 210 808.
Direct evidence of inertially confined fusion ignition appears in the abrupt temperature increase and consequent rapid increase in the thermonuclear burn rate as seen in the reaction history. The Gamma Reaction History (GRH) and Gas Cherenkov Detector (GCD) diagnostics are γ-based Cherenkov detectors that provide high quality measurements of deuterium–tritium fusion γ ray production and are, thus, capable of monitoring the thermonuclear burn rate. Temporal shifts in both peak burn time and burn width have been observed during recent high-yield shots (yields greater than 1017 neutrons) and are essential diagnostic signatures of the ignition process. While the current GRH and GCD detectors are fast enough to sense the changes of reaction history due to alpha heating, they do not have enough dynamic range to capture the onset of alpha heating. The next generation of instrumentation, GRH-15m, is proposed to increase the yield-rate coverage to measure the onset of alpha-heating.
This report includes the first set of MOSS data collected in the context of fusion-relevant plasma-material interactions. The data collected show that the surface curvature and stress irreversibly change during low-energy deuterium ion irradiation. This holds for both tungsten thin films, bulk tungsten, and single-crystal silicon. Meanwhile, argon irradiation did not produce this same effect. The surface morphology of the exposed samples was virtually identical regardless of gas species. There were, however, some surface chemistry differences between argon-exposed and deuterium-exposed samples, mainly concerning the oxide states of the W4f XPS signal. The changes in surface chemistry could be due to hydrogen dynamics on the tungsten surfaces, though further investigation will be required with in-situ irradiation and analysis. The changes in surface stress could also not be explained by temperature alone, as the temperature reached during the ion irradiations was too low to induce non-elastic strain. With this in mind, deuterium retention in the materials could change the surface curvature and stress. To confirm this, further work will be needed to quantify the change in curvature and relate it to ion flux, fluence, and deuterium retention. Unfortunately, when this work was conducted, the experimental setup for thermal desorption spectroscopy was not operating correctly, and that work could not be included in this dissertation. Despite experimental difficulties, this work shows the first results of a promising technique to study deuterium dynamics in tungsten-based reconstituted thin films.
The ability to provide fusion burn control without requiring physical access through the first wall and fuel breeding blankets, would be vital for any future, magnetically confined fusion power reactor. A multi-sensor, fusion fuel cycle exhaust, neutral gas analysis system on JET, capable of delivering real time data, and accessing only the sub-divertor region, provides an excellent example of such capability. Optimized for and operated during the deuterium–tritium experimental campaigns 2 and 3 (DTE2, DTE3), it is proving valuable for planning to explore fusion reactor burn control in ITER with a comparable diagnostic system called the Diagnostic Residual Gas Analyzer (DRGA). This paper aims to show feasibility of developing model-based controllers for ITER and next generation, reactor-relevant devices, by building both on the empirical experience in JET-DTE2, and on the already emerging experience on developing such models specifically for ITER. The paper begins with a specific use-case from JET-DTE2, pertaining to the observed sensitivity of the fusion neutron yield on the concentration of isotopic helium-3 ( 3 He), with data from one of the high-performance DT shots exhibited with emphasis on the 3 He measurement via the sub-divertor. Then, a first model is developed and then explored with simulations that aim to discover how well the controllers in the model react to either insufficient levels of 3 He or excessive levels of 3 He. The simulations then explore potential impact from a delay in the measurement (or the response) that would be comparable to the ∼1 s, conductance limited response for the ITER DRGA system, currently in its final design. The simulations show that control is feasible, and that its effectiveness is not significantly impacted by such delay.
In indirect-drive implosions, the final core hot spot energy and pressure and, hence, neutron yield attainable in 1D increase with increasing laser peak power and, hence, radiation drive temperature at the fixed capsule and Hohlraum size. Here we present simple analytic scalings validated by 1D simulations that quantify the improvement in performance and use this to explain existing data and simulation trends. Extrapolating to the 500 TW National Ignition Facility peak power limit in a low gas-fill 5.4 mm diameter Hohlraum based on existing high adiabat implosion data at 400 TW, 1.3 MJ and 1 × 10 16 yield, we find that a 2–3 × 10 17 yield (0.5–0.7 MJ) is plausible using only 1.8 MJ of laser energy. Based on existing data varying deuterium–tritium (DT) fuel thickness and dopant areal density, further improvements should be possible by increasing DT fuel areal density, and hence confinement time and yield amplification.