SEARCH · Engineering Papers

Results for “Plasma kinetic 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 541 records · Page 30

All sources NASA NTRS DOE OSTI

NASA NTRS2017

Kinetic Scale Structure of Low-Frequency Waves and Fluctuations

The dissipation of solar wind turbulence at kinetic scales is believed to be important for the heating of the corona and for accelerating the wind. The linear Vlasov kinetic theory is a useful tool for identifying various wave modes, including kinetic Alfven, fast magnetosonic/whistler, and ion-acoustic (or kinetic slow), and their possible roles in the dissipation. However, the kinetic mode structure in the vicinity of ion-cyclotron modes is not clearly understood. The present paper aims to further elucidate the structure of these low-frequency waves by introducing discrete particle effects through hybrid simulations and Klimontovich formalism of spontaneous emission theory. The theory and simulation of spontaneously emitted low-frequency fluctuations are employed to identify and distinguish the detailed mode structures associated with ion-Bernstein modes versus quasi-modes. The spontaneous emission theory and simulation also confirm the findings of the Vlasov theory in that the kinetic Alfven waves can be defined over a wide range of frequencies, including the proton cyclotron frequency and its harmonics, especially for high-beta plasmas. This implies that these low-frequency modes may play predominant roles even in the fully kinetic description of kinetic scale turbulence and dissipation despite the fact that cyclotron harmonic and Bernstein modes may also play important roles in wave-particle interactions.

waves↗

DOE OSTI2020

Physics-guided machine learning approaches to predict the ideal stability properties of fusion plasmas

One of the biggest challenges to achieve the goal of producing fusion energy in tokamak devices is the necessity of avoiding disruptions of the plasma current due to instabilities. The disruption event characterization and forecasting (DECAF) framework has been developed in this purpose, integrating physics models of many causal events that can lead to a disruption. Two different machine learning approaches are proposed to improve the ideal magnetohydrodynamic (MHD) no-wall limit component of the kinetic stability model included in DECAF. First, a random forest regressor (RFR), was adopted to reproduce the DCON computed change in plasma potential energy without wall effects for a large database of equilibria from the national spherical torus experiment (NSTX). This tree-based method provides an analysis of the importance of each input feature, giving an insight into the underlying physics phenomena. Secondly, a fully-connected neural network has been trained on sets of calculations with the DCON code, to get an improved closed form equation of the no-wall β limit as a function of the relevant plasma parameters indicated by the RFR. The neural network has been guided by physics theory of ideal MHD in its extension outside the domain of the NSTX experimental data. The estimated value has been incorporated into the DECAF kinetic stability model and tested against a set of experimentally stable and unstable discharges. Moreover, the neural network results were used to simulate a real-time stability assessment using only quantities available in real-time. Finally, the portability of the model was investigated, showing encouraging results by testing the NSTX-trained algorithm on the mega ampere spherical tokamak (MAST).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Drift-ordered fluid vorticity equation with energy consistency

Although drift-ordered fluid models are widely applied in tokamak edge turbulence simulations, the models used are acknowledged not to conserve energy or even electrical charge. The present paper aims to remove many of the existing pitfalls in drift-fluid models, however, with the objective of finding a solution simple enough to be implemented in numerical applications. Our main result is an improved version of the drift-Braginskii equations involving a generalized vorticity function. In the new drift-Braginskii system, the quasi-neutrality condition translates into a transport equation for a generalized vorticity, expressed in conservation form, and related to the total mass-weighted circulation. It is found that kinetic energy conservation can be achieved if the polarization flow is defined recursively. In conclusion, the resulting model conserves the kinetic energy associated with E×B and diamagnetic flows and retains the associated perpendicular kinetic energy flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2021

Reconnection and particle acceleration in 3D current sheet evolution in moderately-magnetized astrophysical pair plasma

