Search NASA⌕ Search

SEARCH · Search NASA

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

An experimental study of the existence regions and non-linear interactions of drift wave and Kelvin–Helmholtz instabilities in a linear magnetized plasma

Experimental observations of the intrinsic excitation and non-linear interactions of drift wave (DW) and Kelvin–Helmholtz (KH) instabilities in a linear magnetized plasma column are presented. The experiments are carried out in the inverse mirror plasma experimental device (IMPED)—a cylindrical, magnetized, linear plasma machine designed to study low-frequency waves and instabilities in plasma. A novel feature of IMPED is the ability to control plasma profiles, such as the density n(r)⁠, electron temperature T e (r)⁠, and plasma potential V p (r) by varying the ratio Rm of the magnetic field in the main chamber to that in the source chamber. At high values of Rm, higher-density gradient scale length promotes the drift wave (DW) instability while lower Rm value results in a higher radial electric field, inducing a sheared poloidal flow that enhances the dominance of the Kelvin–Helmholtz (KH) mode. The background and fluctuating plasma parameters are characterized using various configurations of multiple in situ electric probes at different spatial locations to quantify the local gradients that excite the low-frequency primary instabilities. Statistical, spectral, and bispectral analysis of the density and potential signals help identify these modes in terms of wave number, frequency, phase, and amplitude and also delineate the nature of their non-linear interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Shock Stabilization of the HLLC Riemann Solver for the Carbuncle Instability

The HLLC approximate Riemann solver improves upon the HLL Riemann solver by resolving contact discontinuities. This is a particularly desirable property for multi-material codes in which problems usually contain material interfaces. However, the HLLC solver is known to suffer from the carbuncle phenomenon, a numerical instability most apparent at grid-aligned shocks in multi-dimensional simulations. Many problems of interest, including high energy-density physics applications, require the accurate resolution of both material interfaces and hydrodynamic shocks. A variety of methods have been developed to cure this instability, with varying degrees of complexity. The objective of this work is to describe a simple approach to modify the HLLC Riemann solver and prevent the carbuncle instability. Here, the method is then demonstrated for assorted two-dimensional test problems known to exhibit the shock instability. The performance of the new solver is compared with that of the standard HLL and HLLC Riemann solvers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Local Lagrangian reduced-order modeling for the Rayleigh-Taylor instability by solution manifold decomposition

The Rayleigh-Taylor instability is a classical hydrodynamic instability of great interest in various disciplines of science and engineering, including astrophysics, atmospheric sciences and climate, geophysics, and fusion energy. Analytical methods cannot be applied to explain the long-time behavior of the Rayleigh-Taylor instability, and therefore, numerical simulation of the full problem is required. However, in order to capture the growth of amplitude of perturbations accurately, both the spatial and temporal discretizations need to be extremely fine for traditional numerical methods, and long-time simulation may become prohibitively expensive. In this paper, we propose efficient reduced order model techniques to accelerate the simulation of the Rayleigh-Taylor instability in compressible gas dynamics. Here, we introduce a general framework for decomposing the solution manifold to construct the temporal domain partition and temporally-local reduced order model construction with varying Atwood number. We propose two practical approaches in this framework, namely decomposition by physical time and by penetration distance. Numerical results are presented to examine the performance of the proposed approaches.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Structure formation in turbulence as an instability of effective quantum plasma

Structure formation in turbulence can be understood as an instability of “plasma” formed by fluctuations serving as effective particles. These “particles” are quantumlike in the sense that their wavelengths are non-negligible compared to the sizes of background coherent structures. The corresponding “kinetic equation” describes the Wigner matrix of the turbulent field, and the coherent structures serve as collective fields. This formalism is usually applied to manifestly quantumlike or scalar waves. Here, we show how to systematically extend it to more complex systems using compressible Navier–Stokes turbulence as an example. In this case, the fluctuation Hamiltonian is a five-dimensional matrix operator and diverse modulational modes are present. As an illustration, we calculate these modes for a sinusoidal shear flow and find two modulational instabilities. One of them is specific to supersonic flows, and the other one is a Kelvin–Helmholtz-type instability that is a generalization of the known zonostrophic instability. Our calculations are readily extendable to other types of turbulence, for example, magnetohydrodynamic turbulence in plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Kelvin-Helmholtz instability as one of the key features for fast and efficient emulsification by hydrodynamic cavitation

