The phenomenon of coagulation occurs in a wide range of applications, e.g., in physics (aggregation of particles, growth of gas bubbles), meteorology (merging of drops in atmospheric clouds, aerosol transport), chemistry (reacting polymers, soot formation) and astrophysics (formation of stars and planets).
A common tool for the mathematical treatment of the subject are nonlinear partial integrodifferential equations. In particular, the time evolution of the average concentration of objects of a given size is described by Smoluchowski's coagulation equation and its spatially inhomogeneous extensions. The close connections between the deterministic equations and stochastic interacting particle systems provide the basis both for analytical studies and for numerical approximations.
Contribution of the Insitute
Many processes in chemical engineering have been successfully modelled by stochastic particle systems in recent years. For example, the formation of soot particles in combustion systems is a highly complex process, which needs rather sophisticated modelling of the resulting soot particle populations. The stochastic approach allows for very detailed models of individual particles to be included with little additional impact on computation times.
In collaboration with the Department of Chemical Engineering and Biotechnology (University of Cambridge, UK), a mathematical model for such processes has been developed. In this model, each particle is described by its constituent primary particles and the connectivity between these primaries (see Figure). Each primary in turn has internal variables that describe its chemical composition.
Stochastic weighted algorithms were applied to this multidimensional particle model. The weighting mechanisms are designed in such a way that various physical processes changing individual particles (such as growth, or other surface reactions) can be treated by the algorithms. A spatially resolved stochastic particle model for coagulationadvection problems in bounded domains with inception at the boundary was studied. Convergence of the particle system to a deterministic limit was proved. It was demonstrated that stochastic weighted algorithms can be successfully used with realistic coagulation kernels and high dimensional particle models.
Stochastic particle methods depend on several discretization parameters, due to the finite number of particles, the decoupling of transport and interaction (time step), and the spatial delocalization of the interaction (cell size). The dependence on the cell size of the steady state solution obtained in the infiniteparticlenumber limit was studied.
Publications
Articles in Refereed Journals

A. Mielke, R.I.A. Patterson, M.A. Peletier, D.R.M. Renger, Nonequilibrium thermodynamical principles for chemical reactions with massaction kinetics, SIAM Journal on Applied Mathematics, 77 (2017), pp. 15621585, DOI 10.20347/WIAS.PREPRINT.2165 .
Abstract
We study stochastic interacting particle systems that model chemical reaction networks on the micro scale, converging to the macroscopic Reaction Rate Equation. One abstraction level higher, we study the ensemble of such particle systems, converging to the corresponding Liouville transport equation. For both systems, we calculate the corresponding large deviations and show that under the condition of detailed balance, the large deviations induce a nonlinear relation between thermodynamic fluxes and free energy driving force. 
R.I.A. Patterson, Properties of the solutions of delocalised coagulation and inception problems with outflow boundaries, Journal of Evolution Equations, 16 (2016), pp. 261291.
Abstract
Well posedness is established for a family of equations modelling particle populations undergoing delocalised coagulation, advection, inflow and outflow in a externally specified velocity field. Very general particle types are allowed while the spatial domain is a bounded region of $d$dimensional space for which every point lies on exactly one streamline associated with the velocity field. The problem is formulated as a semilinear ODE in the Banach space of bounded measures on particle position and type space. A local Lipschitz property is established in total variation norm for the propagators (generalised semigroups) associated with the problem and used to construct a Picard iteration that establishes local existence and global uniqueness for any initial condition. The unique weak solution is shown further to be a differentiable or at least bounded variation strong solution under smoothness assumptions on the parameters of the coagulation interaction. In the case of one spatial dimension strong differentiability is established even for coagulation parameters with a particular bounded variation structure in space. This one dimensional extension establishes the convergence of the simulation processes studied in [Patterson, Stoch. Anal. Appl. 31, 2013] to a unique and differentiable limit. 
R.I.A. Patterson, W. Wagner, Cell size error in stochastic particle methods for coagulation equations with advection, SIAM Journal on Numerical Analysis, 52 (2014), pp. 424442.
Abstract
The paper studies the approximation error in stochastic particle methods for spatially inhomogeneous population balance equations. The model includes advection, coagulation and inception. Sufficient conditions for second order approximation with respect to the spatial discretization parameter (cell size) are provided. Examples are given, where only first order approximation is observed. 
M. Sander, R.I.A. Patterson, A. Braumann, A. Raj, M. Kraft, Developing the PAHPP soot particle model using process informatics and uncertainty propagation, Proceedings of the Combustion Institute, 33 (2011), pp. 675683.
Abstract
n this work we present the new PAHPP soot model and use a data collaboration approach to determine some of its parameters. The model describes the formation, growth and oxidation of soot in laminar premixed flames. Soot particles are modelled as aggregates containing primary particles, which are built from polycyclic aromatic hydrocarbons (PAHs), the main building blocks of a primary particle (PP). The connectivity of the primary particles is stored and used to determine the rounding of the soot particles due to surface growth and condensation processes. Two neighbouring primary particles are replaced by one if the coalescence level between the two primary particles reaches a threshold. The model contains, like most of the other models, free parameters that are unknown a priori. The experimental premixed flame data from Zhao et al. [B. Zhao, Z. Yang, Z. Li, M.V. Johnston, H. Wang, Proc. Combust. Inst. 30 (2) (2005) 1441?1448] have been used to estimate the smoothing factor of soot particles, the growth factor of PAHs within particles and the soot density using a low discrepancy series method with a subsequent response surface optimisation. The optimised particle size distributions show good agreement with the experimental ones. The importance of a standardised data mining system in order to optimise models is underlined. 
R.I.A. Patterson, W. Wagner, M. Kraft, Stochastic weighted particle methods for population balance equations, Journal of Computational Physics, 230 (2011), pp. 74567472.
Abstract
A class of stochastic algorithms for the numerical treatment of population balance equations is introduced. The algorithms are based on systems of weighted particles, in which coagulation events are modelled by a weight transfer that keeps the number of computational particles constant. The weighting mechanisms are designed in such a way that physical processes changing individual particles (such as growth, or other surface reactions) can be conveniently treated by the algorithms. Numerical experiments are performed for complex laminar premixed flame systems. Two members of the class of stochastic weighted particle methods are compared to each other and to a direct simulation algorithm. One weighted algorithm is shown to be consistently better than the other with respect to the statistical noise generated. Finally, run times to achieve fixed error tolerances for a real flame system are measured and the better weighted algorithm is found to be up to three times faster than the direct simulation algorithm. 
M. Celnik, R.I.A. Patterson, M. Kraft, W. Wagner, A predictorcorrector algorithm for the coupling of stiff ODEs to a particle population balance, Journal of Computational Physics, 228 (2009), pp. 27582769.

