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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Transition characteristics and electron kinetics in microhollow cathode discharges

We demonstrate the transition characteristics and electron kinetics of microdischarges in a microgap with a cathode having microstructures using two-dimensional particle-in-cell/Monte Carlo collision (2d3v) simulations. It is found that the microdischarge is inside the hollow cavity at lower pressures, forming microhollow cathode discharges (MHCDs), while the dominant discharge moves outside the cavity at higher pressures, locating above the cathode rectangular protrusion tip. The spatial distributions of the microdischarge parameters (e.g., electric potential and electron density) are presented by capturing the transition characteristics. The electron kinetics of the microdischarges are examined based on the moment analysis of the electron Boltzmann equation. The collisional and non-collisional components of the electron power absorption are decomposed, which were seldomly investigated for microdischarges previously. The results provide fundamental insights into MHCD formation with a structured electrode at varying pressure conditions, which could be beneficial for the design and fabrication of microplasma devices in practical applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

The evolution of the ion diffusion region during collisionless magnetic reconnection in a force-free current sheet

Two-dimensional particle-in-cell simulation is performed to investigate magnetic reconnection in a force-free current sheet. The results show that the evolution of the ion diffusion region has two different phases. In the first phase, the electrons flow toward the X line along one pair of separatrices and away from the X line along the other pair of separatrices. Therefore, in the ion diffusion region, a distorted quadrupole structure of the out-of-plane magnetic field is formed, which is similar to that of a typical guide field reconnection in the Harris current sheet. In the second phase, the electrons move toward the X line along the separatrices and then flow away from the X line at the inner side of the separatrices. In the ion diffusion region, the out-of-plane magnetic field exhibits a characteristic quadrupole pattern with a good symmetry, which is similar to that of antiparallel reconnection in the Harris current sheet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

On the electron dynamics during island coalescence in asymmetric magnetic reconnection

We present an analysis of the electron dynamics during rapid island merging in asymmetric magnetic reconnection. We consider a doubly periodic system with two asymmetric transitions. The upper layer is an asymmetric Harris sheet of finite width perturbed initially to promote a single reconnection site. The lower layer is a tangential discontinuity that promotes the formation of many X-points, separated by rapidly merging islands. Across both layers, the magnetic field and the density have a strong jump, but the pressure is held constant. Our analysis focuses on the consequences of electron energization during island coalescence. We focus first on the parallel and perpendicular components of the electron temperature to establish the presence of possible anisotropies and non-gyrotropies. Thanks to the direct comparison between the two different layers simulated, we can distinguish three main types of behavior characteristic of three different regions of interest. The first type represents the regions where traditional asymmetric reconnections take place without involving island merging. The second type of regions instead shows reconnection events between two merging islands. Finally, the third regions identify the regions between two diverging island and where typical signature of reconnection is not observed. Electrons in these latter regions additionally show a flat-top distribution resulting from the saturation of a two-stream instability generated by the two interacting electron beams from the two nearest reconnection points. Finally, the analysis of agyrotropy shows the presence of a distinct double structure laying all over the lower side facing the higher magnetic field region. This structure becomes quadrupolar in the proximity of the regions of the third type. The distinguishing features found for the three types of regions investigated provide clear indicators to the recently launched Magnetospheric Multiscale NASA mission for investigating magnetopause reconnection involving multiple islands.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Global fluid turbulence simulations in the scrape-off layer of a stellarator island divertor

Isothermal fluid turbulence simulations have been performed in the edge and scrape-off layer (SOL) of an analytic stellarator configuration with an island divertor, thereby providing numerical insight into edge turbulence in regions around islands in a stellarator. The steady-state transport follows the a curvature drive that is inverse to the major radius ( $1/R$ ) toward the outboard side, but large fluctuations are present throughout the island divertor region, with the average wavelength of similar size to the island width. The system exhibits a prominent $m=2$ , $n=5$ mode, where m is the poloidal mode number and n is the toroidal mode number, although other modes are present. The amplitude and radial extent of the density fluctuations are similar throughout the edge and SOL, but can decrease near island O-points. The fluctuations exhibit a predominantly positive skewness on the outboard midplane, indicating blob-like perturbations for the transport into the outer SOL. It is determined that a point on the separatrix is generally more correlated with regions outside of the SOL than a nearby reference point which does not lie on the separatrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Particle-in-cell simulation of multipactor discharge on a dielectric in a parallel-plate waveguide

