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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Hybrid Simulations of FRC Merging and Compression

An improved understanding of field-reversed configuration (FRC) merging and stability in high acceleration and compression magnetic fields is needed to speed up the development of the pulsed fusion concept developed at Helion Energy. All previous theoretical and simulation work on FRC merging and compression was performed using two-dimensional (2D) magnetohydrodynamic (MHD) models. The results of novel 2D hybrid simulations (fluid electrons and full-orbit kinetic ions) of FRC merging and compression are presented. Results of kinetic and MHD simulations, computed using the HYM code, are compared and analyzed. In cases without axial magnetic compression, both the MHD and hybrid simulations show a high sensitivity to the initial parameters (i.e. FRC separation, velocity, normalized separatrix radius, and plasma viscosity), showing that FRCs with large elongation and separatrix radius either do not merge or merge partially, forming a doublet FRC. In conclusion, application of a mirror coil field at the FRC ends with increasing strength is shown to lead to fast and complete merging of the FRCs in MHD and kinetic simulations.

FRC↗
DOE OSTI2021

Effect of hydrogen isotope on plasma impedance and thermal pinching induced by single aerosol droplets injected into an inductively coupled plasma

The mass of the hydrogen isotope in micrometer-sized aqueous droplets introduced into an ICP has been found to affect the degree of resulting plasma perturbation. Studied plasma perturbations include local plasma cooling, plasma shrinkage due to thermal pinching, and changes in plasma impedance. The local cooling effect caused by a vaporizing aerosol droplet was found to be similar between H 2 O and D 2 O. This finding is likely a result of the similar heat capacities and vaporization rates of H 2 O and D 2 O droplets and of the comparable thermal conductivities and dissociation kinetics of H 2 O and D 2 O vapor. In contrast, a clear isotope effect was observed for plasma shrinkage and impedance change, with H 2 O causing a considerably stronger effect. Plasma shrinkage, gauged by the reduction in Ar emission at the edge of the plasma, was detected at all observation heights within the load-coil region, with D 2 O droplets causing as much as 30% less of an effect than H 2 O. Likewise, D 2 O droplets induce significantly less change in plasma impedance than H 2 O droplets. Depending on experimental conditions and the probing method, changes in plasma impedance caused by monodisperse D 2 O droplets were only 52%–80% that of H 2 O. The lower influence on plasma shrinkage and impedance change for D 2 O droplets is believed to be related to the much lower thermal conductivity of atomic deuterium, around 70% that of atomic hydrogen (protium). The time required for the plasma impedance to return to its original steady-state value was comparatively long (~10 ms) and similar for H 2 O and D 2 O droplets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

W7-X and the sawtooth instability: towards realistic simulations of current-driven magnetic reconnection

Magnetic reconnection in W7-X is studied by means of global numerical simulations in a series of models of increasing complexity. The magnetic geometry ranges from that of a cylinder to the full three-dimensional field of W7-X, and the equations solved range from ideal magnetohydrodynamics (MHD) to gyrokinetics. We simulate plasmas from the first operation phase with electron cyclotron current drive (ECCD). These are characterized by an equilibrium magnetic field featuring an ECCD-distorted 'humped' profile of the rotational transform ι, with ι = 1 in two radial locations. Such plasmas generally show sawtooth activity, hence motivating the present study. We pay particular attention to the role of equilibrium current density gradients in the destabilization of reconnecting modes. When the equilibrium temperature and density gradients are artificially suppressed (to eliminate the pressure gradient drive), the perturbed electrostatic potential is radially localized between the locations at which ι = 1. This is shown with a purely collisionless gyrokinetic model, in cylindrical geometry. In the real toroidal geometry of W7-X, for a non-ideal MHD model including a uniform resistivity, electron inertia and (numerical) viscosity, the same qualitative behaviour is observed. In particular, even if a resonant (m, n) = (1, -1) perturbation is initialized, the most unstable mode is the (m, n) = (-4, 4), where m and n are the poloidal and toroidal mode numbers, respectively. Other modes are destabilized due to geometric coupling. The growth rate of this instability scales as η 1/3 , where η is the plasma resistivity, thus suggesting that ECCD drives ideal MHD stable W7-X plasmas towards non-ideal marginality. An ideal magnetohydrodynamic analysis confirms the result. A fluid-kinetic hybrid version of the EUTERPE code shows that gyrokinetic ions have a stabilizing effect on these modes. For W7-X relevant collisionalities, the growth rate scales linearly with the electron skin depth, d e . Implications of our results for sawtoothing W7-X operation are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Species Separation and Hydrogen Streaming upon Shock Release from Polystyrene under Inertial Confinement Fusion Conditions

