Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PLASMA DYNAMICS”

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,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Spacecraft surface charging as a function of material properties

Spacecraft material behavior plays a very important role in space missions. Spacecraft immersed in plasma get charged by absorbing plasma particles and by emitting electrons from spacecraft surfaces via photoelectron and secondary electron emission. Spacecraft charging depends heavily on material properties such as work function, secondary electron yield, dielectric constant, and electric conductivity among other. Material properties are typically assumed to be static in charging models. However, it is well known that this is not the case in space. This makes spacecraft charging predictions very challenging. Material properties are well characterized before the spacecraft is put in orbit through characterization in the lab under clean conditions. However, once in space, material properties change due to the harsh and very dynamic space environment. We present a new capability to predict material behavior in space from first-principles modeling. The ongoing effort seeks to couple material models, density functional theory (DFT) and molecular dynamic (MD) codes, with environment models, plasma kinetic codes. This preliminary study will show results of surface charging as a function of material work function, dielectric constant, and conductivity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Interstitial hydrogen enhances the mobility of some grain boundaries in tungsten

Segregation of interstitials at a grain boundary (GB) is known to generally lower its mobility. This phenomenon, called 'solute-drag', has important ramifications on the process of recrystallization and microstructural evolution. Here, we present predictions from molecular dynamics (MD) simulations which demonstrate that interstitial hydrogen in tungsten can in fact increase the mobility of some GBs which exhibit shear coupling. Assuming a disconnection-based mechanism, activation energies and pre-factors for disconnection nucleation are predicted from simulations of shear-coupled motion. In GBs where enhanced mobility is predicted, interstitial H reduces both the activation energy and the pre-factor for disconnection nucleation, thus effectively increasing the mobility. For GBs with diminished mobility, MD predicts that presence of interstitial H reduces the pre-factor and, in some cases, increases the activation energy. The reduction in the activation energy inferred from MD simulations are confirmed by nudged elastic band calculations. Temperature-dependent structural transitions are observed for some GBs, and the effect of interstitial H is found to change with the changes in structure. The effect of interstitial H is predicted to be complex and highly variable, providing some plausible explanations for experimental observations on the recrystallization of tungsten in presence of H-loaded plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Effects of the relative motion of different particles on the wave instability in dusty plasmas

Theoretical study is carried out for the dust acoustic waves in a dusty plasma. The dispersion relation is obtained from the fluid dynamical model. It is found that the wave frequency depends on the electron number density, wave number, and the speed of electrons or ions. The dependencies of the instability on the system parameters are also discussed. It is found that the long wavelength wave is usually unstable, which may be useful in the plasma physics, especially in the magnetic confined fusion system or the charged beam dynamics in the accelerated driven system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

ContourNet: Salient Local Contour Identification for Blob Detection in Fusion Plasma Simulation Data

We present ContourNet, a deep learning approach to identify salient local isocontours as blobs in large-scale 5D gyrokinetic tokamak simulation data. Blobs-regions of high turbulence that run along the edge wall down toward the diverter and can damage the tokamak-are non-well-defined features but have been empirically localized by isocontours in 2D normalized fluctuating density fields. The key of our study is to train ContourNet to follow the empirical rules to detect blobs over the time-varying simulation data. The architecture of ContourNet is a convolutional neural segmentation network: the inputs are the density field and a rasterized isocontour; the output is a set of isocontour encircling blobs. At the training stage, we feed the network with manually identified isocontours and propagated labels. At the inference stage, we extract isocontours from the segmented blob regions. Results show that our approach can achieve both high accuracy and performance, which enables scientists to understand the blob dynamics influencing the confinement of the plasma.

Blob detection↗
DOE OSTI2015

Nonlinear dynamics of the ion Weibel-filamentation instability: An analytical model for the evolution of the plasma and spectral properties

