Search NASA⌕ Search

SEARCH · Search NASA

Results for “Shock 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Detectability of standing accretion shock instabilities activity in supernova neutrino signals

We introduce a novel methodology for establishing the presence of standing accretion shock instabilities (SASI) in the dynamics of a core collapse supernova from the observed neutrino event rate at water- or ice- based neutrino detectors. The methodology uses a likelihood ratio in the frequency domain as a test- statistics; it is also employed to assess the potential to estimate the frequency and the amplitude of the SASI modulations of the neutrino signal. The parameter estimation errors are consistent with the minimum possible errors as evaluated from the inverse of the Fisher information matrix, and close to the theoretical minimum for the SASI amplitude. Using results from a core-collapse simulation of a 15 solar-mass star by Kuroda et al. [Astrophys. J. 851, 62 (2017).] as a test bed for the method, we find that SASI can be identified with high confidence for a distance to the supernova of up to ~6 kpc for IceCube and up to ~3 kpc for a 0.4 Mt mass water Cherenkov detector. This methodology will aid the investigation of a future galactic supernova.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Characterizing a supernova’s standing accretion shock instability with neutrinos and gravitational waves

Here, we perform a novel multi-messenger analysis for the identification and parameter estimation of the Standing Accretion Shock Instability (SASI) in a core collapse supernova with neutrino and gravitational wave (GW) signals. In the neutrino channel, this method performs a likelihood ratio test for the presence of SASI in the frequency domain. For gravitational wave signals we process an event with a modified constrained likelihood method. Using simulated supernova signals, the properties of the Hyper-Kamiokande neutrino detector, and O3 LIGO Interferometric data, we produce the two-dimensional probability density distribution (PDF) of the SASI activity indicator and calculate the probability of detection P D as well as the false identification probability P FI . We discuss the probability to establish the presence of the SASI as a function of the source distance in each observational channel, as well as jointly. Compared to a single-messenger approach, the joint analysis results in P D (at P FI = 0.1) of SASI activities that is larger by up to ≈ 40% for a distance to the supernova of 5 kpc. We also discuss how accurately the frequency and duration of the SASI activity can be estimated in each channel separately. Our methodology is suitable for implementation in a realistic data analysis and a multi-messenger setting.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Modeling Hydrodynamic Instabilities, Shocks, and Radiation Waves in High Energy Density Experiments [Dissertation]

This thesis presents the computational design, modeling, and analysis of three experiments in high-energy density physics (HEDP), all of them concerning fundamental radiation flows. The first experiment is a laboratory astrophysics experiment to investigate the role of the Kelvin-Helmholtz instability (KHI) in the process of galactic filaments supplying gas to galactic halos. The achieved goal was to provide a first study in which the role of the instability is maximal and predict behavior in future iterations of an experiment accessing a more radiative regime where the role of the instability is stifled. This experiment would help answer how certain galaxies are able to grow so rapidly and produce many stars as the KHI limits this process. The second experiment, COAX, is a radiation flow experiment with a novel spectroscopy diagnostic configuration, designed to spatially measure the temperature of a radiation wave as it travels down a doped foam. A key result of this work was the development of a synthetic spectroscopy application and application of modern spectroscopy comparison techniques to provide our first temperature reconstructions from the experimental data. This experimental platform serves as the launching ground for a number of new experiments that vary the basic premise and thus is foundational to our ongoing research. The final experiment is a full integration of modeling, design, and theoretical development for the Radishock experiment. This experimental platform studies the head-on collision of a radiation wave with a counter-propagating shock, and like COAX, uses spectroscopy to diagnose and detect the interaction. My research analyzes the successful shots, indicating aspects of successful detections and suggests improvements to future iterations of the design.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Same-sided successive-shock HED instability experiments

