Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma fusion”

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

LLNL FESP Theory Highlights: November 2024

This talk presented a refined understanding of plasma-facing thermionic emitter cathodes by examining the limitations of the conventional space charge-limited model, which overlooks ion trapping dynamics within the virtual cathode. While conventional theory distinguishes between temperature-limited and space charge-limited regimes based on electron current constraints, recent findings reveal that trapped ions can form a quasi-neutral "second plasma," expanding upstream and interacting with the primary plasma. These effects are amplified in non-planar geometries, such as cylindrical or spherical sheaths, where trapped ions create a positive feedback loop with the upstream plasma. To address these phenomena, a new "Aid-and-Compete" model was developed, offering insights into current enhancement and sheath mode transitions relevant for plasma devices, including emissive cathodes in electric propulsion, fusion applications, thin film deposition systems, and plasma diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

STORM: Scrape-off layer turbulence in tokamak fusion reactors

The scrape-off layer of a tokamak fusion reactor carries the plasma exhaust from the hot core plasma to the material surfaces of the reactor vessel. The heat loads imposed by the exhaust are a critical limit on the performance of fusion power plants. Turbulent transport of the plasma regulates the width of the scrape-off layer plasma and must be modelled to understand the intensity of these heat loads. STORM is a plasma turbulence code capable of simulating three dimensional turbulence across the full scrape-off layer of a tokamak fusion reactor, using a drift reduced, collisional fluid model. STORM uses mostly finite difference schemes, with a staggered grid in the direction parallel to the magnetic field. We describe the model, geometry and initialisation options used by STORM, as well as the numerical methods, which are implemented using the BOUT++ plasma simulation framework. BOUT++ has been enhanced alongside the development of STORM, providing better support for staggered grid methods. We summarise these enhancements, including a detailed explanation of the parallel derivative methods, which underwent a major update for version 4 of BOUT++.

BOUT++↗
DOE OSTI2020

The Materials Plasma Exposure eXperiment: Status of the Physics Basis Together With the Conceptual Design and Plans Forward

The Materials Plasma Exposure eXperiment (MPEX) is a linear plasma device that will address the plasma–material interaction (PMI) science for future fusion reactors and will enable testing of plasma-facing components (PFCs). It is designed as a steady-state device eventually delivering an ion fluence of up to 10 31 m -2 to the target. The device will be designed to handle neutron-activated materials. These capabilities should push the technical readiness level of PFCs up to six for some end-of-life aspects. In order to achieve the relevant plasma conditions, as they are expected in future fusion reactor divertors, MPEX will utilize a novel plasma source concept. This plasma source concept uses a high-power helicon (200 kW, 13.56 MHz), an electron cyclotron heating (ECH) system, which will heat electrons via electron Bernstein wave (EBW) heating (about 400 kW, 70 GHz), and an ion cyclotron heating (ICH) system (400 kW, 6–9 MHz). The physics basis for this plasma source concept including the heating physics and transport is based on experiments on Proto-MPEX. An overview of the experimental results and the physics basis will be given. The physics basis directly translates into functional requirements of MPEX and the conceptual design. The status of the conceptual design of MPEX will be shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Interactions Between Microturbulence and Magnetohydrodynamic Instabilities in Magnetized Plasmas (Final Technical Report)

The overall goal of this project was to develop a better understanding between small-scale microturbulence and large-scale magnetohydrodynamic (MHD) instabilities, particularly magnetic islands and tearing modes, in fusion relevant magnetically confined plasmas. This project advanced understanding through numerical simulations of different model equations. Questions addressed included how the magnitudes, spatial structures, and scalings of turbulent forcings and drags impact island dynamics, as well as what the "region of influence" of a tearing mode on the turbulence dynamics. A key finding of the project is that the microturbulence (and associated transport) exhibits spatial asymmetries for sufficiently large islands. The results have improved our fundamental understanding of plasma dynamics, motivated new experimental investigations, and will help inform future fusion predictive modeling efforts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Sensitivities of time-dependent temperature profile predictions for NSTX with the multi-mode model