We present a predictive model of the nonlinear phase of the Weibel instability induced by two symmetric, counter-streaming ion beams in the non-relativistic regime. This self-consistent model combines the quasilinear kinetic theory of Davidson et al. [Phys. Fluids 15, 317 (1972)] with a simple description of current filament coalescence. It allows us to follow the evolution of the ion parameters up to a stage close to complete isotropization, and is thus of prime interest to understand the dynamics of collisionless shock formation. Its predictions are supported by 2-D and 3-D particle-in-cell simulations of the ion Weibel instability. The derived approximate analytical solutions reveal the various dependencies of the ion relaxation to isotropy. In particular, it is found that the influence of the electron screening can affect the results of simulations using an unphysical electron mass.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Plasma-assisted atomic layer etching of single-crystal diamond

Applications of near-surface nitrogen-vacancy (NV) centers in diamond are often limited by surface defects created during processing. Understanding and controlling plasma-induced surface damage is important for preserving the optical and spin properties of diamond NV centers. We report molecular dynamics simulations of a novel form of plasma-aided atomic layer etching of diamond. In this proposed scheme, the initial surface modification step consists of Ar + ion bombardment. This creates an amorphous layer at the surface, the thickness of which is controlled by the ion energy. Amorphization is known to help smooth the surface, at least locally, and would also serve to sputter clean an initially contaminated surface. A second step impacts the amorphous carbon layer with O + between about 1 and 5 eV. Simulations show that this energy range will remove the amorphous carbon but will not etch the underlying diamond. Though this energy range is not trivial to achieve in a conventional low-temperature plasma, potential methods for creating these low-energy O + ions are discussed. In addition to a-C etching, the O + impacts are predicted to remove any isolated (100) diamond terraces by selectively attacking the edges of the terraces. The proposed atomic layer etching (ALE) approach reverses the conventional ALE sequence in which the surface modification step is usually chemical modification (oxidation in this case), followed by a removal step using Ar + impacts. This plasma ALE procedure is predicted to create a diamond surface that is atomically flat and defect free.

