Integration of Particle-Gas Systems with Stiff Mutual Drag Interaction

Chao-Chin Yang,Anders Johansen

DOI: https://doi.org/10.3847/0067-0049/224/2/39

2016-04-15

Abstract:Numerical simulation of numerous mm/cm-sized particles embedded in a gaseous disk has become an important tool in the study of planet formation and in understanding the dust distribution in observed protoplanetary disks. However, the mutual drag force between the gas and the particles can become so stiff, particularly because of small particles and/or strong local solid concentration, that an explicit integration of this system is computationally formidable. In this work, we consider the integration of the mutual drag force in a system of Eulerian gas and Lagrangian solid particles. Despite the entanglement between the gas and the particles under the particle-mesh construct, we are able to devise a numerical algorithm that effectively decomposes the globally coupled system of equations for the mutual drag force and makes it possible to integrate this system on a cell-by-cell basis, which considerably reduces the computational task required. We use an analytical solution for the temporal evolution of each cell to relieve the time-step constraint posed by the mutual drag force as well as to achieve the highest degree of accuracy. To validate our algorithm, we use an extensive suite of benchmarks with known solutions in one, two, and three dimensions, including the linear growth and the nonlinear saturation of the streaming instability. We demonstrate numerical convergence and satisfactory consistency in all cases. Our algorithm can for example be applied to model the evolution of the streaming instability with mm/cm-sized pebbles at high mass loading, which has important consequences for the formation scenarios of planetesimals.

Instrumentation and Methods for Astrophysics,Earth and Planetary Astrophysics,Computational Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the numerical calculation difficulties encountered when simulating the mutual drag force between millimeter - to - centimeter - sized particles and gas in the protoplanetary disk. Specifically, when the particle size is small or the local solid concentration is high, the mutual drag force between the particles and the gas becomes very strong (i.e., "stiff"), which makes the traditional explicit integration method extremely difficult computationally, because the time step must be very small to maintain numerical stability. This small time step leads to very high computational costs and even makes large - scale simulations impractical.

To solve this problem, the author proposes a new numerical algorithm that can effectively decompose the globally coupled equation system, allowing the mutual drag force to be integrated cell - by - cell, thereby significantly reducing the required computational tasks. In addition, by using the analytical solution to predict the velocity changes of gas and particles in each cell, this algorithm not only alleviates the time - step constraints caused by the mutual drag force but also achieves the highest accuracy.

Specifically, the main contributions of the paper include:
1. **Algorithm Design**: Propose a numerical algorithm that can effectively handle the stiff mutual drag force between the Eulerian gas grid and the Lagrangian solid particles.
2. **Application of Analytical Solution**: Use the analytical solution to predict the velocity changes in each cell, thereby avoiding the strict time - step limitations caused by the drag force.
3. **Verification and Testing**: Verify the numerical convergence and consistency of the algorithm through a series of one - dimensional, two - dimensional, and three - dimensional benchmark tests with known solutions.
4. **Application Prospects**: This algorithm can be applied to simulate the evolution of fluid instabilities (such as fluid instabilities) under high mass loading, which is of great significance for understanding the formation process of kilometer - sized planetesimals in the protoplanetary disk.

In summary, this paper aims to provide an efficient and accurate method for simulating the complex interactions between particles and gas in the protoplanetary disk, especially when the drag force is very strong.