Stochastic modelling of silicon nanoparticle synthesis
Repository URI
Repository DOI
Change log
Authors
Abstract
This thesis presents new methods to study the aerosol synthesis of nano-particles and a new model to simulate the formation of silicon nanoparticles. Population balance modelling is used to model nanoparticle synthesis and a stochastic numerical method is used to solve the governing equations. The population balance models are coupled to chemical kinetic models and offer insight into the fundamental physiochemical processes leading to particle formation.
The first method developed in this work is a new mathematical expression for calculating the rate of Brownian coagulation with stochastic weighted algorithms (SWAs). The new expression permits the solution of the population balance equations with SWAs using a computationally-efficient technique of majorant rates and fictitious jumps. Convergence properties and efficiency of the expression are evaluated using a detailed silica particle model. A sequential-modular algorithm is subsequently presented which solves networks of perfectly stirred reactors with a population balance model using the stochastic method. The algorithm is tested in some simple network configurations, which are used to identify methods through which error in the stochastic solution may be reduced. It is observed that SWAs are useful in preventing accumulation of error in reactor networks. A new model for silicon nanoparticle synthesis is developed. The model includes gas-phase reactions describing silane decomposition, and a detailed multivariate particle model which tracks particle structure and composition. Systematic parameter estimation is used to fit the model to experimental cases. Results indicated that the key challenge in modelling silicon systems is obtaining a correct description of the particle nucleation process. Finally, the silicon model is used in conjunction with the reactor network algorithm to simulate the start-up of a plug-flow reactor. The power of stochastic methods in resolving characteristics of a particle ensemble is highlighted by investigating the number, size, degree of sintering and polydispersity along the length of the reactor.