Search NASASearch

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 73 records · Page 4

Numerical simulation of plasma transport driven by the Io torus

The Rice convection model (RCM) has been modified to a form suitable for Jupiter (RCM-J) to study plasma interchange motion in and near the Io plasma torus. The net result of the interchange is that flux tubes, heavily loaded with torus plasma, are transported outward, to be replaced by tubes containing little low-energy (less than 1 keV) plasma. The process is numerically simulated in terms of time evolution from an initial torus that is longitudinally asymmetric and with gradually decreasing density outward from Io's orbit. In the simulations, the nonlinear stage of the instability characteristically exhibits outreaching fingers of heavily-loaded flux tubes that lengthen at an accelerating rate. The principal finding is that the primary geometrical form of outward transport of torus plasma in Jupiter's magnetosphere is through long, outward-moving fingers of plasma. In the simulations, the fingers mainly form in the active sector of the Io torus (the heavier side of the asymmetric torus), and they are spaced longitudinally roughly 20 deg apart.

Yang, Y. S.

Radio Plasma Imager Simulations and Measurements

The Radio Plasma Imager (RPI) will be the first-of-its kind instrument designed to use radio wave sounding techniques to perform repetitive remote sensing measurements of electron number density (N(sub e)) structures and the dynamics of the magnetosphere and plasmasphere. RPI will fly on the Imager for Magnetopause-to-Aurora Global Exploration (IMAGE) mission to be launched early in the year 2000. The design of the RPI is based on recent advances in radio transmitter and receiver design and modern digital processing techniques perfected for ground-based ionospheric sounding over the last two decades. Free-space electromagnetic waves transmitted by the RPI located in the low density magnetospheric cavity will be reflected at distant plasma cutoffs. The location and characteristics of the plasma at those remote reflection points can then be derived from measurements of the echo amplitude, phase, delay time, frequency, polarization, Doppler shift, and echo direction. The 500 m tip-to-tip X and Y (spin plane) antennas and 20 m boom Z axis antenna on RPI will be used to measures echoes coming from distances of several R(sub E).

Green, J. L.

A simulation model of time-dependent plasma-spacecraft interactions

A plasma simulation code is presented that models the time-dependent plasma properties in the vicinity of a spherical, charged spacecraft. After showing agreement with analytic, steady-state theories and ATS-6 satellite data, the following three problems are treated: (1) transient pulses from photoemission at various emission temperatures and ambient plasma conditions, (2) spacecharge limited emission, and (3) simulated plasma oscillations in the long wavelength limit.

Rothwell, P. L.

High-Density Plasma Reactors: Simulations for Design

The development of improved and more efficient plasma reactors is a costly process for the semiconductor industry. Until five years ago, the Industry made most of its advancements through a trial and error approach. More recently, the role of computational modeling in the design process has increased. Both conventional computational fluid dynamics (CFD) techniques like Navier-Stokes solvers as well as particle simulation methods are used to model plasma reactor flowfields. However, since high-density plasma reactors generally operate at low gas pressures on the order of 1 to 10 mTorr, a particle simulation may be necessary because of the failure of CFD techniques to model rarefaction effects. The direct simulation Monte Carlo method is the most widely accepted and employed particle simulation tool and has previously been used to investigate plasma reactor flowfields. A plasma DSMC code is currently under development at NASA Ames Research Center with its foundation as the object-oriented parallel Cornell DSMC code, MONACO. The present investigation is a follow up of a neutral flow investigation of the effects of process parameters as well as reactor design on etch rate and etch rate uniformity. The previous work concentrated on silicon etch of a chlorine flow in a configuration typical of electron cyclotron resonance (ECR) or helical resonator type reactors. The effects of the plasma on the dissociation chemistry were modeled by making assumptions about the electron temperature and number density. The electrons or ions themselves were not simulated.The present work extends these results by simulating the charged species.The electromagnetic fields are calculated such that power deposition is modeled self-consistently. Electron impact reactions are modeled along with mechanisms for charge exchange. An bipolar diffusion assumption is made whereby electrons remain tied to the ions. However, the velocities of tile electrons are allowed to be modified during collisions and are not confined to a Maxwellian distribution. The interaction between the neutral flow and plasma is examined, and results for etch rate uniformity from the previous research and the present plasma simulations are compared.

Hash, David B.

Simulation of plasma double-layer structures

