Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinetic plasmas”

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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Demonstrations that the Solar Wind Is Not Accelerated by Waves

The present work uses both observations and theoretical considerations to show that hydromagnetic waves cannot produce the acceleration of the fast solar wind and the related heating of the open solar corona. Waves do exist, and can play a role in the differential heating and acceleration of minor ions, but their amplitudes are not sufficient to power the wind, as demonstrated by extrapolation of magnetic spectra from Helios and Ulysses observations. Dissipation mechanisms invoked to circumvent this conclusion cannot be effective for a variety of reasons. In particular, turbulence does not play a strong role in the corona as shown by both observations of coronal striations and theoretical considerations of line-tying to a nonturbulent photosphere, nonlocality of interactions, and the nature of the kinetic dissipation. In the absence of wave heating and acceleration, the chromosphere and transition region become the natural source of open coronal energization. We suggest a variant of the 'velocity filtration' approach in which the emergence and complex churning of the magnetic flux in the chromosphere and transition region continuously and ubiquitously produces the nonthermal distributions required. These particles are then released by magnetic carpet reconnection at a wide range of scales and produce the wind as described in kinetic approaches. Since the carpet reconnection is not the main source of the energization of the plasma, there is no expectation of an observable release of energy in nanoflares.

Roberts, Aaron↗
NASA NTRS2008

Observational and Theoretical Challenges to Wave or Turbulence Accelerations of the Fast Solar Wind

We use both observations and theoretical considerations to show that hydromagnetic waves or turbulence cannot produce the acceleration of the fast solar wind and the related heating of the open solar corona. Waves do exist as shown by Hinode and other observations, and can play a role in the differential heating and acceleration of minor ions but their amplitudes are not sufficient to power the wind, as demonstrated by extrapolation of magnetic spectra from Helios and Ulysses observations. Dissipation mechanisms invoked to circumvent this conclusion cannot be effective for a variety of reasons. In particular, turbulence does not play a strong role in the corona as shown by both eclipse observations of coronal striations and theoretical considerations of line-tying to a nonturbulent photosphere, nonlocality of interactions, and the nature of kinetic dissipation. In the absence of wave heating and acceleration, the chromosphere and transition region become the natural source of open coronal energization. We suggest a variant of the velocity filtration approach in which the emergence and complex churning of the magnetic flux in the chromosphere and transition region continuously and ubiquitously produces the nonthermal distributions required. These particles are then released by magnetic carpet reconnection at a wide range of scales and produce the wind as described in kinetic approaches. Since the carpet reconnection is not the main source of the energization of the plasma, there is no expectation of an observable release of energy in nanoflares.

Roberts, D. Aaron↗
NASA NTRS2014

Si-Ge Nano-Structured with Tungsten Silicide Inclusions

Traditional silicon germanium high temperature thermoelectrics have potential for improvements in figure of merit via nano-structuring with a silicide phase. A second phase of nano-sized silicides can theoretically reduce the lattice component of thermal conductivity without significantly reducing the electrical conductivity. However, experimentally achieving such improvements in line with the theory is complicated by factors such as control of silicide size during sintering, dopant segregation, matrix homogeneity, and sintering kinetics. Samples are prepared using powder metallurgy techniques; including mechanochemical alloying via ball milling and spark plasma sintering for densification. In addition to microstructural development, thermal stability of thermoelectric transport properties are reported, as well as couple and device level characterization.

thermoelectricity↗
DOE OSTI2020

Kinetic beaming in radiative relativistic magnetic reconnection: a mechanism for rapid gamma-ray flares in jets

