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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

A unique control system simulator for the evaluation of pulsed plasma thrusters

Because of the low thrust characteristics of solid-propellant pulsed plasma thrusters and their operational requirement to operate in a vacuum environment, unique and sensitive test techniques are required. A technique evolved for testing and evaluating pulsed plasma thrusters in an open- or closed-loop system mode employs a unique air bearing platform as a single-axis simulator on which the thruster is mounted. The simulator described was developed to evaluate pulsed plasma thrusters in the low micropound range; however, the simulator can be extended to cover the operational range of currently developed millipound thrusters.

Dahlgren, J. B.↗
DOE OSTI2025

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation

Here, this paper formulates a new particle-in-cell method for the Vlasov–Maxwell system. Under the Lorenz gauge condition, Maxwell’s equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The use of potentials for the fields motivates the adoption of a Hamiltonian formulation for particles that employs the generalized (conjugate) momentum. A notable advantage offered by the Hamiltonian formulation is the elimination of time derivatives in the Lorenz gauge formulation that are required by the standard Newton–Lorentz treatment of the particles. This allows the fields to retain the full time-accuracy guaranteed by the field solver. The resulting updates for particles require only knowledge of the fields and their spatial derivatives. An analytical method for constructing these spatial derivatives is presented that exploits the underlying integral solution used in the field solver for the wave equations. Moreover, these derivatives are demonstrated to converge at the same rate as the fields in both time and space. The Method of Lines Transpose field solver we consider in this work is globally first-order accurate in time and high-order accurate in space (e.g., fourth- and fifth-order) and belongs to a larger class of methods which are unconditionally stable, can address geometry, and leverage $\mathcal {O}(N)$ fast summation methods for efficiency. We demonstrate the method on several well-established benchmark problems on bounded domains, including a plasma sheath as well as a relativistic particle beam. The efficacy of the proposed formulation is established by comparing with a second-order accurate finite-difference time-domain method that employs a leapfrog time advance for particles and a charge conserving map suitable for bounded domains. The new method shows mesh-independent numerical heating properties even in cases where the plasma Debye length is smaller than the grid spacing. This is an important feature of the new method for problems defined on bounded domains, because it permits the use of coarser grids in space in the representation of the fields. Such a capability has significant implications for the simulation of plasmas in bounded domains with complex geometry, where the ratio between the largest and smallest cells can vary significantly. The use of high-order spatial approximations in the new method also means that fewer grid points are required in order to achieve a fixed accuracy. Our results also suggest that the new method can be used with fewer simulation particles per cell compared to the benchmark explicit method, which permits further computational savings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

3D simulation of H-mode plasmas with localized divertor impurity injection on Alcator C-Mod using the edge transport code EMC3-EIRENE

Experiments in Alcator C-Mod to assess the level of toroidal asymmetry in divertor conditions resulting from poloidally and toroidally localized extrinsic impurity gas seeding show a weak toroidal peaking (~1.1) in divertor electron temperatures for high-power Enhanced D-alpha H-mode plasmas. This is in contrast to similar experiments in Ohmically heated L-mode plasmas, which showed a clear toroidal modulation in the divertor electron temperature. Modeling of these experiments using the 3D edge transport code EMC3-EIRENE [Y. Feng, et al., J. Nucl. Mater., 241 (1997) 930] qualitatively reproduces these trends, and indicates that the different response in the simulations is due to the ionization location of the injected nitrogen. Low electron temperatures in the private flux region (PFR) in L-mode

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Laser-plasma interactions with a Fourier-Bessel particle-in-cell method

A new spectral particle-in-cell (PIC) method for plasma modeling is presented and discussed. In the proposed scheme, the Fourier-Bessel transform is used to translate the Maxwell equations to the quasi-cylindrical spectral domain. In this domain, the equations are solved analytically in time, and the spatial derivatives are approximated with high accuracy. In contrast to the finite-difference time domain (FDTD) methods, that are used commonly in PIC, the developed method does not produce numerical dispersion and does not involve grid staggering for the electric and magnetic fields. These features are especially valuable in modeling the wakefield acceleration of particles in plasmas. The proposed algorithm is implemented in the code PLARES-PIC, and the test simulations of laser plasma interactions are compared to the ones done with the quasi-cylindrical FDTD PIC code CALDER-CIRC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2001

A Coupled Plasma and Sheath Model for High Density Reactors

