A Survey of Numerical Solutions to the Coagulation Equation
Abstract
We present the results of a systematic survey of numerical solutions to the coagulation equation for a rate coefficient of the form Aij (imu jnu + inu jmu) and monodisperse initial conditions. The results confirm that there are three classes of rate coefficients with qualitatively different solutions. For nu ≤ 1 and lambda = mu + nu ≤ 1, the numerical solution evolves in an orderly fashion and tends toward a self-similar solution at large time t. The properties of the numerical solution in the scaling limit agree with the analytic predictions of van Dongen and Ernst. In particular, for the subset with mu > 0 and lambda < 1, we disagree with Krivitsky and find that the scaling function approaches the analytically predicted power-law behavior at small mass, but in a damped oscillatory fashion that was not known previously. For nu ≤ 1 and lambda > 1, the numerical solution tends toward a self-similar solution as t approaches a finite time t0. The mass spectrum nk develops at t0 a power-law tail nk k-tau at large mass that violates mass conservation, and runaway growth/gelation is expected to start at tcrit = t0 in the limit the initial number of particles n0 -> ∞. The exponent tau is in general less than the analytic prediction (lambda + 3)/2, and t0 = K/[(lambda - 1) n0 A11] with K = 1--2 if lambda > 1.1. For nu > 1, the behaviors of the numerical solution are similar to those found in a previous paper by us. They strongly suggest that there are no self-consistent solutions at any time and that runaway growth is instantaneous in the limit n0 -> ∞. They also indicate that the time tcrit for the onset of runaway growth decreases slowly toward zero with increasing n0.
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.