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 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Numerical simulation of the interaction of the plasma sheet with the lobes of the Earth's magnetotail

Codes involving one and two spatial dimensions and three velocity dimensions were used to model the Earth's magnetotail. It was shown that the magnetotail can become inflated as a consequence of low energy plasma convection toward the neutral plane. The computer study exhibits a conversion of both magnetic field energy and of energy supplied by the convection electric field into particle energy. The numerical simulations suggest that much of the magnetotail substorm morphology may be a simple consequence of an increase, followed by a decrease, in the convection electric field, without the requirement of any magnetospheric size scale plasma instability or other disruptive processes. It is also concluded that the presence of the convection electric field and a continuing replenishment of low energy particles in the magnetotail are both necessary for maintenance of the magnetotail.

Swift, D. W.↗
NASA NTRS2007

Numerical Simulations of Flow Separation Control in Low-Pressure Turbines using Plasma Actuators

A recently introduced phenomenological model to simulate flow control applications using plasma actuators has been further developed and improved in order to expand its use to complicated actuator geometries. The new modeling approach eliminates the requirement of an empirical charge density distribution shape by using the embedded electrode as a source for the charge density. The resulting model is validated against a flat plate experiment with quiescent environment. The modeling approach incorporates the effect of the plasma actuators on the external flow into Navier Stokes computations as a body force vector which is obtained as a product of the net charge density and the electric field. The model solves the Maxwell equation to obtain the electric field due to the applied AC voltage at the electrodes and an additional equation for the charge density distribution representing the plasma density. The new modeling approach solves the charge density equation in the computational domain assuming the embedded electrode as a source therefore automatically generating a charge density distribution on the surface exposed to the flow similar to that observed in the experiments without explicitly specifying an empirical distribution. The model is validated against a flat plate experiment with quiescent environment.

Suzen, Y. B.↗
DOE OSTI2023

Modeling study of divertor particle flow pattern and in–out plasma density asymmetry due to drifts with SOLPS and BOUT++

A study of the effects of drifts on the particle flow pattern and in–out divertor plasma density asymmetry for L-mode and H-mode plasmas is carried out for EAST discharges by the edge plasma transport codes SOLPS and BOUT++. The simulation of L-mode plasmas is done by SOLPS while the simulation of H-mode plasmas is done by BOUT++. The toroidal magnetic field direction for the simulated discharge is artificially reversed in the codes to study the effects of different drift directions on the divertor particle flow pattern and the in–out asymmetry of divertor plasma density. The divertor particle flows induced by diamagnetic and E × B drifts are found to have similar directions in the divertor region for the same discharge. The directions of the flows induced by drifts would be reversed with the reverse of toroidal magnetic field direction. The diamagnetic drift seems to have no effect on the in–out asymmetry of divertor plasma density due to its divergence-free nature. However, the E × B drift could result in a pronounced asymmetry of plasma density between the inner and outer divertor targets. The density in–out asymmetry caused by E × B drift is reversed with the reverse of E × B drift flow direction. Detailed analysis shows that the radial component of the E × B drift flow is the main cause of density asymmetry. The results from the simulation of H-mode plasmas with BOUT++ are similar to those of the L-mode plasmas with SOLPS except that the drift effects seem to be slightly larger in the H-mode plasmas compared to the L-mode plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2008

Fast Plasma Instrument for MMS: Data Compression Simulation Results

