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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Kinetic Simulation of Collisional Magnetized Plasmas with Semi-implicit Time Integration

Plasmas with varying collisionalities occur in many applications, such as tokamak edge regions, where the flows are characterized by significant variations in density and temperature. While a kinetic model is necessary for weakly-collisional high-temperature plasmas, high collisionality in colder regions render the equations numerically stiff due to disparate time scales. In this study, we propose an implicit–explicit algorithm for such cases, where the collisional term is integrated implicitly in time, while the advective term is integrated explicitly in time, thus allowing time step sizes that are comparable to the advective time scales. This partitioning results in a more efficient algorithm than those using explicit time integrators, where the time step sizes are constrained by the stiff collisional time scales. Finally, we implement semi-implicit additive Runge–Kutta methods in COGENT, a high-order finite-volume gyrokinetic code and test the accuracy, convergence, and computational cost of these semi-implicit methods for test cases with highly-collisional plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Premixed Combustion Simulations with a Self-Consistent Plasma Model for Initiation

Combustion simulations of H2-O2 ignition are presented here, with a self-consistent plasma fluid model for ignition initiation. The plasma fluid equations for a nanosecond pulsed discharge are solved and coupled with the governing equations of combustion. The discharge operates with the propagation of cathode directed streamer, with radical species produced at streamer heads. These radical species play an important role in the ignition process. The streamer propagation speeds and radical production rates were found to be sensitive to gas temperature and fuel-oxidizer equivalence ratio. The oxygen radical production rates strongly depend on equivalence ratio and subsequently results in faster ignition of leaner mixtures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Final Technical Report for Grant No. DE-FC02-01ER54648, SciDAC Center for Simulation of Wave-Plasma Interactions

Harnessing the energy that is released in fusion reactions would provide a safe and abundant source of power to meet the growing energy needs of the world population. The next step toward the development of fusion as a practical energy source is the construction of ITER, a device capable of producing and controlling the high performance plasma required for self-sustaining fusion reactions, or “burning” plasma. The input power required to drive the ITER plasma into the burning regime will be supplied primarily with a combination of external power from radio frequency waves in the ion cyclotron range of frequencies and energetic ions from neutral beam injection sources, in addition to internally-generated Ohmic heating from the induced plasma current that also serves to create the magnetic equilibrium for the discharge. Since the success of the ITER project depends critically on the ability to create and maintain burning plasma conditions, it is absolutely necessary to have physics-based models that can accurately simulate the RF processes that affect the dynamical evolution of the ITER discharge. The Center for Simulation of Wave-Plasma Interactions (CSWPI) is a multi-institutional collaboration that has conducted ongoing research aimed at developing: (1) Coupled core-to-edge simulations that will lead to an increased understanding of parasitic losses of the applied RF power in the boundary plasma between the RF antenna and the core plasma; (2) Development of models for core interactions of RF waves with energetic electrons and ions (including fusion alpha particles and fast neutral beam ions) that include a more accurate representation of the particle dynamics in the combined equilibrium and wave fields; and (3) Development of improved algorithms that will take advantage of massively parallel computing platforms at the petascale level and beyond to achieve the needed physics, resolution, and/or statistics to address these issues. This research reported in this final technical report is focused on further verification and validation of simulation models for interactions of RF waves with energetic particles in the core with a long term goal of assessing their readiness for inclusion in integrated modeling frameworks and utilizing the world’s largest supercomputers to quantitatively predict the detailed electric field structure in the vicinity of the launching structures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Reformulation of plasma spectral line broadening with computer simulations

Spectral line broadening by plasmas can be computed by solving the equation of motion for the dipole of the radiating system perturbed by a fluctuating radiator–plasma interaction obtained from computer simulations. These calculations usually described the radiator by a localized charge distribution of negligible size and the plasma by an electric potential plus inhomogeneities at the radiator nucleus. Using an alternative approach, it is shown that starting from the full Coulomb interaction and then introducing approximations leads to expressions that naturally take advantage of readily available data from atomic packages. Applying this alternative approach to screened Coulomb interactions in so-called trivial molecular dynamics to account for neglected plasma correlations leads to a modified multipole expansion compared to the standard method. Issues associated with the ad hoc introduction of the screened Coulomb potential are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Quantum Information for Fusion Energy Sciences (Final Technical Report)

