Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Plasma confinement”

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

Nonlinear control of the minimum safety factor in tokamaks by optimal allocation of spatially moving electron cyclotron current drive

The minimum value of the safety factor profile is related to the magnetohydrodynamic (MHD) stability of the plasma confined in a tokamak. Therefore, active control of the minimum safety factor may mitigate MHD instabilities that can degrade or even terminate plasma confinement. Typically, in most tokamak scenarios, the minimum safety factor evolves spatially with time, i.e., the location at which the safety factor achieves the minimum value changes with time. In addition to the inherent nonlinearities in the minimum safety factor evolution, its spatial variation makes the control design challenging. In particular, complexity in control design may arise from the need for time-dependent nonlinear models that account for spatial variation of the minimum safety factor. Furthermore, the minimum safety factor may drift to locations where the actuator authority is low. The problem of minimum safety factor control with target location tracking and moving electron cyclotron current drive (ECCD) is addressed in this work. A nonlinear time-dependent model that incorporates the spatial variation of the minimum safety factor is presented. A nonlinear controller based on optimal feedback linearization is developed to track a target minimum safety factor. The proposed controller treats the ECCD position as a controllable variable. In other words, the controller prescribes the ECCD position (in addition to the non-inductive powers) in real time based on an optimal criterion that is defined a priori. This work also presents the steps necessary to integrate the minimum safety factor controller with a total energy controller to achieve multiple control objectives simultaneously. In conclusion, the proposed integrated control algorithm is tested using nonlinear simulations in the Control Oriented Transport SIMulator (COTSIM) for a DIII-D tokamak scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

MPEX AI Digital Twins