The paper investigates the oil–water emulsification process inside a micro-venturi channel. More specifically, the possible influence of Kelvin-Helmholtz instability on the emulsification process. High-speed visualizations were conducted inside a square venturi constriction with throat dimensions of 450 µm by 450 µm, both under visible light and X-Rays. We show that cavity shedding caused by the instability results in the formation of several cavity vortices. Their rotation causes the deformation of the oil stream into a distinct wave-like shape, combined with fragmentation into larger drops due to cavitation bubble collapse. Later on, the cavity collapse further disperses the larger drops into a finer emulsion. Thus, it turns out that the Kelvin-Helmholtz instability is similarly characteristic for hydrodynamic cavitation emulsification inside a microchannel as is the Rayleigh-Taylor instability for acoustically driven emulsion formation.

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

Ion-temperature- and density-gradient-driven instabilities and turbulence in Wendelstein 7-X close to the stability threshold

Electrostatic gyrokinetic instabilities and turbulence in the Wendelstein 7-X stellarator are studied. Particular attention is paid to the ion-temperature-gradient (ITG) instability and its character close to marginal stability [Floquet-type turbulence (Zocco et al., Phys. Rev. E, vol. 106, 2022, p. L013202) with no electron temperature gradient]. The flux tube version of the δƒ code stella (Barnes et al., J. Comput. Phys., vol. 391, 2019, pp. 365–380) is used to run linear and nonlinear gyrokinetic simulations with kinetic electrons. The nature of the dominant instability depends on the wavelength perpendicular to the magnetic field, and the results are conveniently displayed in stability diagrams that take this dependence into account. This approach highlights the presence of universal instabilities, which are less unstable but have longer wavelengths than other modes. A quasi-linear estimate of the heat flux suggests they are relevant for transport. Close to the stability threshold, the linear eigenmodes and turbulence form highly extended structures along the computational domain if the magnetic shear is small. Numerical experiments and diagnostics are undertaken to assess the resulting radial localisation of the turbulence, which affects the interaction of the latter with zonal flows. Increasing the amplitude of the magnetic shear (e.g.through current drive) has a stabilising effect on the turbulence and, thus, reduces the nonlinear energy transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Role of ion acoustic instability in magnetic reconnection

We report on a first-principles numerical study of magnetic reconnection in plasmas with different initial ion-to-electron temperature ratios. In cases where this ratio is significantly below unity, we observe intense wave activity in the diffusion region, driven by the ion-acoustic instability. Our analysis shows that the dominant macroscopic effect of this instability is to drive substantial ion heating. In contrast to earlier studies reporting significant anomalous resistivity, we find that anomalous contributions due to the ion-acoustic instability are minimal. These results shed light on the dynamical impact of this instability on reconnection processes, offering new insights into the fundamental physics governing collisionless reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Global Simulation of Electron Cyclotron Harmonic Wave Instability in a Storm-Time Magnetosphere

Electron cyclotron harmonic (ECH) waves are electrostatic emissions between the ECHs and play a dominant role for precipitating energetic electrons in the magnetotail. Statistically, the ECH wave intensity is stronger at nightside and dawnside than at dayside and duskside. In this study, we, for the first time, simulate the global ECH wave evolution during a geomagnetic storm event using Ring current Atmosphere interactions Model with Self-Consistent Magnetic field (RAM-SCB) combined with a linear growth rate solver.We find that the simulation results are generally consistent with the statistical and real-time observations. The ECH wave instability is much stronger at nightside and dawnside, compared to the instability at dayside and duskside. Before a geomagnetic storm (quiet time), the unstable regions of the ECH waves lie beyond $5R_E$ with a weak instability level. During the main phase of a geomagnetic storm, the unstable regions can extend to a lower altitude $(3R_E)$ with a strong instability level. Finally, during the recovery phase, the unstable regions return to $5$ ~ $6R_E$. We also find that the inner boundary of unstable ECH wave regions is coincident with the plasmapause location during the whole geomagnetic storm event.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Evaluation method of beam instability in laser ion source using solenoid

In a laser ion source using solenoid field confinement, it is known that an ion beam becomes unstable in the certain range of a magnetic field. Although it is essential to quantify the instability when discussing the unstable region of the beam, it is difficult to evaluate the beam instability by peak current or the amount of charge because an irregular change of temporal profile occurs in addition to the shot-by-shot fluctuation of amplitude. In this study, I propose the most appropriate method to evaluate the beam instability using the difference from an average waveform. The validity of the new method was evaluated by comparing three evaluation methods (variation of maximum value, variation of integral of waveforms, and the proposed evaluation method) with the experimentally obtained waveforms with the stable and unstable regions of a solenoid field. The proposed method was verified to best represent the beam instability by laser-induced plasma.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Generalized slab universal instability and its appearance in pair plasma

