We examine the nonlinear development of unstable magnetosonic waves driven by a background radiative flux -- the Radiation-Driven Magneto-Acoustic Instability (RMI, a.k.a. the photon bubble instability). The RMI may serve as a persistent source of de
nsity, radiative flux, and magnetic field fluctuations in stably-stratified, optically-thick media. The conditions for instability are present in a variety of astrophysical environments, and do not require the radiation pressure to dominate or the magnetic field to be strong. Here we numerically study the saturation properties of the RMI, covering three orders of magnitude in the relative strength of radiation, magnetic field, and gas energies. Two-dimensional, time-dependent radiation-MHD simulations of local, stably-stratified domains are conducted with Zeus-MP in the optically-thick, highly-conducting limit. Our results confirm the theoretical expectations of Blaes and Socrates (2003) in that the RMI operates even in gas pressure-dominated environments that are weakly magnetized. The saturation amplitude is a monotonically increasing function of the ratio of radiation to gas pressure. Keeping this ratio constant, we find that the saturation amplitude peaks when the magnetic pressure is comparable to the radiation pressure. We discuss the implications of our results for the dynamics of magnetized stellar envelopes, where the RMI should act as a source of sub-photospheric perturbations.
We examine the stability of a standing shock wave within a spherical accretion flow onto a gravitating star, in the context of core-collapse supernova explosions. Our focus is on the effect of nuclear dissociation below the shock on the linear growth
, and non-linear saturation, of non-radial oscillations of the shocked fluid. We combine two-dimensional, time-dependent hydrodynamic simulations using FLASH2.5 with a solution to the linear eigenvalue problem, and demonstrate the consistency of the two approaches. Previous studies of this `Standing Accretion Shock Instability (SASI) have focused either on zero-energy accretion flows without nuclear dissociation, or made use of a detailed finite-temperature nuclear equation of state and included strong neutrino heating. Our main goal in this and subsequent papers is to introduce equations of state of increasing complexity, in order to isolate the various competing effects. In this work we employ an ideal gas equation of state with a constant rate of nuclear dissociation below the shock, and do not include neutrino heating. We find that a negative Bernoulli parameter below the shock significantly lowers the real frequency, growth rate, and saturation amplitude of the SASI. A decrease in the adiabatic index has similar effects. The non-linear development of the instability is characterized by an expansion of the shock driven by turbulent kinetic energy at nearly constant internal energy. Our results also provide further insight into the instability mechanism: the rate of growth of a particular mode is fastest when the radial advection time from the shock to the accretor overlaps with the period of a standing lateral sound wave. The fastest-growing mode can therefore be modified by nuclear dissociation.
