PIC-DSMC Simulations of Plasma Dynamics and Expansion During Vacuum Arc Initiation
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Abstract Not Provided
Abstract not provided.
Abstract Not Provided
Abstract not provided.
Explore the source record for details and available documents.
Abstract not provided.
Explore the source record for details and available documents.
The properties of the turbulence that develops in the outflows of magnetic reconnection have been investigated using self-consistent plasma simulations, in three dimensions. As commonly observed in space plasmas, magnetic reconnection is characterized by the presence of turbulence. In this study, we provide a direct comparison of our simulations with reported observations of reconnection events in the magnetotail, investigating the properties of the electromagnetic field and the energy conversion mechanisms. In particular, simulations show the development of a turbulent cascade consistent with spacecraft observations, statistics of the dissipation mechanisms in the turbulent outflows similar to the ones observed in reconnection jets in the magnetotail, and that the properties of turbulence vary as a function of the distance from the reconnecting X-line.
We develop new numerical schemes for Vlasov–Poisson equations with high-order accuracy. Our methods are based on a spatially monotonicity-preserving (MP) scheme and are modified suitably so that the positivity of the distribution function is also preserved. We adopt an efficient semi-Lagrangian time integration scheme that is more accurate and computationally less expensive than the three-stage TVD Runge–Kutta integration. We apply our spatially fifth- and seventh-order schemes to a suite of simulations of collisionless self-gravitating systems and electrostatic plasma simulations, including linear and nonlinear Landau damping in one dimension and Vlasov–Poisson simulations in a six-dimensional phase space. The high-order schemes achieve a significantly improved accuracy in comparison with the third-order positive-flux-conserved scheme adopted in our previous study. With the semi-Lagrangian time integration, the computational cost of our high-order schemes does not significantly increase, but remains roughly the same as that of the third-order scheme. Vlasov–Poisson simulations on 128{sup 3}×128{sup 3} mesh grids have been successfully performed on a massively parallel computer.
Damped Lyα (DLA) and sub-DLA absorbers in quasar spectra provide the most sensitive tools for measuring the element abundances of distant galaxies. The estimation of abundances from absorption lines depends sensitively on the accuracy of the atomic data used. We have started a project to produce new atomic spectroscopic parameters for optical and UV spectral lines using state-of-the-art computer codes employing a very broad configuration interaction (CI) basis. Here we report our results for Zn ii, an ion used widely in studies of the interstellar medium (ISM) as well as DLAs and sub-DLAs. We report new calculations of many energy levels of Zn ii and the line strengths of the resulting radiative transitions. Our calculations use the CI approach within a numerical Hartree–Fock framework. We use both nonrelativistic and quasi-relativistic one-electron radial orbitals. We have incorporated the results of these atomic calculations into the plasma simulation code Cloudy and applied them to a lab plasma and examples of a DLA and a sub-DLA. Our values of the Zn ii λλ 2026, 2062 oscillator strengths are higher than previous values by 0.10 dex. The Cloudy calculations for representative absorbers with the revised Zn atomic data imply ionization corrections lower than calculated earlier by 0.05 dex. The new results imply that Zn metallicities should be lower by 0.1 dex for DLAs and by 0.13–0.15 dex for sub-DLAs than in past studies. Our results can be applied to other studies of Zn ii in the Galactic and extragalactic ISM.
Fusion edge plasmas can be far from thermal equilibrium and require the use of a non-linear collision operator for accurate numerical simulations. In this article, the non-linear single-species Fokker–Planck–Landau collision operator developed by Yoon and Chang (2014) [9] is generalized to include multiple particle species. The finite volume discretization used in this work naturally yields exact conservation of mass, momentum, and energy. The implementation of this new non-linear Fokker–Planck–Landau operator in the gyrokinetic particle-in-cell codes XGC1 and XGCa is described and results of a verification study are discussed. Finally, the numerical techniques that make our non-linear collision operator viable on high-performance computing systems are described, including specialized load balancing algorithms and nested OpenMP parallelization. The collision operator's good weak and strong scaling behavior are shown.
This work demonstrates that simulations of advanced burning plasma operation scenarios can be successfully parallelized in time using the parareal algorithm. CORSICA - an advanced operation scenario code for tokamak plasmas is used as a test case. This is a unique application since the parareal algorithm has so far been applied to relatively much simpler systems except for the case of turbulence. In the present application, a computational gain of an order of magnitude has been achieved which is extremely promising. A successful implementation of the Parareal algorithm to codes like CORSICA ushers in the possibility of time efficient simulations of ITER plasmas.
This milestone work baselines electromagnetic particle-in-cell capability of the EMPIRE plasma simulation code to model key processes germane to the physics of electrode plasmas arising in magnetically-insulated transmission lines operating at or near 20 MA. This evaluation is done so through the provision of benchmark verification problems designed to exercise the individual and combined physics models on a small-scale surrogate geometry for the final-feed-to-load region of the Z accelerator under representative operating conditions. In this report, we overview our test designs, and present a portfolio of simulation results along with performance assessments which altogether establish state-of-the-art. In particular, two main verification categories are covered this report: (1) Z-relevant desorption physics (Temkin isotherm), and (2) two approaches to simulate electrode plasma creation and dynamics (automatic creation versus self-consistent creation through direct simulation Monte Carlo collisions).