Abstract

A one-dimensional model of clarifier-thickener units in engineering applications can be expressed as a conservation law with a flux that is discontinuous with respect to the spatial variable. This model also includes a singular feed source. In this paper, the clarifier-thickener model studied in a previous paper (Numer. Math. 97 (2004) 25-65) is extended by a singular sink through which material is extracted from the unit. A difficulty is that in contrast to the singular source, the sink term cannot be incorporated into the flux function; rather, the sink is represented by a new non-conservative transport term. To focus on the new analytical difficulties arising due to this non-conservative term, a reduced problem is formulated, which contains the new sink term of the extended clarifier-thickener model, but not the source term and flux discontinuities. The paper is concerned with numerical methods for both models (extended and reduced) and with the well-posedness analysis for the reduced problem. For the reduced problem, a definition of entropy solutions, based on Krǔkov-type entropy functions and fluxes, is provided. Jump conditions are derived and uniqueness of the entropy solution is shown. Existence of an entropy solution is shown by proving convergence of a monotone difference scheme. Two variants of the numerical scheme are introduced. Numerical examples illustrate that all three variants converge to the entropy solution, but introduce different amounts of numerical diffusion. © 2006 Cambridge University Press.