Draney, J. S. [Princeton University, NJ (United St↗
DOE OSTI2018

Investigation of the dynamics of ionization induced injected electrons under the influence of beam loading effects

In laser-driven wakefield, ionization induced injection is an efficient way to inject electrons in the plasma wave. A detailed study on the beam dynamics under the influence of beam loading effects, which can be controlled by the concentration of nitrogen impurity introduced in the hydrogen gas was conducted. For a specific value of this percentage, the final energy of the high-energy electron bunch becomes nearly independent of the trapped positions, thus leading to a small energy dispersion. We also show that the final beam emittance is mainly determined by the injection process.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

The charge imbalance in ultracold plasmas

Ultracold plasmas are regarded as quasineutral but not strictly neutral. The results of charge imbalance in the expansion of ultracold plasmas are reported. The calculations are performed by a full molecular-dynamics simulation. The details of the electron velocity distributions are calculated without the assumption of electron global thermal equilibrium and Boltzmann distribution. Spontaneous evolutions of the charge imbalance from the initial states with perfect neutrality are given in the simulations. The expansion of outer plasma slows down with the charge imbalance. The influences of plasma size and parameters on the charge imbalance are discussed. The radial profiles of electron temperature are given for the first time, and the self-similar expansion can still occur even if there is no global thermal equilibrium. The electron disorder induced heating is also found in the simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1980

Plasmas in Saturn's magnetosphere

The solar wind plasma analyzer on board Pioneer 2 provides first observations of low-energy positive ions in the magnetosphere of Saturn. Measurable intensities of ions within the energy-per-unit charge (E/Q) range 100 eV to 8 keV are present over the planetocentric radial distance range about 4 to 16 R sub S in the dayside magnetosphere. The plasmas are found to be rigidly corotating with the planet out to distances of at least 10 R sub S. At radial distances beyond 10 R sub S, the bulk flows appear to be in the corotation direction but with lesser speeds than those expected from rigid corotation. At radial distances beyond the orbit of Rhea at 8.8 R sub S, the dominant ions are most likely protons and the corresponding typical densities and temperatures are 0.5/cu cm and 1,000,000 K, respectively, with substantial fluctuations. It is concluded that the most likely source of these plasmas in the photodissociation of water frost on the surface of the ring material with subsequent ionization of the products and radially outward diffusion. The presence of this plasma torus is expected to have a large influence on the dynamics of Saturn's magnetosphere since the pressure ratio beta of these plasmas approaches unity at radial distances as close to the planet as 6.5 R sub S. On the basis of these observational evidences it is anticipated that quasi-periodic outward flows of plasma, accompanied with a reconfiguration of the magnetosphere beyond about 6.5 R sub S, will occur in the local night sector in order to relieve the plasma pressure from accretion of plasma from the rings.

Frank, L. A.↗
DOE OSTI2024

Generalized gluon distribution for quarkonium dynamics in strongly coupled N = 4 Yang-Mills theory

We study the generalized gluon distribution that governs the dynamics of quarkonium inside a non-Abelian thermal plasma characterizing its dissociation and recombination rates. This gluon distribution can be written in terms of a correlation function of two chromoelectric fields connected by an adjoint Wilson line. We formulate and calculate this object in N = 4 supersymmetric Yang-Mills theory at strong coupling using the AdS/CFT correspondence, allowing for a nonzero center-of-mass velocity v of the heavy quark pair relative to the medium. The effect of a moving medium on the dynamics of the heavy quark pair is described by the simple substitution T → γ T in agreement with previous calculations of other observables at strong coupling, where T is the temperature of the plasma in its rest frame, and γ = ( 1 − v 2 ) − 1 / 2 is the Lorentz boost factor. Such a velocity dependence can be important when the quarkonium momentum is larger than its mass. Contrary to general expectations for open quantum systems weakly coupled with large thermal environments, the contributions to the transition rates that are usually thought of as the leading ones in Markovian descriptions vanish in this strongly coupled plasma. This calls for new theoretical developments to assess the effects of strongly coupled non-Abelian plasmas on in-medium quarkonium dynamics. Finally, we compare our results with those from weakly coupled QCD, and find that the QCD result moves toward the N = 4 strongly coupled result as the coupling constant is increased within the regime of applicability of perturbation theory. This behavior makes it even more pressing to develop a non-Markovian description of quarkonium in-medium dynamics. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Bursty magnetic reconnection at the Earth's magnetopause triggered by high-speed jets

The impact of high-speed jets—dynamic pressure enhancements in the magnetosheath—on the Earth's magnetopause has been observed to trigger local magnetic reconnection. We perform a three-dimensional hybrid simulation to study the magnetosheath and magnetopause under turbulent conditions using a quasi-radial southward interplanetary magnetic field (IMF). In contrast to quasi-steady reconnection with a strong southward IMF, we show that after the impact of a jet on the magnetopause, the magnetopause moves inward, the current sheet is compressed and intensified and signatures of local magnetic reconnection are observed, showing similarities to spacecraft measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

An equilibrium-preserving discretization for the nonlinear Rosenbluth–Fokker–Planck operator in arbitrary multi-dimensional geometry

Here, the Fokker–Planck collision operator is an advection-diffusion operator which describe dynamical systems such as weakly coupled plasmas, photonics in high temperature environment biological, and even social systems. For plasmas in the continuum, the Fokker–Planck collision operator supports such important physical properties as conservation of number, momentum, and energy, as well as positivity. It also obeys the Boltzmann's H-theorem, i.e., the operator increases the system entropy while simultaneously driving the distribution function towards a Maxwellian. In the discrete, when these properties are not ensured, numerical simulations can either fail catastrophically or suffer from significant numerical pollution. There is strong emphasis in the literature on developing numerical techniques to solve the Fokker–Planck equation while preserving these properties. In this short note, we focus on the analytical equilibrium preserving property, meaning that the Fokker–Planck collision operator vanishes when acting on an analytical Maxwellian distribution function. The equilibrium preservation property is especially important, for example, when one is attempting to capture subtle transport physics. Since transport arises from small Ο(ϵ) corrections to the equilibrium (where ϵ is a small expansion parameter), numerical truncation error present in the equilibrium solution may dominate, overwhelming transport dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Hydrodynamic expansion and near-infrared absorption of x-ray heated aluminum plasmas

We use x-ray pulses from dense argon plasmas at the Z Machine (Sandia National Laboratories) to generate hypersonic aluminum plasmas akin to material ejecta during proposed planetary defense missions, fusion reactor wall excursions, and other high-energy density processes. Near-infrared absorption is used to diagnose the controlled expansion of the plasmas through cylindrical cavities following their generation from x-ray heating of solid aluminum 7075 alloy. The data are compared to multidimensional radiation hydrodynamics simulations utilizing the ALEGRA multiphysics code, accounting for the dynamics of radiation scattering, material phase change, plasma expansion, thermal re-irradiation, and interactions with the cavity and with the infrared beams. To allow for accurate simulation, density functional theory is used to apply the Hagen–Rubens relation for the far-infrared and is adjoined with a detailed configuration accounting model using the Propaceos code, producing opacities spanning 10 −1 –10 4 eV photon energy for aluminum 7075 alloy, and in comparison with pure aluminum. The model is found to agree with experimental data in the higher-fluence regime when the Hagen–Rubens relation is applied. The ejected material, which is observed to travel up to 55 km/s, is comprised of a strongly ionized, non-LTE plasma front at ∼10 eV temperature followed by a weakly ionized LTE gas at higher density. The present findings lend some confidence to the broad-range equation of state and infrared opacity models for weakly ionized aluminum plasmas while demonstrating an approach to their future refinement, with potential application to astrophysical plasmas and other extreme processes.

Adiabatic process↗
DOE OSTI2014

Coupling strength in Coulomb and Yukawa one-component plasmas

In a non-ideal classical Coulomb one-component plasma (OCP), all thermodynamic properties are known to depend only on a single parameter—the coupling parameter Γ. In contrast, if the pair interaction is screened by background charges (Yukawa OCP) the thermodynamic state depends, in addition, on the range of the interaction via the screening parameter κ. How to determine in this case an effective coupling parameter has been a matter of intensive debate. Here we propose a consistent approach for defining and measuring the coupling strength in Coulomb and Yukawa OCPs based on a fundamental structural quantity, the radial pair distribution function (RPDF). The RPDF is often accessible in experiments by direct observation or indirectly through the static structure factor. Alternatively, it is directly computed in theoretical models or simulations. Our approach is based on the observation that the build-up of correlation from a weakly coupled system proceeds in two steps: First, a monotonically increasing volume around each particle becomes devoid of other particles (correlation hole), and second (upon further increase of the coupling), a shell structure emerges around each particle giving rise to growing peaks of the RPDF. Using molecular dynamics simulation, we present a systematic study for the dependence of these features of the RPDF on Γ and κ and derive a simple expression for the effective coupling parameter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Compression dynamics of quasi-spherical wire arrays with different linear mass profiles

Results of experimental studies of the implosion of quasi-spherical wire (or metalized fiber) arrays are presented. The goal of the experiments was to achieve synchronous three-dimensional compression of the plasma produced in different regions of a quasi-spherical array into its geometrical center. To search for optimal synchronization conditions, quasi-spherical arrays with different initial profiles of the linear mass were used. The following dependences of the linear mass on the poloidal angle were used: m{sub l}(θ) ∝ sin{sup –1}θ and m{sub l}(θ) ∝ sin{sup –2}θ. The compression dynamics of such arrays was compared with that of quasi-spherical arrays without linear mass profiling, m{sub l}(θ) = const. To verify the experimental data, the spatiotemporal dynamics of plasma compression in quasi-spherical arrays was studied using various diagnostics. The experiments on three-dimensional implosion of quasi-spherical arrays made it possible to study how the frozen-in magnetic field of the discharge current penetrates into the array. By measuring the magnetic field in the plasma of a quasi-spherical array, information is obtained on the processes of plasma production and formation of plasma flows from the wire/fiber regions with and without an additionally deposited mass. It is found that penetration of the magnetic flux depends on the initial linear mass profile m{sub l}(θ) of the quasi-spherical array. From space-resolved spectral measurements and frame imaging of plasma X-ray emission, information is obtained on the dimensions and shape of the X-ray source formed during the implosion of a quasi-spherical array. The intensity of this source is estimated and compared with that of the Z-pinch formed during the implosion of a cylindrical array.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Gyrokinetic simulations of momentum flux parasitic to free-energy transfer

Ion Landau damping interacts with a portion of the E × B drift to cause a nondiffusive outward flux of co-current toroidal angular momentum. Quantitative evaluation of this momentum flux requires nonlinear simulations to determine f L , the fraction of fluctuation free energy that passes through ion Landau damping, in fully developed turbulence. Nonlinear gyrokinetic simulations with the GKW code confirm the presence of the systematic symmetry-breaking momentum flux. For simulations with adiabatic electrons, f L scales inversely with the ion temperature gradient, because only the ion curvature drift can transfer free energy to the electrostatic potential. Although kinetic electrons should, in principle, relax this restriction, the ion Landau damping measured in collisionless kinetic-electron simulations remained at low levels comparable with ion-curvature-drift transfer, except when magnetic shear $\hat{s}$ was strong. A set of simulations scanning the electron pitch-angle scattering rate showed only a weak variation of f L with the electron collisionality. Still, collisional-electron simulations with electron temperature greater than ion temperature unambiguously showed electron-curvature-drift transfer supporting ion Landau damping, leading to a corresponding enhancement of the symmetry-breaking momentum flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Multi-dimensional modeling of non-equilibrium plasma for automotive applications

While spark-ignition (SI) engine technology is aggressively moving towards challenging (dilute and boosted) combustion regimes, advanced ignition technologies generating non-equilibrium types of plasma are being considered by the automotive industry as a potential replacement for the conventional spark-plug technology. However, there are currently no models that can describe the low-temperature plasma (LTP) ignition process in computational fluid dynamics (CFD) codes that are typically used in the multi-dimensional engine modeling community. A key question for the engine modelers that are trying to describe the non-equilibrium ignition physics concerns the plasma characteristics. A key challenge is also represented by the plasma formation timescale (nanoseconds) that can hardly be resolved within a full engine cycle simulation. This paper reports on the multi-dimensional modeling of LTP generated by a nanopulsed high-voltage discharge in a pin-to-pin electrode configuration, and evaluates the effects of ambient pressure on the post-discharge results. It is shown that a nanopulsed delivery can result in a LTP or a mode transition to an arc-like event depending on the mixture properties in the gap between the two electrodes. It is further shown that numerical predictions of the LTP-to-arc transition point as a function of ambient pressure along with qualitative distributions of important radicals and temperature closely match recent experiments with identical geometry and discharge characteristics. As a result, the high-fidelity modeling effort described in this paper can be used as a key tool to deliver in-depth understanding of non-equilibrium plasmas for automotive applications and lead to the future implementation of dedicated models for ignition technologies that are based on non-equilibrium plasma physics.

Scarcelli, Riccardo↗
DOE OSTI2020

Particlization in fluid dynamical simulations of heavy-ion collisions: The IS3D module

The IS3D particlization module simulates the emission of hadrons from heavy-ion collisions via Monte-Carlo sampling of the Cooper–Frye formula which converts fluid dynamical information into local phase-space distributions for hadrons. The code package includes multiple choices for the non-equilibrium correction to these distribution functions: the 14-moment approximation, first-order Chapman–Enskog expansion, and two types of modified equilibrium distributions. This makes it possible to explore to what extent heavy-ion experimental data are sensitive to different choices for δf n , presently the main source of theoretical uncertainty in the particlization stage. Here, we validate our particle sampler with a high degree of precision by generating several million hadron emission events from a longitudinally boost-invariant hypersurface and comparing the event-averaged particle spectra and space–time distributions to the Cooper–Frye formula.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