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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Simulations and observations of plasma depletion, ion composition, and airglow emissions in two auroral ionospheric depletion experiments

Observations of plasma depletion, ion composition modification, and airglow emissions in the Waterhole experiments are presented. The detailed ion chemistry and airglow emission processes related to the ionospheric hole formation in the experiment are examined, and observations are compared with computer simulation results. The latter indicate that the overall depletion rates in different parts of the depletion region are governed by different parameters.

Yau, A. W.↗
DOE OSTI2022

Verification of a fully kinetic ion model for electromagnetic simulations of high-frequency waves in toroidal geometry

For the research of high-frequency electromagnetic waves in tokamaks, an electromagnetic simulation model, in which the ion dynamics is described by a six-dimensional Vlasov equation and the electron dynamics is described by a drift kinetic equation, is formulated and implemented in the global gyrokinetic toroidal code (GTC). Analytic dispersion relations are derived in reduced systems and compared with various theories to verify the model. Linear simulations of a generalized ion Bernstein wave and ion cyclotron emission are verified by comparing the GTC simulation results with analytic dispersion relation theory and magnetoacoustic cyclotron instability theory, respectively, in cylindrical geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2003

Numerical Simulation of Waves Driven by Plasma Currents Generated by Low-Frequency Alfven Waves in a Multi-Ion Plasma

When multi-ion plasma consisting of heavy and light ions is permeated by a lowfrequency Alfien (LFA) wave, the EXB and the polarization drifts of the different ion species and the electrons could be quite different. The relative drifts between the charged-particle species drive waves, which energize the plasma. Using 2.5-D particle-in-cell simulations, we study this process of wave generation and its nonlinear consequences in terms of acceleration and heating plasma. Specifically we study the situation for LFA wave frequency being lower than the heavyion cyclotron frequency in a multi-ion plasma. We impose such a wave to the plasma assuming that its wavelength is much larger than that of the waves generated by the relative drifts. For better understanding, the LFA-wave driven simulations are augmented by those driven by initialized ion beams.

Singh, Nagendra↗
DOE OSTI2017

HIS-TORIC: extending core ICRF wave simulation to include realistic SOL plasmas

Here, this paper presents a novel approach to incorporating an arbitrarily shaped edge scrape-off-layer (SOL) plasma and an ion cyclotron range of frequency (ICRF) antenna structure into existing core ICRF wave simulation models. We partition the entire computation domain into two sub-domains: a core and an edge region. Simulations in each domain are performed separately with appropriate numerical solvers. For the core, the TORIC ICRF solver (Brambilla 1999 Plasma Phys. Control. Fusion 41 1) was modified to impose an essential (Dirichlet) boundary condition at its interface with the edge domain. In the edge, a finite element method is used to solve a cold collisional plasma model. The domains are then joined together using the continuity boundary condition for the tangential electric and magnetic fields at their interfaces (Hybrid Integration of SOL to TORIC: HIS-TORIC). The model developed here was tested using an ICRH H minority heating scenario on the Alcator C-Mod tokamak (Hutchinson et al 1994 Phys. Plasmas 1 1511). The simulated pattern of core wave propagation agrees well with a standard TORIC simulation. This approach opens the possibility of using a realistic diverted SOL plasma and a complicated 3D RF antenna together with a rigorous hot core plasma model, while requiring only minimal modification to existing RF codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Numerical simulation of an atmospheric pressure plasma jet with coaxial shielding gas

A two-dimensional ( r, z ) numerical simulation of the discharge characteristics of an atmospheric pressure plasma jet (APPJ), with coaxial shielding gas, was performed. The helium working gas flowed in a central capillary tube, engulfing a needle electrode powered by 13.7 MHz radio frequency sinusoidal voltage. The N 2 shielding gas flowed in the annular space of a coaxial tube. These gases emerged, in laminar flow, in a 78%N 2 -21%O 2 -1%Ar dry air ambient. The characteristics of the APPJ with shielding gas were compared to those of the APPJ without shielding gas. The nitrogen shielding gas hindered the diffusion of oxygen and argon from the ambient air into the helium jet. With the shielding gas present, more nitrogen penetrated into the helium core, causing a shorter plasma ‘plume’. The flow rates of the working and shielding gas, critically affected the gas temperature, and in turn the discharge characteristics. For a He flow of 2 standard liters per minute (slm), switching on the nitrogen shielding gas flow (at 4.5 slm) reduced the on-axis O 2 and Ar mole fractions from 3.9 × 10 - 4 to 6.8 × 10 - 5 and from 1.9 × 10 - 5 to 3.3 × 10 - 6 , respectively, at an axial distance of 3 mm downstream of the nozzle. The radial profiles of the mole fractions of the ambient gases were monotonically and strongly decreasing towards the system axis, for short axial distances from the nozzle (~1 mm), but became progressively flatter at longer distances from the nozzle (3 mm and 5 mm). Simulation predictions captured the salient features of experimental data of ambient species mole fractions in the plasma jet, and the 706 nm optical emission intensity profiles of the He 3 3 S excited state.