Shock release from inertial confinement fusion (ICF) shells poses a great challenge to single-fluid hydrodynamic equations, especially for describing materials composed of different ion species. This has been evidenced by a recent experiment [Haberberger et al., Phys. Rev. Lett. 123, 235001 (2019)], in which low-density plasmas ( 10 19 to 10 20 cm – 3 ) are measured to move far ahead of what radiation-hydrodynamic simulations predict. To understand such experimental observations, in this work we have performed large-scale nonequilibrium molecular-dynamics simulations of shock release in polystyrene (CH) at experimental conditions. These simulations revealed that upon shock releasing from the back surface of a CH foil, hydrogen can stream out of the bulk of the foil due to its mass being lighter than carbon. This released hydrogen, exhibiting a much broader velocity distribution than carbon, forms low-density plasmas moving in nearly constant velocities ahead of the in-flight shell, which is in quantitative agreement with the experimental measurements. Such kinetic effect of species separation is currently missing in single-fluid radiation-hydrodynamics codes for ICF simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2008

Simulating Sources of Superstorm Plasmas

We evaluated the contributions to magnetospheric pressure (ring current) of the solar wind, polar wind, auroral wind, and plasmaspheric wind, with the surprising result that the main phase pressure is dominated by plasmaspheric protons. We used global simulation fields from the LFM single fluid ideal MHD model. We embedded the Comprehensive Ring Current Model within it, driven by the LFM transpolar potential, and supplied with plasmas at its boundary including solar wind protons, polar wind protons, auroral wind O+, and plasmaspheric protons. We included auroral outflows and acceleration driven by the LFM ionospheric boundary condition, including parallel ion acceleration driven by upward currents. Our plasmasphere model runs within the CRCM and is driven by it. Ionospheric sources were treated using our Global Ion Kinetics code based on full equations of motion. This treatment neglects inertial loading and pressure exerted by the ionospheric plasmas, and will be superceded by multifluid simulations that include those effects. However, these simulations provide new insights into the respective role of ionospheric sources in storm-time magnetospheric dynamics.

Fok, Mei-Ching↗
NASA NTRS1991

Plasma convection and ion beam generation in the plasma sheet boundary layer

Because of the dawn-dusk electric field E(dd), plasma in the magnetotail convects from the lobe toward the central plasma sheet (CPS). In the absence of space or velocity diffusion due to plasma turbulence, convection would yield a steady state distribution function f = V exp (-2/3) g(v exp 2 V exp 2/3), where V is the flux tube volume. Starting with such a distribution function and a plasma beta which varies from beta greater than 1 in the CPS to beta much smaller than 1 in the lobe, the evolution of the ion distribution function was studied considering the combined effects of ion diffusion by kinetic Alfven waves (KAW) in the ULF frequency range (1-10 mHz) and convection due to E(dd) x B drift in the plasma sheet boundary layer (PSBL) and outer central plasma sheet (OCPS). The results show that, during the early stages after launching the KAWs, a beamlike ion distribution forms in the PSBL and at the same time the plasma density and temperature decrease in the OCPS. Following this stage, ions in the beams convect toward the CPS resulting in an increase of the plasma temperature in the OCPS.

Moghaddam-Taaheri, E.↗
NASA NTRS1995

The sheath/ionosphere boundary layer at Venus