Magnetospheric Multiscale (MMS) mission will study small-scale reconnection structures and their rapid motions from closely spaced platforms using instruments capable of high angular, energy, and time resolution measurements. To meet these requirements, the Fast Plasma Instrument (FPI) consists of eight (8) identical half top-hat electron sensors and eights (8) identical ion sensors and an Instrument Data Processing Unit (IDPU). The sensors (electron or ion) are grouped into pairs whose 6 deg x 180 deg fields-of-view (FOV) are set 90 deg apart. Each sensor is equipped with electrostatic aperture steering to allow the sensor to scan a 45 deg x 180 deg fan about its nominal viewing (0 deg deflection) direction. Each pair of sensors, known as the Dual Electron Spectrometer (DES) and the Dual Ion Spectrometer (DIS), occupies a quadrant on the MMS spacecraft and the combination of the eight electron/ion sensors, employing aperture steering, image the full-sky every 30-ms (electrons) and 150-ms (ions), respectively. To probe the results in the DES complement of a given spacecraft generating 6.5-Mbs(exp -1) of electron data while the DIS generates 1.1-Mbs(exp -1) of ion data yielding an FPI total data rate of 6.6-MBs(exp -1). The FPI electron/ion data is collected by the IDPU then transmitted to the Central Data Instrument Processor (CIDP) on the spacecraft for science interest ranking. Only data sequences that contain the greatest amount of temporal/spatial structure will be intelligently down-linked by the spacecraft. Currently, the FPI data rate allocation to the CIDP is 1.5-Mbs(exp -1). Consequently, the FPI-IDPU must employ data/image compression to meet this CIDP telemetry allocation. Here, we present simulations of the CCSDS 122.0-B-1 algorithm-based compression of the FPI-DES electron data. Compression analysis is based upon a seed of re-processed Cluster/PEACE electron measurements. Topics to be discussed include: review of compression algorithm; data quality; data formatting/organization; and, implications for data/matrix pruning. To conclude a presentation of the base-lined FPI data compression approach is provided.

Barrie, A.↗
DOE OSTI2023

Benchmarking magnetised three-wave coupling for laser backscattering: analytic solutions and kinetic simulations

Understanding magnetised laser–plasma interactions is important for controlling magneto-inertial fusion experiments and developing magnetically assisted radiation and particle sources. For nanosecond pulses at non-relativistic intensities, interactions are dominated by coherent three-wave interactions, whose nonlinear coupling coefficients became known only recently when waves propagate at oblique angles with the magnetic field. In this paper, backscattering coupling coefficients predicted by warm-fluid theory are benchmarked using particle-in-cell simulations in one spatial dimension, and excellent agreements are found for a wide range of plasma temperatures, magnetic field strengths and laser propagation angles, when the interactions are mediated by electron-dominant hybrid waves. Systematic comparisons between theory and simulations are made possible by a rigorous protocol. On the theory side, the initial boundary value problem of linearised three-wave equations is solved, and the transient-time solutions allow the effects of growth and damping to be distinguished. On the simulation side, parameters are carefully chosen and calibration runs are performed to ensure that comparisons are well controlled. Fitting simulation data to analytical solutions yields numerical growth rates that match theory predictions within error bars. Although warm-fluid theory is found to be valid for a wide parameter range, genuine kinetic effects have also been observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Thermodynamic and Kinetic Properties of Shocks in Two-Dimensional Yukawa Systems [Thermodynamic and Kinetic Properties of Shocks in 2D Yukawa Systems]

Particle-level simulations of shocked plasmas are carried out to examine kinetic properties not captured by hydrodynamic models. In particular, molecular dynamics simulations of 2D Yukawa plasmas with variable couplings and screening lengths are used to examine shock features unique to plasmas, including the presence of dispersive shock structures for weak shocks. A phase-space analysis reveals several kinetic properties, including anisotropic velocity distributions, non-Maxwellian tails, and the presence of fast particles ahead of the shock, even for moderately low Mach numbers. As a result, we also examine the thermodynamics (Rankine-Hugoniot relations) of recent experiments and find no anomalies in their equations of state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Dataset, Code, and Models for Training Deep Learning Potentials for Low Temperature Plasma-Surface Interactions

This repository contains datasets, training scripts, and finished models, and test simulations used in the development of DeepREBO— a machine-learned interatomic potential trained to emulate the REBO2 empirical potential. The data was generated to study deep potential development for simulations of plasma-surface interactions. It uses an active learning framework, starting from a minimal dataset and iteratively expanding it. Included are those generated datasets, the trained models, and simulations used to evaluate the performance of the training process. This resource supports reproducibility and provides a reference framework for training deep potentials in plasma-surface interaction studies.

active learning↗
DOE OSTI2024

Temporal evolution of the light emitted by a thin, laser-ionized plasma source