Physics↗
DOE OSTI2015

The effects of dilution on turbulence and transport in C-Mod ohmic plasmas and comparisons with gyrokinetic simulations

Main ion dilution has been predicted by gyrokinetic simulations to have a significant effect on ion thermal transport in C-Mod ohmic plasmas. This effect was verified experimentally with a specific set of experiments on C-Mod in which ohmic deuterium plasmas across the linear ohmic confinement (LOC) through the saturated ohmic confinement (SOC) regimes were diluted by seeding with nitrogen gas (Z=7) injection. The seeding was observed to increase the normalized ion temperature gradients (ITGs) by up to 30% without a corresponding increase in the gyrobohm normalized ion energy flux, indicating a change in either the stiffness or the critical ion temperature gradient associated with ITG turbulence. The seeding also reversed the direction of the intrinsic toroidal rotation in plasmas slightly above the normal intrinsic rotation reversal critical density. GYRO simulations of the seeded and unseeded plasmas show that the seeding affected both the critical gradient and the stiffness. For plasmas in the LOC regime, the dilution primarily increased the critical gradient, while for plasmas in the SOC regime the dilution primarily decreased the stiffness. At r/a = 0.8, where the experimental fluxes were above marginal stability, local GYRO predicted and experimental energy fluxes agreed, except for Q i in the SOC regime where GYRO under-predicted the experimental energy flux. At r/a = 0.6, where the experimental fluxes were close to marginally stable, local GYRO predicted ITG modes to be strongly unstable and are responsible for both Q i and Q e (with Q i > Q e ), as opposed to the experiment where Q i < Q e . In contrast, global GYRO in this region predicted the ITG modes to be closer to marginal stability, and accurately predict the experimental Q i when the T i profile is modified within experimental uncertainties. Further, the fact that Q e is always less than Q i in the r/a = 0.6 simulations with k θ ρ s ≤1 indicates that high-k electron temperature gradient driven (ETG) modes must be included in future simulations and may be responsible for the electron energy transport in this case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry

A brute-force, long-time gyrokinetic simulation of plasma profile evolution in magnetic fusion devices is not desirable due to large computational resource requirements and a possible accumulation of numerical error. The equation-free projective integration method of Keverekidis et al. [Commun. Math. Sci. 1(4), 715–762 (2003)] is one of the outstanding candidates in projecting micro-scale simulations to a longer timescale. However, its application to tokamak plasma has not been fruitful due to the appearance of spurious transient oscillations in the lifting process, which are present when the kinetic simulations are initialized with a simplified model distribution function and which make the kinetic simulations to deviate from the desired paths. In this work, a kinetically informed lifting algorithm is added to the equation-free projective integration method, which is then verified in the electrostatic gyrokinetic particle-in-cell code XGCa [R. Hager and C. S. Chang, Phys. Plasmas 23, 042503 (2016)] for a neoclassical ion heat transport problem with adiabatic electrons. This new lifting operator is demonstrated to control spurious transients, enabling an over four-times reduction in the overall computing time in the time-evolution of the ion temperature profile in an axisymmetric toroidal plasma. Further reduction in the computing time is found to be limited due to the stability properties of the linear least squares projective integrator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Accuracy of the explicit energy-conserving particle-in-cell method for under-resolved simulations of capacitively coupled plasma discharges

The traditional explicit electrostatic momentum-conserving particle-in-cell algorithm requires strict resolution of the electron Debye length to deliver numerical stability and accuracy. The explicit electrostatic energy-conserving particle-in-cell algorithm alleviates this constraint with minimal modification to the traditional algorithm, retaining its simplicity, ease of parallelization, and acceleration on modern supercomputing architectures. In this article, we apply the algorithm to model a one-dimensional radio frequency capacitively coupled plasma discharge relevant to industrial applications. The energy-conserving approach closely matches the results from the momentum-conserving algorithm and retains accuracy even for cell sizes up to 8 times the electron Debye length. For even larger cells, the algorithm loses accuracy due to poor resolution of steep gradients within the radio frequency sheath. Accuracy can be recovered by adopting a non-uniform grid, which resolves the sheath and allows for cell sizes up to 32 times the electron Debye length in the quasi-neutral bulk of the discharge. The effect is an up to 8 times reduction in the number of required simulation cells, an improvement that can compound in higher-dimensional simulations. We therefore consider the explicit energy-conserving algorithm as a promising approach to significantly reduce the computational cost of full-scale device simulations and a pathway to delivering kinetic simulation capabilities of use to industry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2012

