Search NASA⌕ Search

SEARCH · Search NASA

Results for “Plasma waves”

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 OSTI2025

Direct evidence of the effect of a moderate external magnetic field on stimulated Raman scattering in the kinetic regime

We present results from an experiment carried out at the OMEGA-EP laser facility that investigated the effect of a perpendicular magnetic field on stimulated Raman scattering (SRS) and report the first direct measurement of magnetic mitigation of SRS. A 13-T magnetic field generated by pulsed-power coils was imposed on a gas jet plasma, and a novel three-picket interaction beam was used to explore SRS reflectivity for several plasma conditions within single shots. The time-resolved backscattered light shows that SRS was mitigated by the external 13-T magnetic field in the kinetic regime (kλ D ∼ 0.3, where k is the electron plasma wave's wavenumber and λ D is electron Debye length) and at a density of n e /n cr ∼ 0.10. On the other hand, we also measured an enhancement of SRS reflectivity at lower density (n e /n cr < 0.08). We discuss experimental results in the context of magnetohydrodynamic and particle-in-cell simulations that scan SRS dynamics for a variety of plasma conditions. While the experimental evidence of SRS mitigation validates prior work on the kinetic simulation of SRS in an external magnetic field, we also find that other mechanisms such as SRS rescatter can lead to an enhancement of measured SRS reflectivity in simulations of parameters relevant to our experiment.

Gas jet↗
DOE OSTI2026

Construction of generalized quasilinear diffusion coefficient using neural networks with physical restrictions

The quasilinear diffusion coefficient (D QL ) derived from our machine learning framework shows comparable trends with the ground truth D QL obtained from GENRAY-CQL3D simulations. Additionally, for the strong absorption cases, the radial current drive profiles generated using the D QL from our model exhibit consistent behavior with those obtained from the original simulation. These findings indicate the potential of our surrogate modeling approach with physical restrictions to replicate key wave–plasma interaction characteristics while reducing computational costs. Traditionally, calculating D QL for wave–particle interactions relies on computationally intensive wave simulations coupled with Fokker–Planck solvers. To address this challenge, we developed a machine learning-based surrogate model with physical restrictions derived from cold plasma theory and bounce-averaged damping effects. First, we establish the propagation domain of Lower Hybrid Waves in the (N∥, ρ) space by identifying the accessibility limit and determining the upper and lower bounds of N∥ using the Potential Power Deposition (PPD) method. Subsequently, leveraging a database constructed using Latin hypercube sampling alongside the underlying physical restrictions (e.g. PPD), machine learning methods including U-Net and Recurrent Neural Networks are employed to design a physics-restricted machine learning framework capable of reconstructing D QL .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Quasilinear theory: the lost ponderomotive effects and why they matter

Quasilinear theory (QLT) has been used for modeling wave–plasma interactions for decades but remains largely heuristic. Plasma inhomogeneity, ponderomotive effects, microscopic fluctuations, and collisions are not easily accommodated from first principles in QLT, and typically are ignored entirely, due to the limitations of the standard Fourier–Laplace global-mode approach. This results in inconsistencies, for example, violation of the action conservation for nonresonant waves. However, these issues can be avoided, and the theory can be substantially generalized and corrected, if QLT is formulated using more suitable analytical tools, particularly, the Weyl symbol calculus. Here, an attempt is made to deliver an accessible review of this modern formulation, provide intuitive calculations for special cases, and elaborate on the connection with the ‘oscillation-center QLT’ originally proposed by Dewar (Phys Fluids 16:1102, 1973). A Fokker–Planck equation for a ‘dressed’ distribution is derived from the Klimontovich equation and captures quasilinear diffusion, ponderomotive forces, and interactions with background fields for a generic Hamiltonian, so many known formulations of QLT for specific plasma models become corollaries of a single unifying theory. Also, waves are allowed to be off-shell (not constrained by a dispersion relation), which allows them to accommodate microscopic fluctuations. This leads to a collision integral of the Balescu–Lenard type that has all the usual properties but is not restricted to any specific plasma model. For on-shell waves, a generalized version of the classic oscillation-center QLT is obtained. Finally, combined with the wave-kinetic equation, this formulation not only conserves particles, momentum, and energy, like the classic QLT but also reinstates the action conservation for nonresonant waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Transversely pumped laser driven particle accelerator

We present a novel acceleration scheme capable of accelerating electrons and ions in an underdense plasma. Transversely Pumped Acceleration (TPA) uses multiple arrays of counter-propagating laser beamlets that focus onto a central acceleration axis. Tuning the injection timing and the spacing between the adjacent beamlets allows for precise control over the position and velocity of the intersection point of the counter-propagating beam arrays. This results in an accelerating structure that propagates orthogonal to the direction of laser propagation. We present the theory that sets the injection timing of the incoming pulses to accelerate electrons and ions with a tunable phase velocity plasma wave. Simulation results are also presented which demonstrate 1.12 GeV proton beams accelerated in 3.6 mm of plasma and electron acceleration gradients on the order of 1 TeV/m in a scheme that circumvents dephasing. This work has potential applications as a compact accelerator for medical physics and high energy physics colliders.

Laser-produced plasmas↗
DOE OSTI2025

Analysis of runaway electron driven whistler wave instability experiments

Data acquired on the DIII-D tokamak were analyzed. The data are from experiments that were conducted to study an instability that is driven unstable by intense populations of electrons with MeV energies that are known as runaway electrons. The instability is a type of plasma wave called a whistler wave that occurs at frequencies above the ion cyclotron frequency but well below the electron cyclotron frequency. The waves were measured by magnetic fluctuation coils that are embedded in the DIII-D vacuum vessel wall. After upgrades to this diagnostic were completed, new experiments were conducted on July 13, 2020 in order to measure the toroidal mode number of the whistler waves and to extend the frequency of the detected waves. Through the use of mixers, instability between 600-700 MHz was detected. (The initial experiments only measured up to 200 MHz.) Analysis of the data was led by Hari Choudhury, a PhD student at Columbia University. Mr. Choudhury has submitted two papers for publication that include contributions by UC Irvine (UCI) Professor Heidbrink and his graduated PhD student Genevieve DeGrandchamp: “Detailed Characterization of Runaway Electron Driven Whistler Waves in Low-Density DIII-D Discharges” and “First Demonstration of Resonant Pitch-Angle Scattering of Relativistic Electrons by Externally-Launched Helicon Waves.” The first paper, which has been submitted to Physics of Plasmas, has significant contributions to both the data and the interpretation by UCI scientists. In contrast, UCI contributions to the second paper, which has been submitted to Physical Review Letters, are relatively minor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Implementation of a Mesh refinement algorithm into the quasi-static PIC code QuickPIC

Plasma-based acceleration (PBA) has emerged as a promising candidate for the accelerator technology used to build a future linear collider and/or an advanced light source. In PBA, a trailing or witness particle beam is accelerated in the plasma wave wakefield (WF) created by a laser or particle beam driver. The WF is often nonlinear and involves the crossing of plasma particle trajectories in real space and thus particle-in-cell methods are used. The distance over which the drive beam evolves is several orders of magnitude larger than the wake wavelength. This large disparity in length scales is amenable to the quasi-static approach. Three-dimensional (3D), quasi-static (QS), particle-in-cell (PIC) codes, e.g., QuickPIC, have been shown to provide high fidelity simulation capability with 2-4 orders of magnitude speedup over 3D fully explicit PIC codes. In PBA, the witness beam needs to be matched to the focusing forces of the WF to reduce the emittance growth. In some linear collider designs, the matched spot size of the witness beam can be 2 to 3 orders of magnitude smaller than the spot size (and wavelength) of the wakefield. Such an additional disparity in length scales is ideal for mesh refinement where the WF within the witness beam is described on a finer mesh than the rest of the WF. A mesh refinement scheme is described that has been implemented into the 3D QS PIC code, QuickPIC. Very fine (high) resolution is used in a small spatial region that includes the witness beam and progressively coarser resolutions in the rest of the simulation domain. A fast multigrid Poisson solver has been implemented for the field solve on the refined meshes and a Fast Fourier Transform (FFT) based Poisson solver is used for the coarse mesh. The code has been parallelized with both MPI and OpenMP, and the parallel scalability has also been improved by using pipelining. A preliminary adaptive mesh refinement technique is described to optimize the computational time for simulations with an evolving witness beam size. Several test problems are used to verify that the mesh refinement algorithm provides accurate results. Additionally, the results are benchmarked against highly resolved simulations exhibiting near-azimuthal symmetry, performed using QPAD—a novel hybrid QS PIC code that uses a PIC description in the coordinates (r, ct – z) and a gridless description in the azimuthal angle, Φ.

Linear collider↗
DOE OSTI2026

Effects of density turbulence on helicon wave propagation in the core plasmas

Radio frequency wave propagation can be significantly affected by density irregularities, such as filaments in the scrape-off layer or instabilities in the core plasma. In this study, we examine the impact of edge turbulence on helicon wave propagation using the Petra-M simulation tool. To analyze the effect of edge turbulence, we utilize a realistic background plasma derived from XGC simulations, which includes spatial density fluctuations at the edge. This focus is particularly relevant for the DIII-D configuration, characterized by edge density turbulence of the core plasmas in a wide pedestal QH-mode. We focus on helicon wave propagation in the core plasma since the slow mode cannot propagate into the core plasma when the density is higher than the lower hybrid resonance, where a wide pedestal QH-mode occurs. The simulation results indicate that edge density fluctuations have a substantial impact on helicon wave coupling due to scattering in the core. We specifically demonstrate that the toroidal mode number and the level of density fluctuations are significant factors influencing scattering. While we have minimized the excitation of slow modes from the antenna, we still observe mode-converted slow modes resulting from the incoming helicon waves. The insights derived from these simulations will inform upcoming tokamak experiments regarding helicon antenna coupling in long pulse scenarios.

Kim, Eun-Hwa [Princeton Plasma Physics Laboratory ↗
DOE OSTI2026

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Overview of recent results from the ST40 compact high-field spherical tokamak

ST40 is a compact, high-field ($B_{\mathrm{T0}}\unicode{x2A7D} 2.1\,\,\,\textrm{T}$) spherical tokamak (ST) with a mission to expand the physics and technology basis for the ST route to commercial fusion. The ST40 research programme covers confinement and stability; solenoid-free start-up; high-performance operating scenarios; and plasma exhaust. In 2022, ST40 obtained central deuterium ion temperatures of $9.6 \pm 0.4\ \textrm{keV}$, demonstrating for the first time that pilot plant relevant ion temperatures can be reached in a compact, high-field ST. Analysis of these high-ion temperature plasmas is presented, including a summary of confinement, transport and microstability characteristics, and energetic particle instabilities. Recent scenario development activities have focused on establishing diverted H-mode plasmas across a range of toroidal fields and plasma currents, along with scenarios with high non-inductive current fractions. In future operations, beginning in 2025, a 1 MW dual frequency (104/137 GHz) electron cyclotron (EC) system will be installed to enable the study of EC and electron Bernstein wave plasma start-up and current drive. Predictive modelling of the potential performance of these systems is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Thermal Weibel instability induced magnetic fields co-exist with linear wakes in laser-ionized plasmas

When a moderately intense, few-picosecond-long laser pulse ionizes gas to produce an underdense plasma column, a linear relativistic plasma wave or wake can be excited by the self-modulation instability that may prove useful for multi-bunch acceleration of externally injected electrons or positrons to high energies in a short distance. At the same time, due to the anisotropic temperature distributions of the ionized plasma electrons, the Weibel instability can self-generate magnetic fields throughout such a plasma on a few picoseconds timescale that can persist even longer than the lifetime of the wake. In the present paper, we first show using simulations that both these effects do indeed co-exist in space and time in the plasma. Using our simulations, we make preliminary estimates of the contribution to the transverse emittance growth of an externally injected beam due to the Weibel magnetic fields in a few-millimeter-long plasma. We then present the results of an experiment that has allowed us to measure the spatiotemporal evolution of the magnetic fields using an ultrashort relativistic electron probe beam. Both the topology and the lifetime of the Weibel instability induced magnetic fields in the experiment are in reasonable agreement with the fields induced by the Weibel instability in the simulations.

43 PARTICLE ACCELERATORS↗
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

Photon acceleration of high-intensity vector vortex beams into the extreme ultraviolet

Extreme ultraviolet (XUV) light sources allow for the probing of bound electron dynamics on attosecond scales, interrogation of high-energy-density matter, and access to novel regimes of strong-field quantum electrodynamics. Despite the importance of these applications, coherent XUV sources remain relatively rare, and those that do exist are limited in their peak intensity and spatio-polarization structure. Here, we demonstrate that photon acceleration of an optical vector vortex pulse in the moving density gradient of an electron beam–driven plasma wave can produce a high-intensity, tunable-wavelength XUV pulse with the same vector vortex structure as the original pulse. Quasi-3D, boosted-frame particlein- cell simulations show the transition of optical vector vortex pulses with 800-nm wavelengths and intensities below 10 18 W/cm 2 to XUV vector vortex pulses with 36-nm wavelengths and intensities exceeding 10 20 W/cm 2 over a distance of 1.2 cm. The XUV pulses have sub-femtosecond durations and nearly flat phase fronts. The production of such high-quality, high-intensity XUV vector vortex pulses could expand the utility of XUV light as a diagnostic and driver of novel light–matter interactions.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Gate-Tunable Single Terahertz Meta-Atom Ultrastrong Light-Matter Coupling

We study the electrical tunability of ultrastrong light-matter interactions between a single terahertz circuit-based complementary split ring resonator (cSRR) and a two-dimensional electron gas. For this purpose, transmission spectroscopy measurements are performed under the influence of a strong magnetic field at different set points for the electric gate bias. The resulting Landau polariton dispersion depends on the applied electric bias, as the gating technique confines the electrons in-plane down to extremely subwavelength dimensions as small as d = 410 nm. This confinement allows for the excitation of standing plasma waves at zero magnetic field and an effective tunability of the electron number coupled to the THz resonator. This allows the normalized coupling strength to be tuned in situ from η = 0.46 down to η = 0.18. This is the first demonstration of terahertz far-field spectroscopy of an electrically tunable interaction between a single terahertz resonator and electrons in a GaAs quantum well heterostructure.

Landau levels↗
DOE OSTI2025

Laser wakefield acceleration driven by a discrete flying focus

Laser wakefield acceleration (LWFA) may enable the next generation of TeV-scale lepton colliders. Reaching such energies will likely require multiple LWFA stages to overcome limitations on the energy gain achievable in a single stage. The use of stages, however, introduces challenges such as alignment, adiabatic matching between stages, and a lower average accelerating gradient. Here, we propose a discrete flying focus that can deliver higher energy gain in a single stage, thereby reducing the number of stages required for a target energy. A sequence of laser pulses with staggered focal points and delays drives a plasma wave in which an electron beam experiences a near-constant accelerating gradient over distances beyond those attainable with a conventional pulse. Simulations demonstrate that a discrete flying focus with a total energy of 150 J can transfer 40 GeV per electron to a 50-pC beam in a single 30-cm stage, corresponding to 50 dephasing lengths.

laser wakefield acceleration↗
DOE OSTI2025

GloPW

This model predicts global plasma wave activity, with a focus on electromagnetic ion cyclotron (EMIC) waves in Earth’s magnetosphere.

Noh, Sung Jun↗
DOE OSTI2024

Effects of wave damping and finite perpendicular scale on three-dimensional Alfvén wave parametric decay in low-beta plasmas

Shear Alfvén wave parametric decay instability (PDI) provides a potential path toward significant wave dissipation and plasma heating. However, fundamental questions regarding how PDI is excited in a realistic three-dimensional (3D) open system and how the finite perpendicular wave scale—as found in both laboratory and space plasmas—affects the excitation remain poorly understood. Here, we present the first 3D, open-boundary, hybrid kinetic-fluid simulations of kinetic Alfvén wave PDI in low-beta plasmas. Key findings are that the PDI excitation is strongly limited by the wave damping present, including electron–ion collisional damping (represented by a constant resistivity) and geometrical attenuation associated with the finite-scale Alfvén wave, and ion Landau damping of the child acoustic wave. The perpendicular wave scale alone, however, plays no discernible role: waves of different perpendicular scales exhibit similar instability excitation as long as the magnitude of the parallel ponderomotive force remains unchanged. These findings are corroborated by theoretical analysis and estimates. This new understanding of 3D kinetic Alfvén wave PDI physics is essential for laboratory study of the basic plasma process and may also aid future evaluation of the relevance/role of PDI in low-beta space plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Modeling thermal radiation waves in silica plasmas for the Mooncat NIF experiment

The Mooncat experiment on the National Ignition Facility uses a laser-driven hohlraum to create a thermal radiation wave in a titanium-doped silica plasma. The titanium dopant enables absorption spectroscopy measurements to infer the temperature of the wave as it propagates. This measurement can be used to constrain multi-physics simulation codes to better understand when simulations do not match an experiment. In this paper, we present radiation-hydrodynamics simulations of the thermal radiation wave in the first full-platform shots of the Mooncat experiment. We examine the important parameters of the simulation, focusing on the radiation temperature source, the material model of the silica plasma as it pertains to radiation transport, and lateral leakage through a beryllium tube enclosing the silica. We compare different simulation modeling strategies to an analytic model of diffusive radiation transport and find that the simulation agrees with the analytic model when it is sufficiently simplified. These simulations show how radiation energy couples to matter to develop a shock wave in a radiative heat wave, an important topic in astrophysics and nuclear fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Shock wave formation in radiative plasmas

The temporal evolution of weak shocks in radiative media is theoretically investigated in this work. The structure of radiative shocks has traditionally been studied in a stationary framework. Their systematic classification is complex because layers of optically thick and thin regions alternate to form a radiatively-driven precursor and a temperature-relaxation layer, between which the hydrodynamic shock is embedded. In this work, we analyze the formation of weak shocks when two radiative plasmas with different pressures are put in contact. Applying a reductive perturbative method yields a Burgers-type equation that governs the temporal evolution of the perturbed variables including the radiation field. The conditions upon which optically thick and thin solutions exist have been derived and expressed as a function of the shock strength and Boltzmann number. Below a certain Boltzmann number threshold, weak shocks always become optically thick asymptotically in time, while thin solutions appear as transitory structures. The existence of an optically thin regime is related to the presence of an overdense layer in the compressed material. Scaling laws for the characteristic formation time and shock width are provided for each regime. The theoretical analysis is supported by FLASH simulations, and a comprehensive testcase has been designed to benchmark radiative hydrodynamic codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