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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ~1/$\sqrt{N_p}$ with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly nonequilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Temporal parallelization of edge plasma simulations using the parareal algorithm and the SOLPS code

We show that numerical modelling of edge plasma physics may be successfully parallelized in time. The parareal algorithm has been employed for this purpose and the SOLPS code package coupling the B2.5 finite-volume fluid plasma solver with the kinetic Monte-Carlo neutral code Eirene has been used as a test bed. The complex dynamics of the plasma and neutrals in the scrape-off layer (SOL) region makes this a unique application. It is demonstrated that a significant computational gain (more than an order of magnitude) may be obtained with this technique. The use of the IPS framework for event-based parareal implementation optimizes resource utilization and has been shown to significantly contribute to the computational gain.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Code Verification of a Warm Electron Diode Using the EMPIRE Plasma Simulation Code.

This report documents recent code verification exercises for a warm electron diode problem using the EMPIRE plasma code. This computationally expensive test was performed three times, including two different code versions and two different time integration algorithms, and the resulting code responses were analyzed for convergence to the analytical solution and orders-of-convergence using the StREEQ numerical error estimation tool. Significant deviations from the exact solution and expected orders-of-convergence, as well as changes in code behavior over time and due to choice of time integration algorithm were observed, illustrating the need to fix fundamental code issues as well as additional code verification testing before the code can be relied upon to accurately solve critical problems within its application space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

UCSD Performance in the Edge Plasma Simulation (EPSI) Project. Final Report

This report contains a final report on the activities of UC San Diego PI G.R. Tynan and his collaborators as part of the EPSI Project, that was led by Dr. C.S. Chang, from PPPL. As a part of our work, we carried out several experiments on the ALCATOR C-­MOD tokamak device, aimed at unraveling the “trigger” or cause of the spontaneous transition from low-­mode confinement (L-­mode) to high confinement (H-­mode) that is universally observed in tokamak devices, and is planned for use in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

A pairwise nuclear fusion algorithm for weighted particle-in-cell plasma simulations

In this work, a pairwise nuclear fusion algorithm for arbitrarily weighted macroparticles in a particle-in-cell simulation is described. The method is benchmarked in situations with like-particles, D(d,n) 3 He, unlike-particles, D(t,n) 4 He, thermonuclear plasmas, beam-target fusion, and for large difference in macroparticle weights. Studies of the required number of macroparticles in thermonuclear plasmas show that 100–1000 macroparticles are required to achieve repeatability of yields around 10%, likewise 10 4 –10 5 for 1%, depending on the fusion interaction and ion temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1972

Development of RF plasma simulations of in-reactor tests of small models of the nuclear light bulb fuel region

Experiments were conducted to develop test configurations and technology necessary to simulate the thermal environment and fuel region expected to exist in in-reactor tests of small models of nuclear light bulb configurations. Particular emphasis was directed at rf plasma tests of approximately full-scale models of an in-reactor cell suitable for tests in Los Alamos Scientific Laboratory's Nuclear Furnace. The in-reactor tests will involve vortex-stabilized fissioning uranium plasmas of approximately 200-kW power, 500-atm pressure and equivalent black-body radiating temperatures between 3220 and 3510 K.

Roman, W. C.↗
NASA NTRS2018

Investigations of Spot-to-plume Mode Transition in a Hollow Cathode Discharge Using 2-D Axisymmetric Plasma Simulations

After decades of laboratory investigations that provided invaluable measurements and insight, the physics behind the transition from spot to plume modes in hollow cathodes remains one of the longest standing theoretical problems in electric propulsion. This has prohibited the development of ab initio models that allow for the prediction of the transition across different cathodes and operating conditions. Since the beginning of its development over a decade ago, simulations with the 2-D axisymmetric Orificed Cathode (OrCa2D) code have helped elucidate a wide range of processes in hollow cathode discharges. However, the code has never been used to investigate the onset of plume mode. We present results from the first OrCa2D simulations of a 25-A LaB6 cathode for a range of flow rates (5-20 sccm) in which transition from spot to plume modes is known to occur. The cathode in this study was one of the two technologies considered for the 12.5 kW Hall Effect Rocket with Magnetic Shielding (HERMeS) and operates nominally at 21 A and 14.8 sccm. The simulations capture the characteristic rise of the peak-to-peak amplitude in the keeper voltage oscillations and underscore the significance of the plume neutral gas in the transition. The plasma inside the cathode is found to be relatively quiescent throughout the transition, in agreement with previous experimental observations. The computed keeper voltage fluctuations at low flow rates (<8 sccm) are found to be driven by oscillations of the same frequency in the plasma plume with the following main characteristics: (1) they are of low frequency (<10 kHz), and associated with small longitudinal motion in the direction of the applied magnetic field, (2) they occur in a region of the plume where the neutral gas provided by the cathode has been fully depleted, and (3) they have a (small) wave velocity of about 100 m/s, which is at least ~10 smaller than the drift, thermal and acoustic speeds of the ions. At 8 sccm, when the transition to the large-amplitude oscillations begins, the ionization frequency in the neutral-depleted plume region ranges ~2-100 kHz. The frequency of the oscillations in the plasma (and keeper voltage) is found to be equal to the ionization frequency (~5 kHz) at the center of this region. The findings suggest that the transition to plume mode is driven by ionization processes in the near-plume of the cathode, in line with previous conjectures that were based solely on laboratory observations.

