Draft:Super Droplet Method algorithm

Submission declined on 20 May 2026 by Ldm1954 (talk). This draft is not written from a neutral point of view. Wikipedia articles must be written neutrally in a formal, impersonal, and dispassionate way. They should not read like a blog post, advertisement, or fan page. Rewrite the draft to remove: promotional language: see Words to watch; personal commentary: opinions or direct addresses to the reader; informal language. Instead, only summarize in your own words a range of independent, reliable, published sources that discuss the subject. If you would like to continue working on the submission, click on the "Edit" tab at the top of the window. If you have not resolved the issues listed above, your draft will be declined again and potentially deleted. If you need extra help, please ask us a question at the AfC Help Desk or get live help from experienced editors. Please do not remove reviewer comments or this notice until the submission is accepted. Where to get help If you need help editing or submitting your draft, please ask us a question at the AfC Help Desk or get live help from experienced editors. These venues are only for help with editing and the submission process, not to get reviews. If you need feedback on your draft, or if the review is taking a lot of time, you can try asking for help on the talk page of a relevant WikiProject. Some WikiProjects are more active than others so a speedy reply is not guaranteed. How to improve a draft Wikipedia:Contributing to Wikipedia – a basic overview on how to edit Wikipedia. Help:Wikitext – how to use the markup Help:Referencing for beginners – how to include references Wikipedia:Article development – how to develop your article Wikipedia:Writing better articles – how to improve your article Wikipedia:Verifiability – make sure your article includes reliable third-party sources You can also browse Wikipedia:Featured articles and Wikipedia:Good articles to find examples of Wikipedia's best writing on topics similar to your proposed article. Improving your odds of a speedy review To improve your odds of a faster review, tag your draft with relevant WikiProject tags using the button below. This will let reviewers know a new draft has been submitted in their area of interest. For instance, if you wrote about a female astronomer, you would want to add the Biography, Astronomy, and Women scientists tags. Add tags to your draft Editor resources Easy tools: Citation bot (help) | Advanced: Fix bare URLs Declined by Ldm1954 3 months ago. Last edited by Citation bot 34 days ago. Reviewer: Inform author. Resubmit Please note that if the issues are not fixed, the draft will be declined again.

Submission declined on 26 September 2025 by Ldm1954 (talk). This draft is not written from a neutral point of view. Wikipedia articles must be written neutrally in a formal, impersonal, and dispassionate way. They should not read like a blog post, advertisement, or fan page. Rewrite the draft to remove: promotional language: see Words to watch; personal commentary: opinions or direct addresses to the reader; informal language. Instead, only summarize in your own words a range of independent, reliable, published sources that discuss the subject. Declined by Ldm1954 10 months ago.

• Comment: Making a few minor changes and resubmitting does not solve the issues. If this is done again I suggest rejection. Ldm1954 (talk) 12:02, 20 May 2026 (UTC)

• Comment: Currently too much of this is a how-to/textbook, and we do not have those, see WP:NOTTEXTBOOK. Look at some other pages such as Broyden's method, QM/MM, Anderson acceleration and others. I thnk the whole #SDM Monte-Carlo Algorithm section should be heavily trimmed, perhaps removed. Ldm1954 (talk) 13:19, 26 September 2025 (UTC)

In mathematical modeling of aerosols, clouds and precipitation, Super Droplet Method (SDM) is a Monte-Carlo algorithm for representing collisions and coalescence of particles in atmospheric fluid dynamics simulations. The algorithm and its name was introduced in a 2007 arXiv e-print by Shin-ichiro Shima et al...^[1] (preceded by a 2006 patent application^[2] and followed by 2008 RIMS Kôkyûroku^[3] and 2009 QJRMS^[4] papers).

SDM scheme is a probabilistic alternative to the deterministic model of the process embodied in the Smoluchowski coagulation equations. Among the key characteristics of SDM is that it is not subject to the "curse of dimensionality" that hampers application of other methods when multiple particle attributes need to be resolved in a simulation^[7]. The algorithm is embarrassingly parallel, has linear time complexity, constant state vector size (number conservation of simulated particles during collisions) and zero numerical diffusion.

