Quantitative analysis of multi-particle interactions during particle breakage: A discrete non-linear population balance framework

Abstract Discrete linear population balance models have been widely used for comminution processes and particle breakage tests that involve a series of breakage events. The fundamental tenet of these models is the independence of particle breakage from the surrounding particle population. On the other hand, ample experimental data have shown that multi-particle interactions among particles of different sizes at the particle ensemble scale lead to a population-dependent breakage probability. In an attempt to analyze and model the multi-particle interactions, we formulate a time- or space-discrete and particle size-discrete non-linear population balance model using a recently introduced non-linear kinetics framework for rate-based comminution processes. The discrete model is developed as a limiting case and presented in equivalent matrix and index notation. The practical applicability, accuracy, and stability aspects of the model are discussed. With an empirical, flexible, non-linear functional accounting for multi-particle interactions explicitly, the proposed model can serve as a framework for quantitatively analyzing multi-particle interactions. Numerical simulations with a simplified discrete model have provided significant insight into the retardation effect of fines on the breakage probability of the coarser particles experimentally observed in dense-phase particulate systems such as particle bed compression tests.

[1]  Colin Webb,et al.  On predicting roller milling performance: Part I: the breakage equation , 2001 .

[2]  Ecevit Bilgili,et al.  Formulation of a non-linear framework for population balance modeling of batch grinding : Beyond first-order kinetics , 2006 .

[3]  Wolfgang Peukert,et al.  Characterisation of Grinding‐Relevant Particle Properties by Inverting a Population Balance Model , 2002 .

[4]  Benjamin Epstein,et al.  Logarithmico-Normal Distribution in Breakage of Solids , 1948 .

[5]  Lawrence F. Shampine,et al.  The MATLAB ODE Suite , 1997, SIAM J. Sci. Comput..

[6]  W. Wildeboer,et al.  Identification of the grinding mechanisms and their origin in a stirred ball mill using population balances , 2009 .

[7]  Peng Wang,et al.  Mathematical modeling of fluid energy milling based on a stochastic approach , 2010 .

[8]  J. Dodds,et al.  Identification of the fragmentation mechanisms in wet-phase fine grinding in a stirred bead mill , 1997 .

[9]  D. Fuerstenau,et al.  Effect of fine particles on the kinetics and energetics of grinding coarse particles , 1991 .

[10]  Ecevit Bilgili,et al.  Nano-milling of pigment agglomerates using a wet stirred media mill: Elucidation of the kinetics and breakage mechanisms , 2006 .

[11]  Rajesh N. Dave,et al.  Identification of the breakage rate and distribution parameters in a non-linear population balance model for batch milling , 2011 .

[12]  Leonard G. Austin,et al.  Methods for determination of breakage distribution parameters , 1972 .

[13]  J. B. Riggs An Introduction to Numerical Methods for Chemical Engineers , 1988 .

[14]  Edmond P. Byrne,et al.  Mathematical modelling of granola breakage during pipe pneumatic conveying , 2011 .

[15]  P. Moin Fundamentals of Engineering Numerical Analysis , 2001 .

[16]  Leonard G. Austin,et al.  An analysis of fine dry grinding in ball mills , 1981 .

[17]  C. L. Prasher,et al.  Crushing and Grinding Process Handbook , 1987 .

[18]  Henri Berthiaux Analysis of grinding processes by Markov chains , 2000 .

[19]  Ecevit Bilgili,et al.  Numerical Simulation of Open‐Circuit Continuous Mills Using a Non‐Linear Population Balance Framework: Incorporation of Non‐First‐Order Effects , 2005 .

[20]  Brian Scarlett,et al.  Population balances for particulate processes - a volume approach. , 2002 .

[21]  Leonard G. Austin,et al.  The back-calculation of specific rates of breakage and non-normalized breakage distribution parameters from batch grinding data , 1977 .

[22]  G. Tanovic,et al.  Predicting the size and compositional distributions of wheat flour stocks following first break roller milling using the breakage matrix approach , 2006 .

[23]  Ecevit Bilgili,et al.  Population balance modeling of non-linear effects in milling processes , 2005 .

[24]  L. G. Austin,et al.  The kinetics and shape factors of ultrafine dry grinding in a laboratory tumbling ball mill , 1990 .

[25]  Fracture kinetics of particle bed comminution — ramifications for fines production and mill optimization , 1999 .

[26]  D. Verkoeijen The measurement and modeling of granulation processes: A turbo view , 2001 .

[27]  Anton P. J. Middelberg,et al.  A mathematical model for Escherichia coli debris size reduction during high pressure homogenisation based on grinding theory , 1997 .

[28]  M. C. Williams,et al.  Problems in population balance modeling of wet grinding , 1992 .

[29]  E. Bilgili,et al.  Nonlinear effects in particulate processes , 2005 .

[30]  Lawrence F. Shampine,et al.  Solving ODEs with MATLAB , 2002 .

[31]  B. Epstein The mathematical description of certain breakage mechanisms leading to the logarithmico-normal distribution , 1947 .

[32]  THE ANALYSIS OF REPEATED BREAKAGE EVENTS AS AN EQUIVALENT RATE PROCESS , 1990 .

[33]  Dominique Bertrand,et al.  Evaluation of the transition matrix for comminuting pea seeds in an impact mill using a linear neural network , 1999 .

[34]  T. G. Callcott,et al.  A matrix analysis of processes involving particle assemblies , 1956, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[35]  Leonard G. Austin,et al.  Introduction to the mathematical description of grinding as a rate process , 1971 .