All magnetically confined plasma fusion power plant concepts (Tokamak, Spherical Tokamak, Stellarator, Mirror, ...) must exhaust the heat and plasma from the core confinement region to the material walls. The primary channel for this exhaust is through a plasma divertor which directs plasma along open magnetic field lines to a material target. The Material Plasma Exposure eXperiment (MPEX) illustrated in Figure 1, is a high-power, steady-state linear plasma device designed to produce the plasma material interaction (PMI) conditions of the divertor of future magnetic confinement fusion power plants: energy flux 20MW/m 2 , ion fluence 1031/m 2 , pulse duration 106 sec. These goals of plasma exposure in MPEX are well beyond those achieved in magnetic fusion experimental devices. Successfully achieving these high power steady state conditions for long pulses requires operational control of the heating and particle sources and the plasma flux to the walls and target. The MPEX AI Hot Spot Controller, proposed in this project, will help achieve the operational milestones of MPEX. The MPEX device will begin commissioning at the end of FY26. A smaller proto-MPEX was operated for 14,666 plasma discharges and will resume operation in September of 2025 as proto-MPEX-lite, with reduced capability, to test a new window for the Helicon plasma source. The proto-MPEX data has undergone surrogate modeling with machine learning methods (R. Archibald, 2022 IEEE International Conference on Big Data). This proto-MPEX data will be used to begin development of the AI digital twins described in this white paper. The scientific mission of MPEX is to qualify materials of different composition for use in the high energy and plasma flux conditions of a fusion power plant. The materials exposed in MPEX will in some cases be exposed to high neutron fluxes at other ORNL facilities to measure the changes to their PMI properties. The targets exposed in MPEX will be transported under vacuum to a Surface Analysis Station (SAS). The SAS will be equipped with the following diagnostics: Focused Ion Beam (FIB) for trench milling, 100-400 angstrom resolution scanning electron microscope (SEM), surface mapping x-ray spectrometer, high resolution camera, and a future upgrade to a laser induced breakdown spectroscopy quadruple mass spectrometer (LIBS-QMS). The MPEX experiments will generate diverse pre- and post-exposure measurement data of detailed material properties down to the crystal grain level in 3D for post-exposure assessment of PMI damage (e.g. cracking, melting, erosion and redeposition of the material). Physics models for the PMI, and how the material composition and manufacturing impact its performance under high energy plasma exposure, need to be validated with MPEX data to guide the selection of new candidate materials. Our vision for the MPEX AI Digital Twins project is to supply experimental and physics model simulation data to train Artificial Intelligence (AI) models for data processing, analysis, operational control, PMI and materials simulation to maximize the scientific output of the MPEX device. Ultimately, an AI digital twin of MPEX material assessment metrics for tested and synthetic material types with simulated PMI will be trained by the AI Modeling Teams on the experimental and physics simulation data submitted to the American Science Cloud by this project. A purely empirical search for the best material is inefficient given the finite number of samples that can be tested on MPEX. In order to expand the material properties database for training the MPEX Material Assessment AI Digital Twin, and to gain physics understanding of the PMI processes, physics models of the material properties and PMI processes are required. The physics simulations provide detailed simulation data, like impact angles for plasma ions, sputtering yields, transport of the ionized sputtered target material in the plasma, and redeposition locations. This simulation data expands the measurement data for deeper physics understanding. The experimental data is essential to validate the PMI and material structure simulation models. The validated models can then be used to generate new simulation data of MPEX material assessments for synthetic material compositions that have not been exposed in MPEX. These predictive simulations, plus the whole experimental dataset, will be used to train the MPEX Material Assessment AI Digital Twin allowing a rapid generative AI search for new materials with reduced PMI damage by interpolating the domain of the training set. These new optimum materials can be simulated with the physics codes and/or tested in MPEX. The ability of AI neural networks to interpolate multi-dimensional parameter spaces and generate virtual data is exploited for a more efficient search for optimum materials. The advent of the Transformational AI Models Consortium (TAIMC) is an opportunity to engage with state of the art private and public AI developers to achieve the goals of the AI digital twins and AI accelerated physics models proposed in this project. Our partners at ORNL from the Advance Scientific Computing Research (ASCR) organization will collaborate in accelerating the integrated plasma material interaction simulation framework. This simulation framework will provide a platform for generating simulation data across a range of physical fidelities, including hybrid methods that produce multi-fidelity results. This data will be leveraged for AI model development, both for generation of surrogates and the automation of simulation campaigns. A part of the research below will include collaborative efforts with the TAIMC to (i) adapt data storage approaches to ensure AI-readiness, (ii) provide a protypical exemplar to inform and exercise constructed workflows, and (iii) generate and share data, using the TAIMC unified AI data standard, for foundational models that will be trained from multiple sources across the DOE complex. We will also collaborate with the TAIMC, as well as the planned AI modeling teams, to develop approaches for reducing the cost of data generation. These include tailored multi-fidelity approaches as well as fine-tuning strategies to augment general, large-scale foundational models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Analysis of the linear and nonlinear stability of Alfven eigenmodes and fish-bones in JET DT discharges: mode identification and shear flows generation

The plasma in future nuclear fusion reactors will be heated by neutral beam injectors (NBIs) and high frequency electromagnetic waves as well as fusion born alpha particles. Energetic particles (EPs), with energies up to two orders of magnitude larger than the thermal plasma, can trigger EP driven modes and induce harmful EP losses, reducing the plasma heating efficiency and the economical viability of the reactor. The present study is dedicated to analyze the Alfven Eigenmode (AE) activity in JET D–T discharges, the closest experiment to reactor-like operation performed until now. There, EP driven modes are induced by the combined effect of tangential NBIs and ion cyclotron resonance heating (ICRH) driven EP. Linear and nonlinear simulations are performed with the gyro-fluid FAR3d code to analyze the AE activity observed in the discharge 99896. The linear simulations reproduce the unstable n = 3 to 5 toroidal AEs (TAE) at the inner plasma region observed in the experiment, triggered by highly energetic passing deuterium populations injected by the tangential NBIs, further accelerated by the effect of the ICRH up to 1 MeV. In addition, fish-bones triggered by energetic trapped hydrogen induced by the ICRH are also reproduced. On the other hand, the alpha particles density is too small to destabilize AEs in the experiment. Nonetheless, increasing artificially the alpha density by one order of magnitude, an n = 1 beta induced AE can be destabilized in the inner plasma region. Nonlinear simulations indicate the generation of zonal structures during the AE/fish-bone saturation phase. TAE and fish-bones causes a rather weak increase of the passing D and trapped H EP (around 2%), respectively. Shear flows and zonal currents are generated during the saturation of TAE and fish-bones. Nonlinear simulations performed for D–T and pure deuterium thermal plasma indicate AE/fish-bone activity is weaker and shear flows are less intense in the pure deuterium case, trends consistent with the experimental observations that also indicates a deterioration of the thermal plasma confinement. Therefore, both numerical studies and experimental evidence indicate the generation of shear flows by AE/fish-bones could be connected with an improvement of the thermal plasma confinement.