The Multi-Mode Model (MMM) for turbulent transport was applied to a large set of well-analyzed discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate its sensitivities to a wide range of plasma conditions. MMM calculations were performed for hundreds of milliseconds in each discharge by performing time-dependent predictive simulations with the 1.5D tokamak integrated modeling code TRANSP. A closely related study (Lestz et al 2025 Plasma Phys. Control. Fusion 67 105029) concluded that MMM predicted electron and ion temperature profiles that were in reasonable agreement with NSTX observations, generally outperforming a different reduced transport model, TGLF. This finding motivates the more thorough investigation of the characteristics of the MMM predictions conducted in this work. The simulations with MMM have electron energy transport dominated by electron temperature gradient modes in the examined discharges with relatively low plasma β (ratio of kinetic plasma pressure to magnetic field pressure) and high collisionality, transitioning to a mixture of different modes for higher β and lower collisionality. The thermal ion diffusivity predicted by MMM is much smaller than the neoclassical contribution, in line with previous experimental analysis of NSTX. Nonetheless, the electron and ion temperature profiles are coupled via collisional energy exchange and thus sensitive to which transport channels are predicted. The time-dependent simulations with MMM are robust to the simulation start time, converging to remarkably similar temperature profiles later during the discharge. MMM typically overpredicts confinement relative to NSTX observations, leading to the prediction of overly steep temperature profiles. Plasmas with spatially broader temperature profiles, higher plasma β, and longer energy confinement times tend to be predicted by MMM with better agreement with the experiment. As a result, these findings provide useful context for understanding the regime-dependent tendencies of MMM in anticipation of self-consistent, time-dependent predictive simulations of NSTX-U discharges with these same modeling tools.

MMM↗
DOE OSTI2025

Conference Report on the 8th International Symposium on Liquid Metals Applications for Fusion (ISLA-8)

The International Symposium on Liquid Metals Applications for Fusion (ISLA) aims to assemble scientists and engineers engaged in research on lithium and liquid metal applications for fusion devices, facilitating discussions on recent advancements and challenges in an open forum to support the development of viable fusion reactors. The 8th International Symposium on Liquid Metals Applications for Fusion (ISLA-8) was organized by the Institute of Plasma Physics, Chinese Academy of Sciences, from 8 to 12 September 2024, in Hefei, China. The symposium was attended by over 70 participants, marking one of the highest attendance figures in the series. A total of 68 presentations were delivered, including 59 presented on-site, while the remainder were conducted online. The participants represented 10 countries, namely China, Japan, the Netherlands, Russia, the USA, the Czech Republic, Italy, Thailand, Spain, and Germany. The symposium covered 10 topics, structured into 13 sessions. Additionally, an opening session provided an overview of the current symposium, while a closing session summarized reports from each session chair.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Informing ARC divertor design and plasma facing material selection through integrated modeling

This INFUSE 2023 project between UCLA and Commonwealth Fusion Systems used computer modeling to test whether tungsten materials can survive in CFS's ARC fusion reactor divertor. The team simulated plasma conditions and material responses, finding that tungsten-rhenium alloys resist grain growth better than pure tungsten, and that hydrogen buildup depends more on particle flux than temperature. The work helps CFS design durable plasma-facing components for their fusion power plant.

36 MATERIALS SCIENCE↗
DOE OSTI2025

First result of boronization assisted by the ICWC on EAST with full metal wall

