Search NASA⌕ Search

SEARCH · Search NASA

Results for “PRINCETON PLASMA PHYSICS LABORATORY”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Assessing plasma face component thermal response to rotating 3D magnetic fields for SPARC tokamak

Thermal response simulations of plasma-facing components (PFCs) in the SPARC tokamak, performed with the HEAT code, show that three-dimensional (3D) heat loads resulting from stationary n=1 perturbations require highly radiative scenarios, with up to 95% of the power crossing the separatrix (P SOL ) being radiated, to maintain PFC temperatures within acceptable operational limits, whereas the application of slowly rotating 3D fields substantially reduces the thermal loads. The HEAT module, developed to predict heat loads from non-axisymmetric plasmas, is extended to model time-dependent heat flux patterns generated by rotating 3D fields, and a comprehensive thermal analysis is performed on PFCs subjected to both the maximum and minimum power loads, as well as to rotating heat flux distributions, to evaluate the temperature evolution for varying perturbation amplitudes and rotation frequencies. The extension of this analysis to 3D fields with toroidal mode number n=2 shows that this configuration leads to weaker localized heat flux peaks relative to the n=1 case, enabling safe operation with less than 80% of the power radiated when static 3D fields of low amplitude are applied, while using slowly rotating fields at higher amplitudes. These results indicate that n=2 perturbations are generally less detrimental to divertor power exhaust, emphasizing the strong dependence of divertor power exhaust on the characteristics of the applied 3D fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Visualization techniques for the gyrokinetic tokamak simulation code

Gyrokinetic simulations of plasma microturbulence in tokamaks are challenging to visualize because the compute grid follows the magnetic field lines that spiral around the torus. We have overcome this challenge by developing three new approaches that improve visualization of gyrokinetics. Our techniques work directly with the topology of magnetic flux surfaces where the simulation stores variables in concentric rings on poloidal planes (vertical cross sections of the torus). Our visualization preview step triangulates each consecutive pair of rings to display the data on a poloidal plane. The second visualization technique follows spiral field lines around the torus and constructs polygons to visualize a flux surface. Third, the poloidal triangles are connected between planes to form prisms that compose a 3-D model of the entire torus. The visualization workflow produces detailed geometry that matches the high resolution, irregular compute grid for every time step. The surface and solid models are displayed in scientific visualization programs to effectively explore and communicate the results, including fluctuation of electron density, ion temperature, and electrostatic potential. Highly detailed renderings verify plasma behavior along magnetic field lines over time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Enabling integrated AI control on DIII-D: a control system design with state-of-the-art experiments

We present the design and application of a general algorithm for Prediction And Control using MAchiNe learning (PACMAN) in DIII-D. Machine learning (ML)-based predictors and controllers have shown great promise in achieving regimes in which traditional controllers fail, such as tearing mode (TM) free scenarios, ELM-free scenarios and stable advanced tokamak conditions. The architecture presented here was deployed on DIII-D to facilitate the end-to-end implementation of advanced control experiments, from diagnostic processing to final actuation commands. This paper describes the detailed design of the algorithm and explains the motivation behind each design point. We also describe several successful ML control experiments in DIII-D using this algorithm, including a reinforcement learning controller targeting advanced non-inductive plasmas, a wide-pedestal quiescent H-mode ELM predictor, an Alfvén Eigenmode controller, a Model Predictive Control plasma profile controller and a state-machine TM predictor-controller. There is also discussion on guiding principles for real-time ML controller design and implementation.

machine learning↗
DOE OSTI2025

Anomalous resistivity and electron heating by lower hybrid drift waves inside reconnecting current sheets