In this work, a generalized linear dispersion relation of electromagnetic slab universal modes is derived, taking into account arbitrary ion charge state, electron finite Larmor radius (FLR) effects, and Debye shielding effects. As a consequence, it provides more accurate predictions and is applicable to a wider range of plasmas. We find that electron FLR effects have a weakly stabilizing effect on the slab universal instability in electron–ion plasma, while Debye shielding strongly stabilizes this instability when λ D approaches ρ i (λ D is the Debye length and ρi is the ion gyroradius). In particular, we examine the stability criterion for this instability in electron–positron pair plasmas and find that the instability persists in this simplest plasma system as long as the pair plasma number density exceeds the critical value n c = B 2 /(8πm e c 2 ).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Thermomagnetic instability of plasma composition gradients

In this paper, we show that, under Braginskii magneto-hydrodynamics, anti-parallel gradients in an average ion charge state and electron temperature can be unstable to the growth of self-generated magnetic fields. The instability is analogous to the field-generating thermomagnetic instability, although it is driven by the collisional thermal force magnetic source term rather than the Biermann battery term. The gradient in ion charge state causes a gradient in collisionality, which couples with temperature perturbations to create a self-generated magnetic field. This magnetic field deflects the electron heat flux in a way that reinforces the temperature perturbation. The derived linearized growth rate, typically on hydrodynamic timescales, includes resistive and thermal smoothing. It increases with large ion composition gradients and electron heat flux, conditions typical of the hohlraum walls or contaminant mix jets in inertial confinement fusion implosions. However, extended magneto-hydrodynamic simulations indicate that the instability is usually dominated and stabilized by nonlinear Nernst advection, in a similar manner to the standard thermomagnetic instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Vortex-sheet modeling of hydrodynamic instabilities produced by an oblique shock interacting with a perturbed interface in the HED regime

In this work, we consider hydrodynamic instabilities produced by the interaction of an oblique shock with a perturbed material interface under high-energy-density (HED) conditions. During this interaction, a baroclinic torque is generated along the interface due to the misalignment between the density and pressure gradients, thus leading to perturbation growth. Our objective is to understand the competition between the impulsive acceleration due to the normal component of the shock velocity, which drives the Richtmyer–Meshkov instability, and the shear flow across the interface due to the tangential component of the shock velocity, which drives the Kelvin–Helmholtz instability, as well as its relation to perturbation growth. Since the vorticity resulting from the shock-interface interaction is confined to the interface, we describe the perturbation growth using a two-dimensional vortex-sheet model. We demonstrate the ability of the vortex-sheet model to reproduce roll-up dynamics for non-zero Atwood numbers by comparing to past laser-driven HED experiments. We determine the dependence of the interface dynamics on the tilt angle and propose a time scaling for the perturbation growth at early time. Eventually, this scaling will serve as a platform for the design of future experiments. This study is the first attempt to incorporate into a vortex-sheet model the time-dependent interface decompression and the deceleration (as well as the corresponding Rayleigh–Taylor instability) arising from laser turn-off.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Simulations of the shock-driven Kelvin–Helmholtz instability in inclined gas curtains

In this paper, we present simulation results for the two-dimensional, shock-driven Kelvin–Helmholtz instability. Simulations are performed with a Mach 2.0 shock propagating through a finite-thickness curtain of gas inclined at an angle α0=30° with respect to the shock plane. After the passage of the shock, the gas curtain is accelerated along its axis. A perturbation develops due to shock reflection near the lower wall, and a Kelvin–Helmholtz instability forms near the vertical center of the curtain. This is the first known numerical reproduction of these phenomena that have previously been observed in experiments with an inclined cylindrical gas column. The effects of varying Mach number and column width were explored in detail to complement experimental data. Additionally, the dependence of the Kelvin–Helmholtz wavelength on Mach number closely matches the relationship observed in experiments. This supports the notion that the observed instability is effectively two-dimensional and inviscid (like classical Kelvin–Helmholtz). The growth rate of the perturbations in the gas curtain was also found to be similar for different Mach numbers. The perturbation at the curtain foot, previously unreported in experiments, was found to have a similar relationship to Mach number as the Kelvin–Helmholtz instability. Both perturbation wavelengths are found to be proportional to layer width. Simulations were performed with the fast interfaces and transport in the atmosphere, an exascale ready, graphics processing unit-accelerated compressible flow solver developed at the University of New Mexico.

42 ENGINEERING↗
DOE OSTI2021

Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in planar geometry are derived in the small Atwood number (Boussinesq) limit. The models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar (heavy-fluid mass fraction) variance S and its dissipation rate Χ modeled either differentially or algebraically. The models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing, respectively. Mixing layer growth parameters and other physical observables relevant to each instability are obtained explicitly as functions of the model coefficients. The turbulent fields are also expressed in terms of the model coefficients, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. The model calibration methodology is described and discussed. Expressions for a subset of the various physical observables are used to calibrate each of the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to these values of the observables and to specific canonical Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz turbulent flows. A calibrated four-equation model is then used to reconstruct the mean and turbulent fields, and late-time turbulent equation budgets for each instability-induced flow across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗
DOE OSTI2024