Boron (B), a low-Z (atomic number) material, has been widely utilized in wall conditioning to improve plasma performance in fusion devices. In 2023, boronization was successfully conducted on EAST featuring an ITER-like tungsten divertor and fully metallic first wall. The process employed predischarge coating with carborane (C 2 B 10 H 12 ) as the working material, assisted by ion cyclotron wall conditioning (ICWC). After one time 12 g boronization, it was found the thickness of B film was approximately 120 nm. Post-boronization observations indicated that substantial hydrogen (H) release during initial plasma discharges compared with the consumed W/B wall, attributed to H co-deposition during the ICWC-boronization processing, which led to uncontrollable divertor neutral pressure and plasma density. The H/(H + D) ratio demonstrated a gradual reduction from ∼85% to 30% over more than 1850 s of deuterium plasma, with a cumulative injected energy of 2325 MJ. The B coating significantly enhanced the stored energy in plasma and improved confinement performance. The stored energy in plasma showed an increase of about 20%, primarily due to a reduction in impurity radiation, including oxygen (O) and heavy impurities such as tungsten (W), iron (Fe), and copper (Cu). The effective ion charge (Z eff ) decreased from 2.3 to 2.0. Following ICWC-boronization, the line-integrated radiation profile decreased by nearly 35% in the plasma core, plasma density and electron temperature exhibited an increase of ∼7% and 12% due to enhanced wall fueling and reduced impurity radiation. The lifetime of boronization, as evaluated by the line emissions from boron and other impurity radiation, was about 1700 s of deuterium plasma, with a cumulative injected energy of 2125 MJ on EAST. These findings provide significant insights for evaluating ICWC-boronization applicability in ITER with full W wall structure.

EAST↗
DOE OSTI2021

Phase-space sensitivity (weight functions) of 3 MeV proton diagnostics

Abstract The weight function W ( X ) is a diagnostic sensitivity to phase-space variables X that relates the measured signal C to the distribution function F ( X ) through the equation C = ∫ W ( X ) F ( X ) d X . In the present work, an algorithm to calculate W for a diagnostic that measures 3 MeV protons produced in d(d,p)t fusion reactions between a fast ion and a thermal deuteron is developed. The emitted protons escape the tokamak on curved orbits and are detected. These curved orbits constitute effective diagnostic ‘sightlines’. The presented algorithm accounts for the complications associated with these curved sightlines. An initial calculation of time-reversed proton orbits computes effective solid angles and sightlines for the relevant range of incident proton velocity vectors. These precomputed orbits are then used within the framework of FIDASIM (2020 Plasma Phys. Cont. Fusion 62 105008) to calculate the reactivity averaged over the thermal distribution of the ‘target’ deuterons and the probability that a fast ion of specified energy and pitch has a gyroangle that is consistent with the kinematic equations along each of the sightlines. Comparisons with analytic formulas and with independent calculations for the Mega Amp Spherical Tokamak 3 MeV proton diagnostic verify the algorithm.

Physics↗
DOE OSTI2021

Fiber-optic silicon Fabry–Perot interferometric bolometer with improved detection limit for magnetic confinement fusion

Fiber-optic bolometers (FOBs) intended for plasma radiation measurement in magnetically confined fusion have been previously developed using a silicon pillar that functions as both a Fabry–Perot interferometer (FPI) for temperature measurement and an absorber for the radiation. We report an FOB design that can significantly improve the detection sensitivity over earlier designs by engineering the absorber of the FOB. Our design uses the fact that, compared with a silicon pillar, a gold film with the same x-ray absorption thickness will show a much higher temperature rise from a given power density of the radiation. Therefore, the responsivity of an FOB can be improved by attaching a large gold disk to the silicon FPI as the absorber. We have developed a fabrication method for FOBs of such design and obtained an FOB with a 4-µm-thick, 0.6-mm-diameter gold disk attached to a 200-µm-diameter, 100-µm-thick silicon FPI. We have characterized the noise, responsivity, response time, and noise-equivalent power density (NEPD) and compared these with the earlier design where the absorber is mainly the silicon FPI itself. The experimental result suggests that the FOB with the gold disk achieves a responsivity of ∼2.8 mK/(W/m2) and a noise-equivalent-power-density of 0.13 W/m2, which are, respectively, more than nine times larger and six times smaller compared to the FOB using a previous design. Improved NEPD and good absorption over a broad frequency range will make the FOB more attractive for applications in magnetic-confinement fusion devices.

Uddin, Nezam↗
DOE OSTI2026

Suprathermal electron transport studies in radio-frequency-heated tokamak plasmas