Electrostatic plasma double layers are numerically simulated by means of a magnetized 2 1/2 dimensional particle in cell method. The investigation of planar double layers indicates that these one dimensional potential structures are susceptible to periodic disruption by instabilities in the low potential plasmas. Only a slight increase in the double layer thickness with an increase in its obliqueness to the magnetic field is observed. Weak magnetization results in the double layer electric field alignment of accelerated particles and strong magnetization results in their magnetic field alignment. The numerical simulations of spatially periodic two dimensional double layers also exhibit cyclical instability. A morphological invariance in two dimensional double layers with respect to the degree of magnetization implies that the potential structures scale with Debye lengths rather than with gyroradii. Electron beam excited electrostatic electron cyclotron waves and (ion beam driven) solitary waves are present in the plasmas adjacent to the double layers.

Borovsky, J. E.

Plasma Heating Simulation in the VASIMR System

The paper describes the recent development in the simulation of the ion-cyclotron acceleration of the plasma in the VASIMR experiment. The modeling is done using an improved EMIR code for RF field calculation together with particle trajectory code for plasma transport calculat ion. The simulation results correlate with experimental data on the p lasma loading and predict higher ICRH performance for a higher density plasma target. These simulations assist in optimizing the ICRF anten na so as to achieve higher VASIMR efficiency.

Ilin, Andrew V.

Numerical Simulations of Plasma Based Flow Control Applications

A mathematical model was developed to simulate flow control applications using plasma actuators. The effects of the plasma actuators on the external flow are incorporated into Navier Stokes computations as a body force vector. In order to compute this body force vector, the model solves two additional equations: one for the electric field due to the applied AC voltage at the electrodes and the other for the charge density representing the ionized air. The model is calibrated against an experiment having plasma-driven flow in a quiescent environment and is then applied to simulate a low pressure turbine flow with large flow separation. The effects of the plasma actuator on control of flow separation are demonstrated numerically.

Suzen, Y. B.

Plasma Jet Simulations Using a Generalized Ohm's Law

Plasma jets are important physical phenomena in astrophysics and plasma propulsion devices. A currently proposed dual jet plasma propulsion device to be used for ISS experiments strongly resembles a coronal loop and further draws a parallel between these physical systems [1]. To study plasma jets we use numerical methods that solve the compressible MHD equations using the generalized Ohm s law [2]. Here, we will discuss the crucial underlying physics of these systems along with the numerical procedures we utilize to study them. Recent results from our numerical experiments will be presented and discussed.

Ebersohn, Frans

Investigation of Drift Effects in UEDGE Simulations of NSTX-U Edge Plasma With Lithium Divertors

Lithium is a low-Z material, and lithium-based plasma-facing components (PFCs) are planned for the National Spherical Torus Experiment Upgrade (NSTX-U) to explore potential benefits for divertor power exhaust and core plasma management. NSTX-U is a medium-sized spherical tokamak with up to 12 MW of auxiliary heating, capable of generating reactor-relevant plasma conditions. This work presents boundary plasma simulations for NSTX-U with lithium PFCs using the UEDGE code, incorporating full magnetic and 𝐄 ×𝐁 drift physics. The simulations show that drifts strongly influence heat and particle transport: they enhance convective transport, broaden the scrape-off layer heat-flux width 𝜆 𝑞 , and reduce the anomalous heat diffusivity 𝜒 required to reproduce predicted SOL heat-flux width. 𝐄 ×𝐁 drifts provide poloidal transport, while ∇𝐵 (which includes both gradB and curvature) drifts provide radial heat and particle transport. Lithium transport is also affected by drifts, with lithium ions migrating from the outer divertor to the inner divertor through the private flux region (PFR) following the 𝐄 ×𝐁 drifts flow, lowering upstream impurity lithium densities. UEDGE is self-consistently coupled with the Wall-Li model to study plasma lithium PFC interactions depending on the local lithium sourcing based on local plasma conditions and lithium surface temperature. In these simulations, lithium evaporation shows a vapor-shielding effect that reduces divertor heat flux and increases radiative losses once surface temperatures exceed 450°C. This research work provides a first step toward self-consistent modeling of lithium PFCs in NSTX-U, demonstrating the impact of drift-driven plasma transport in SOL and divertor regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Current-driven double layers and the auroral plasma

Simulation results on current-driven double layers are presented. The role of the Buneman instability in the formation of double layers is clearly identified. The recurring formation and motion of double layers produce a temporal variation of the electric fields similar to that observed in the auroral plasma by the S3-3 satellite.

Singh, N.

Direct measurements of plasma characteristics in space-simulation beam-plasma interactions

