Search NASA⌕ Search

SEARCH · Search NASA

Results for “Instability waves”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Hall interchange instability as a seed for helical magneto-Rayleigh–Taylor instabilities in magnetized liner inertial fusion Z-Pinches scaled from Z-Machine parameters to a next generation pulsed power facility

Magnetized liner inertial fusion (MagLIF) is a magneto-inertial-fusion concept that is studied on the 20-MA, 100-ns rise time Z Pulsed Power Facility at Sandia National Laboratories. Given the relative success of the platform, there is a wide interest in studying the scaled performance of this concept at a next-generation pulsed-power facility that may produce peak currents upward of 60 MA. An important aspect that requires more research is the instability dynamics of the imploding MagLIF liner, specifically how instabilities are initially seeded. It has been shown in magnetized 1-MA thin-foil liner Z-pinch implosion simulations that a Hall interchange instability (HII) effect can provide an independent seeding mechanism for helical magneto-Rayleigh–Taylor instabilities. Here in this paper, we explore this instability at higher peak currents for MagLIF using 2D discontinuous Galerkin PERSEUS simulations, an extended magneto-hydrodynamics code, which includes Hall physics. Our simulations of scaled MagLIF loads show that the growth rate of the HII is invariant to the peak current, suggesting that studies at 20-MA are directly relevant to 60-MA class machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Parameterized anomalous transport model for current-carrying collisionless plasmas in pulsed power inertial confinement fusion

Current delivery in pulsed power inertial confinement fusion is influenced by collisionless current-carrying microturbulent plasmas, which are sourced from electrode surfaces. In this setting, the lower hybrid drift instability—triggered by plasma acceleration—is a leading candidate driver of difficult-to-predict momentum and energy transport. To characterize the nonlinear state of the microturbulent plasma, a parameterized anomalous transport model is developed for the instability, with analytic formulas for anomalous collision frequency, resistivity, and species heating rates. The formulas are expressed in terms of linear-theory variables and four dimensionless parameters that characterize the macroscopic plasma state. The model is built on linear theory analysis, power law analysis, and quasilinear theory analysis, and is validated using a series of nonlinear continuum kinetic Vlasov–Poisson simulations. The theoretical and computational investigation demonstrates that the anomalous collision frequency associated with the instability can be reliably approximated, within about a factor of five or better, by the unscaled linear theory growth rate of the fastest-growing wavenumber mode. This finding enables efficient calculation of anomalous resistivity and species heating rates over a wide range of plasma conditions, resulting in improved predictive capabilities.

Complex functions↗
DOE OSTI2020

Valley magnetism, nematicity, and density wave orders in twisted bilayer graphene

We analyze density wave and Pomeranchuk orders in twisted bilayer graphene. This complements our earlier analysis of the pairing instabilities. We assume that near half filling of either conduction or valence band, the Fermi level is close to Van Hove points, where the density of states diverges, and study potential instabilities in the particle-hole channel within a patch model with two valley degrees of freedom. The hexagonal symmetry of twisted bilayer graphene allows for either six or twelve Van Hove points. We consider both cases and find the same two leading candidates for particle-hole order. One is an SU(2)-breaking spin state with ferromagnetism within a valley. A subleading intervalley hopping induces antiferromagnetism between the valleys. The same state has also been obtained in strong-coupling approaches, indicating that this order is robust. Here, the other is a mixed state with 120° complex spin order and orthogonal complex charge order. In addition, we find a weaker but still attractive interaction in nematic channels, and discuss the type of a nematic order.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Multiscale Richtmyer-Meshkov instability experiments to isolate the strain rate dependence of strength

