Core-collapse supernovae span a wide range of energies, from much less than to much greater than the binding energy of the progenitor star. As a result, the shock wave generated from a supernova explosion can have a wide range of Mach numbers. In this paper, we investigate the propagation of shocks with arbitrary initial strengths in polytropic stellar envelopes using a suite of spherically symmetric hydrodynamic simulations. We interpret these results using the three known self-similar solutions for this problem: the Sedov-Taylor blastwave describes an infinitely strong shock and the self-similar solutions from Coughlin et al. (2018b) (Paper I) and Coughlin et al. (2019) (Paper II) describe a weak and infinitely weak shock (the latter being a rarefaction wave). We find that shocks, no matter their initial strengths, evolve toward either the infinitely strong or infinitely weak self-similar solutions at sufficiently late times. For a given density profile, a single function characterizes the long-term evolution of a shock's radius and strength. However, shocks with strengths near the self-similar solution for a weak shock (from Paper I) evolve extremely slowly with time. Therefore, the self-similar solutions for infinitely strong and infinitely weak shocks are not likely to be realized in low-energy stellar explosions, which will instead retain memory of the shock strength initiated in the stellar interior.

Early-type Wolf-Rayet(WR) stellar models harbor a super-Eddington layer in their outer envelopes due to a prominent iron opacity bump. In the past few decades, one-dimensional hydrostatic and time-steady hydrodynamic models have suggested a variety of WR responses to a super-Eddington force including envelope inflation and optically thick winds. In this paper, I study these responses using semianalytical estimates and WR models from both MESA and Ro & Matzner; four conclusions are present. First, early-type WR stars do not harbor inflated envelopes because they have either strong winds or insufficient luminosities. Second, the condition for an opacity bump to harbor a sonic point is expressible as a minimum mass-loss rate, $\dot{M}_{\rm sp}(L_*)$. In agreement with Grassitelli et al. and Ro, the majority of galactic early-type WR stars can harbor sonic points at the iron opacity bump. However, about half of those in the Large Magellanic Cloud cannot given typical wind parameters. Third, WR winds driven by the iron opacity bump must have mass-loss rates that exceed a global minimum of $10^{-5.8}-10^{-6}M_\odot\,{\rm yr}^{-1}$. Lastly, the observed early-type WR distribution follows a simple mass-loss relation derived here if the radiation-to-gas pressure ratio is approximately $p_r/p_g \sim 145$ in the wind; a value consistent with studies by Grafener et al. and Nakauchi & Saio.

Coughlin et al. (2018) (Paper I) derived and analyzed a new regime of self-similarity that describes weak shocks (Mach number of order unity) in the gravitational field of a point mass. These solutions are relevant to low energy explosions, including failed supernovae. In this paper, we develop a formalism for analyzing the stability of shocks to radial perturbations, and we demonstrate that the self-similar solutions of Paper I are extremely weakly unstable to such radial perturbations. Specifically, we show that perturbations to the shock velocity and post-shock fluid quantities (the velocity, density, and pressure) grow with time as $t^{\alpha}$, where $\alpha \le 0.12$, implying that the ten-folding timescale of such perturbations is roughly ten orders of magnitude in time. We confirm these predictions by performing high-resolution, time-dependent numerical simulations. Using the same formalism, we also show that the Sedov-Taylor blastwave is trivially stable to radial perturbations provided that the self-similar, Sedov-Taylor solutions extend to the origin, and we derive simple expressions for the perturbations to the post-shock velocity, density, and pressure. Finally, we show that there is a third, self-similar solution (in addition to the the solutions in Paper I and the Sedov-Taylor solution) to the fluid equations that describes a rarefaction wave, i.e., an outward-propagating sound wave of infinitesimal amplitude. We interpret the stability of shock propagation in light of these three distinct self-similar solutions.

We present solutions for the self-similar propagation of a shock wave in a hydrostatic, adiabatic medium with a point mass gravitational field. In contrast to the well-known, Sedov-Taylor blastwave, these solutions apply to the case when the shock Mach number is of order a few, and the energy of the shocked fluid is not conserved but self-consistently modified by the binding energy of the ambient medium that is swept up by the passage of the shock. Furthermore, we show that there is one solution (for a given ambient density profile) that smoothly passes through a sonic point in the post-shock flow and results in accretion onto the central object; in analogy with the Bondi problem, we propose that these solutions are the ones that are most relevant in astrophysical environments. We apply these accreting models to failed supernovae, in which neutron star formation does not unbind the envelope, but a weak shock is still generated in the outer layers of the star from neutrino-induced mass loss. We find excellent agreement between the predictions of our self-similar, shock propagation model and numerical simulations of the collapse of a yellow supergiant; the self-similar solutions reproduce the overall scaling of the shock speed, the time and space-dependent evolution of the velocity, density, and pressure behind the shock, and the accretion rate onto the black hole. Our results have important implications for the fallback and ejection of material in failed supernovae.