At Venus the interaction of the shocked solar wind and cold planetary ions takes place in the dayside mantle. The shocked solar wind is a warm, drifting Maxwellian plasma whereas the planetary plasma is cold; the plasma in the mantle is strongly magnetized. The coexistence of these two populations is unstable, and it leads to wave excitations that organize the energy and momentum exchange between the shocked solar wind and the plasma of planetary origin. The source of the free energy is the solar wind. The intensive wave activity seen in the 100 Hz channel of the wave instrument onboard the Pioneer-Venus Orbiter in the dayside mantle region of Venus can be identified as almost electrostatic VLF waves excited by the kinetic branch of the modified two-stream lower hybrid instability. The waves interact with the particles, and the planetary plasma is heated and accelerated outside the ionosphere, close to its upper boundary. This way solar wind scavenges the ionosphere, and planetary ions leave the planetary magnetosphere. A portion of the wave energy is capable of penetrating directly into the ionosphere and heating it.

Szego, K.↗
DOE OSTI2025

Interactions of a Microdischarge with a Premixed Flame (Final Technical Report)

In this project we introduced an experimental platform that facilitates the study of the two-way interaction of plasma and flame, and the exploration of the impact of the spatial location of the discharge on its kinetics, in a controlled environment. The experimental platform is a mesoscale CH4/air premixed flame burner that also acts as a dielectric barrier discharge (DBD) operated by nanosecond pulsed voltage. A unique feature of this setup is that the DBD utilizes transparent electrodes, made of Indium Tin Oxide (ITO), to unlock optical access along the line-of-sight parallel to the dominant electric field. This platform is suitable to study both the forward problem (impact of plasma-on-flame), which manifests itself as a modified flame speed as well as flame wrinkling; and the backward problem (impact of flame-on-plasma), which includes plasma regime transitions, between microdischarges and uniform modes, driven by the flame passage. Prior work at MIT has focused on electrical measurements and direct imaging studies of the plasma-flame interaction through the ITO electrode. This project, which was awarded experimental runtime at the Low-temperature Plasma Research Facility (PRF) at Sandia National Laboratories, quantifies the flame-plasma interaction through temporally resolved measurement of key radicals (oxygen atom, O) produced by the nanosecond repetitively pused discharge (NRPD) in an individual microdischarge, using two-photon absorption laser-induced fluorescence (TALIF).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Particle drift in the Earth's plasma sheet

We generalize the derivation of the average gradient/curvature-drift for a flux tube filled with an isotropic distribution of particles at specified kinetic energy. The present treatment is restricted to a two-dimensional magnetic field with zero electric field, but it includes all chaotic and Speiser orbits, which do not correspond to the simple picture of gradient/curvature drift. We assume that particles are evenly distributed throughout the regions of phase space allowed by their energy and canonical momentum. This assumption is closely related but not exactly equivalent to the assumption of isotropic pitch-angle distribution. Our derivation assumes that the maximum Larmor radius is small compared to the scale length for equatorial variations in the flux tube volume, but it does not involve any restrictions on the curvature of the field line. The resulting expression for the drift rate is valid for situations where the particle drift velocity is comparable to the thermal speed in some regions. The apparent implication of this generalized treatment is that the existence of very complex non-adiabatic particle trajectories in the plasma sheet may not invalidate previous estimates of the average rate of particle drift out the sides of the tail, estimates that were made under the assumption of simple guiding-center drifts.

Wolf, R. A.↗
DOE OSTI2022

Electron Weibel instability induced magnetic fields in optical-field ionized plasmas

Generation and amplification of magnetic fields in plasmas is a long-standing topic that is of great interest to both plasma and space physics. The electron Weibel instability is a well-known mechanism responsible for self-generating magnetic fields in plasmas with temperature anisotropy and has been extensively investigated in both theory and simulations, yet experimental verification of this instability has been challenging. Recently, we demonstrated a new experimental platform that enables controlled initialization of highly nonthermal and/or anisotropic plasma electron velocity distributions via optical-field ionization. Using an external electron probe bunch from a linear accelerator, the onset, saturation, and decay of the self-generated magnetic fields due to electron Weibel instability were measured for the first time to our knowledge. In this paper, we will first present experimental results on time-resolved measurements of the Weibel magnetic fields in non-relativistic plasmas produced by Ti:Sapphire laser pulses (0.8 μm) and then discuss the feasibility of extending the study to a quasi-relativistic regime by using intense CO 2 (e.g., 9.2 μm) lasers to produce much hotter plasmas.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Kinetic simulation study of magnetized collisionless shock formation on a terawatt laser system