We present an experimental and simulation-based investigation of the temporal evolution of light emission from a thin, laser-ionized helium plasma source. We demonstrate an analytic model to calculate the approximate scaling of the time-integrated, on-axis light emission with the initial plasma density and temperature, supported by the experiment, which enhances the understanding of plasma light measurement for plasma wakefield accelerator (PWFA) plasma sources. Our model simulates the plasma density and temperature using a split-step Fourier code and a particle-in-cell code. A fluid simulation is then used to model the plasma and neutral density, and the electron temperature as a function of time and position. We then show the numerical results of the space-and-time-resolved light emission and that collisional excitation is the dominant source of light emission. We validate our model by measuring the light emitted by a laser-ionized plasma using a novel statistical method capable of resolving the nanosecond-scale temporal dynamics of the plasma light using a cost-effective camera with microsecond-scale timing jitter. This method is ideal for deployment in the high radiation environment of a particle accelerator that precludes the use of expensive nanosecond-gated cameras. Our results show that our models can effectively simulate the dynamics of a thin, laser-ionized plasma source. In addition, this work provides a detailed understanding of the plasma light measurement, which is one of the few diagnostic signals available for the direct measurement of PWFA plasma sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

A simulation study of interactions of Space-Shuttle generated electron beams with ambient plasma and neutral gas

The object was to conduct large scale simulations of electron beams injected into space. The study of active injection of electron beams from spacecraft is important since it provides valuable insight into beam-plasma interactions and the development of current systems in the ionosphere. However, the beam injection itself is not simple, being constrained by the ability of the spacecraft to draw return current from the ambient plasma. The generation of these return currents is dependent on several factors, including the density of the ambient plasma relative to the beam density, the presence of neutrals around the spacecraft, the configuration of the spacecraft, and the motion of the spacecraft through the plasma. Two dimensional particle simulations with collisional processes included are used to show how these different and often coupled processes can be utilized to enhance beam propagation from the spacecraft. To understand the radical expansion of mechanism of an electron beam from a highly charged spacecraft, two dimensional particle in cell simulations were conducted for a high density electron beam injected parallel to magnetic fields from an isolated equipotential conductor into a cold background plasma. The simulations indicate that charge buildup at the beam stagnation point causes the beam to expand radially to the beam electron gyroradius.

Source record↗
DOE OSTI2025

Global Gyrokinetic Simulations of Isotope Effects under Ambipolar Electric Fields and Advances Toward Whole-Volume Modeling

We review global gyrokinetic simulation studies on plasma transport in the Large Helical Device using XGC-S. XGC-S is an extended version of X-point Gyrokinetic Code for stellarators and has been progressively verified throughout the code development process. Verification tests of neoclassical transport successfully demonstrate the generation of an ambipolar electric field due to ripple-trapped particles. We perform quasi-linear analyses of the ion temperature gradient mode under the influence of the ambipolar electric field. The results reveal that the ambipolar electric field and the heavy hydrogen component in mixed isotope plasmas can lead to the favorable isotope effect observed in recent deuterium experiments. We also present recent efforts in code development toward whole-volume simulations, including the helical divertor region. A mesh generation scheme based on field-line tracing and the construction of curved surfaces perpendicular to the magnetic field would be promising for global field calculations in the whole-volume simulations.

Basic Plasma Phenomena and Gas Discharges↗
NASA NTRS1993

Studies on equatorial shock formation during plasmaspheric refilling

During the grant period starting August 1, 1992, our major effort has been on examining the presence of equatorially trapped hot plasma on plasmaspheric refilling. We performed one-dimensional PIC simulations of cold plasmas expanding into a hot plasma, consisting of hot anisotropic ions and warm isotropic electrons, trapped in a region of minimum magnetic field. Simulations showed that the electric potential barrier built up by the anisotropy of the hot ion population facilitates in the formation of electrostatic shocks when the cold ion beams begin to come into contact with the hot plasma. The shock formation occurs even when the cold ion beams are highly supersonic with respect to the ion-acoustic speed. This finding is interesting because equatorial shock formation during the early stage of plasmaspheric refilling has been debated over about two decades. In the past ion-ion instability has been invoked as the main mechanism for the coupling between the cold ion beams approaching the equator from the conjugate ionspheres. This coupling occurs when the beams are sufficiently slow; the beam velocity being less than three times the ion-acoustic speed. In the presence of hot plasma, the beams slow down by the potential barrier. The slowing down and the reflection process lead to the formation of the electrostatic shock even for highly supersonic ion beams. The mixing of hot and cold plasma was also studied.

