Abstract

This paper analyses a continuum model of the sedimentation of polydisperse suspensions, where the volume concentrations of the solids species are the sought unknowns. This leads to systems of conservation laws which are non-genuinely nonlinear ("non-convex") in the sense that they are only piecewise genuinely nonlinear. The solution of a Riemann problem is represented as a concatenation of elementary waves, where the linear degeneracy requires using the Liu entropy condition. After a general description of the composition of elementary waves, an algorithm for the computational construction is given. The solution of the Riemann problem is then utilized as a problem-adapted building block of a front tracking method. For the model of bidisperse suspensions, Riemann problems are classified by the location of the left and right state in the phase space of unknown variables. The front tracking method is applied to solve the initial value problem of a first-order hyperbolic 2×2 system of conservation laws describing the settling of a bidisperse suspension. The solution obtained by front tracking is compared with results obtained by an experiment and a finite difference scheme. © 2007 WILEY-VCH Verlag GmbH & Co. KGaA.