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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Measurements and Simulations of Surface Dielectric Barrier Discharges Used as Plasma Actuators

This report is a Ph.D. dissertation performed under NRA cooperative agreement and submitted as part of the final report. Asymmetric surface dielectric barrier discharges (DBDs) have shown promise for use as aerodynamic actuators for active flow control. In this project we studied DBD actuators experimentally and numerically. Our DBDs used a symmetric triangular high voltage waveform to generate plasma in atmospheric pressure air. Time-averaged measurements indicated that the induced force of a single barrier actuator design (one electrode insulated from the plasma) can be increased exponentially above the results of previous studies by decreasing both the length and thickness of the electrode exposed to the plasma. This increased force may allow these devices to control flow separation in a wider range of flow environments. Experiments using an intensified digital camera to examine the plasma on time scales of a few nanoseconds showed that, in addition to the previously-observed filamentary and jet-like plasma structures, discharges with very thin exposed electrodes exhibited a weak but constant plasma immediately adjacent to those electrodes. In double-barrier actuators (both electrodes insulated), decreasing the diameter of the narrower electrode lead to increasing forces, and recorded images showed the simultaneous existence of both filamentary and jet-like plasma structures. The development and application of a time-dependent, two-dimensional computational fluid plasma model has aided in understanding the detailed physics of surface DBDs at all-time scales. For simulated single-barrier discharges, the model qualitatively reproduced the filamentary and jet-like micro-discharge structures. The model was somewhat successful in reproducing the observed characteristics of double-barrier actuators. For both actuator geometries, the model indicated that the majority of the forces induced on the neutral gas occur in between micro-discharges as the plasmas decay.

Hoskinson, Alan R.↗
DOE OSTI2018

Design and Analysis of an Actively Cooled Window for a High-Power Helicon Plasma Source

The material's plasma exposure experiment (MPEX) is a steady-state, linear plasma facility proposed to provide plasma conditions that simulate the first wall and divertor of the next-step fusion devices. MPEX will use an actively cooled, high-power helicon antenna as a plasma source. The antenna is located outside of the vacuum region, and during operation, a significant percentage of the antenna power is deposited on the plasma facing surface of the antenna window. The window thus must be actively cooled so that thermal stresses do not become excessive. The design presented includes two cylindrical windows with an annular water passage in-between in order to allow active cooling. A ceramic-to-metal seal is proposed in order to produce a water-to-vacuum seal that will not outgas and produce impurities. Experimental results from an uncooled window on a short-pulse prototype experiment are modeled using finite-element (FE) results in order to correlate with the heat flux profile. The design of the actively cooled window assembly is presented with computational fluid dynamics and FE analyses to confirm that the hardware can function with a 200-kW antenna within performance limits.

Lumsdaine, Arnold↗
NASA NTRS1986

Numerical simulation of spacecraft charging by impact-induced plasmas during a cometary flyby

A numerical model is developed for the interaction of a cometary probe, such as Giotto, with its environment, i.e., dust particles and gas. The spacecraft was set on a course to pass the comet at a velocity of 69 km/sec, so a chance existed that a potential field would form around the spacecraft and block lower energy particles from reaching the spacecraft instruments. The motion of electrons and ions is traced as a function of time to examine the evolution of the electric field, electric potential and the total space charge distributions on the surface of the spacecraft and its environment. Account is taken of the density of the particles and gas molecules at various distances from the comet, the collision energies involved, and the Giotto geometry. A solution is defined for the Poisson equation to describe the evolution of the plasma around Giotto, including the effects of ion collisions with the Al bumper protecting the spacecraft. The simulation predicts formation of an ion wake behind Giotto and the evolution of a positive potential on the order of 10 V around the spacecraft, i.e., sufficient for a positive potential barrier near the surface of the spacecraft.

Thiemann, H.↗
NASA NTRS2009

Combined Modeling of Acceleration, Transport, and Hydrodynamic Response in Solar Flares: The Numerical Model - 1