An original 2D3V (two-dimensional in coordinate space and three-dimensional in velocity space) particle-in-cell code has been developed for simulation of multipactor discharge on a dielectric in a parallelplate metal waveguide with allowance for secondary electron emission (SEE) from the dielectric surface and waveguide walls, finite temperature of secondary electrons, electron space charge, and elastic and inelastic scattering of electrons from the dielectric and metal surfaces. The code allows one to simulate all stages of the multipactor discharge, from the onset of the electron avalanche to saturation. It is shown that the threshold for the excitation of a single-surface multipactor on a dielectric placed in a low-profile waveguide with absorbing walls increases as compared to that in the case of an unbounded dielectric surface due to escape of electrons onto the waveguide walls. It is found that, depending on the microwave field amplitude and the SEE characteristics of the waveguide walls, the multipactor may operate in two modes. In the first mode, which takes place at relatively low microwave amplitudes, a single-surface multipactor develops only on the dielectric, the surface of which acquires a positively potential with respect to the waveguide walls. In the second mode, which occurs at sufficiently high microwave intensities, a single-surface multipactor on the dielectric and a two-surface multipactor between the waveguide walls operate simultaneously. In this case, both the dielectric surface and the interwall space acquire a negative potential. It is shown that electron scattering from the dielectric surface and waveguide walls results in the appearance of high-energy tails in the electron distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Effects of a Guide Field on the Larmor Electric Field and Upstream Electron Temperature Anisotropy in Collisionless Asymmetric Magnetic Reconnection

We perform the first study of the properties of the Larmor electric field (LEF) in collisionless asymmetric magnetic reconnection in the presence of an out-of-plane (guide) magnetic field for different sets of representative upstream parameters at Earth’s dayside magnetopause with an ion temperature greater than the electron temperature (the ion-to-electron temperature ratio fixed at 2) using two-dimensional particle-in-cell simulations. We show that the LEF does persist in the presence of a guide field. We study how the LEF thickness and strength change as a function of guide field and the magnetospheric temperature and reconnecting magnetic field strength. We find that the thickness of the LEF structure decreases, while its magnitude increases when a guide field is added to the reconnecting magnetic field. The added guide field makes the Larmor radius smaller, so the scaling with the magnetospheric ion Larmor radius is similar to that reported for the case without a guide field. Note, however, that the physics causing the LEF is not well understood, so future work in other parameter regimes is needed to fully predict the LEF for arbitrary conditions. We also find that a previously reported upstream electron temperature anisotropy arises in the vicinity of the LEF region both with and without a guide field. We argue that the generation of the anisotropy is linked to the existence of the LEF. The LEF can be used in combination with the electron temperature anisotropy as a signature to effectively identify dayside reconnection sites in observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Diffusive shock re-acceleration

Here, we have performed two-dimensional hybrid simulations of non-relativistic collisionless shocks in the presence of pre-existing energetic particles (‘seeds’); such a study applies, for instance, to the re-acceleration of galactic cosmic rays (CRs) in supernova remnant (SNR) shocks and solar wind energetic particles in heliospheric shocks. Energetic particles can be effectively reflected and accelerated regardless of shock inclination via a process that we call diffusive shock re-acceleration. We determine that re-accelerated seeds can drive the streaming instability in the shock upstream and produce effective magnetic field amplification. This can later trigger the injection of thermal protons even at oblique shocks that ordinarily cannot inject thermal particles. We characterize the current in reflected seeds, finding that it tends to a universal value $J\simeq en_{\text{CR}}v_{\text{sh}}$ , where $en_{\text{CR}}$ is the seed charge density and $v_{\text{sh}}$ is the shock velocity. When applying our results to SNRs, we find that the re-acceleration of galactic CRs can excite the Bell instability to nonlinear levels in less than ${\sim}10~\text{yr}$ , thereby providing a minimum level of magnetic field amplification for any SNR shock. Finally, we discuss the relevance of diffusive shock re-acceleration also for other environments, such as heliospheric shocks, galactic superbubbles and clusters of galaxies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

He{sup 2+} HEATING VIA PARAMETRIC INSTABILITIES OF PARALLEL PROPAGATING ALFVÉN WAVES WITH AN INCOHERENT SPECTRUM