Theoretical analysis of Richtmyer-Meshkov instability (RMI) experiments for solid strength shows that the strain rate for a given shock should be inversely proportional to the length scale of the sine wave perturbations when η 0 k , the nondimensional amplitude to wavelength ratio, is held fixed. To isolate the effect of strain rate on strength, free-surface RMI specimens of annealed copper were prepared with three perturbation regions with the same η 0 k but different length scales, characterized by the wavelength λ varying by a factor of 4.9 from 65 to 130 to 320 µ m . Three such targets with different fixed η 0 k ′ s were impacted to a shock pressure of 25 GPa, and the instability evolution was measured with photon Doppler velocimetry. Strengths estimated by comparing hydrocode simulation to the data increased from 700 to 1200 MPa as λ decreased. The different η 0 k targets exercised increasing amounts of plastic strain yet showed no evidence of strain hardening. Physical regime sensitivity analysis determined that for 320 − 65 µ m wavelength perturbations, the effective strain rates increased from 8.7 × 10 6 to 3.3 × 10 7 s − 1 , a factor of 3.8. Thus, the predicted strain rate scaling was mostly achieved but slightly suppressed by increased strength at higher rates. The RMI strength estimates were plotted against constitutive testing data on copper from the literature to show striking evidence of the strength upturn at higher strain rates. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2020

Beam-ion Studies in NSTX and NSTX-U (Final Technical Report)

The confinement of high energy "fast ions" is crucial for the success of magnetic fusion as a practical energy source. The research performed by UC Irvine on NSTX and NSTX-U provided new information about this important topic in the configuration known as a "spherical tokamak." Injected neutral beams provided the fast ions (also known as "beam ions"). There were five overarching goals of the research. One objective was to measure the confinement of beam ions in order to ascertain if they were accomplishing their desired purpose of transferring their energy to the bulk plasma. A second goal was to better understand instabilities that are driven unstable by the fast ions. A related goal was to use the similarities and differences between NSTX and the DIII-D conventional tokamak to better understand the fast-ion driven instabilities. A fourth goal was to develop new instruments to measure fast ions and techniques that facilitate analysis of the data. The fifth objective was to measure and understand acceleration of beam ions by RF waves. Much was accomplished in all of these five areas. In the first, it was shown that the large magnetic moment of the beam ions did not harm confinement but instabilities known as the "sawtooth" and "long-lived mode" do. In the second area, much attention was devoted to the fast-ion instabilities known as "Alfvén eigenmodes." In particular, Alfvén eigenmode "avalanches" can cause ~ 30% of the fast ions to be lost in a single explosive burst. The comparative studies with DIII-D showed that an instability discovered on NSTX is also important in conventional tokamaks. In the fourth category, several instruments were developed, a number of effects that complicate interpretation of the data were understood, and progress toward a new method to infer the fast-ion distribution function from the data was made. In the fifth category, the measured profile of accelerated fast ions was initially much broader than theoretical predictions but subsequent improvements in the theoretical modeling achieved good agreement with the data.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

High spatial resolution and contrast radiography of hydrodynamic instabilities at the National Ignition Facility

In this work, we are developing techniques for studying the Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities in a planar geometry at high-energy-densities at the National Ignition Facility (NIF). In particular, through the improvement of experimental imaging quality, we are progressing toward the study of the turbulent regime of the mixing regions in capsule implosion experiments for inertial confinement fusion, which requires few micrometers resolution. Using 60 NIF beams, a solid shock tube is driven launching a shock wave that crosses the interface between a dense and a light material pre-machined in the target to obtain sinusoidal ripples, which results in RM and RT instabilities that are imaged using the NIF Crystal Backlighter Imager. High-quality images were obtained with a mean resolution of 7 μm and improved contrast. While the obtained resolution does not allow the observation of the smallest scale of the “turbulent” energy spectrum, the generated image encompasses 63% of the total flow energy, a 50% improvement over previous studies, which is observed for the first time a roll-up feature in a high energy density-type RT experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Hydrostatic and Non‐Hydrostatic Baroclinic Instability in the Dynamical Core of the DOE Global Climate Model

