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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

NUMERICAL SIMULATIONS OF PLASMA DYNAMICS IN THE VICINITY OF A RETRACTING FLUX TUBE

In a previous paper, we presented an analytical, zero- β model for supra-arcade downflows in which a retracting flux tube deforms the surrounding magnetic field, constricting the flow of plasma along affected field lines and, in some cases, forcing the plasma to exhibit collimated shocks. Here we present a numerical simulation based on the same model construction—a retracting flux tube is treated as a rigid boundary around which the plasma is forced to flow and the magnetic field and plasma evolve according to the governing equations of magnetohydrodynamics. We find that the collimated shocks described in our previous study are recovered for plasma β in the range of 0 ≤ β ≲ 1, while for 1 ≲ β the behavior is similar to the simpler hydrodynamic case, with classical bow shocks forming when the acoustic Mach number approaches or exceeds unity. Furthermore, we find that while the plasma β is important for identifying the various types of behaviors, more important still is the Alfvén Mach number, which, if large, implies that the bulk kinetic energy of the fluid exceeds the internal energy of the magnetic field, thereby leading to the formation of unconfined, fast-mode magnetosonic shocks, even in the limit of small β .

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1990

Plasma contactor research, 1989

The characteristics of double layers observed by researchers investigating magnetospheric phenomena are contrasted to those observed in plasma contacting experiments. Experiments in the electron collection mode of the plasma contacting process were performed and the results confirm a simple model of this process for current levels ranging to 3 A. Experimental results were also obtained in a study of the process of electron emission from a hollow cathode plasma contactor. High energy ions are observed coming from the cathode in addition to the electrons and a phenomenological model that suggests a mechanism by which this could occur is presented. Experimental results showing the effects of the design parameters of the ambient plasma simulator on the plasma potential, electron temperature, electron density and plasma noise levels induced in plasma contacting experiments are presented. A preferred simulator design is selected on the basis of these results.

Williams, John D.↗
DOE OSTI2017

Fast low-to-high confinement mode bifurcation dynamics in a tokamak edge plasma gyrokinetic simulation

Transport barrier formation and its relation to sheared flows in fluids and plasmas are of fundamental interest in various natural and laboratory observations and of critical importance in achieving an economical energy production in a magnetic fusion device. Here we report the first observation of an edge transport barrier formation event in an electrostatic gyrokinetic simulation carried out in a realistic diverted tokamak edge geometry under strong forcing by a high rate of heat deposition. Here, the results show that turbulent Reynolds-stress-driven sheared E x B flows act in concert with neoclassical orbit loss to quench turbulent transport and form a transport barrier just inside the last closed magnetic flux surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Numerical simulation of plasma processes driven by transverse ion heating

The plasma processes driven by transverse ion heating in a diverging flux tube are investigated with numerical simulation. The heating is found to drive a host of plasma processes, in addition to the well-known phenomenon of ion conics. The downward electric field near the reverse shock generates a doublestreaming situation consisting of two upflowing ion populations with different average flow velocities. The electric field in the reverse shock region is modulated by the ion-ion instability driven by the multistreaming ions. The oscillating fields in this region have the possibility of heating electrons. These results from the simulations are compared with results from a previous study based on a hydrodynamical model. Effects of spatial resolutions provided by simulations on the evolution of the plasma are discussed.

Singh, Nagendra↗
DOE OSTI2018

Community Project for Accelerator Science and Simulation-3

The UCLA Plasma Simulation Group is a major partner of the “Community Petascale Project for Accelerator Science and Simulation”. This is the technical progress report for the UCLA effort. The UCLA effort is organized into HEP and ASCR activities.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

A generative artificial intelligence framework for long-time plasma turbulence simulations

Generative deep learning techniques are employed in a novel framework for the construction of surrogate models capturing the spatiotemporal dynamics of 2D plasma turbulence. The proposed Generative Artificial Intelligence Turbulence (GAIT) framework enables the acceleration of turbulence simulations for long-time transport studies. GAIT leverages a convolutional variational auto-encoder and a recurrent neural network to generate new turbulence data from existing simulations, extending the time horizon of transport studies with minimal computational cost. The application of the GAIT framework to plasma turbulence using the Hasegawa–Wakatani (HW) model is presented, evaluating its performance via various analyses. Very good agreement is found between the GAIT and the HW models in the spatiotemporal Fourier and Proper Orthogonal Decomposition spectra, the flow topology characterized by the Okubo–Weiss parameter, and the time autocorrelation function of turbulent fluctuations. Excellent agreement has also been obtained in the probability distribution function of particle displacements and the effective turbulent diffusivity. In-depth analyses of the latent space of turbulent states, choice of hyperparameters and alternative deep learning models for the time prediction are presented. Our results highlight the potential of Artificial Intelligence-based surrogate models to overcome the computational challenges in turbulence simulation, which can be extended to other situations such as geophysical fluid dynamics.

