Search NASA⌕ Search

SEARCH · Search NASA

Results for “instability”

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

Peeling-ballooning modes in spherical tokamaks: Multi-branch instabilities and effects beyond ideal MHD

A number of important physics effects on the stability of relatively high-n (n is the toroidal mode number) peeling-ballooning modes (PBMs) are investigated based on an equilibrium reconstructed from a NSTX discharge, utilizing extended magnetohydrodynamic (MHD) eigenvalue solvers. For a given toroidal mode number n, multiple branches of instabilities are computed, with the total number of unstable branches roughly linearly scaling with n. Most of the unstable branches are located in the plasma core region, but edge-localized branches, i.e., PBMs, are also identified at higher n-numbers. For the single-fluid-wise most unstable PBM with n = 19⁠, stabilizing/destabilizing effects due to various physics beyond ideal MHD are systematically investigated. Plasma toroidal flow is found to be weakly stabilizing. Local flow shear is generally stabilizing as well, with the degree of stabilization depending on the initial growth rate (without flow shear) of the mode. The plasma resistivity can strongly destabilize the PBM within the single-fluid framework. Anisotropic thermal transport, strong parallel sound wave damping, as well as two-fluid effects are all stabilizing to the mode. In particular, diamagnetic stabilization (within the two-fluid model) is found to be very strong for this mode.

Linear stability analysis↗
DOE OSTI2026

Effects of external magnetic fields on the multi-mode Rayleigh–Taylor instability

The magneto-Rayleigh–Taylor instability (mRTI) is a key process in inertial confinement fusion and is thought to be widespread in the interstellar medium, where it can concentrate plasma into discrete structures. We present resistive MHD simulations of the nonlinear evolution of multi-mode mRTI in both two and three dimensions, examining the effects of uniform external magnetic fields oriented either parallel or perpendicular to the initial interface. In both 2D and 3D, weak parallel fields enhance mixing-zone growth, whereas stronger fields suppress it. For perpendicular fields, growth is initially inhibited but becomes enhanced at later times. These behaviors arise from magnetic tension, which modifies flow anisotropy, buoyancy, drag, and vortex dynamics. The interplay of these mechanisms governs the distinct ways in which magnetic fields influence mRTI evolution.

Flow instabilities↗
DOE OSTI2025

Dependence of ion-cyclotron range of frequencies instabilities on species mix and fast-ion distribution: II. Time evolution

Dedicated DIII-D experiments studied instabilities below the ion cyclotron frequency f ci and at cyclotron harmonics by injecting thirteen different neutral beam populations into plasmas with different thermal compositions of H and D at five different values of toroidal field. The resulting database of 2529 conditions is described in Part 1: Stability. This paper focuses on the temporal evolution. For some beam geometries, the ion cyclotron emission (ICE) is much stronger at the beginning of injection, suggesting the importance of ‘bump-on-tail’ energy drive. Sawtooth crashes cause the emission to cease for some beam types and emission bands, but the opposite effect occurs for others. When unstable, for virtually all conditions and frequency bands, the emission bursts in a cycle reminiscent of predator-prey oscillations. Evidence for nonlinear interaction between ICE harmonics is observed. Strongly unstable modes typically appear within 1 ms from the start of beam injection and often decay within 0.1 ms when the beam turns off, but both the initial rise and final decay of the mode power exhibit a wide variety of behavior.

electromagnetic ion cyclotron instabilities↗
DOE OSTI2025

Electron-ion collider transverse instabilities due to the crab cavity fundamental impedance

The electron-ion collider crab Cavities will present a significant transverse impedance to the beam. The low-level radio frequency system will have to reduce the crab cavity impedance to prevent transverse instabilities, while regulating the crabbing voltage and minimizing the radio frequency noise levels injected into the beam. This work presents an estimate of the crab cavity impedance in the presence of feedback and the resulting stability margins using simplified stability criteria. Generalized stability criteria are also presented in this work and are used to more accurately estimate stability margins.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Evolution of radiative thermal instability in a confined medium