The preferential heating of heavy ions in the solar corona and solar wind has been a long-standing hot topic. In this paper we use a one-dimensional hybrid simulation model to investigate the heating of He{sup 2+} particles during the parametric instabilities of parallel propagating Alfvén waves with an incoherent spectrum. The evolution of the parametric instabilities has two stages and involves the heavy ion heating during the entire evolution. In the first stage, the density fluctuations are generated by the modulation of the pump Alfvén waves with a spectrum, which then results in rapid coupling with the pump Alfvén waves and the cascade of the magnetic fluctuations. In the second stage, each pump Alfvén wave decays into a forward density mode and a backward daughter Alfvén mode, which is similar to that of a monochromatic pump Alfvén wave. In both stages the perpendicular heating of He{sup 2+} particles occurs. This is caused by the cyclotron resonance between He{sup 2+} particles and the high-frequency magnetic fluctuations, whereas the Landau resonance between He{sup 2+} particles and the density fluctuations leads to the parallel heating of He{sup 2+} particles. The influence of the drift velocity between the protons and the He{sup 2+} particles on the heating of He{sup 2+} particles is also discussed in this paper.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Experimental room temperature hohlraum performance study on the National Ignition Facility [Experimental evidence for improved performance in room temperature hohlraums on the National Ignition Facility]

Room temperature or “warm” (273 K) indirect drive hohlraum experiments have been conducted on the National Ignition Facility with laser energies up to 1.26 MJ and compared to similar cryogenic or “cryo” (~20 K) experiments. Warm experiments use neopentane (C 5 H 12 ) as the low pressure hohlraum fill gas instead of helium, and propane (C 3 H 8 ) to replace the cryogenic DT or DHe3 capsule fill. The increased average Z of the hohlraum fill leads to increased inverse bremsstrahlung absorption and an overall hotter hohlraum plasma in simulations. The cross beam energy transfer (CBET) from outer laser beams (pointed toward the laser entrance hole) to inner beams (pointed at the equator) was inferred indirectly from measurements of Stimulated Raman Scattering (SRS). These experiments show that a similar hot spot self-emission shape can be produced with less CBET in warm hohlraums. The measured inner cone SRS reflectivity (as a fraction of incident power neglecting CBET) is ~2.5× less in warm than cryo shots with similar hot spot shapes, due to a less need for CBET. The measured outer-beam stimulated the Brillouin scattering power that was higher in the warm shots, leading to a ceiling on power to avoid the optics damage. These measurements also show that the CBET induced by the flow where the beams cross can be effectively mitigated by a 1.5 Å wavelength shift between the inner and outer beams. A smaller scale direct comparison indicates that warm shots give a more prolate implosion than cryo shots with the same wavelength shift and pulse shape. Lastly, the peak radiation temperature was found to be between 5 and 7 eV higher in the warm than the corresponding cryo experiments after accounting for differences in backscatter.

42 ENGINEERING↗
DOE OSTI2018

Simulations of the Stable Expansion of a Monoelement Three-Component Explosive-Emission Plasma

The results of numerical simulations of the stable expansion stage of a micro explosive plasma plume are presented. A three-liquid hydrodynamic model was used with corrections allowing its application to the consistent description of the transfer of particles in a monoelement plasma with ions with two different charges and of the electron transfer in a changing plasma–collector vacuum gap. The evolution of concentration profiles, the drift velocity, and electron and ion temperatures in plasma are analyzed depending on conditions of a collector. It is shown in the model that the velocity separation of ions with different charges occurs efficiently in the absence of an external electric field, whereas in a current-carrying plasma, the ions with different charges move with almost identical velocities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Transitioning from File-Based HPC Workflows to Streaming Data Pipelines with openPMD and ADIOS2

This paper aims to create a transition path from file-based IO to streaming-based workflows for scientific applications in an HPC environment. By using the openPMP-api, traditional workflows limited by filesystem bottlenecks can be overcome and flexibly extended for in situ analysis. The openPMD-api is a library for the description of scientific data according to the Open Standard for Particle-Mesh Data (openPMD). Its approach towards recent challenges posed by hardware heterogeneity lies in the decoupling of data description in domain sciences, such as plasma physics simulations, from concrete implementations in hardware and IO. The streaming backend is provided by the ADIOS2 framework, developed at Oak Ridge National Laboratory. This paper surveys two openPMD-based loosely-coupled setups to demonstrate flexible applicability and to evaluate performance. In loose coupling, as opposed to tight coupling, two (or more) applications are executed separately, e.g. in individual MPI contexts, yet cooperate by exchanging data. This way, a streaming-based workflow allows for standalone codes instead of tightly-coupled plugins, using a unified streaming-aware API and leveraging high-speed communication infrastructure available in modern compute clusters for massive data exchange. We determine new challenges in resource allocation and in the need of strategies for a flexible data distribution, demonstrating their influence on efficiency and scaling on the Summit compute system. The presented setups show the potential for a more flexible use of compute resources brought by streaming IO as well as the ability to increase throughput by avoiding filesystem bottlenecks.

