Agglomeration and breakage of nanoparticles in stirred media mills-a comparison of different methods and models

The increasing industrial demand for nanoparticles challenges the application of stirred media mills to grind in the sub-micron size range. It was shown recently [Mende et al., 2003. Mechanical production and stabilization of submicron particles in stirred media mills. Powder Technology 132, 64-73] that the grinding behavior of particles in the sub-micron size range in stirred media mills and the minimum achievable particle size is strongly influenced by the suspension stability and thus the agglomeration behavior of the suspension. Therefore, an appropriate modeling of the process must include a superposition of the two opposing processes in the mill i.e., breakage and agglomeration which can be done by means of population balance models. Modeling must now include the influence of colloidal surface forces and hydrodynamic forces on particle aggregation and breakup. The superposition of the population balance models for agglomeration and grinding with the appropriate kernels leads to a system of partial differential equations, which can be solved in various ways numerically. Here a modified h-p Galerkin algorithm which is implemented in the commercially available software package PARSIVAL developed by CiT (CiT GmbH, Rastede, Germany) and the moment methodology according to [Diemer and Olsen, 2002a. A moment methodology for coagulation and breakage problems: Part I-analytical solution of the steady-state population balance. Chemical Engineering Science 57 (12), 2193-2209; Diemer and Olsen, 2002b. A moment methodology for coagulation and breakage problems: Part II-moment models and distribution reconstruction. Chemical Engineering Science 57 (12), 2211-2288] are used and compared to explicit data on alumina. This includes a comparison of the derived particle size distributions, moments and its accuracy depending on the starting particle size distribution and the used agglomeration and breakage kernels. Finally, the computational effort of both methods in comparison to the prior mentioned parameters is evaluated in terms of practical application.