Search NASA⌕ Search

SEARCH · Search NASA

Results for “Explosive instabilities”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

SNEWPY: A Data Pipeline from Supernova Simulations to Neutrino Signals

Current neutrino detectors will observe hundreds to thousands of neutrinos from a Galactic supernovae, and future detectors will increase this yield by an order of magnitude or more. With such a data set comes the potential for a huge increase in our understanding of the explosions of massive stars, nuclear physics under extreme conditions, and the properties of the neutrino. However, there is currently a large gap between supernova simulations and the corresponding signals in neutrino detectors, which will make any comparison between theory and observation very difficult. SNEWPY is an open-source software package which bridges this gap. The SNEWPY code can interface with supernova simulation data to generate from the model either a time series of neutrino spectral fluences at Earth, or the total time-integrated spectral fluence. Data from several hundred simulations of core-collapse, thermonuclear, and pair-instability supernovae is included in the package. This output may then be used by an event generator such as sntools or an event rate calculator such as SNOwGLoBES. Additional routines in the SNEWPY package automate the processing of the generated data through the SNOwGLoBES software and collate its output into the observable channels of each detector. In this paper we describe the contents of the package, the physics behind SNEWPY, the organization of the code, and provide examples of how to make use of its capabilities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

SNEWPY: A Data Pipeline from Supernova Simulations to Neutrino Signals

Current neutrino detectors will observe hundreds to thousands of neutrinos from Galactic supernovae, and future detectors will increase this yield by an order of magnitude or more. With such a data set comes the potential for a huge increase in our understanding of the explosions of massive stars, nuclear physics under extreme conditions, and the properties of the neutrino. However, there is currently a large gap between supernova simulations and the corresponding signals in neutrino detectors, which will make any comparison between theory and observation very difficult. SNEWPY is an open-source software package that bridges this gap. The SNEWPY code can interface with supernova simulation data to generate from the model either a time series of neutrino spectral fluences at Earth, or the total time-integrated spectral fluence. Data from several hundred simulations of core-collapse, thermonuclear, and pair-instability supernovae is included in the package. This output may then be used by an event generator such as sntools or an event rate calculator such as the SuperNova Observatories with General Long Baseline Experiment Simulator (SNOwGLoBES). Additional routines in the SNEWPY package automate the processing of the generated data through the SNOwGLoBES software and collate its output into the observable channels of each detector. In this paper we describe the contents of the package, the physics behind SNEWPY, the organization of the code, and provide examples of how to make use of its capabilities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Three-dimensional Core-collapse Supernova Simulations with 160 Isotopic Species Evolved to Shock Breakout

We present three-dimensional simulations of core-collapse supernovae using the FLASH code that follow the progression of the explosion to the stellar surface, starting from neutrino radiation hydrodynamic simulations of the neutrino-driven phase performed with the Chimera code. We consider a 9.6 M ⊙ zero-metallicity progenitor starting from both 2D and 3D Chimera models and a 10 M ⊙ solar-metallicity progenitor starting from a 2D Chimera model, all simulated until shock breakout in 3D while tracking 160 nuclear species. The relative velocity difference between the supernova shock and the metal-rich Rayleigh–Taylor (R-T) "bullets" determines how the metal-rich ejecta evolves as it propagates through the density profile of the progenitor and dictates the final morphology of the explosion. We find maximum 56 Ni velocities of ~1950 and ~1750 km s -1 at shock breakout from 2D and 3D 9.6 M ⊙ Chimera models, respectively, due to the bullets' ability to penetrate the He/H shell. When mapping from 2D, we find that the development of higher-velocity structures is suppressed when the 2D Chimera model and 3D FLASH model meshes are aligned. The development of faster-growing spherical-bubble structures, as opposed to the slower-growing toroidal structure imposed by axisymmetry, allows for interaction of the bullets with the shock and seeds further R-T instabilities at the He/H interface. We see similar effects in the 10 M ⊙ model, which achieves maximum 56 Ni velocities of ~2500 km s -1 at shock breakout.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The Three-dimensional Collapse of a Rapidly Rotating 16 M ⊙ Star