M. Celnik, R.I.A. Patterson, M. Kraft, W. Wagner, Coupling a stochastic soot population balance to gasphase chemistry using operator splitting, , 148 (2007), pp. 158176.

N. Morgan, C. Wells, M. Goodson, M. Kraft, W. Wagner, A new numerical approach for the simulation of the growth of inorganic nanoparticles, Journal of Computational Physics, 211 (2006), pp. 638658.

R.I.A. Patterson, J. Singh, M. Balthasar, M. Kraft, W. Wagner, Extending stochastic soot simulation to higher pressures, Combustion and Flame, 145 (2006), pp. 638642.

C. Wells, N. Morgan, M. Kraft, W. Wagner, A new method for calculating the diameters of partiallysintered nanoparticles and its effect on simulated particle properties, , 61 (2006), pp. 158166.

M. Kraft, S. Mosbach, W. Wagner, Using a web module to teach stochastic modelling, Chemical Engineering Education, 39 (2005), pp. 244248.

N. Morgan, W. Wells, M. Kraft, W. Wagner, Modelling nanoparticle dynamics: Coagulation, sintering, particle inception and surface growth, Combustion Theory and Modelling, 9 (2005), pp. 449461.

J. Singh, M. Balthasar, M. Kraft, W. Wagner, Stochastic modeling of soot particle size and age distributions in laminar premixed flames, Proceedings of the Combustion Institute, 30 (2005), pp. 14571465.

C. Lecot, W. Wagner, A quasiMonte Carlo scheme for Smoluchowski's coagulation equation, Mathematics of Computation, 73 (2004), pp. 19531966.

M. Kraft, W. Wagner, An improved stochastic algorithm for temperaturedependent homogeneous gas phase reactions, Journal of Computational Physics, 185 (2003), pp. 139157.

M. Kraft, W. Wagner, Numerical study of a stochastic particle method for homogeneous gas phase reactions, Computers & Mathematics with Applications. An International Journal, 45 (2003), pp. 329349.

W. Wagner, Stochastic, analytic and numerical aspects of coagulation processes, Mathematics and Computers in Simulation, 62 (2003), pp. 265275.
Contributions to Collected Editions

D. Knees, R. Kornhuber, Ch. Kraus, A. Mielke, J. Sprekels, C3  Phase transformation and separation in solids, in: MATHEON  Mathematics for Key Technologies, M. Grötschel, D. Hömberg, J. Sprekels, V. Mehrmann ET AL., eds., 1 of EMS Series in Industrial and Applied Mathematics, European Mathematical Society Publishing House, Zurich, 2014, pp. 189203.
Preprints, Reports, Technical Reports