AE↗
DOE OSTI2024

Verification of electromagnetic simulation capabilities in global gyrokinetic particle-in-cell code GTS

Recently, the numerical scheme presented by Mishchenko et al. enabled explicit gyrokinetic simulations of low-frequency electromagnetic instabilities in tokamaks at experimentally relevant values of plasma β⁠. This scheme resolved the long-standing cancellation problem that previously hindered gyrokinetic particle-in-cell code simulations of magnetohydrodynamic phenomena with inherently small parallel electric fields. Moreover, the scheme did not employ approximations that eliminate critical tearing-type instabilities. Here, we report on the implementation of this numerical scheme in the global gyrokinetic particle-in-cell code GTS. This implementation allows for a more complete and accurate picture of interaction between small scale turbulence and MHD modes in tokamaks. Additionally, we present a comprehensive set of verification simulations of numerous electromagnetic instabilities relevant to present-day tokamaks. These simulations encompass the kinetic ballooning mode, the internal kink mode, the tearing mode, the micro-tearing mode, and the toroidal Alfven eigenmode destabilized by energetic ions, which are all instrumental in understanding tokamak physics. We will also showcase the preliminary nonlinear simulations of kinetic ballooning instabilities and (2,1) island formation due to tearing mode instability. These simulations validate the accuracy of the scheme implementation and pave the way for studying how these instabilities affect plasma confinement and performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Multimachine analysis of intermittent fluctuations in the scrape-off layer of magnetically confined fusion plasmas

Far scrape-off layer fluctuations, also known as blobs/filaments, are compared for the first time across six tokamaks: Alcator C-Mod, DIII-D, TCV, KSTAR, MAST, and MAST-U, under the context of stochastic modeling in Ohmic and low confinement mode scenarios. In this model, which well describes these fluctuations, measurements from Langmuir probes are modeled as a superposition of uncorrelated, two-sided exponential temporal pulses. The key parameters of this model are the mean waiting times between consecutive pulses, the pulse duration times and the mean of the pulse amplitudes. The study examines the influence of line-averaged density, plasma current, magnetic field and aspect ratio. The mean amplitudes and pulse waiting times increase with the line-averaged density across devices, while the pulse duration times remain constant. Whereas the variation in waiting times is gradual across the density scan, the mean amplitudes increase abruptly as the density profile broadens and flattens for Alcator C-Mod. The blob waiting times decrease with increasing plasma current and the toroidal magnetic field, while they increase with aspect ratio for conventional tokamaks. The duration times follow the same trends as the waiting times in these parameters but vary more slowly, resulting in the trends in the waiting times dominating the intermittency level of the fluctuations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Mode transitions and spoke structures in E × B Penning discharge

