Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Instabilities”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Evolution of the electron cyclotron drift instability in two-dimensions

The Electron Cyclotron Drift Instability driven by the electron E × B drift in partially magnetized plasmas is investigated with highly resolved particle-in-cell simulations. The emphasis is on two-dimensional effects involving the parallel dynamics along the magnetic field in a finite length plasma with dielectric walls. It is found that the instability develops as a sequence of growing cyclotron harmonics demonstrating wave breaking and complex nonlinear interactions, being particularly pronounced in ion density fluctuations at short wavelengths. At the same time, nonlinear evolution of fluctuations of the ion and electron density, as well as the anomalous electron current, shows cascade toward long wavelengths. Tendency to generate long wavelength components is most clearly observed in the spectra of the electron density and the anomalous current fluctuations. An intense but slowly growing mode with a distinct eigen-mode structure along the magnetic field develops at a later nonlinear stage enhancing the tendency toward long wavelength condensation. Furthermore, the latter mode having a finite wavelength along the magnetic field is identified as the Modified Two-Stream Instability (MTSI). It is shown that the MTSI mode results in strong parallel heating of electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Amplification due to two-stream instability of self-electric and magnetic fields of an ion beam propagating in background plasma

Propagation of charged particle beams in background plasma as a method of space charge neutralization has been shown to achieve a high degree of charge and current neutralization and therefore enables nearly ballistic propagation and focusing of charged particle beams. Correspondingly, the use of plasmas for propagation of charged particle beams has important applications for transport and focusing of intense particle beams in inertial fusion and high energy density laboratory plasma physics. However, the streaming of beam ions through a background plasma can lead to the development of two-stream instability between the beam ions and the plasma electrons. The beam electric and magnetic fields enhanced by the two-stream instability can lead to defocusing of the ion beam. Using particle-in-cell simulations, we study the scaling of the instability-driven self-electromagnetic fields and consequent defocusing forces with the background plasma density and beam ion mass. In conclusion, we identify plasma parameters where the defocusing forces can be reduced.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Development of new platforms for hydrodynamic instability and asymmetry measurements in deceleration phase of indirectly driven implosions on NIF

Hydrodynamic instabilities and asymmetries are a major obstacle in the quest to achieve ignition at the National Ignition Facility (NIF) as they cause pre-existing capsule perturbations to grow and ultimately quench the fusion burn in experiments. This paper reviews the development of two new experimental techniques to measure high-mode instabilities and low-mode asymmetries in the deceleration phase of indirect drive inertial confinement fusion implosions. In the first innovative technique, self-emission from the hot spot was enhanced with an argon dopant to “self-backlight” the shell in-flight, imaging the perturbations in the shell near peak velocity. Experiments with pre-imposed two-dimensional perturbations showed hydrodynamic instability growth of up to 7000× in areal density. These experiments discovered unexpected three-dimensional structures originating from the capsule support structures. These new 3-D structures became one of the primary concerns for the indirect drive ICF program that requires their origin to be understood and their impact mitigated. In a second complementary technique, the inner surface of the decelerating shell was visualized in implosions using x-ray emission of a high-Z dopant added to the inner surface of the capsule. With this technique, low mode asymmetry and high mode perturbations, including perturbations seeded by the gas fill tube and capsule support structure, were quantified near peak compression. Using this doping method, the role of perturbations and radiative losses from high atomic number materials on neutron yield was quantified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Kelvin-Helmholtz instability is the result of parity-time symmetry breaking

Parity-Time (PT)-symmetry is being actively investigated as a fundamental property of observables in quantum physics. We show that the governing equations of the classical two-fluid interaction and the incompressible fluid system are PT-symmetric, and the well-known Kelvin-Helmholtz instability is the result of spontaneous PT-symmetry breaking. It is expected that all classical conservative systems governed by Newton's law admit PT-symmetry, and the spontaneous breaking thereof is a generic mechanism for classical instabilities. Here discovering the PT-symmetry of systems in fluid dynamics and plasma physics and identifying the PT-symmetry breaking responsible for instabilities enable new techniques to classical physics and enrich the physics of PT-symmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Removing instabilities in the hierarchical equations of motion: Exact and approximate projection approaches