Thermally bistable fluid tends to self-organize into clouds of hot and cold material, which are internally uniform and separated by thin conduction fronts. The evolution of these clouds has been studied for isobaric systems, but when pressure is instead treated as a dynamical quantity and allowed to evolve self-consistently, fundamentally different dynamics appear. Such a treatment is necessary in some laboratory plasmas, whose volume is constrained but whose pressure can vary. Here, solutions are derived for the evolution of clouds, accounting for pressure variation and interactions between conduction fronts. Additional stable configurations and secondary instabilities are derived, which may be relevant to fusion plasmas and to the study of photoionized plasma in the laboratory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Electromagnetic Inhomogeneous-Energy-Density-Driven Instability in the UCLA Large Plasma Device (EM-IEDDI in LAPD) (Final Technical Report)

Using the Large Plasma Device (LAPD), the electromagnetic signatures associated with the inhomogeneous-energy-density-driven (IEDD) instability were targeted for documentation in a lab experiment. Exploitable features of LAPD are the extreme-length, high-beta, low-collisionality, and excellent reproducibility. The signatures of interest were mode amplitude, mode structure, frequency spectrum, phase velocity, and excitation threshold. The abrupt inhomogeneity at the interface between two coaxial plasma columns, that have disparate space potential values, were expected to generate a strong, localized, transverse, dc electric field and form a strong ExB-flow inhomogeneity that was expected to trigger broadband, intermediate-frequency waves that propagate many gyroradii beyond the interfacial source region. WVU was unable to destabilize the IEDDI mode in LAPD. Experimental limitation of three different parameters conspire to prevent IEDDI-favorable conditions even though Alfven-favorable conditions are met. Electron-neutral collision frequency is too high, Structuring of the radial profile of the radial electric field is negligible and the control of that structure is minimal, and the Magnetic field is too large for some purposes and too small for other purposes. WVU designed a new antenna and relied on externally-driving the IEDDI (instead of spontaneous destabilization). The antenna operation failed for unforeseen reasons and lack of electrician time and availability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Collaborative Research: Parametric Instabilities of Alfven Waves in Low‐beta Plasmas

Alfven waves are of fundamental importance in magnetized plasmas. This project aims to advance our understanding of the basic physics of the parametric instabilities of Alfven waves using advanced computer simulations. The results will have implications for important processes in the heliosphere involving Alfven waves, such as coronal heating and solar wind heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Global Simulations of Gravitational Instability in Protostellar Disks with Full Radiation Transport. I. Stochastic Fragmentation with Optical-depth-dependent Rate and Universal Fragment Mass

Fragmentation in a gravitationally unstable accretion disk can be an important pathway for forming stellar/planetary companions. To characterize quantitatively the condition for and outcome of fragmentation under realistic thermodynamics, we perform global 3D simulations of gravitationally unstable disks at various cooling rates and cooling types, including the first global simulations of gravitational instability that employ full radiation transport. We find that fragmentation is a stochastic process, with the fragment generation rate per disk area p frag showing an exponential dependence on the parameter β ≡ Ω K t cool , where Ω K is the Keplerian rotation frequency, and t cool is the average cooling timescale. Compared to a prescribed constant β, radiative cooling in the optically thin/thick regime makes p frag decrease slower/faster in β; the critical β corresponding to ∼1 fragment per orbit is ≈3, 5, and 2 for constant β, optically thin, and optically thick cooling, respectively. The distribution function of the initial fragment mass is remarkably insensitive to disk thermodynamics. Regardless of cooling rate and optical depth, the typical initial fragment mass is m frag ≈ 40M tot h 3 , with M tot being the total (star+disk) mass and h = H/R being the disk aspect ratio. Applying this result to typical Class 0/I protostellar disks, we find m frag ∼ 20M J , suggesting that fragmentation more likely forms brown dwarfs. Given the finite width of the m frag distribution, forming massive planets is also possible.

gravitational instability↗
DOE OSTI2025

Global Simulations of Gravitational Instability in Protostellar Disks with Full Radiation Transport. II. Locality of Gravitoturbulence, Clumpy Spirals, and Implications for Observable Substructure

