A representative particle approach to coagulation and fragmentation of dust aggregates and fluid droplets

TL;DR

This paper introduces a new Monte-Carlo algorithm for modeling dust particle coagulation and fragmentation, capable of including complex grain properties and suitable for integration into hydrodynamic simulations.

Contribution

The paper presents a novel Monte-Carlo method that efficiently models aggregation and fragmentation of dust particles with evolving properties, improving upon existing techniques.

Findings

Reproduces analytic solutions of simplified kernels

Matches numerical results for Brownian motion

Performs well with moderate swarm numbers

Abstract

Context: There is increasing need for good algorithms for modeling the aggregation and fragmentation of solid particles (dust grains, dust aggregates, boulders) in various astrophysical settings, including protoplanetary disks, planetary- and sub-stellar atmospheres and dense molecular cloud cores. Here we describe a new algorithm that combines advantages of various standard methods into one. Aims: The aim is to develop a method that 1) can solve for aggregation and fragmentation, 2) can easily include the effect and evolution of grain properties such as compactness, composition, etc., and 3) can be built as a coagulation/fragmentation module into a hydrodynamics simulations. Methods: We develop a Monte-Carlo method in which we follow the 'life' of a limited number of representative particles. Each of these particles is associated with a certain fraction of the total dust mass and thereby represents a large number of true particles which all are assumed to have the same properties as their representative particle. Under the assumption that the total number of true particles vastly exceeds the number of representative particles, the chance of a representative particle colliding with another representative particle is negligibly small, and we therefore ignore this possibility. This now makes it possible to employ a statistical approach to the evolution of the representative particles. Results: The method reproduces the known analytic solutions of simplified coagulation kernels, and compares well to numerical results for Brownian motion using other methods. For reasonably well-behaved kernels it produces good results even for moderate number of swarms.