Acceleration and transport of high-energy particles and fluid dynamics of atmospheric plasma are interrelated aspects of solar flares, but for convenience and simplicity they were artificially separated in the past. We present here self consistently combined Fokker-Planck modeling of particles and hydrodynamic simulation of flare plasma. Energetic electrons are modeled with the Stanford unified code of acceleration, transport, and radiation, while plasma is modeled with the Naval Research Laboratory flux tube code. We calculated the collisional heating rate directly from the particle transport code, which is more accurate than those in previous studies based on approximate analytical solutions. We repeated the simulation of Mariska et al. with an injection of power law, downward-beamed electrons using the new heating rate. For this case, a -10% difference was found from their old result. We also used a more realistic spectrum of injected electrons provided by the stochastic acceleration model, which has a smooth transition from a quasi-thermal background at low energies to a non thermal tail at high energies. The inclusion of low-energy electrons results in relatively more heating in the corona (versus chromosphere) and thus a larger downward heat conduction flux. The interplay of electron heating, conduction, and radiative loss leads to stronger chromospheric evaporation than obtained in previous studies, which had a deficit in low-energy electrons due to an arbitrarily assumed low-energy cutoff. The energy and spatial distributions of energetic electrons and bremsstrahlung photons bear signatures of the changing density distribution caused by chromospheric evaporation. In particular, the density jump at the evaporation front gives rise to enhanced emission, which, in principle, can be imaged by X-ray telescopes. This model can be applied to investigate a variety of high-energy processes in solar, space, and astrophysical plasmas.

STOCHASTIC ACCELERATION MODEL↗
NASA NTRS1972

Collisionless plasma flow over a conducting sphere.

Laboratory simulation study of spacecraft-space plasma interaction using a plasma wind tunnel. Some preliminary results on a collisionless plasma flow pattern over a conducting sphere are presented. A physical explanation is given for the observed plasma flow behavior.

Stone, N. H.↗
DOE OSTI2016

Spontaneous focusing of plasma flow in a weak perpendicular magnetic field

Structure formation of high-beta plasma flow in a perpendicular magnetic field is investigated in the ion kinetic regime by a fully kinetic particle-in-cell simulation. We demonstrate that directional plasma flow is spontaneously focused to form a sharp density structure. The primary focusing process comes from field-aligned electron inflow associated with the whistler mode and plasma confinement due to a self-generated magnetic field. The resulting concave magnetic field lines modulate ion gyration to cause a secondary focusing process with significant plasma concentration. Required conditions for these processes are determined by a dimensionless parameter α ≡ β{sub i0}(Δ{sub W}/ρ{sub i0}), where β{sub i0}, Δ{sub W}, and ρ{sub i0} denote the plasma kinetic beta, window size, and ion gyration radius, respectively. The focusing process is apparent for small α, whereas diamagnetic expansion is dominant for large α. This condition describes a transition between diamagnetic cavity formation and the focusing process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Coupled EM-PIC/Radiation Transport Simulations of HERMES Courtyard Experiments

A suite of coupled computational models for simulating the radiation, plasma, and electromagnetic (EM) environment in the High-Energy Radiation Megavolt Electron Source (HERMES) courtyard has been developed. In principle, this provides a predictive forward-simulation capability based solely on measured upstream anode and cathode current waveforms in the Magnetically Insulated Transmission Line (MITL). First, 2D R-Z ElectroMagnetic Particle-in-Cell (EM-PIC) simulations model the MITL and diode to compute a history of all electrons incident on the converter. Next, radiation transport simulations use these electrons as a source to compute the time-dependent dose rate and volumetric electron production in the courtyard. Finally, the radiation transport output is used as sources for EM-PIC simulations of the courtyard to com- pute electromagnetic responses. This suite has been applied to the November 2016 trials, shots 10268-10313. Modeling and experiment differ in significant ways. This is just the first iteration of a long process to improve the agreement, as outlined in the summary.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Modeling of a 43 GHz microplasma within a photonic crystal

