On the Courant-Friedrichs-Lewy condition for numerical solvers of the coagulation equation

Abstract

Evolving the size distribution of solid aggregates challenges simulations of young stellar objects. Among other difficulties, generic formulae for stability conditions of explicit solvers provide severe constrains when integrating the coagulation equation for astrophysical objects. Recent numerical experiments have recently reported that these generic conditions may be much too stringent. By analysing the coagulation equation in the Laplace space, we explain why this is indeed the case and provide a novel stability condition which avoids time over-sampling.

6 pages, 4 figures, Accepted for publication in MNRAS

References in corpus (8)

• Coagulation, fragmentation and radial motion of solid particles in protoplanetary disks
• Gas- and dust evolution in protoplanetary disks
• Forming Planetesimals in Solar and Extrasolar Nebulae
• Including Dust Coagulation in Hydrodynamic Models of Protoplanetary Disks: Dust Evolution in the Vicinity of a Jupiter-mass Planet
• Grand challenges in protoplanetary disc modelling
• Modeling dust growth in protoplanetary disks: The breakthrough case
• Grain growth for astrophysics with Discontinuous Galerkin schemes
• High-order WENO scheme for Polymerization-type equations