Inertial confinement fusion (ICF) and high-energy density (HED) physics experiments experience complicated forcing for instability growth and mix due to the ubiquitous presence of multiple shocks interacting with perturbations on multiple material interfaces. One common driver of instability growth is successive shocks from the same direction. However, there is a severe lack of analytic work and modeling validation for same-sided successive shocks since they are extremely difficult to achieve with conventional (non-HED) drivers. Successive shocks access a large instability parameter space; idealized fluid theory [K. O. Mikaelian, Phys. Rev. A 31, 410 (1985)] predicts 15 different interface evolution scenarios for a sinusoidal perturbation. Growth becomes more complex for multi-mode, compressible HED systems. The Mshock campaign is the first experiment in any fluid regime to probe a wide portion of successive shock parameter space. This is enabled by our development of a hybrid direct/indirect drive platform capable of creating independently controllable successive shocks on the National Ignition Facility. These experiments have delivered the first data capable of rigorously challenging our models and their ability to accurately capture Richtmyer–Meshkov growth under successive shocks. Single-mode and two-mode experiments have successfully demonstrated the ability to access and control the various growth scenarios of the shocked interface, including re-inversion, freeze out, and continued growth. Simulations and theoretical modeling are shown to accurately capture the experimental observations in the linear growth phase, giving us confidence in our ICF/HED design codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Visualization of shocked material instabilities using a fast-framing camera and XFEL four-pulse train

Many questions regarding dynamic materials could be answered by using time-resolved ultra-fast imaging techniques to characterize the physical and chemical behavior of materials in extreme conditions and their evolution on the nanosecond scale. In this work, we perform multi-frame phase-contrast imaging (PCI) of micro-voids in low density polymers under laser-driven shock compression. At the Matter in Extreme Conditions (MEC) Instrument at the Linac Coherent Light Source (LCLS), we used a train of four x-ray free electron laser (XFEL) pulses to probe the evolution of the samples. To visualize the void and shock wave interaction, here, we deployed the Icarus V2 detector to record up to four XFEL pulses, separated by 1-3 nanoseconds. In this work, we image elastic waves interacting with the micro-voids at a pressure of several GPa. Monitoring how the material’s heterogeneities, like micro-voids, dictate its response to a compressive wave is important for benchmarking the performances of inertial confinement fusion energy materials. For the first time in a single sample, we have combined an ultrafast x-ray framing camera and four XFEL pulse train to create an ultrafast movie of micro-void evolution under laser-driven shock compression. Eventually, we hope this technique will resolve the material density as it evolves dynamically under laser shock compression.

fusion↗
DOE OSTI2023

Intricate structure of the plasma Rayleigh–Taylor instability in shock tubes

Spikes and bubbles grow on unstable interfaces that are accelerated in high-energy-density conditions. If a shock propagates ahead of the interface, the plasma can be heated to extreme conditions where conduction and radiation fluxes influence the hydrodynamics. For example, a National Ignition Facility experiment found reduced single-mode nonlinear mixed-width growth in conditions scaled from a supernova explosion. We present high-resolution two-dimensional radiation hydrodynamic simulations with the Flash code that quantitatively reproduce the experiment. Radiative fluxes are primarily responsible for ablating the spike and removing the mushroom caps. The ablated plasma increases the mixed mass and forms a low-density halo with spikes forming in both directions. This is considerably more complex than classical instability. The halo is sensitive to ablative physics, so radiographing it may aid in the verification of energy transport modeling. Although ablation changes the spike shape, it has little effect on the overall mixed width for these parameters. This is because ablation enhances the bubble velocity but it has the opposite effect on the spike. The radiation transport instead suppresses the growth via increasing the shocked foam density, thus decreasing the Atwood number. Furthermore, a terminal velocity model including the rarefaction expansion agrees with the experimental mixed-width growth.

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 OSTI2025

Stable low diffusion flux splitting schemes on unstructured meshes