Poeschel, Franz↗
DOE OSTI2018

Energy levels and radiative rates for transitions in Cr-like Kr XIII, Tc XX and Xe XXXI

Energy levels, radiative transition probabilities, Landé g-factors, and radiative lifetimes are reported for 144 levels of the 3p{sup 6}3d{sup 6} and 3p{sup 5}3d{sup 7} configurations of the Cr-like ions Kr XIII, Tc XX and Xe XXXI. Since to our knowledge no prior data for energy levels or A-values, experimental or theoretical, are available in the literature for these ions, calculations have been performed with three independent atomic structure codes, namely GRASP, FAC and a code implementing the QR approach, with increasing configuration interaction (CI) basis with each code. This is to establish the accuracy of the calculated results so that data can be confidently applied in the modeling of plasmas. Our recommended data are determined within the QR approach, because not only has the largest CI basis been included in this calculation, the reliability has already been established in our earlier work for other Cr-like ions, for which some experimental and theoretical data were available for comparison. Results for radiative rates (A-values) and transition wavelengths (λ) are listed for all strong E1, E2 and M1 transitions among the 144 levels of each ion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Calculations for energies, transition rates, and lifetimes in Al-like Kr XXIV

Using the second-order many-body perturbation theory (MBPT) method, a complete and accurate data set of excitation energies, lifetimes, wavelengths, and electric dipole (E1), magnetic dipole (M1), electric quadrupole (E2), and magnetic quadrupole (M2) line strengths, transition rates, and oscillator strengths for the lowest 880 levels arising from the 3l{sup 3} (0≤l≤2), 3l{sup 2}4l{sup ′} (0≤l≤2, 0≤l{sup ′}≤3), 3s{sup 2}5l (0≤l≤4), 3p{sup 2}5l (0≤l≤1), and 3s3p5l (0≤l≤4) configurations in Al-like Kr XXIV is provided. Comparisons are made with available experimental and theoretical results. Our calculated energies are expected to be accurate enough to facilitate identifications of observed lines involving the n=4,5 levels. The complete data set is also useful for modeling and diagnosing fusion plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

An efficient transport solver for tokamak plasmas

A simple approach to efficiently solve a coupled set of 1-D diffusion-type transport equations with a stiff transport model for tokamak plasmas is presented based on the 4th order accurate Interpolated Differential Operator scheme along with a nonlinear iteration method derived from a root-finding algorithm. Here, numerical tests using the Trapped Gyro-Landau-Fluid model show that the presented high order method provides an accurate transport solution using a small number of grid points with robust nonlinear convergence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Multidisciplinary benchmarks of a conservative spectral solver for the nonlinear Boltzmann equation

The Boltzmann equation describes the evolution of the phase-space probability distribution of classical particles under binary collisions. Approximations to it underlie the basis for several scholarly fields, including aerodynamics and plasma physics. While these approximations are appropriate in their respective domains, they can be violated in niche but diverse applications which require direct numerical solution of the original nonlinear Boltzmann equation. An expanded implementation of the Galerkin–Petrov conservative spectral algorithm is employed to study a wide variety of physical problems. Enabled by distributed precomputation, solutions of the spatially homogeneous Boltzmann equation can be achieved in seconds on modern personal hardware, while spatially-inhomogeneous problems are solvable in minutes. Here, several benchmarks are presented focusing on accuracy compared to both analytic theoretical predictions and other Boltzmann solvers. These benchmarks span several physical domains including weakly ionized plasma, gaseous fluids, and atomic-plasma interaction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

A hyperbolic-equation system approach for magnetized electron fluids in quasi-neutral plasmas

A new approach using a hyperbolic-equation system (HES) is proposed to solve for the electron fluids in quasi-neutral plasmas. The HES approach avoids treatments of cross-diffusion terms which cause numerical instabilities in conventional approaches using an elliptic equation (EE). A test calculation reveals that the HES approach can robustly solve problems of strong magnetic confinement by using an upwind method. The computation time of the HES approach is compared with that of the EE approach in terms of the size of the problem and the strength of magnetic confinement. The results indicate that the HES approach can be used to solve problems in a simple structured mesh without increasing computational time compared to the EE approach and that it features fast convergence in conditions of strong magnetic confinement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