ABSTRACT Rapid gamma-ray flares pose an astrophysical puzzle, requiring mechanisms both to accelerate energetic particles and to produce fast observed variability. These dual requirements may be satisfied by collisionless relativistic magnetic reconnection. On the one hand, relativistic reconnection can energize gamma-ray emitting electrons. On the other hand, as previous kinetic simulations have shown, the reconnection acceleration mechanism preferentially focuses high energy particles – and their emitted photons – into beams, which may create rapid blips in flux as they cross a telescope’s line of sight. Using a series of 2D pair-plasma particle-in-cell simulations, we explicitly demonstrate the critical role played by radiative (specifically inverse Compton) cooling in mediating the observable signatures of this ‘kinetic beaming’ effect. Only in our efficiently cooled simulations do we measure kinetic beaming beyond one light crossing time of the reconnection layer. We find a correlation between the cooling strength and the photon energy range across which persistent kinetic beaming occurs: stronger cooling coincides with a wider range of beamed photon energies. We also apply our results to rapid gamma-ray flares in flat-spectrum radio quasars, suggesting that a paradigm of radiatively efficient kinetic beaming constrains relevant emission models. In particular, beaming-produced variability may be more easily realized in two-zone (e.g. spine-sheath) set-ups, with Compton seed photons originating in the jet itself, rather than in one-zone external Compton scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A Generalized Grain-Scale Model for the Non-Plasma and Plasma-Assisted Hydrogen Direct Reduction of Iron Ore

Direct Reduction of Iron ore using hydrogen (H-DRI) is a promising pathway towards efficient steelmaking and accurate predictive models are a necessity for scale-up and optimization of this technology. However, accurate models of this process remain limited because existing models oversimplify grain-scale phenomena, such as nonlinearity inside grain, self-sufficient porosity, surface reactions, and the role of plasma species. These phenomena are important for flash steelmaking and plasma-assisted H-DRI processes. To address this need, we present a phenomenological model for simulating H-DRI at the scale of a single micron-sized grain of the iron ore. We call this the Transient Reactive Grain Model (TRGM). TRGM incorporates key physical process: gas species transport, a chemical kinetics of material conversion, nanopore structural evolution and, adsorption-desorption surface kinetics at the reactive nanopore surface. The important contribution of this work is that the model provides a dependence on different reductant species, specifically hydrogen atoms versus molecules, so that role of hydrogen plasma reduction can be clarified compared to the use of pure hydrogen gas reduction. TRGM predictions agree well with experimental data for both molecular H2 reduction of Fe2O3 and plasma hydrogen reduction of Fe3O4. Results reveal species concentration gradients with a diffuse reaction zone, and enhanced hydrogen diffusion at the grain outer surface due to evolving porosity. These findings challenge common assumptions in existing models, including sharp reaction fronts, quasi-steady diffusion and kinetics, and the neglect of surface chemistry. As a generalized grain-scale model for H-DRI processes, TRGM has practical applications in flash steelmaking and in-flight reduction using both molecular and plasma hydrogen.

08 HYDROGEN↗
DOE OSTI2021

Fermi acceleration in relativistic collisionless plasma shocks correlates with anisotropic energy gains

Collisionless shocks generated by two colliding relativistic electron-positron plasma shells are studied using particle-in-cell simulations. Shocks are mediated by the Weibel instability (WI), and the kinetic energy of the fastest accelerated particles is found to be anisotropically modified by WI-induced electric fields. Specifically, we show that all particles interacting with the shock bifurcate into two groups based on their final relativistic Lorentz factor γ: slow (γ < γ bf ) and fast (γ < γ bf ), where γ bf is the bifurcation Lorentz factor that was found to be approximately twice the initial (upstream) Lorentz factor γ0. We have found that the kinetic energies of the slow particles are equally affected by the longitudinal and transverse components of the shock electric field, whereas the fast particles are primarily accelerated by the transverse field component.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Plasma-Catalyst Reactivity Control of Surface Nitrogen Species through Plasma-Temperature-Programmed Hydrogenation to Ammonia