Observations and Simulations of Formation of Broad Plasma Depletions Through Merging Process

Broad plasma depletions in the equatorial ionosphere near dawn are region in which the plasma density is reduced by 1-3 orders of magnitude over thousands of kilometers in longitude. This phenomenon is observed repeatedly by the Communication/Navigation Outage Forecasting System (C/NOFS) satellite during deep solar minimum. The plasma flow inside the depletion region can be strongly upward. The possible causal mechanism for the formation of broad plasma depletions is that the broad depletions result from merging of multiple equatorial plasma bubbles. The purpose of this study is to demonstrate the feasibility of the merging mechanism with new observations and simulations. We present C/NOFS observations for two cases. A series of plasma bubbles is first detected by C/NOFS over a longitudinal range of 3300-3800 km around midnight. Each of the individual bubbles has a typical width of approx 100 km in longitude, and the upward ion drift velocity inside the bubbles is 200-400 m/s. The plasma bubbles rotate with the Earth to the dawn sector and become broad plasma depletions. The observations clearly show the evolution from multiple plasma bubbles to broad depletions. Large upward plasma flow occurs inside the depletion region over 3800 km in longitude and exists for approx 5 h. We also present the numerical simulations of bubble merging with the physics-based low-latitude ionospheric model. It is found that two separate plasma bubbles join together and form a single, wider bubble. The simulations show that the merging process of plasma bubbles can indeed occur in incompressible ionospheric plasma. The simulation results support the merging mechanism for the formation of broad plasma depletions.

plasma bubbles↗
DOE OSTI2015

Trapping of intense light in hollow shell

A small hollow shell for trapping laser light is proposed. Two-dimensional particle-in-cell simulation shows that under appropriate laser and plasma conditions a part of the radiation fields of an intense short laser pulse can enter the cavity of a small shell through an over-critical density plasma in an adjacent guide channel and become trapped. The trapped light evolves into a circulating radial wave pattern until its energy is dissipated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Verification of gyrokinetic particle simulations of neoclassical tearing modes in fusion plasmas

Abstract The ability to simulate neoclassical tearing modes (NTMs) in the gyrokinetic toroidal code (GTC) has been developed and verified, in which ions are treated with a gyrokinetic model and electrons are treated as a resistive fluid. The simulation results demonstrate that the neoclassical bootstrap current effect can destabilize an otherwise stable classical tearing mode. In the cylindrical geometry, GTC simulations in the magnetohydrodynamic limit show quantitative agreement with the modified Rutherford theory, both in terms of the scaling law in the small island limit and in terms of the saturation level and pressure flattening effect in the large island limit. The toroidal effects are slightly destabilizing for the NTM, while the kinetic effects of thermal ions are stabilizing for the NTM and increase its excitation threshold.

Physics↗
DOE OSTI2020

Interaction of ultra relativistic ${e}^{-}{e}^{+}$ fireball beam with plasma

Ab initio simulations of the propagation in a plasma of a soon to be available relativistic electron–positron beam or fireball beam provide an effective mean for the study of microphysics relevant to astrophysical scenarios. We show that the current filamentation instability associated with some of these scenarios reaches saturation after only 10 cm of propagation in a typical laboratory plasma with a density~1017 cm-3 . The different regimes of the instability, from the purely transverse to the mixed mode filamentation, can be accessed by varying the background plasma density. The instability generates large local plasma gradients, intense transverse magnetic fields, and enhanced emission of radiation. We suggest that these effects may be observed experimentally for the first time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Developing the science and technology for the Material Plasma Exposure eXperiment