Artificial intelligence↗
DOE OSTI2023

An extended Vlasov–Fokker–Planck approach for kinetic simulations of laser plasmas

Vlasov–Fokker–Planck simulation codes occupy an important niche in modeling laser-produced plasmas, since they are well suited to studying the effect of collisions on electron kinetic phenomena, especially energy transport. One of the most important elements of energy transport is the absorption of laser light by the plasma; however, simulating this in detail requires resolving oscillations of the laser light, whose characteristic timescale is orders of magnitude shorter than the simulation time needed to study transport physics. For this reason, most Vlasov–Fokker–Planck codes used to study electron transport in laser plasmas rely on simplified models of the laser–plasma coupling. Their underlying assumptions nominally preclude their use for modeling laser light having short-scale structure in space or time, such as broadband lasers. In this work, we derive a more general computational framework suitable for arbitrarily structured laser fields. Furthermore, our approach is based on an extended set of Vlasov–Fokker–Planck equations that separately solve for the low- and high-frequency plasma response. We implement these extended Vlasov–Fokker–Planck equations in the spherical harmonic code K2 and demonstrate the performance of the method on several laser absorption test problems, with particular attention to the judicious selection of time steps, time integrators, and spherical harmonic truncation, according to the intensity and spectrum of the laser light under consideration. Comparison with the widely used Langdon absorption operator shows the Langdon operator performs remarkably well for predicting laser heating in the simple cases considered here, even in situations that would seem to violate its underlying assumptions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Quantum signal processing for simulating cold plasma waves

Numerical modeling of radio-frequency waves in plasma with sufficiently high spatial and temporal resolution remains challenging even with modern computers. However, such simulations can be sped up using quantum computers in the future. In this work, we propose how to do such modeling for cold plasma waves, in particular, for an X wave propagating in an inhomogeneous one-dimensional plasma. The wave system is represented in the form of a vector Schrödinger equation with a Hermitian Hamiltonian. Block encoding is used to represent the Hamiltonian through unitary operations that can be implemented on a quantum computer. To perform the modeling, we apply the so-called quantum signal processing algorithm and construct the corresponding circuit. Quantum simulations with this circuit are emulated on a classical computer, and the results show agreement with traditional classical calculations. We also discuss how our quantum circuit scales with the resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Importance of Resolving the Spectral Support of Beam-plasma Instabilities in Simulations

Many astrophysical plasmas are prone to beam-plasma instabilities. For relativistic and dilute beams, the spectral support of the beam-plasma instabilities is narrow, i.e., the linearly unstable modes that grow with rates comparable to the maximum growth rate occupy a narrow range of wavenumbers. This places stringent requirements on the box-sizes when simulating the evolution of the instabilities. We identify the implied lower limits on the box size imposed by the longitudinal beam plasma instability, i.e., typically the most stringent condition required to correctly capture the linear evolution of the instabilities in multidimensional simulations. We find that sizes many orders of magnitude larger than the resonant wavelength are typically required. Using one-dimensional particle-in-cell simulations, we show that the failure to sufficiently resolve the spectral support of the longitudinal instability yields slower growth and lower levels of saturation, potentially leading to erroneous physical conclusion.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1984

Turbulent heating of heavy ions on auroral field lines

Electrostatic ion cyclotron turbulence and the associated acceleration of ions on auroral field lines are investigated analytically and by plasma simulations. The auroral plasma is assumed to consist of drifting electrons and stationary hydrogen and oxygen ions. It is found that for a given critical drift, the maximum ion perpendicular heating is generally larger for oxygen ions than for hydrogen ions. Simulation results show that unless the oxygen ions are a minority species, oxygen transverse heating generally exceeds that of hydrogen ions. Theory and numerical simulations are in good agreement.

Ashour-Abdalla, M.↗
DOE OSTI2018

Extension of Landau-fluid closure to weakly collisional plasma regime

