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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Accelerating kinetic plasma simulations with machine-learning-generated initial conditions

Computational models of plasma technologies often solve for the system operating conditions by time-stepping an initial value problem to a quasi-steady solution. However, the strongly nonlinear and multi-timescale nature of plasma dynamics often necessitate millions, or even hundreds of millions, of steps to reach convergence, reducing the effectiveness of these simulations for computer-aided engineering. We consider acceleration of kinetic plasma simulations via data-driven machine-learning-generated initial conditions, which initialize the simulations close to their final quasi-steady-state, thereby reducing the number of steps to reach convergence. Three machine-learning models are developed to predict the density and ion kinetic profiles of capacitively coupled plasma discharges relevant to the microelectronics industry. The models are trained on kinetic simulations over a range of device operating frequencies and pressures. Best performance was observed when simulations were initialized with ion kinetic profiles generated by a convolutional neural network, reducing the mean number of steps to reach convergence by 17.1× when compared to initialization with a zero-dimensional global model. We also outline a workflow for continuous data-driven model improvement and simulation speedup, with the aim of generating sufficient data for full device digital twins.

Artificial neural networks↗

DOE OSTI2026

Characterization of core neutrals using passive measurements of the D α spectrum near the X -point in the DIII-D tokamak

New spectroscopic measurements of deuterium Balmer-α emission are used to infer the spatial distribution of neutral particles near the X-point in a diverted high confinement mode plasma. The charge exchange neutral spectroscopy (CENS) diagnostic—recently installed on the DIII-D tokamak—uses 15 lines of sight extending from the edge of the confined region (ρ ≈ 0.7) to the X-point in lower single-null plasmas. Each CENS line of sight is spectrally resolved around the deuterium Balmer-α wavelength (6561 Å) to measure the Doppler shift, line broadening, and intensity of Dα emission. Thermal line broadening is used to identify emission from neutrals undergoing charge-exchange (CX) with high temperature ions in the confined plasma. This spectral information allows a more accurate determination of the neutral density deeper inside the confined plasma compared to traditional filter-based diagnostics. We present two methods of analyzing CENS measurements: (1) fitting the spectrum from each view-chord separately, and (2) a tomographic method for inverting the neutral density over a 2D region of space using all CENS views collectively. The neutral density profile is found to decay exponentially in the radial direction with two scale-lengths. In the pedestal region the neutral density decays at a rate approximately equal to the local mean-free-path for CX collisions, $L_{n_\mathrm{D0}}\approx \lambda_\text{CX}$. Further inside the plasma the neutral density decays at a rate equal to the mean-free-path for ionization, $L_{n_\mathrm{D0}}\approx \lambda_\text{inz.}$. The separatrix $n_\mathrm{D0}$ value is found to be approximately $2\times10^{15}$ m−3. Based on results from the 2D inversion, the density of neutrals is found to increase along the separatrix approaching the X-point by a factor of $\approx 5\times$ over the region covered by the CENS diagnostic.

X-point neutrals↗

DOE OSTI2025

Active Divertor Heat Flux Control using Impurity Powder Dropper

Divertor plasma-facing components (PFCs) in a tokamak are typically designed to withstand average steady-state heat loads of about 5–10 MW/m 2 , a limit that applies to both solid and liquid lithium (LL) PFCs. Exceeding these design values can result in surface damage to tungsten PFCs or excessive lithium (Li) evaporation in liquid lithium divertor (LLD) PFCs. Since exceeding the divertor heat load limits has serious consequences, it is therefore prudent to develop a tool to reduce the divertor heat load and bring the heat load to within the design limit without affecting the plasma performance. Active low Z impurity injection such as Li has been suggested as a potential solution to mitigate excess heat flux as suggested previously, given that non-coronal radiation can be quite large ~ 20–30 MJ per mole of injected Li. Li is considered desirable for reducing the edge neutral recycling helping to improve plasma energy confinement. In this paper, we model the Impurity Power Dropper (IPD) to investigate its potential of divertor heat flux control. The IPD is typically located at the top of the tokamak device and uses a vertical drift tube of a few meters. In the 2 m drift tube case, the IPD powder is accelerated to ~ 6 m/sec before reaching the plasma with the upper divertor configuration, matching the condition for the in-board side pellet injection case. By modeling the IPD geometry we determined the IPD powder deposition profile, and thus the non-coronal radiation and ionization profiles in time as well. From the enhanced radiation power loss, it is therefore possible to reduce the divertor heat load using the divertor simulation code. In conclusion, the IPD divertor heat flux control can be tested in the facilities with IPD including ST-40, DIII-D, EAST, WEST and NSTX-U.

