Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA SIMULATION”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Computational study of nonlinear plasma waves. I - Simulation model and monochromatic wave propagation. II - Sideband instability and satellite growth

A hybrid plasma simulation model is described and applied to the study of electrostatic wave propagation in a one-dimensional Maxwellian plasma with periodic boundary conditions. The model employs a cloud-in-cell scheme which can drastically reduce the fluctuations in particle simulation models and greatly ease the computational difficulties of the Vlasov equation approach. A grid in velocity space is introduced and the particles are represented by points in the x-v phase space. The model is tested first in the absence of an applied signal and then in the presence of a small-amplitude perturbation. The method is also used to study propagation of an essentially monochromatic plane wave. Results on amplitude oscillation and nonlinear frequency shift are compared with available theories.

Matsuda, Y.↗
DOE OSTI2015

Computational hydrodynamics and optical performance of inductively-coupled plasma adaptive lenses

This study addresses the optical performance of a plasma adaptive lens for aero-optical applications by using both axisymmetric and three-dimensional numerical simulations. Plasma adaptive lenses are based on the effects of free electrons on the phase velocity of incident light, which, in theory, can be used as a phase-conjugation mechanism. A closed cylindrical chamber filled with Argon plasma is used as a model lens into which a beam of light is launched. The plasma is sustained by applying a radio-frequency electric current through a coil that envelops the chamber. Four different operating conditions, ranging from low to high powers and induction frequencies, are employed in the simulations. The numerical simulations reveal complex hydrodynamic phenomena related to buoyant and electromagnetic laminar transport, which generate, respectively, large recirculating cells and wall-normal compression stresses in the form of local stagnation-point flows. In the axisymmetric simulations, the plasma motion is coupled with near-wall axial striations in the electron-density field, some of which propagate in the form of low-frequency traveling disturbances adjacent to vortical quadrupoles that are reminiscent of Taylor-Görtler flow structures in centrifugally unstable flows. Although the refractive-index fields obtained from axisymmetric simulations lead to smooth beam wavefronts, they are found to be unstable to azimuthal disturbances in three of the four three-dimensional cases considered. The azimuthal striations are optically detrimental, since they produce high-order angular aberrations that account for most of the beam wavefront error. A fourth case is computed at high input power and high induction frequency, which displays the best optical properties among all the three-dimensional simulations considered. In particular, the increase in induction frequency prevents local thermalization and leads to an axisymmetric distribution of electrons even after introduction of spatial disturbances. The results highlight the importance of accounting for spatial effects in the numerical computations when optical analyses of plasma lenses are pursued in this range of operating conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Impact of spatiotemporal smoothing on the two-plasmon–decay instability

Higher levels of hot electrons from the two-plasmon–decay instability are observed when smoothing by spectral dispersion (SSD) is turned off in directly driven inertial confinement fusion experiments at the Omega Laser Facility. This finding is explained using a hot-spot model based on speckle statistics and simulation results from the laser–plasma simulation environment. The model accurately reproduces the relative increase in hot-electron activity at two different drive intensities, although it slightly overestimates the absolute number of hot electrons in all cases. Extrapolating from the current ≈ 360-GHz system while adhering to the logic of the hot-spot model suggests that larger SSD bandwidth should significantly mitigate hot-electron generation, and legacy 1-THz OMEGA experiments appear to support this conclusion. These results demonstrate that it is essential to account for laser speckles and spatiotemporal smoothing to obtain quantitative agreement with experiments. A compilation of hot-electron data from the past two decades reveals several other important points: (1) many prior experiments are more easily understood using recent results from multibeam absolute instability theory; and (2) experiments with ignition-scale conditions produce less hot electrons compared to OMEGA spherical experiments for a given vacuum overlapped intensity, which is a promising result for validating performance predictions based on hydrodynamic scaling relations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Ionization dynamics of intense laser-produced argon plasmas revealed by NLTE modeling