We present a coupled plasma and collisionless; sheath model for the simulation of high density plasma processing reactors. Due to inefficiencies in numerical schemes and the resulting computational burden, a coupled multidimensional plasma and sheath simulation has not been possible model for gas mixtures and high density reactors of practical interest. In this work we demonstrate that with a fully implicit algorithm and a refined computational mesh, a self-consistent plasma and sheath simulation is feasible. We discuss the details of the model equations, the importance of ion inertia, and the resulting sheath profiles for argon and chlorine plasmas. We find that at low operating pressures (10-30 mTorr), the charge separation occurs only within a 0.5 mm layer near the surface in a 300 mm inductively coupled plasma etch reactor. A unified model eliminates the use of off-line or loosely coupled sheath models with simplifying assumptions which generally lead to uncertainties in ion flux and sheath electrical properties.

Deepak, Bose↗
DOE OSTI2020

Plasma decay in hydrocarbons and hydrocarbon- and H 2 O-containing mixtures excited by high-voltage nanosecond discharge at elevated gas temperatures

Plasma decay after a high-voltage nanosecond discharge was experimentally and numerically studied in pure hydrocarbons (C 2 H 6 and C 3 H 8 ), and H 2 O:N 2 and C 3 H 8 :O 2 mixtures for pressures in the range 2–4 Torr and gas temperatures from 300 to 600 K. In a stoichiometric C 3 H 8 :O 2 mixture, plasma decay was also studied in a repetitively pulsed discharge for varing numbers of discharge pulses (varying degrees of fuel oxidation). The rate of plasma decay was determined from the temporal evolution of electron density measured using the microwave interferometer. It was observed that gas heating to 600 K leads to a decrease in the rate of plasma decay in all cases. The effect of heating on plasma decay was most profound in the H 2 O:N 2 mixture (after a single discharge pulse) and in the C 3 H 8 :O 2 mixture for high degrees of fuel oxidation. A kinetic scheme was developed to numerically simulate the plasma decay in hydrocarbons and combustible mixtures. Numerical analysis showed that, under the conditions studied, plasma decay was controlled by dissociative electron recombination with simple molecular and cluster ions. Gas heating led to a decrease in the rate of the electron-ion recombination and the rate of conversion of molecular ions to cluster ions. As a result, the gas temperature increase caused a decrease in the fraction of cluster ions for which the recombination coefficients are an order of magnitude higher than the recombination coefficients for molecular ions. The influence of gas heating on the decrease of the amount of cluster ions was more important when the ion composition was dominated by hydrated H 3 O + (H 2 O) k ions. The rates of the formation of these ions are extremely sensitive to any variations in gas temperature. A result, in agreement with our observations, gas heating led to an anomalous decrease in the rates of plasma decay in the H 2 O:N 2 mixture, as well as in the C 3 H 8 :O 2 mixture when H 2 O molecules were produced due to fuel oxidation.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2015

THIN CURRENT SHEETS AND ASSOCIATED ELECTRON HEATING IN TURBULENT SPACE PLASMA

Intermittent structures, such as thin current sheets, are abundant in turbulent plasmas. Numerical simulations indicate that such current sheets are important sites of energy dissipation and particle heating occurring at kinetic scales. However, direct evidence of dissipation and associated heating within current sheets is scarce. Here, we show a new statistical study of local electron heating within proton-scale current sheets by using high-resolution spacecraft data. Current sheets are detected using the Partial Variance of Increments (PVI) method which identifies regions of strong intermittency. We find that strong electron heating occurs in high PVI (>3) current sheets while no significant heating occurs in low PVI cases (<3), indicating that the former are dominant for energy dissipation. Current sheets corresponding to very high PVI (>5) show the strongest heating and most of the time are consistent with ongoing magnetic reconnection. This suggests that reconnection is important for electron heating and dissipation at kinetic scales in turbulent plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Electrochemical Corrosion Behavior of Plasma-Sprayed FeCrNiMoCBSi Amorphous/Nanocrystalline Coatings in Simulated Seawater Medium

