Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userInstabilities of Spiral Shocks I: Onset of Wiggle Instability and its Mechanism
95
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We found that loosely wound spiral shocks in an isothermal gas disk caused by a non-axisymmetric potential are hydrodynamically unstable, if the shocks are strong enough. High resolution, global hydrodynamical simulations using three different numerical schemes, i.e. AUSM, CIP, and SPH, show similarly that trailing spiral shocks with the pitch angle of larger than ~10 deg wiggle, and clumps are developed in the shock-compressed layer. The numerical simulations also show clear wave crests that are associated with ripples of the spiral shocks. The spiral shocks tend to be more unstable in a rigidly rotating disk than in a flat rotation. This instability could be an origin of the secondary structures of spiral arms, i.e. the spurs/fins, observed in spiral galaxies. In spite of this local instability, the global spiral morphology of the gas is maintained over many rotational periods. The Kelvin-Helmholtz (K-H) instability in a shear layer behind the shock is a possible mechanism for the wiggle instability. The Richardson criterion for the K-H stability is expressed as a function of the Mach number, the pitch angle, and strength of the background spiral potential. The criterion suggests that spiral shocks with smaller pitch angles and smaller Mach numbers would be more stable, and this is consistent with the numerical results.
rate research
Read More
Gas in disk galaxies interacts nonlinearly with an underlying stellar spiral potential to form galactic spiral shocks. While numerical simulations typically show that spiral shocks are unstable to wiggle instability (WI) even in the absence of magnetic fields and self-gravity, its physical nature has remained uncertain. To clarify the mechanism behind the WI, we conduct a normal-mode linear stability analysis as well as nonlinear simulations assuming that the disk is isothermal and infinitesimally thin. We find that the WI is physical, originating from the generation of potential vorticity at a deformed shock front, rather than Kelvin-Helmholtz instabilities as previously thought. Since gas in galaxy rotation periodically passes through the shocks multiple times, the potential vorticity can accumulate successively, setting up a normal mode that grows exponentially with time. Eigenfunctions of the WI decay exponentially downstream from the shock front. Both shock compression of acoustic waves and a discontinuity of shear across the shock stabilize the WI. The wavelength and growth time of the WI depend on the arm strength quite sensitively. When the stellar-arm forcing is moderate at 5%, the wavelength of the most unstable mode is about 0.07 times the arm-to-arm spacing, with the growth rate comparable to the orbital angular frequency, which is found to be in good agreement with the results of numerical simulations.
We argue that self-excited instabilities are the cause of spiral patterns in simulations of unperturbed stellar discs. In previous papers, we have found that spiral patterns were caused by a few concurrent waves, which we claimed were modes. The superposition of a few steadily rotating waves inevitably causes the appearance of the disc to change continuously, and creates the kind of shearing spiral patterns that have been widely reported. Although we have found that individual modes last for relatively few rotations, spiral activity persists because fresh instabilities appear, which we suspected were excited by the changes to the disc caused by previous disturbances. Here we confirm our suspicion by demonstrating that scattering at either of the Lindblad resonances seeds a new groove-type instability. With this logical gap closed, our understanding of the behaviour in the simulations is almost complete. We believe that our robust mechanism is a major cause of spiral patterns in the old stellar discs of galaxies, including the Milky Way where we have previously reported evidence for resonance scattering in the recently released Gaia data.
Using one-dimensional hydrodynamic simulations including interstellar heating, cooling, and thermal conduction, we investigate nonlinear evolution of gas flow across galactic spiral arms. We model the gas as a non-self-gravitating, unmagnetized fluid, and follow its interaction with a stellar spiral potential in a local frame comoving with the stellar pattern. Initially uniform gas rapidly separates into warm and cold phases as a result of thermal instability (TI), and also forms a quasi-steady shock that prompts phase transitions. After saturation, the flow follows a recurring cycle: warm and cold phases in the interarm region are shocked and immediately cool to become a denser cold medium in the arm; post-shock expansion reduces the mean density to the unstable regime in the transition zone and TI subsequently mediates evolution back into warm and cold interarm phases. For our standard model with n_0 = 2 cm^-3, the gas resides in the dense arm, thermally-unstable transition zone, and interarm region for 14%, 22%, 64% of the arm-to-arm crossing time. These regions occupy 1%, 16%, and 83% of the arm-to-arm distance, respectively. Gas at intermediate temperatures represents ~25-30% of the total mass, similar to the fractions estimated from HI observations. Despite transient features and multiphase structure, the time-averaged shock profiles can be matched to that of a diffusive isothermal medium with temperature 1,000 K and particle mean free path of l_0 = 100 pc. Finally, we quantify numerical conductivity associated with translational motion of phase-separated gas on the grid, and show that convergence of numerical results requires the numerical conductivity to be comparable to or smaller than the physical conductivity. (Abridged)
In this paper we study the feathering substructures along spiral arms by considering the perturbational gas response to a spiral shock. Feathers are density fluctuations that jut out from the spiral arm to the inter-arm region at pitch angles given by the quantum numbers of the doubly-periodic structure. In a localized asymptotic approximation, related to the shearing sheet except that the inhomogeneities occur in space rather than in time, we derive the linearized perturbation equations for a razor-thin disk with turbulent interstellar gas, frozen-in magnetic field, and gaseous self-gravity. Apart from the modal quantum numbers, the individual normal modes of the system depend on seven dimensionless quantities that characterize the underlying time-independent axisymmetric state plus its steady, nonlinear, two-armed spiral-shock (TASS) response to a hypothesized background density-wave supported by the disk stars of the galaxy. We show that some of these normal modes have positive growth rates. Their over-density contours in the post-shock region are very reminiscent of observed feathering substructures in full magnetohydrodynamic (MHD) simulations. The feathering substructures are parasitic instabilities intrinsic to the system; thus, their study not only provides potential diagnostics for important parameters that characterize the interstellar medium of external galaxies, but also yields a deeper understanding of the basic mechanism that drives the formation of the giant molecular clouds (GMCs) and the OB stars that outline observed grand-design spirals.
We present a study of the spiral responses in a stable disc galaxy model to co-orbiting perturbing masses that are evenly spaced around rings. The amplitudes of the responses, or wakes, are proportional to the masses of the perturbations, and we find that the response to a low-mass ring disperses when it is removed -- behaviour that is predicted by linear theory. Higher mass rings cause nonlinear changes through scattering at the major resonances, provoking instabilities that were absent before the scattering took place. The separate wake patterns from two rings orbiting at differing frequencies, produce a net response that is an apparently shearing spiral. When the rings have low mass, the evolution of the simulation is both qualitatively and quantitatively reproduced by linear superposition of the two separate responses. We argue that apparently shearing transient spirals in simulations result from the superposition of two or more steadily rotating patterns, each of which is best accounted for as a normal mode of the non-smooth disc.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