Shock instabilities are shown to manifest in modern low-diffusion flux-vector splitting (FVS) schemes when used on unstructured meshes, or situations where shocks do not align with the mesh lines. These instabilities occur irrespective of the Mach number of the shock. Three types of dissipative mechanisms that suppress these instabilities are presented. These mechanisms are carefully designed in order to affect only problematic regions of the flux-splittings. The AUSM + and LDFSS schemes are stabilized using the proposed modifications. It is shown that the added dissipation improves the shock behavior of AUSM and LDFSS on unstructured meshes. It is also shown that the AUSM + -up scheme is prone to the “carbuncle” instability, a specific type of shock instability, when used on unstructured meshes. The modifications proposed in this work do not lead to carbuncle instabilities for the problems considered here. Furthermore, the modified schemes are shown to satisfy certain properties that are crucial for accurate shear layer computations, such as stationary contact preservation. Using benchmark problems, it is demonstrated that despite the diffusion added for stabilization, these schemes are not overly diffusive. Furthermore, due to these advantages, the modified FVS schemes presented here are promising candidates for high-speed compressible flow computations on unstructured meshes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Magnetic Field Amplification by a Plasma Cavitation Instability in Relativistic Shock Precursors

Abstract Plasma streaming instabilities play an important role in magnetic field amplification and particle acceleration in relativistic shocks and their environments. However, in the far shock precursor region where accelerated particles constitute a highly relativistic and dilute beam, streaming instabilities typically become inefficient and operate at very small scales when compared to the gyroradii of the beam particles. We report on a plasma cavitation instability that is driven by dilute relativistic beams and can increase both the magnetic field strength and coherence scale by orders of magnitude to reach near-equipartition values with the beam energy density. This instability grows after the development of the Weibel instability and is associated with the asymmetric response of background leptons and ions to the beam current. The resulting net inductive electric field drives a strong energy asymmetry between positively and negatively charged beam species. Large-scale particle-in-cell simulations are used to verify analytical predictions for the growth and saturation level of the instability and indicate that it is robust over a wide range of conditions, including those associated with pair-loaded plasmas. These results can have important implications for the magnetization and structure of shocks in gamma-ray bursts, and more generally for magnetic field amplification and asymmetric scattering of relativistic charged particles in plasma astrophysical environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Isolation and phase-space energization analysis of the instabilities in collisionless shocks

We analyse the generation of kinetic instabilities and their effect on the energization of ions in non-relativistic, oblique collisionless shocks using a 3D-3V (three spatial with three velocity components) simulation by dHybridR, a hybrid particle-in-cell code. At sufficiently high Mach number, quasi-perpendicular and oblique shocks can experience rippling of the shock surface caused by kinetic instabilities arising from free energy in the ion velocity distribution due to the combination of the incoming ion beam and the population of ions reflected at the shock front. To understand the role of the ripple on particle energization, we devise a new instability isolation method to identify the unstable modes underlying the ripple and interpret the results in terms of the governing kinetic instability. We generate velocity-space signatures using the field–particle correlation technique to look at energy transfer in phase space from the isolated instability driving the shock ripple, providing a viewpoint on the different dynamics of distinct populations of ions in phase space. Together, the field–particle correlation technique and our new instability isolation method provide a unique viewpoint on the different dynamics of distinct populations of ions in phase space and allow us to completely characterize the energetics of the collisionless shock under investigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

High speed PLIF study of the Richtmyer–Meshkov instability upon re-shock

The Richtmyer–Meshkov instability (RMI) of a twice-shocked gas interface is researched using high-speed planar-laser induced fluorescence (PLIF) in the Wisconsin Shock Tube Laboratory’s vertical shock tube. The initial condition (IC) is a shear layer with broadband diffuse perturbations at the interface between a helium-acetone mixture and argon. This IC is accelerated by a shock of nominal strength M = 1.8, and then accelerated again by the transmitted shock that reflects off the end wall of the tube. An estimate of the light gas mole fraction is extracted from high-speed imaging using an iterative process that accounts for the nonlinear temperature dependence of the acetone’s fluorescence quantum yield (FQY) and absorption cross-section. A vorticity deposition model for the initial growth rate after reshock is compared with the Mikaelian model for re-shock. Previously used in literature, the number of generations is shown to naturally arise from a normalisation of the scalar transport equation. A self-similar analysis is then further performed using the mole fraction data to explore the evolution of the RMI after reshock and the higher order moments of the light gas mole fraction are compared with a proposed model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Parametric Study of the SASI Comparing General Relativistic and Nonrelativistic Treatments