Nonthermal plasma activation of N 2 can facilitate nitrogen adsorption on metal catalysts at low bulk temperatures and atmospheric pressure. Here, we apply a plasma-assisted temperature-programmed reaction (plasma-TPRxn) for ammonia (NH 3 ) synthesis using sequential exposure of a silica-supported metal catalyst to N 2 plasma followed by thermal hydrogen treatment while ramping the temperature to decouple the plasma activation of N 2 from surface catalyzed hydrogenation steps. This approach eliminates the effects from bulk plasma phase reactions, thereby allowing for direct interrogation of plasma activated nitrogen on the active metal surfaces. We confirm previously reported spectroscopic observations that show plasma-generated surface nitrogen can be converted to NH 3 through surface catalyzed pathways. Further, we demonstrate that the ammonia desorption peak temperature is sensitive to metal, with Pt desorbing NH 3 at the lowest temperature. Unsteady state microkinetic models of desorption kinetics as a function of initial N coverage and metal recover observed trends in NH 3 desorption temperatures and confirm that observed results reflect hydrogenation of plasma-induced N accommodation at each surface. In total, we show that the hydrogenation ability of the catalyst after plasma activation of N 2 is responsible for the reactivity trends observed in plasma-assisted NH 3 synthesis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2022

Hybrid simulation of energetic particles interacting with magnetohydrodynamics using a slow manifold algorithm and GPU acceleration

The hybrid method combining particle-in-cell and magnetohydrodynamics can be used to study the interaction between energetic particles and global plasma modes. In this paper we introduce the M3D-C1-K code, which is developed based on the M3D-C1 finite element code solving the magnetohydrodynamics equations, with a newly developed kinetic module simulating energetic particles. The particle pushing is done using a new algorithm by applying the Boris pusher to the classical Pauli particles to simulate the slow-manifold of particle orbits, with long-term accuracy and fidelity. The particle pushing can be accelerated using GPUs with a significant speedup. The moments of the particles are calculated using the δƒ method, and are coupled into the magnetohydrodynamics simulation through pressure or current coupling schemes. Several linear simulations of magnetohydrodynamics modes driven by energetic particles have been conducted using M3D-C1-K with the δƒ method, including fishbone, toroidal Alfvén eigenmodes and reversed shear Alfvén eigenmodes. Good agreement with previous results from other eigenvalue, kinetic and hybrid codes have been achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Perpendicular magnetofluid theory for magnetically confined plasmas in the collisionless limit

Note that the magnetically confined plasmas in thermonuclear fusion turn to be collisionless. Some particles are trapped by the inhomogeneous equilibrium field or by the waves, while some others keep circulating. One cannot expect the charged particles to be collectively localized in space by collisions along the magnetic field lines. Consequently, the fluid parameters relating to the parallel motion, such as parallel fluid velocity, become inadequate. In this paper, the perpendic- ular magnetofluid theory for magnetically confined plasmas in the collisionless limit is developed. It revises Braginskii’s transport theory by releasing the collisional dominance assumption. In the perpendicular direction, the particle spatial localization is solo resulted from the strong magnetic field, while the finite Larmor radius (FLR) effects are taken into account in the next order. In the parallel direction, the particle mobility feature is fully retained by constructing only the per- pendicular moments. Therefore, besides the FLR effects till order O($k$ $^{2}_{⊥}$ $ρ$ $^{2}_{i}$ ), the Landau damping, the trapped particle effects, and other parallel kinetic effects are all kept in the current formalism. Here, k ⊥ is the perpendicular wave number and ρ i is the ion Larmor radius. Here, the theory can be used both for the analytical and numerical studies of magnetically confined plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2001

Time-dependent oral absorption models

The plasma concentration-time profiles following oral administration of drugs are often irregular and cannot be interpreted easily with conventional models based on first- or zero-order absorption kinetics and lag time. Six new models were developed using a time-dependent absorption rate coefficient, ka(t), wherein the time dependency was varied to account for the dynamic processes such as changes in fluid absorption or secretion, in absorption surface area, and in motility with time, in the gastrointestinal tract. In the present study, the plasma concentration profiles of propranolol obtained in human subjects following oral dosing were analyzed using the newly derived models based on mass balance and compared with the conventional models. Nonlinear regression analysis indicated that the conventional compartment model including lag time (CLAG model) could not predict the rapid initial increase in plasma concentration after dosing and the predicted Cmax values were much lower than that observed. On the other hand, all models with the time-dependent absorption rate coefficient, ka(t), were superior to the CLAG model in predicting plasma concentration profiles. Based on Akaike's Information Criterion (AIC), the fluid absorption model without lag time (FA model) exhibited the best overall fit to the data. The two-phase model including lag time, TPLAG model was also found to be a good model judging from the values of sum of squares. This model also described the irregular profiles of plasma concentration with time and frequently predicted Cmax values satisfactorily. A comparison of the absorption rate profiles also suggested that the TPLAG model is better at prediction of irregular absorption kinetics than the FA model. In conclusion, the incorporation of a time-dependent absorption rate coefficient ka(t) allows the prediction of nonlinear absorption characteristics in a more reliable manner.

