Abstract:We analyze the linear stability of a stalled accretion shock in a perfect gas with a parametrized cooling function L ~ rho^{beta-alpha} P^alpha. The instability is dominated by the l=1 mode if the shock radius exceeds 2-3 times the accretor radius, depending on the parameters of the cooling function. The growth rate and oscillation period are comparable to those observed in the numerical simulations of Blondin & Mezzacappa (2006). The instability mechanism is analyzed by separately measuring the efficiencies of the purely acoustic cycle and the advective-acoustic cycle. These efficiencies are estimated directly from the eigenspectrum, and also through a WKB analysis in the high frequency limit. Both methods prove that the advective-acoustic cycle is unstable, and that the purely acoustic cycle is stable. Extrapolating these results to low frequency leads us to interpret the dominant mode as an advective-acoustic instability, different from the purely acoustic interpretation of Blondin & Mezzacappa (2006). A simplified characterization of the instability is proposed, based on an advective-acoustic cycle between the shock and the radius r_nabla where the velocity gradients of the stationary flow are strongest. The importance of the coupling region in this mechanism calls for a better understanding of the conditions for an efficient advective-acoustic coupling in a decelerated, nonadiabatic flow, in order to extend these results to core-collapse supernovae.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the linear stability problem of the stalled accretion shock, especially in the context of core - collapse supernovae. Specifically, the author analyzes the linear stability of the stalled accretion shock with a parameterized cooling function \( L \propto \rho^{\beta-\alpha} P^{\alpha} \) in a perfect gas. ### Main problems 1. **Determine the instability mechanism**: The author attempts to understand and distinguish the roles of the acoustic cycle and the advective - acoustic cycle in the instability mechanism. 2. **Verify the numerical simulation results**: Through linear stability analysis, confirm whether the dominant \( l = 1 \) mode observed in the numerical simulation by Blondin & Mezzacappa (2006) indeed corresponds to the most unstable eigenmode. 3. **Explain the influence of different cooling functions**: Study the influence of different cooling parameters (such as \( \alpha \) and \( \beta \)) on instability and characteristic frequencies, especially whether the mode \( l = 1 \) is always the most unstable when the shock radius is large. ### Key findings - **Instability of the advective - acoustic cycle**: Through direct measurement of the characteristic spectrum and WKB analysis, it is proved that the advective - acoustic cycle is unstable, while the pure acoustic cycle is stable. - **Dominance of the mode \( l = 1 \)**: When the shock radius exceeds 2 - 3 times the radius of the accretor, the mode \( l = 1 \) is dominant. - **Consistency with numerical simulation**: The growth rate and oscillation period of the characteristic frequencies are consistent with the numerical simulation results of Blondin & Mezzacappa (2006). ### Formula presentation The cooling function is defined as: \[ L \propto \rho^{\beta-\alpha} P^{\alpha} \] The sound speed \( c \) is defined as: \[ c^2 \equiv \frac{\gamma P}{\rho} \] The differential system of the eigenvalue problem includes the following equations: \[ \frac{\partial f}{\partial r} = \delta\left(\frac{L}{\rho v}\right) + i \omega \frac{v}{1 - M^2} \left( h - \frac{f}{c^2} + \left[\gamma - 1+\frac{1}{M^2}\right] \frac{\delta S}{\gamma} \right) \] \[ \frac{\partial h}{\partial r} = i \frac{\delta K}{\omega r^2 v} + i \omega \frac{1}{v(1 - M^2)} \left( \mu^2 \frac{f}{c^2} - M^2 h - \delta S \right) \] \[ \frac{\partial \delta S}{\partial r} = i \omega \frac{\delta S}{v} + \delta\left(\frac{L}{P v}\right) \] \[ \frac{\partial \delta K}{\partial r} = i \omega \frac{\delta K}{v} + l(l + 1) \delta\left(\frac{L}{\rho v}\right) \] where: \[ \mu^2 \equiv 1-\frac{l(l + 1)c^2}{\omega^2 r^2(1 - M^2)} \] These formulas are used to describe the differential system of the perturbed flow.

Instability of a stalled accretion shock: evidence for the advective-acoustic cycle

Non-radial instabilities of isothermal Bondi accretion with a shock: vortical-acoustic cycle vs post-shock acceleration

Multidimensional supernova simulations with approximative neutrino transport. II. Convection and the advective-acoustic cycle in the supernova core

Analytic insight into the physics of SASI I. Shock instability in a non-rotating stellar core

High Mach-number Bondi--Hoyle--Lyttleton flow around a small accretor

Non-Axisymmetric g-Mode and p-Mode Instability in a Hydrodynamic Thin Accretion Disk

A Shallow Water Analogue of the Standing Accretion Shock Instability: Experimental Demonstration and Two-Dimensional Model

Stable and Unstable Accretion Flows with Angular Momentum near a Point Mass

The Local Stability of Accretion Disks with Advection

Radiative cooling instability in 1D colliding flows

HYDRODYNAMICS OF CORE-COLLAPSE SUPERNOVAE AT THE TRANSITION TO EXPLOSION. I. SPHERICAL SYMMETRY

Adiabatic Perturbations in Homologous Conventional Polytropic Core Collapses of a Spherical Star

A New Mechanism for Core-Collapse Supernova Explosions

Nonlinear Evolution of Instabilities Between Dust and Sound Waves

A shared accretion instability for black holes and neutron stars

Perturbation Analysis of a General Polytropic Homologously Collapsing Stellar Core

Three-dimensional Hydrodynamic Instabilities in Stellar Core Collapses

The Radial-Azimuthal Instability of Gas-Pressure-Dominated Accretion Disk with Advection

Thermal instabilities in accretion disks II: Numerical Experiments for the Goldreich-Schubert-Fricke Instability and the Convective Overstability in disks around young stars

Weak Shock Propagation with Accretion. III. A Numerical Study On Shock Propagation & Stability