The simulation of plasma dynamics is a critical area of Fusion Energy Sciences (FES) due to it’s usefulness in predicting, controlling, and confining plasmas in the context of potential fusion reactors. The simulation of plasmas is a computationally difficult problem in both classical and quantum physics, motivating investigation into the potential of quantum computers to simulate these systems. This project took several concrete steps towards this goal by developing tools for improving the control, characterization, and calibration of quantum gates on a superconducting quantum computer, developing error suppression and mitigation tools to reduce errors on the quantum computer, and utilizing these advancements to simulate reduced models of plasma dynamics on the quantum computer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Multispecies density peaking in gyrokinetic turbulence simulations of low collisionality Alcator C-Mod plasmas

Peaked density profiles in low-collisionality AUG and JET H-mode plasmas are probably caused by a turbulently driven particle pinch, and Alcator C-Mod experiments confirmed that collisionality is a critical parameter. Density peaking in reactors could produce a number of important effects, some beneficial, such as enhanced fusion power and transport of fuel ions from the edge to the core, while others are undesirable, such as lower beta limits, reduced radiation from the plasma edge, and consequently higher divertor heat loads. Fundamental understanding of the pinch will enable planning to optimize these impacts. We show that density peaking is predicted by nonlinear gyrokinetic turbulence simulations based on measured profile data from low collisionality H-mode plasma in Alcator C-Mod. Multiple ion species are included to determine whether hydrogenic density peaking has an isotope dependence or is influenced by typical levels of low-Z impurities, and whether impurity density peaking depends on the species. We find that the deuterium density profile is slightly more peaked than that of hydrogen, and that experimentally relevant levels of boron have no appreciable effect on hydrogenic density peaking. The ratio of density at r/a = 0.44 to that at r/a = 0.74 is 1.2 for the majority D and minority H ions (and for electrons), and increases with impurity Z: 1.1 for helium, 1.15 for boron, 1.3 for neon, 1.4 for argon, and 1.5 for molybdenum. The ion temperature profile is varied to match better the predicted heat flux with the experimental transport analysis, but the resulting factor of two change in heat transport has only a weak effect on the predicted density peaking.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Effects of a magnetic field on hot electron generation from laser-plasma instabilities

We performed two-dimensional simulations of turbulent laser-plasma instabilities in the presence and absence of external magnetic fields using the Laser Plasma Simulation Environment (LPSE) code. The results demonstrate that, in the presence of a magnetic field, the transition from ballistic to gyrating electron motion enhances the energy transfer from electron plasma waves to the electron population. Although stronger magnetic fields produce a larger population of hot electrons, these electrons tend to remain confined near the quarter-critical density, where the instabilities also localize, thereby reducing the potential for hot electron transport deeper into the target. Additionally, we present a scaling analysis that quantifies hot electron generation as a function of plasma electron temperature, density scale length, and applied magnetic field strength. These findings may have important applications for controlling hot electron flux and mitigating preheat in inertial confinement fusion targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Applying the Gyrokinetic Formulation of Magnetized Plasma Dynamics to Z-pinch plasmas

This project adapts and applies the gyrokinetic Eulerian code COGENT to an emerging area of highly-magnetized Z-pinch plasma systems. Simulation efforts for these plasmas currently use either magneto-hydrodynamic (MHD) fluid codes, which suffer from fidelity being restricted to narrow plasma conditions, or fully kinetic particle-in-cell (PIC) methods, which suffer from the run-time needed to resolve the very small spatial and temporal scales required for the plasma conditions. The progress with achieving detailed theoretical understanding of Z-pinch plasma systems can be greatly expedited using some combination of gyrokinetic and drift-ordered MHD models, developed in the tokamak community and implemented in the COGENT code. These formulations are applicable to a much larger parameter space than standard MHD models, while providing a substantial speed-up in run-time over fully kinetic codes. In particular, the gyrokinetic formalism is successfully applied to a flow-stabilized Z pinch (FSZP) system – a promising fusion concept where a sheared axial flow is used to provide stability against large-scale ideal-MHD modes. Results of gyrokinetic simulations with COGENT agree well with those from fully kinetic simulations, while requiring ~500x less run-time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1978