Pre-supernova (SN) outbursts from massive stars may be driven by hydrodynamical wave energy emerging from the core of the progenitor star during late nuclear burning phases. Here, we examine the effects of wave heating in stars containing little or no hydrogen, i.e., progenitors of type IIb/Ib SNe. Because there is no massive hydrogen envelope, wave energy is thermalized near the stellar surface where the overlying atmospheric mass is small but the optical depth is large. Wave energy can thus unbind this material, driving an optically thick, super-Eddington wind. Using 1D hydodynamic MESA simulations of $\sim \! 5 \, M_\odot$ He stars, we find that wave heating can drive pre-SN outbursts composed of a dense wind whose mass loss rate can exceed $\sim \! 0.1 \, M_\odot/{\rm yr}$. The wind terminal velocities are a few $100 \, {\rm km}/{\rm s}$, and outburst luminosities can reach $\sim \! 10^6 \, L_\odot$. Wave-driven outbursts may be linked with observed or inferred pre-SN outbursts of type Ibn/transitional/transformational SNe, and pre-SN wave-driven mass loss is a good candidate to produce these types of SNe. However, we also show that non-linear wave breaking in the core of the star may prevent such outbursts in stars with thick convective helium-burning shells. Hence, only a limited subset of SN progenitors are likely to experience wave-driven pre-SN outbursts.

Wave-driven outflows and non-disruptive explosions have been implicated in pre-supernova outbursts, supernova impostors, LBV eruptions, and some narrow-line and superluminous supernovae. To model these events, we investigate the dynamics of stars set in motion by strong acoustic pulses and wave trains, focusing here on nonlinear wave propagation, shock formation, and an early phase of the development of a weak shock. We identify the shock formation radius, showing that a heuristic estimate based on crossing characteristics matches an exact expansion around the wave front and verifying both with numerical experiments. Our general analytical condition for shock formation applies to one-dimensional motions within any static environment, including both eruptions and implosions, and can easily be extended to non-stationary flows. We also consider the early phase of shock energy dissipation. We find that waves of super-Eddington acoustic luminosity always create shocks, rather than damping by radiative diffusion. Therefore, shock formation is integral to super-Eddington outbursts.

Hydrostatic models of Wolf-Rayet stars typically contain low-density outer envelopes that inflate the stellar radii by a factor of several and are capped by a denser shell of gas. Inflated envelopes and density inversions are hallmarks of envelopes that become super-Eddington as they cross the iron-group opacity peak, but these features disappear when mass loss is sufficiently rapid. We re-examine the structures of steady, spherically symmetric wind solutions that cross a sonic point at high optical depth, identifying the physical mechanism by which outflow affects the stellar structure, and provide an improved analytical estimate for the critical mass loss rate above which extended structures are erased. Weak-flow solutions below this limit resemble hydrostatic stars even in supersonic zones; however, we infer that these fail to successfully launch optically thick winds. Wolf-Rayet envelopes will therefore likely correspond to the strong, compact solutions. We also find that wind solutions with negligible gas pressure are stably stratified at and below the sonic point. This implies that convection is not the source of variability in Wolf-Rayet stars, as has been suggested; but, acoustic instabilities provide an alternative explanation. Our solutions are limited to high optical depths by our neglect of Doppler enhancements to the opacity, and do not account for acoustic instabilities at high Eddington factors; yet they provide useful insights into Wolf-Rayet stellar structures.

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.

In a non-spherical stellar explosion, non-radial motions become important near the stellar surface. For realistic deviations from spherical symmetry, non-radial flow dramatically alters the dynamics and emission of shock emergence on a significant fraction of the surface. The breakout flash is stifled, ejecta speeds are limited, and matter is cast sideways. Non-radial ejection allows for collisions outside the star, which may engender a new type of transient. Strongly oblique breakouts are most easily produced in compact stellar progenitors, such as white dwarfs and stripped-envelope core collapse supernovae. We study the shock structure and post-shock acceleration using conservation laws, a similarity analysis, and an approximate theory for oblique shocks. The shock is likely to extend vertically from the stellar surface, then kink before joining a deep asymptotic solution. Outflow from the region crossed by an oblique shock is probably unsteady and may affect the surface ahead of the main shock. We comment on the implications for several notable explosions in which the non-spherical dynamics described in this paper are likely to play an important role. We also briefly consider relativistic and superluminal pattern speeds.

We examine the dynamics of accelerating normal shocks in stratified planar atmospheres, providing accurate fitting formulae for the scaling index relating shock velocity to the initial density and for the post-shock acceleration factor as functions of the polytropic and adiabatic indices which parameterize the problem. In the limit of a uniform initial atmosphere there are analytical formulae for these quantities. In the opposite limit of a very steep density gradient the solutions match the outcome of shock acceleration in exponential atmospheres.

The existence of an extended hot gaseous corona surrounding clusters, groups and massive galaxies is well established by observational evidence and predicted by current theories of galaxy formation. When a small galaxy collides with a larger one, their coronae are the first to interact, producing disturbances that remove gas from the smaller system and settle it into the corona of the larger one. For a Milky-Way-size galaxy merging into a low-mass group, ram pressure stripping and the Kelvin-Helmholtz instability are the most relevant of these disturbances. We argue that the turbulence generated by the latter mixes the material of both coronae in the wake of the orbiting satellite creating a "warm phase" mixture with a cooling time a factor of several shorter than that of the ambient intragroup gas. We reach this conclusion using analytic estimates, as well as adiabatic and dissipative high resolution numerical simulations of a spherical corona subject to the ablation process of a constant velocity wind with uniform density and temperature. Although this is a preliminary analysis, our results are promising and we speculate that the mixture could potentially trigger in situ star formation and/or be accreted into the central galaxy as a cold gas flow resulting in a new mode of star formation in galaxy groups and clusters.