Nonlinear Evolution of Instabilities Between Dust and Sound Waves
Abstract
We study the non-linear evolution of the acoustic 'Resonant Drag Instability' (RDI) using numerical simulations. The acoustic RDI is excited in a dust-gas mixture when dust grains stream through gas, interacting with sound waves to cause a linear instability. We study this process in a periodic box by accelerating neutral dust with an external driving force. The instability grows as predicted by linear theory, eventually breaking into turbulence and saturating. As in linear theory, the non-linear behavior is characterized by three regimes - high, intermediate, and low wavenumbers - the boundary between which is determined by the dust-gas coupling strength and the dust-to-gas mass ratio. The high and intermediate wavenumber regimes behave similarly to one another, with large dust-to-gas ratio fluctuations while the gas remains largely incompressible. The saturated state is highly anisotropic: dust is concentrated in filaments, jets, or plumes along the direction of acceleration, with turbulent vortex-like structures rapidly forming and dissipating in the perpendicular directions. The low-wavenumber regime exhibits large fluctuations in gas and dust density, but the dust and gas remain more strongly coupled in coherent 'fronts' perpendicular to the acceleration. These behaviors are qualitatively different from those of dust 'passively' driven by external hydrodynamic turbulence, with no back-reaction force from dust onto gas. The virulent nature of these instabilities has interesting implications for dust-driven winds in a variety of astrophysical systems, including around cool-stars, in dusty torii around active-galactic-nuclei, and in and around giant molecular clouds.
Additional Information
We would like to thank Stefania Moroianu for a number of helpful discussions. Support for PFH, JS, & ERM was provided by an Alfred P. Sloan Research Fellowship, NSF Collaborative Research Grant #1715847 and CAREER grant #1455342, and NASA grants NNX15AT06G, JPL 1589742, 17-ATP17-0214. Numerical calculations were run on the Caltech compute cluster "Wheeler," allocations from XSEDE TG-AST130039 and PRAC NSF.1713353 supported by the NSF, and NASA HEC SMD-16-7592.Attached Files
Submitted - 1810.08214.pdf
Files
Name | Size | Download all |
---|---|---|
md5:aec892de5c634cd3e1a397d3a516cffb
|
10.8 MB | Preview Download |
Additional details
- Eprint ID
- 92722
- Resolver ID
- CaltechAUTHORS:20190206-105610955
- Alfred P. Sloan Foundation
- NSF
- AST-1715847
- NSF
- AST-1455342
- NASA
- NNX15AT06G
- JPL
- 1589742
- JPL
- 17-ATP17-0214
- NSF
- TG-AST130039
- NSF
- OAC-1713353
- NASA
- SMD-16-7592
- Created
-
2019-02-08Created from EPrint's datestamp field
- Updated
-
2023-06-02Created from EPrint's last_modified field
- Caltech groups
- TAPIR, Astronomy Department