Perpendicular, magnetized, collisionless shocks in hydrogen and neon plasmas are studied with 2-D particle-in-cell simulations for parameters accessible to experiments on OMEGA EP. The simulations are performed with realistic ion-electron mass ratios by which the relative importance of different micro-instabilities can be accurately captured. The dispersion relation of the modified two-stream instability (MTSI), the main dissipation mechanism for these shocks, is used to find suitable parameters for upcoming experiments. Simulations show that magnetized collisionless shocks can be readily formed within a few tenths of an ion gyro-period in both hydrogen and neon gases, with a background magnetic field of 50 T, achievable using the magneto-inertial fusion electrical discharge system (MIFEDS). Here, a portion of the ions are reflected to the upstream region and accelerated in both shock normal and tangential directions, indicating the formation of a supercritical shock. Shock front reformation is seen in longer time 1-D simulations. The results show that the formation time and width of these shocks are determined by MTSI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electron cyclotron drift instability and anomalous transport: two-fluid moment theory and modeling

In the presence of a strong electric field perpendicular to the magnetic field, the electron cross-field (E × B) flow relative to the unmagnetized ions can cause the so-called electron cyclotron drift instability (ECDI) due to resonances of the ion acoustic mode and the electron cyclotron harmonics. This occurs in, for example, collisionless shock ramps in space, and in E × B discharge devices such as Hall thrusters. A prominent feature of ECDI is its capability to induce an electron flow parallel to the background >E field at a speed greatly exceeding predictions by classical collision theory. Such anomalous transport is important due to its role in particle thermalization at space shocks, and in causing plasma flows towards the walls of E × B devices, leading to unfavorable erosion and performance degradation, etc. The development of ECDI and anomalous transport is often considered requiring a fully kinetic treatment. In this work, however, we demonstrate that a reduced variant of this instability, and more importantly, the associated anomalous transport, can be treated self-consistently in a collisionless two-fluid framework without any adjustable collision parameter. By treating both electron and ion species on an equal footing, the free energy due to the inter-species velocity shear allows the growth of an anomalous electron flow parallel to the background E field. We will first present linear analyses of the instability in the two-fluid five- and ten-moment models, and compare them against the fully-kinetic theory. At low temperatures, the two-fluid models predict the fastest-growing mode in good agreement with the kinetic result. Also, by including more (> = 10) moments, secondary (and possibly higher) unstable branches can be recovered. The dependence of the instability on ion-to-electron mass ratio, plasma temperature, and background B field strength is also thoroughly explored. We then carry out direct numerical simulations of the cross-field setup using the five-moment model. The development of the instability, as well as the anomalous transport, is confirmed and in excellent agreement with theoretical predictions. The force balance properties are also studied using the five-moment simulation data. Here this work casts new insights into the nature of ECDI and the associated anomalous transport and demonstrates the potential of the two-fluid moment model in efficient modeling of E × B plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The plasma–sheath transition and Bohm criterion in a high recycling divertor

The high recycling regime of a divertor is characterized by high plasma particle fluxes and low temperature at the target, where a strong hydrogen recirculation loop exists. Atomic processes in the high recycling regime, such as ion–neutral friction and radiation, can affect the plasma momentum and energy transport in the sheath transition region. Here, in this work, the plasma–sheath transition near a high recycling wall is investigated. The Bohm speed, which constraints the ion exit flow speed, is evaluated from a transport model that accounts for the effect of the anisotropic transport and atomic collisions in the transition layer. A first principles kinetic code vector particle-in-cell with the atomic collision package is used to investigate a 1D self-consistent slab plasma with a high recycling boundary for the tungsten and carbon divertors. The results demonstrate the accuracy of the Bohm speed model in predicting the ion exit flow speed in the transition region, as well as the reduction of the Bohm speed due to the ion–neutral friction. The effect of different wall materials, tungsten, and carbon, on the Bohm speed and near-wall plasma profile is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2011

Modeling the Inner Magnetosphere: Radiation Belts, Ring Current, and Composition

The space environment is a complex system defined by regions of differing length scales, characteristic energies, and physical processes. It is often difficult, or impossible, to treat all aspects of the space environment relative to a particular problem with a single model. In our studies, we utilize several models working in tandem to examine this highly interconnected system. The methodology and results will be presented for three focused topics: 1) Rapid radiation belt electron enhancements, 2) Ring current study of Energetic Neutral Atoms (ENAs), Dst, and plasma composition, and 3) Examination of the outflow of ionospheric ions. In the first study, we use a coupled MHD magnetosphere - kinetic radiation belt model to explain recent Akebono/RDM observations of greater than 2.5 MeV radiation belt electron enhancements occurring on timescales of less than a few hours. In the second study, we present initial results of a ring current study using a newly coupled kinetic ring current model with an MHD magnetosphere model. Results of a dst study for four geomagnetic events are shown. Moreover, direct comparison with TWINS ENA images are used to infer the role that composition plays in the ring current. In the final study, we directly model the transport of plasma from the ionosphere to the magnetosphere. We especially focus on the role of photoelectrons and and wave-particle interactions. The modeling methodology for each of these studies will be detailed along with the results.

Glocer, Alex↗
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.↗
DOE OSTI2024

Experimental validation of a collision-radiation dataset for molecular hydrogen in plasmas

Quantitative spectroscopy of molecular hydrogen has generated substantial demand, leading to the accumulation of diverse elementary process data encompassing radiative transitions, electron-impact transitions, predissociations, and quenching. However, their rates currently available are still sparse, and there are inconsistencies among those proposed by different authors. In this study, we demonstrate an experimental validation of such a molecular dataset by composing a collisional-radiative model (CRM) for molecular hydrogen and comparing experimentally obtained vibronic populations across multiple levels. From the population kinetics of molecular hydrogen, the importance of each elementary process in various parameter space is studied. In low-density plasmas (electron density ne≲1017 m−3) the excitation rates from the ground states and radiative decay rates, both of which have been reported previously, determine the excited state population. The inconsistency in the excitation rates affects the population distribution the most significantly in this parameter space. However, in higher density plasmas (ne≳1018 m−3), the excitation rates from excited states become important, which have never been reported in the literature, and may need to be approximated in some way. In order to validate these molecular datasets and approximated rates, we carried out experimental observations for two different hydrogen plasmas; a low-density radio frequency heated plasma (ne≈1016 m−3) and the Large Helical Device (LHD) divertor plasma (ne≳1018 m−3). The visible emission lines from EF1Σg+, HH¯1Σg+, D1Πu±, GK1Σg+, I1Πg±, J1Δg±, h3Σg+, e3Σu+, d3Πu±,g3Σg+, i3Πg±, and j3Δg± states were observed simultaneously and their population distributions were obtained from their intensities. We compared the observed population distributions with the CRM prediction, in particular the CRM with the rates compiled by Janev et al., Miles et al., and those calculated with the molecular convergent close-coupling (MCCC) method. The MCCC prediction gives the best agreement with the experiment, particularly for the emission from the low-density plasma. However, the population distribution in the LHD divertor shows a worse agreement with the CRM than those from low-density plasma, indicating the necessity of the precise excitation rates from excited states. We also found that the rates for the electron attachment is inconsistent with experimental results. This requires further investigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Second-order spectral line shift comparisons

The second-order spectral line width formulae from the projection operator and kinetic theory methods were recently compared. It was shown that a systematic expansion of the projection operator width expression including initial correlations formally agrees with the second-order kinetic theory result. It is now shown that the second-order dynamic shifts are also formally the same. The static shifts, however, differ due to an ad hoc treatment of electron-electron correlations in the projection operator method. The approximation is necessary in order to screen the radiator-electron interactions. The differences, however, are expected to be small. Finally, the results suggest using the rigorous and more compact second-order width and shift expressions from the kinetic theory method as the starting point for spectral line shape calculations. At line center, however, the projection operator second-order expression for the width and shift simplifies and reduces to the kinetic theory result.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