The ionization dynamics and transient behavior of under-dense plasma irradiated by an intense laser are investigated. We report two significant effects in the ionization behavior: (1) a surprisingly large delay in ionization response and (2) a stepwise ionization process which involves collisional and laser-driven photoionization (LDP) processes. Ionization induced by intense lasers can exhibit delayed responses due to rapid changes in conditions, particularly when atomic transition processes occur more slowly than the relevant time scales. Furthermore, modeling reveals that the two-step ionization process—collisional excitation followed by LDP—plays an important role in this ionization delay, with collisional excitation acting as the bottleneck. Even low-energy photons (∼3.5 eV) can predominantly ionize plasmas, challenging the conventional belief that such energies are insufficient to overcome the binding energy of bound electrons. These findings underscore the necessity of including such processes into plasma simulations for various laser-plasma experiments.

Collisional excitation↗
DOE OSTI2024

Simulating nonlinear optical processes on a superconducting quantum device

Simulating plasma physics on quantum computers is difficult because most problems of interest are nonlinear, but quantum computers are not naturally suitable for nonlinear operations. In weakly nonlinear regimes, plasma problems can be modelled as wave–wave interactions. In this paper, we develop a quantization approach to convert nonlinear wave–wave interaction problems to Hamiltonian simulation problems. We demonstrate our approach using two qubits on a superconducting device. Unlike a photonic device, a superconducting device does not naturally have the desired interactions in its native Hamiltonian. Nevertheless, Hamiltonian simulations can still be performed by decomposing required unitary operations into native gates. To improve experimental results, we employ a range of error-mitigation techniques. Apart from readout error mitigation, we use randomized compilation to transform undiagnosed coherent errors into well-behaved stochastic Pauli channels. Moreover, to compensate for stochastic noise, we rescale exponentially decaying probability amplitudes using rates measured from cycle benchmarking. We carefully consider how different choices of product-formula algorithms affect the overall error and show how a trade-off can be made to best utilize limited quantum resources. This study provides an example of how plasma problems may be solved on near-term quantum computing platforms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

KINETIC SIMULATIONS OF THE LOWEST-ORDER UNSTABLE MODE OF RELATIVISTIC MAGNETOSTATIC EQUILIBRIA

We present the results of particle-in-cell numerical pair plasma simulations of relativistic two-dimensional magnetostatic equilibria known as the “Arnold–Beltrami–Childress” fields. In particular, we focus on the lowest-order unstable configuration consisting of two minima and two maxima of the magnetic vector potential. Breaking of the initial symmetry leads to exponential growth of the electric energy and to the formation of two current layers, which is consistent with the picture of “X-point collapse” first described by Syrovatskii. Magnetic reconnection within the layers heats a fraction of particles to very high energies. After the saturation of the linear instability, the current layers are disrupted and the system evolves chaotically, diffusing the particle energies in a stochastic second-order Fermi process, leading to the formation of power-law energy distributions. The power-law slopes harden with the increasing mean magnetization, but they are significantly softer than those produced in simulations initiated from Harris-type layers. The maximum particle energy is proportional to the mean magnetization, which is attributed partly to the increase of the effective electric field and partly to the increase of the acceleration timescale. We describe in detail the evolving structure of the dynamical current layers and report on the conservation of magnetic helicity. These results can be applied to highly magnetized astrophysical environments, where ideal plasma instabilities trigger rapid magnetic dissipation with efficient particle acceleration and flares of high-energy radiation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Hermes-3

Hermes-3 is a plasma simulation model built on BOUT++, developed as part of research studies of the edge and divertor regions of magnetically confined fusion plasmas such as tokamaks. The main purpose of this model is multi-species simulation of plasmas containing a mixture of deuterium, tritium, helium and other plasma and neutral gas species.

Dudson, BenjaminD↗
DOE OSTI2022

hPIC2: A hardware-accelerated, hybrid particle-in-cell code for dynamic plasma-material interactions

The exascale era of high performance computing promises to bring the field of computational plasma physics ever closer to the goal of accurate multiscale modeling. Such computers will rely on hardware acceleration to offload work to dedicated components, notably general-purpose graphics processing units (GPUs). However, devices from different manufacturers require software to be written with different parallel programming models, greatly increasing the code maintenance burden of applications designed to perform on more than one such device. hPIC2 is a hybrid plasma simulation code developed with the Kokkos performance portability framework to target the architectures that will drive exascale computing for the foreseeable future. As a hybrid simulation code, hPIC2 investigates the simultaneous use of various plasma models on the same domain, at the same time. hPIC2 also optionally couples to RustBCA, which accurately models ion-material interactions using the binary collision approximation (BCA) method. In conclusion, hPIC2 therefore achieves scalable performance on a variety of computing architectures when simulating complex and diverse plasmas, particularly near plasma-material interfaces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

