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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Bubble explosion induced melt pool instabilities in electron beam melting of aluminum alloy

Electron beam melting (EBM) is an additive manufacturing technology that can process materials and manufacture components otherwise impossible or uneconomical. However, defects, including porosity and surface irregularities, are widely reported in EBM-built components, and their formation mechanisms are not fully understood. Here, using in-situ high-speed synchrotron X-ray imaging, we reveal that bubble explosions in Al6061 during EBM induce melt pool instabilities contributing to defect formation. The melt pool and keyhole evolve through three stages: (1) initial formation of a melt pool, (2) subsurface bubble formation and explosion, and (3) periodic keyhole oscillation. During scanning, periodic bubble explosions can eject molten liquid as spatters and disturb the vapor depression and melt pool, contributing to surface humping, that may trigger lack-of-fusion defects in subsequent layers. The physical insights we report could provide guidance for EBM machine development, process innovation, alloy design and model development.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Phase instability-coupled fracture behavior in garnet LLZO solid electrolytes: a machine learning-enabled atomistic study

Fracture in the garnet-type solid electrolyte Li 7 La 3 Zr 2 O 12 (LLZO) poses a critical threat to both the performance and safety of solid-state batteries. To unravel the coupled chemomechanical processes that govern fracture in LLZO under external loading, we carry out large-scale molecular-dynamics simulations with a validated machine-learning force field. Our results demonstrate that triaxial stresses at crack flanks trigger a localized cubic-to-tetragonal phase transformation, which is accompanied by Li-ion rearrangement. The emergent tetragonal domains feature lattice contraction normal to the fracture plane, imposing coherent misfit strains that provide an additional driving force for further crack propagation. Crucially, introducing Li deficiencies stabilizes the cubic phase, postponing the phase transition and thereby delaying fracture initiation. These findings highlight the role of intrinsic phase instability in dictating LLZO's fracture resistance and its critical connection to local Li concentration. This chemomechanical coupling points toward targeted strategies to enhance the mechanical robustness of garnet electrolytes, including tuning Li content, ensuring dopant homogeneity, and refining processing protocols.

Materials science↗
DOE OSTI2024

Eigenmode twist by energy flux during fast-ion driven instabilities

Twisting of the spatial structure of eigenmodes by energy fluxes generated due to the spatial channeling (SC)—the fast-ion energy and momentum transfer across the magnetic field by destabilized modes—is considered. It is revealed that there is a correlation between the direction of the radial wave energy flux, orientation of the mode twist (MT), and the direction of mode rotation. The mode twist parameter is introduced and relations connecting it with the energy flux are established. It is shown the energy flux transforms zeros of the radial profile of the mode amplitude into minima (i.e., zeros disappear). It is found that the radial group velocity and phase velocity of reversed shear Alfvén eigenmodes (RSAEs) have opposite directions. These findings can be used for diagnostics of the energy fluxes during fast-ion driven instabilities and understanding whether the SC degrades or improves plasma performance. Specific calculations are carried out for fast magnetoacoustic modes, FMMs (known also as compressional Alfvén eigenmodes, CAEs), and Alfvén eigenmodes, AEs. A DIII-D experiment where RSAEs were observed is considered. It is concluded that an outward SC took place in this experiment. Peculiarities of various modes are discussed, which may explain why twisting of AEs, but not of FMMs, was observed experimentally in the DIII-D and NSTX tokamaks.

Kolesnichenko, Ya. I. (ORCID:0000000323219681)↗
DOE OSTI2025

Shock tube simulations for the three-layer Richtmyer–Meshkov instability with single-mode and multimode perturbations

While the canonical two-component, single-mode Richtmyer–Meshkov instability (RMI) has been extensively studied, relatively less work has focused on the effects of an additional intermediate-density middle layer. This work investigates such three-material RMI configurations at two Atwood number scenarios using the ARES hydrodynamics code. After validation against previous experimental and computational studies, setups corresponding to recent three-layer shock tube experiments are simulated. Cases with both single-mode and multimode perturbations are studied to quantify mixing across the interface between the materials with highest and intermediate density. In particular, this work is able to comprehensibly examine differences between two- and three-dimensional setups for the single-mode and multimode problems. Observations from previous two-layer investigations still apply in the three-layer setup, but over the time horizons considered, there appears to be insufficient nonlinear mode coupling to create significant differences between two- and three-dimensional simulations following the first passage of a shock. Finally, additional reshock simulations have additional nonlinear growth that does result in expected differences between two- and three-dimensional cases in this three-layer setup, but significant differences do not manifest during the time horizon studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Axisymmetric hydrodynamics in numerical relativity: treating coordinate singularity, artificial heating and modeling MHD instabilities

