Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fusion Reactor”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Generation of shear flows induced by AE / EPM in LHD plasma

The generation of shear flows (SFs) by Alfven Eigenmodes (AEs) and energetic particle modes (EPMs) have important effects on the operation of future nuclear fusion reactors, because SFs regulate the saturation of the AEs/EPMs, the transport of EPs and thermal plasma, as well as the formation of transport barriers among other consequences. The aim of this study is the analysis of SFs generation during the saturation phase of AEs and EPMs in LHD plasma. Experiments performed in the 23rd and 24th LHD experimental campaigns are dedicated to explore the destabilization of AEs/EPMs in discharges with different heating patterns, thermal plasma and magnetic field configurations. In particular, the shots 176490 and 179697 show the destabilization of MHD bursts and energetic-ion-driven resistive interchange modes (EIC), respectively. Charge exchange spectroscopy measurements in both discharges indicate that the generation of SFs by AE/EPM is uncorrelated with the perturbation induced by the neutral beam injector (NBI). Nonlinear simulations performed using the gyro-fluid code FAR3d show the generation of zonal structures, especially SFs, induced during the saturation phase of Toroidal Alfven Eigenmodes (TAEs) triggered in the MHD burst as well as by the 1/1 EIC in the bursting phase. The simulations indicate that SFs are caused by the radial electric fields powered by energy transfers from the unstable AE/EPM towards the thermal plasma. The strongest SFs are measured during the EIC bursting phase once the 1/1 EPM overlaps with nearby resonances at the plasma periphery. Likewise, the largest SFs during the MHD burst are observed once TAEs radially overlap in the inner-middle plasma region.

AE↗
DOE OSTI2025

Turbulence suppression at extreme plasma densities on DIII-D and EAST

Recent high-poloidal-beta (high-βP) experiments on DIII-D and EAST have made coordinated breakthroughs for high confinement quality at high density near the Greenwald limit. Density gradient amplification of turbulence suppression at high βP can explain both of these achievements. Experiments on DIII-D have achieved Greenwald fraction (fGr = line-averaged density/Greenwald density) above 1 simultaneously with normalized energy confinement (H98y2) around 1.5, as required in fusion reactor designs but never before verified in tokamak experiments with the divertor configuration. A synergy between increased H98y2 and fGr is observed with strong gas puffing, due to the build-up of an internal transport barrier at large radius in the temperature and density channels. Transport simulations reveal that the favorable trend of reduced turbulent energy transport at higher density is only expected when increasing the density gradient at high local safety factor and high β, thus at high βP to ensure strong α-stabilization. These conditions are crucial to many conceptual designs for steady-state reactors. New experiments on EAST have nearly doubled the ion temperature at fGr ∼ 0.9, consistent with predict-first modeling results based on the same physics revealed from the DIII-D analysis. All previous EAST long-pulse H-modes have Ti ≪ Te near plasma axis. Transport modeling indicates that the profiles are limited by ion-temperature-gradient modes at mid-radius. The modeling also suggested potential solutions, including reducing magnetic shear, enhancing density gradients, and higher impurity concentration. Following this guidance, EAST experiments directly show a strong enhancement of Ti achieved with a combination of a second plasma current ramp-up, a density gradient increase, and a Zeff perturbation by a short pulse (100 ms) of impurity injection, as predicted by the earlier modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Techno-Economic Analysis of Gas-Liquid Contactors for Tritium Extraction from Lead-Lithium

To enable a sustainable fuel cycle, any deuterium-tritium fusion reactor must breed its tritium fuel onsite. Lead-lithium (PbLi), a eutectic metal, is a leading liquid breeder material for tritium generation. One challenge with PbLi blanket technology is the extraction of tritium from the molten eutectic. Three technologies are the focus of worldwide research: the vacuum permeator, the vacuum sieve tray, and the gas-liquid contactor (GLC). The present work offers a methodology for designing, sizing, optimizing, and costing a trickle-bed GLC. Here, we analyzed tritium extraction from PbLi using MELODIE experimental data by applying traditional packed bed mass transfer efficiency models along with supplementary models, like film theory. Our analysis revealed that traditional packed bed mass transfer efficiency models do not fit the MELODIE loop experimental data. Moreover, uncertainty in PbLi solubility resulted in a 325-fold increase in required gas flow rates when comparing identical packing heights. The film theory liquid mass transfer coefficient, Delt-Olujic wettability model, and Reiter tritium solubility values fit the MELODIE data best and were used both in the design and to conduct the economic analysis. Techno-economic analysis of the GLC was performed to evaluate three design sizes, all achieving a minimum extraction efficiency of 90 [%] for a total tritium extraction of 31 [kg/yr].