We present numerical results from a parameter study of the standing accretion shock instability (SASI), investigating the impact of general relativity (GR) on the dynamics. Using GR hydrodynamics with GR gravity, and nonrelativistic (NR) hydrodynamics with Newtonian gravity, in an idealized model setting, we vary the initial radius of the shock, and by varying its mass and radius in concert, the proto-neutron star compactness. We investigate four compactnesses expected in a post-bounce core-collapse supernova (CCSN). We find that GR leads to a longer SASI oscillation period, with ratios between the GR and NR cases as large as 1.29 for the highest-compactness suite. We also find that GR leads to a slower SASI growth rate, with ratios between the GR and NR cases as low as 0.47 for the highest-compactness suite. We discuss implications of our results for CCSN simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Core collapse supernova gravitational wave emission for progenitors of 9.6, 15, and 25 M ⊙

Here, we present gravitational wave emission predictions based on three core collapse supernova simulations corresponding to three different progenitor masses. The masses span a large range, between 9.6 and 25M ⊙ , are all initially nonrotating, and are of two metallicities: zero and solar. We compute both the temporal evolution of the gravitational wave strains for both the plus and the cross polarizations, as well as their spectral decomposition and characteristic strains. The temporal evolution of our zero metallicity 9.6M ⊙ progenitor model is distinct from the temporal evolution of our solar metallicity 15M ⊙ progenitor model and our zero metallicity 25M ⊙ progenitor model. In the former case, the high-frequency gravitational wave emission is largely confined to a brief time period ~75 ms after bounce, whereas in the latter two cases high-frequency emission does not commence until ~125 ms after bounce or later. The excitation mechanisms of the high-frequency emission in all three cases correspond to proto-neutron star convection and accretion onto the proto-neutron star from the convective gain layer above it, with the former playing the dominant role for most of the evolution. The low-frequency emission in all three models exhibits very similar behavior. At frequencies below ~250 Hz, gravitational waves are emitted by neutrino-driven convection and the standing accretion shock instability (SASI). This emission extends throughout the simulations when a gain region is present. In all three models, explosion is observed at ~125, ~500, and ~250 ms after bounce in the 9.6, 15, and 25M ⊙ progenitor models, respectively. At these times, the low-frequency gravitational wave emission is joined by very low-frequency emission, below ~10 Hz. These very low-frequency episodes are the result of explosion and begin at the above designated explosion times in each of our models. Our characteristic strains tell us that, in principle, all three gravitational wave signals would be detectable by current-generation detectors for a supernova at a distance of 10 kpc. However, our 9.6M ⊙ progenitor model is a significantly weaker source of gravitational waves, with strain amplitudes approximately 5–10 times less than in our other two models. The characteristic strain for this model tells us that such a supernova would be detectable only within a much more narrow frequency range around the maximum sensitivity of today’s detectors. Finally, in our 9.6M ⊙ progenitor model, we see very high-frequency gravitational radiation, extending up to ~2000 Hz. This feature results from the interaction of shock- and deleptonization-induced convection with perturbations introduced in the progenitor by nuclear burning during core collapse. While unique to the 9.6M ⊙ progenitor model analyzed here, this very high-frequency emission may, in fact, be a generic feature of the predictions for the gravitational wave emission from all core collapse supernova models when simulations are performed with three-dimensional progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Core-collapse supernova neutrino emission and detection informed by state-of-the-art three-dimensional numerical models