Active lithium injection↗

DOE OSTI2025

Gyrokinetic simulations of the effects of magnetic islands on microturbulence in KSTAR

Gyrokinetic simulations are utilized to study effects of magnetic islands (MIs) on the ion temperature gradient (ITG) turbulence in the KSTAR tokamak with resonant magnetic perturbations. Simulations show that the transport is controlled by the nonlinear interactions between the ITG turbulence and self-generated vortex flows and zonal flows, leading to an anisotropic structure of fluctuation and transport on the poloidal plane and in the toroidal direction. MIs greatly enhance turbulent transport of both particle and heat. The turbulent transport exhibits variations in the toroidal direction, with transport through the resonant layer near the island X-point being enhanced when the X-point is located at the outer mid-plane. A quantitative agreement is shown between simulations and KSTAR experiments in terms of time frequency and perpendicular wavevector spectrum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Parallelized real-time physics codes for plasma control on DIII-D

A real-time safe multi-threading library was developed on the DIII-D plasma control system to optimize the real-time TORBEAM and real-time STRIDE physics codes. These physics codes are crucial for future fusion power plant operation as they provide information about electron cyclotron wave propagation and heating as well as inform about ideal plasma stability limits. The real-time TORBEAM code executed consistently in under 20 ms while the real-time STRIDE code computes in 100 ms. The multi-threading library developed in this work can be applied to other real-time physics-based codes that will be crucial for the next generation of fusion devices.

DIII-D↗

DOE OSTI2025

Machine-Learned Force Field for Molecular Dynamics Simulations of Nonequilibrium Ammonia Synthesis on Iron Catalysts

Ammonia (NH 3 ) is one of the most important industrial chemicals. The conventional NH 3 synthesis method-the Haber–Bosch process-converts atmospheric nitrogen (N 2 ) into NH 3 using H 2 with an iron (Fe) catalyst. However, this process requires high pressures (100–200 atm) and temperatures (700–800 K) near thermal equilibrium. Recently, Fe-based nanocatalysts have been reported to produce promising NH 3 yields under atmospheric pressures and temperature-modulated nonequilibrium conditions. Understanding the mechanism of nonequilibrium catalysis with programmed temperature variation could help to optimize this fully electrified and less energy-intensive process. Although reactive molecular dynamics (RMD) simulations can be a useful tool to model nonequilibrium catalytic processes, they require the development of accurate force fields (i.e., interatomic potentials). Here, we present a machine-learned (ML) force field within the Deep Potential MD (DPMD) framework, trained using periodic density functional theory (DFT) calculations, to model NH 3 synthesis on Fe catalysts with various surface adsorbates such as *N, *H, *N 2 , *H 2 , *NH, *NH 2 , and *NH 3 . Here, we generated the DFT data from static models of elementary reactions on the most stable (110) surface of body-centered cubic Fe, which then were augmented by data from constant number of particles–volume–temperature (NVT) DFT-MD trajectories at various temperatures. Finally, we utilized the fully optimized ML force field to investigate reaction dynamics at an Fe(110) surface at linearly increasing temperatures using NVT-DPMD simulations. Our simulations indicate that pulsed temperature ramping could prove favorable for NH3 synthesis. For example, we conducted ramping under multiple sets of conditions: (i) from 900 to 1200 K over periods of 0.1–0.3 ns for Fe surfaces precovered with N or NH along with H; and (ii) from 300 to 600 K over 0.1–0.3 ns for Fe surfaces precovered with NH 3 . While our simulations so far are limited to short time scales (very rapid heating), these observations shed light on the mechanism of the high NH 3 synthesis rate achieved in a novel temperature-modulated nonequilibrium catalytic reactor using pulsed heating and cooling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

From Weibel seeds to dynamo beyond pair-plasmas

Bridging the spatiotemporal scales of magnetic seed field generation and subsequent dynamo amplification in the weakly collisional intracluster medium presents an extreme numerical challenge. We perform collisionless turbulence simulations with initially unmagnetised electrons that capture both magnetic seed generation via the electron Weibel instability and the ensuing dynamo amplification. Going beyond existing pair-plasma studies, we use an ion-to-electron mass ratio of 100 for which we find electron and ion dynamics are sufficiently decoupled. These simulations are enabled by the 10-moment collisionless fluid solver of Gkeyll, which evolves the full pressure tensor for all species. The electron heat-flux closure regulates pressure isotropisation and effectively sets the magnetic Reynolds number. We investigate how the strength of the closure influences the transition between a regime reminiscent of previous kinetic pair-plasma simulations and a regime exhibiting dynamo behaviour qualitatively similar to magnetohydrodynamics.