This contribution reports on the study of suprathermal electron radial transport enhancement in the presence of electron-heating radio-frequency waves, such as electron-cyclotron waves or lower-hybrid waves. Since suprathermal electrons emit hard X-rays from Bremsstrahlung, experimental data from hard X-ray measurements are used to follow their dynamics. In particular, electron-cyclotron power modulation experiments performed in TCV show that this transport depends on the total electron-cyclotron power. These results are corroborated by experimentally constrained Fokker-Planck simulations [J. Cazabonne et al, Plasma Phys. Control. Fusion 65, 104001 (2023)]. Similar experiments have been performed in WEST with lower-hybrid waves, and are reported for the first time. Preliminary analyses show potential evidence of suprathermal electron transport, but no clear dependence on wave power nor photon energy at this stage.

Cazabonne, Jean [CEA, IRFM, Saint Paul Lez Durance↗
DOE OSTI2025

Assessing time-dependent temperature profile predictions using reduced transport models for high performing NSTX plasmas

Time-dependent, predictive simulations were performed with the 1.5D tokamak integrated modeling code TRANSP on a large set of well-analyzed, high performing discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate how well modern reduced transport models can reproduce experimentally observed temperature profiles in spherical tokamaks. Overall, it is found that simulations using the Multi-Mode Model (MMM) more consistently agree with the NSTX observations than those using the Trapped Gyro-Landau Fluid (TGLF) model, despite TGLF requiring orders of magnitude greater computational cost. When considering all examined discharges, MMM has median overpredictions of electron temperature (T e ) and ion temperature (T i ) profiles of 28% and 27%, respectively, relative to the experiment. TGLF overpredicts T e by 46%, with much larger variance than MMM, and underpredicts T i by 25%. As the ratio of kinetic to magnetic field pressure (β) is increased across NSTX discharges, TGLF predicts lower T e and significant flattening of the T i profile, conflicting with NSTX observations. When using an electrostatic version of TGLF, both T e and T i are substantially overpredicted, underscoring the importance of electromagnetic turbulence in the high β spherical tokamak regime. Additionally, calculations with neural net surrogate models for TGLF were performed outside of TRANSP with a time slice flux matching transport solver, finding better agreement with experiment than the TRANSP simulations, highlighting the impact of different transport solvers and simulation techniques. Furthermore, the reasonable agreement with experiment of temperature profiles predicted by MMM motivates a more detailed examination of the sensitivities of the TRANSP simulations with MMM to different NSTX plasma regimes in a companion paper (Lestz et al Plasma Phys. Control. Fusion 67 105030), in preparation for self-consistent, time-dependent predictive modeling of NSTX-U scenarios.

Lestz, J. B. [General Atomics, San Diego, CA (Unit↗
DOE OSTI2024

The role of RF-induced E×B flows in the mitigation of scrape-off-layer convective transport during ion cyclotron resonance heating

Abstract While multiple experiments have reported a decrease in intermittent fluctuations in the far Scrape-Off-Layer (SOL) during ion cyclotron resonance heating (Antar et al 2010 Phys. Rev. Lett. 105 165001, Li et al 2022 Nucl. Eng. Technol. 54 207–19, Antar et al 2012 Nucl. Fusion 52 103005), the physical mechanism behind this observation has not been fully established yet. In this work, we demonstrate, for the first time, a direct correlation between the amplitude of RF-induced E×B flows and turbulence suppression in the far SOL. Using the Gas Puff Imaging (GPI) diagnostic on Alcator C-Mod, we show again that Ion Cyclotron Range of Frequencies can significantly alter the flow in the SOL and introduce a shear layer in regions magnetically connected close to the antenna (Cziegler et al 2012 Plasma Phys. Control. Fusion 54 105019). With the 4-strap field-aligned antenna operated in dipole phasing, the ratio of the power coupled by the central two straps to the power coupled by the outer two straps was varied. The resulting RF-induced radial electric field magnitude thus varied, and we show that the impact on the far SOL turbulence correlates with the modified E×B velocity. We then apply a newly-developed blob tracking algorithm (Han et al 2022 Sci. Rep. 12 18142) to higher-resolution GPI videos in order to directly observe the process of blob shearing by RF-induced E×B flows. We show that the radially sheared poloidal flows act as a transport barrier by stretching, stopping, and destroying filaments, which is consistent with the observed difference in turbulence statistics.

Diab, R. (ORCID:0000000342065868)↗
DOE OSTI2022

Beam Spray Thresholds in ICF-Relevant Plasmas

Beam spray measurements suggest thresholds that are a factor of ≈2 to 15× less than expected based on the filamentation figure of merit often quoted in the literature. In this moderate-intensity regime, the relevant mechanism is forward stimulated Brillouin scattering (FSBS). Both weak ion acoustic wave damping and thermal enhancement of ion acoustic waves contribute to the low thresholds. FSBS imparts a redshift to the transmitted beam. Here, regarding the specific possibility of beam spray occurring outside the laser entrance holes of an indirectly driven hohlraum, this shift may be the most concerning feature owing to the high sensitivity of crossed-beam energy transfer to the interacting beam wavelengths in the subsequent overlap region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

On numerical errors to the fields surrounding a relativistically moving particle in PIC codes

The particle-in-cell (PIC) method is widely used to model the self-consistent interaction between discrete particles and electromagnetic fields. It has been successfully applied to problems across plasma physics including plasma based acceleration, inertial confinement fusion, magnetically confined fusion, space physics, astrophysics, high energy density plasmas. In many cases the physics involves how relativistic particles (those with high relativistic γ factors) are generated and interact with plasmas. However, when relativistic particles stream across the grid, both in “vacuum” and in plasma, many numerical issues may arise which can lead to unphysical results. We present a detailed analysis of how discretized Maxwell solvers used in PIC codes can lead to numerical errors to the fields that surround particles that move at relativistic speeds across the grid. Expressions for the axial electric field as integrals in k space are presented that reveal two types of errors. The first arises from errors to the numerator of the integrand and leads to unphysical fields that are antisymmetric about the particle. Furthermore, the second arises from errors to the denominator of the integrand and lead to Cherenkov like radiation in “vacuum”. These fields are not anti-symmetric, extend behind the particle, and cause the particle to accelerate or decelerate depending on the solver and parameters. The unphysical fields are studied in detail for two representative solvers - the Yee solver and the FFT based solver. Although the Cherenkov fields are absent, the space charge fields are still present in the fundamental Brillouin zone for the FFT based solvers. In addition, the Cherenkov fields are present in higher order zones for the FFT based solvers. Comparison between the analytical solutions and PIC simulation results are presented. A solution for eliminating these unphysical fields by modifying the k operator in the axial direction is also presented. Using a customized finite difference solver, this solution was successfully implemented into OSIRIS [1]. Results from the customized solver are also presented. Additionally, this solution will be useful for a beam of particles that all move in one direction with a small angular divergence.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Thomson Scattering with Gain

Thomson-scattering signals can be significantly modified by convective gains associated with the stimulated Raman scattering and stimulated Brillouin scattering instabilities as the scattered light traverses the probe beam. Gain results in amplification and narrowing of the red-shifted (i.e., Stokes) scattered light, whereas blue-shifted (anti-Stokes) features are depleted and broadened, making the relative prominence of these features a key signature of gain. Not accounting for instability growth has been a ubiquitous source of error in previous Thomson-scattering measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Prediction of electric and magnetic fields from spectral data using machine learning algorithms for Doppler-free saturation spectroscopy diagnostics

The prediction of electric and magnetic field amplitudes from atomic spectral data is critical for plasma control in fusion devices such as tokamaks. Conventional approaches that rely on physics-based models are computationally expensive and unsuitable for real-time applications. In this work, we develop and benchmark three machine learning algorithms—simulation-based inference (SBI), fully connected neural networks (FCNN), and histogram-based gradient boosting regression (GBR-Hist)—to infer field intensities directly from Doppler-free saturation spectroscopy (DFSS) spectra. Synthetic datasets of spectra were generated using the EZSSS code and evaluated both with and without added Poisson noise to mimic experimental conditions. We find that SBI achieves the highest accuracy and robustness, FCNN provides a strong balance of accuracy and computational efficiency for real-time applications, and GBR-Hist offers the fastest inference but is more sensitive to noise. Furthermore, these results demonstrate the potential of machine learning to accelerate DFSS analysis and enhance its utility for plasma diagnostics and control.

Doppler-free saturation spectroscopy↗