ABSTRACT Based on our recent three-dimensional core-collapse supernova (CCSN) simulations including both exploding and non-exploding models, we study the detailed neutrino signals in representative terrestrial neutrino observatories, namely Super-Kamiokande (Hyper-Kamiokande), DUNE, JUNO, and IceCube. We find that the physical origin of difference in the neutrino signals between 1D and 3D is mainly proto-neutron-star convection. We study the temporal and angular variations of the neutrino signals and discuss the detectability of the time variations driven by the spiral standing accretion shock instability (spiral SASI) when it emerges for non-exploding models. In addition, we determine that there can be a large angular asymmetry in the event rate (${\gtrsim} 50 {{\ \rm per\ cent}}$), but the time-integrated signal has a relatively modest asymmetry (${\lesssim} 20 {{\ \rm per\ cent}}$). Both features are associated with the lepton-number emission self-sustained asymmetry and the spiral SASI. Moreover, our analysis suggests that there is an interesting correlation between the total neutrino energy (TONE) and the cumulative number of neutrino events in each detector, a correlation that can facilitate data analyses of real observations. We demonstrate the retrieval of neutrino energy spectra for all flavours of neutrino by applying a novel spectrum reconstruction technique to the data from multiple detectors. We find that this new method is capable of estimating the TONE within the error of ∼20 per cent if the distance to the CCSN is ≲6 kpc.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Black hole to breakout: 3D GRMHD simulations of collapsar jets reveal a wide range of transients

We present a suite of the first 3D GRMHD collapsar simulations, which extend from the self-consistent jet launching by an accreting Kerr black hole (BH) to the breakout from the star. We identify three types of outflows, depending on the angular momentum, l, of the collapsing material and the magnetic field, B, on the BH horizon: (i) subrelativistic outflow (low l and high B), (ii) stationary accretion shock instability (SASI; high l and low B), (iii) relativistic jets (high l and high B). In the absence of jets, free-fall of the stellar envelope provides a good estimate for the BH accretion rate. Jets can substantially suppress the accretion rate, and their duration can be limited by the magnetization profile in the star. We find that progenitors with large (steep) inner density power-law indices (≳ 2), face extreme challenges as gamma-ray burst (GRB) progenitors due to excessive luminosity, global time evolution in the light curve throughout the burst and short breakout times, inconsistent with observations. Our results suggest that the wide variety of observed explosion appearances (supernova/supernova + GRB/low-luminosity GRBs) and the characteristics of the emitting relativistic outflows (luminosity and duration) can be naturally explained by the differences in the progenitor structure. Our simulations reveal several important jet features: (i) strong magnetic dissipation inside the star, resulting in weakly magnetized jets by breakout that may have significant photospheric emission and (ii) spontaneous emergence of tilted accretion disc-jet flows, even in the absence of any tilt in the progenitor.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2024

Feasibility of an experiment on clumping induced by the Crow instability along a shocked cylinder

The growth of three-dimensional perturbations subject to the Crow instability along a vortex dipole resulting from the passage of a shock wave through a heavy gaseous cylinder is examined numerically. A linear stability analysis is performed based on geometric parameters extracted from two-dimensional simulations to determine the range of unstable wavenumbers, which is found to extend from 0.0 to 1.3 when normalized by the core separation distance. The analysis is then verified by comparison to three-dimensional simulations, which clearly show the development of the instability and the pinch-off of the vortex dipole into isolated vortex rings, which manifest as clumps of the original cylinder material. A scaling law is developed to determine the relevant spatiotemporal scales of the instability development, which is then used to assess the feasibility of a high-energy-density experiment visualizing clump formation. Specifically, a shocked cylinder with an initial diameter of 100 μm consisting of a perturbation of approximate wavelength and amplitude of 600 and 10 μm, respectively, is expected to form clumps resulting from the Crow instability approximately 40 ns after it is shocked, with dynamics which can be readily visualized on the Omega EP laser facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