Polydisperse streaming instability - II. Methods for solving the linear stability problem

doi:10.1093/mnras/stab111

2021-01-19

Sijme-Jan Paardekooper *

Colin P. McNally

Francesco Lovascio

DOI：10.1093/mnras/stab111

原文链接

原文求助

摘要

En 中文

Occurring in protoplanetary discs composed of dust and gas, streaming instabilities are a favoured mechanism to drive the formation of planetesimals. The polydispserse streaming instability is a generalization of the streaming instability to a continuum of dust sizes. This second paper in the series provides a more in-depth derivation of the governing equations and presents novel numerical methods for solving the associated linear stability problem. In addition to the direct discretization of the eigenproblem at second order introduced in the previous paper, a new technique based on numerically reducing the system of integral equations to a complex polynomial combined with root finding is found to yield accurate results at much lower computational cost. A related method for counting roots of the dispersion relation inside a contour without locating those roots is also demonstrated. Applications of these methods show they can reproduce and exceed the accuracy of previous results in the literature, and new benchmark results are provided. Implementations of the methods described are made available in an accompanying python package psitools.

Keyword:

hydrodynamics

instabilities

methods: numerical

planets and satellites: formation