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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Coupled 2D MHD and runaway electron fluid simulations of SPARC disruptions

Runaway electrons (REs) generated during disruption events in tokamaks can carry mega-Ampère level currents, potentially causing damage to plasma-facing components. Understanding RE evolution during disruption events is important for evaluating strategies to mitigate RE damage. Using two-dimensional toroidally symmetric magnetohydrodynamic (MHD) simulations in M3D-C1, which incorporates a fluid RE model evolved self-consistently with the bulk MHD fluid, we examine the seeding and avalanching of REs during disruptions in the SPARC tokamak – a compact, high-field, high-current device designed to achieve a fusion gain Q > 2 in deuterium–tritium plasmas. The M3D-C1 simulations of unmitigated disruptions demonstrate RE plateau formation and peaking of the final current density, which agree well with the results of lower-fidelity reduced RE fluid models. This work provides the first systematic comparison and benchmarking of different primary sources, including activated tritium beta decay and Compton scattering, in SPARC disruption simulations with self-consistent MHD and RE coupling.

Compton scattering↗
DOE OSTI2025

Investigation of boron powder flow rates on real-time wall

The limit of boron flow rates for real-time conditioning of the first walls has been systematically investigated in the Experimental Advanced Superconducting Tokamak (EAST) with a full metal wall. Initially, solid boron injection demonstrated effective control over carbon impurities and deuterium recycling on the basis of pre-discharge boronization. A minimum flow rate, identified between 1.0 mg/s and 2.0 mg/s, was necessary for actively improving wall conditions under specific plasma operating scenarios, with this effect progressively enhancing as boron flow rates increased. Additionally, a maximum flow rate, estimated between 3.5 mg/s and 8.0 mg/s, was identified for these plasma conditions. When boron flow rates exceeded this maximum, boron-induced fueling effects influenced the plasma line-averaged density, and at excessively high flow rates, plasma disruption was observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Extended-MHD simulations of disruption mitigation via massive gas injection in SPARC

Recent developments to the M3D-C1 code enable higher fidelity modeling of disruptions, and can be applied in the design verification of reactor-scale tokamaks. Among these new capabilities is a method to mesh conducting vessel structures such as coils and passive plates, packing of the toroidal mesh around gas injectors, as well as anisotropic resistivity inside the vessel structures. We present extended-magnetohydrodynamic (MHD) simulations of disruption mitigation via massive gas injection (MGI) in SPARC. The goal of this study is to inform the disruption mitigation layout of SPARC and aid in the design of an effective gas injector configuration. Fully three-dimensional simulations with M3D-C1 are carried out for various injector configurations with the primary goal of determining the effect of different MGI parameters on heat loads and vessel forces. The simulations include a model for impurity ionization, recombination, advection and radiation, as well as spatially resolved conducting structures around the plasma. A localized mixture of deuterium and neon with a small toroidal and poloidal width is injected in up to six locations. We demonstrate that M3D-C1 can model a rapid shutdown via MGI using narrow and more realistic gas plumes than in previous simulations. As a result of the q = 1 surface in the SPARC baseline case a sawtooth is observed early in the simulations. Despite the sawtooth and the onset of edge MHD instabilities, the impurity distribution remains localized around the injector locations, but enables a radiative shutdown of the plasma. We find that using the maximum of six gas injectors results in a lower peaking factor and leads to a more even distribution of radiation toroidally than using two injectors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Advanced Simulation of ITER Core X-ray Crystal Spectroscopy

X-Ray Simulation Analysis (XRSA) is an analytical ray-tracing mixed code developed specifically for the ITER Core X-Ray Crystal Spectroscopy (XRCS-Core) diagnostic, which employs a dual-reflection configuration incorporating multiple pre-reflectors made of Highly Oriented Pyrolytic Graphite (HOPG) and spherically curved analyzing crystals. The ITER XRCS-Core is designed for high spectral resolution measurement in specific wavelength ranges, including narrow bands around 1.354 Å for W 64+ , 2.19 Å for Xe 51+ , and 2.555 Å for Xe 44+ and Xe 47+ , enabling diagnostic capability across a broad electron temperature range in the ITER plasma. XRSA facilitates efficient simulation of the spectral performance of this complex X-ray spectroscopic system. Recent updates to the XRSA code have incorporated two critical effects: auto-focusing, which specifically applies to HOPG, and polarization. These two effects are particularly important in the dual-reflection configuration used in the ITER XRCS-Core system to provide more accurate modeling results. Here, simulations conducted with the updated code demonstrate that polarization has a substantial impact on the performance of the dual-reflection system. Additionally, the combined influence of polarization and system layout introduces performance variations across channels through the same crystal.