Fusion Fuel Cycle↗
DOE OSTI2023

Influence of fatigue precracking and specimen size on Master Curve fracture toughness measurements of EUROFER97 and F82H steels

EUROFER97 and F82H steels are two leading reduced-activation ferritic-martensitic (RAFM) steels for fusion first wall and blanket applications. Exposure to the harsh environment of fusion reactors can result in severe degradation of fracture toughness. Thus, the post-irradiation evaluation of fracture toughness is critical for understanding the material behavior. Due to the space constraints of irradiation facilities and challenges in controlling a uniform irradiation condition for large size specimens, the development of small specimen test techniques (SSTT) is indispensable to evaluate the performance of irradiated materials. In this study, we evaluated specimen size effects on the Master Curve fracture toughness of EUROFER97 and F82H steels. A wide variety of specimens, including 0.5 T compact tension (C(T)) specimens, 0.16 T mini-compact tension (miniC(T)) specimens, and 1.65 mm miniature bend bar specimens, were tested. The testing methodology was based on the Master Curve method in the ASTM E1921 standard. No specimen size effect was observed in 0.5 T C(T) and 0.16 T miniC(T) specimens on the Master Curve reference temperature T0, while 1.65 mm miniature bend bar specimens yielded a higher T 0Q . A strong effect of fatigue precracking on T 0 for 0.5 T C(T) and 0.16 T miniC(T) specimens was observed, such that testing on specimens with skewed fatigue precrack fronts resulted in lower T 0 than for specimens with ASTM standard qualified straight fatigue precrack fronts. The results highlight the importance of experimental quality control in developing SSTT for Master Curve fracture toughness testing. Lastly, we also evaluated and provided recommendations on the minimum number of specimens needed for each specimen type for yielding reliable T 0Q values.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Optimization of lithium vapor box divertor evaporator location on NSTX-U using SOLPS-ITER

Commercial fusion reactors will be faced with extremely high divertor target heat fluxes that will require mitigation. Simulations of detachment in an NSTX-U scenario projected to have 92 MW/m 2 unmitigated peak target heat flux are presented, which reaches sub-10 MW/m 2 target heat flux using a highly dissipating lithium vapor box divertor design. The lithium vapor box is a detached divertor design which employs lithium vapor evaporation and condensation to contain lithium below the X-point. Previous SOLPS modeling has indicated a lithium vapor box can reduce the heat flux down to 10 MW/m 2 via simultaneous evaporation from the Private Flux Region (PFR) and the Common Flux Region (CFR) sides of the vapor box. It is found here that PFR evaporation has improved access to the separatrix leading to significantly more efficient power dissipation than CFR evaporation. Simulations of target evaporation with an evaporation distribution that is self-consistent with the temperature of a Capillary Porous System with Fast flowing liquid lithium (CPSF) could reach n Li / n e ~ 0.025-0.030 at the Last Closed Flux Surface (LCFS) depending on the liquid metal flow speeds and lithium sputtering yield, while PFR-side evaporation can reach acceptable heat fluxes with n Li / n e ~ 0.038 at the LCFS. However, PFR evaporator performance can be improved if the target is allowed to be hot enough such that it reflects lithium, reaching n Li / n e ~ 0.028 and reducing required lithium evaporation. Ultimately PFR evaporation and target evaporation are found to have similar ability to produce acceptable heat flux solutions with minimal upstream concentration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

MX precipitate behavior in an irradiated advanced Fe-9Cr steel: Self-ion irradiation effects on phase stability