Two-dimensional axisymmetric simulations of binary neutron star (BNS) merger remnant are a cheap alternative to 3D simulations. To maintain realism for secular timescales, simulations must avoid accumulated errors from drifts in conserved quantities and artificial heating, and they must model turbulent transport in a way that remains plausible throughout the evolution. It is also crucial to avoid numerical artifacts due to the polar coordinate axis singularity. Methods that behave well near the axis often break flux-conservative form of the hydrodynamic equations, resulting in significant drifts in conserved quantities. We present a flux-conservative scheme that maintains smoothness near the axis without sacrificing conservative formulation of the equations or incurring drifts in conserved global quantities. We compare the numerical performance of different treatments of the hydrodynamic equations when evolving a hypermassive neutron star resembling the remnant of a BNS merger. These simulations demonstrate that the new scheme combines the axis smoothness of non-conservative methods with the mass and angular momentum conservation of other conservative methods on $\sim 10^2$ ms timescales of viscous and neutrino-driven evolution. Because fluid profiles remain smooth in the remnant interior, it is possible to remove artificial heating by evolving the entropy density. We show how physical heating and cooling terms can be easily calculated from source terms of the conservative evolution variables and demonstrate our implementation. Finally, we discuss and implement improvements to the effective viscosity scheme to better model the effect of magnetohydrodynamic instabilities as the remnant evolves.

axisymmetry↗
DOE OSTI2025

First observations of edge instabilities in strongly shaped negative triangularity plasmas on DIII-D

Throughout the 2023 DIII-D negative triangularity (NT) campaign, unexpected magnetic fluctuations indicative of an edge localized instability were observed in the majority of strongly shaped (–0.55 ≤ δ ≤ –0.45) discharges. Moreover, these same fluctuations were found to be absent in weaker NT plasma shapes. Based on a database study, three separate mode behaviors are observed based primarily on the edge safety factor and normalized plasma beta β N . The main regime is the bursty state, in which short (~5 ms) magnetic bursts occur every 5–10 ms. Each burst drives a small reduction in the edge electron temperature (T e ) as measured by electron cyclotron emissions (ECE), which then rapidly recovers between the bursts. The quiescent regime appears to have governing mechanics similar to the bursty regime, but with multiple unique mode signatures observed in the magnetics. These additional modes greatly reduce the bursty nature of the plasma. Finally, a saturated state is observed in a number of discharges, particularly at early time slices. The saturated regime is often short-lived (maximum 500 ms), but can have a large impact on the dynamics of the plasma edge region, causing significant drops in electron temperature as observed in the ECE diagnostic. 2D ECE imaging (ECE-I) is utilized, showing large kink-like perturbations localized near the last closed flux surface that extend into the SOL indicative of an external kink/peeling behavior. Additionally, evidence of phase mixing and observed tearing-parity perturbations localized near the edge rational surfaces is observed in the quiescent regime plasmas. As a result, the various modes do not appear to have a large impact on global confinement but may play a role in regulating the transport of the NT edge and the observed edge gradients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Infernal instabilities in negative-triangularity plasmas with negative central shear

Abstract A systematic numerical investigation is carried out to understand magnetohydrodynamic stability of the ideal infernal-kink instability in tokamak plasmas with both negative triangularity (neg-D) shaping and negative central shear for the equilibrium safety factor profile. The latter is motivated by the desire to form the internal transport barrier in the neg-D configuration, which is known to have difficulty in forming the edge transport barrier. The infernal-kink mode is generally found to be more unstable in neg-D plasmas as compared to their positive D-shaped (pos-D) counterpart. This is mainly due to less favorable (or even unfavorable) average magnetic curvature near the radial location of the minimum safety factor ( q min ) as compared to the pos-D configuration. The larger Shafranov shift associated with the neg-D shape helps the mode stabilization but is not sufficient to overcome the destabilizing effect due to bad curvature. Strong poloidal mode coupling due to plasma shaping (toroidicity, elongation, triangularity, etc.) helps explain the slight shift with respect to that predicted by the analytic theory of the peak location of the computed mode growth versus q min .