The hierarchical equations of motion (HEOM) provide a numerically exact approach for computing the reduced dynamics of a quantum system linearly coupled to a bath. We have found that HEOM contains temperature-dependent instabilities that grow exponentially in time. In the case of continuous-bath models, these instabilities may be delayed to later times by increasing the hierarchy dimension; however, for systems coupled to discrete, nondispersive modes, increasing the hierarchy dimension does little to alleviate the problem. We show that these instabilities can also be removed completely at a potentially much lower cost via projection onto the space of stable eigenmodes; furthermore, we find that for discrete-bath models at zero temperature, the remaining projected dynamics computed with few hierarchy levels are essentially identical to the exact dynamics that otherwise might require an intractably large number of hierarchy levels for convergence. Recognizing that computation of the eigenmodes might be prohibitive, e.g., for large or strongly coupled models, we present a Prony filtration algorithm that may be useful as an alternative for accomplishing this projection when diagonalization is too costly. We present results demonstrating the efficacy of HEOM projected via diagonalization and Prony filtration. Finally, we discuss issues associated with the non-normality of HEOM.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Instability-enhanced transport in low temperature magnetized plasma

It is shown that the transport in low temperature, collisional, bound plasma is enhanced by instabilities at a high magnetic field. While the magnetic field confines the electrons in a stable plasma, the instability completely destroys the confinement such that the transport becomes independent of the magnetic field in the highly magnetized limit. Finally, an analytical expression of the instability-enhanced collision frequency is proposed, based on a magnetic field independent edge-to-center density ratio.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Multi-stream instability of a single long electron bunch in a storage ring

A multistream instability is observed experimentally during the longitudinal expansion of an electron beam in a storage ring. The instability is observed when the beam expands in length to several times the circumference of the ring so that portions of the beam overlap. While portions of the beam overlap in physical space, they form multiple streams and remain separate in velocity space. The streams become unstable as their number increases and their separation in velocity decreases. An analytical theory predicts the onset of the instability, in agreement with simulations and measurements, over a wide range of peak line-charge densities (10.3 pC/m–1.8 nC/m) and bunch lengths. Finally, this work extends previous calculations to include the dynamic nonlinear elongation of the bunch, with a given initial length, and defines an onset criterion for the filament velocity separation for the beam to be unstable.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Scalings pertaining to current sheet disruption mediated by the plasmoid instability

Analytic scaling relations are derived for a phenomenological model of the plasmoid instability in an evolving current sheet, including the effects of reconnection outflow. Two scenarios are considered, where the plasmoid instability can be triggered either by an injected initial perturbation or by the natural noise of the system (here referred to as the system noise). The two scenarios lead to different scaling relations because the initial noise decays when the linear growth of the plasmoid instability is not sufficiently fast to overcome the advection loss caused by the reconnection outflow, whereas the system noise represents the lowest level of fluctuations in the system. The leading order approximation for the current sheet width at disruption takes the form of a power law multiplied by a logarithmic factor, and from that, the scaling relations for the wavenumber and the linear growth rate of the dominant mode are obtained. When the effects of the outflow are neglected, the scaling relations agree, up to the leading order approximation, with previously derived scaling relations based on a principle of least time. As a result, the analytical scaling relations are validated with numerical solutions of the model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Multi-mode ionization instability induced striations in RF driven He/H 2 O atmospheric pressure plasma (APP) discharges