astrophysical plasmas↗

DOE OSTI2025

Azimuthal and axial structures in 3D particle-in-cell simulation of Penning discharge

We report the results of a 3D particle-in-cell simulation of cylindrical Penning discharge with the axial magnetic field in the so-called reflex configuration, where the cathode and anti-cathode are biased to the same negative potential. The discharge is supported by thermal electron emission from the cathode. Electron and ion collisions, including ionization, are fully accounted for. The emphasis is on a specific regime in which the plasma potential at the center of the discharge is positive with respect to the chamber walls, serving as an anode. Spatial and temporal scales of the observed azimuthal and axial fluctuations and structures are characterized. It is suggested that azimuthal structures are caused by the dissipative gradient-drift instability. We find that the axial fluctuations related to the plasma-beam instabilities are weakly correlated with the azimuthal perturbations of the density, so that the azimuthal modes rotate as a whole and do not show any axial shear. The behavior of the electric potential is more involved, demonstrating intermittent standing wave and propagating structures in the axial direction that modulate the electron transport, producing the standing wave pattern (along the z-direction) in the radial electron flux.

3D PIC simulation of discharge↗

DOE OSTI2026

Full-wave modeling of RF waves in fusion plasmas with finite element method: Progress in past decades and its future role

This paper reviews the progress in computing radio frequency (RF) wave fields in fusion plasmas, specifically focusing on simulations utilizing the finite element method (FEM) over the past few decades. Computing RF wave fields in fusion plasmas presents unique challenges due to large simulation domains, complex antenna structures, non-local dielectric properties, and wide ranges of spatial scales. It highlights key developments and outlines future directions, primarily addressing waves in the ion cyclotron (IC) to lower hybrid (LH) frequency range. We begin with briefly revisiting earlier developments before the widespread availability of modern computer-aided engineering (CAE) software based on FEM. This historical perspective illuminates early progress and the physics difficulties that motivated ongoing work within the community. Modern wave simulations for RF antennas based on FEM are characterized by the use of detailed 3D antenna model geometry generated from engineering CAD software and localized wave dielectric model. Significant advancements have also been made in improving physics models to include phenomena such as RF sheath rectification and wave scattering. FEM-based RF simulations have also been applied to compute wave propagation in the core region, where the inclusion of non-local dielectric response is crucial. This is a challenging goal, and several promising approaches have been proposed in this area. Additionally, RF simulation development initiatives based on open-source libraries have gained popularity, demonstrating scalability and flexibility in extending physics models. This paper will discuss the advantages and disadvantages of using such a publicly available FEM library.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Plasma modification through boron particulate injection in the full tungsten environment of WEST

Recent experiments have confirmed the compatibility of extended boron particulate injections with high performance plasma discharges in the full tungsten (W) environment of WEST. Utilizing an impurity powder dropper (IPD) equipped with boron (B) powders a series of extended experimental programs providing controlled injections have quantified plasma response to varying levels of injection rate and total injection quantity. Calibration of injection quantities confirmed through post-situ testing of the IPD and cross-correlated with both high-speed camera illumination and spectroscopic measurement have allowed for the first time a fine scale determination of the effects of powder introduction on plasma performance. Plasma enhancement, consistent with turbulence reduction through profile modification, has been observed with sustained increases in the stored energy (WMHD), by 18%, electron temperature (T e ) by 35%, and neutron rate (N n ) by up to 200%, all of which scale positively with increasing powder injection rates. These injections have also resulted in both prompt and extended reductions in native impurity content, decreases in post injection radiated power, and strong decreases in divertor deuterium signatures signifying a reduction in recycling suggesting enhanced boron layer formation which provides a reduction of source terms and leads to enhanced gettering of main ion and impurity sources.