Spiral perturbations in a gravitationally unstable accretion disk regulate disk evolution through angular momentum transport and heating and provide an observational signature of gravitational instability (GI). We use global 3D simulations to systematically characterize and understand these spiral perturbations. The spiral perturbations and the resulting transport are overall insensitive to the cooling type, with the exception that radiative cooling, especially in the optically thick regime, reduces the amplitude of temperature perturbations. Spiral perturbations are localized around corotation, allowing transport to be approximated by a local α viscosity to zeroth order in aspect ratio (H/R), but only after averaging over multiple orbits in time and/or multiple scale heights in space. Meanwhile, large-amplitude perturbations from strong gravitoturbulence can cause $\mathcal{O}$(α 1/2 ) deviation in the cooling rate of the disk. We develop empirical prescriptions for the angular momentum transport, heating, and cooling in a gravitoturbulent disk that capture the deviation from a viscous, unperturbed disk to first order in H/R and α 1/2 . The spiral perturbations in saturated gravitoturbulence are clumpy, with dense clumps forming through the nonlinear coupling between multiple modes at different m. Observationally, the clumpy gravitoturbulence produced by saturated GI can be mistaken with observational noise or embedded companions, especially under finite resolution. Meanwhile, grand-design spirals with m-fold symmetry may be uncommon among disks in saturated gravitoturbulence, and we speculate that they may instead be a signature of recently triggered or decaying GI.

Gravitational instability↗
DOE OSTI2025

Mechanistic Insights into the Surface Instabilities of TiNb 2 O 7, a High‐Power Li‐Ion Anode

TiNb 2 O 7 (TNO) is a promising Li-ion battery anode for high-power applications, such as implantable medical devices and heavy-duty equipment. Hailed as being safe due to its elevated operating potential near 1.6 V, TNO has long been assumed to be highly stable in the carbonate-based electrolytes used in Li-ion batteries. Herein, all mechanisms occurring at the surface of both TNO and Nd-doped TNO are identified, and both materials in fact show significant gassing. CO 2 is even released at open circuit conditions, demonstrating the poor chemical stability of the material in the electrolyte even prior to battery operation. Such extreme instability is a critical safety concern. In addition, it was found that Ti dissolves from the surface of TNO particles at low voltage (below 1.4 V vs Li), and in fact deposits on the counter electrode. Ti further inside TNO particles then diffuses to the Ti-poor surface during discharge. Partial carbon-coating as a mitigating measure has also been tested and found to exacerbate these processes. The findings identify novel reactions occurring within TNO, and clearly highlight the need to stabilize the surfaces of TNO in order to prevent such aggressive deterioration at the surface of the particles.

25 ENERGY STORAGE↗
DOE OSTI2024

Heterogeneous Solid Electrolyte Interphase Interactions Dictate Interface Instability in Sodium Metal Electrodes

Sodium (Na) metal batteries have attracted recent attention due to their low cost and high abundance of Na. However, the advancement of Na metal batteries is impeded due to key challenges such as dendrite growth, solid electrolyte interphase (SEI) fracture, and low Coulombic efficiency. This study examines the coupled electro‐chemo‐mechanical interactions governing the electrodeposition stability and morphological evolution at the Na/electrolyte interface. The SEI heterogeneities influence transport and reaction kinetics leading to the formation of current and stress hotspots during Na plating. Further, it is demonstrated that the heterogeneity‐induced Na metal evolution and its influence on the stress distribution critically affect the mechanical overpotential, contributing to a faster SEI failure. The analysis reveals three distinct failure mechanisms—mechanical, transport, and kinetic—that govern the onset of instabilities at the interface. Finally, a comprehensive comparative study of SEI failure in Na and lithium (Li) metal anodes illustrates that the electrochemical and mechanical characteristics of the SEI are crucial in tailoring the anode morphology and interface stability. This work delineates mechanistic stability regimes cognizant of the SEI attributes and underlying failure modes and offers important guidelines for the design of artificial SEI layers for stable Na metal electrodes.

25 ENERGY STORAGE↗
DOE OSTI2024