Previously [E. Kawamura et al., Plasma Sources Sci. Technol. 25, 054009 (2016) and E. Kawamura et al., J. Phys. D: Appl. Phys. 50, 145204 (2017)], one dimensional (1D) particle-in-cell (PIC) simulations of 1 to 4mm gap, He/2%H 2 O atmospheric pressure plasma discharges showed an ionization instability resulting in bulk striations. Assuming that the ionization rate coefficient K iz is related to the root mean square electric field E by K iz &prop; E q , a striation theory showed that q<0 is a necessary condition for the instability. A local calculation yielded q>0, implying that nonlocal electron kinetics are required for the instability. Wider gaps can fit in a wider range of wavelengths, resulting in multi-mode striations. Previously, we assumed a single mean q value for each discharge, and did not calculate q for each mode separately. Here, we develop a wavelength resolved global striation model and apply it to 1D PIC simulations of 4mm gap discharges with input currents J rf = 0.04 to 0.30A/cm 2 . We first examine a base case at 0.23A/cm 2 and observe a mixture of unstable modes within a window of wavelengths λ. At shorter λ, the modes are suppressed by diffusion. At longer λ, we observe a transition to locality in which q becomes less negative with increasing λ, approaching its local positive value and stabilizing the modes. The window of unstable modes shifts to shorter λ with increasing J rf , causing the modes to be suppressed by diffusion at higher J rf . At lower J rf , the decrease in bulk plasma density with decreasing J rf suppresses the striations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Laser-plasma instabilities in long scale-length plasmas relevant to shock-ignition

We present 2D particle-in-cell simulations of laser-plasma instabilities (LPI) performed using conditions relevant to experiments on the OMEGA and NIF laser facilities. The laser intensity used is 2×10 15 W cm -2 , which allows comparison with previous experiments and simulations investigating shock ignition, though this is a relatively low intensity for the scheme. We find that the large convective gain of the instabilities leads to dynamics that are largely controlled by pump depletion and to scattering or absorption by LPIs predominantly occurring away from the quarter-critical density. The two-plasmon decay (TPD) instability is dominant for OMEGA-type conditions and remains important at the NIF-scale. In both cases, most absorption from TPD occurs near its Landau cutoff density, and we examine its nonlinear dynamics in this region, focusing on the ion-acoustic waves responsible for saturation. For NIF-type conditions, stimulated Raman scattering (SRS) and stimulated Brillouin scattering divert significant fractions of laser energy, and we find that SRS backscatter can occur with high gain at low densities due to kinetic effects. An in-depth analysis of hot-electron production is performed to identify the main sources and their characteristic electron temperature. For both simulations presented, the overall hot-electron temperature lies between 30 and 35 keV. This is relatively low and will likely increase shock pressure; however, the distributions also contain a significant number of high-energy (E k >100 keV) electrons that would likely cause unacceptable preheat. We, therefore, suggest some strategies that may be utilized to minimize this high-energy component. I. INTROD

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Alfvénic modes excited by the kink instability in PHASMA

Magnetic flux ropes have been successfully created with plasma guns in the newly commissioned PHAse Space MApping (PHASMA) experiment. The flux ropes exhibit the expected m = 1 kink instability. The observed threshold current for the onset of this kink instability is half of the Kruskal–Shafranov current limit, consistent with predictions for the non-line tied boundary condition of PHASMA. The helicity, paramagnetism, and growth rate of the observed magnetic fluctuations are also consistent with kink instability predictions. The observed fluctuation frequency appears to be a superposition of a real frequency due to a Doppler shift of the kink mode arising from plasma flow (~2kHz) and a contribution from a wave mode (~5 kHz). The dispersion of the wave mode is consistent with an Alfvén wave. Distinct from most previ-ous laboratory studies of flux ropes, the working gas in PHASMA is argon. Thus, the ion cyclotron frequency in PHASMA is quite low and the frequency of the Alfvénic mode plateaus at ~0.5 of the ion gyro frequency with increasing background magnetic field strength.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Studying the Richtmyer–Meshkov instability in convergent geometry under high energy density conditions using the Decel platform