Inside an electron diffusion region of laboratory reconnection experiments, the quasi-electrostatic lower hybrid drift wave (ES-LHDW) is observed when a significant guide field component is present. Through direct measurement of the anomalous drag term and quasilinear analysis, it is shown that ES-LHDW can account for approximately 20% of the mean reconnection electric field in a case with moderate guide field. This value exceeds the contribution from classical resistivity, which is around 10%. The effects of the Lorentz force term, often neglected for electrostatic waves, are crucial for the observed correlation between electric field and density fluctuations. Anomalous electron heating by the perturbed current and resistivity (2.6 MW/m 3 ) also surpasses the classical Ohmic heating, which is about 2.0 MW/m 3 . For the case with a high guide field, significantly higher local electron temperatures were observed during periods of strong ES-LHDW activity. A statistical analysis further supports electron heating by LHDW, showing a larger increase in electron temperature with a high guide field. Finally, data from the Magnetospheric Multiscale mission provide evidence of Landau damping of ES-LHDW, suggesting that ES-LHDW may contribute to the generation of nonthermal electrons along the direction parallel to the magnetic field.

Magnetic reconnection↗
DOE OSTI2026

Permanent magnets for ELM suppression in tokamaks: feasibility and operational compatibility

Permanent magnets can provide static magnetic fields without power supplies or feed lines, offering a simpler alternative to conventional electromagnets in tokamaks. This work examines permanent magnet arrays (PMAs) for edge localized mode (ELM) control in double-null (DN) configurations, where traditional RMP coils exhibit limited effectiveness. Linear plasma response calculations using IPEC assess high-field-side (HFS) permanent magnet placement, benchmarking performance against DIII-D low-field-side (LFS) internal coils (I-coils) while evaluating engineering constraints. Modeling results demonstrate that HFS permanent magnet configurations generate HFS plasma response amplitudes more than 5 times larger than conventional I-coil systems. While other resonant response metrics show comparable or reduced response relative to I-coils, the combination of conventional RMP systems and permanent magnets is expected to provide enhanced performance. Operational impact assessments of persistent magnetic fields, including sideband field effects and startup/ramp-up compatibility, reveal acceptable performance within established operational boundaries. This analysis establishes PMAs as a technically viable approach for DN ELM control with reductions in system complexity, motivating further experimental validation on existing tokamak facilities.

3D magnetic perturbation↗
DOE OSTI2026

Variations in lithium vapor cave performance predictions due to radial transport and recycling assumptions

The lithium vapor cave is a detached divertor design that uses a single private flux region baffle to contain a dense cloud of lithium vapor to dissipate heat flux. Plasma flows are created via fuel gas puffing in order to minimize lithium contamination of the main plasma. Significant modeling using the 2D edge code SOLPS-ITER has already been performed, predicting that sufficient target heat flux reductions (q$^{max}_{Target}$ MW m −2 ) with ( $n$ Li /$n$ e ) LCFS < 0.05 is possible in a case with 90 MW m −2 unmitigated heat flux in NSTX-U. Low heat flux and low upstream concentration was found with a variety of combinations of divertor geometries, target recycling coefficients, upstream plasma parameters, lithium evaporation locations and deuterium fueling locations, with variations in performance found for each design choice. However, the most universal uncertainty of SOLPS-ITER simulations has until now remained unaddressed systematically, namely the cross-field anomalous particle and heat diffusivities. This article aims to bound the uncertainty in lithium concentration prediction as a result of the assumed deuterium cross-field transport and recycling. For the simulations presented here, a factor of 2.1 increase in the upstream lithium density prediction across a factor of four decrease to the assumed deuterium radial particle diffusivity is found. This result is compared across different assumed deuterium recycling coefficients, known to reduce with lithium injection. Upstream lithium density is found to vary by a factor of 2.4 across a feasible range of recycling coefficients at similar $n$$^{OMP,sep}_{e}$.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Numerical modeling of impurity powder injection in W7-X

In this work, we present numerical simulation results of powder injection in W7-X using the EMC3-EIRENE and DIS codes. First, we model powder injection experiments performed in W7-X with the Probe Mounted Powder Injector. The simulation results qualitatively agree with visible imaging measurements. Secondly, we perform predictive simulations to guide the installation of an Impurity Powder Dropper in W7-X, allowing to choose in between several available non-vertical ports to maximize the amount of powder penetrating into the plasma, as well as the verticality of the port, to minimize sticking of the powders in the in-vessel stainless steel guiding tube. Port AEM41 is selected as the best candidate for IPD installation. The robustness of the simulation results has been verified for different plasma densities, powder materials and sizes, powder friction coefficient and changes in the plasma flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Directional finite difference method for directly solving 3D gyrokinetic field equations with enhanced accuracy