The electron density (n e ) and argon metastable density (1 s 5 ) of a 43GHz microplasma are obtained using the zero-dimensional plasma kinetics solver (ZDPlasKin®) for an incident power of 160 and 1000mW and argon pressures of 10–600Torr (1.3 ×10 3 to 8.0 × 10 4 Pa). These simulations are compared with previously published experimental data. To determine the self-consistent electric field in the modeled plasma, the three-dimensional millimeter wave fields are computed as a function of electron density using ANSYS EM19.2, HFSS®. This electromagnetic field model is coupled to ZDPlasKin such that any increase in the simulated plasma density correctly attenuates the simulated electric field within the microplasma. The electron density is found to be sensitive to argon gas temperature, so a two-zone temperature model was needed to obtain agreement with experimental measurements. The temperature in the central core of the microplasma was determined by a previous experimental study. That temperature was used as an input to the model for the simulation of volume recombination losses. The outer regions of the microplasma are assumed to be in equilibrium with the walls (300K). This second temperature was used in the model to determine diffusion losses. The modeled electron and metastable densities are of the order 10 20 and 10 18 m –3 , respectively. This is in good agreement with those measured experimentally as long as the two-zone temperature model is used. In the absence of a hot gas core, the modeled three-body recombination rates are excessive and the simulation severely under-predicts the electron density and over-estimates the metastable density. We conclude that the millimeter wave microplasma has a hot core (2500K at 600 Torr) that rarifies the argon gas and effectively reduces three-body recombination. Furthermore, this allows one to achieve high electron density on the order of 10 20 m –3 with only 100mW of wave power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Laser-plasma instabilities in long scale-length plasmas relevant to shock-ignition

We present 2D particle-in-cell simulations of laser-plasma instabilities (LPI) performed using conditions relevant to experiments on the OMEGA and NIF laser facilities. The laser intensity used is 2×10 15 W cm -2 , which allows comparison with previous experiments and simulations investigating shock ignition, though this is a relatively low intensity for the scheme. We find that the large convective gain of the instabilities leads to dynamics that are largely controlled by pump depletion and to scattering or absorption by LPIs predominantly occurring away from the quarter-critical density. The two-plasmon decay (TPD) instability is dominant for OMEGA-type conditions and remains important at the NIF-scale. In both cases, most absorption from TPD occurs near its Landau cutoff density, and we examine its nonlinear dynamics in this region, focusing on the ion-acoustic waves responsible for saturation. For NIF-type conditions, stimulated Raman scattering (SRS) and stimulated Brillouin scattering divert significant fractions of laser energy, and we find that SRS backscatter can occur with high gain at low densities due to kinetic effects. An in-depth analysis of hot-electron production is performed to identify the main sources and their characteristic electron temperature. For both simulations presented, the overall hot-electron temperature lies between 30 and 35 keV. This is relatively low and will likely increase shock pressure; however, the distributions also contain a significant number of high-energy (E k >100 keV) electrons that would likely cause unacceptable preheat. We, therefore, suggest some strategies that may be utilized to minimize this high-energy component. I. INTROD

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Classification of Orbits in Poincare Maps Using Machine Learning

Poincare plots, also called Poincare maps, are used by plasma physicists to understand the behavior of magnetically confined plasma in numerical simulations of a tokamak. These plots are created by the intersection of field lines with a two-dimensional poloidal plane that is perpendicular to the axis of the torus representing the tokamak. A plot is composed of multiple orbits, each created by a different field line as it goes around the torus. Each orbit can have one of four distinct shapes, or classes, that indicate changes in the topology of the magnetic fields confining the plasma. Given the (x, y) coordinates of the points that form an orbit, the analysis task is to assign a class to the orbit, a task that appears ideally suited for a machine learning approach. In this paper, we describe how we overcame two major challenges in solving this problem - creating a high-quality training set, with few mislabeled orbits, and converting the coordinates of the points into features that are discriminating, despite the variation within the orbits of a class and the apparent similarities between orbits of different classes. Our automated approach is not only more objective and accurate than visual classification, but is also less tedious, making it easier for plasma physicists to analyze the topology of magnetic fields from numerical simulations of the tokamak.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Analytical and numerical studies of dark current in radiofrequency structures for short-pulse high-gradient acceleration

