A numerical scheme for the early steps of nucleation-aggregation models

In the formation of large clusters out of small particles, the initializing step is called the nucleation, and consists in the spontaneous reaction of agents which aggregate into small and stable polymers called nuclei. After this early step, the polymers are involved in a number of reactions such as polymerization, fragmentation and coalescence. Since there may be several orders of magnitude between the size of a particle and the size of an aggregate, building efficient numerical schemes to capture accurately the kinetics of the reaction is a delicate step of key importance. In this article, we propose a conservative scheme, based on finite volume methods on an adaptive grid, which is capable of simulating well the early steps of the reaction as well as the later chain reactions.