Singh, Nagendra↗
DOE OSTI2024

Nonlinear gyrokinetic modelling of high confinement negative triangularity plasmas

Abstract Nonlinear gyrokinetic simulations correctly predict particle as well as ion and electron energy fluxes of high confinement plasmas with a negative triangularity cross sectional shape, showing that core transport in these plasmas is well described by standard gyrokinetic models. Experimentally inferred power balance fluxes are mostly reproduced within one standard deviation across a wide portion of the minor radius. Experimental conditions are reproduced by ion scale simulations, without the need to include density and temperature profile curvature effects. The experimental case is used as baseline to predict that the non-dimensional confinement scaling in negative triangularity plasmas increases strongly with plasma current while slightly degrading at increasing normalized pressure and decreasing collisionality. Recent experiments showed that low toroidal rotation negatively impacts confinement; consistent with the experiment, simulations predict that low rotational shear significantly affects confinement unless the plasma effective charge is maintained above a minimum level. Core confinement is predicted to significantly degrade in low aspect ratio devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Numerical simulation to study the transient self focusing of laser beam in plasma

In this paper, we present the numerical simulation for the coupled system of equations governing the dynamics of laser and Ion Acoustic Wave (IAW) in a collisionless plasma, when the coupling between the waves is through ponderomotive non-linearity. The nonlinear evolution of the laser beam is studied when the pump laser is perturbed by a periodic perturbation. By changing the perturbation wave number, we have studied its effect on the nonlinear evolution pattern of laser beam. In order to have a physical insight into the nonlinear dynamics of laser beam evolution in time and space, we have studied the laser and IAW spectra containing spatial harmonics. The magnitude of these harmonics changes with time and leads to time dependent localization of laser beam in spatial domain. The nonlinear dynamics of this localization is investigated in detail by using simulation and a semi-analytical model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Assessing time-dependent temperature profile predictions using reduced transport models for high performing NSTX plasmas

Time-dependent, predictive simulations were performed with the 1.5D tokamak integrated modeling code TRANSP on a large set of well-analyzed, high performing discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate how well modern reduced transport models can reproduce experimentally observed temperature profiles in spherical tokamaks. Overall, it is found that simulations using the Multi-Mode Model (MMM) more consistently agree with the NSTX observations than those using the Trapped Gyro-Landau Fluid (TGLF) model, despite TGLF requiring orders of magnitude greater computational cost. When considering all examined discharges, MMM has median overpredictions of electron temperature (T e ) and ion temperature (T i ) profiles of 28% and 27%, respectively, relative to the experiment. TGLF overpredicts T e by 46%, with much larger variance than MMM, and underpredicts T i by 25%. As the ratio of kinetic to magnetic field pressure (β) is increased across NSTX discharges, TGLF predicts lower T e and significant flattening of the T i profile, conflicting with NSTX observations. When using an electrostatic version of TGLF, both T e and T i are substantially overpredicted, underscoring the importance of electromagnetic turbulence in the high β spherical tokamak regime. Additionally, calculations with neural net surrogate models for TGLF were performed outside of TRANSP with a time slice flux matching transport solver, finding better agreement with experiment than the TRANSP simulations, highlighting the impact of different transport solvers and simulation techniques. Furthermore, the reasonable agreement with experiment of temperature profiles predicted by MMM motivates a more detailed examination of the sensitivities of the TRANSP simulations with MMM to different NSTX plasma regimes in a companion paper (Lestz et al Plasma Phys. Control. Fusion 67 105030), in preparation for self-consistent, time-dependent predictive modeling of NSTX-U scenarios.