Magnetic reconnection, a plasma process converting magnetic energy to particle kinetic energy, is often invoked to explain magnetic energy releases powering high-energy flares in astrophysical sources including pulsar wind nebulae and black hole jets. Reconnection is usually seen as the (essentially 2D) nonlinear evolution of the tearing instability disrupting a thin current sheet. To test how this process operates in 3D, we conduct a comprehensive particle-in-cell simulation study comparing 2D and 3D evolution of long, thin current sheets in moderately-magnetized, collisionless, relativistically-hot electron- positron plasma, and find dramatic differences. We first systematically characterize this process in 2D, where classic, hierarchical plasmoid-chain reconnection determines energy release, and explore a wide range of initial configurations, guide magnetic field strengths, and system sizes. We then show that 3D simulations of similar configurations exhibit a diversity of behaviours, including some where energy release is determined by the nonlinear relativistic drift-kink instability. Thus, 3D current-sheet evolution is not always fundamentally classical reconnection with perturbing 3D effects, but, rather, a complex interplay of multiple linear and nonlinear instabilities whose relative importance depends sensitively on the ambient plasma, minor configuration details, and even stochastic events. It often yields slower but longer-lasting and ultimately greater magnetic energy release than in 2D. Intriguingly, nonthermal particle acceleration is astonishingly robust, depending on the upstream magnetization and guide field, but otherwise yielding similar particle energy spectra in 2D and 3D. Though the variety of underlying current-sheet behaviours is interesting, the similarities in overall energy release and particle spectra may be more remarkable.

Gregory R. Werner↗

DOE OSTI2020

Onset of fast magnetic reconnection and particle energization in laboratory and space plasmas

The onset of magnetic reconnection in space, astrophysical and laboratory plasmas is reviewed discussing results from theory, numerical simulations and observations. After a brief introduction on magnetic reconnection and approach to the question of onset, in this paper we first discuss recent theoretical models and numerical simulations, followed by observations of reconnection and its effects in space and astrophysical plasmas from satellites and ground-based detectors, as well as measurements of reconnection in laboratory plasma experiments. Mechanisms allowing reconnection spanning from collisional resistivity to kinetic effects as well as partial ionization are described, providing a description valid over a wide range of plasma parameters, and therefore applicable in principle to many different astrophysical and laboratory environments. Finally, we summarize the implications of reconnection onset physics for plasma dynamics throughout the Universe and illustrate how capturing the dynamics correctly is important to understanding particle acceleration. The goal of this review is to give a view on the present status of this topic and future interesting investigations, offering a unified approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

MARS-F/K modeling of plasma response and fast ion losses due to RMP in KSTAR

Abstract The toroidal single-fluid magnetohydrodynamic (MHD) code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) and the MHD-kinetic hybrid code MARS-K (Liu et al 2008 Phys. Plasmas 15 112503) are utilized to study the plasma response to the n = 1 ( n is the toroidal mode number) resonant magnetic perturbation (RMP), applied to suppress the type-I edge localized mode (ELM) in a KSTAR discharge. Both the resistive-rotating and ideal-static plasma models identify strong screening of the resonant radial field harmonics of the applied RMP due to the plasma response, and predict a strong edge-peeling response of the plasma which is consistent with the optimal ELM control coil current configuration adopted in experiment. The RMP-induced radial displacement of the plasma, computed by the resistive-rotating plasma model, agrees reasonably well with that reconstructed from the measured data in the plasma core. Taking into account the drift kinetic response of fast ions, MARS-K hybrid modeling also finds quantitative agreement of the plasma core fluid pressure perturbation with experiment. Based on the MARS-F computed plasma response, a guiding-center orbit-tracing simulation finds about 0.3% of fast ion losses due the n = 1 RMP in the KSTAR ELM control experiment considered. Most losses are associated with counter-current fast ions located near the plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Spatiotemporal Energy Cascade in Three-Dimensional Magnetohydrodynamic Turbulence

We present a new scale decomposition method to investigate turbulence in wavenumber-frequency space. Using 3D magnetohydrodynamic turbulence simulations, we show that magnetic fluctuations with time scales longer than the nonlinear time exhibit an inverse cascade toward even smaller frequencies. Low frequency magnetic fluctuations support turbulence, acting as an energy reservoir that is converted into plasma kinetic energy, the latter cascading toward large wavenumbers and frequencies, where it is dissipated. Our results shed new light on the spatiotemporal properties of turbulence, potentially explaining the origin and role of low frequency turbulent fluctuations in the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Fluctuation dynamo in a weakly collisional plasma