R.I.A. Patterson, D.R.M. Renger, Dynamical large deviations of countable reaction networks under a weak reversibility condition, Preprint no. 2273, WIAS, Berlin, 2016, DOI 10.20347/WIAS.PREPRINT.2273 .
Abstract, PDF (343 kByte)
A dynamic large deviations principle for a countable reaction network including coagulationfragmentation models is proved. The rate function is represented as the infimal cost of the reaction fluxes and a minimiser for this variational problem is shown to exist. A weak reversibility condition is used to control the boundary behaviour and to guarantee a representation for the optimal fluxes via a Lagrange multiplier that can be used to construct the changes of measure used in standard tilting arguments. Reflecting the pure jump nature of the approximating processes, their paths are treated as elements of a BV function space. 
A. Mielke, R.I.A. Patterson, M.A. Peletier, D.R.M. Renger, Nonequilibrium thermodynamical principles for chemical reactions with massaction kinetics, Preprint no. 2165, WIAS, Berlin, 2015, DOI 10.20347/WIAS.PREPRINT.2165 .
Abstract, PDF (363 kByte)
We study stochastic interacting particle systems that model chemical reaction networks on the micro scale, converging to the macroscopic Reaction Rate Equation. One abstraction level higher, we study the ensemble of such particle systems, converging to the corresponding Liouville transport equation. For both systems, we calculate the corresponding large deviations and show that under the condition of detailed balance, the large deviations induce a nonlinear relation between thermodynamic fluxes and free energy driving force.
Talks, Poster

R.I.A. Patterson, Monte Carlo simulation of soot, King Abdullah University of Science and Technology (KAUST), Clean Combustion Research Center, Thuwal, Saudi Arabia, January 11, 2016.

R.I.A. Patterson, Population balance simulation, University of Cambridge, Department for Chemical Engineering and Biotechnology, UK, May 5, 2016.

R.I.A. Patterson, Simulations of flame generated particles, Advances in Uncertainty Quantification Methods, Algorithms and Applications (UQAW 2016), January 5  10, 2016, King Abdullah University of Science and Technology (KAUST), Thuwal, Saudi Arabia, January 5, 2016.

D.R.M. Renger, Large deviations for reacting particle systems: The empirical and ensemble processes, Workshop ``Interplay of Analysis and Probability in Applied Mathematics'', July 26  August 1, 2015, Mathematisches Forschungsinstitut Oberwolfach, Oberwolfach, July 30, 2015.

D.R.M. Renger, The empirical process of reacting particles: Large deviations and thermodynamic principles, Minisymposium ``Real World Phenomena Explained by Microscopic Particle Models'' of the 8th International Congress on Industrial and Applied Mathematics (ICIAM 2015), August 8  22, 2015, International Council for Industrial and Applied Mathematics, Beijing, China, August 10, 2015.

R.I.A. Patterson, Approximation errors for Smoluchowski simulations, 10 th IMACS Seminar on Monte Carlo Methods, July 6  10, 2015, Johannes Kepler Universität Linz, Austria, July 7, 2015.

R.I.A. Patterson, Particle systems, kinetic equations and their simulation, 8th International Congress on Industrial and Applied Mathematics, ICIAM 2015, August 8  15, 2015, CNCC  China National Convention Center, Beijing, China.

R.I.A. Patterson, Uniqueness and regularity for coagulationadvection problems, Workshop on Theory and Numerics of Kinetic Equations, June 1  4, 2015, Universität Saarbrücken, June 2, 2015.

R.I.A. Patterson, Convergence of coagulationadvection simulations, Workshop on Particle Transport with Emphasis on Stochastics, November 6  7, 2014, Aarhus University, Denmark, November 6, 2014.

R.I.A. Patterson, Monte Carlo simulation of nanoparticle formation, University of Technology Eindhoven, Institute for Complex Molecular Systems, Netherlands, September 5, 2013.

R.I.A. Patterson, Stochastic methods for particle coagulation problems in flows with boundaries, 5th Workshop ``Theory and Numerics of Kinetic Equations'', May 13  15, 2013, Universität des Saarlandes, Saarbrücken, May 14, 2013.

W. Wagner, Stochastic particle methods for population balance equations, 5th Workshop ``Theory and Numerics of Kinetic Equations'', May 13  15, 2013, Universität des Saarlandes, Saarbrücken, May 13, 2013.

W. Wagner, Stochastic particle methods for coagulation problems, 28th International Symposium on Rarefied Gas Dynamics, July 9  13, 2012, University of Zaragoza, Spain, July 9, 2012.