I report on the three-dimensional (3D) hydrodynamic evolution of a rapidly rotating 16 M ⊙ star to iron core collapse. For the first time, I follow the 3D evolution of the angular momentum (AM) distribution in the iron core and convective shell burning regions for the final 10 minutes up to and including gravitational instability and core collapse. In 3D, convective regions show efficient AM transport that leads to an AM profile that differs in shape and magnitude from MESA within a few shell convective turnover timescales. For different progenitor models, such as those with tightly coupled Si/O convective shells, efficient AM transport in 3D simulations could lead to a significantly different AM distribution in the stellar interior affecting estimates of the natal neutron star or black hole spin. The results suggest that 3D AM transport in convective and rotating shell burning regions are critical components in models of massive stars and could qualitatively alter the explosion outcome and inferred compact remnant properties.

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↗
NASA NTRS1987

Supershells and propagating star formation

Stellar winds and repeated supernovae from an OB association will create a cavity of coronal gas in the interstellar medium, with radius greater than 100 pc, surrounded by a dense, expanding shell of cool interstellar gas. If the association has a typical initial mass function, its supernovae explosions will inject energy into the supershell at a nearly constant rate for about 50 Myr. The supershell loses its interior pressure and enters the snowplow phase when radiative cooling becomes important or when the shell bursts through the gas disk of a galaxy, typically after a few times 10 Myr and with a radius of 100-300 pc. At approximately the same time, the supershell becomes gravitationally unstable, forming giant molecular clouds which are sites for new star formation. There is widespread evidence for supershells in the Galaxy and other spiral and irregular galaxies from 21-cm emission-line surveys, optical emission-line surveys, and studies of supernova remnants. The gravitational instability of the supershells provides a physical mechanism for induced star formation and may account for bursts of star formation, especially in irregular galaxies.

Mccray, Richard↗
NASA NTRS2008

How Very Massive Metal-Free Stars Start Cosmological Reionization

The initial conditions and relevant physics for the formation of the earliest galaxies are well specified in the concordance cosmology. Using ab initio cosmological Eulerian adaptive mesh refinement radiation hydrodynamical calculations, we discuss how very massive stars start the process of cosmological reionization. The models include nonequilibrium primordial gas chemistry and cooling processes and accurate radiation transport in the case B approximation using adaptively ray-traced photon packages, retaining the time derivative in the transport equation. Supernova feedback is modeled by thermal explosions triggered at parsec scales. All calculations resolve the local Jeans length by at least 16 grid cells at all times and as such cover a spatial dynamic range of approx.10(exp 6). These first sources of reionization are highly intermittent and anisotropic and first photoionize the small-scale voids surrounding the halos they form in, rather than the dense filaments they are embedded in. As the merging objects form larger, dwarf-sized galaxies, the escape fraction of UV radiation decreases and the H II regions only break out on some sides of the galaxies, making them even more anisotropic. In three cases, SN blast waves induce star formation in overdense regions that were formed earlier from ionization front instabilities. These stars form tens of parsecs away from the center of their parent DM halo. Approximately five ionizing photons are needed per sustained ionization when star formation in 10(exp 6) stellar Mass halos is dominant in the calculation. As the halos become larger than approx.10(exp 7) Stellar Mass, the ionizing photon escape fraction decreases, which in turn increases the number of photons per ionization to 15-50, in calculations with stellar feedback only. Radiative feedback decreases clumping factors by 25% when compared to simulations without star formation and increases the average temperature of ionized gas to values between 3000 and 10,000 K.

Wise, John H.↗
NASA NTRS1995

Observations in the vicinity of substorm onset: Implications for the substrom process

Multi-instrument data sets from the ground and satellites at both low and high altitude have provided new results concerning substorm onset and its source region in the magnetosphere. Twenty-six out of 37 substorm onset events showed evidence of azimuthally spaced auroral forms (AAFs) prior to the explosive poleward motion associated with optical substorm onset. AAFs can span 8 hours of local time prior to onset and generally propagate eastward in the morning sector. Onset itself is, however, more localized spanning only about 1 hour local time. AAF onset occur during time periods when the solar wind pressure is relatively high. AAFs brighten in conjunction with substorm onset leading to the conclusion that they are a growth phase activity casually related to substorm onset. Precursor activity associated with these AAFs is also seen near geosynchronous orbit altitude and examples show the relationship between the various instrumental definitions of substorm onset. The implied mode number (30 to 135) derived from this work is inconsistent with cavity mode resonances but is consistent with a modified flute/ballooning instability which requires azimuthal pressure gradients. The extended source region and the distance to the open-closed field line region constrain reconnection theory and local mechanisms for substorm onset. It is demonstrated that multiple onset substorms can exist for which localized dipolarizations and the Pi 2 occur simultaneously with tail stretching existing elsewhere. These pseudobreakups can be initiated by auroral streamers which originate at the most poleward set of arc systems and drift to the more equatorward main UV oval. Observations are presented of these AAFs in conjunction with low- and high-altitutde particle and magnetic field data. These place the activations at the interface between dipolar and taillike field lines probably near the peak in the cross-tail current. These onsets are put in the context of a new scenario for substorm morphology which employs individual modules which operate independently or couple together. This allows particular substorm events to be more accurately described and investigated.