Computer graphics↗
DOE OSTI2026

Physics of Ultra-High Energy Density Relativistic Plasmas from Ordered Nanostructures

This award funded a four-year program at Colorado State University on the physics of ultra-high energy density plasmas generated by focusing femtosecond, >1021 Wcm-2 laser pulses onto ordered nanowire arrays. We used the ALEPH petawatt-class Ti:Sa laser with its high contrast 400 nm second harmonic beamline to irradiate arrays of Ni and TiO2 nanowires made by template assisted electrodeposition and atomic layer deposition over a wide range of average densities. We met two of the three original project goals: radial ion acceleration via doppler shifted K shell emission, and time and spectrally resolved K shell measurements. The main experimental advance was integrating the Lawrence Livermore National Laboratory (LLNL) TREX ultrafast x-ray streak camera with a suite of spherically bent quartz Bragg crystal spectrometers built using Princeton Plasma Physics Laboratory (PPPL) methodology. This is the first TREX-on-bent-quartz setup fielded at a femtosecond, >1021 Wcm-2 facility. It is now a permanent diagnostic at ALEPH and is offered to the LaserNetUS user community. The program produced three peer reviewed journal articles that acknowledge this award. The platform will continue at the upgraded CSU ATLAS facility, where the ALEPH laser is being upgraded to 2 PW.

Hollinger, Reed [Colorado State University]↗
DOE OSTI2026

Decomposition of methane diluted with inert gas in an RF discharge cell

Decomposition of methane using non-thermal plasmas is an attractive route for producing hydrogen-rich gases and valuable carbon nanomaterials. Understanding how plasma discharge modes influence methane decomposition in optimizing plasma-assisted chemical conversion remains unexplored. This study explores the coupling between the discharge structure and product selectivity in RF capacitively coupled discharges operating in methane/inert gas mixtures in the pressure range of 2–3 torr. Under our experimental conditions the discharge exhibits mode transitions from uniform to striated in Ar and Kr and from diffuse to contracted in Ar and Kr with <5% CH 4 . The discharges in He and Ne remained uniform under our operating conditions, and their mixtures with CH 4 remained diffuse. A 0-d model for Ar/CH 4 discharge established a threshold for contraction while also asserting the importance of Ar$^{*}_{\textrm{m}}$ in the dissociation and ionization processes. The highest degree of methane decomposition, >99.7% with the main products of acetylene and graphitized solid carbon was achieved in the contracted discharge mode for both Kr or Ar with ⩽5% CH 4 . We demonstrate that contraction can play a crucial role in the effective decomposition of methane with value-added products and that both the electronic and thermal properties of plasma gas are responsible for this effect.

RF plasma↗
DOE OSTI2025

Three-Dimensional Helical-Rotating Plasma Structures in Beam-Generated Partially Magnetized Plasmas

Azimuthal structures emerging in beam-generated partially magnetized plasmas are investigated using three-dimensional particle in cell-Monte Carlo collision simulations. Two distinct instability regimes are identified at low pressures. Furthermore, when the gas pressure is sufficiently high, quasineutrality is attained and 2D spiral-arm structures form as a result of the development of a lower-hybrid instability, resulting in enhanced cross-field transport. At lower pressures, quasineutrality is not achieved and a 3D helical-rotating plasma structure forms due to development of the diocotron instability. Analytical formulas are proposed for the critical threshold pressure between these regimes and for the rotation frequency of the helical structures. Preliminary experimental verification is provided.

Penning discharges↗
DOE OSTI2025

Bootstrap current modeling in M3D-C1