Monte Carlo simulation of the classical two-dimensional one component plasma

Monte Carlo simulation, lattice dynamics in the harmonic approximation, and solution of the hypernetted chain equation were used to study the classical two-dimensional one component plasma. The system consists of a single species of charged particles immersed in a uniform neutralizing background. The particles interact via a l/r potential, where r is the two dimensional separation. Equations of state were calculated for both the liquid and solid phases. Results of calculation of the thermodynamic functions and one and two particle correlation functions are presented.

Gann, R. C.↗
DOE OSTI2021

Geometric electrostatic particle-in-cell algorithm on unstructured meshes

We present a geometric particle-in-cell (PIC) algorithm on unstructured meshes for studying electrostatic perturbations with frequency lower than electron gyrofrequency in magnetized plasmas. In this method, ions are treated as fully kinetic particles and electrons are described by the adiabatic response. The PIC method is derived from a discrete variational principle on unstructured meshes. To preserve the geometric structure of the system, the discrete variational principle requires that the electric field is interpolated using Whitney 1-forms, the charge is deposited using Whitney 0-forms and the electric field is computed by discrete exterior calculus. The algorithm has been applied to study the ion Bernstein wave (IBW) in two-dimensional magnetized plasmas. The simulated dispersion relations of the IBW in a rectangular region agree well with theoretical results. In a two-dimensional circular region with fixed boundary condition, the spectrum and eigenmode structures of the IBW are obtained from simulations. We compare the energy conservation property of the geometric PIC algorithm derived from the discrete variational principle with that of previous PIC methods on unstructured meshes. The comparison shows that the new PIC algorithm significantly improves the energy conservation property.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Simulation of electron Bernstein waves using FullWave with a 2D non-local hot plasma model

Hot plasma wave simulation capability is expanded in the FullWave code by updating the hybrid iterative solver in the code with a semi-implicit time stepping method. The new approach is used to simulate Electron Bernstein Wave (EBW) heating in over-dense spherical tokamak plasmas. The code’s hybrid iterative solver circumvents the prohibitive memory cost of direct methods by combining a time evolution of Maxwell’s equations with frequency-domain relaxation, while the conductivity kernel, calculated via 3D particle tracking, captures the essential non-local wave–particle interactions. One-dimensional EBW simulations verify the algorithm’s accuracy by demonstrating mode conversion from X-mode wave to EBW at the upper hybrid resonance and a strong cyclotron damping near the plasma core. Two-dimensional simulation reproduces the predicted short EBW wavelength and quantitatively matches the hot-plasma dispersion relation. This study demonstrates the fidelity of the hybrid solver for the electron cyclotron frequency range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Impact of collisionality on turbulence in the edge of tokamak plasma using 3D global simulations

Collisionality is one of the key parameters in determining turbulent transport in the plasma edge, regulating phenomena such as 'shoulder formation', separation of scale lengths in the scrape-off layer (SOL), turbulence damping and zonal flow dynamics. Understanding its role is therefore of primary importance for future reactors like ITER. Getting reliable predictions and a better characterization of plasma flow properties when varying collisionality remains, however, a critical challenge for the simulations. This paper focuses on the impact of varying collisionality in a non-isothermal three-dimensional fluid model of the plasma edge. A high field side limited configuration encompassing open and closed magnetic field lines with parameters typical of a medium-sized tokamak is considered. The present model can consistently account for the variations of collisionality and its impact on both the parallel resistivity η ∥ and the ion and electron parallel thermal conductivities χ ∥e,i . Details on mean flow and turbulence properties are given. Changing collisionality leads to significant changes in the flow properties both on the mean and fluctuating quantities. In particular, lowering collisionality decreases the size of coherent structures, the fluctuation levels of turbulence, and steepens the density and temperature equilibrium profiles around the separatrix leading to a global reduction of the turbulent transport. The SOL width is observed to increase with collisionality, eventually resulting in the disappearance of the scale lengths separation between near and far SOL, consistently with previous experimental observations. At low collisionality, where the presence of narrow feature is well-established, a contribution of heat conduction increases up to compete with heat convection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Nonlinear coupling of whistler waves to oblique electrostatic turbulence enabled by cold plasma