ZERNIPAX: A fast and accurate Zernike polynomial calculator in Python

Zernike polynomials serve as an orthogonal basis on the unit disc, and have proven to be effective in optics simulations, astrophysics, and more recently in plasma simulations. Unlike Bessel functions, Zernike polynomials are inherently finite and smooth at the disc center (r=0), ensuring continuous differentiability along the axis. This property makes them particularly suitable for simulations, requiring no additional handling at the origin. We developed ZERNIPAX, an open-source Python package capable of utilizing CPU/GPUs, leveraging Google's JAX package and available on GitHub as well as the Python software repository PyPI. Furthermore, our implementation of the recursion relation between Jacobi polynomials significantly improves computation time compared to alternative methods by use of parallel computing while still performing more accurately for high-mode numbers.

Astrophysics↗
DOE OSTI2017

Magnetohydrodynamic simulation study of plasma jets and plasma-surface contact in coaxial plasma accelerators

Recent experiments by Loebner et al. [IEEE Trans. Plasma Sci. 44, 1534 (2016)] studied the effect of a hypervelocity jet emanating from a coaxial plasma accelerator incident on target surfaces in an effort to mimic the transient loading created during edge localized mode disruption events in fusion plasmas. In this study, we present a magnetohydrodynamic (MHD) numerical model to simulate plasma jet formation and plasma-surface contact in this coaxial plasma accelerator experiment. The MHD system of equations is spatially discretized using a cell-centered finite volume formulation. The temporal discretization is performed using a fully implicit backward Euler scheme and the resultant stiff system of nonlinear equations is solved using the Newton method. The numerical model is employed to obtain some key insights into the physical processes responsible for the generation of extreme stagnation conditions on the target surfaces. Simulations of the plume (without the target plate) are performed to isolate and study phenomena such as the magnetic pinch effect that is responsible for launching pressure pulses into the jet free stream. The simulations also yield insights into the incipient conditions responsible for producing the pinch, such as the formation of conductive channels. The jet-target impact studies indicate the existence of two distinct stages involved in the plasma-surface interaction. A fast transient stage characterized by a thin normal shock transitions into a pseudo-steady stage that exhibits an extended oblique shock structure. A quadratic scaling of the pinch and stagnation conditions with the total current discharged between the electrodes is in qualitative agreement with the results obtained in the experiments. Finally, this also illustrates the dominant contribution of the magnetic pressure term in determining the magnitude of the quantities of interest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Cross-beam energy transfer in direct-drive ICF. II. Theory and simulation of mitigation through increased laser bandwidth

The response of the cross-beam energy transfer instability (CBET) to laser bandwidth is investigated through a combination of theory and simulation. Existing linear theory is generalized to treat broadband lasers, demonstrating that CBET is most effectively suppressed when the bandwidth exceeds the ion-acoustic wave (IAW) frequency. It is shown that for such bandwidths, reverse (seed to pump) transfer becomes possible, which reduces the net energy transfer rapidly as bandwidth is increased. The CBET gain exponent in this regime scales with bandwidth (Δω) as Δω−3 for Gaussian or Lorentzian laser spectra with different scalings possible for other spectra. Comparison of our theory with linearized fluid and particle-in-cell simulations, performed with the laser-plasma simulation environment (LPSE) and vector particle in cell (VPIC) codes, respectively, finds that the model is accurate in the absence of nonlinear processes. However, linear analysis also finds that the IAW energy density scales as Δω−1, implying that nonlinear effects may be more difficult to control than the CBET scaling would suggest. Indeed, nonlinear effects are found to be present in VPIC simulations with high-intensity lasers, despite minimal apparent CBET. Nonlinear processes in the VPIC cases include particle trapping, the two-ion wave decay, and ion wave self-focusing. In some high intensity VPIC cases, these effects lead to net energy transfer from seed to pump and increases to backscatter stimulated Brillouin scattering reflectivities. Finally, for a given bandwidth, we show that improved control of nonlinear processes can be achieved via smoothing by spectral dispersion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Hypercube technology

