A Robust Numerical Solution of the Stochastic Collection–Breakup Equation for Warm Rain

The focus of this paper is on the numerical solution of the stochastic collection equation–stochastic breakup equation (SCE–SBE) describing the evolution of raindrop spectra in warm rain. The drop size distribution (DSD) is discretized using the fixed-pivot scheme proposed by Kumar and Ramkrishna, and new discrete equations for solving collision breakup are presented. The model is evaluated using established coalescence and breakup parameterizations (kernels) available in the literature, and in that regard this paper provides a substantial review of the relevant science. The challenges posed by the need to achieve stable and accurate numerical solutions of the SCE–SBE are examined in detail. In particular, this paper focuses on the impact of varying the shape of the initial DSD on the equilibrium solution of the SCE–SBE for a wide range of rain rates and breakup kernels. The results show that, although there is no dependence of the equilibrium DSD on initial conditions for the same rain rate and breakup kernel, there is large variation in the time that it takes to reach steady state. This result suggests that, in coupled simulations of in-cloud motions and microphysics and for short time scales (30 min) for which transient conditions prevail, the equilibrium DSD may not be attainable except for very heavy rainfall. Furthermore, simulations for the same initial conditions show a strong dependence of the dynamic evolution of the DSD on the breakup parameterization. The implication of this result is that, before the debate on the uniqueness of the shape of the equilibrium DSD can be settled, there is critical need for fundamental research including laboratory experiments to improve understanding of collisional mechanisms in DSD evolution. 1