Lunsford, R. [Princeton Plasma Physics Laboratory ↗

DOE OSTI2026

Measurements and scaling of ion propulsion impulse during driven magnetic reconnection

Impulse scaling during magnetic reconnection, the magnetic energy conversion to kinetic energy, via direct Mach probe measurements in Magnetic Reconnection Experiment is examined. Ion exhaust velocity and impulse scalings with reconnecting magnetic field during the push phase of driven reconnection are presented. The outflows and impulse measurements are compared with global MHD simulations. Both measurements and simulations reveal a favorable scaling, greater than linear, of impulse with reconnecting field. These scaling results establish that magnetic reconnection could be utilized for plasma propulsion.

Ebrahimi, Fatima [Princeton Plasma Physics Laborat↗

DOE OSTI2025

Systematic analysis of magnetic equilibrium reconstruction with eddy currents on LTX- β

Magnetic equilibrium reconstruction in the Lithium Tokamak Experiment-Beta (LTX-β) is complicated by strong eddy currents and toroidal asymmetries arising from its segmented, close-fitting conducting shell. In the earlier experiment LTX, these three-dimensional (3D) conductor effects significantly distorted diagnostic signals and challenged conventional axisymmetric reconstruction methods. In this work, we demonstrate that accurate plasma equilibria can be recovered on LTX-β by incorporating a small number of dominant eddy current modes derived from realistic conductor models. Using the open-source $\tt{TokaMaker}$ Grad–Shafranov solver, we reconstruct equilibria across a systematically selected set of LTX-β discharges and validate them against the legacy $\tt{PSI-Tri}$ hybrid 2D-3D code. Our results show that fully 2D $\tt{TokaMaker}$ reconstructions achieve significantly improved agreement with flux loop and Mirnov probe measurements, reducing total chi-squared fitting errors, especially during startup. Finally, we develop new hybrid 2D-3D $\tt{TokaMaker}$ reconstructions by integrating the $\tt{ThinCurr}$ 3D eddy current model, which we optimize for significant further reduction of chi-squared in most scenarios. These findings underscore the importance of realistic wall-current modeling in short-pulse tokamaks, and establish a physics-based reconstruction framework that is extensible to devices with complex passive structures and 3D wall interactions.

Grad–Shafranov solver↗

DOE OSTI2025

A study of resistive peeling–ballooning modes across spherical tokamaks

We investigate how non-ideal-magnetohydrodynamics (MHD) effects, in particular plasma resistivity, impact the peeling–ballooning stability thresholds in spherical tokamaks. This analysis follows the discovery of resistive kink-peeling modes in ELMing National Spherical Torus Experiment (NSTX) discharges. In the present study we extend this modeling to ELMing pulses in the Mega Ampere Spherical Tokamak (MAST) and MAST—Upgrade (MAST-U), where we find a clear resistive scaling for peeling–ballooning modes. While in NSTX ideal-MHD predicts stability for ELMing discharges, in MAST-U we find that the plasma is slightly unstable to peeling–ballooning modes, but is fully stabilized once diamagnetic effects are considered in terms of a growth rate normalization. A resistive power law scaling is calculated for these modes on MAST-U, which lies in between that of tearing modes and resistive interchange modes. A comparison between M3D-C1 and NIMROD shows reasonable agreement for this scaling. Resistivity destabilizes the modes and the peeling–ballooning unstable domain is considerably expanded in both, MAST and MAST-U. In addition to the MAST/-U pulses we also analyze resistive PB stability in a NSTX-similarity discharge on DIII-D. While having a different aspect ratio from NSTX, this discharge uses NSTX-like shaping parameters, toroidal field and plasma current. By considering these discharges alongside NSTX cases, we identify conditions influencing the onset of resistive peeling–ballooning modes. Furthermore, our findings indicate that magnetic shear in the pedestal region is closely linked to the emergence of resistive edge modes.

MAST-U↗

DOE OSTI2025

Electromagnetic Analysis of ITER Electron Cyclotron Emission Components and Support Structures for the Preliminary Design Review

The ITER electron cyclotron emission (ECE) diagnostic system is located at Design Shield Module (DSM) 2, Equatorial Port (EP) 9 to measure electron temperature profile and electron temperature fluctuations and also to assess nonthermal electron distributions via the oblique view. Therefore, ECE has both radial and oblique views with two couples of mirrors for two different optical views. This ECE diagnostic system shall be exposed to significant power due to unabsorbed electron cyclotron heating power in the plasma. It shall also receive large electromagnetic (EM) loads up to 100 MN/m 3 force density due to the eddy currents generated by the short 16-ms transient plasma disruption. The global EM models of EP12 and EP11 with the worst case of plasma disruption Major Disruption Downward Exponential 16 ms Category III (MD_DW_EXP16MS_CATIII) have been provided by ITER. However, the support structure for these critical components, including four mirrors, is still under the design and development phase, especially during the period of preliminary design reviews (PDRs). For those small in-vessel components that are not modeled explicitly for the PDRs, we can extract the magnetic field (B) data and flux variations (dB/dt) at those locations from the global model and then use the B and dB/dt method to calculate the eddy currents, EM forces, and torques on those components in the local submodels. In conclusion, this paper will present the detailed PDR EM analysis results of the ECE components and support structures.

Electromagnetic↗

DOE OSTI2026

Time-dependent scenario modeling for the ST-E1 fusion power plant

ST-E1 is a low aspect ratio fusion power plant being designed by Tokamak Energy targeting 1.5 GW of fusion power. Characterization of the ST-E1 flat-top scenario is described elsewhere McNamara et al (2026 Nucl. Fusion 66 086008); here we focus on addressing the question of how to ramp-up the ST-E1 plasma from an initial state following breakdown and flux-surface formation to the target flat-top state. Being low-aspect ratio, the available solenoid flux of ST-E1 is limited. Therefore, particular consideration is placed on developing ramp-up scenarios that predominantly use inductive flux provided by external vertical field coils. Through time-dependent modeling with METIS, we show that this is possible when the ramp-up is performed at relatively high plasma density: although auxiliary current drive efficiency is reduced, this is significantly outweighed by (1) higher electron-ion collisional equilibration, (2) higher fusion power once ions become sufficiently hot, (3) higher poloidal beta for increased vertical-field flux, and (4) potentially favorable exhaust compatibilities. Ultimately, we show the target ST-E1 flat-top performance can be reached after a ramp-up period lasting 150 s using less than 40 Vs of solenoid flux (with vertical field providing ∼ 90 Vs of flux). The sensitivity to model assumptions are presented, with the general observation that deleterious effects can be mitigated through minor alterations of the auxiliary power temporal waveform and/or total auxiliary power level. The impact of a solenoid and the auxiliary power mix (electron cyclotron heating only versus electron and ion cyclotron heating) on the ST-E1 ramp-up success are also discussed in appendices. On this latter topic, we show that the effect of direct-ion heating during ramp-up is obscured by the uncertainty in the pedestal dynamics, identifying a clear line of future work required to make a definite decision on the ST-E1 auxiliary power mix.

current ramp-up↗

DOE OSTI2025

In situ studies of a molten metal anode ablation in a nearly atmospheric pressure DC arc

A DC arc with a meltable metal anode in a near-atmospheric pressure hydrocarbon gas is an emerging method for producing single-walled carbon nanotubes (SWCNTs). In these systems, evaporation of the molten metal anode determines the formation of catalyst seed particles needed for SWCNT growth, and therefore, should be monitored, controlled, and optimized. Evaluating the anode ablation rate by weighing the anode before and after a synthesis run is unfeasible due to anode carburization in the hydrocarbon atmosphere. To overcome this, we implemented a high-speed, 2D, 2-color pyrometry for reliable temperature measurements of the molten anode in a DC arc. The obtained temperature fields were used to calculate the anode ablation rates. Results showed the importance of resolving the arc and molten pool dynamics, as well as addressing the issue of reflections. Furthermore, significant changes in ablation rates were revealed upon addition of CH 4 , which must be considered when scaling up the production of SWCNTs.

Marangoni↗

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 OSTI2026

Difference in neoclassical edge flows between strongly negative and positive triangularities in the XGC gyrokinetic simulation

The neoclassical baseline study of a strongly negative triangularity (NT) plasma and the corresponding positive triangularity (PT) plasma is performed using the edge-specialized, total-f gyrokinetic code (X-point Gyrokinetic Code, XGC). A DIII-D-like plasma is used, based on the negative triangularity discharge of DIII-D #193 793. An artificial positive triangularity (PT) equilibrium has been constructed to compare the edge rotation physics at the same triangularity strength but with opposite sign, while keeping the same elongation and other geometric parameters. Carbon(6+) ions are added to the deuterium plasma at an experimentally relevant level. By using the experimental profile of carbon toroidal rotation profile as an input, XGC finds that the deuteron rotation is significantly different from the carbon rotation at the inboard and outboard midplanes, mostly caused by the difference in the Pfirsch–Schlüter rotation. More importantly, a significant difference in the X-point orbit loss physics, thus the rotation source, is found between the positive and negative triangularity equilibrium models. However, it is also found that the agreement between the present neoclassical simulation and the experimental NT data is validated only within the middle of the pedestal slope, indicating the importance of edge turbulence. Finally, this study could establish a baseline for the multiphysics, multiscale studies that include turbulence of negative triangularity plasmas.

Deuterium↗