The JPL designed MARKIII hypercube supercomputer has been in application service since June 1988 and has had successful application to a broad problem set including electromagnetic scattering, discrete event simulation, plasma transport, matrix algorithms, neural network simulation, image processing, and graphics. Currently, problems that are not homogeneous are being attempted, and, through this involvement with real world applications, the software is evolving to handle the heterogeneous class problems efficiently.

Parker, Jay W.↗
DOE OSTI2016

Vlasov simulations of electron-ion collision effects on damping of electron plasma waves

Collisional effects can play an essential role in the dynamics of plasma waves by setting a minimum damping rate and by interfering with wave-particle resonances. Kinetic simulations of the effects of electron-ion pitch angle scattering on Electron Plasma Waves (EPWs) are presented here. In particular, the effects of such collisions on the frequency and damping of small-amplitude EPWs for a range of collision rates and wave phase velocities are computed and compared with theory. Both the Vlasov simulations and linear kinetic theory find the direct contribution of electron-ion collisions to wave damping significantly reduced from that obtained through linearized fluid theory. To our knowledge, this simple result has not been published before. Simulations have been carried out using a grid-based (Vlasov) approach, based on a high-order conservative finite difference method for discretizing the Fokker-Planck equation describing the evolution of the electron distribution function. Details of the implementation of the collision operator within this framework are presented. Such a grid-based approach, which is not subject to numerical noise, is of particular interest for the accurate measurements of the wave damping rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Numerical Improvements in Magnetohydrodynamic, Pulsed Power Simulations of Near-Target Plasmas

Magnetohydrodynamic (MHD) simulations of pulsed power experiments frequently result in unphysical runaway heating when extended to low values of mass density. Traditionally, this has been addressed by the use of conductivity floors below which the plasma is given an arbitrarily low conductivity value. Here, a low-density treatment is presented that allows for low-density material to carry current while maintaining a stable temperature. This treatment is implemented in the Ares multiphysics code. Numerical modifications to the conventional vacuum treatment include an energy-conserving density floor, a modified averaging procedure to determine the thermal conductivity at the edge of adjacent zones, and a modified averaging procedure for determining the thermal conductivity in mixed zones that contain material from differing regions. Additionally, anomalous resistivity and Bohm diffusion models are implemented as simplified models for microphysics-induced enhancement of collisional transport. Importantly, the advantage of these various improvements are illustrated in a simple 1-D pulsed power target where the combined changes result in stable temperatures within the lower density regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2000

Cloudy 94 and Applications to Quasar Emission Line Regions

This review discusses the most recent developments of the plasma simulation code Cloudy and its application to the, emission-line regions of quasars. The longterm goal is to develop the tools needed to determine the chemical composition of the emitting gas and the luminosity of the central engine for any emission line source. Emission lines and the underlying thermal continuum are formed in plasmas that are far from thermodynamic equilibrium. Their thermal and ionization states are the result of a balance of a vast set of microphysical processes. Once produced, radiation must, propagate out of the (usually) optically thick source. No analytic solutions are possible, and recourse to numerical simulations is necessary. I am developing the large-scale plasma simulation code Cloudy as an investigative tool for this work, much as an observer might build a spectrometer. This review describes the current version of Cloudy, version 94. It describes improvements made since the, release of the previous version, C90. The major recent, application has been the development of the "Locally Optimally-Emitting Cloud" (LOC) model of AGN emission line regions. Powerful selection effects, introduced by the atomic physics and line formation process, permit individual lines to form most efficiently only near certain selected parameters. These selection effects, together with the presence of gas with a wide range of conditions, are enough to reproduce the spectrum of a typical quasar with little dependence on details. The spectrum actually carries little information to the identity of the emitters. I view this as a major step forward since it provides a method to handle accidental details at the source, so that we can concentrate on essential information such as the luminosity or chemical composition of the quasar.

Ferland, Gary J.↗