Antidiffusive Lagrangian‐remap schemes for models of polydisperse sedimentation

Abstract

Listen

Translate article

One‐dimensional models of gravity‐driven sedimentation of polydisperse suspensions with particles that belong to N size classes give rise to systems of N strongly coupled, nonlinear first‐order conservation laws for the local solids volume fractions. As the eigenvalues and eigenvectors of the flux Jacobian have no closed algebraic form, characteristic‐wise numerical schemes for these models become involved. Alternative simple schemes for this model directly utilize the velocity functions and are based on splitting the system of conservation laws into two different first‐order quasi‐linear systems, which are solved successively for each time iteration, namely, the Lagrangian and remap steps (so‐called Lagrangian‐remap [LR] schemes). This approach was advanced in (Bürger, Chalons, and Villada, SIAM J Sci Comput 35 (2013), B1341–B1368) for a multiclass Lighthill–Whitham‐Richards traffic model with nonnegative velocities. By incorporating recent antidiffusive techniques for transport equations a new version of these Lagrangian‐antidiffusive remap (L‐AR) schemes for the polydisperse sedimentation model is constructed. These L‐AR schemes are supported by a partial analysis for N = 1. They are total variation diminishing under a suitable CFL condition and therefore converge to a weak solution. Numerical examples illustrate that these schemes, including a more accurate version based on MUSCL extrapolation, are competitive in accuracy and efficiency with several existing schemes. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1109–1136, 2016