Bootstrap current plays a crucial role in the equilibrium of magnetically confined plasmas, particularly in quasi-symmetric stellarators and in tokamaks, where it can represent bulk of the electric current density. Accurate modeling of this current is essential for understanding the magnetohydrodynamic (MHD) equilibrium and stability of these configurations. This study expands the modeling capabilities of M3D-C1, an extended-MHD code, by implementing self-consistent physics models for bootstrap current. It employs two analytical frameworks: a generalized Sauter model (Sauter et al. 1999 Phys. Plasmas vol. 6, no. 7, pp. 2834–2839), and a revised Sauter-like model (Redl et al. 2021 Phys. Plasmas vol. 28, no. 2, pp. 022502). The isomorphism described by Landreman et al. (2022 Phys. Rev. Lett. vol. 128, pp. 035001) is employed to apply these models to quasi-symmetric stellarators. The implementation in M3D-C1 is benchmarked against neoclassical codes, including NEO, XGCa and SFINCS, showing excellent agreement. These improvements allow M3D-C1 to self-consistently calculate the neoclassical contributions to plasma current in axisymmetric and quasi-symmetric configurations, providing a more accurate representation of the plasma behavior in these configurations. A workflow for evaluating the neoclassical transport using SFINCS with arbitrary toroidal equilibria calculated using M3D-C1 is also presented. This workflow enables a quantitative evaluation of the error in the Sauter-like model in cases that deviate from axi- or quasi-symmetry (e.g. through the development of an MHD instability).

fusion plasma↗
NASA NTRSDate not supplied

Nuclear Fusion Space Propulsion Research, Experimentation, Theory Development, and Systems Analysis Efforts Led by the NASA Glenn Research Center (1994-2004)

This review paper summarizes work performed from 1994 to 2004 by a several interrelated government, academic, and industry teams led by the NASA Glenn Research Center. The nuclear fusion space propulsion system concept was predicated on a spherical torus reactor, which enabled manned missions to the outer planets in less than one year. Moderate thrust levels (1,000’s lbf)from direct nuclear fusion exhaust plasma via a magnetic nozzle enabled high thrust-to-weight. An entire vehicle conceptual design, including an artificial gravity crew habitat, was created by the NASA Glenn Research Center. The proof of concept experiment test article and facility upgrade was performed at the Ohio State University which also included staff from the Ohio Aerospace Institute and Science Applications International Corporation.The governing equations for the plasma physics theory of magnetic nozzle operation were derived by the Los Alamos National Laboratory. A preliminary investigation of a proof of concept test utilizing Coaxial Helicity Ejection as a means to supply plasma for propulsion at the National Spherical Torus Experiment reactor was outlined by the Princeton Plasma Physics Laboratory. An industry standard on nuclear fusion propulsion conceptual design was created by two AIAA teams. Despite extremely modest funding levels, significant progress was made advancing the state of the art.The result was a coordinated conceptual, theoretical, and experimental design effort to guide fusion space propulsion development.

Nuclear fusion↗
DOE OSTI2025

Investigation of resonant layer response in electron viscosity regime

We present a supplementary study of previous work in Waybright and Park [Phys. Plasmas 31, 022502 (2024)] which demonstrates a substantial effect of electron viscosity on the resonant layer response to non-axisymmetric magnetic perturbations. A main refinement is to include a curl element of electron viscosity in the generalized Ohm's law. The refinement reveals a resonant layer response in the electron viscosity regime corresponding to slowly rotating and highly viscous plasmas.

Asymptotic analysis↗
DOE OSTI2025

Core-edge integrated predictive studies of ST40 and NSTX plasmas with the scrape-off layer box model

The ability to model the interplay between the core and edge of tokamak plasmas is crucial to designing both the plasma operating scenario of a fusion pilot plant and the design of the tokamak itself. Scrape-off-layer (SOL) models that are tailored to integrated scenario modeling need to have fast turn-around time and minimal computational burden to enable wide parameter-space coverage for design scoping. The SOL 0-D Box model is a reduced SOL model based on global power and particle balance that captures the essential physics of SOL transport with little computational cost. The usage of the 0-D Box model in core-edge coupled simulations has been demonstrated in both interpretive and predictive modes on a variety of devices. This paper presents a sensitivity study of the 0-D Box model to the input SOL heat-flux width for an ST40 plasma. This study demonstrates that accurate prediction of this width is crucial to predicting global performance parameters of a plasma scenario, such as energy confinement time and flux consumption. We also present an extension of the Box model to 1-D to allow for parallel variation of plasma parameters along the magnetic field lines. The 1-D Box model is then compared with SOLPS-ITER simulations of an NSTX plasma. Advantages and limitations of the Box model are discussed, and future directions are outlined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Strong gradient neoclassical transport in the plateau regime