The gyrokinetic (GK) field equation is a three-dimensional (3D) elliptic equation, but it is often simplified to a set of two-dimensional (2D) equations by assuming that the field does not vary along a specific direction. However, this simplification can introduce inevitable 0th-order numerical errors, as nonlinear mode coupling in toroidal geometry can produce undesirable harmonic modes that violate the assumption. In this work, we propose a novel directional finite difference method (FDM) with a local coordinate transformation to better resolve the target field of interest. The directional FDM can accurately solve 3D GK field equations without simplifications, which can overcome the limitations of conventional methods. The accuracy and efficiency of different FDMs are analyzed in great detail for a variety of geometries, from simple 2D Cartesian coordinates to realistic 3D curvilinear coordinates. The 0th-order numerical errors of simplified 2D GK equations were found to be more problematic for low-harmonic modes and low aspect ratio geometries such as spherical tokamaks. On the other hand, the directional 3D FDM can accurately resolve a much wider range of harmonic modes aligned to the direction of interest, including the low-harmonic modes. In conclusion, we demonstrate that the directional 3D FDM is a highly effective algorithm for solving the 3D GK field equations, achieving accuracy improvements of 10 to 100 times or more, particularly for low-harmonic modes in spherical tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Commissioning and operation of a real-time multi-point Thomson scattering evaluation system at Wendelstein 7-X

A real-time multi-point Thomson scattering (rtMPTS) evaluation system has recently been installed and commissioned at the Wendelstein 7-X (W7-X) stellarator during the 2024 -2025 operational phases. This system has demonstrated successful evaluation of Thomson scattered spectra on millisecond timescales. We discuss initial results from this system, consisting of a comparison of different evaluation techniques, assessment of system evaluation times, and proof-of-principle control experiments carried out with the W7-X poloidal correlation reflectometer (PCR) system. This diagnostic represents the first step towards the development of a density profile controller at W7-X, which can be leveraged to support steady-state operation.

