O(N^2) fragmentation algorithm
arXiv:1911.02043 · doi:10.3847/1538-4365/ab7b71
Abstract
Collisional fragmentation is a ubiquitous phenomenon arising in a variety of astrophysical systems, from asteroid belts to debris and protoplanetary disks. Numerical studies of fragmentation typically rely on discretizing the size distribution of colliding objects into a large number N of bins in mass space, usually logarithmically spaced. A standard approach for redistributing the debris produced in collisions into the corresponding mass bins results in O(N^3) calculation, which leads to significant computational overhead when N is large. Here we formulate a more efficient explicit O(N^2) fragmentation algorithm, which works when the size spectrum of fragments produced in an individual collision has a self-similar shape with only a single characteristic mass scale (which can have arbitrary dependence on the energy and masses of colliding objects). Fragment size spectra used in existing fragmentation codes typically possess this property. We also show that our O(N^2) approach can be easily extended to work with non-self-similar fragment size distributions, for which we provide a worked example. This algorithm offers a substantial speedup of fragmentation calculations for large N> 100, even over the implicit methods, making it an attractive tool for studying collisionally evolving systems.
arXiv comment9 pages, 5 figures, submitted to AAS Journals
References in corpus (18)
- Coagulation, fragmentation and radial motion of solid particles in protoplanetary disks
- Gas- and dust evolution in protoplanetary disks
- The outcome of protoplanetary dust growth: pebbles, boulders, or planetesimals? I. Mapping the zoo of laboratory collision experiments
- Planetesimal formation by sweep-up: How the bouncing barrier can be beneficial to growth
- Long-Term Collisional Evolution of Debris Disks
- Accretion in the Early Kuiper Belt II. Fragmentation
- Steady-State Size Distributions for Collisional Populations: Analytical Solution with Size-Dependent Strength
- Numerical Modeling of the Coagulation and Porosity Evolution of Dust Aggregates
- The Size Distribution of Kuiper Belt Objects
- Prospects for Detection of Catastrophic Collisions in Debris Disks
- From dust to planetesimals: an improved model for collisional growth in protoplanetary disks
- Including Dust Coagulation in Hydrodynamic Models of Protoplanetary Disks: Dust Evolution in the Vicinity of a Jupiter-mass Planet
- Redistribution of CO at the Location of the CO Ice Line in evolving Gas and Dust Disks
- Collisional and Thermal Emission Models of Debris Disks: Towards Planetesimal Population Properties
- Inhibited Coagulation of Micron-size Dust Due to the Electrostatic Barrier
- Effects of Ringed Structures and Dust Size Growth on Millimeter Observations of Protoplanetary Disks
- Breakthrough revisited: investigating the requirements for growth of dust beyond the bouncing barrier
- Non-Power Law Behavior in Fragmentation Cascades