The “Decel” platform at Sandia National Laboratories investigates the Richtmyer–Meshkov instability (RMI) in converging geometry under high energy density conditions [Knapp et al., Phys. Plasmas 27, 092707 (2020)]. In Decel, the Z machine magnetically implodes a cylindrical beryllium liner filled with liquid deuterium, launching a converging shock toward an on-axis beryllium rod machined with sinusoidal perturbations. The passage of the shock deposits vorticity along the Be/D2 interface, causing the perturbations to grow. Here, we present platform improvements along with recent experimental results. To improve the stability of the imploding liner to the magneto Rayleigh–Taylor instability, we modified its acceleration history by shortening the Z electrical current pulse. Next, we introduce a “split rod” configuration that allows two axial modes to be fielded simultaneously in different axial locations along the rod, doubling our data per experiment. We then demonstrate that asymmetric slots in the return current structure modify the magnetic drive pressure on the surface of the liner, advancing the evolution on one side of the rod by multiple ns compared to its 180° counterpart. This effectively enables two snapshots of the instability at different stages of evolution per radiograph with small deviations of the cross-sectional profile of the rod from the circular. Using this platform, we acquired RMI data at 272 and 157 μm wavelengths during the single shock stage. Finally, we demonstrate the utility of these data for benchmarking simulations by comparing calculations using ALEGRA MHD and RageRunner.

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 OSTI2023

Transition from ITG to MTM linear instabilities near pedestals of high density plasmas

Here, investigation of linear gyrokinetic ion-scale modes ( k θ ρ s = 0.3) finds that a transition from ion temperature gradient to microtearing mode (MTM) dominance occurs as the density is increased near the pedestal region of a parameterized DIII-D sized tokamak. H-modes profile densities, temperatures, and equilibria are parameterized utilizing the OMFIT PRO_create module. With these profiles, linear gyrokinetic ion-scale instabilities are predicted with CGYRO. This transition (nMTM) has a weak dependence on radial location in the region near the top of the pedestal ( ρ = 0.7 – 0.9), which allows simulating single radii to examine the approximate scaling of nMTM with global parameters. The critical nMTM is found to scale with plasma current. Additionally, increasing the minor radius by decreasing the aspect ratio and increasing the major radius are found to reduce nMTM. However, any relationship between nMTM and density limit physics remains unclear as nMTM increases relative to the Greenwald density with larger minor radius and with larger magnetic field, suggesting that the transport due to MTM may be less important for a reactor. Additionally, nMTM is sensitive to the pedestal temperature, the local electron and ion gradients, the ratio of ion to electron temperature T i / T e, and the current profile. MTMs are predicted to be the dominant instability in the core at similar Greenwald fractions for DIII-D, NSTX, and NSTX-U H-mode experiments, supporting the results of the parameterized study. Additionally, MTMs continue to be the dominant linear instability in a DIII-D L-mode after an H–L transition as the plasma approaches a density limit disruption despite the large change in plasma profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Electrostatic upper-hybrid mode instability driven by a ring electron distribution

Quasi electrostatic fluctuations in the upper-hybrid frequency range are commonly detected in the planetary magnetospheric environment. The origin of such phenomena may relate to the instability driven by a loss-cone feature associated with the electrons populating the dipole-like magnetic field. The present paper carries out a one-dimensional electrostatic particle-in-cell simulation accompanied by a reduced quasilinear kinetic theoretical analysis to investigate the dynamics of the upper-hybrid mode instability driven by an initial ring electron distribution function, which is a form of loss-cone distribution. A favorable comparison is found between the two approaches, which shows that the reduced quasilinear theory, which is grounded in the concept of a model of the particle distribution function that is assumed to maintain a fixed mathematical form except that the macroscopic parameters that define the distribution are allowed to evolve in time, can be an effective tool in the study of plasma instabilities, especially if it is guided by and validated against the more rigorous simulation result.

Physics↗
DOE OSTI2023

Three-dimensional shear-flow instability saturation via stable modes