The radiation instability of thermally stable nanocrystalline platinum gold

Here, recent experimentally validated alloy design theories have demonstrated nanocrystalline binary alloys that are stable against thermally induced grain growth. An open question is whether such thermal stability also translates to stability under irradiation. In this study, we investigate the response to heavy ion irradiation of a nanocrystalline platinum gold alloy that is known to be thermally stable from previous studies. Heavy ion irradiation was conducted at both room temperature and elevated temperatures on films of nanocrystalline platinum and platinum gold. Using scanning/transmission electron microscopy equipped with energy-dispersive spectroscopy and automated crystallographic orientation mapping, we observe substantial grain growth in the irradiated area compared to the controlled area beyond the range of heavy ions, as well as compositional redistribution under these conditions, and discuss mechanisms underpinning this instability. These findings highlight that grain boundary stability against one external stimulus, such as heat, does not always translate into grain boundary stability under other stimuli, such as displacement damage.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Further Development of the Tamped Richtmyer-Meshkov Instability Method and Application to Molybdenum Dynamic Strength Calibration and Tabulation

The high pressure and high strain rate dynamic strength of Mo is experimentally and computationally investigated in the 3–20 GPa stress, 50–600 C temperature, and 10 5 –10 6 /s strain rate regimes using a modified tamped Richtmyer-Meshkov instability (RMI) method. Modifications to the tamped RMI method include a method to determine loading states during strain, a new strength calibration function based on interface shape, and a robust uncertainty quantification method. These modifications improve fidelity of the tamped RMI method, allowing evaluation of the compensating effects of pressure hardening, strain rate hardening, strain hardening, and thermal softening. The new calibration function based on interface shape is not limited to sinusoidal corrugations and could be applied to additional interface shapes. Plate impact experiments are performed at Argonne National Laboratory’s Advanced Photon Source’s Dynamic Compression Sector operated by Washington State University (DCS), driving a planar shock front through a corrugated Mo-D 2 O or Mo-C 8 F 18 interface, forcing the corrugation to significantly deform. The extent of interfacial deformation, RMI growth, is experimentally observed using X-ray phase contrast imaging at the DCS. RMI jet lengths and jet shapes are extracted from the experimental radiographs, then used to calibrate numerical simulations performed with the Sandia National Laboratories (SNL) hydrocode CTH. Mo yield strength, Y, as a function of shock pressure, P, strain rate, $\dot{\varepsilon }$, accumulated strain, ϵ, relative volumetric compression, RD , and temperature, T , is determined for each impact experiment and presented. The calibrated Mo yield strength values range 1.2–1.8 GPa, with strength generally decreasing as the impact stress increases. This trend is likely caused by thermal softening or strain localization. The tabular yield strength versus loading condition data presented in this paper can be used to fit complex strength models.