A novel FeCrNiMoCBSi coating was fabricated on F22 marine steel by plasma spraying. The coating microstructure and composition were characterized by scanning electron microscope, transmission electron microscope, x-ray diffraction, differential scanning calorimetry, and x-ray photoelectron spectroscopy. The results showed that the coating consisted of amorphous and nanocrystalline phases. Additionally, the corrosion resistance in simulated seawater was evaluated by electrochemical corrosion methods. It was found that the coating exhibited an excellent resistance to pitting corrosion. For long-time immersion tests, the corrosion resistance was gradually decreased with the increasing through-porosity due to the permeation of the corrosive solution into the coating. The coating lost its efficacy through three stages, namely the protection stage, the permeation stage, and the failure stage.

36 MATERIALS SCIENCE↗
DOE OSTI2015

THE PLASMA ENVIRONMENT IN COMETS OVER A WIDE RANGE OF HELIOCENTRIC DISTANCES: APPLICATION TO COMET C/2006 P1 (MCNAUGHT)

On 2007 January 12, comet C/2006 P1 (McNaught) passed its perihelion at 0.17 AU. Abundant remote observations offer plenty of information on the neutral composition and neutral velocities within 1 million kilometers of the comet nucleus. In early February, the Ulysses spacecraft made an in situ measurement of the ion composition, plasma velocity, and magnetic field when passing through the distant ion tail and the ambient solar wind. The measurement by Ulysses was made when the comet was at around 0.8 AU. With the constraints provided by remote and in situ observations, we simulated the plasma environment of Comet C/2006 P1 (McNaught) using a multi-species comet MHD model over a wide range of heliocentric distances from 0.17 to 1.75 AU. The solar wind interaction of the comet at various locations is characterized and typical subsolar standoff distances of the bow shock and contact surface are presented and compared to analytic solutions. We find the variation in the bow shock standoff distances at different heliocentric distances is smaller than the contact surface. In addition, we modified the multi-species model for the case when the comet was at 0.7 AU and achieved comparable water group ion abundances, proton densities, plasma velocities, and plasma temperatures to the Ulysses/SWICS and SWOOPS observations. We discuss the dominating chemical reactions throughout the comet-solar wind interaction region and demonstrate the link between the ion composition near the comet and in the distant tail as measured by Ulysses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Final Report for Grant DE-FG02-93ER54215

The research supported by this grant led to many important publications in Physical Review Letters, Physics of Plasmas and Nuclear Fusion. Two invited talks at the American Physical Society Division of Plasma Physics, oral presentations at the Sherwood Fusion Theory Conference and MHD workshops. The most important achievements include: MHD equilibrium code FLOW for tokamak equilibria with arbitrary poloidal and toroidal flow; Two fluid tokamak equilibria with arbitrary flows; 2D Compressible MHD code to simulate tokamak plasma evolution with poloidal flow and the MHD edge pedestal; Resistive Wall Mode fluid theory in the presence of flow and differentially rotating (liquid) walls; Kinetic theory of the Resistive Wall Mode; Generalized kinetic energy principle for Resistive Wall Mode stability and the RWM kinetic stability code MISK; Kinetic theory of the m=1 mode for ITER including trapped thermal ions and electrons and fishbone oscillations driven by tangential beam injection; Two-fluid Lawson ignition criterion for tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1995

A magnetospheric magnetic field model with flexible current systems driven by independent physical parameters

A tilt-dependent magnetic field model of the Earth's magnetosphere with variable magnetopause standoff distance is presented. Flexible analytic representations for the ring and cross-tail currents, each composed of the elements derived from the Tsyganenko and Usmanov (1982) model, are combined with the fully shielded vacuum dipole configurations of Voigt (1981). Although the current sheet does not warp in the y-z plane, changes in the shape and position of the neutral sheet with dipole tilt are consistent with both MHD equilibrium theory and observations. In addition, there is good agreement with observed Delta B profiles and the average equatorial contours of magnetic field magnitude. While the dipole field is rigorously shielded within the defined magnetopause, the ring and cross-tails currents are not similarly confined, consequently, the model's region of validity is limited to the inner magnetosphere. The model depends on four independent external parameters. We present a simple but limited method of simulating several substorm related magnetic field changes associated with the disrupion of the near-Earth cross-tail current sheet and collapse of the midnight magnetotail field region. This feature further facilitates the generation of magnetic field configuration time sequences useful in plasma convection simulations of real magnetospheric events.

Hilmer, Robert V.↗
NASA NTRS2005