Abstract The dynamical core of the Department of Energy global climate model is used to understand the role of non‐hydrostatic dynamics in the simulation of dry and moist baroclinic waves. To reduce computational cost, the Diabatic Acceleration and REscaling approach is adopted. Scale analysis and numerical simulations suggest that the model solution is not distorted by the change of spherical metric terms due to the change of Earth radius, neither are the inter‐scale interactions strongly altered as indicated by the analysis of spectral flux. Compared with hydrostatic simulations, the onset of baroclinic instability is delayed under non‐hydrostatic dynamics as the associated weaker vertical motions tend to increase the critical wavelength and narrow the range of waves that can be baroclinically unstable. During the development of baroclinic waves, non‐hydrostatic dynamics tends to induce vertical motions in the upper troposphere, accelerating the eastward propagation of upper‐level ridges through their impact on local vorticity tendency. These processes reduce the westward tilt of the vertical ridge axes and suppress the conversion of mean flow available potential energy to eddy kinetic energy, leading to weaker baroclinic eddies than in hydrostatic simulations. The contrasts between hydrostatic and non‐hydrostatic settings hold in supplemental experiments that vary the background flow Rossby number, the amount of water vapor content and vertical resolution. We also find that mesoscale and smaller‐scale activities can be considerably under‐represented when the vertical resolution is limited to that typically used in global climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A note on the instability and pattern formation of shrinkage cracks in viscoplastic soils

In this work we present a theoretical study on the conditions for the onset of cracks, as well as the corresponding pattern formation, in saturated viscoplastic soils under isotropic loading (extension). The type of stress applied is left unspecified, to cover a variety of loadings including shrinkage due to desiccation, isotropic thermal expansion, mechanical loading and so forth. By treating the saturated soil as rigid viscoplastic, we obtain a 2D extension of the Cnoidal Waves equations (Veveakis and Regenauer-Lieb, 2015). By numerically solving the corresponding boundary value problem, we retrieve conditions for the onset of cracking instability in 2D loading, and identify the characteristic spacing between cracks to be a length scale combining all the hydro-mechanical parameters of the problem. Finally, we show that in a rectangular slab of clay under isotropic extension, patterns of triangular, rectangular and hexagonal cracks can tessellate the domain, with the hexagonal pattern being the energetically favored, as it minimizes the free energy of the system.

58 GEOSCIENCES↗
DOE OSTI2020

Gravitational-wave signal of a core-collapse supernova explosion of a 15M ⊙ star

In this paper we report on the gravitational-wave signal computed in the context of a three-dimensional simulation of a core-collapse supernova explosion of a 15 M ⊙ star. The simulation was performed with our neutrino hydrodynamics code chimera . We detail the gravitational wave strains as a function of time, for both polarizations, and discuss their physical origins. We also present the corresponding spectral signatures. Gravitational wave emission in our model has two key features: low-frequency emission (less than 200 Hz) emanates from the gain layer as a result of neutrino-driven convection and the standing accretion shock instability (SASI), and high-frequency emission (greater than 600 Hz) emanates from the proto–neutron star due to Ledoux convection within it. The high-frequency emission dominates the gravitational wave emission in our model and emanates largely from the convective layer itself, not from the convectively stable layer above it, due to convective overshoot. Moreover, the low-frequency emission emanates from the gain layer itself, not from the proto–neutron star, due to accretion onto it. We provide evidence of the SASI in our model and demonstrate that the peak of our low-frequency gravitational wave emission spectrum corresponds to it. Given its origin in the gain layer, we classify the SASI emission in our model as p-mode emission and assign a purely acoustic origin, not a vortical–acoustic origin, to it. We compare the results of our three-dimensional model analysis with those obtained from the model’s two-dimensional counterpart and find a significant reduction in the strain amplitudes in the former case, as well as significant reductions in all related quantities. Our dominant proto–neutron star gravitational wave emission is not well characterized by emission from surface g modes, complicating the relationship between peak frequencies observed and the mass and radius of the proto–neutron star expressed by analytic estimates under the assumption of surface g-mode emission. We present our frequency normalized characteristic strain along with the sensitivity curves of current- and next-generation gravitational wave detectors. This simple analysis indicates that the spectrum of gravitational wave emission between approximately 20 Hz and approximately 1 kHz, stemming from neutrino-driven convection, the SASI, accretion onto the proto–neutron star, and proto–neutron star convection, will be accessible for a Galactic event.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Topological and nematic superconductivity mediated by ferro-SU(4) fluctuations in twisted bilayer graphene