In an effort to optimize Fe-9Cr reduced activation ferritic/martensitic (RAFM) steels and to inform the design and operation of fusion reactors, this work represents the first in a series of cohesive studies dedicated to the evolution of MX-TiC precipitates under accelerated single and dual ion irradiations. This study investigates CNA9, a simplified Fe-9Cr RAFM steel featuring initial MX-TiC precipitate densities of (2.3±0.3)×10²¹ m⁻³. This material was subjected to single self-ion irradiation at damage levels ranging from 1 to 100 displacements per atom (dpa) over a temperature range of 300 to 600°C, with a nominal dose rate of 7×10⁻⁴ dpa/s. Irradiation-induced coarsening was observed, as evidenced by statistically significant increases in mean diameter sizes, at 15 dpa at both 500°C and 600°C, whereas no coarsening was noted at 300°C or 400°C. Further, complete dissolution of precipitates occurred at damage levels of 50 and 100 dpa across the two temperatures tested (300°C and 500°C) while no significant changes were observed at any doses below 15 dpa at 500°C. Experimentally parameterized recoil resolution modeling suggests that the observed radiation stability of MX-TiC precipitates is intricately linked to diffusional changes of solutes resulting from the co-evolution of microstructural features within the experiments. The findings align with current theoretical perspectives on radiation-induced precipitate stability in complex alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

One-dimensional simulation and validation of divertor detachment induced through nitrogen seeding on HL-2A