Particle acceleration, magnetic field generation, and emission in relativistic pair jets

Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Recent simulations show that the Weibel instability created by relativistic pair jets is responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet propagating through an ambient plasma with and without initial magnetic fields. The growth rates of the Weibel instability depends on the distribution of pair jets. The Weibel instability created in the collisionless shock accelerates particles perpendicular and parallel to the jet propagation direction. This instability is also responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron s transverse deflection behind the jet head. The jitter radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nishikawa, K.-I.↗
DOE OSTI2015

A comparative study of capacitively coupled HBr/He, HBr/Ar plasmas for etching applications: Numerical investigation by fluid model

Fluid model has been applied to perform a comparative study of hydrogen bromide (HBr)/He and HBr/Ar capacitively coupled plasma discharges that are being used for anisotropic etching process. This model has been used to identify the most dominant species in HBr based plasmas. Our simulation results show that the neutral species like H and Br, which are the key player in chemical etching, have bell shape distribution, while ions like HBr{sup +}, Br{sup +}, which play a dominant rule in the physical etching, have double humped distribution and show peaks near electrodes. It was found that the dilution of HBr by Ar and/or He results in an increase in electron density and electron temperature, which results in more ionization and dissociation and hence higher densities of neutral and charged species can be achieved. The ratio of positive ion flux to the neutral flux increases with an increase in additive gas fraction. Compare to HBr/He plasma, the HBr/Ar plasma shows a maximum change in the ion density and flux and hence the etching rate can be considered in the ion-assisted and in the ion-flux etch regime in HBr/Ar discharge. The densities of electron and other dominant species in HBr/Ar plasma are higher than those of HBr/He plasma. The densities and fluxes of the active neutrals and positive ions for etching and subsequently chemical etching versus physical sputtering in HBr/Ar and HBr/He plasmas discharge can be controlled by tuning gas mixture ratio and the desire etching can be achieved.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Validation of NSFsim as a Grad-Shafranov equilibrium solver at DIII-D

Plasma shape is a significant factor that must be considered for any Fusion Pilot Plant (FPP) as it has significant consequences for plasma stability and core confinement. A new simulator, NSFsim, has been developed based on a historically successful code, DINA [1], offering tools to simulate both transport and plasma shape. Specifically, NSFsim is a free boundary equilibrium and transport solver and has been configured to match the properties of the DIII-D tokamak. This paper is focused on validating the Grad-Shafranov (GS) solver of NSFsim by analyzing its ability to recreate the plasma shape, the poloidal flux distribution, and the measurements of the simulated diagnostic signals originating from flux loops and magnetic probes in DIII-D. Five different plasma shapes are simulated to show the robustness of NSFsim to different plasma conditions; these shapes are Lower Single Null (LSN), Upper Single Null (USN), Double Null (DN), Inner Wall Limited (IWL), and Negative Triangularity (NT). The NSFsim results are compared against real measured signals, magnetic profile fits from EFIT [2], and another plasma equilibrium simulator, GSevolve [3]. EFIT reconstructions of shots are readily available at DIII-D, but GSevolve was manually ran by us to provide simulation data to compare against.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Kinetic Physics of Homogeneous Turbulence in Collisionless Plasmas

We present results of large scale, fully kinetic 3D simulations of decaying plasma turbulence that span a range of scales from tens of ion gyroradii to the electron Debye length. A significant fraction of the dissipated energy is transferred to ions, in contrast to the scenario where the dominant dissipation is due to Landau damping of wave-like fluctuations. Intense current sheets are formed in the simulation with thicknesses over a range of scales approximately between electron and ion inertial scales. The rate of ion heating is well correlated with intermittency level. Geometrical analysis of the current sheets suggest that magnetic reconnection is active in the present simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

BOUT++ nonlinear simulation for a comparative study with the measured 2D ELM structures in the KSTAR H-mode plasma

Here, time-dependent nonlinear simulation using the BOUT++ 3-field model is performed for a study of edge localized mode (ELM) dynamics in KSTAR H-mode plasmas. The simulated results of the ELM evolution are directly compared with the 2D ELM structures, measured by electron cyclotron emission imaging systems, from the initial growth to crash phase. The simulation results are qualitatively in good agreement with the observations. Specifically, it is notable in simulation and experiment that higher-n modes in a quasi-steady state phase evolve into lower-n modes close to the ELM-crash in time. It indicates that lower-n modes may be responsible for the onset of the crash. A postanalysis confirms that nonlinear mode coupling between toroidal harmonics may drive the lower-n mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