Recent laboratory investigations of space-simulated electron-beam-plasma interactions are discussed. The plasma has been characterized with respect to its spatial distributions in density, temperature, and density fluctuation power spectra. The latter results have been further classified in terms of amplitude and spectral distributions. The overall results include: (1) detailed radial profiles of plasma density as a function of beam parameters; (2) the association of low-frequency large amplitude fluctuations with ion acoustic, ion cyclotron, and drift-wave modes; (3) the identification and spatial mapping of suprathermal electrons in the beam core; and (4) the experimental definition of a density-dependent criterion for the ignition of the beam-plasma-discharge.

Szuszczewicz, E. P.

Simulations of plasma current induced by toroidal field ramping down on tokamaks

In tokamak experiments, the BT ramping down can drive off-axis and parallel inductive current. This approach leads to a decrease in q95 and an increase in normalized beta, βN. The off-axis inductive current also broadens the current profile. Typically, the 1.5D transport code is used to simulate the time evolution of the plasma current profile, which is based on the flux-surface-averaged Faraday’s law. The ONETWO code is one of those transport codes. However, this code cannot simulate the situation of the evolving toroidal field BT. This article proposes a new time-dependent model to take into account the BT ramping down situation. A modified formula of flux-surface-averaged Faraday’s law was derived to consider the effect of the BT ramping rate on the current evolution, and it was implemented in the ONETWO code. Then, the modified ONETWO code was used to simulate the BT ramping down experiment on DIII-D. The simulation result of the plasma current evolution with BT ramping down shows a broader current profile with smaller ohmic current induced by the poloidal field, compared to that without BT ramping down.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Scale Separation Effects on Simulations of Plasma Turbulence

Abstract Understanding plasma turbulence requires a synthesis of experiments, observations, theory, and simulations. In the case of kinetic plasmas such as the solar wind, the lack of collisions renders the fluid closures such as viscosity meaningless and one needs to resort to higher-order fluid models or kinetic models. Typically, the computational expense in such models is managed by simulating artificial values of certain parameters such as the ratio of the Alfvén speed to the speed of light ( v A / c ) or the relative mass ratio of ions and electrons ( m i / m e ). Although, typically care is taken to use values as close as possible to realistic values within the computational constraints, these artificial values could potentially introduce unphysical effects. These unphysical effects could be significant at sub-ion scales, where kinetic effects are the most important. In this paper, we use the 10-moment fluid model in the Gkeyll framework to perform controlled numerical experiments, systematically varying the ion–electron mass ratio from a small value down to the realistic proton–electron mass ratio. We show that the unphysical mass ratio has a significant effect on the kinetic range dynamics as well as the heating of both plasma species. The dissipative process for both ions and electrons becomes more compressive in nature, although the ions remain nearly incompressible in all cases. The electrons move from being dominated by incompressive viscous-like heating/dissipation to very compressive heating/dissipation dominated by compressions/rarefactions. While the heating change is significant for the electrons, a mass ratio of m i / m e ∼ 250 captures the asymptotic behavior of electron heating.

79 ASTRONOMY AND ASTROPHYSICS

Numerical simulation of plasma insulator interactions in space. Part 1: The self consistent calculation

A computer program is being developed to simulate the interaction of a plasma with a conducting disk partially covered by an insulator. Initial runs consider only charge sticking to the dielectric. Results indicate that the current density drawn by the hole in the dielectric increases approximately linearly with voltage for conductor voltages between 5 volts and 250 volts.

Nonnast, J. H.

Simulation of plasma and neutral transport in PISCES-RF using SOLPS-ITER

In this research, the fluid plasma transport code SOLPS-ITER is applied and validated against experimental data from the plasma interaction surface component experimental station (PISCES)-RF linear plasma device to establish a physics basis for plasma and neutral transport in its two magnetic field (B-field) geometry setups-(1)the cusp and (2) non-cusp or linear B-field. The main focus of this study is to understand (1) radial plasma transport (2) heat and particle loads on the upstream dump and downstream target plate, and (3) the physics of plasma-neutral interactions in PISCES-RF. The simulation setup adheres to typical PISCES-RF experimental conditions, with a 2D helicon power deposition profile as an input heating source. SOLPS-ITER simulations reproduce experimental conditions with the Bohm diffusion model for both B-field configurations of the PISCES-RF experiment. Major energy loss channels include neutral radiation and power deposited on the wall and dump plate, with only 1% of the input power reaching the target. The ionization front is well confined near the dump plate due to the heating and puffing regions. Additionally, SOLPS-ITER simulation results are also found to be in very good agreement with the particle-in-cell calculations using the code-PICOS++ which supports the validity of SOLPS in low collisionality regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY