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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Hybrid gyrokinetic ion/fluid electron simulation of toroidal tearing modes

The effects of toroidicity and kinetic ions on the resistive tearing mode are systematically studied with the gyrokinetic particle-in-cell simulation code GEM and compared with analytic theory. A new field solver in toroidal geometry has been developed for the simulation of low-n (n = 1, 2) modes in tokamaks. It is found that the toroidal effect significantly reduces the growth rate of the tearing mode. The toroidal effect can also increase the radial width of the tearing mode and change the scaling between the radial mode width and resistivity due to the toroidal pressure term in the electron continuity equation. The kinetic effects of ions can decrease the growth rate of the tearing mode. The plasma flux-surface shaping is found to have significant effect on the tearing mode.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Design elements and first data from a new Doppler backscattering system on the MAST-U spherical tokamak

A new Doppler backscattering (DBS) system has been installed and tested on the MAST-U spherical tokamak. It utilizes eight simultaneous fixed frequency probe beams (32.5, 35, 37.5, 40, 42.5, 45, 47.5, and 50 GHz). These frequencies provide a range of radial positions from the edge plasma to the core depending on plasma conditions. The system utilizes a combination of novel features to provide remote control of the probed density wavenumber, the launched polarization (X vs O-mode), and the angle of the launched DBS to match the magnetic field pitch angle. The range of accessible density turbulence wavenumbers (kθ) is reasonably large with normalized wavenumbers kθρs ranging from ≤0.5 to 9 (ion sound gyroradius ρs = 1 cm). This wavenumber range is relevant to a variety of instabilities believed to be important in establishing plasma transport (e.g., ion temperature gradient, trapped electron, electron temperature gradient, micro-tearing, kinetic ballooning modes). The system is specifically designed to address the requirement of density fluctuation wavevector alignment which can significantly reduce the SNR if not accounted for.

Instruments & Instrumentation↗
DOE OSTI2021

Novel Relaxation Time Approximation to the Relativistic Boltzmann Equation

In this work, we show that the widely used relaxation time approximation to the relativistic Boltzmann equation contains basic flaws, being incompatible with micro- and macroscopic conservation laws if the relaxation time depends on energy or general matching conditions are applied. We propose a new approximation that fixes such fundamental issues and maintains the basic properties of the linearized Boltzmann collision operator. We show how this correction affects transport coefficients, such as the bulk viscosity and particle diffusion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Kinetic Simulations of Instabilities and Particle Acceleration in Cylindrical Magnetized Relativistic Jets

Relativistic magnetized jets, such as those from AGN, GRBs, and XRBs, are susceptible to current- and pressure-driven MHD instabilities that can lead to particle acceleration and nonthermal radiation. Here, we investigate the development of these instabilities through 3D kinetic simulations of cylindrically symmetric equilibria involving toroidal magnetic fields with electron–positron pair plasma. Generalizing recent treatments by Alves et al. and Davelaar et al., we consider a range of initial structures in which the force due to toroidal magnetic field is balanced by a combination of forces due to axial magnetic field and gas pressure. We argue that the particle energy limit identified by Alves et al. is due to the finite duration of the fast magnetic dissipation phase. We find a rather minor role of electric fields parallel to the local magnetic fields in particle acceleration. In all investigated cases, a kink mode arises in the central core region with a growth timescale consistent with the predictions of linearized MHD models. In the case of a gas-pressure-balanced (Z-pinch) profile, we identify a weak local pinch mode well outside the jet core. We argue that pressure-driven modes are important for relativistic jets, in regions where sufficient gas pressure is produced by other dissipation mechanisms.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗
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↗
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↗
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↗
DOE OSTI2020

Driven electrostatic phase space vortices in a 1D weakly dissipative Vlasov–Poisson system

In this paper, the effect of collisions on driven electrostatic phase space vortices is analyzed by means of Eulerian simulation for two different collision models. It was demonstrated recently that in the absence of collisions, at late times, steady state phase space vortices manifest to form a plateau in the resonant region of the particle velocity distribution function, due to trapping of particles supporting multiextrema giant phase space vortices (PSVs). In the presence of collisions, over long time, this multiextrema plateau are found to smooth out, since collisions drive the velocity distribution toward Maxwellian, irrespective of how weak the collisions are as long as they are non-zero. In these conditions, kinetic processes and collisionality are found to be in competition, and the evolution of the plasma is found, therefore, to be a result of nontrivial combination of these two effects. An attempt has been made by means of numerical simulations to study the effect of weak collisionality on the electrostatic driven phase space vortices with two types of collision operators: (1) Bhatnagar–Gross–Krook (Krook) collision operator, where the colliding particles can be treated as isolated pairs and, (2) Fokker–Planck (FP) type collision operator (Zakharov–Karpman) in one dimension, where many weak collisions lead to particle diffusion in velocity space. It is shown that depending on the collision model used, the nature of smoothing in velocity space of giant PSVs results in qualitatively very different phase space structures. However, irrespective of the collision model used, excess density fractions over 10% are retained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Spatially dependent modeling and simulation of runaway electron mitigation in DIII-D

New simulations with the Kinetic Orbit Runaway electron (RE) Code (KORC) show RE deconfinement losses to the wall during plasma scrape off are the primary current dissipation mechanism in DIII-D experiments with high-Z impurity injection, and not collisional slowing down. The majority of simulations also exhibit an increase in the RE beam energy due to acceleration by the induced toroidal electric field, even while the RE beam current is decreasing. In this study, KORC integrates RE orbits using the relativistic guiding center equations of motion and incorporates time-sequenced, experimental reconstructions of the magnetic and electric fields and line integrated electron density to construct spatiotemporal models of electron and partially ionized impurity transport in the companion plasma. Comparisons of experimental current evolution and KORC results demonstrate the importance of including Coulomb collisions with partially ionized impurity physics, initial RE energy, pitch angle, and spatial distributions, and spatiotemporal electron and partially ionized impurity transport. The research presented here provides an initial quantification of the efficacy of RE mitigation via injected impurities and identification of the critical role played by loss of confinement due to plasma scrape off on the inner wall as compared to the relatively slow collisional damping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