Approximating the Becker-Doring cluster equations

Dugald B. Duncan, Ali R. Soheili

Research output: Contribution to journalArticlepeer-review

24 Citations (Scopus)


The Becker-Doring equations model the dynamics of coagulation and fragmentation of clusters of identical particles. The model is an infinite system of ordinary differential equations (ODEs) which specify the rates of change of the concentrations of r-particle clusters. For numerical computation the system is truncated at clusters of a finite size, but this might have to be prohibitively large to capture the metastable behaviour of the system. In this work we derive and investigate the properties of approximations of the Becker-Doring equations which aim to capture the metastable behaviour of the problem with a much reduced system of equations. The schemes are based on a piecewise constant flux approximation, a Galerkin method using a discrete inner product, a discretization of a PDE that in turn approximates the Becker-Doring equations and a lumped coefficients model. We establish a posteriori error estimates for three of these schemes and report on the results of a set of numerical experiments involving metastable behaviour in the solution. Three of the schemes give accurate results on a reduced set of equations, and the PDE scheme appears to be more efficient than the others for this problem.

Original languageEnglish
Pages (from-to)1-29
Number of pages29
JournalApplied Numerical Mathematics
Issue number1-2
Publication statusPublished - Apr 2001


Dive into the research topics of 'Approximating the Becker-Doring cluster equations'. Together they form a unique fingerprint.

Cite this