Voorhees, T. J. [Sandia National Lab. (SNL-CA), Li↗
DOE OSTI2026

Oxide nucleation via threshold-driven volume instability in Ni oxidation

Oxide nucleation dictates the onset of metal oxidation, yet the atomistic pathways bridging initial oxygen adsorption to bulk phase transformation remain elusive. Here, we utilize in-situ atomic-resolution imaging to directly capture the nucleation of nickel oxide, revealing a threshold-driven, cooperative transformation of Ni into NiO. Our observations identify a distinct incubation period during which oxygen progressively accumulates in subsurface layers. Once a critical concentration and penetration depth are reached, a collective lattice reconfiguration is triggered, abruptly converting multiple Ni layers into NiO. Atomistic simulations corroborate this mechanism, identifying a cooperative lattice instability induced by subsurface oxygen saturation. These results establish subsurface oxygen incorporation as the missing mechanistic link between surface adsorption and three-dimensional oxide formation, providing an atomistic framework to understand and control reactive phase transformations in materials synthesis, catalysis, and degradation.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Experiments and simulations on the Richtmyer-Meshkov instability with a thin intermediate layer

Experiments and simulations of the Richtmyer-Meshkov instability (RMI) in two- and three-layer configurations are presented. The two-layer case utilizes a light-over-heavy configuration and consists of air as the light gas and sulfur hexafluoride (SF 6 ) as the heavy gas. The three-layer case utilizes a light-intermediate-heavy configuration, with helium (He) as the light gas, air as the intermediate gas, and SF 6 as the heavy gas. Statistically significant differences in the mixing layer width of the lower interface are not observed between the two cases. This differs from the experiments of Schalles et al., where a small, though statistically significant, difference in mixing layer growth was observed between the two- and three-layer cases with a nominally two-dimensional, single mode perturbation. Notably, the perturbations on the lower interface in the present work do not grow large enough to significantly interact with the upper interface during the duration of the experiments. In conclusion, this suggests that the differences in mixing layer growth observed by Schalles et al. may be due to interactions of the perturbations on one interface with the other interface rather than being inherent to the three-layer problem.

42 ENGINEERING↗
DOE OSTI2024

The influence of the shock-to-reshock time on the Richtmyer–Meshkov instability in reshock

Experiments on the Richtmyer–Meshkov instability (RMI) in a dual driver vertical shock tube (DDVST) are described. An initially planar, stably stratified membraneless interface is formed by flowing air from above and sulfur hexafluoride from below the interface location using the method of Jones & Jacobs ( Phys. Fluids , vol. 9, issue 1997, 1997, pp. 3078–3085). A random three-dimensional, multi-modal initial perturbation is imposed by vertically oscillating the gas column to produce Faraday waves. The DDVST design generates two shock waves, one originating above and one below the interface, with these shocks having independently controllable strengths and interface arrival times. The shock waves have nominal strengths of $M_L=1.17$ and $M_H=1.18$ for the shock wave originating in the light and heavy gas, respectively, with these strengths chosen to result in arrested bulk interface motion following reshock. The influence of the length of the shock-to-reshock time, as well as the order of shock arrival, on the post-reshock RMI is examined. The mixing layer width grows according to $h\propto t^\theta$ , where $\theta _H=0.36\pm 0.018$ (95 %) and $\theta _L=0.38\pm 0.02$ (95 %) for heavy and light shock first experiments, respectively, indicating no strong dependence on the order of shock wave arrival. Volume integrated specific turbulent kinetic energy (TKE) in the mixing layer versus time is found to decay according to $E_{tot}/\bar {\rho }\propto t^p$ with $p_H=-0.823\pm 0.06$ (95 %) and $p_L=-1.061\pm 0.032$ (95 %) for heavy and light shock first experiments, respectively. Notably, the 95 % confidence intervals do not overlap. Analysis on the influence of the shock-to-reshock time on turbulent length scales, transition criteria, spectra and mixing layer anisotropy are also presented.

Ferguson, Kevin (ORCID:0000000166537482)↗
DOE OSTI2026

Pressure-Induced Chemical Bonding Effects on Lattice and Magnetic Instabilities in Antiferromagnetic Insulating CaMn 2 Sb 2

Exotic quantum phenomena often emerge near an electronic delocalization transition (EDT) from an antiferromagnetic insulating phase to a strongly correlated metallic state under pressure. We report the pressure-induced structural and magnetic evolution of the antiferromagnetic insulator CaMn 2 Sb 2 . Single-crystal X-ray diffraction reveals a first-order phase transition near 5.4 GPa from a trigonal P-3m1 structure to a monoclinic P2 1 /m phase accompanied by a ∼7% volume collapse. Residual electron density analysis at intermediate pressures reveals charge localization along Mn–Sb chains, signaling electronic instability preceding the structural transition. Bonding analysis indicates anisotropic Mn–Sb orbital reconfiguration under pressure, driving a distorted square-pyramidal geometry. Neutron scattering confirms the transition and identifies a pressure-induced incommensurate magnetic order, distinct from the ambient antiferromagnetic state. In the monoclinic phase, zigzag Mn chains exhibit antiferromagnetic coupling along the ac-plane, enabled by enhanced orbital overlap. These results establish CaMn 2 Sb 2 as a model system for studying the coupling of structural distortion, charge redistribution, and magnetic order in layered Mn pnictides under pressure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