Divertor detachment is a promising method to solve the power exhaust problem in tokamak devices or even in future magnetic fusion reactors. In this work, a detachment experiment (HL-2A shot #38008) with mixed gas seeding (60% nitrogen and 40% deuterium) is simulated using the SD1D module in BOUT++. In the process from attachment to detachment, the target electron temperature and the target ion saturation current in simulations are found to be consistent with the experimental results measured by Langmuir probes on the target plate. In order to understand the underlying detachment mechanism on HL-2A, this work analyses the role of different particle species in the cases with different seeding rates. It shows that the plasma density varies little and the density of neutrals (D and D 2 ) slightly decreases with the increase of nitrogen seeding rate, such that plasma–neutral interactions cannot effectively reduce plasma energy and plasma momentum in the divertor. The case with a high seeding rate predicts that increasing seeding rate cannot reach a target temperature lower than 2.5 eV, which is the required temperature for strong plasma–molecule interactions. Thus, the plasma–molecule interactions may not be important in the divertor during nitrogen seeding. Here this work also studies the parallel forces including (1) force due to the parallel electric field, (2) friction forces, (3) ion- and electron-thermal forces, and (4) collisional reactions (e.g. charge exchange recombination and ionisation). It is found that the friction force between nitrogen ions and other particle species (primarily D + ) is the dominant force towards the target, while ion- and electron-thermal forces are the dominant force pushing nitrogen ions back to upstream. Parallel forces determine the parallel distribution of nitrogen impurities, and therefore decide the region of nitrogen radiation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Material Plasma Exposure eXperiment (MPEX) High Heat Flux Bellows for Component Alignment in High Microwave Environment

The Material Plasma Exposure eXperiment (MPEX) device is a linear plasma device developed to perform plasma material interaction experiments under the conditions prototypic of a fusion reactor divertor. MPEX has multiple systems that must be precisely aligned to the plasma axis, including an electron cyclotron heating system that emits up to 400 kW of microwave power into the vacuum vessel. Five distinct systems require precise alignment on the MPEX device, thus requiring four bellows, all of which are adjacent to the plasma at a relatively high heat flux of approximately 47 kW/m2 and microwave power regions. The MPEX high heat flux bellows (HHFB) is designed to deliver 6 degrees of freedom positioning. The HHFB includes titanium–zirconium–molybdenum (TZM) inserts that are brazed into a Glidcop AL-15 body using a high-temperature braze alloy, thus blocking direct line of sight to an edge-welded bellows from the plasma and microwave screen to block microwaves from the bellows. Further, a custom ConFlat knife edge is machined into the Glidcop AL-15, so the vacuum flanges do not need a braze or weld joint on the vacuum interface. Fingerstock or copper mesh is used to restrict microwave power from entering the interstitial space between the water-cooled Glidcop AL-15 body and edge-welded bellows. Glidcop AL-15 was selected as the material of choice for the water-cooled body because it can maintain mechanical integrity at elevated temperatures, and it also allows for a high-temperature braze. TZM was selected for its machinability and compatibility with vacuum and plasma requirements. A test article is also being considered for a similarly shaped component, the MPEX limiter, that will demonstrate the integrity of the braze joint under high thermal load. Results from this testing will be extrapolated to deduce the lifetime and integrity of the HHFB design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Innovative laser-based methods for monitoring fuel retention in ITER

This paper addresses the challenge of tritium inventory management in ITER and future fusion reactors, highlighting the importance of accurate tritium measurement and its spatial distribution within the vacuum vessel. Given ITER’s operational constraints, especially the limit on tritium retention, precise measurement is essential for both safety and regulatory compliance. To tackle these questions, the paper presents the T-monitor diagnostic system developed by Forschungszentrum Jülich, which uses Laser-Induced Desorption (LID) in combination with Diagnostic Residual Gas Analysis (DRGA) to measure hydrogen isotope concentrations on the surface of divertor tiles. The system integrates a high-power laser, advanced optical components, and a Fast Scanning Mirror Unit (FSMU) for accurate laser spot positioning with rapid response. Designed to measure in situ tritium retention, the diagnostic provides high-resolution spatial mapping, vital for evaluating detritiation strategies. The laser heating process increases the divertor surface temperature to 1600 K within the laser spot, promoting hydrogen isotope desorption. Accurate measurements require the precise control of laser parameters, including pulse duration and spot size, with a target relative accuracy of 20%. The optical design includes both in-vessel and ex-vessel components, such as durable high-reflectivity mirrors made of gold and copper, selected not only for their infrared performance but also for their transmission of visible wavelengths for observation purposes. To protect optical components from contamination, a pneumatic shutter is used.

Hydrogen isotopes↗
DOE OSTI2022

Benchmarking a 9Cr-2WVTa Reduced Activation Ferritic Martensitic Steel Fabricated via Additive Manufacturing

Reduced activation ferritic (RAF) martensitic steels are promising candidates for the first wall of fusion reactors. However, current manufacturing capabilities call for these components to be made by welding wrought plates. This limits design freedom and necessitates the use of post-weld heat treatments (PWHT) in accordance with the boiler and pressure vessel code. Additive manufacturing (AM) can offer a unique solution to solve this challenge by leveraging the layer-wise deposition strategy to come up with temper bead deposition techniques to eliminate post-processing heat treatments (PPHT). However, it is necessary to benchmark the properties of RAF steels fabricated by AM with their wrought counterparts to identify the process-structure-property correlation, which is the goal of this study. The study demonstrates that while tensile properties at room temperature and high temperatures are satisfactory, the as fabricated and samples after PPHT have significant heterogeneity in tensile elongation. This has been attributed to the presence of discontinuities in the build. The as-fabricated samples have an average tensile strength of 1190 + 12 MPa and an average elongation of 15 + 5% at room temperature and 658 ± 20 MPa ultimate tensile strength (UTS) and 14 ± 7% at 600 °C. After the post-weld heat treatment, mechanical properties decrease to around 600–650 MPa and an elongation between 20–25% at room temperature to 300 MPa UTS and 25–28% elongation at 600 °C. The characterization of microstructures at various length scales demonstrates that the as-fabricated structure has a significant fraction of delta ferrite in a lath martensitic matrix. No precipitates could be identified in the as-fabricated structure. PPHT led to a decrease in the area fraction of delta ferrite and precipitation of M 23 C 6 and MX. Detailed characterization clearly demonstrates that the lack of precipitates in the as-fabricated structure could be due to the slow tempering response of the alloy. Finally, the needs to develop new alloys to achieve the objectives stated above are articulated.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Density drop at the divertor target in the prototype material plasma exposure eXperiment (Proto-MPEX)

The steady-state linear device “Material Plasma Exposure eXperiment” (MPEX) is currently under construction at Oak Ridge National Laboratory with the goal of enabling Plasma-Material Interaction studies at future fusion reactor relevant plasma conditions. In this work, a newly in-house developed hybrid Particle-In-Cell code-PICOS++ is applied to understand the experimental results obtained from the prototype of MPEX referred to as the “Proto-MPEX” during its helicon-only and helicon with ion cyclotron resonance heating (ICRH) experiments. This study explains the physics of the experimentally observed plasma density-drop at the divertor target in Proto-MPEX device during ICRH. In contrast to previous work on ICRH in MPEX, this study demonstrates that the mirror force plays a central role in the Proto-MPEX plasma transport during ICRH, which has new features not previously explored. Force balance analyses reveal that the temperature anisotropy produced by ICRH leads to a significant increase in the mirror force downstream of the resonance where the magnetic field is diverging. This force accelerates ions toward the target and leads to a drop in plasma density to ensure conservation of particle flux. Simulations with ICRH where the magnetic field divergence downstream of the resonance has been removed, do not produce plasma acceleration nor density drop at the target despite efficient ion heating at the resonance. Moreover, simulation results demonstrate that for a given ICRH power, lowering the source rate produces ions with increased perpendicular energy which interact with the mirror force to produce higher plasma acceleration which increases the strength of the density-drop at the target. The strength of the density drop appears to reach an asymptotic limit at a certain threshold ICRH power. Simulations show that this threshold power increases with increasing particle source rate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Turbulence link to L-mode, I-mode, and H-mode confinement in the DIII-D tokamak

Understanding the physics of low-confinement (L-), improved-confinement (I-), and high-confinement (H-) modes is critical for fusion reactors. The finding herein reports observations of two types of turbulence coexisting near the L-mode edge, one magnetohydrodynamic (MHD)-like and another micro-tearing mode (MTM)-like, linked to the H-mode and I-mode confinement in the DIII-D tokamak. Ion-scale magnetic and density turbulence is measured using a Faraday-effect radial-interferometer-polarimeter and beam-emission-spectroscopy (BES). Broadband turbulence spectra of up to ~600 kHz are observed in two discharges where transitions between L-mode, I-mode, and H-mode occurs. Turbulence is found to be inversely correlated with confinement, meaning lower turbulence power at higher confinement. Distinctively, the high-frequency (HF, >~100 kHz) magnetic turbulence power changes by the most (55%) during transitions primarily involving energy confinement change, whereas the low-frequency (LF, <~100 kHz) magnetic and density turbulence power changes by the most (80%) during transitions primarily involving particle confinement change. The LF turbulence amplitude oscillates with and leads to deuterium-alpha emission oscillations before an H-mode. These results imply that HF turbulence mainly affects energy confinement whereas LF turbulence can affect particle confinement. The magnetic and density turbulence exhibits coherence up to 0.6 and cross-phase magnitude close to π/2 in most cases, suggesting they have a common origin in both the LF and HF ranges. BES suggests that LF turbulence resides at the edge (ρ = 0.95) and HF turbulence can be at the outer core (ρ = 0.8) or edge (ρ = 0.95). Comparisons of measurements, theory, and gyrokinetic simulations suggest that HF turbulence is MTM-like in all cases, whereas LF turbulence is more consistent with MHD-like modes and the exact instability might change during transitions—except that a drift-wave origin is possible in a low collisionality H-mode. These results suggest that the H-mode involves suppressed MHD-like turbulence, whereas the I-mode mitigates MTM-like turbulence along with largely unchanged MHD-like turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Taming harmful bursts and heat flux in high-confinement tokamak plasmas

A major challenge in tokamak fusion research is first-wall erosion caused by steady heat loads and sudden energy bursts known as edge-localized modes. Divertor detachment reduces steady-state heat flux, while resonant magnetic perturbations can suppress these instabilities. However, integrating the two has been difficult because they require conflicting operating conditions. Here we demonstrate simultaneous achievement of resonant magnetic perturbations mitigated small edge-localized modes and impurity seeded partial divertor detachment in plasmas with an ITER-similar shape on the DIII-D tokamak. Experiments and simulations show that resonant magnetic perturbations facilitate detachment by redistributing particles, lowering the core density and increasing the scrape-off layer density, thereby reducing the amount of injected gas required. Cooling-gas injection eliminates the secondary heat-flux peak created by three-dimensional magnetic lobes, while edge cooling weakens the plasma response to the applied magnetic fields. These advances illustrate a viable pathway for integrating edge stability control with power exhaust in future fusion reactors.

magnetically confined plasmas↗
DOE OSTI2021

Predicted Heat Flux Performance of Actively Cooled Tungsten-Armored Graphitic Foam Monoblocks

Tungsten (W)–armored graphitic foam monoblocks were developed for applications requiring high-Z plasma-facing material in long-pulse fusion experiments and ultimately deuterium-tritium fusion reactors. The monoblocks are an integrated material system combining the advantages of a chemical vapor deposited (CVD) W coating with a high-conductivity graphitic foam. The W is a high-melting-point, high-Z material with low tritium retention. The graphitic foam coupled to a swirl tube serves as a high-thermal-conductivity heat sink that cannot melt, although it can sublime at much higher temperatures than copper melts. Together, they comprise a robust plasma-facing component (PFC) weighing roughly 5% of an all-W component or 17% of a traditional W-coated copper heat sink.A single-channel mock-up consisting of four graphitic foam monoblocks equipped with a water-cooled swirl tube was fabricated for eventual testing in the 60-kW, EB-60, rastered electron beam at the Applied Research Laboratory of The Pennsylvania State University. Two monoblocks have a thin 50-μm-thick coating of pure W chemically vapor deposited over NbC and pure Nb interlayers. Two others have a 2-mm-thick pure W coating CVD on graphitic monoblocks using the same interlayers. The mock-up will be cooled with available 10 m/s, 0.7 MPa water with a 22°C inlet temperature and subjected to varying uniform heat loads up to 20 MW/m 2 . It is equipped with type-K thermocouples at various depths, and calibrated infrared thermography and spot pyrometry will be used to characterize the heated surface. Real-time water calorimetry will be used to ascertain the absorbed steady-state power and infer the heat flux during testing. Since testing cannot be done under prototypic divertor flow conditions, it is necessary to predict the thermal response of this novel PFC system and investigate the power sharing between radiation and convection at divertor heat flux levels and its inherent ability to avoid critical heat flux. Results are reported for predictions obtained from computational fluid dynamics models up to 30 MW/m 2 of steady-state uniform heat flux. Leading-edge heat loads of 30 MW/m 2 on a 2-mm-wide side strip were also investigated to ascertain if coating delamination is likely.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Understanding the formation of a low-pressure pedestal in the presence of a strong internal transport barrier in DIII-D high β p plasmas

As a promising scenario for fusion reactors, the high poloidal-beta (β P ) scenario is characterized by a strong large radius internal transport barrier (ITB), which significantly enhances the overall confinement quality and the bootstrap current fraction for fully non-inductive operation. It is frequently observed that in the presence of a strong ITB, the pedestal height is lower and is accompanied by small edge localized modes (ELMs), which further improves the compatibility of a high performance core with an edge solution. A mechanism for the formation of the low pedestal is proposed in this paper. It is found that the strong ITB creates an off-axis bootstrap current to clamp the local safety factor q, and thus the magnetic shear in the outer core/pedestal region is increased. Gyrokinetic simulations with the CGYRO code show that the higher magnetic shear brings the experimental profiles into the range where the growth rate of drift-wave instabilities and thus transport is higher, and therefore a lower pedestal gradient is expected. Here, the combination of low pedestal and high magnetic shear further enhances the turbulent transport across the whole pedestal, consistent with power balance analysis. Such a positive feedback mechanism ultimately results in a lower pressure pedestal as observed in experiments. Under such a low pedestal, linear simulations with BOUT++ predict the growth rates of peeling–ballooning modes to be lower across the whole toroidal mode number spectra, and the nonlinear BOUT++ simulation exhibits lower saturated fluctuation intensity as well, consistent with the experimentally observed lower ELM size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Progress of Divertor Heat and Particle Flux Control in EAST for Advanced Steady-State Operation in the Last 10 Years

Active control of the excessively high heat and particle fluxes on the divertor target plates is of fundamental importance to the steady state operation of tokamaks, especially for fusion reactors. A series of experiments have been carried out on this critical issue to relieve the plasma-wall interactions in the experimental advanced superconducting tokamak (EAST) in the last ten years, not only contributing to the long pulse operation of EAST itself, but also providing physical understandings and potential techniques to the next-generation devices like ITER. Here we have characterized the power deposition pattern and broadened the divertor footprint width effectively. The plasma-wetted area is actively handled using either 3-dimentional edge magnetic topology or advanced plasma equilibrium, thereby peak heat flux around the strike point is reduced. Active control of detachment or radiation compatible with core plasma performance has progressed significantly in very recent years, with a series of active feedback control modules developed and utilized successfully, based on the divertor physics advances with both experiments and simulation. The upper divertor of EAST was upgraded from graphite to active water-cooling ITER-like tungsten in 2014, exhibiting much enhanced heat removal capability. As for the particle exhaust including both fueling and impurity particles, in addition to wall conditioning and impurity source control, the efficiency of particle flux exhaust is optimized by making full use of the divertor closure and the plasma drifts in both scrape-off layer and divertor volume. These heat and particle exhaust advances contribute greatly to a series of EAST achievements like H-mode operation over 100 s. A brief near-term plan on the integrated control of divertor plasma-wall interactions in long-time scale will also be introduced, aiming to provide favorable divertor operation solution for ITER and CFETR.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Quantum Information for Fusion Energy Sciences (Final Technical Report)

The simulation of plasma dynamics is a critical area of Fusion Energy Sciences (FES) due to it’s usefulness in predicting, controlling, and confining plasmas in the context of potential fusion reactors. The simulation of plasmas is a computationally difficult problem in both classical and quantum physics, motivating investigation into the potential of quantum computers to simulate these systems. This project took several concrete steps towards this goal by developing tools for improving the control, characterization, and calibration of quantum gates on a superconducting quantum computer, developing error suppression and mitigation tools to reduce errors on the quantum computer, and utilizing these advancements to simulate reduced models of plasma dynamics on the quantum computer. In order to efficiently simulate plasma physics, an optimal control method which synthesizes, directly at the pulse level, any quantum gate on qubit and qutrit systems was developed. Using four superconducting transmon quantum processors at Rigetti and LLNL, it was demonstrated that any arbitrary quantum gate on qubits and qutrits could be implemented with high fidelity, leading to a significantly reduced length of a gate sequence. A problem of interest in FES is the nonlinear optical process of laser pulse compression within a plasma. Since quantum physics is linear, simulating nonlinear operations is not naturally feasible on a quantum computer, however it is possible to simulated a quantized version of the nonlinear process. A quantization approach to convert nonlinear wave-wave interaction problems to Hamiltonian simulation problems was developed and demonstrated using two qubits on a Rigetti device. In this experiment, a number of error suppression and mitigation techniques were investigated to determine how best to utilize the finite quantum resources. This study provides an example of how plasma problems may be solved on near-term, noisy quantum computing platforms and identified a promising set of techniques. Building on the insights of these experiments, the investigation turned to linear electron-plasma wave physics. A connection was identified between a local one-dimensional lattice spin model and linear wave phenomena, allowing a plasma physics problem to be efficiently mapped to the quantum computer. In this framework, reflection and transmission of plasma waves at a sharp boundary was studied, as well as the propagation of waves through an inhomogeneous plasma medium. In addition to the suite of error suppression and mitigation techniques developed, this experiment introduced the use of a digital-analog gate scheme designed to efficiently simulate the plasma Hamiltonian. With hardware available at the conclusion of the project, simulation at the scale of 9 qubits and 15 timesteps (60 entangling layers) was achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A Sparse-Grid Probabilistic Scheme for Approximation of the Runaway Probability of Electrons in Fusion Tokamak Simulation

Runaway electrons (RE) generated during magnetic disruptions present a major threat to the safe operation of plasma nuclear fusion reactors. A critical aspect of understanding RE dynamics is to calculate the runaway probability, i.e., the probability that an electron in the phase space will runaway on, or before, a prescribed time. Such probability can be obtained by solving the adjoint equation of the underlying Fokker-Planck equation that controls the electron dynamics. In this effort, we present a sparse-grid probabilistic scheme for computing the runaway probability. The key ingredient of our approach is to represent the solution of the adjoint equation as a conditional expectation, such that discretizing the differential operator reduces to the approximation of a set of integrals. Adaptive sparse grid interpolation is utilized to approximate the map from the phase space to the runaway probability. The main novelties of this effort are the integration of the sparse-grid method into the probabilistic numerical scheme for computing escape probability, and the application of the proposed method in computing RE probabilities. Two numerical examples are given to illustrate that the proposed method can achieve O(Δt) convergence, and that the local anisotropic adaptive refinement strategy (M. Stoyanov, Adaptive sparse grid construction in a context of local anisotropy and multiple hierarchical parents. In: Sparse Grids and Applications-Miami 2016, Springer, Berlin, 2018, pp. 175–199) can effectively handle the sharp transition layer between the runaway and non-runaway regions.

Yang, Minglei↗