High-gradient acceleration is a key research area that could enable compact linear accelerators for future colliders, light sources, and other applications. In the pursuit of high-gradient operation, rf breakdown limits the attainable accelerating gradient in normal-conducting rf structures. Recent experiments at the Argonne Wakefield Accelerator suggest a promising approach: using short rf pulses with durations of a few nanoseconds. Experimental studies show that these 𝒪⁡(1⁢ ⁢ns) rf pulses can mitigate breakdown limitations, resulting in higher gradients. For example, an electric field of nearly 400 MV/m was achieved in an 𝑋-band photoemission gun driven by 6-ns-long rf pulses, with rapid rf conditioning and low dark current observed. Despite these promising results, the short-pulse regime remains an underexplored parameter space, and rf breakdown physics under nanosecond-long pulses requires further investigation. In this paper, we present analytical and numerical simulations of dark current dynamics in accelerating cavities operating in the short-pulse regime. We study breakdown-associated processes spanning different time scales, including field emission, multipacting, and plasma formation, using simulations of the 𝑋-band photogun cavities. The results reveal the advantages of using short rf pulses to reduce dark current and mitigate rf breakdown, offering a path toward a new class of compact accelerators with enhanced performance and reduced susceptibility to breakdown.

Linear accelerators↗
DOE OSTI2021

Modelling dust transport in DIII-D with DTOKS-Upgrade

Comprehensive upgrades to the dust transport code Dust in TOKamaks (DTOKS) that extend the plasma-dust interaction model are presented and compared with recent measurements of dust transport in DIII-D. Simulations incorporating variation in physical properties of graphite dust with temperature and size in a stationary plasma background suggest a substantial decrease in lifetimes due principally to thermal expansion. The trajectories of 53 dust grains identified from analysis of visible camera data taken across two similar shots were used to measure the dust particle velocity distributions. Dust tracks terminated mostly at the outer divertor strike point having a mean observation time of 2.1 ± 0.4 ms. Stochastic modelling of 200 graphite dust particles in the DIII-D tokamak performed with DTOKS-U using plasma simulations generated by OEDGE found similar behaviour, with particles ablating rapidly after acquiring a positive charge in the region close to the outer strike point, creating an acute source of neutral carbon atoms. The simulated mean lifetime, 11 ± 2 ms, showed approximate agreement with experimental observation when corrected by accounting for dust visibility and ignoring the longest trajectories 1.5 ± 0.2 ms. Synthetic diagnostic data generated from coupling the results of DTOKS-U with the visualisation software Calcam offers a powerful new tool for validation of simulations and predictive calculations of dust dynamics.

Physics↗
DOE OSTI2025

Simulations of Radiation Hydrodynamics in Strongly Coupled Plasmas (Final Report)

Radiation propagates throughout the Universe and is often the only information we have from distant astronomical objects. Additionally, radiation plays a key role in systems ranging from the interior of stars to the establishment of shocks at the edges of forming galaxies to neutron star mergers. Here radiative processes, absorption, emission, and scattering, are fundamental to transport of momentum and energy, as well as, the overall evolution of these and other astrophysical systems. In some cases, radiation-dominant systems are also sufficiently dense and ionized, such that, the ions are strongly coupled (ie electrostatic energy >> thermal energy or Γ >> 1), which will affect radiative processes and the transport of radiation. During this project, we developed a computational ray-tracing algorithm to study radiation transport through a simulated neutron star capsule implosion under conditions relevant to Omega-60 laser experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A two-dimensional numerical study of ion-acoustic turbulence