Turbulence in three dimensions (3D) supports vortex stretching that has long been known to accomplish energy transfer to small scales. Moreover, net energy transfer from large-scale, forced, unstable flow-gradients to smaller scales is achieved by gradient-flattening instability. Despite such enforcement of energy transfer to small scales, it is shown here that the shear-flow-instability-supplied 3D-fluctuation energy is largely inverse-transferred from the fluctuation to the mean-flow gradient, and such inverse transfer is more efficient for turbulent fluctuations in 3D than in two dimensions (2D). The transfer is due to linearly stable eigenmodes that are excited nonlinearly. The stable modes, thus, reduce both the nonlinear energy cascade to small scales and the viscous dissipation rate. The vortex-tube stretching is also suppressed. Up-gradient momentum transport by the stable modes counters the instability-driven down-gradient transport, which also is more effective in 3D than in 2D (≈70% vs ≈50%). From unstable modes, these stable modes nonlinearly receive energy via zero-frequency fluctuations that vary only in the direction orthogonal to the plane of 2D shear flow. The more widely occurring 3D turbulence is thus inherently different from the commonly studied 2D turbulence, despite both saturating via stable modes.

Mechanics↗
DOE OSTI2024

Analysis of fusion fuel conditions and instability development in shocklessly compressed MagLIF implosion simulations

Magnetized liner inertial fusion (MagLIF) implosions on the Z accelerator have almost exclusively been driven by ∼100-ns rise time current pulses. The rise time is selected to be as short as achievable on Z partially to minimize the time during which deleterious implosion instabilities can develop. Modifying the shape of the current pulse could provide benefits for MagLIF, including more efficient compression of the fusion fuel and the magnetic flux inside the liner cavity. Quasi-isentropic compression of the liner prevents formation of shocks in the liner material and reduces the amount of entropy generation within the liner. This allows for more final compression of the liner and fuel assembly. We present results from one-dimensional (1D) radiation-magnetohydrodynamic (rad-MHD) simulations comparing thermonuclear fuel conditions in MagLIF implosions driven with two different current pulses: a ∼100-ns rise time, ∼21.5 MA peak current “short pulse” and a ∼200-ns rise time, ∼21.5 MA peak current “shockless” pulse. We also quantify and compare the instability development in three-dimensional (3D) MHD implosion simulations driven by these two different pulse shapes. Our 1D simulations indicate that the shocklessly compressed MagLIF implosion performs better than the short pulse driven implosion with a >50% higher thermonuclear neutron yield, and 3D simulations indicate comparable implosion instability development, suggesting that pulse shaping could enable improvements to MagLIF performance on Z without compromising implosion stability.

Shipley, G. A. (ORCID:0000000205364001)↗
DOE OSTI2024

Investigating instabilities in magnetized low-pressure capacitively coupled RF plasma using particle-in-cell (PIC) simulations

The effect of a uniform magnetic field on particle transport in low-pressure radio frequency (RF) capacitively coupled plasma (CCP) has been studied using a particle-in-cell model. Three distinct regimes of plasma behavior can be identified as a function of the magnetic field. In the first regime at low magnetic fields, asymmetric plasma profiles are observed within the CCP chamber due to the effect of $\overrightarrow{E}$ × $\overrightarrow{B}$ drift. As the magnetic field increases, instabilities develop and form self-organized spoke-shaped structures that are distinctly seen within the bulk plasma closer to the sheath. In this second regime, the spoke-shaped coherent structures rotate inside the plasma chamber in the -$\overrightarrow{E}$ × $\overrightarrow{B}$ direction, where $\overrightarrow{E}$ × $\overrightarrow{B}$ are the DC electric and magnetic field vectors, respectively, and the DC electric field exists in the sheath and pre-sheath regions. The spoke rotation frequency is in the megahertz range. As the magnetic field strength increases further, the rotating coherent spokes continue to exist near the sheath. The coherent structures are, however, accompanied by new small-scale incoherent structures originating and moving within the bulk plasma region away from the sheath. This is the third regime of plasma behavior. The threshold values of the magnetic field between these regimes were found not to vary with changing plasma reactor geometry (e.g., area ratio between ground and powered electrodes) or the use of an external capacitor between the RF-powered electrode and the RF source. The threshold values of the magnetic field between these regimes shift toward higher values with increasing gas pressure. Analysis of the results indicates that the rotating structures are due to the lower hybrid instability driven by density gradients and electron-neutral collisions. This paper provides guidance on the upper limit of the magnetic field for instability-free operation in low-pressure CCP-based semiconductor deposition and etch systems that use the external magnetic field for plasma uniformity control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