Elphinstone, R. D.↗
DOE OSTI2026

Reformulating a Gurson-based dynamic damage model and demonstrating improved predictive power and numerical robustness

The original Tepla (TEnsile PLAsticity) ductile damage model, based on the Gurson yield surface, has long been used to model damage evolution and material failure under dynamic loading. Unfortunately, Tepla suffered from mesh sensitivity, numerical instability, and limited predictive capability. Here, we theoretically reformulate Tepla to address these issues. We especially focus on the prediction of porosity, which is the key state variable used for modeling ductile damage, as compared to more easily measured surface velocities, which are at best an indirect measure of damage. Key model changes include separating the viscosity during volumetric void growth from underlying shear strength behavior and switching to an iterative bisection solver. The new Tepla is then calibrated on incipient spall experiments on half-hard copper and tantalum, which demonstrate its ability to simultaneously fit the model to recovered porosity distributions and measured surface velocities, a stringent test. Improved numerical behavior, such as greatly reduced mesh sensitivity, is also shown in those simulations. Finally, the new Tepla model is applied to several high-explosive loaded, sweeping wave experiments, showing the ability of the model to predict behavior on tests with significantly different loading conditions and histories than the calibration data.

36 MATERIALS SCIENCE↗
NASA NTRS1987

Kinetic aspects of tail dynamics - Theory and simulation

Kinetic theories relevant to the geomagnetic tail are reviewed. The topics discussed include kinetic instabilities, simulations, and current-sheet particle acceleration. Tearing mode and reconnection theories are emphasized. Kinetic treatment is appropriate for these topics since the tail plasma is collisionless. Fluid calculations are appropriate when stochastic processes dominate and for studies where long wavelengths are important. However, fluid treatments of tearing modes and reconnection require a finite resistivity in the diffusion region. Thus, although 'anomalous resistivity' can be guessed or in some cases calculated, ideally the kinetic treatment is often to be preferred. Particle motion and acceleration in the current sheet can give rise to beam-like distributions in the plasma-sheet boundary layer. Studies of current-sheet particle motion have also been used as the basis for 'kinetic' tail equilibrium models. Furthermore, quite recently current-sheet particle motion is used directly in Coroniti's explosive tail reconnection model. The 'inertial conductivity' from the equilibrium models provides the 'dissipation' necessary for reconnection.

Speiser, T. W.↗
DOE OSTI2022

Pulsating one-dimensional detonation in ammonia-hydrogen–air mixtures

In this study, one-dimensional detonations in ammonia/hydrogen-air mixtures are numerically investigated by solving the fully compressible Navier-Stokes equations with detailed chemistry. Pulsating instabilities with single-mode are observed during the detonation wave propagation, accompanied by periodic coupling and decoupling of the lead shock wave and the reactive front. The ratio between driver pressure and initial pressure determines the overdrive degree and thus the oscillatory mode of detonation for a premixture with certain composition. The effects of hydrogen dilution and mixture equivalence ratio on pulsating detonations are also examined under a constant driver pressure. The growing hydrogen fraction in fuel blends significantly increases the oscillation frequency. In addition, the pulsating detonation frequency rises with increasing equivalence ratio under fuel-lean conditions, peaks under stoichiometric conditions, and falls under fuel-rich conditions as the equivalence ratio increases further. Evolutions of reactants, main intermediate radicals, and products are analysed in both fuel-lean and fuel-rich conditions. A chemical explosive mode analysis further confirms the highly-autoignitive nature of the mixture in the induction zone between reaction front and shock front where thermal diffusion plays a negligible role.

08 HYDROGEN↗
NASA NTRS2003

How Can The SN-GRB Time Delay Be Measured?