Non-NASA Center↗
DOE OSTI2022

Simulation of runaway electron production with CQL3D coupled to NIMROD

Abstract A coupling between two distinctly different codes—one magnetohydrodynamic (MHD) and another kinetic—is achieved and applied for simulation of runaway electron (RE) production. The 3D initial value MHD code NIMROD simulates a DIII-D pure neon shattered pellet injection plasma quench including the propagation and ablation of the fragments, ionization and recombination of the impurities, and the radiated and transported energies. The field data from NIMROD is then used by the bounce-averaged Fokker–Planck Collisional QuasiLinear 3D (CQL3D) kinetic code to simulate the production of REs and their radial transport. The coupling procedure involves mapping of data between different grids and adjustment of the NIMROD toroidal electric field when REs appear. It is shown that without the radial transport, a large RE current is generated, up to 30% of the pre-pellet ohmic current. However, when the radial transport is included in CQL3D, the RE current is reduced to undetectable level, consistent with experiment. Various forms of the radial diffusion are surveyed to determine conditions when the fast electrons would not have time to be accelerated to relativistic energies before they are lost to chamber wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Time-resolved X-ray imaging of the current filamentation instability in solid-density plasmas

The streaming of energetic charged particles can magnetize astrophysical and laboratory plasmas via the current filamentation instability. Despite its importance, the experimental characterization of this instability has remained a challenge. Here, we report an experiment combining a high-intensity optical laser with a high-brightness X-ray free electron laser that successfully images the instability in solid-density plasmas with 200nm spatial and 50fs temporal resolution. We characterize the development of µm-scale filamentary structures and their evolution over tens of picoseconds through a non-linear merging process. The measured plasma density modulations and long merging time reveal the critical importance of space-charge effects and ion motion on this electron-driven instability. Supporting theoretical analysis and kinetic simulations help distinguish the relative role of space-charge and resistive effects. Our findings indicate that magnetic fields on the order of 10megagauss are produced, with important implications for transport and radiation emission of energetic particles in plasmas.

Schoenwaelder, Christopher↗
NASA NTRS1987

MPD thruster plasma instability studies

Results of an ongoing experimental/theoretical effort to understand the influence of plasma waves and instabilities upon coaxial plasma thruster performance and longevity are discussed. An experiment which permits measurement of the plasma dispersion relation plus a suite of plasma parameters for each one millisecond quasi-steady discharge is described. This investigation is guided by a linearized Boltzmann-Poisson kinetic numerical model which can accommodate the effects of collisions, inhomogeneities and the magnetic field. Preliminary tests reveal the presence of the low-frequency branch of an ion cyclotron wave. Although the stability of these waves still has to be established through spatial growth rate measurements, it is well known that they are prone to 'current-driven' instabilities. 'Current-driven' instabilities in the collisional tensor conducting plasma of the MPD thruster are analytically discussed and their relevance to the MPD plasma is illustrated by estimating a critical total current at which their onset occurs.

Kelly, A. J.↗
DOE OSTI2021

Computing the Paschen curve for argon with speed-limited particle-in-cell simulation