The turbulent amplification of cosmic magnetic fields depends upon the material properties of the host plasma. In many hot, dilute astrophysical systems, such as the intracluster medium (ICM) of galaxy clusters, the rarity of particle–particle collisions allows departures from local thermodynamic equilibrium. These departures – pressure anisotropies – exert anisotropic viscous stresses on the plasma motions that inhibit their ability to stretch magnetic-field lines. We present an extensive numerical study of the fluctuation dynamo in a weakly collisional plasma using magnetohydrodynamic (MHD) equations endowed with a field-parallel viscous (Braginskii) stress. When the stress is limited to values consistent with a pressure anisotropy regulated by firehose and mirror instabilities, the Braginskii-MHD dynamo largely resembles its MHD counterpart, particularly when the magnetic field is dynamically weak. If instead the parallel viscous stress is left unabated – a situation relevant to recent kinetic simulations of the fluctuation dynamo and, we argue, to the early stages of the dynamo in a magnetized ICM – the dynamo changes its character, amplifying the magnetic field while exhibiting many characteristics reminiscent of the saturated state of the large-Prandtl-number ( ${Pm}\gtrsim {1}$ ) MHD dynamo. We construct an analytic model for the Braginskii-MHD dynamo in this regime, which successfully matches simulated dynamo growth rates and magnetic-energy spectra. A prediction of this model, confirmed by our numerical simulations, is that a Braginskii-MHD plasma without pressure-anisotropy limiters will not support a dynamo if the ratio of perpendicular and parallel viscosities is too small. This ratio reflects the relative allowed rates of field-line stretching and mixing, the latter of which promotes resistive dissipation of the magnetic field. Finally, in all cases that do exhibit a viable dynamo, the generated magnetic field is organized into folds that persist into the saturated state and bias the chaotic flow to acquire a scale-dependent spectral anisotropy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Collisionless conduction in a high-beta plasma: a collision operator for whistler turbulence

The regulation of electron heat transport in high-β, weakly collisional, magnetized plasma is investigated. A temperature gradient oriented along a mean magnetic field can induce a kinetic heat-flux-driven whistler instability (HWI), which back-reacts on the transport by scattering electrons and impeding their flow. Previous analytical and numerical studies have shown that the heat flux for the saturated HWI scales as β$^{-1}_{e}$. These numerical studies, however, had limited scale separation and consequently large fluctuation amplitudes, which calls into question their relevance at astrophysical scales. To this end, we perform a series of particle-in-cell simulations of the HWI across a range of β$_e$ and temperature-gradient length scales under two different physical set-ups. The saturated heat flux in all of our simulations follows the expected β$^{-1}_{e}$ scaling, supporting the robustness of the result. We also use our simulation results to develop and implement several methods to construct an effective collision operator for whistler turbulence. The results point to an issue with the standard quasi-linear explanation of HWI saturation, which is analogous to the well-known 90° scattering problem in the cosmic-ray community. Despite this limitation, the methods developed here can serve as a blueprint for future work seeking to characterize the effective collisionality caused by kinetic instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

How transient effects can change the inferred plasma temperatures when using K-shell sulfur line ratios in a time dependent iron sulfide plasma

This paper describes how the steady state atomic kinetics approximation can underestimate the electron temperature determined from K-shell lines in ps-time-scale transient plasmas. In particular, we model the temperature determination of solid FeS targets used in opacity experiments at the Orion laser facility from the ratio of sulfur He-α to Ly-α lines. Such experiments use short-pulse lasers to heat a thin microdot of FeS buried in a plastic target to temperatures of more than 1 keV and densities of approximately 1–2 g/cm 3 . Using atomic kinetics calculations based on a temperature history from a radiation hydrodynamic simulation of the target evolution, the peak temperature inferred from the sulfur line ratios is 1.29 keV at 3.1 ps as compared with the input peak temperature of 1.41 keV at 2.0 ps. There is a time lag of 1.2 ps at the peak, and an overall 0.5 ps time lag in the temporal history of the temperature as the plasma cools over the next 10 ps.

Rilloraza, Francisco↗

DOE OSTI2024

Fast ion relaxation in ITER mediated by Alfvén instabilities

Abstract We address the critical issue for future burning plasmas of whether high-energy fusion products or auxiliary heating-beam ions will be confined for a sufficiently long time to compensate for thermal plasma energy losses. This issue can be mitigated by one of the most deleterious collective phenomena—the instability of low, sub-cyclotron frequency Alfvén eigenmodes (AEs), such as toroidicity-induced AEs and reversed-shear AEs in the ITER steady-state scenario. Using a revised quasi-linear (QL) theory applied to energetic particle (EP) relaxation in the presence of AEs, we find that the AE instabilities can affect both neutral beam ions and alpha particles, although the resulting fast ion transport is expected to be modest if classical particle slowing down is assumed. On the other hand, the QL theory predicts that the AE amplitudes will be enhanced by the background microturbulence, although this topic remains outside our scope due to the significant numerical effort required to evaluate these effects. We report our results for EP relaxation dynamics obtained utilizing several tools: (i) a comprehensive linear stability study of the sub-cyclotron Alfvénic spectrum as computed by ideal magnetohydrodynamic NOVA simulations for the AE eigenproblem, (ii) drift kinetic NOVA-C calculations for wave–particle interaction and AE growth/damping rates, and (iii) predictive QL modeling coupled with the global transport code TRANSP to assess the EP relaxation on the equilibrium timescale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2008