Two-dimensional particle-in-cell simulations in the (radial-azimuthal) plane perpendicular to the axial direction of a cylindrical $E$ x $B$ Penning discharge are presented. The low-pressure discharge is self-consistently supported by plasma ionization from the electron beam injected axially, along the direction of the external magnetic field. It is shown that with the increasing strength of the external magnetic field, the discharge undergoes a sequence of transitions between several azimuthal modes. Azimuthal m > 1 spiral arm structures are excited at low magnetic field values as plasma confinement improves and the radial density profile becomes peaked. With a larger field, spiral arms with m > 1 are replaced by the m = 1 spoke mode, most clearly seen in plasma density. A transition from spiral arms to the spoke regime occurs when the plasma potential in the center changes from weakly positive (or zero) to negative. Further increase in the magnetic field results in a well-developed m = 1 spoke mode with additional small-scale higher-frequency m > 1 structures inside and around the spoke. It is shown that while ionization and collisions affect some characteristics of the observed fluctuations, the basic features of the spoke and m > 1 spiral structure remained similar without ionization. The role of energy conservation in small-scale high-frequency modes and spoke dynamics is discussed. It is demonstrated that in regimes with the m = 1 spoke mode, additional m = 4 harmonics of the ion and electron fluxes to the wall appear due to the square boundary. The frequency of the m = 1 mode is weakly affected by the geometry of the boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Final Report: Performance-structure characterization to improve REBCO Fusion conductor production at SuperPower

Several companies have used ReBCO tapes to make magnets operating at 20 K, at magnetic fields up to 20T suitable for plasma confining for controlled fusion applications. The transport properties of ReBCO conductors should be optimized for such conditions. NHMFL is one of the few laboratories where such tests can be accomplished. Recently, SuperPower Inc. came up with a new product, REBCO MOCVD HM tape, optimized for plasma confining applications. These HM tapes are designed to have increased critical currents at intermediate temperatures of 4.2K

36 MATERIALS SCIENCE↗
DOE OSTI2024

NSTX-U liquid metal core-edge facility (LMCE)

NSTX-U/LMCE will provide a unique and world-leading research facility to address the primary challenge to delivering economic and timely magnetic fusion energy, namely the need to develop a power and particle exhaust and first-wall system that can withstand very high edge heat fluxes, maximize energy confinement, and avoid the production of large masses of solid eroded first-wall material. The NSTX-U/LMCE facility will assess the ability of liquid metals (LMs) – especially liquid lithium – to provide a new boundary condition for magnetic fusion systems, to extend the lifetime of the plasma facing components (PFCs) and improve core plasma confinement. Such capability is needed to establish the basis for next-step fusion facilities including fusion pilot plants, and to maintain U.S. world leadership in core-edge integration research. NSTX-U/LMCE will leverage the ability to generate very high divertor perpendicular heat flux q⊥ ~ 100MW/m 2 , extensive diagnostics, and liquid-metal-applicable infrastructure of NSTX-U. NSTX-U/LMCE will provide access to a high-confinement plasma core with majority self-driven plasma current, the flexibility to test a range of liquid metal divertor concepts, access to a range of separatrix collisionalities (from high to very low), and the ability to controllably vary the first-wall temperature to vary the plasma- wall interaction physics on liquid lithium components. Further, NSTX-U/LMCE will utilize more reactor-relevant high-Z refractory-metal PFC substrates. With these capabilities the NSTX-U/LMCE facility will explore the full continuum of core-edge solutions ranging from high core radiated power, to conditions with radiative losses concentrated in the scrape-off layer (SOL), and ultimately low recycling conditions. The low collisionality SOL that may be accessible in the low recycling regime is relatively unexplored and will require a kinetic treatment of the edge, which can be addressed theoretically, and with experiments in LTX-β. Additional smaller-scale preparatory R&D facilities will be required to reduce the risk of premature technical/engineering failure of liquid metal systems implemented in NSTX-U. The NSTX-U/LMCE facility aligns very well with recommendations in the FESAC Long-Range Plan and NASEM Pilot Plant reports and the Bold Decadal Vision, will be unique in the world program throughout the next decade, and is garnering private company interest in utilizing NSTX-U/LMCE for development of LM PFCs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Global gyro-kinetic ion temperature gradient and trapped electron mode turbulence modelling in X -point geometry in negative and positive triangularity