Upon inclusion of collisions, the speed-limited particle-in-cell (SLPIC) simulation method successfully computed the Paschen curve for argon. The 1D3V simulations modeled an electron cascade across an argon-filled capacitor, including electron-neutral ionization, electron-neutral elastic collisions, electron-neutral excitation, and ion-induced secondary electron emission. In electrical breakdown, the timescale difference between ion and electron motion makes traditional PIC methods computationally slow. To decrease this timescale difference and speed up computation, here we used SLPIC, a time-domain algorithm that limits the speed of the fastest electrons in the simulation. The SLPIC algorithm facilitates a straightforward, fully kinetic treatment of dynamics and collisions. SLPIC was as accurate as PIC, but ran up to 200 times faster. SLPIC accurately computed the Paschen curve for argon over three orders of magnitude in pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1987

Emerging flux, magnetic reconnection, plasma turbulence and waves in the transition zone of the solar atmosphere

Ultraviolet spectroscopy of the solar transition zone, especially Spacelab 2 results, is reviewed. Emerging magnetic flux is considered as the energy source of the solar corona and the solar wind. The conversion of magnetic into kinetic energy is facilitated in the transition zone because of unique conditions. Radiation losses are sufficient to create fast instabilities. Observed nonthermal velocities in spectra of the transition zone indicate the existence of strong plasma oscillations and turbulence. The small filling factor of the transition zone indicates a highly filamentary structure of density, magnetic field strength, and currents, which results in rapid reconnection time scales. Particle acceleration and the heating of the corona are placed in the areas of strong plasma oscillations of the transition zone.

Brueckner, G. E.↗
DOE OSTI2022

Global gyrokinetic simulations of ASDEX Upgrade up to the transport timescale with GENE–Tango

An accurate description of turbulence up to the transport timescale is essential for predicting core plasma profiles and enabling reliable calculations for designing advanced scenarios and future devices. Here, we exploit the gap separation between turbulence and transport timescales and couple the global gyrokinetic code GENE to the transport-solver Tango, including kinetic electrons, collisions, realistic geometries, toroidal rotation and electromagnetic effects for the first time. This approach overcomes gyrokinetic codes' limitations and enables high-fidelity profile calculations in experimentally relevant plasma conditions, significantly reducing the computational cost. We present numerical results of GENE–Tango for two ASDEX Upgrade discharges, one of which exhibits a pronounced peaking of the ion temperature profile not reproduced by TGLF–ASTRA. We show that GENE–Tango can correctly capture the ion temperature peaking observed in the experiment. By retaining different physical effects in the GENE simulations, e.g., collisions, toroidal rotation and electromagnetic effects, we show that the ion temperature profile's peaking can be linked to electromagnetic effects of submarginal (stable) KBM modes. Furthermore, based on these results, the expected GENE–Tango speedup for the ITER standard scenario is larger than two orders of magnitude compared to a single gyrokinetic simulation up to the transport timescale, possibly making first-principles ITER simulations feasible on current computing resources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Effects of emitting surfaces and trapped ions on the sheath physics and current flow in multidimensional plasma systems

Recent one-dimensional simulations of planar sheaths with strong electron emission have shown that trapping of charge-exchange ions causes transitions from space-charge limited (SCL) to inverse sheaths. However, multidimensional emitting sheath phenomena with collisions remained unexplored, due in part to high computational cost. As such, we developed a novel continuum kinetic code to study the sheath physics, current flow and potential distributions in two-dimensional unmagnetized configurations with emitting surfaces. For small negatively biased thermionic cathodes in a plasma, the cathode sheath can exist in an equilibrium SCL state. The SCL sheath carries an immense density of trapped ions, neutralized by thermoelectrons, within the potential well of the virtual cathode. For further increases of emitted flux, the trapped ion cloud expands in space. The trapped ion space charge causes an increase of thermionic current far beyond the saturation limit predicted by conventional collisionless SCL sheath models without ion trapping. For sufficiently strong emission, the trapped ion cloud consumes the entire 2D plasma domain, forming a mode with globally confined ions and an inverse sheath at the cathode. In situations where the emitted flux is fixed and the bias is swept (e.g. emissive probe), the trapped ions cause a large thermionic current to escape for all biases below the plasma potential. Strong suppression of the thermionic emission, required for the probe to float, only occurs when the probe is above the plasma potential.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