Both collisionless and collisional Landau fluid (LF) closure have been developed and implemented in two-fluid plasma simulation framework BOUT++, taking full advantage of the fast non-Fourier method (Dimits, 2014). The spectral range of good fit can be conveniently extended to larger domain by adding more Lorentzians. The scaling for upper limit of wavenumber resolved in collisionless LF operator and the fitting coefficients for the collisional LF operator, the later is extended to weakly collisional regime (ν$^{*}_{e}$ ~ 0.01), are presented in this paper. The flux limited expression q$^{FL}_{∥e}$ is reduced to free streaming expression q$^{FS}_{∥e}$ in collisionless limit and gives classical Spitzer–Härm result q$^{SH}_{∥e}$ in collisional limit as we expected. The collisional LF closure recovers the classical Spitzer–Härm result q$^{SH}_{∥e}$ in collisional limit and is same as the collisionless LF closure in collisionless limit. The test cases also obviously reflect the nonlocal effects from LF closures and responding to different boundary conditions in open surface region. Our report will make it more applicable and possible to include the Landau damping kinetic effect in the fluid simulation of tokamak plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2021

Terminator Double Layer Explorer (TerDLE): Examining the Near-Moon Lunar Wake

As the solar wind flows by the Moon, an antisunward-directed low-density wake forms as the plasma expands to fill in the trailing void in the plasma flow. Analytical modeling and modern plasma simulations suggest that plasma quasi-neutrality could possibly be broken close to the terminator obstruction as solar wind electrons expand into the wake ahead of the ions, leading to the formation of a standing (time-stationary) double layer. The objective of the Terminator Double Layer Explorer is to extend the fundamental understanding of the plasma expansion into the trailing near-vacuum wake region by (1) identifying any plasma expansion density anomalies at low altitudes near the terminator wake initiation region, (2) assessing the highly variable solar wind's effect on the low-altitude wake region, and (3) determining if plasma neutrality is maintained or lost during passages through the low-altitude expansion region. The mission concept uses a propulsion-driven CubeSat with ion spectrometer and plasma wave system in elliptical orbit about the Moon with periselene near the terminator. Over the course of the mission, the periselene decreases, placing the CubeSat ever closer to the terminator wake initiation location and the possible nonneutral region.

lunar science↗
DOE OSTI2022

Simulation of Plasma Emission in Magnetized Plasmas

The recent Parker Solar Probe observations of type III radio bursts show that the effects of the finite background magnetic field can be an important factor in the interpretation of data. In the present paper, the effects of the background magnetic field on the plasma-emission process, which is believed to be the main emission mechanism for solar coronal and interplanetary type III radio bursts, are investigated by means of the particle-in-cell simulation method. The effects of the ambient magnetic field are systematically surveyed by varying the ratio of plasma frequency to electron gyrofrequency. The present study shows that for a sufficiently strong ambient magnetic field, the wave–particle interaction processes lead to a highly field-aligned longitudinal mode excitation and anisotropic electron velocity distribution function, accompanied by a significantly enhanced plasma emission at the second-harmonic plasma frequency. For such a case, the polarization of the harmonic emission is almost entirely in the sense of extraordinary mode. On the other hand, for moderate strengths of the ambient magnetic field, the interpretation of the simulation result is less clear. The underlying nonlinear-mode coupling processes indicate that to properly understand and interpret the simulation results requires sophisticated analyses involving interactions among magnetized plasma normal modes, including the two transverse modes of the magneto-active plasma, namely, the extraordinary and ordinary modes, as well as electron-cyclotron-whistler, plasma oscillation, and upper-hybrid modes. At present, a nonlinear theory suitable for quantitatively analyzing such complex-mode coupling processes in magnetized plasmas is incomplete, which calls for further theoretical research, but the present simulation results could provide a guide for future theoretical efforts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Synthetic diagnostics platform for fusion plasmas (invited)

A Synthetic Diagnostics Platform (SDP) for fusion plasmas has been developed which provides state of the art synthetic reflectometry, beam emission spectroscopy, and Electron Cyclotron Emission (ECE) diagnostics. Interfaces to the plasma simulation codes GTC, XGC-1, GTS, and M3D-C{sup 1} are provided, enabling detailed validation of these codes. In this paper, we give an overview of SDP’s capabilities, and introduce the synthetic diagnostic modules. A recently developed synthetic ECE Imaging module which self-consistently includes refraction, diffraction, emission, and absorption effects is discussed in detail. Its capabilities are demonstrated on two model plasmas. The importance of synthetic diagnostics in validation is shown by applying the SDP to M3D-C{sup 1} output and comparing it with measurements from an edge harmonic oscillation mode on DIII-D.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1988

An MHD simulation of plasma flow past Io - Alfven and slow mode perturbations

The flow of plasma past Io is studied using a time dependent, three-dimensional magnetohydrodynamic (MHD) simulation. In addition to observing the 'Alfven wing', a standing Alfven wave perturbation expected from analytic theory, it is found that the other MHD modes contribute important perturbations. Standing slow mode perturbations also are present in the flow.

Linker, Jon A.↗