Comparative modelling of ion temperature gradient (ITG)/trapped electron mode (TEM) turbulence in negative (NT) and positive (PT) triangularity plasma shapes was done using the nonlinear global full-f gyrokinetic particle code JOREK-GK in the realistic X-point tokamak geometry including the Scrape Off Layer (SOL) for TCV and DIII-D parameters. A comparison of JOREK-GK code with the gyrokinetic codes GS2 and GENE was done using NT/PT triangularity TCV L-modes parameters showing good agreement between codes in linear growth rates and clear beneficial effect of NT as compared to PT. Global non-linear modelling of the ITG/TEM saturated turbulence for realistic DIII-D NT pulses was done and compared with numerically constructed PT equilibrium with the same plasma profiles. Existence of longer correlation length of density fluctuations in PT compared to NT was demonstrated. Stronger and more sheared zonal flows are generated via Reynolds stress in NT compared to PT. These factors are stabilizing for TEM/ITG turbulence in NT and lead to smaller heat fluxes and heat conductivities in NT compared to PT configuration. Weak dependence of plasma confinement on collisionality and plasma rotation was found in modelling of DIII-D NT shots similar to the experiment. The confinement scaling with normalized ion gyro-radius ρ* was estimated both for NT and PT. Bohm-like scaling was obtained in both configurations, however with better confinement for NT compared to PT which could be favourable factor for reactor size machines.

ITG↗
DOE OSTI2025

Evolution of intermittent filaments in the scrape-off layer of NSTX

