The Expulsion of Stellar Envelopes in Core-Collapse Supernovae

We examine the relation between presupernova stellar structure and the distribution of ejecta in core-collapse supernovae, assuming adiabatic, spherically symmetric flow. We develop a simple yet accurate formula for the blastwave shock velocity, and demonstrate that the entire final density distribution can be approximated with simple models for the final pressure distribution, along with the approximate shock-deposited entropy, in a way that matches the results of simulations. We find that the distribution of density in a star's ejecta depends on whether its outer envelope is radiative or convective, and if convective, on the composition structure of the star; simple approximate forms are presented for red and blue supergiant ejecta. Our models are most accurate for the high-velocity ejecta from the periphery of a star, where the shock dynamics are predictable. We present formulae for the final density distribution of this material, for both radiative and efficiently convective envelopes. These formulae limit to the well-known planar, self-similar solutions for mass shells approaching the stellar surface. But, the assumption of adiabatic flow fails at low optical depth, so this planar limit need not be attained. Formulae are given for the observable properties of the X-ray burst accompanying shock emergence, and their dependence on the parameters of the explosion. Motivated by the relativistic expansion recently inferred by Kulkarni et al. (1998) for the synchrotron shell around SN1998bw, we estimate the criterion for relativistic mass ejection and the rest mass of relativistic ejecta.Comment: 57 pages, 10 eps figures, aaspp4, submitted to Ap

Oblique Shock Breakout in Supernovae and Gamma-Ray Bursts: II. Numerical Solutions For Non-Relativistic Pattern Speeds

Non-spherical explosions develop non-radial flows as the pattern of shock emergence progresses across the stellar surface. In supernovae these flows can limit ejecta speeds, stifle shock breakout emission, and cause collisions outside the star. Similar phenomena occur in stellar and planetary collisions, tidal disruption events, accretion-induced collapses, and propagating detonations. We present two-dimensional, nested-grid Athena simulations of non-radial shock emergence in a frame comoving with the breakout pattern, focusing on the adiabatic, non-relativistic limit in a plane stratified envelope. We set boundary conditions using a known self-similar solution and explore the role of box size and resolution on the result. The shock front curves toward the stellar surface, and exhibits a kink from which weak discontinuities originate. Flow around the point of shock emergence is neither perfectly steady nor self-similar. Waves and vortices, which are not predominantly due to grid effects, emanate from this region. The post-shock flow is deflected along the stellar surface, and its pressure disturbs the stellar atmosphere upstream of the emerging shock. We use the numerical results and their analytical limits to predict the effects of radiation transfer and gravity, which are not included in our simulations.Comment: 15 pages, 12 figures, submitted to Ap