The connection between SNe and GRBs, launched by SN 1998bw / GRB 980425 and clinched by SN 2003dh / GRB 030329-with the two GRBs differing by a factor of approximately 50000 in luminosity-so far suggests a rough upper limit of approximately 1-2 days for the delay between SN and GRB. Only four SNe have had nonnegligible coverage in close coincidence with the initial explosion, near the W shock breakout: two Qpe II, and two Type IC, SN 1999ex and SN 1998bw. For the latter, only a hint of the minimum between the UV maximum and the radioactivity bump served to help constrain the interval between SN and GRB. Swift GRB alerts may provide the opportunity to study many SNe through the UV breakout phase: GRB 980425 look dikes -apparently nearby, low- luminosity, soft-spectrum, long-lag GRBs-accounted for half of BATSE bursts near threshold, and may dominate the Swift yield near threshold, since it has sensitivity to lower energies than did BATSE. The SN to GRB delay timescale should be better constrained by prompt UV/optical observations alerted by these bursts. Definitive delay measurements may be obtained if long-lag bursters are truly nearby: The SNe/GRBs could emit gravitational radiation detectable by LIGO-II if robust non-axisymmetric bar instabilities develop during core collapse, and/or neutrino emission may be detectable as suggested by Meszaros et al.

Norris, J. P.↗
DOE OSTI2024

Modification of Jet Velocities in an Explosively Loaded Copper Target with a Conical Cavity

Here, in this work, the design and execution of an experiment with the goal of demonstrating control over the evolution of a copper jet is described. Simulations show that when using simple multi-material buffers placed between a copper target with a conical cavity and a cylinder of high-explosive, a variety of jetting behaviors occur based on material placement, including both jet velocity augmentation and mitigation. A parameter sweep was performed to determine optimal buffer designs in two configurations. Experiments using the optimal buffer designs verified the effectiveness of the buffers at altering jet velocities. Similar trends were shown between the experimental results and the modeling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2013

An Update on the Effectiveness of Postmission Disposal in LEO

The commonly‐adopted orbital debris mitigation measures were developed to reduce the growth of the future debris population. A major component in debris mitigation is post-mission disposal (PMD). The key PMD element for LEO satellites is the 25‐year rule. It is intended to limit the long‐term presence of rocket bodies (R/Bs) and spacecraft (S/C), as well as mission-related debris, in the environment. The effectiveness of PMD has been demonstrated and documented since the development of mitigation measures began in the 1990s. This paper summarizes an updated study, based on the current environment, using the NASA LEGEND model. The study focused on the > or = 10 cm population in LEO. The historical simulation covered 1957 through 2011 and followed the recorded launches and known breakup events. The future projection was carried out for 200 years. An eight‐year launch traffic, 2004 - 2011, was repeated during the projection period. An eight‐year mission lifetime was assumed for future S/C. No stationkeeping and no collision avoidance maneuver were implemented. Only objects 10 cm and larger were included in collision consideration. No explosion was allowed for R/Bs and S/C launched after 2011. The 25‐year PMD rule success rates were set at 0%, 10%, 50%, 75%, and 95%, respectively, for the 5 study scenarios. Results of the simulations were analyzed to quantify the differences among the different compliance rates. As expected, the 0% PMD projection followed a rapid and non‐linear increase in the next 200 years. The LEO population, on average, more than tripled at the end of the simulations. With a 50% compliance of the 25‐year rule, the population growth was reduced approximately by half. However, even with a 95% compliance of the 25‐year rule, the LEO debris population would still increase by an average of more than 50% in 200 years. These simulation results provide an updated assessment of the effectiveness of the 25‐year rule. It is the first and the most cost‐effective defense against future population growth. In addition, the results also confirm the instability of the LEO population and lay the foundation for the need to consider environment remediation in the future.

Liou, J.-C.↗
NASA NTRS1985

Plasma waves associated with the AMPTE artificial comet

Numerous plasma wave effects were detected by the AMPTE/IRM spacecraft during the artificial comet experiment on December 27, 1984. As the barium ion cloud produced by the explosion expanded over the spacecraft, emissions at the electron plasma frequency and ion plasma frequency provided a determination of the local electron density. The electron density in the diamagnetic cavity produced by the ion cloud reached a peak of more than 5 x 10 to the 5th per cu cm, then decayed smoothly as the cloud expanded, varying approximately as t exp-2. As the cloud began to move due to interactions with the solar wind, a region of compressed plasma was encountered on the upstream side of the diamagnetic cavity. The peak electron density in the compression region was about 1.5 x 10 to the 4th per cu cm. Later, a very intense (140 mVolt/m) broadband burst of electrostatic noise was encountered on the sunward side of the compression region. This noise has characteristics very similar to noise observed in the earth's bow shock, and is believed to be a shocklike interaction produced by an ion beam-plasma instability between the nearly stationary barium ions and the streaming solar wind protons.

Gurnett, D. A.↗
NASA NTRS1993