Linear plasma generators are cost effective facilities to simulate divertor plasma conditions of present and future fusion reactors. They are used to address important R&D gaps in the science of plasma material interactions and towards viable plasma facing components for fusion reactors. Next generation plasma generators have to be able to access the plasma conditions expected on the divertor targets in ITER and future devices. The steady-state linear plasma device MPEX will address this regime with electron temperatures of 1–10 eV and electron densities of $10^{21}{\text{}}\!-\!10^{20}$ ${\rm m}^{-3}$. The resulting heat fluxes are about 10 MW ${\rm m}^{-2}$ . MPEX is designed to deliver those plasma conditions with a novel Radio Frequency plasma source able to produce high density plasmas and heat electron and ions separately with electron Bernstein wave (EBW) heating and ion cyclotron resonance heating with a total installed power of 800 kW. The linear device Proto-MPEX, forerunner of MPEX consisting of 12 water-cooled copper coils, has been operational since May 2014. Its helicon antenna (100 kW, 13.56 MHz) and EC heating systems (200 kW, 28 GHz) have been commissioned and 14 MW ${\rm m}^{-2}$ was delivered on target. Furthermore, electron temperatures of about 20 eV have been achieved in combined helicon and ECH heating schemes at low electron densities. Overdense heating with EBW was achieved at low heating powers. The operational space of the density production by the helicon antenna was pushed up to $1.1 \times 10^{20}$ ${\rm m}^{-3}$ at high magnetic fields of 1.0 T at the target. Finally, the experimental results from Proto-MPEX will be used for code validation to enable predictions of the source and heating performance for MPEX. MPEX, in its last phase, will be capable to expose neutron-irradiated samples. In this concept, targets will be irradiated in ORNL's High Flux Isotope Reactor and then subsequently exposed to fusion reactor relevant plasmas in MPEX.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1987

Real time animation of space plasma phenomena

In pursuit of real time animation of computer simulated space plasma phenomena, the code was rewritten for the Massively Parallel Processor (MPP). The program creates a dynamic representation of the global bowshock which is based on actual spacecraft data and designed for three dimensional graphic output. This output consists of time slice sequences which make up the frames of the animation. With the MPP, 16384, 512 or 4 frames can be calculated simultaneously depending upon which characteristic is being computed. The run time was greatly reduced which promotes the rapid sequence of images and makes real time animation a foreseeable goal. The addition of more complex phenomenology in the constructed computer images is now possible and work proceeds to generate these images.

Jordan, K. F.↗
DOE OSTI2015

WarpVisit

WarpVisit is an insitu simulation application that integrates the Warp laser plasma accelerator simulation framework with Visit a parallel visualization application. WarpVisit is written in python and supports interactive or live mode where user can connect to Warp with the Visit GUI and batch mode for batch for non-interactive use on high-performance computing resources.

Loring, Burlen↗
DOE OSTI2017

Quasi-optical simulation of the electron cyclotron plasma heating in a mirror magnetic trap

The resonance microwave plasma heating in a large-scale open magnetic trap is simulated taking into account all the basic wave effects during the propagation of short-wavelength wave beams (diffraction, dispersion, and aberration) within the framework of the consistent quasi-optical approximation of Maxwell’s equations. The quasi-optical method is generalized to the case of inhomogeneous media with absorption dispersion, a new form of the quasi-optical equation is obtained, the efficient method for numerical integration is found, and simulation results are verified on the GDT facility (Novosibirsk).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2005

NASA GRC and MSFC Space-Plasma Arc Testing Procedures

Tests of arcing and current collection in simulated space plasma conditions have been performed at the NASA Glenn Research Center (GRC) in Cleveland, Ohio, for over 30 years and at the Marshall Space flight Center (MSFC) for almost as long. During this period, proper test conditions for accurate and meaningful space simulation have been worked out, comparisons with actual space performance in spaceflight tests and with real operational satellites have been made, and NASA has achieved our own internal standards for test protocols. It is the purpose of this paper to communicate the test conditions, test procedures, and types of analysis used at NASA GRC and MSFC to the space environmental testing community at large, to help with international space-plasma arcing testing standardization. To be discussed are: 1. Neutral pressures, neutral gases, and vacuum chamber sizes. 2. Electron and ion densities, plasma uniformity, sample sizes, and Debye lengths. 3. Biasing samples versus self-generated voltages. Floating samples versus grounded. 4. Power supplies and current limits. Isolation of samples from power supplies during arcs. Arc circuits. Capacitance during biased arc-threshold tests. Capacitance during sustained arcing and damage tests. Arc detection. Preventing sustained discharges during testing. 5. Real array or structure samples versus idealized samples. 6. Validity of LEO tests for GEO samples. 7. Extracting arc threshold information from arc rate versus voltage tests. 8 . Snapover and current collection at positive sample bias. Glows at positive bias. Kapton pyrolization. 9. Trigger arc thresholds. Sustained arc thresholds. Paschen discharge during sustained arcing. 10. Testing for Paschen discharge thresholds. Testing for dielectric breakdown thresholds. Testing for tether arcing. 11. Testing in very dense plasmas (ie thruster plumes). 12. Arc mitigation strategies. Charging mitigation strategies. Models. 13. Analysis of test results. Finally, the necessity of testing will be emphasized, not to the exclusion of modeling, but as part of a complete strategy for determining when and if arcs will occur, and preventing them from occurring in space.

Ferguson, Dale C.a↗