Onset of kinetic effects on Rayleigh–Taylor instability: Advective–diffusive asymmetry

In nature and engineering applications, the Rayleigh–Taylor instability (RTI) occurs over a wide range of Atwood, Reynolds, Mach, and Knudsen numbers. At low Atwood, Mach, and Knudsen numbers, the classic advective instability causes quasi-symmetric bubble and spike growth on the two sides of the interface. However, recent findings suggest that at high degrees of rarefaction, advective effects are suppressed and molecular diffusion leads to planar growth of the density fronts on either side of the interface. This study aims to investigate the flow physics of the transition from advective to diffusive behavior, focusing on the onset of the kinetic effects. Using the gas kinetic methodology, RTI is simulated over a range of Knudsen and Mach numbers in the transition regime. The simulation results reveal the various stages of transformation from advective instability to diffusive transport. For the first time, the study demonstrates the existence of a Knudsen–Mach parameter regime where the bubble side exhibits advective instability, while the other side shows a planar density front due to molecular diffusion, rather than the canonical advective spike shape. The dominance of different mechanisms on the two sides of the interface leads to the advective–diffusive asymmetry. In conclusion, the findings of this study can lead to a more comprehensive understanding of RTI over a wide range of Mach and Knudsen numbers.

42 ENGINEERING↗
DOE OSTI2025

A future of inertial confinement fusion without laser-plasma instabilities

From the beginning of inertial confinement fusion (ICF) research, laser-plasma instabilities excited by narrowband lasers have limited the hydrodynamic design space of all laser-based approaches to inertial fusion energy (IFE). With advances in broadband laser technologies, the next generation of ICF drivers will likely have large bandwidth that is engineered to mitigate laser-plasma instabilities, thereby expanding the hydrodynamic design space to include both robust high yields (>200−MJ) with large-energy laser systems (>4 MJ) and high gains (>10) with moderate-energy laser facilities (<2 MJ). State-of-the-art simulations indicate that laser bandwidths of a few percent are required to mitigate instabilities for IFE-relevant conditions. To test these models and demonstrate that high-bandwidth lasers mitigate laser-plasma instabilities, the Fourth-generation Laser for Ultrabroadband eXperiments will be used with the OMEGA Laser System. The goal is to provide the community with the confidence to invest in a multiple-beam high-bandwidth laser facility that will demonstrate the necessary ablation pressures for robust direct-drive ignition without detrimental levels of hot electrons that degrade the implosion performance. It is important to recognize that the highest performing ICF implosions will likely never be completely LPI free because of the significant advantages to maximizing the laser intensity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Cooling and instabilities in colliding flows

ABSTRACT Collisional self-interactions occurring in protostellar jets give rise to strong shocks, the structure of which can be affected by radiative cooling within the flow. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in both one and three dimensions with a power-law cooling function. The characteristic length and time-scales for cooling are temperature dependent and thus may vary as shocked gas cools. When the cooling length decreases sufficiently and rapidly, the system becomes unstable to the radiative shock instability, which produces oscillations in the position of the shock front; these oscillations can be seen in both the one- and three-dimensional cases. Our simulations show no evidence of the density clumping characteristic of a thermal instability, even when the cooling function meets the expected criteria. In the three-dimensional case, the nonlinear thin shell instability (NTSI) is found to dominate when the cooling length is sufficiently small. When the flows are subjected to the radiative shock instability, oscillations in the size of the cooling region allow NTSI to occur at larger cooling lengths, though larger cooling lengths delay the onset of NTSI by increasing the oscillation period.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Interplay between the non-resonant streaming instability and self-generated pressure anisotropies

The non-thermal particles escaping from collision-less shocks into the surrounding medium can trigger a non-resonant streaming instability that converts parts of their drift kinetic energy into large amplitude magnetic field perturbations, and promote the confinement and acceleration of high energy cosmic rays. We present simulations of the instability using an hybrid-Particle-in-Cell approach including Monte Carlo collisions, and demonstrate that the development of the non-resonant mode is associated with important ion pressure anisotropies in the background plasma. Depending on the initial conditions, the anisotropies may act on the instability by lowering its growth and trigger secondary micro-instabilities. Introducing collisions with neutrals yield a strong reduction of the magnetic field amplification as predicted by linear fluid theory. In contrast, Coulomb collisions in fully ionized plasmas are found to mitigate the self-generated pressure anisotropies and promote the growth of the magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