W. Wagner, Fragmentation equations and stochastic models, Workshop ``Asymptotic Methods for Dissipative Particle Systems'', May 18  22, 2009, University of California, Institute for Pure and Applied Mathematics, Los Angeles, USA, May 18, 2009.

W. Wagner, Fragmentation models, Workshop on Theory and Applications of Classical and Quantum Kinetic Theory, June 21  26, 2009, Banff International Research Station for Mathematical Innovation and Discovery, Canada, June 23, 2009.

W. Wagner, Introduction to Markov jump processes, Scuola Superiore di Catania, Laboratory for Complex Systems, Italy, October 27, 2009.

W. Wagner, Random and deterministic fragmentation models, Seventh IMACS Seminar on Monte Carlo Methods (MCM2009), September 6  11, 2009, Université Libre de Bruxelles, Brussels, Belgium, September 7, 2009.

W. Wagner, Explosion properties of random fragmentation models, Workshop ``Coagulation et Fragmentation Stochastiques'', April 15  18, 2008, Université Paris VI, Laboratoire de Probabilités, France, April 16, 2008.

W. Wagner, Explosion properties of random fragmentation models, Università di Catania, Dipartimento di Matematica e Informatica, Italy, May 8, 2008.

W. Wagner, Stochastic models for coagulationfragmentation processes, Workshop on CoagulationFragmentation Processes: Theory and Application, July 4  8, 2005, University of Edinburgh, UK, July 5, 2005.

W. Wagner, Modelling nanoparticle dynamics: Coagulation, sintering, particle inception and surface growth, MC$^2$QMC 2004, June 7  10, 2004, JuanlesPins, France, June 7, 2004.

W. Wagner, Numerical simulation of coagulation dynamics using quasirandom numbers, MC$^2$QMC 2004, June 7  10, 2004, JuanlesPins, France, June 7, 2004.

W. Wagner, A quasiMonte Carlo simulation of cluster coagulation, IV IMACS Seminar on Monte Carlo Methods, September 15  19, 2003, Berlin, September 18, 2003.
External Preprints

W.J. Menz, R.I.A. Patterson, W. Wagner, M. Kraft, Application of stochastic weighted algorithms to a multidimensional silica particle model, Preprint no. 120, University of Cambridge, Cambridge Center for Computational Chemical Engineering, 2012.

A. Braumann, M. Kraft, W. Wagner, Numerical study of a stochastic particle algorithm solving a multidimensional population balance model for high shear granulation, Technical report no. 93, c4ePreprint Series, Cambridge Centre for Computational Chemical Engineering, University of Cambridge, Department of Chemical Engineering, 2010.
Abstract
This paper is concerned with computational aspects of a multidimensional population balance model of a wet granulation process. Wet granulation is a manufacturing method to form composite particles, granules, from small particles and binders. A detailed numerical study of a stochastic particle algorithm for the solution of a fivedimensional population balance model for wet granulation is presented. Each particle consists of two types of solids (containing pores) and of external and internal liquid (located in the pores). Several transformations of particles are considered, including coalescence, compaction and breakage. A convergence study is performed with respect to the parameter that determines the number of numerical particles. Averaged properties of the system are computed. In addition, the ensemble is subdivided into practically relevant size classes and analysed with respect to the amount of mass and the particle porosity in each class. These results illustrate the importance of the multidimensional approach. Finally, the kinetic equation corresponding to the stochastic model is discussed. 
M.S. Celnik, R.I.A. Patterson, M. Kraft, W. Wagner, A predictorcorrector algorithm for the coupling of stiff ODEs to a particle population balance, Preprint no. 58, Cambridge Center for Computational Chemical Engineering, 2008.
Abstract
In this paper a novel predictorcorrector algorithm is presented for the solution of coupled gasphaseparticulate systems. The emphasis of this work is the study of soot formation, but the concepts can be applied to other systems. This algorithm couples a stiff ODE solver to a Monte Carlo population balance solver. Such coupling has been achieved previously for similar systems using a Strang operator splitting algorithm, however, that algorithm demonstrated several numerical issues which resulted in a high computational cost to acquire adequate precision. In particular a sourcesink instability was identified whereby a largemagnitude source term present in the ODE system was competing with a similarly sized sink term in the population balance. This instability required that the splittingstep size was very small in order to keep numerical error sufficiently low. A predictorcorrector algorithm has been formulated to negate this instability. An additional efficiency is gained with this algorithm as a principal computational cost of the Strang splitting algorithm is removed: the requirement to reinitialise the ODE solver every splittingstep. The numerical convergence of the new algorithm is demonstrated, and its efficiency is compared to that of the Strang splitting algorithm. Substantial computation time savings are demonstrated, which allow a fixed error in three studied system functionals to be achieved with an orderofmagnitude reduction in computation time.