We propose an SU(4) spin-valley-fermion model to investigate the superconducting instabilities of twisted bilayer graphene (TBG). In this approach, bosonic fluctuations associated with an emergent SU(4) symmetry, corresponding to combined rotations in valley and spin spaces, couple to the low-energy fermions that comprise the flat bands. These fluctuations are peaked at zero wave vector, reflecting the “ferromagnetic-like” SU(4) ground state recently found in strong-coupling solutions of microscopic models for TBG. Focusing on electronic states related to symmetry-imposed points of the Fermi surface, dubbed here “valley hot-spots” and “van Hove hot-spots,” we find that the coupling to the itinerant electrons partially lifts the huge degeneracy of the ferro-SU(4) ground state manifold, favoring intervalley order, spin-valley coupled order, ferromagnetic order, spin-current order, and valley-polarized order, depending on details of the band structure. Here, these fluctuations, in turn, promote attractive pairing interactions in a variety of closely competing channels, including a nodeless f-wave state, a nodal i-wave state, and topological d+id and p+ip states with unusual Chern numbers 2 and 4, respectively. Nematic superconductivity, although not realized as a primary instability of the system, appears as a consequence of the near-degeneracy of superconducting order parameters that transform as one-dimensional and two-dimensional irreducible representations of the point group D 6 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Interpreting and Stabilizing Machine-Learning Parametrizations of Convection

Neural networks are a promising technique for parameterizing subgrid-scale physics (e.g., moist atmospheric convection) in coarse-resolution climate models, but their lack of interpretability and reliability prevents widespread adoption. For instance, it is not fully understood why neural network parameterizations often cause dramatic instability when coupled to atmospheric fluid dynamics. This paper introduces tools for interpreting their behavior that are customized to the parameterization task. First, we assess the nonlinear sensitivity of a neural network to lower-tropospheric stability and the midtropospheric moisture, two widely studied controls of moist convection. Second, we couple the linearized response functions of these neural networks to simplified gravity wave dynamics, and analytically diagnose the corresponding phase speeds, growth rates, wavelengths, and spatial structures. To demonstrate their versatility, these techniques are tested on two sets of neural networks, one trained with a super-parameterized version of the Community Atmosphere Model (SPCAM) and the second with a near-global cloud-resolving model (GCRM). Additionally, even though the SPCAM simulation has a warmer climate than the cloud-resolving model, both neural networks predict stronger heating/drying in moist and unstable environments, which is consistent with observations. Moreover, the spectral analysis can predict that instability occurs when GCMs are coupled to networks that support gravity waves that are unstable and have phase speeds larger than 5 m s -1 . In contrast, standing unstable modes do not cause catastrophic instability. Using these tools, differences between the SPCAM-trained versus GCRM-trained neural networks are analyzed, and strategies to incrementally improve both of their coupled online performance unveiled.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Multi-scale Interaction Mechanism for Edge-Localized-Mode Suppression in the Tokamak Edge

A central challenge in fusion energy is reconciling the high-confinement mode required for reactor performance with the intense intermittent relaxation events it produces, known as edge-localized modes. These instabilities arise in the steep pressure pedestal at the plasma edge when magnetohydrodynamic thresholds are crossed, inflicting damaging heat loads on reactor components. Here, we show that multiscale interactions between microscopic turbulence and macroscopic magnetohydrodynamic modes provide encouraging prospects for self-organized edge-localized modes regulation. Using direct quantitative measurements of multiscale modes, eddy dynamics, and turbulent flux, we show that small-scale electron drift wave turbulence actively scatters the large-scale peeling-ballooning modes. This scattering decorrelates the pressure and velocity fields of the instability, so arresting its growth. Our modeling and theoretical analysis confirm this suppression mechanism is effective even when conventional linear stability thresholds are exceeded. This work establishes a nonlinear principle for edge-localized modes stability, revealing how ambient micro-turbulence can be leveraged to maintain a macro-stable, high-performance pedestal for future fusion reactors.

