A difference scheme for a degenerating convection-diffusion-reaction system modelling continuous sedimentation

¹ CI 2MA and Departamento de Ingeniería Matemática, Facultad de Ciencias Físicas y Matemáticas, Universidad de Concepción, Casilla 160-C, Concepción, Chile
² Centre for Mathematical Sciences, Lund University, P.O. Box 118, S-221 00 Lund, Sweden

^* Corresponding author: diehl@maths.lth.se

Received: 18 October 2016
Revised: 14 August 2017
Accepted: 17 August 2017

Abstract

Continuously operated settling tanks are used for the gravity separation of solid-liquid suspensions in several industries. Mathematical models of these units form a topic for well-posedness and numerical analysis even in one space dimension due to the spatially discontinuous coefficients of the underlying strongly degenerate parabolic, nonlinear model partial differential equation (PDE). Such a model is extended to describe the sedimentation of multi-component particles that react with several soluble components of the liquid phase. The fundamental balance equations contain the mass percentages of the components of the solid and liquid phases. The equations are reformulated as a system of nonlinear PDEs that can be solved consecutively in each time step by an explicit numerical scheme. This scheme combines a difference scheme for conservation laws with discontinuous flux with an approach of numerical percentage propagation for multi-component flows. The main result is an invariant-region property, which implies that physically relevant numerical solutions are produced. Simulations of denitrification in secondary settling tanks in wastewater treatment illustrate the model and its discretization.