How Ionospheric Ions Populate the Magnetosphere during a Magnetic Storm

Ionospheric oxygen ions have been observed throughout the magnetosphere, from the plasma sheet to the ring current region. I t has been found that the O+ /H+ density ratio in the magnetosphere increases with geomagnetic activity and varies with storm phases. During the magnetic storm in late September to earIy October 2002, Cluster was orbiting in the plasma sheet and ring current regions. At prestorm time, Cluster observed high H+ density and low O+ density in the plasma sheet and lobes. During the storm main phase, 0+ density has increased by 10 times over the pre-storm level. Strong field-aligned beams of O+ were observed in the lobes. O+ fluxes were significantly reduced in the central plasma sheet during the storm recovery. However, 0+ was still evident on the boundaries of the plasma sheet and in the lobes. In order to interpret the Cluster observations and to understand how O+ ions populate the magnetosphere during a magnetic storm, we model the storm in early October 2002 using our global ion kinetic simulation (GIK). We use the LFN global simulation model to produce electric and magnetic fields in the outer magnetosphere, the Strangeway outflow scaling with Delcourt ion trajectories to include ionospheric outflows, and the Fok inner magnetospheric model for the plasmaspheric and ring current response to all particle populations. We find that the observed composition features are qualitatively reproduced by the simulations, with some quantitative differences that point to future improvements in the models.

Fok, Mei-Ching↗

DOE OSTI2025

Simulation Models for Exploring Magnetic Reconnection

Simulations have played a critical role in the advancement of our knowledge of magnetic reconnection. However, due to the inherently multiscale nature of reconnection, it is impossible to simulate all physics at all scales. For this reason, a wide range of simulation methods have been crafted to study particular aspects and consequences of magnetic reconnection. This article reviews many of these methods, laying out critical assumptions, numerical techniques, and giving examples of scientific results. Plasma models described include magnetohydrodynamics (MHD), Hall MHD, Hybrid, kinetic particle-in-cell (PIC), kinetic Vlasov, Fluid models with embedded PIC, Fluid models with direct feedback from energetic populations, and the Rice Convection Model (RCM).

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Similarity properties in capacitive radio frequency plasmas with nonlinear collision processes

Abstract Similarity laws (SL) are essential for correlating the characteristics of plasmas at different dimensional scales, which have been validated for radio frequency (rf) discharges at low pressure but under limited conditions. In this work, we evaluate the effects of nonlinear collisions (e.g. stepwise ionization) on the similarity properties in capacitive rf discharges in argon across a wide range of pressure regimes via fully kinetic particle-in-cell simulations. The SL scalings of fundamental discharge parameters, e.g. the electron power absorption and electron energy probability function, are examined with and without nonlinear collisions, respectively. Without the nonlinear collisions, the similarity scalings are found to be rigorously valid. When the nonlinear collisions are considered, the similarity properties in rf discharges still exist approximately, which indicates that the violations caused by the nonlinear collisions are still minor in the studied cases. The reasons for the effectiveness of SL scalings with nonlinear collision processes are also discussed. The results from this study confirmed the validity of similarity transformations with more complex reaction kinetics, which may promote the applicability of the SL scalings for the design and fabrications of plasma devices.

Physics↗

DOE OSTI2026

Mathematical Models and Numerical Methods for High-Fidelity Simulation of Ignition of Reactive Mixtures by Nanosecond Plasma Discharges in Realistic Configurations

We present a newly developed framework for the numerical simulation of ignition of reactive mixtures using single or repeated nanosecond discharge pulses. The framework builds upon the AMReX library, using the existing compressible solver PeleC and low-Mach solver PeleLMeX and allowing for adaptive mesh refinement, complex geometries, and execution on next-generation high-performance computing (HPC) systems. High-fidelity elementary models are adopted for weakly-ionised plasma discharges with significant energy deposition, consistent with nanosecond discharge pulses, and then implemented in the solver. The treatment of non-thermal electrons and charged species, thermodynamics of non-equilbrium species, plasma kinetics, limiting time scales, and boundary conditions for charged species are discussed and addressed for computational efficiency. The framework is demonstrated for three relevant applications: single and multi-pulse discharges in air, single pulse ignition of an ethylene/air mixture, and a three-dimensional plasma discharge in air with temperature stratification. The successful application of the framework demonstrates the feasibility of high-fidelity simulation of ignition of air/hydrocarbon mixtures in three-dimensions with multiple discharge pulses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Symmetric Set of Transport Coefficients for Collisional Magnetized Plasma

