Analysis of electron analog simulation of ion containment in magnetic fields
Physical model for plasma containment simulation in magnetic fields
Kaufman, H. R.↗
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.
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.↗
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↗
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 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↗
Field-aligned plasma flow in MHD simulations of magnetotail reconnection and the formation of boundary layers
The outflow from a reconnection region in a realistic magnetotail geometry is studied using a two-dimensional time-dependent, compressible, resistive MHD code. Two cases are emphasized: (1) the evolution of near-earth reconnection, which grows in the form of an internal unstable mode after gradual externally forced changes have initiated an anomalous dissipation process, and (2) the evolution of more distant reconnection under influence of a nonuniform inflow that forces reconnection to occur at a given location in the distant tail. In both cases, it is demonstrated that plasma flow is primarily parallel to the magnetic field in regions away from the localized area of reconnection and outside a narrow central layer.
Birn, J.↗
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) in Huntsville, Alabama, 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 Debuy 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. 5. Arc circuits. Capacitance during biased arc-threshold tests. Capacitance during sustained arcing and damage tests. Arc detection. Prevention sustained discharges during testing. 6. Real array or structure samples versus idealized samples. 7. Validity of LEO tests for GEO samples. 8. Extracting arc threshold information from arc rate versus voltage tests. 9. Snapover and current collection at positive sample bias. Glows at positive bias. Kapon (R) pyrolisis. 10. Trigger arc thresholds. Sustained arc thresholds. Paschen discharge during sustained arcing. 11. Testing for Paschen discharge threshold. Testing for dielectric breakdown thresholds. Testing for tether arcing. 12. Testing in very dense plasmas (ie thruster plumes). 13. Arc mitigation strategies. Charging mitigation strategies. Models. 14. 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.↗
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) in Huntsville, Alabama, 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. Discussed herein are neutral gas conditions, plasma densities and uniformity, vacuum chamber sizes, sample sizes and Debye lengths, biasing samples versus self-generated voltages, floating samples versus grounded samples, test electrical conditions, arc detection, preventing sustained discharges during testing, real samples versus idealized samples, validity of LEO tests for GEO samples, extracting arc threshold information from arc rate versus voltage tests, snapover, current collection, and glows at positive sample bias, Kapton pyrolysis, thresholds for trigger arcs, sustained arcs, dielectric breakdown and Paschen discharge, tether arcing and testing in very dense plasmas (i.e. thruster plumes), arc mitigation strategies, charging mitigation strategies, models, and 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.↗
When should PIC simulations be applied to atmospheric pressure plasmas? Impact of correlation heating
Abstract Molecular dynamics simulations are used to test when the particle-in-cell (PIC) method applies to atmospheric pressure plasmas. It is found that PIC applies only when the plasma density and macroparticle weight are sufficiently small because of two effects associated with correlation heating. The first is the physical effect of disorder-induced heating (DIH). This occurs if the plasma density is large enough that a species (typically ions) is strongly correlated in the sense that the Coulomb coupling parameter exceeds one. In this situation, DIH causes ions to rapidly heat following ionization. PIC is not well suited to capture DIH because doing so requires using a macroparticle weight of one and a grid that well resolves the physical interparticle spacing. These criteria render PIC intractable for macroscale domains. The second effect is a numerical error due to artificial correlation heating (ACH). ACH is like DIH in that it is caused by the Coulomb repulsion between particles, but differs in that it is a numerical effect caused by a macroparticle weight larger than one. Like DIH, it is associated with strong correlations. However, here the macroparticle coupling strength is found to scale as Γ w 2 / 3 , where Γ is the physical coupling strength and w is the macroparticle weight. So even if the physical coupling strength of a species is small, as is expected for electrons in atmospheric pressure plasmas, a sufficiently large macroparticle weight can cause the macroparticles to be strongly coupled and therefore heat due to ACH. Furthermore, it is shown that simulations in reduced dimensions exacerbate these issues.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Influence of Ionization and Beam Quality on Interaction of TW-Peak CO2 Laser with Hydrogen Plasma
3D numerical simulations of the interaction of a powerful CO2 laser with hydrogen jets demonstrating the role of ionization and laser beam quality are presented. Simulations are performed in support of the plasma wakefield accelerator experiments being conducted at the BNL Accelerator Test Facility (ATF). The CO2 laser at BNL ATF has several potential advantages for laser wakefield acceleration compared to widely used solid-state lasers. SPACE, a parallel relativistic Particle-in-Cell code, developed at SBU and BNL, has been used in these studies. A novelty of the code is its set of efficient atomic physics algorithms that compute ionization and recombination rates on the grid and transfer them to particles. The primary goal of the initial BNL experiments was to characterize the plasma density by measuring the sidebands in the spectrum of the probe laser. Simulations, that resolve hydrogen ionization and laser spectra, help explain several trends that were observed in the experiments.
97 MATHEMATICS AND COMPUTING↗
Enabling co-simulation of tokamak plant models and plasma control systems
A system for connecting the Plasma Control System and a model of the tokamak Plant in closed loop co-simulation for plasma control development has been in routine use at DIII-D for more than 20 years and at other fusion labs that use variants of the DIII-D PCS for approximately the last decade. Here, co-simulation refers to the simultaneous execution of two independent codes with the exchange of data - Plant actuator commands and tokamak diagnostic data - between them during execution. Interest in this type of PCS-Plant simulation technology has also been growing recently at other fusion facilities. In fact, use of such closed loop control simulations is assumed to play an even larger role in the development of both the ITER Plasma Control System (PCS) and the experimental operation of the ITER device, where they will be used to support verification/validation of the PCS and also for ITER pulse schedule development and validation. We describe the key use cases that motivate the co-simulation capability and the features that must be provided by the Plasma Control System to support it. These features could be provided by the PCS itself or by a model of the PCS. If the PCS itself is chosen to provide them, there are requirements imposed on its architecture. If a PCS model is chosen, there are requirements imposed on the initial implementation of this simulation as well as long-term consequences for its continued development and maintenance. We describe these issues for each use case and discuss the relative merits of the two choices. Several examples are given illustrating uses of the co-simulation method to address problems of plasma control during the operation of DIII-D and of other devices that use the DIII-D PCS.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Simulation Tool for Dielectric Barrier Discharge Plasma Actuators
Traditional approaches for active flow separation control using dielectric barrier discharge (DBD) plasma actuators are limited to relatively low speed flows and atmospheric conditions. This results in low feasibility of the DBDs for aerospace applications. For active flow control at turbine blades, fixed wings, and rotary wings and on hypersonic vehicles, DBD plasma actuators must perform at a wide range of conditions, including rarified flows and combustion mixtures. An efficient, comprehensive, physically based DBD simulation tool can optimize DBD plasma actuators for different operation conditions. Researchers are developing a DBD plasma actuator simulation tool for a wide range of ambient gas pressures. The tool will treat DBD using either kinetic, fluid, or hybrid models, depending on the DBD operational condition.
Likhanskii, Alexander↗
An implementation of a high-order generalized finite difference method for solving the time-harmonic cold plasma wave equation in toroidal geometry
A high-order physics-informed meshless finite difference numerical technique is introduced for solving the time-harmonic cold plasma wave equation in toroidal geometries, presenting a novel application of the generalized finite difference (GFD) method to plasma wave simulations. The algorithm employs an irregular distribution of computational points, with local point density informed by the shortest wavelength derived from the cold plasma dispersion relation. Numerical stability and robustness are addressed using regularization techniques. The algorithm, implemented for two spatial dimensions, solves for the wave electric field and is demonstrated to achieve convergence rates of $\mathcal{O}$($\mathcal{h}$ $\mathcal{P}$ ). Verification tests reproduce plane wave solutions, and example simulations of ion cyclotron resonance heating and electron cyclotron resonance heating demonstrate its capability, approaching realistic tokamak plasma scenarios. This work contributes to laying a foundation for the GFD method to be used in more sophisticated, optimized, and physically realistic full-wave simulations in time-harmonic plasma wave research.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Improved analytic modeling of the linear Rayleigh–Taylor instability with plasma transport
The Rayleigh–Taylor (R–T) instability with transport effects has long been studied analytically in the neutral fluid regime. How this picture may differ in plasmas has been given less consideration, however. Leveraging previous numerical simulations with plasma viscosity and diffusion of R–T at a binary plasma interface, we show how plasma transport and kinetic effects alter the linear R–T dispersion relation. Our insight informs a much improved analytical dispersion relation, which better matches the simulation data than previously established models. Additionally, we demonstrate the importance of employing asymptotically correct viscosity coefficients in fluid codes for accurately capturing plasma R–T growth.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Gyrokinetic GENE simulations of DIII-D near-edge L-mode plasmas
We present gyrokinetic simulations with the GENE code addressing the near-edge region of an L-mode plasma in the DIII-D tokamak. At radial position ρ = 0.80, simulations with the ion temperature gradient (ITG) increased by 40% above the nominal value give electron and ion heat fluxes that are in simultaneous agreement with the experiment. This gradient increase is consistent with the combined statistical and systematic uncertainty σ of the charge exchange recombination spectroscopy measurements at the 1.6σ level. Multiscale simulations are carried out with a realistic mass ratio and geometry for the first time in the near-edge. These multiscale simulations suggest that the highly unstable ion temperature gradient (ITG) modes of the flux-matched ion-scale simulations suppress electron-scale transport, such that ion-scale simulations are sufficient at this location. At radial position ρ = 0.90, nonlinear simulations show a hybrid state of ITG and trapped electron modes, which was not expected from linear simulations. The nonlinear simulations reproduce the total experimental heat flux with the inclusion of E × B shear effects and an increase in the electron temperature gradient by ∼23%. This gradient increase is compatible with the combined statistical and systematic uncertainty of the Thomson scattering data at the 1.3σ level. These results are consistent with previous findings that gyrokinetic simulations are able to reproduce the experimental heat fluxes by varying input parameters close to their experimental uncertainties, pushing the validation frontier closer to the edge region.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Modeling transient edge plasma transport with dynamic recycling
The work presents numerical simulation studies of the role that dynamic plasma recycling on the main wall and divertor target surfaces plays in transient edge plasma transport phenomena, such as edge localized modes (ELMs). The studies are performed by coupling the edge plasma transport code UEDGE [Rognlien et al., J. Nucl. Mater. 196–198, 347 (1992)] and the wall reaction–diffusion transport code FACE [Smirnov et al., Fusion Sci. Technol. 71, 75 (2017)]. The two-dimensional, time-dependent, two-way coupling of the codes, in a realistic tokamak geometry, is accomplished using the Integrated Plasma Simulator framework [Elwasif et al., in 18th Euromicro Conference on Parallel, Distributed and Network-Based Processing (PDP 2010), Pisa, Italy (IEEE, 2010), pp. 419–427] for all modeled material plasma boundaries. The simulations show that dynamic plasma recycling has substantially different characteristics on the main wall and on the divertor plates. It is demonstrated that during an ELM cycle the outer wall can dynamically absorb and release a number of particles comparable to that expelled by the ELM from the core plasma, by far exceeding the dynamic retention capacity of the divertor surfaces. The resulting evolution of the edge and divertor plasma conditions during an ELM cycle is analyzed.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Gyrokinetic study of the impact of the electron to ion heating ratio on the turbulent diffusion of highly charged impurities
A gyrokinetic study based on numerical and analytical calculations is presented, which computes the dependence of the turbulent diffusion of highly charged impurities on the ratio of the electron to the ion heat flux of the plasma. Nonlinear simulations show that the size of the turbulent diffusion of heavy impurities can vary by one order of magnitude with fixed total heat flux and is an extremely sensitive function of the electron to ion heat flux ratio. Numerical linear calculations are found to reproduce the nonlinear results. Thereby, a quasi-linear analytical approach is used to explain the origin of this dependence.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Numerical simulation study of impurity B transport during real-time B powder injection in EAST
To understand the beneficial effects of impurity boron (B) on the plasma, the transport of B species in both atomic and ionic states in the edge plasma and the ensuing evolution of the background plasma are simulated during real-time B powder injection by using two-dimensional fluid code SOLPS-ITER. The initial background plasma profiles for the simulation are reconstructed based on an EAST upper-single null discharge (shot No. 70601). With the injection rate of B atoms, 4.0 × 10 21 atoms/s, the simulation shows that effective charge number is less than 1.1 at the inner boundary R - R sep = -5.2 cm (background plasma density, 5.1 × 10 19 m -3 ) at the outer mid-plane indicating that the B impurity does not have much influence on the energy confinement. The present work finds that the impurity B exhibits very specific features in the edge plasma, mostly present in the narrow region just outside the separatrix. The particular distribution pattern can reduce the heat flux greatly, leading to the divertor detachment. The divertor detachment can be achieved easily under an injection rate of B powder in a reasonable range. This study suggests that the impurity B has a promising potential application in future tokamaks.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