A GENERALIZED KINETIC MODEL OF SEDIMENTATION OF POLYDISPERSE SUSPENSIONS WITH A CONTINUOUS PARTICLE SIZE DISTRIBUTION

Polydisperse suspensions with particles of a finite number N of size classes have been widely studied in laboratory experiments. However, in most real-world applications the particle sizes are distributed continuously. In this paper, a well-studied one-dimensional kinematic model for batch sedimentation of polydisperse suspensions of small equal-density spheres is extended to suspensions with a continuous particle size distribution. For this purpose, the phase density function Φ = Φ(t, x, ξ), where ξ ∈ [0, 1] is the normalized squared size of the particles, is introduced, whose integral with respect to ξ on an interval [ξ1, ξ2] is equivalent to the volume fraction at (t, x) occupied by particles of that size range. Combining the Masliyah–Lockett–Bassoon (MLB) model for the solid-fluid relative velocity for each solids species with the concept of phase density function yields a scalar, first-order equation for Φ, namely the equation of the generalized kinetic theory. Three numerical schemes for the solution of this equation are introduced, and a numerical example and an L1 error study show that one of these schemes introduces less numerical diffusion and less spurious oscillations near discontinuities than the others. Several numerical examples illustrate the simulated behavior of this kind of suspensions. Numerical results also illustrate the solution of an eigenvalue problem associated with the equation of the generalized kinetic theory.

[1]  Raimund Bürger,et al.  Neumann Problems for Quasi-linear Parabolic Systems Modeling Polydisperse Suspensions , 2006, SIAM J. Math. Anal..

[2]  S. Wong,et al.  Hyperbolicity and kinematic waves of a class of multi-population partial differential equations , 2006, European Journal of Applied Mathematics.

[3]  Monika Bargieł,et al.  An extension of the particle-based approach to simulating the sedimentation of polydisperse suspensions , 2006 .

[4]  Alice J. Kozakevicius,et al.  Adaptive multiresolution WENO schemes for multi-species kinematic flow models , 2007, J. Comput. Phys..

[5]  Robert H. Davis,et al.  Classification of concentrated suspensions using inclined settlers , 1996 .

[6]  Andrea Tosin,et al.  Mathematical modeling of vehicular traffic: a discrete kinetic theory approach , 2007 .

[7]  P. Rosin.  The Laws Governing the Fineness of Powdered Coal , 1933 .

[8]  Bandaru V. Ramarao,et al.  On instabilities arising during sedimentation of two-component mixtures of solids , 1984, Journal of Fluid Mechanics.

[9]  Jacob H. Masliyah,et al.  ]Hindered settling in a multi-species particle system , 1979 .

[10]  Shamsul Qamar,et al.  Numerical solution of population balance equations for nucleation, growth and aggregation processes , 2007, Comput. Chem. Eng..

[11]  G. Anestis,et al.  Sediment composition due to settling of particles of different sizes , 1985 .

[12]  Chi-Wang Shu,et al.  A weighted essentially non-oscillatory numerical scheme for a multi-class Lighthill-Whitham-Richards traffic flow model , 2003 .

[13]  Rinaldo M. Colombo,et al.  An $n$-populations model for traffic flow , 2003, European Journal of Applied Mathematics.

[14]  Monika Bargieł,et al.  Stability of tridisperse suspensions , 2007 .

[15]  H. P. Greenspan,et al.  On hindered settling of particles of different sizes , 1982 .

[16]  Shamsul Qamar,et al.  Adaptive high-resolution schemes for multidimensional population balances in crystallization processes , 2007, Comput. Chem. Eng..

[17]  Raimund Bürger,et al.  Strongly Degenerate Parabolic-Hyperbolic Systems Modeling Polydisperse Sedimentation with Compression , 2003, SIAM J. Appl. Math..

[18]  R. A. Ford,et al.  Simulation of Sedimentation of Polydisperse Suspensions: A Particle-Based Approach , 2005 .

[19]  Shamsul Qamar,et al.  A comparative study of high resolution schemes for solving population balances in crystallization , 2006, Comput. Chem. Eng..

[20]  Raimund Bürger,et al.  Model Equations and Instability Regions for the Sedimentation of Polydisperse Suspensions of Spheres , 2002 .

[21]  Nicola Bellomo,et al.  ON THE MATHEMATICAL THEORY OF VEHICULAR TRAFFIC FLOW I: FLUID DYNAMIC AND KINETIC MODELLING , 2002 .

[22]  Francis Filbet,et al.  Numerical Simulation of the Smoluchowski Coagulation Equation , 2004, SIAM J. Sci. Comput..

[23]  Sohrab Rohani,et al.  Dynamic and steady-state sedimentation of polydisperse suspension and prediction of outlets particle-size distribution , 2004 .

[24]  Stefan Heinrich,et al.  Improved accuracy and convergence of discretized population balance for aggregation: The cell average technique , 2006 .

[25]  B. Perthame,et al.  Derivation of hyperbolic models for chemosensitive movement , 2005, Journal of mathematical biology.

[26]  Nicola Bellomo,et al.  Generalized kinetic (Boltzmann) models: mathematical structures and applications , 2002 .

[27]  Z. Ha,et al.  Settling Velocities of Polydisperse Concentrated Suspensions , 2008 .

[28]  M. J. Lockett,et al.  Sedimentation of Binary Particle Mixtures , 1979 .

[29]  Sze Chun Wong,et al.  A weighted essentially non-oscillatory numerical scheme for a multi-class traffic flow model on an inhomogeneous highway , 2006, J. Comput. Phys..

[30]  Elena De Angelis,et al.  Modelling complex systems in applied sciences; methods and tools of the mathematical kinetic theory for active particles , 2006, Math. Comput. Model..

[31]  J. D. Towers,et al.  Closed-form and finite difference solutions to a population balance model of grinding mills , 2005 .

[32]  Raimund Bürger,et al.  A family of numerical schemes for kinematic flows with discontinuous flux , 2008 .

[33]  R. Bürger,et al.  Central difference solutions of the kinematic model of settling of polydisperse suspensions and three-dimensional particle-scale simulations , 2001 .

[34]  Raimund Bürger,et al.  Numerical simulation of the settling of polydisperse suspensions of spheres , 2000 .

[35]  Monika Bargieł,et al.  Simulation of sedimentation and fluidization of polydisperse suspensions via a Markov model , 2006 .

[36]  Modeling of sedimentation of polydisperse spherical beads with a broad size distribution , 2003 .

[37]  R. Colombo,et al.  Measure valued solutions to conservation laws motivated by traffic modelling , 2006, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[38]  E. Tadmor,et al.  New High-Resolution Central Schemes for Nonlinear Conservation Laws and Convection—Diffusion Equations , 2000 .

[39]  N. Bellomo,et al.  MULTICELLULAR BIOLOGICAL GROWING SYSTEMS: HYPERBOLIC LIMITS TOWARDS MACROSCOPIC DESCRIPTION , 2007 .

[40]  Raimund Bürger,et al.  On Riemann problems and front tracking for a model of sedimentation of polydisperse suspensions , 2007 .

[41]  G. Batchelor,et al.  Structure formation in bidisperse sedimentation , 1986, Journal of Fluid Mechanics.

[42]  Robert H. Davis,et al.  Spreading of the interface at the top of a slightly polydisperse sedimenting suspension , 1988, Journal of Fluid Mechanics.

[43]  Monika Bargieł,et al.  A five-parameter Markov model for simulating the paths of sedimenting particles , 2007 .

[44]  Krishnaswamy Nandakumar,et al.  Gravity separation of bidisperse suspensions: Light and heavy particle species , 1987 .