Mohammed, Ameer [Princeton Plasma Physics Laborato↗
DOE OSTI2025

The spherical tokamak advanced reactor (STAR) fusion power plant design

Scientific and technical advancements have been made that improve fusion’s prospects to provide a new energy source, showing enhanced plasma confinement conditions with plasma temperatures reaching or exceeding 100 million degrees. Overshadowing this progress is the challenge involved in developing an economically viable fusion power plant design. Many proposed next-step DEMO and pilot plant designs are extensions of existing physics-focused experimental devices defined to understand and control plasma operations to achieve and sustain a fusion reaction. Transitioning scientific and technical advancements into a functional power plant requires a dedicated focus on architectural designs that integrate diverse technologies, while optimizing physics conditions, with a focus on economic viability. This holistic approach is essential in turning the promise of fusion energy into a reality. The Spherical Tokamak Advanced Reactor (STAR) is a fusion power plant conceptual design with the architectural focus that strives to balance physics, engineering, and cost considerations. In conclusion, it has been set up to introduce relevant physics, engineering and concept features that an intermediate pilot plant might follow, with the goal of meeting system performances and economic requirements that lead to a commercially competitive fusion power plant.

Blanket segmentation↗
DOE OSTI2025

Active tungsten expulsion in ELM-absent H-mode plasmas via on-demand ELM triggering with lithium granule injection

Lithium granules gravitationally injected into the upper X-point region demonstrated on-demand edge-localized modes (ELM) triggering in otherwise ELM-suppressed H-mode plasmas on the Experimental Advanced Superconducting Tokamak. Sub-mm lithium granules dropped into enhanced D-alpha H-mode plasmas achieved a high triggering efficiency, while enabling ELM frequencies from several to hundreds of hertz. Core radiation from heavy impurities, dominated by W, was reduced by up to 60%, and the normalized energy confinement increased by up to 30%. At low injection frequencies, ELMs of substantially reduced size compared to spontaneous type-I ELMs were observed. At high injection frequencies, a transition to a mixed ELM phase occurred, characterized by intermittent larger ELMs and suppression of the quasi-coherent mode, achieving the most significant W reduction and energy confinement improvement. These results highlight a promising pathway for active W control via controlled, small ELMs in long-pulse, high-performance scenarios.

EDA H-mode↗
DOE OSTI2026

Effects of negative triangularity on microinstabilities in a low-recycling lithium-wall spherical tokamak

In this work, we present a linear gyrokinetic study of the impact of negative triangularity (NT) on microinstabilities in the Lithium Tokamak eXperiment-β (LTX-β), a low-recycling spherical tokamak with liquid lithium plasma-facing components that produce flat electron-temperature profiles [Elliott et al., IEEE Trans. Plasma Sci. 48, 1382 (2020)]. While NT is widely recognized as a stabilizing mechanism and often associated with improved confinement in conventional tokamaks, this study reveals that its effect is not universally stabilizing in the parameter regime of LTX-β and is shown to be highly sensitive to local equilibria. Using local linear simulations with the GS2 code [Kotschenreuther et al., Comput. Phys. Commun. 88, 128 (1995)] at ρ=0.3, 0.5, and 0.8 for two representative discharges (#103955 and #109355), and employing the Miller equilibrium model to isolate shaping effects, we find that NT can transition from stabilizing to destabilizing depending on magnetic shear, safety factor, and electron-temperature gradient. In shot #103955, NT reduces growth rates across radii, with strongest stabilization at the edge, whereas in shot #109355, it is stabilizing in the core but destabilizing at mid-radius and edge under experimental conditions. Parametric scans show that flattening the electron-temperature profile, increasing magnetic shear, and reducing the safety factor recover NT stabilization. These results demonstrate that NT stabilization is tunable rather than intrinsic and requires coordinated control of magnetic geometry and gradient drive.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Planar laser-induced fluorescence system for spatiotemporal ion velocity distribution function measurements

In this work, we present a planar laser-induced fluorescence (PLIF) system for measuring two-dimensional (2D), spatiotemporally resolved ion velocity distribution functions (IVDFs). A continuous-wave tunable diode laser produces a laser sheet that irradiates the plasma, and the resulting fluorescence is captured by an intensified CCD (ICCD) camera. Fluorescence images recorded at varying laser wavelengths are converted into 2D IVDFs using the Doppler shift principle. The developed diagnostic is implemented in an electron beam generated E × B plasma with a bulk plasma density of $\sim\!\!{10^{10}}{\text{c}}{{\text{m}}^{ - 3}}$ . The developed diagnostic is validated against a conventional single-point LIF method using photomultiplier tube-based detection, while significantly reducing the total measurement time by the number of spatial positions measured. The time-resolving capability of this diagnostic is tested by oscillating the plasma between two nominal operating modes with different density profiles and triggering the ICCD camera by the externally driven plasma oscillation. The measured 2D IVDF maps reveal several signatures of ion dynamics in this plasma source, including radially outflowing ions across the electric field and anomalous ion heating at the periphery, consistent with recent kinetic simulations and theoretical studies. A possible correlation between these ion kinetic features and rotating spoke structures is discussed.

E×B plasma↗
DOE OSTI2026

Synthetic modeling of soft x-ray emissivity for magnetic island analysis in LTX-β

We present a synthetic soft x-ray (SXR) forward-modeling framework that characterizes the emissivity structure of rotating magnetic islands in the Lithium Tokamak eXperiment-β across space and time. Magnetic islands associated with tearing modes produce modulations in line-integrated SXR brightness. Traditional tomographic methods struggle to resolve these structures in devices with limited sightlines, resulting in an under-determined inversion problem, or otherwise necessitate equilibrium reconstruction data, preventing their use in active control. Our approach removes this challenge by combining a fully three-dimensional ray-tracing model with a time-dependent emissivity prescription of a $m/n$ = 2/1 magnetic island. The model is validated through helical island geometry, rotation measured by magnetic diagnostics, and equilibrium constraints from PSI-Tri reconstructions. Synthetic brightness signals are generated for each photodiode sightline from a modeled emissivity profile and directly compared with experimental data from a tangential midplane SXR array. By fitting the synthetic diagnostic output to the observed brightness across time, we simultaneously infer key geometric parameters-such as island radial location, island width, and rotation frequency-without requiring full tomographic reconstruction, all with an average absolute deviation of less than 3%. This work demonstrates that forward modeling with a single tangential array can extract key island parameters in a spherical tokamak, providing a computationally efficient alternative to conventional SXR tomography and providing a pathway toward real-time magnetic-island characterization in future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Laboratory Observation of Transition from Collisional Slow to Collisionless Fast Reconnection

Temporal transition of an externally driven antiparallel asymmetric magnetic reconnection from collisional slow to collisionless fast regime is observed in a laboratory plasma for the first time. Here, this transition is initiated when the two-fluid Hall effect begins to dominate over collisional effects at the X point, characterized by the ratio of electron-ion collision mean free path to current sheet thickness exceeding unity. Prior to the transition, an enhanced reconnection electric field develops on the low-density side where electrons are heated both ohmically and by large-amplitude lower-hybrid drift waves (LHDW) before collisionality at the X-point drop significantly. These LHDWs generate anomalous resistivity accounting for 30% of the reconnection electric field. The observed time evolution is consistent with the hybrid collisionless-collisional reconnection scenario relevant to onset of asymmetric reconnection in natural plasmas.

Shi, Peiyun [Princeton Plasma Physics Laboratory (↗
DOE OSTI2025

Error field measurements and correction on MUSE permanent magnet stellarator

The magnetic field topology and error fields are measured in the MUSE permanent magnet stellarator with an electron beam and fluorescent rod. The effective magnetization of the magnets is determined by measuring rational flux surface locations. Electron magnetic drift effects must be considered self-consistently due to the small magnetic shear on MUSE. Error fields are characterized by flux surface shape deviation and magnetic island chain width and phase. Error fields are then corrected by adjusting the permanent magnet (PM) holders as rigid bodies. Resonant error fields are on the order of or smaller than 3 x 10 -6 of the toroidal field after correction and surface shape matches well with the design. Resonant components of the quasi-symmetry deviation are estimated, using a novel magnetic drift measurement, to be in the 10 −4 range. The PM approach is validated as a low-cost and precise way to build optimized stellarators.

electron beam mapping↗
DOE OSTI2025

3D modeling of n = 1 RMP driven heat fluxes on the SPARC tokamak PFCs using HEAT

3D heat flux calculations at the lower outer divertor plate of SPARC using the HEAT code show that 3D fields generated from error field correction coils can lead to enhanced peak heat fluxes up to 15 times larger compared to the axisymmetric case. Previously employed to simulate axisymmetric heat flux on 3D plasma facing components, the HEAT code can now predict 3D heat flux generated by non-axisymmetric plasmas. This is achieved via a new HEAT module which leverages the 3D field line tracing capabilities of MAFOT starting from an M3D-C1 (MHD resistive code) perturbed equilibria. The resulting heat flux is assigned using the magnetic footprint and the heat flux layer model, an extension of the 2D heat flux model also known as the Eich, to 3D non-axisymmetric plasmas. For SPARC, the new capabilities of HEAT are used to calculate the 3D heat loads resulting from n = 1 perturbation fields (with n indicating the toroidal periodicity) applied through a toroidal array of six picture frame coils with different amplitude. The comparison with the unperturbed case shows significant changes in shape and intensity of the heat flux profile. The results show that the application of n = 1 3D field leads to a localized enhancement of the heat flux peak, influenced by the wetted area impacted by the magnetic footprint, and the appearance of a secondary heat flux peak, whose intensity depends on amplitude of the applied 3D field and toroidal location.

3D heat flux↗