Filamentary structures naturally arise from background turbulence in the scrape-off layer (SOL) of plasmas, leading to significant particle and heat transport that can degrade overall plasma confinement. This enhanced transport can contribute to unacceptably high heat loads on plasma-facing components. As such, understanding the physics of SOL plasma filaments is critical for predicting and mitigating their effects in future fusion devices. On the National Spherical Torus Experiment (NSTX), plasma filaments—commonly referred to as “blobs”—were investigated using the gas-puff imaging (GPI) diagnostic in the edge and SOL regions. The analysis involved identifying, segmenting, and tracking the characteristic contours of the blobs in each frame of the GPI video sequences. Their evolution was characterized through shape descriptors, velocity, and angular velocity derived from their contour coordinates. The results indicate that as the blob area increases, their shapes become more concave and less circular, suggesting reduced structural stability in larger blobs. This result aligns with previous theoretical results where it was shown that larger blobs are more susceptible to instabilities [Krasheninnikov et al., J. Plasma Phys. 74, 679–717 (2008) and D'Ippolito et al., Phys. Plasmas 18, 060501 (2011)]. A positive correlation was observed between radial velocity and radial position, suggesting radially outward acceleration of the filaments, potentially driven by decreasing viscous drag toward the far SOL. Interestingly, blobs in background SOL turbulence exhibited minimal spinning in contrast to filaments originating from edge localized modes, which show substantial rotation during their paths [Lampert et al., Phys. Plasmas 29, 102502 (2022)]. Statistical analysis of the solidity and total curvature shape descriptors, along with their temporal evolution, revealed relatively broad, near-Gaussian distributions. This suggests that blob morphology is strongly influenced by stochastic turbulent processes in the surrounding plasma environment. Blob parameters were also compared with bulk plasma and radial profile measurements. Notable trends were found between blob rotation and poloidal velocity with collisionality and line-integrated density. These findings contribute to a deeper understanding of blob dynamics and provide valuable insights for refining SOL turbulence models.

Covariance and correlation↗
DOE OSTI2024

Validation of SOLPS-ITER and EDGE2D-EIRENE simulations for H, D, and T JET ITER-like wall low-confinement mode plasmas

Both experiments and simulations with SOLPS-ITER and EDGE2D-EIRENE show that the onset of detachment for the low-field side (LFS) divertor – defined here as the line-averaged upstream density (〈n e 〉 edge ) at which the plasma flux to the LFS target (I LFS-plate ) starts to decrease with increasing 〈n e 〉 edge – is independent of the isotope mass. However, there are three major simulation-experiment discrepancies: (i) the absolute values of I LFS-plate and the electron density (n e ) in the LFS divertor at the onset of detachment are significantly lower in simulations, i.e., approximately a factor of 2 for I LFS-plate and a factor of 3-4 for n e ; (ii) the degree of detachment – defined here as the difference between I LFS-plate at the onset of detachment and at an 〈n e 〉 edge value close to the density limit – is smaller in simulations compared to experiments; and (iii) the experimentally observed larger degree of detachment for D and T plasmas compared to H plasmas cannot be clearly distinguished from the simulation results. There are strong indications that discrepancy (i) is to a large extent caused by neglecting Lyman-opacity effects in our simulations. The simulations predict a similar net volumetric recombination source for all isotopes due to the fact that molecule-activated recombination (MAR) compensates for the reduced electron–ion recombination (EIR) for H, whereas MAR is negligible for D and T. This similar net volumetric recombination source for all isotopes leads to an isotope-independent degree of detachment in simulations. An analysis of the Balmer-α and Balmer-$γ$ emission confirms the underestimate of MAR in simulations (especially for D and T) for the JET metallic wall, which was previously observed for devices with a carbon wall. The underestimate of MAR is an important cause for discrepancy (ii) and the fact that there is a stronger underestimate of MAR for D and T than for H explains discrepancy (iii). Extending the plasma grid to the vessel wall increases I LFS-plate and n e at the onset of detachment by 25%, and the EIR source increases by 80% in detached conditions. Hence, while the extended grid results are closer to the experimental observations, the previously described qualitative discrepancies still persist.

Balmer emission↗
DOE OSTI2025

Gyrokinetic PIC Study on RMP Affected Neoclassical Transport in Toroidal Plasmas

In magnetically confined fusion plasmas, the breaking of ‘magnetic flux-surfaces’ due to resonant magnetic perturbations (RMPs) can generate magnetic islands and alter field topology to significantly impact plasma confinement and transport. Here, this work investigates the effect of magnetic islands on neoclassical radial energy transport within the core plasma of an analytic circular tokamak using the XGC-S global gyrokinetic particle-in-cell code. Findings from our simulations revealed substantial enhancements in electron neoclassical radial energy diffusivity in and around the islands, in addition to a newly observed two-peak structure at the O/X-points and outer island boundary in the electron diffusivity profile.

RMP↗
DOE OSTI2024

Nonlinear gyrokinetic modelling of high confinement negative triangularity plasmas

Abstract Nonlinear gyrokinetic simulations correctly predict particle as well as ion and electron energy fluxes of high confinement plasmas with a negative triangularity cross sectional shape, showing that core transport in these plasmas is well described by standard gyrokinetic models. Experimentally inferred power balance fluxes are mostly reproduced within one standard deviation across a wide portion of the minor radius. Experimental conditions are reproduced by ion scale simulations, without the need to include density and temperature profile curvature effects. The experimental case is used as baseline to predict that the non-dimensional confinement scaling in negative triangularity plasmas increases strongly with plasma current while slightly degrading at increasing normalized pressure and decreasing collisionality. Recent experiments showed that low toroidal rotation negatively impacts confinement; consistent with the experiment, simulations predict that low rotational shear significantly affects confinement unless the plasma effective charge is maintained above a minimum level. Core confinement is predicted to significantly degrade in low aspect ratio devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Efficient analysis of magnetic field line behavior in toroidal plasmas

The confinement of plasmas in tokamaks and stellarators depends on magnetic field lines lying in nested toroidal surfaces. The transition near the plasma edge away from the lines lying in magnetic surfaces defines properties of divertors. The transition in time defines properties of disruptions. Divertor design and disruption analyses require a detailed understanding of these transitions. The use of a Fourier transform coupled with a Gaussian window function allows far more information to be extracted about these transitions using far shorter field line integrations than can be obtained using traditional methods based on Poincaré plots. The physics of divertors and disruptions is reviewed to clarify why the type of information that can be gained from more efficient methods of analysis is of central importance to the fusion program based on magnetic confinement.

Boozer, Allen H. [Columbia Univ., New York, NY (Un↗
DOE OSTI2024

Localized phase contrast imaging at the Wendelstein 7-X stellarator

In its basic form, phase contrast imaging (PCI) provides line-integrated measurements of electron density fluctuations in plasmas. As turbulent fluctuations in magnetically confined plasmas have wave vectors almost perpendicular to the background magnetic field, the signals scattered by fluctuations from different parts of the PCI line-of-sight (LoS) are spatially separated in focal planes of the plasma. This allows localized PCI measurements by placing a mask in such a plane, to only permit signals from specific parts of the LoS to reach the PCI detectors. The present paper describes modeling and design of localization masks for the PCI system at the Wendelstein 7-X (W7-X) stellarator as well as the first results obtained using the masks in the recent long-pulse W7-X experimental campaign. During this project, we have extended the theory describing the mask response within the Fraunhofer diffraction model. As a novel development, we show from first principles that the mask response is determined by the fraction of power of the scattered beam spots that passes the mask. These insights have been used to select the W7-X mask design, consisting of a circular cutout, allowing the unscattered beam spot to pass the mask, with wedges covering a fixed angular range outside the central cutout. In the recent W7-X experimental campaign, the masks have verified the location of the main turbulence features observed by the PCI system and provided new information about the location of short-wavelength magnetohydrodynamic modes.

Physics↗
DOE OSTI2025

Enhanced Collisional Losses from a Magnetic Mirror Using the Lenard-Bernstein Collision Operator

Collisions are crucial in governing particle and energy transport in plasmas confined in a magnetic mirror trap. Modern gyrokinetic codes model transport in magnetic mirrors, but some utilize approximate model collision operators. This study focuses on a Pastukhov-style method of images calculation of particle and energy confinement times using a Lenard-Bernstein model collision operator. Prior work on parallel particle and energy balances used a different Fokker-Planck plasma collision operator. The method must be extended in non-trivial ways to study the Lenard-Bernstein operator. To assess the effectiveness of our approach, we compare our results with a modern finite element solver. Our findings reveal that the particle confinement time scales like a exp( a 2 ) using the Lenard-Bernstein operator, in contrast to the more accurate scaling that the Coulomb collision operator would yield a 2 exp( a 2 ), where a 2 is approximately proportional to the ambipolar potential. We propose that codes solving for collisional losses in magnetic mirrors utilizing the Lenard-Bernstein or Dougherty collision operator scale their collision frequency of any electrostatically confined species. This study illuminates the collision operator’s intricate role in the Pastukhov-style method of images calculation of collisional confinement.

fusion plasma↗
DOE OSTI2025

ThinCurr: An open-source 3D thin-wall eddy current modeling code for the analysis of large-scale systems of conducting structures

In this paper we present a new thin-wall eddy current modeling code, ThinCurr, for studying inductively-coupled currents in 3D conducting structures -- with primary application focused on the interaction between currents flowing in coils, plasma, and conducting structures of magnetically-confined plasma devices. The code utilizes a boundary finite element method on an unstructured, triangular grid to accurately capture device structures. The new code, part of the broader Open FUSION Toolkit, is open-source and designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. Scalability to large models is enabled through use of hierarchical off-diagonal low-rank compression of the inductance matrix, which is otherwise dense. Ease of handling large models of complicated geometry is further supported by automatic determination of supplemental elements through a greedy homology approach. Here, a detailed description of the numerical methods of the code and verification of the implementation of those methods using cross-code comparisons against the VALEN code and Ansys commercial analysis software is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A scaling law of the neutral penetration length and Balmer- α wing shape in high-temperature plasmas

Hydrogen atoms penetrating deep inside high-temperature magnetically confined plasmas by repetitive charge-exchange collisions result in a particle source of the plasma, which affects the plasma transport significantly. In this paper, we present an approximate solution of the fluid equations for neutral transport and an analytical representation of the neutral penetration length, in a simplified plasma geometry. This analysis predicts a power-law decay in the Balmer-α line wings which reflects the velocity distribution of the neutral atoms, with the power-law exponent analytically represented as well. These scaling laws are compared with a simple Monte–Carlo simulation and spectroscopic observations of Large Helical Device plasmas. Since the Balmer-α line wings are experimentally accessible, our formulation opens the possibility to quickly estimate the neutral penetration length from spectroscopic observations.

neutral opacity↗