Zhang, X. M.↗
DOE OSTI2025

Shafranov-shift destabilization of ballooning-type micro-instabilities

Contrary to common wisdom, we argue that the Shafranov shift is destabilizing for the ambient ballooning-type instabilities, which account for most of the transport flux in tokamak core plasmas. Higher Shafranov shift indeed reduces the magnetic drift frequency ω d , i.e. improves the bad curvature mildly, around the outboard midplane. However, this improvement is limited to a very narrow spatial region while ω d is increased over the remaining poloidal space. The eigenfunction averaged $\langle$ω d $\rangle$ is effectively enhanced due to the finite mode width of the eigenfunction and is thus destabilizing to the ballooning-type mode, as demonstrated by gyrokinetic simulations with the CGYRO code using local Miller equilibrium geometry. The predicted nonlinear flux also increases with Shafranov shift, consistent with linear simulations. The reduced transport model TGLF can capture the physics reasonably well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Measurements of the polarization of several instabilities in the DIII-D tokamak

Recently, a method to infer the polarization of modes with frequencies much less than the ion cyclotron frequency was published [X.D. Du et al., Phys. Rev. Lett. 132 (2024) 215101]. The method uses measurements of electron temperature and density fluctuations δTₑ and δnₑ at the same spatial position to infer the local ratio of “acoustic polarization,” |δϕ ∥ |/(|δϕ ∥ |+|δψ|), where δϕ ∥ is the effective parallel potential and δψ is related to the parallel magnetic vector potential A ∥ . This paper summarizes key formulas, with emphasis on their range of validity, and elaborates on the workflow required to infer the acoustic polarization from experimental data. The drift-acoustic polarization of ellipticity-induced, toroidicity-induced, and reversed shear Alfvén eigenmodes is nearly zero, as expected for modes with predominately shear-Alfvénic polarization. The polarization of beta-induced Alfvén eigenmodes contains an acoustic component that increases with poloidal wave number. “Low frequency modes,” (instabilities that appear transiently when the minimum of the safety factor qₘᵢₙ passes through rational values) have large and highly variable acoustic polarization. In both experiment and simulation, fishbones have non-zero acoustic polarization that increases as the mode chirps down in frequency.

Alfven eigenmode↗
DOE OSTI2026

An overview of 3D field optimization for control of transport and edge instabilities on KSTAR

An international team from several laboratories and universities has made key advances over the last few years in the control of plasma transport and edge instabilities with applied 3D fields in the KSTAR tokamak to optimize long pulse operation scenarios. This overview begins with the optimization of both core and edge resonant magnetic perturbations (RMP) to improve fast ion confinement to avoid excessive limiter heat loads due to fast ion losses and successful modeling of the experimental results. Integrated and advanced plasma control techniques with machine learning (ML) and adaptive control were then used to optimize the 3D field spectrum in real-time to control edge localized modes (ELMs) while avoiding core locked modes that could disrupt the plasma. Accelerating the offline model of 3D fields with a surrogate ML model can optimize ELM suppression in the edge while limiting the impact of the applied RMP fields deeper in the plasma core in real-time. In addition, the impact of the 3D fields on the divertor heat load has been modeled and compared with experimental measurements. An analysis of a multi-machine database including KSTAR has been performed to better understand the metrics for the observed RMP thresholds for ELM suppression and the resulting plasma performance. Predictive modeling of the operational space for ELM suppression and density pumpout due to RMP has shown the importance of magnetic islands in the plasma edge and their impact on plasma turbulence. This research has culminated in the development of successful long pulse operational scenarios on KSTAR while attempting to overcome challenges of the new tungsten divertor.

3D fields↗
DOE OSTI2025

Atomic Coherence of 2 Minutes and Instability of 1.5 × 10 −18 at 1 s in a Wannier-Stark Lattice Clock

We explore the limits of atomic coherence and measurement precision in a 87 Sr optical lattice clock. We perform a detailed characterization of key effects, including lattice Raman scattering and atomic collisions in a shallow lattice configuration, determining a 174(28) s 3 𝑃 0 clock state lifetime. Investigation of atomic coherence across a range of lattice depths and atomic densities reveals decoherence mechanisms related to photon scattering and atomic interaction. At a reduced density, we observe a coherence time of 118(9) s, approaching the fundamental limit set by spontaneous emission. Guided by this coherence understanding, we demonstrate a clock instability for an atomic ensemble of 1.5 × 10 −18 at 1 s in fractional frequency units. Our results are important for further advancing the state of the art of an optical lattice clock for fundamental physics applications.