Strong gradient regions in tokamaks such as the pedestal or internal transport barriers are regions of reduced turbulence where neoclassical transport can play a dominant role. In pedestals, gradient lengths comparable to the ion poloidal gyroradius have been measured. Standard neoclassical theory can miss important strong gradient effects in these regions because it assumes that the gradient length scales of density, temperature and potential are larger than the ion poloidal gyroradius. We extend plateau regime neoclassical theory into regions of gradients of the order of the ion poloidal gyroradius to capture strong gradient effects on transport processes in the pedestal and internal transport barriers. The fundamental idea behind our new framework is to keep a scale separation between the orbit widths and the gradient length scales by performing a large aspect ratio expansion. In the plateau regime, strong gradients cause poloidal variation that is in–out as well as up–down asymmetric. We study two different test cases assuming either radial force balance or the absence of turbulence and show that strong gradient effects can enhance or reduce standard neoclassical theory predictions in the plateau regime in strong gradient regions.

fusion plasma↗
DOE OSTI2025

Properties of Lower Hybrid Drift Waves and Energy Transfer Near and Inside the Magnetic Reconnection Electron Diffusion Regions

Here, we investigate properties of lower hybrid drift waves (LHDWs) near and inside the electron diffusion regions in 17 magnetopause and 9 magnetotail reconnection events using Magnetospheric MultiScale (MMS) mission observations. Our analysis show that LHDW type depend on the electron beta, as electron beta increases LHDWs become more electromagnetic in nature. The energy transfer from electromagnetic fields to particles is higher in electrostatic LHDWs and it is largely in parallel direction with respect to the local magnetic field. Linear dispersion analysis shows that electrostatic LHDWs are perpendicular propagating while electromagnetic waves propagate in oblique directions and the normalized wavenumber of all LHDW types falls within 0.5–0.8 range. A simple estimate on the LHDW nonlinear saturation suggests a possibly important roles played by these waves in supporting the reconnection electric field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Structured light approaches in laser-based plasma diagnostics

There is a growing demand for plasma diagnostics suitable for industrial plasma reactors employed in semiconductor nanofabrication, especially relevant to microelectronics and quantum information systems. Such reactors typically have limited optical access and pose considerable diagnostic challenges, including intense background emission, significant thermal loads, and contamination of optical viewports. In this study, we outline research into structured light techniques (laser beams with tailored spatial, temporal, or phase characteristics) that effectively overcome these issues using laser-induced fluorescence (LIF) as an example. The focus of presented diagnostics is on ion kinetics analysis within an industrial plasma source, although this approach is broadly applicable to other plasma systems and diagnostic contexts. We present a confocal LIF implementation using an axicon-generated Bessel annular beam, achieving spatial resolutions of approximately 5 mm at a focal distance of 300 mm, with potential improvements to about 1 mm. This approach matches conventional orthogonal LIF performance but requires only one optical port. Wavelength-modulation LIF employs nonlinear laser wavelength tuning to measure spectral line derivatives, suppressing background emission and enhancing details of spectral line shape. Additionally, we present new results on applying vortex beams (laser beams carrying orbital angular momentum, OAM) for LIF measurements in an industrial plasma device. These measurements enable simultaneous axial and tangential velocity determination using a single laser beam and have been tested with xenon ion transition. Initial quantification of results was performed. Together, these structured-light approaches provide robust, background-resilient, multi-dimensional diagnostics for complex plasma environments.