Lestz, J. B. [General Atomics, San Diego, CA (Unit↗
DOE OSTI2023

Photon Doppler Velocimetry to Spatially Resolve Plasma Density in a Power Flow Gap

The understanding of power flow plasmas is important as we look towards next generation pulsed power (NGPP) as current losses could prohibit the goals of that facility. Therefore, it is important to have accurate diagnostics of the plasma parameters on the current machines, which can be used to help inform and improve simulations. Having these plasma parameters will help validate models and simulations to provide confidence when they are expanded to conditions relevant to NGPP. One important plasma parameter that can be measured is the electron density, which can be measured by photonic Doppler velocimetry (PDV). A PDV system has several key advantages over other interferometers by measuring relatively low densities (> 1 × 10 15 cm -2 ) with both spatial and temporal resolution. Experiments were performed on the Mykonos pulsed power machine, which is a 1 MA sub scale machine in which recent platforms have been developed to explore current densities relevant to the inner magnetically insulated transmission line (MITL) on the Z machine. Experiments were performed on two different platforms, the thin foil platform and the Mykonos parallel plate platform (MP3). In addition, a combination of both single-point and multi-point measurements were used. The single-point measurements proved to be very promising, providing a clear increase in density at about 70 ns into the current rise on thin foil experiments up to about 5 × 10 17 cm -3 before the probe stopped providing signal. While we did also see returns from multi-point measurements on both platforms, the signals were not as easy to interpret due to strong background effects. However, they do show initial promise for this diagnostic to measure density at several points across a 1 mm gap. These measurements provide insights in how to improve the diagnostic so that it can provide useful information on power flow relevant experiments.

47 OTHER INSTRUMENTATION↗
NASA NTRS1988

Simulation of electrostatic turbulence in the plasma sheet boundary layer with electron currents and bean-shaped ion beams

Plasma data from ISEE-1 show the presence of electron currents as well as energetic ion beams in the plasma sheet boundary layer. Broadband electrostatic noise and low-frequency electromagnetic bursts are detected in the plasma sheet boundary layer, especially in the presence of strong ion flows, currents, and steep spacial gradients in the fluxes of few-keV electrons and ions. Particle simulations have been performed to investigate electrostatic turbulence driven by a cold electron beam and/or ion beams with a bean-shaped velocity distribution. The simulation results show that the counterstreaming ion beams as well as the counterstreaming of the cold electron beam and the ion beam excite ion acoustic waves with a given Doppler-shifted real frequency. However, the effect of the bean-shaped ion velocity distributions reduces the growth rates of ion acoustic instability. The simulation results also show that the slowing down of the ion bean is larger at the larger perpendicular velocity. The wave spectra of the electric fields at some points of the simulations show turbulence generated by growing waves.

Nishikawa, K.-I.↗
DOE OSTI2016

Efficient production of high-energy nonthermal particles during magnetic reconnection in a magnetically dominated ion-electron plasma

Magnetic reconnection is a leading mechanism for dissipating magnetic energy and accelerating nonthermal particles in Poynting-flux-dominated flows. We investigate nonthermal particle acceleration during magnetic reconnection in a magnetically dominated ion–electron plasma using fully kinetic simulations. For an ion–electron plasma with a total magnetization of ${\sigma }_{0}={B}^{2}/(4\pi n({m}_{i}+{m}_{e}){c}^{2})$, the magnetization for each species is ${\sigma }_{i}\sim {\sigma }_{0}$ and ${\sigma }_{e}\sim ({m}_{i}/{m}_{e}){\sigma }_{0}$, respectively. We have studied the magnetically dominated regime by varying σe = 103–105 with initial ion and electron temperatures ${T}_{i}={T}_{e}=5-20{m}_{e}{c}^{2}$ and mass ratio ${m}_{i}/{m}_{e}=1-1836$. Our results demonstrate that reconnection quickly establishes power-law energy distributions for both electrons and ions within several (2–3) light-crossing times. For the cases with periodic boundary conditions, the power-law index is $1\lt s\lt 2$ for both electrons and ions. The hard spectra limit the power-law energies for electrons and ions to be ${\gamma }_{{be}}\sim {\sigma }_{e}$ and ${\gamma }_{{bi}}\sim {\sigma }_{i}$, respectively. The main acceleration mechanism is a Fermi-like acceleration through the drift motions of charged particles. When comparing the spectra for electrons and ions in momentum space, the spectral indices sp are identical as predicted in Fermi acceleration. We also find that the bulk flow can carry a significant amount of energy during the simulations. Finally, we discuss the implication of this study in the context of Poynting-flux dominated jets and pulsar winds, especially the applications for explaining nonthermal high-energy emissions.

79 ASTRONOMY AND ASTROPHYSICS↗