Particle-based approaches, including SDM, for simulating the dynamics of sizes of particles in clouds are considered as one of three main cloud microphysics modelling approaches, along with moment-resolved and bin-resolved models^[8]. Unlike bin- or moment-resolved methods which constitute Eulerian formulation, particle-resolved models are based on Lagrangian formulation in terms of dynamics in particle attribute space.

The term "super-droplet" approach/method has been used either in reference to the particular Monte-Carlo algorithm (even if not used to model clouds^[9]), or more broadly in reference to the particle-based approach for modeling of atmospheric clouds (even if neglecting coalescence processes^[10]). Applications of particle-based methods in atmospheric modelling^[11][12], including for cloud microphysics modelling^[13], and including Monte-Carlo techniques and super-particle nomenclature^[14], predate SDM. Analogous Monte-Carlo particle-based methods (or particle swarm methods) have also been used for modelling accretion in astrophysical context^[15]

The canonical version of SDM follows the "all or nothing"^[16] algorithm presented by Shima et al.^[1]. The SDM collision algorithm is a Gillespie algorithm (see also Kinetic Monte Carlo and Dynamic Monte Carlo methods) in the sense that the expected number of droplet collisions which is represented by superdroplets over a given time in a certain volume matches the expected number of real droplet collisions, although the variance can be substantially higher.

The algorithm formulation presented herein follows the notation of Shima et al.^[1], but uses pseudocode rather than mathematical notation. The code snippets are valid Python (used to generate the animation above) and refer to NumPy and SciPy components.

SDM models evolution of a particulate system composed of ${\displaystyle n_{s}}$ super droplets, with ${\displaystyle i}$-th super droplet representing a multiplicity ${\displaystyle \xi _{i}}$ of particles. Each super particle carries a set of attributes, among which there are extensive attributes, such as mass ${\displaystyle m_{i}}$, as well as auxiliary attributes which do not change upon aggregation, such as position in space.

The considered particle-laden volume is split into cells. Cell volumes should be small enough to consider them well mixed, and hence to assume that any particle in the cell can collide with any other with a probability dependent only on particle extensive attributes. In subsequent formulation of the algorithm, a single cell of volume ${\displaystyle \Delta v}$ is considered.

SDM allows for arbitrary initial sampling of the particle attribute distribution. Among the commonly used methods, there is inverse transform sampling which yields uniform initial multiplicities across the population. Such constant-multiplicity sampling implies evaluation of the quantile function at random locations chosen uniformly from the 0...1 range as in the listing below:

import numpy as np
from scipy import stats

def sample(*,
    rng: np.random.Generator,
    dist: stats.rv_continuous,
    norm: float,
    n_s: int,
):
    state = np.empty(
        shape=n_s,
        dtype=[('ξ', np.int64), ('m', np.float64)]
    )
    state['ξ'] =  round(norm / n_s)
    state['m'] =  dist.ppf(rng.uniform(0, 1, n_s))
    return state

Several other sampling methods were discussed and applied in studies employing SDM^[17][18][19]

Simulation using SDM involves time-stepping loop that advances the system state by ${\displaystyle \Delta t}$ in each step. SDM is a Monte Carlo algorithm and each step employs random number generation. As exemplified in the basic serial implementation of an SDM step given in the listing below, each involves:

• permuting the super particles to select a random set of n_pair non-overlapping pairs among all n_s super particles in the considered volume ${\displaystyle \Delta v}$;
• shuffling n_pair random numbers from a uniform distribution (stored into array φ)
• computing the p_ratio of all possible particle unordered pairs to the number of considered candidate pairs n_pair;
• enumerating across all candidate pairs with α numbering the pair, and j,k numbering the super droplets within a pair;
• evaluating the probability p_α of collision within the super-droplet pair α as a product of the coagulation kernel, the timestep-to-cell-volume ratio Δt/Δv, the p_ratio and ξ_j where j is chosen to point to the super droplet of higher multiplicity in the pair;
• evaluating the γ factor, which for p_α<1 is a Boolean flag resulting from comparison of the probability of coagulation with the random number from the 0...1 range; while in the case of p_α≥1, γ is augmented by the number of "certain" coagulations to represent repeated collision events within a single step (note that depending on the value of the ξ_j/ξ_k ratio, the number of repeated events may need be truncated incurring a collision-event deficit);
• updating particle attributes ξ_j, ξ_k, m_j and m_k.

from collections import abc

def sdm(*,
    rng: np.random.Generator,
    ξ: abc.MutableSequence[int],
    m: abc.MutableSequence[float],
    kern: abc.Callable[[float, float], float],
    Δt: float,
    Δv: float,
):
    """ SDM step assuming non-zero multiplicities """
    n_s = len(ξ)
    n_pair = n_s // 2
    
    pairs = rng.permutation(n_s)[: 2 * n_pair]
    φ = rng.uniform(0, 1, n_pair)

    p_ratio = n_s * (n_s - 1) / 2 / n_pair
    for α, (j, k) in enumerate(pairs.reshape(-1, 2)):
        p_jk = kern(m[j], m[k]) * Δt / Δv
        if ξ[j] < ξ[k]:
            j, k = k, j
        p_α = ξ[j] * p_ratio * p_jk
        γ = p_α // 1 + ((p_α - p_α // 1) > φ[α])
        if γ != 0:
            γ = min(γ, (ξ[j] / ξ[k]) // 1)
            if ξ[j] - γ * ξ[k] > 0:
                ξ[j] -= γ * ξ[k]
                m[k] += γ * m[j]
            else:
                ξ[j] = ξ[k] // 2
                ξ[k] -= ξ[j]
                m[k] += γ * m[j]
                m[j] = m[k]

The above algorithm can yield the following collision scenarios (see animations above, symbol m denotes an arbitrary mass unit, symbol M denotes the mass attribute of the super-droplet of respective color):

- A: single-collision event (γ=1 )

- B: multiple-collision event (γ > 1 and ξ[j] ≠ γξ[k])

- C: equal-ratio collision event (γ ≥ 1 and ξ[j] = γξ[k]), note: an equal-ratio collision event may also occur within scenario A

The SDM coagulation algorithm does not feature data dependencies across candidate pairs, hence is particularly well suited for SIMT concurrent operation. Implementations of SDM leverage both CPU and GPU multi-threading.^[20]

The SDM algorithm has been applied in diverse modelling studies, covering topics such as:

• mixed-phase clouds^[19]
• representation of collisional breakup^[21]
• aqueous chemical reactions involving soluble aerosol material^[22]
• in-cloud electro-coalescence^[23]
• role of clouds in nuclear fallout dynamics^[24]
• geo-engineering (cloud brightening)^[25] (the reference uses particle-based microphysics representation, but collisions are represented in a different way than SDM)
• simulations of processes occurring in cloud-chamber laboratories^[26][27]
• interaction of cosmic rays with clouds^[28]
• dynamics of fogs^[29]

(links point to Git repositories)

• reusable SDM libraries and packages:
  • libcloudph++ (C++/CUDA)
  • PySDM (Python/Numba/CUDA)
  • CLEO (C++/Kokkos), also integrated into ICON model.
  • Droplets.jl (Julia)
• SDM implementations integrated in other software:
  • SCALE-SDM (Fortran)
  • Pencil Code (Fortran)
  • PALM LES (Fortran)
  • LCM1D (Python)
  • superdroplet (Cython/Numba/C++11/Fortran 2008/Julia)
  • NTLP (FORTRAN)
  • LacmoPy (Python/Numba)
  • McSnow (FORTRAN)
  • Particula (Python)