atomic, optical, & lattice clocks↗
DOE OSTI2025

Low viscosity of solid MgO at high pressures and strain rates measured using the laser-driven Richtmyer-Meshkov instability

Solids are often assumed to behave as viscous fluids under high-strain rates. This behavior has been studied experimentally in metals but largely unexplored in brittle ceramic materials. Here, in this study, we present a technique for measuring the viscosity of MgO using time-resolved velocimetry to track the growth rate of a perturbation caused by the Richtmyer-Meshkov instability at the OMEGA EP laser facility. To interpret the results, we use an in-house Eulerian hydrocode to simulate our experiments and model the plastic deformation of solid MgO as a viscous fluid. Results indicate that MgO has a surprisingly low upper bound to its effective viscosity of ∼10 2 ⁢ Pa s at 175 ± 15⁢ GPa, ∼3500⁢K, and 10 6 −10 7 s −1 strain rate.

Perez, Tyler M. [Johns Hopkins Univ., Baltimore, M↗
DOE OSTI2024

Complete quasilinear model for the acceleration-driven lower hybrid drift instability and a computational assessment of its validity

A complete quasilinear model is derived for the electrostatic acceleration-driven lower hybrid drift instability in a uniform two-species low-beta plasma in which current is perpendicular to the background magnetic field. The model consists of coupled nonlinear velocity space diffusion equations for the volume-averaged ion and electron distribution functions. Each species' diffusion coefficient depends on a time-evolving spectral density of the electric-field energy per unit volume and a time-evolving dispersion relation. The dispersion relation is expressed analytically in integral form without the use of asymptotic limits and applies to arbitrary distribution functions, so long as they can be expressed as a function of one velocity coordinate, e.g., f⁡(vy) or f⁡(v⊥). The quasilinear model conserves energy and is complete in that it fully describes the evolution of the distribution functions, including resonant and nonresonant particle-wave interactions, while accounting for distribution-function-dependent mixed-complex frequencies. Further, the quasilinear diffusion model is solved numerically and self-consistently using a Crank-Nicolson temporal discretization and a second-order finite-volume velocity-space discretization. Numerical solutions are compared to nonlinear fourth-order accurate continuum kinetic Vlasov-Poisson simulations. Evolution of electric-field energy, growth rates, distribution functions, and diffusion coefficients are shown to be in agreement with Vlasov simulations. The quasilinear model is shown to predict anomalous transport terms, like resistivity and heating, to within a factor of order unity. Discrepancies between the quasilinear model and Vlasov simulations are assessed and attributed primarily to lack of damping in the quasilinear description and to the use of unperturbed-orbit susceptibilities in the linear theory dispersion relation. The results illuminate the predictive accuracy of the quasilinear model, place approximate bounds on its validity, and provide much needed vetting of quasilinear theory's ability to predict the nonlinear state of a microturbulent plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Absence of E 2 g Nematic Instability and Dominant A 1 g Response in the Kagome Metal CsV 3 Sb 5

Ever since the discovery of the charge density wave (CDW) transition in the kagome metal CsV 3 Sb 5 , the nature of its symmetry breaking has been under intense debate. While evidence suggests that the rotational symmetry is already broken at the CDW transition temperature ( T CDW ), an additional electronic nematic instability well below T CDW has been reported based on the diverging elastoresistivity coefficient in the anisotropic channel ( m E 2 g ). Verifying the existence of a nematic transition below T CDW is not only critical for establishing the correct description of the CDW order parameter, but also important for understanding low-temperature superconductivity. Here, we report elastoresistivity measurements of CsV 3 Sb 5 using three different techniques probing both isotropic and anisotropic symmetry channels. Contrary to previous reports, we find the anisotropic elastoresistivity coefficient m E 2 g is temperature independent, except for a step jump at T CDW . The absence of nematic fluctuations is further substantiated by measurements of the elastocaloric effect, which show no enhancement associated with nematic susceptibility. On the other hand, the symmetric elastoresistivity coefficient m A 1 g increases below T CDW , reaching a peak value of 90 at T * = 20 K . Our results strongly indicate that the phase transition at T * is not nematic in nature and the previously reported diverging elastoresistivity is due to the contamination from the A 1 g channel. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Modeling frequency instability in high-quality resonant experiments

Modern resonant sensing tools can achieve increasingly high quality factors, which correspond to extremely narrow linewidths. In such systems, time variation of the resonator’s natural frequency can potentially impact its ability to accumulate power and its resulting sensitivity. One such example is the Dark SRF experiment, which utilizes superconducting radio frequency (SRF) cavities with quality factors of Q ∼ 10 10 . Microscopic deformations of the cavity lead to stochastic jittering of its resonant frequency with amplitude 20 times its linewidth. Naively, one may expect this to lead to a large suppression in accumulated power. In this work, we study in detail the effects of frequency instability on high-quality resonant systems, utilizing the Dark SRF experiment as a case study. We show that the timescale of jittering is crucial to determining its effect on power accumulation. Namely, when the resonant frequency varies sufficiently quickly, the system accumulates power as if there were no jittering at all. This implies that the sensitivity of a jittering resonator is comparable to that of a stable resonator. In the case of Dark SRF, we find that jittering only induces a ∼ 10 % loss in power. Our results allow the dark-photon exclusion bound from Dark SRF’s pathfinder run to be refined, leading to a constraint that is an order of magnitude stronger than previously reported (corresponding to a signal-to-noise ratio that is four orders of magnitude larger). This result represents the world-leading constraint on dark photons over a wide range of masses below 6 μ eV and translates to the best laboratory-based limits on the photon mass m γ < 2.9 × 10 - 48 g .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Unraveling Grain Boundary Instability in Dense Proton-Conducting Oxides

The long-term stability of protonic ceramic electrolysis cell (PCEC) materials under high-steam operating conditions remains a critical barrier to device commercialization. Here, we investigate the fundamental degradation mechanisms of dense BaCe 0.7 Zr 0.1 Y 0.1 Yb 0.1 O 3-δ (BCZYYb) electrolytes operated at 550 °C, 50% H 2 O in air. Over 1,000 h, the total electrolyte conductivity decreases by 11.1%, driven primarily by a >130% increase in grain-boundary resistivity. Post-mortem analyses reveal that damage is localized to near-surface grain boundaries extending ∼50 μm into the dense electrolyte pellet. This surface localization indicates that degradation is likely to be severe in thin, device-level electrolytes. Degradation is primarily attributed to chemo-mechanical grain-boundary weakening arising from hydration-induced chemical expansion, culminating in the formation of intergranular cracks oriented parallel to the pellet surface. These internal cracks subsequently react with steam and/or CO 2 , leading to the formation of nanoscale insulating phases, including Ba(OH) 2 , nanocrystalline BaCO 3 , and amorphous Ce/Zr/Y/Yb-containing oxides or hydroxycarbonates. After an initial degradation period of approximately 200 h, the overall conductivity stabilizes. Incorporating NiO sintering aids reduces grain-boundary density by an order of magnitude under identical sintering conditions. Although addition of NiO increases the initial resistivity by >160% at 550 °C, it substantially suppresses grain-boundary instability and mitigates chemical degradation. These findings underscore the urgent need for chemical and/or physical stabilization of BCZYYb electrolytes and offer design guidelines to enable durable, high-performance PCECs.

08 HYDROGEN↗
DOE OSTI2024

Reconstructing Richtmyer–Meshkov instabilities from noisy radiographs using low dimensional features and attention-based neural networks

We develop an ML-based approach for density reconstruction based on transformer neural networks. This approach is demonstrated in the setting of ICF-like double shell hydrodynamic simulations wherein the parameters related to material properties and initial conditions are varied. The new method can robustly recover the complex topologies given by the Richtmyer-Meshkoff instability (RMI) from a sequence of hydrodynamic features derived from radiographic images corrupted with blur, scatter, and noise. A noise model is developed to characterize errors in extracting features from synthetic radiographs of the simulated density field. The key component of the network is a transformer encoder that acts on a sequence of features extracted from noisy radiographs. This encoder includes numerous self-attention layers that act to learn temporal dependencies in the input sequences and increase the expressiveness of the model. This approach is shown to exhibit an excellent ability to accurately recover the RMI growth rates, despite the gas-metal interface being greatly obscured by radiographic noise. Our approach can be applied in a broad array of fields involving shock physics and material science.

47 OTHER INSTRUMENTATION↗