Braginskii extended magnetohydrodynamics is used to model transport in collisional astrophysical and high energy density plasmas. We show that commonly used approximations to the α ⊥ and β ⊥ transport coefficients [e.g., Epperlein and Haines, Phys. Fluids 29 , 1029 (1986) ] have a subtle inaccuracy that causes significant artificial magnetic dissipation and discontinuities. This is because magnetic transport actually relies on β ∥ - β ⊥ and α ⊥ - α ∥ , rather than α ⊥ and β ⊥ themselves. We provide fit functions that rectify this problem and thus resolve the discrepancies with kinetic simulations in the literature. Furthermore, then implemented in the g orgon code, they reduce the predicted density asymmetry amplitude at laser ablation fronts. Recognizing the importance of α ⊥ - α ∥ and β ∥ - β ⊥ , we recast the set of coefficients. This makes explicit the symmetry of the magnetic and thermal transport, as well as the symmetry of the coefficients themselves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Reconnection and particle acceleration in three-dimensional current sheet evolution in moderately magnetized astrophysical pair plasma

Magnetic reconnection, a plasma process converting magnetic energy to particle kinetic energy, is often invoked to explain magnetic energy releases powering high-energy flares in astrophysical sources including pulsar wind nebulae and black hole jets. Reconnection is usually seen as the (essentially two-dimensional) nonlinear evolution of the tearing instability disrupting a thin current sheet. To test how this process operates in three dimensions, we conduct a comprehensive particle-in-cell simulation study comparing two- and three-dimensional evolution of long, thin current sheets in moderately magnetized, collisionless, relativistically hot electron–positron plasma, and find dramatic differences. We first systematically characterize this process in two dimensions, where classic, hierarchical plasmoid-chain reconnection determines energy release, and explore a wide range of initial configurations, guide magnetic field strengths and system sizes. We then show that three-dimensional (3-D) simulations of similar configurations exhibit a diversity of behaviours, including some where energy release is determined by the nonlinear relativistic drift-kink instability. Thus, 3-D current sheet evolution is not always fundamentally classical reconnection with perturbing 3-D effects but, rather, a complex interplay of multiple linear and nonlinear instabilities whose relative importance depends sensitively on the ambient plasma, minor configuration details and even stochastic events. It often yields slower but longer-lasting and ultimately greater magnetic energy release than in two dimensions. Intriguingly, non-thermal particle acceleration is astonishingly robust, depending on the upstream magnetization and guide field, but otherwise yielding similar particle energy spectra in two and three dimensions. Although the variety of underlying current sheet behaviours is interesting, the similarities in overall energy release and particle spectra may be more remarkable.

Physics↗

NASA NTRS2005

Magnetohydrodynamic Power Generation in the Laboratory Simulated Martian Entry Plasma

This paper addresses the magnetohydrodynamic (MHD) conversion of the energy released during the planetary entry phase of an interplanetary vehicle trajectory. The effect of MHD conversion is multi-fold. It reduces and redirects heat transferred to the vehicle, and regenerates the dissipated energy in reusable and transportable form. A vehicle on an interplanetary mission carries about 10,000 kWh of kinetic energy per ton of its mass. This energy is dissipated into heat during the planetary atmospheric entry phase. For instance, the kinetic energy of Mars Pathfinder was about 4220 kWh. Based on the loss in velocity, Mars Pathfinder lost about 92.5% of that energy during the plasma-sustaining entry phase that is approximately 3900 kWh. An ideal MHD generator, distributed over the probe surface of Mars Pathfinder could convert more than 2000 kWh of this energy loss into electrical energy, which correspond to more than 50% of the kinetic energy loss. That means that the heat transferred to the probe surface can be reduced by at least 50% if the converted energy is adequately stored, or re-radiated, or directly used. Therefore, MHD conversion could act not only as the power generating, but also as the cooling process. In this paper we describe results of preliminary experiments with light and microwave emitters powered by model magnetohydrodynamic generators and discuss method for direct use of converted energy.

Vuskovic, L.↗