experimental nuclear physics↗
DOE OSTI2020

Robust d x 2 - y 2 -wave superconductivity of infinite-layer nickelates

Motivated by the recent observation of superconductivity in strontium-doped NdNiO 2 , we study the superconducting instabilities in this system from various vantage points. Starting with first-principles calculations, we construct two distinct tight-binding models, a simpler single-orbital as well as a three-orbital model, both of which capture the key low-energy degrees of freedom to varying degrees of accuracy. We study superconductivity in both models using the random phase approximation. We then analyze the problem at stronger coupling, and study the dominant pairing instability in the associated t-J model limit. Finally, in all instances, the dominant pairing tendency is in the d x 2 -y 2 channel, analogous to the cuprate superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Unravelling ultralow thermal conductivity in perovskite Cs2AgBiBr6: dominant wave-like phonon tunnelling and strong anharmonicity

Abstract Understanding the lattice dynamics and heat transport physics in the lead-free halide double perovskites remains an outstanding challenge due to their lattice dynamical instability and strong anharmonicity. In this work, we investigate the microscopic mechanisms of anharmonic lattice dynamics and thermal transport in lead-free halide double perovskite Cs 2 AgBiBr 6 from first principles. We combine self-consistent phonon calculations with bubble diagram correction and a unified theory of lattice thermal transport that considers both the particle-like phonon propagation and wave-like tunnelling of phonons. An ultra-low thermal conductivity at room temperature (~0.21 Wm −1 K −1 ) is predicted with weak temperature dependence( ~ T −0.34 ), in sharp contrast to the conventional ~T −1 dependence. Particularly, the vibrational properties of Cs 2 AgBiBr 6 are featured by strong anharmonicity and wave-like tunnelling of phonons. Anharmonic phonon renormalization from both the cubic and quartic anharmonicities are found essential in precisely predicting the phase transition temperature in Cs 2 AgBiBr 6 while the negative phonon energy shifts induced by cubic anharmonicity has a significant influence on particle-like phonon propagation. Further, the contribution of the wave-like tunnelling to the total thermal conductivity surpasses that of the particle-like propagation above around 310 K, indicating the breakdown of the phonon gas picture conventionally used in the Peierls-Boltzmann Transport Equation. Importantly, further including four-phonon scatterings is required in achieving the dominance of wave-like tunnelling, as compared to the dominant particle-like propagation channel when considering only three-phonon scatterings. Our work highlights the importance of lattice anharmonicity and wave-like tunnelling of phonons in the thermal transport in lead-free halide double perovskites.

Chemistry↗
DOE OSTI2025

Generalizing the compressible pairwise interaction extended point-particle model

Ejecta physics plays an important role in material interfaces that are impacted by a strong shock wave. When a shock impacts a rough surface of solid material and melts it, the Richtmyer–Meshkov instability grows perturbations on the surface, which can eject particles. After release, the ejecta travel through the post-shock compressible flow. To accurately simulate a large number of ejecta particles, an Euler–Lagrange approach is preferred, which requires modeling the subgrid-scale physics involved with fluid–particle interactions. We generalize the previous work from Hsiao et al. (2023) to consider systems of moving particles subject to any loading shock. The following improvements were made: (1) Particles are allowed to move relative to each other (2) Non-planar shocks are accounted for along with allowing for variable shock speeds. As a result, the generalized algorithm was tested with particle-resolved simulations for canonical test cases. The results of these tests are discussed and analyzed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A point-particle-based hydride shell-shedding model for ejecta particle transport in reactive environments