High-resolution mid-infrared spectra of Co II, Ni I, and Fe II in SN 1987A

Ground-based infrared observations of SN 1987A on day 612 after the explosion have yielded resolved line profiles of Co II, Ni I, Fe II at 10.52, 11.31, and 17.94 micron, respectively. The spectra were taken at a resolving power of about 1000 with an array grating spectrometer on the 4 m telescope of Cerro Tololo Inter-American Observatory. Based on the observed line intensities we have estimated the minimum mass of each ion: M(Co II) = (6.0 +/- 1.8) x 10 exp -5 solar mass; M(Ni I) = (1.1 +/- 0.1) x 10 exp -3 solar mass; and M(Fe II) = (8.0 +/- 1.5) x 10 exp -3 solar mass. From these we infer total masses for cobalt, nickel, and iron in the ejecta. The nickel and iron line profiles are markedly asymmetric. We interpret these as arising from two components, one centered on the stellar rest velocity with an approximately 3250 km/s full width, and the second at about +1200 km/s with an approximately 1100 km/s full width. The asymmetry may represent a large-scale fracturing of the ejecta by Rayleigh-Taylor instabilities.

Jennings, D. E.↗
DOE OSTI2022

Spatially resolved gas and stellar kinematics in compact starburst galaxies

The kinematics of galaxies provide valuable insights into their physics and assembly history. Kinematics are governed not only by the gravitational potential, but also by merger events and stellar feedback processes such as stellar winds and supernova explosions. We aim to identify what governs the kinematics in a sample of SDSS-selected nearby starburst galaxies, by obtaining spatially resolved measurements of the gas and stellar kinematics. We obtained near-infrared integral-field K-band spectroscopy with VLT/SINFONI for 15 compact starburst galaxies. We derived the integrated as well as spatially resolved stellar and gas kinematics. The stellar kinematics were derived from the CO absorption bands, and Paα and Brγ emission lines were used for the gas kinematics. Based on the integrated spectra, we find that the majority of galaxies have gas and stellar velocity dispersion that are comparable. A spatially resolved comparison shows that the six galaxies that deviate show evidence for a bulge or stellar feedback. Two galaxies are identified as mergers based on their double-peaked emission lines. In our sample, we find a negative correlation between the ratio of the rotational velocity over the velocity dispersion (v rot /σ) and the star formation rate surface density. We propose a scenario where the global kinematics of the galaxies are determined by gravitational instabilities that affect both the stars and gas. This process could be driven by mergers or accretion events. Effects of stellar feedback on the ionised gas are more localised and detected only in the spatially resolved analysis. The mass derived from the velocity dispersion provides a reliable mass even if the galaxy cannot be spatially resolved. The technique used in this paper is applicable to galaxies at low and high redshift with the next generation of infrared-focussed telescopes (JWST and ELT).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

General relativistic pulsations of ultra-massive ZZ Ceti stars

ABSTRACT Ultra-massive white dwarf stars are currently being discovered at a considerable rate, thanks to surveys such as the Gaia space mission. These dense and compact stellar remnants likely play a major role in Type Ia supernova explosions. It is possible to probe the interiors of ultra-massive white dwarfs through asteroseismology. In the case of the most massive white dwarfs, general relativity could affect their structure and pulsations substantially. In this work, we present results of relativistic pulsation calculations employing relativistic ultra-massive ONe-core white dwarf models with hydrogen-rich atmospheres and masses ranging from 1.29 to $1.369 \ \mathrm{M}_{\odot }$ with the aim of assessing the impact of general relativity on the adiabatic gravity (g)-mode period spectrum of very high mass ZZ Ceti stars. Employing the relativistic Cowling approximation for the pulsation analysis, we find that the critical buoyancy (Brunt–Väisälä) and acoustic (Lamb) frequencies are larger for the relativistic case, compared to the Newtonian case, due to the relativistic white dwarf models having smaller radii and higher gravities for a fixed stellar mass. In addition, the g-mode periods are shorter in the relativistic case than those in the Newtonian computations, with relative differences of up to ∼$50$ per cent for the highest mass models ($1.369 \ \mathrm{M}_{\odot }$) and for effective temperatures typical of the ZZ Ceti instability strip. Hence, the effects of general relativity on the structure, evolution, and pulsations of white dwarfs with masses larger than ∼$1.29 \ \mathrm{M}_{\odot }$ cannot be ignored in the asteroseismological analysis of ultra-massive ZZ Ceti stars.

79 ASTRONOMY AND ASTROPHYSICS↗