Engineering Papers⌕ Search

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.

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

At least 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Measurements and scaling of ion propulsion impulse during driven magnetic reconnection

Impulse scaling during magnetic reconnection, the magnetic energy conversion to kinetic energy, via direct Mach probe measurements in Magnetic Reconnection Experiment is examined. Ion exhaust velocity and impulse scalings with reconnecting magnetic field during the push phase of driven reconnection are presented. The outflows and impulse measurements are compared with global MHD simulations. Both measurements and simulations reveal a favorable scaling, greater than linear, of impulse with reconnecting field. These scaling results establish that magnetic reconnection could be utilized for plasma propulsion.

Ebrahimi, Fatima [Princeton Plasma Physics Laborat↗
DOE OSTI2024

Performance prediction applying different reduced turbulence models to the SMART tokamak

The SMall Aspect Ratio Tokamak (SMART) is currently being commissioned at the University of Seville and will be able to compare the performance of positive and negative triangularity plasmas at low aspect ratio. Predictive simulations have been performed for different machine scenarios and heating schemes using the TRANSP code. The objectives of these simulations are to predict the parameters expected in positive triangularity plasmas, to guide diagnostic development, and to validate transport models. Several reduced turbulence models have been used to predict electron and ion temperatures for the operational phase 2. All models provide similar results from approximately mid-radius to the separatrix but important discrepancies are found in the core region. These positive triangularity results are compared with experiments from a similar size machine like GLOBUS-M2. The multi-mode model (MMM) shows the best agreement. Simulations with different boundary conditions have been performed and no strong differences have been observed between them. The impact of neutral beam injection (NBI) on the predicted profiles has also been addressed. Rotation reduces turbulence levels so higher temperatures are achieved when included in the simulations. Studying the different contributions to the thermal diffusivities, it is observed that electron temperature gradient (ETG) turbulence dominates at the plasma core while micro-tearing modes (MTM) dominate at the edge in the electron channel. In the ion channel, the neoclassical contribution is dominant at the core and at the very edge while the Weiland component, which includes ion temperature gradient mode (ITG), trapped electron mode (TEM), kinetic ballooning mode (KBM), peeling mode (PM) and collisionless and collision dominated magnetohydrodynamic (MHD) modes governs the mid-radius region. For phase 3, two plasmas with different electron densities have been studied. The case with lower density matches well a specific discharge of GLOBUS-M2. The higher density plasma shows high performance with β N ≈ 3.8.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Pre-equilibrium photons from the early stages of heavy-ion collisions

Abstract We use QCD kinetic theory to compute photon production in the chemically equilibrating Quark-Gluon Plasma created in the early stages of high-energy heavy-ion collisions. We do a detailed comparison of pre-equilibrium photon rates to the thermal photon production. We show that the photon spectrum radiated from a hydrodynamic attractor evolution satisfies a simple scaling form in terms of the specific shear viscosityη/sand entropy densitydS/dζ∼ (Tτ 1/3 ) 3/2 ∞ . We confirm the analytical predictions with numerical kinetic theory simulations. We use the extracted scaling function to compute the pre-equilibrium photon contribution in$$ \sqrt{s_{NN}} $$ s NN = 2.76 TeV 0–20% PbPb collisions. We demonstrate that our matching procedure allows for a smooth switching from pre-equilibrium kinetic to thermal hydrodynamic photon production. Finally, our publicly available implementation can be straightforwardly added to existing heavy ion models.

Physics↗
DOE OSTI2025

WEST L-mode record long pulses guided by predictions using Integrated Modeling

A new record was set on the WEST Tokamak, designed to operate long duration plasmas in a tungsten (W) environment, with an injected energy of 1.15 GJ and a plasma duration 364s. Scenario development was supported by integrated modeling using the High Fidelity Plasma Simulator (HFPS), the European IMAS-coupled version of JETTO/JINTRAC, which integrates physics-driven modules into a unified framework. In particular, a reduced model for Lower-Hybrid heating and Current-Drive (LHCD) and the quasi-linear turbulent transport model TGLF are crucial for long pulses predictions up to the Last Closed Flux Surface (LCFS). Using this workflow, a 100 s reference discharge was modeled and plasma kinetic profiles and loop voltage were quantitatively well matched. In preparation for the recent long duration experiments, non-inductive current-drive actuators (I P , n e , P LHCD ) were varied to determine the operational domain going towards fully non-inductive discharges. In particular, decreasing the plasma current is shown to ease the access to such conditions, with a careful monitoring of (n e , P LHCD ) to avoid machine limitations. In addition, post-prediction experiments conducted within the investigated parameter range validated the predicted dependencies and were shown to be in quantitative agreement. Exploratory work on the use of ECCD for MHD stability purpose is also introduced.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Fermi-type particle acceleration from magnetic reconnection at the termination shock of a relativistic striped wind

An oblique-rotating pulsar generates a relativistic striped wind in a pulsar wind nebula (PWN). The termination shock of the PWN compresses the Poynting-flux-dominated flow and drives magnetic reconnection. By carrying out particle-in-cell simulations of the termination shock of the PWN, we study the shock structure as well as the energy conversion processes and particle acceleration mechanisms. With the recent advances in the numerical methods, we extend the simulations to the ultrarelativistic regime with a bulk Lorentz factor of up to γ 0 = 10 6 . Magnetic reconnection at the termination shock is highly efficient at converting magnetic energy to particle kinetic energy and accelerating particles to high energies. Similar to earlier studies, we find that the resulting energy spectra crucially depend on λ/d e (λ is the wavelength of the striped wind and de is the relativistic plasma skin depth). When λ/d e is large (λ ≳ 40d e ), the downstream particle spectra form a power-law distribution in the magnetically dominated relativistic wind regime. By analyzing particle trajectories and statistical quantities relevant to particle energization, we find that Fermi-type mechanism dominates the particle acceleration and power-law formation. We find that the results for particle acceleration are scalable as γ 0 and σ 0 increase to large values. The maximum energy for electrons and positrons can reach hundreds of TeV if the wind has a bulk Lorentz factor of γ0 ≈ 106 and a magnetization parameter of σ 0 = 10, which can explain the recent observations of high-energy gamma rays from PWNe

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2014

Bidirectional Energy Cascades and the Origin of Kinetic Alfvenic and Whistler Turbulence in the Solar Wind

The observed steep kinetic scale turbulence spectrum in the solar wind raises the question of how that turbulence originates. Observations of keV energetic electrons during solar quiet time suggest them as a possible source of free energy to drive kinetic turbulence. Using particle-in-cell simulations, we explore how the free energy released by an electron two-stream instability drives Weibel-like electromagnetic waves that excite wave-wave interactions. Consequently, both kinetic Alfvénic and whistler turbulence are excited that evolve through inverse and forward magnetic energy cascades.

plasmas↗
DOE OSTI2024

Statistical data analysis of x-ray spectroscopy data enabled by neural network accelerated Bayesian inference

Bayesian inference applied to x-ray spectroscopy data analysis enables uncertainty quantification necessary to rigorously test theoretical models. However, when comparing to data, detailed atomic physics and radiation transfer calculations of x-ray emission from non-uniform plasma conditions are typically too slow to be performed in line with statistical sampling methods, such as Markov Chain Monte Carlo sampling. Furthermore, differences in transition energies and x-ray opacities often make direct comparisons between simulated and measured spectra unreliable. Here, we present a spectral decomposition method that allows for corrections to line positions and bound–bound opacities to best fit experimental data, with the goal of providing quantitative feedback to improve the underlying theoretical models and guide future experiments. In this work, we use a neural network (NN) surrogate model to replace spectral calculations of isobaric hot-spots created in Kr-doped implosions at the National Ignition Facility. The NN was trained on calculations of x-ray spectra using an isobaric hot-spot model post-processed with Cretin, a multi-species atomic kinetics and radiation code. The speedup provided by the NN model to generate x-ray emission spectra enables statistical analysis of parameterized models with sufficient detail to accurately represent the physical system and extract the plasma parameters of interest.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

The role of electron-scale turbulence in the JET tokamak: experiments and modelling

Dedicated electron heat transport experiments have been carried out in L- and H-mode Deuterium plasmas of the JET-ILW tokamak to identify the amount of electron heat carried by electron-scale electron temperature gradient (ETG) modes. Ion cyclotron resonance heating at different positions has been used to probe the response of the electron temperature inverse gradient length R/L Te to changes in electron heat flux qe, while different amounts of neutral beam heating allowed to scan the ratio of ion to electron temperature T e /T i , which is a key parameter for the onset of ETGs. Results indicate a steepening of the normalized qe vs R/L Te curve above R/LT e ~ 8 for T e /T i ≤ 1, suggestive of the ETG onset. Ion-scale gyro-kinetic (GK) simulations match the ion heat flux and the low-R/L Te part of the qe curve, but do not reproduce such steepening at high R/L Te . Multi-scale GK simulations covering both ion and electron scales and including one impurity bundling light and heavy species indicate an ETG contribution only for R/L Te values larger than the experimental ones. Sensitivity studies of such result are difficult to achieve due to limitation in numerical resources. The quasi-linear TGLF model has been used for sensitivity studies. With the same bundled impurity as the GK multi-scale, TGLF shows the qe steepening at much larger R/L Te values than in experiment, but when using the real mix of light impurities neglecting the heavy impurities, TGLF gets closer to the experimental results. Profile simulations with TGLF including both light and heavy impurities show over-prediction of T e profiles and in some cases also of density, but good T i predictions, confirming issues with the model electron stiffness for these plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

On the local measurement of electric currents and magnetic fields using Thomson scattering in Weibel-unstable plasmas

Here, we demonstrate the capability of the Thomson Scattering (TS) diagnostic to measure locally the microscopic electron and ion currents in counter-streaming plasmas unstable to the Weibel or current-filamentation instability. Synthetic TS spectra are calculated with particle distribution functions obtained from particle-in-cell simulations and used to accurately reproduce the simulated currents. We show that this technique allows accurate local measurements of the magnetic field, thus opening the way for the complete experimental characterization of the growth rate, saturation, and nonlinear dynamics of electromagnetic instabilities in plasmas. We illustrate the application of this diagnostic to experimental TS data, which yields local measurements of the magnetic field in Weibel-unstable plasmas and indicates that the magnetic energy density reaches ~1% of the kinetic energy density of the flows, in agreement with previous numerical studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1981

Effects of auroral-particle anisotropies and mirror forces on high-latitude electric fields

It is noted that, for most of the mechanisms for the strong electric fields that characterize the narrow regions in which there is acceleration and precipitation of ring current and/or plasma-sheet plasma, certain effects must be taken into account in simulations of auroral electric fields. The effects are those of auroral particle anisotropy, of mirror forces due to the inhomogeneous geomagnetic field, of auroral electron backscatter by the atmosphere, and of electron trapping by the combination of magnetic mirroring and electrostatic forces. What is more, the effects of the very strong perpendicular electric field must also be taken into account in a kinetic description of the Poisson equation in order to achieve a unified theory of the auroral electrostatic structure. Progress in these areas during the past few years is reviewed. It is shown that particle anisotropies and mirror forces can account for some basic electrostatic features of the quiet arc, while additional effects may be occurring in strong events in which the parallel potential drop is more than about 10 kV.

Chiu, Y. T.↗
DOE OSTI2023

Probing strong-field QED in beam-plasma collisions

Ongoing progress in laser and accelerator technology opens new possibilities in high-field science, notably to investigate the largely unexplored strong-field quantum electrodynamics (SFQED) regime where electron-positron pairs can be created directly from light-matter or even light-vacuum interactions. Laserless strategies such as beam-beam collisions have also been proposed to access the nonperturbative limit of SFQED. Here we report on a concept to probe SFQED by harnessing the interaction between a high-charge, ultrarelativistic electron beam and a solid conducting target. When impinging onto the target surface, the beam self fields are reflected, partly or fully, depending on the beam shape; in the rest frame of the beam electrons, these fields can exceed the Schwinger field, thus triggering SFQED effects such as quantum nonlinear inverse Compton scattering and nonlinear Breit-Wheeler electron-positron pair creation. Through reduced modeling and kinetic numerical simulations, we show that this single-beam setup can achieve interaction conditions similar to those envisioned in beam-beam collisions, but in a simpler and more controllable way owing to the automatic overlap of the beam and driving fields. This scheme thus eases the way to precision studies of SFQED and is also a promising milestone towards laserless studies of nonperturbative SFQED.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Molecular dynamics simulation of hydrodynamic transport coefficients in plasmas

Molecular dynamics simulations are used to calculate transport coefficients in a two-component plasma interacting through a repulsive Coulomb potential. The thermal conductivity, electrical conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity are computed using the Green–Kubo formalism over a broad range of Coulomb coupling strength, 0.01 ≤ Γ ≤ 140. Emphasis is placed on testing standard results of the Chapman–Enskog solution in the weakly coupled regime (Γ ≪ 1) using these first-principles simulations. As expected, the results show good agreement for Γ ≲ 0.1. However, this agreement is only possible if careful attention is paid to the definitions of linear constitutive relations in each of the theoretical models, a point that is often overlooked. For example, the standard Green–Kubo expression for thermal conductivity is a linear combination of thermal conductivity, electrothermal, and thermoelectric coefficients computed in the Chapman–Enskog formalism. Meaningful results for electrical conductivity are obtained over the full range of coupling strengths explored, but it is shown that potential and virial components of the other transport coefficients diverge in the strongly coupled regime (Γ ≫ 1). In this regime, only the kinetic components of the transport coefficients are meaningful for a classical plasma.

Electrical conductivity↗
NASA NTRS1997

Physics of Boundaries and their Interactions in Space Plasmas

This report describes the work done by SciberNet, Inc. during the month of October. We are working on the further refinement of the model used in our large-scale hybrid simulations of the magnetopause. Specifically, we are experimenting with several ways of modeling the effects of cold magnetospheric ions into our simulations. In addition, we are preparing two presentations for the upcoming Fall AGU highlighting the results of these simulations. We have also made progress in our development of a new kinetic linear code which we are using to study the linear properties of the Kelvin-Helmholtz instability at the magnetopause. We have extended the code from the electrostatic limit to the fully electromagnetic regime and are currently in the process of debugging and testing the code. Finally, we have made several test runs with our 2-D hybrid code for the magnetopause. The inflow-outflow boundary conditions are working properly. However, there are issues related to the setup and evolution of the original equilibrium that we are still trying to resolve. Finally, we are preparing several presentations for the upcoming Fall AGU.

Omidi, Nojan↗
DOE OSTI2024

Complete quasilinear model for the acceleration-driven lower hybrid drift instability and a computational assessment of its validity

A complete quasilinear model is derived for the electrostatic acceleration-driven lower hybrid drift instability in a uniform two-species low-beta plasma in which current is perpendicular to the background magnetic field. The model consists of coupled nonlinear velocity space diffusion equations for the volume-averaged ion and electron distribution functions. Each species' diffusion coefficient depends on a time-evolving spectral density of the electric-field energy per unit volume and a time-evolving dispersion relation. The dispersion relation is expressed analytically in integral form without the use of asymptotic limits and applies to arbitrary distribution functions, so long as they can be expressed as a function of one velocity coordinate, e.g., f⁡(vy) or f⁡(v⊥). The quasilinear model conserves energy and is complete in that it fully describes the evolution of the distribution functions, including resonant and nonresonant particle-wave interactions, while accounting for distribution-function-dependent mixed-complex frequencies. Further, the quasilinear diffusion model is solved numerically and self-consistently using a Crank-Nicolson temporal discretization and a second-order finite-volume velocity-space discretization. Numerical solutions are compared to nonlinear fourth-order accurate continuum kinetic Vlasov-Poisson simulations. Evolution of electric-field energy, growth rates, distribution functions, and diffusion coefficients are shown to be in agreement with Vlasov simulations. The quasilinear model is shown to predict anomalous transport terms, like resistivity and heating, to within a factor of order unity. Discrepancies between the quasilinear model and Vlasov simulations are assessed and attributed primarily to lack of damping in the quasilinear description and to the use of unperturbed-orbit susceptibilities in the linear theory dispersion relation. The results illuminate the predictive accuracy of the quasilinear model, place approximate bounds on its validity, and provide much needed vetting of quasilinear theory's ability to predict the nonlinear state of a microturbulent plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Deterministic scattering of relativistic electron beams by off-resonant circularly polarized electromagnetic waves

A mechanism is presented whereby relativistic electron beams localized in phase space are deterministically scattered by coherent circularly polarized electromagnetic waves without stochastic processes. It is shown via an exact single-particle analysis that the condition for maximal scattering is an off-resonant condition, contrary to previous kinetic analyses that predict maximal diffusion or interaction at exact resonance or its harmonics. The mechanism, verified by single-particle simulations, enables a fast, nonlinear redistribution of the beam particles. A possible application of this mechanism to runaway electron suppression is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A self-consistent hybrid model of kinetic striations in low-current argon discharges

A self-consistent hybrid model of standing and moving striations was developed for low-current DC discharges in noble gases. We introduced the concept of surface diffusion in phase space ( r , u ) (where u denotes the electron kinetic energy) described by a tensor diffusion in the nonlocal Fokker–Planck kinetic equation for electrons in the collisional plasma. Electrons diffuse along surfaces of constant total energy ε = u - eφ ( r ) between energy jumps in inelastic collisions with atoms. Numerical solutions of the 1d1u kinetic equation for electrons were obtained by two methods and coupled to ion transport and Poisson solver. We studied the dynamics of striation formation in Townsend and glow discharges in argon gas at low discharge currents using a two-level excitation-ionization model and a ‘full-chemistry’ model, which includes stepwise and Penning ionization. Standing striations appeared in Townsend and glow discharges at low currents, and moving striations were obtained for the discharge currents exceeding a critical value. These waves originate at the anode and propagate towards the cathode. We have seen two types of moving striations with the two-level and full-chemistry models, which resemble the s and p striations previously observed in the experiments. Simulations indicate that processes in the anode region could control moving striations in the positive column plasma. The developed model helps clarify the nature of standing and moving striations in DC discharges of noble gases at low discharge currents and low gas pressures.

Physics↗
DOE OSTI2024

Pushing the Boundaries of Relativistic Fluid Dynamics: From Particle Colliders to Neutron Star Mergers

This project aimed at establishing connections between fluid dynamics, kinetic theory, and string theory to develop new tools to characterize the far-from-equilibrium properties of relativistic systems, which are instrumental in understanding the emergence of liquid-like behavior of the quark-gluon plasma formed in heavy-ion collisions. This project investigated the fundamental properties of viscous fluids in the full nonlinear regime to pave the way for the systematic inclusion of viscous effects in future numerical simulations of neutron star mergers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Plasma irregularities caused by cycloid bunching of the CRRES G-2 barium release

The Combined Release and Radiation Effects Satellite (CRRES) spacecraft carried a number of barium thermite canisters for release into the upper atmosphere. The barium release labeled G-2 showed evidence of curved irregularities not aligned with the ambient magnetic field B. The newly discovered curved structures can be explained by a process called cycloid bunching. Cycloid bunching occurs when plasma is created by photoionization of a neutral cloud injected at high velocity perpendicular to B. If the injection velocity is much larger than the expansion speed of the cloud, the ion trail will form a cycloid that has irregularities spaced by the product of the perpendicular injection speed and the ion gyroperiod, Images of the solar-illuminated barium ions are compared with the results of a three-dimensional kinetic simulation. Cycloid bunching is shown to be responsible for the rapid generation of both curved and field-aligned irregularities in the CRRES G-2 experiment.

Bernhardt, P. A.↗