Polk, James E.↗
DOE OSTI2021

Challenges and opportunities in verification and validation of low temperature plasma simulations and experiments

This paper describes the verification and validation (V&V) framework developed for the stochastic Particle-in-Cell, Direct Simulation Monte Carlo code Aleph. An ideal framework for V&V from the viewpoint of the authors is described where a physics problem is defined, and relevant physics models and parameters to the defined problem are assessed and captured in a Phenomena Identification and Ranking Table (PIRT). Furthermore, numerous V&V examples guided by the PIRT for a simple gas discharge are shown to demonstrate the V&V process applied to a real-world simulation tool with the overall goal to demonstrably increase the confidence in the results for the simulation tool and its predictive capability. Although many examples are provided here to demonstrate elements of the framework, the primary goal of this work is to introduce this framework and not to provide a fully complete implementation, which would be a much longer document. Comparisons and contrasts are made to more usual approaches to V&V, and techniques new to the low-temperature plasma community are introduced. Specific challenges relating to the sufficiency of available data (e.g., cross sections), the limits of ad hoc validation approaches, the additional difficulty of utilizing a stochastic simulation tool, and the extreme cost of formal validation are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A corrected method for Coulomb scattering in arbitrarily weighted particle-in-cell plasma simulations

A physically and statistically accurate method to perform binary, pairwise Coulomb collisions for arbitrarily-weighted macroparticles in particle-in-cell simulations is described. The method is shown to conserve global energy and momentum on average and to produce physically accurate scattering frequencies. This is justified using the statistically expected values of the method and through a series of numerical tests showing that the method is valid independently of macroparticle weight. This novel method is required due to an error in previous work on binary scattering, initially by Nanbu and Yonemura [J. Comp. Phys. 145, 639 (1998)] and applied to subsequent work, that causes inaccuracies in the scattering frequencies for unequally weighted simulations; a problem that occurs both for scattering within the group and between different groups.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2019

Plasma simulations for the assessment of pole erosion in the Magnetically Shielded Miniature Hall Thruster (MaSMi)

Numerical simulations have been conducted to assess the erosion rates at the pole covers of the magnetically shielded miniature (MaSMi) Hall thruster. MaSMi is part of the Ascendant Sub-kW Transcelestial Electric Propulsion System (ASTRAEUS), whose design requirements demand a xenon throughput of 100 kg, thus enabling deep space missions using SmallSats. The results of our simulations show that the proposed thickness of the graphite pole covers will not be completely eroded at end of life. We find that the erosion rates are larger near the corners of the covers, which are directly exposed to the plasma in the acceleration channel.

Chaplin, Vernon H.↗
DOE OSTI2020

ContourNet: Salient Local Contour Identification for Blob Detection in Fusion Plasma Simulation Data

We present ContourNet, a deep learning approach to identify salient local isocontours as blobs in large-scale 5D gyrokinetic tokamak simulation data. Blobs-regions of high turbulence that run along the edge wall down toward the diverter and can damage the tokamak-are non-well-defined features but have been empirically localized by isocontours in 2D normalized fluctuating density fields. The key of our study is to train ContourNet to follow the empirical rules to detect blobs over the time-varying simulation data. The architecture of ContourNet is a convolutional neural segmentation network: the inputs are the density field and a rasterized isocontour; the output is a set of isocontour encircling blobs. At the training stage, we feed the network with manually identified isocontours and propagated labels. At the inference stage, we extract isocontours from the segmented blob regions. Results show that our approach can achieve both high accuracy and performance, which enables scientists to understand the blob dynamics influencing the confinement of the plasma.

Blob detection↗