A shock wave passing over a rough or perturbed metal surface will induce a limiting case of Richtmyer–Meshkov instability and will cause small particles to eject from the surface and transport into the surrounding medium. These particles are known as ejecta and can be either solid or liquid in nature. Recent experiments have shown that liquid cerium ejecta clouds exhibit unexpected non-monotonic acceleration behaviors as well as temperature plateaus after a brief temperature rise if they are transporting in a chemically reactive, hydrogen-based medium while they act as expected in an inert medium. This work details a point-particle model developed for reactive cerium ejecta transport, which attempts to account for these new physics through the behavior of a developing solid hydride shell, which is believed to form as a product of the reaction. The overall model incorporates the effects of the reaction on the particle properties as well as the effects of potential shedding of the shell into sub-micrometer scale flakes and potential phase change of the hydride if the ejecta particles reach the melt point of the hydride layer. The model is tested by performing simulations of the original motivating experiments and comparing quantities, such as ejected mass, velocimetry, and temperature profiles, against the experimental data. While the model is able to capture many general features of the observed anomalies, some inaccuracies still exist. These point to both missing physics in the model (such as a deuterium adsorption mechanism on the hydride layer) as well as a lack of knowledge of certain material properties (such as the strength of cerium hydride to determine dynamic fracture thicknesses) needed to fully reduce the uncertainties in the model by up to an order of magnitude and perform a true attempt at model validation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Townsend to glow discharge transition for a nanosecond pulse plasma in helium: space charge formation and resulting electric field dynamics

Stark polarization spectroscopy is used to investigate the temporal evolution of the electric field distribution in the cathode region of a nanosecond pulsed discharge in helium at 120 Torr. The measurements are performed on the He I transition at 492.19 nm, during the early stages of the discharge formation. The experimental results are compared with the predictions of a 1D fluid model. Time-resolved ICCD images show that the discharge develops as a diffuse, cathode-directed ionization wave with a Townsend-like feature before transitioning into a glow-like structure. Near anode instabilities characterized by filament formation were observed near the high voltage electrode. Within 30 ns, a reduction of the sheath thickness to about 250 μm is observed, coinciding with a gradual increase of the discharge current and proportional increase in electric field at the cathode. Additionally, the cathode electric field corresponding to this sheath with a thickness of 250 μm is about 40 kV cm -1 . A subsequent steep increase of the discharge current leads to a further reduction of the sheath width. The electric field evolution as obtained by the fluid model is in excellent agreement with the measurements and shows that an enhanced ionization near the cathode is causing the space charge formation responsible for the increase in electric field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Role of electromagnetic gauge-field fluctuations in the selection between chiral and nematic superconductivity

Motivated by the observation of nematic superconductivity in several systems, we revisit the problem of the leading pairing instability of two-component unconventional superconductors on the triangular lattice—such as (p x , p y )-wave and (d x 2 –y 2 , d xy )-wave superconductors. Such a system has two possible superconducting states: the chiral state (e.g., p + ip or d + id), which breaks time-reversal symmetry, and the nematic state (e.g., p + p or d + d), which breaks the threefold rotational symmetry of the lattice. Weak-coupling calculations generally favor the chiral over the nematic superconducting state, raising the question of what mechanism can stabilize the latter. Here, we show that the electromagnetic field fluctuations can play a crucial role in selecting between these two states. Specifically, we derive and analyze the effective free energy for the two-component superconducting order parameter after integrating out the gauge-field fluctuations, which is formally justified if the spatial order parameter fluctuations can be neglected. A nonanalytic cubic term arises, as in the case of a conventional s -wave superconductor. However, unlike the latter, the cubic term depends on the relative phase and on the relative amplitudes between the two order parameter components, in such a way that it generally favors the nematic state. This result is a direct consequence of the fact that the stiffness of the superconducting order parameter is not isotropic. Competition with the quartic term, which favors the chiral state, leads to a renormalized phase diagram in which the nematic state displaces the chiral state over a wide region in the parameter space. Here, we analyze the stability of the fluctuation-induced nematic phase, generalize our results to tetragonal lattices, and discuss their applicability to candidate nematic superconductors, including twisted bilayer graphene.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