We investigate the linear and nonlinear evolution of the current-driven ion-acoustic instability in a collisionless plasma via two-dimensional (2-D) Vlasov–Poisson numerical simulations. We initialise the system in a stable state and gradually drive it towards instability with an imposed, weak external electric field, thus avoiding physically unrealisable super-critical initial conditions. A comprehensive analysis of the nonlinear evolution of ion-acoustic turbulence (IAT) is presented, including the detailed characteristics of the evolution of the particles’ distribution functions, (2-D) wave spectrum and the resulting anomalous resistivity. Our findings reveal the dominance of 2-D quasi-linear effects around saturation, with nonlinear effects, such as particle trapping and nonlinear frequency shifts, becoming pronounced during the later stages of the system's nonlinear evolution. Remarkably, the Kadomtsev–Petviashvili (KP) spectrum is observed immediately after the saturation of the instability. Another crucial and noteworthy result is that no steady saturated nonlinear state is ever reached: strong ion heating suppresses the instability, which implies that the anomalous resistivity associated with IAT is transient and short-lived, challenging earlier theoretical results. Towards the conclusion of the simulation, electron-acoustic waves are triggered by the formation of a double layer and strong modifications to the particle distribution induced by IAT.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Enhanced electron mixing and heating in 3-D asymmetric reconnection at the Earth's magnetopause

Here, electron heating and mixing during asymmetric reconnection are studied with a 3-D kinetic simulation that matches plasma parameters from Magnetospheric Multiscale (MMS) spacecraft observations of a magnetopause diffusion region. The mixing and heating are strongly enhanced across the magnetospheric separatrix compared to a 2-D simulation. The transport of particles across the separatrix in 3-D is attributed to lower hybrid drift turbulence excited at the steep density gradient near the magnetopause. In the 3-D simulation (and not the 2-D simulation), the electron temperature parallel to the magnetic field within the mixing layer is significantly higher than its upstream value in agreement with the MMS observations.

58 GEOSCIENCES↗
DOE OSTI2017

Correlation and transport properties for mixtures at constant pressure and temperature

Transport properties of mixtures of elements in the dense plasma regime play an important role in natural astrophysical and experimental systems, e.g., inertial confinement fusion. In this paper, we present a series of orbital-free molecular dynamics simulations on dense plasma mixtures with comparison to a global pseudo ion in jellium model. Hydrogen is mixed with elements of increasingly high atomic number (lithium, carbon, aluminum, copper, and silver) at a fixed temperature of 100 eV and constant pressure set by pure hydrogen at 2g/cm 3 , namely, 370 Mbars. We compute ionic transport coefficients, such as self-diffusion, mutual diffusion, and viscosity for various concentrations. Small concentrations of the heavy atoms significantly change the density of the plasma and decrease the transport coefficients. The structure of the mixture evidences a strong Coulomb coupling between heavy ions and the appearance of a broad correlation peak at short distances between hydrogen atoms. Finally, the concept of an effective one component plasma is used to quantify the overcorrelation of the light element induced by the admixture of a heavy element.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Study of Plasma Flow Modes in Imploding Nested Arrays

Results from experimental studies of implosion of nested wire and fiber arrays at currents of up to 4 МА at the Angara-5-1 facility are presented. Depending on the ratio between the radii of the inner and outer arrays, different modes of the plasma flow in the space between the inner and outer arrays were implemented: the sub-Alfvénic (V{sub r} < V{sub А}) and super-Alfvénic (V{sub r} > V{sub А}) modes and a mode with the formation of the transition shock wave (SW) region between the cascades. By varying the material of the outer array (tungsten wires or kapron fibers), it is shown that the plasma flow mode between the inner and outer arrays depends on the ratio between the plasma production rates ṁ{sub in} /ṁ{sub out} in the inner and outer arrays. The obtained experimental results are compared with the results of one-dimensional MHD simulation of the plasma flow between the arrays. Stable implosion of the inner array plasma was observed in experiments with combined nested arrays consisting of a fiber outer array and a tungsten inner array. The growth rates of magnetic Rayleigh−Taylor (MRT) instability in the inner array plasma at different numbers of fibers in the outer array and different ratios between the radii of the inner and outer arrays are compared. Suppression of MRT instability during the implosion of the inner array plasma results in the formation of a stable compact Z-pinch and generation of a soft X-ray pulse. A possible scenario of interaction between the plasmas of the inner and outer arrays is offered. The stability of the inner array plasma in the stage of final compression depends on the character of interaction of plasma jets from the outer array with the magnetic field of the inner array.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