Kinetic simulations and theory demonstrate that whistler waves can excite oblique, short-wavelength fluctuations through secondary drift instabilities if a population of sufficiently cold plasma is present. The excited modes lead to heating of the cold populations and damping of the primary whistler waves. The instability threshold depends on the density and temperature of the cold population and can be relatively small if the temperature of the cold population is sufficiently low. Furthermore, this mechanism may thus play a significant role in controlling amplitude of whistlers in the regions of the Earth's magnetosphere where cold background plasma of sufficient density is present.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Modulations of the plasma uniformity by low frequency sources in a large-area dual frequency inductively coupled plasma based on fluid simulations

As the wafer size increases, dual frequency (DF) inductively coupled plasma (ICP) sources have been proposed as an effective method to achieve large-area uniform plasma processing. A two-dimensional (2D) self-consistent fluid model, combined with an electromagnetic module, has been employed to investigate the influence of the low frequency (LF) source on the plasma radial uniformity in an argon DF discharge. When the DF antenna current is fixed at 10 A, the bulk plasma density decreases significantly with the LF due to the less efficient heating, and the best radial uniformity is obtained at 3.39 MHz. As the LF decreases to 2.26 MHz, the plasma density is characterized by an edge-high profile, and meanwhile the maximum of the electron temperature appears below the outer two-turn coil. Moreover, the axial ion flux at 3.39 MHz is rather uniform in the center region except at the radial edge of the substrate, where a higher ion flux is observed. When the inner five-turn coil frequency is fixed at 2.26 MHz, the plasma density profiles shift from edge-high over uniform to center-high as the LF coil current increases from 6 A to 18 A, and the best plasma uniformity is obtained at 14 A. In addition, the maximum of the electron temperature becomes lower with a second peak appears at the radial position of r = 9 cm at 18 A.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A parallel-kinetic-perpendicular-moment model for magnetised plasmas

We describe a new model for the study of weakly collisional, magnetised plasmas derived from exploiting the separation of the dynamics parallel and perpendicular to the magnetic field. This unique system of equations retains the particle dynamics parallel to the magnetic field while approximating the perpendicular dynamics through a spectral expansion in the perpendicular degrees of freedom, analogous to moment-based fluid approaches. In so doing, a hybrid approach is obtained that is computationally efficient enough to allow for larger-scale modelling of plasma systems while eliminating a source of difficulty in deriving fluid equations applicable to magnetised plasmas. We connect this system of equations to historical asymptotic models and discuss advantages and disadvantages of this approach, including the extension of this parallel-kinetic-perpendicular moment beyond the typical region of validity of these more traditional asymptotic models. This paper forms the first of a multi-part series on this new model, covering the theory and derivation, alongside demonstration benchmarks of this approach that include shocks and magnetic reconnection.

astrophysical plasmas↗
DOE OSTI2026

Role of ion acoustic instability in magnetic reconnection

We report on a first-principles numerical study of magnetic reconnection in plasmas with different initial ion-to-electron temperature ratios. In cases where this ratio is significantly below unity, we observe intense wave activity in the diffusion region, driven by the ion-acoustic instability. Our analysis shows that the dominant macroscopic effect of this instability is to drive substantial ion heating. In contrast to earlier studies reporting significant anomalous resistivity, we find that anomalous contributions due to the ion-acoustic instability are minimal. These results shed light on the dynamical impact of this instability on reconnection processes, offering new insights into the fundamental physics governing collisionless reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