Romadanov, Ivan [Princeton Plasma Physics Laborato↗
DOE OSTI2026

Avoidance of disruptions on KSTAR due to vertical displacement events via novel real-time stability assessment

Disruption avoidance via the DECAF approach has been achieved on KSTAR using a novel real-time vertical stability assessment and a multiactuator feedback control strategy. The development of disruption avoidance strategies with reactor-relevant reliability is an urgent activity, enabling future fusion power plants. The stability metric employed is based on a new formulation of a vertical force gradient balance metric evaluated across the poloidal cross section of the plasma, with parameters tuned using historical data. Evaluation of this metric on a validation set of 400 recent KSTAR shots indicates >82% of Vertical displacement events can be avoided via feedback control. Essential to its calculation is the two-dimensional toroidal current density distribution in the plasma. Measurement of this profile faster than fully-converged equilibrium reconstructions can deliver is found to improve forecaster performance and is achieved with a surrogate model that takes as input magnetic diagnostic measurements and outputs the current profile on a basis comprising the top principal components of historical current profiles (from past equilibrium reconstructions). This method solves the non-uniqueness problem typically faced when reconstructing current profiles directly from diagnostics, while improving computational time and accuracy. On average, profiles produced by this model reach coefficients of determination of >0.99 with respect to those from equilibrium reconstructions. The avoidance actuators employed include poloidal field coils and an electron cyclotron current drive system. The multiactuator approach, as shown in this first demonstration, allows disruption avoidance while minimizing impact to operational performance. This ability, along with its flexibility and speed, makes this new approach an attractive option for avoiding these types of disruptions in reactors.

Tobin, Matthew [Columbia Univ., New York, NY (Unit↗
DOE OSTI2025

The importance of accounting for the Tolman correction to surface tension for nucleation and growth modeling of Fe clusters

Gibbs free energies of clusters are required for predictive modeling of cluster growth during condensation of a cooling vapor. Here, we present a straightforward method of calculating free energies of cluster formation using the data from molecular dynamics (MD) simulations. We apply this method to iron clusters having from 2 to 100 atoms. The energies obtained are verified by comparing to an MD-simulated equilibrium cluster size distribution in a sub-saturated vapor. We show that these free energies differ significantly from those obtained with a commonly used spherical cluster approximation, which relies on a surface tension coefficient of a flat surface, as it is used in the classical nucleation theory (CNT). We show that the spherical cluster approximation in CNT can be improved by using a cluster-size-dependent Tolman correction for the surface tension. The Tolman length and effective surface tension values were derived for iron clusters, and they significantly differ from the commonly used experimentally measured values. This improved approximation does not account for geometric magic number effects responsible for spikes and troughs in densities of neighbor cluster sizes. Nonetheless, it allows to more accurately model cluster formation from a cooling vapor. It better reproduces the condensation timeline, overall shape of the cluster size distribution, average cluster size, and the distribution width. In contrast, using a constant surface tension coefficient (as done in CNT) resulted in incorrect condensation dynamics and cluster size distributions. The analytical expression for cluster nucleation rate from CNT was updated to account for the size-dependence of cluster surface tension.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Importance of $\delta B_{\|}$ on ETG stability, turbulence, and transport in NSTX

This study employs electron-scale gyrokinetic simulations to investigate the electron temperature gradient (ETG) driven instabilities, turbulence, and transport in the pedestal region of the National Spherical Torus Experiment, comparing non-lithiated (narrow pedestal) and lithiated (wide pedestal) scenarios. Our findings reveal that, in the non-lithiated case, a branch of strongly unstable ETG modes exhibiting finite parallel magnetic field fluctuations ($\delta B_{\parallel} \neq 0$) emerges at the pedestal top and upper density pedestal region. This branch is uncovered only when $\delta B_{\parallel}$ is retained in the simulations and is associated with substantial electrostatic electron heat flux. This region of strong ETG transport corresponds to the only region in the plasma where the pressure gradient is far below the critical gradient for kinetic ballooning modes. We investigated the origin of this finite $\delta B_{\parallel}$ ETG branch by analyzing the gyrokinetic field equations. Nonlinear saturation is also analyzed and contrasted for simulations with and without $\delta B_{\parallel}$. In contrast with the nonlithiated case, ETG modes in the lithiated case produce substantial transport in the steep gradient region, but are negligible at the pedestal top.

ETG↗
DOE OSTI2026

Recent progress in atomic-scale controlled plasma processing

Atomic-scale control in plasma processing is becoming increasingly critical for fabricating of advanced semiconductor devices, particularly as the industry shifts toward three-dimensional (3D) architectures and high-aspect-ratio (HAR) structures. This review presents a comprehensive overview of recent developments in atomic-scale controlled plasma processes, organized along two key directions: the hierarchical structure of plasma–surface interactions and the generational evolution of atomic layer processing (ALP) technologies. We examined the gas phase, where molecular design enables selective generation of ions and radicals; the boundary layer, where transport phenomena govern species delivery into nanoscale features, and the surface, where temperature-dependent reactions and cyclic processing determine etching selectivity and precision. Building on this foundation, we outline five generations of ALP—from thermal atomic layer deposition to transport-aware, temporally and structurally decoupled processes—highlighting the increasing sophistication of process control. The review further explores the transition from empirical recipe development to science-based, data-driven methodologies. By integrating quantum-chemical modeling, advanced diagnostics, and machine learning, we demonstrated how predictive models can link plasma species composition to process outcomes, enabling autonomous and adaptive control strategies. Finally, this review discusses the broader societal implications of plasma process innovation through the E4 quartet: energy and resource efficiency, environmental sustainability, evolutionary advancement, and educational promotion. These principles guide the development of sustainable and intelligent atomic-scale manufacturing technologies that are not only technically advanced but also socially responsible.

Ishikawa, Kenji [Nagoya Univ. (Japan)] (ORCID:0000↗