A convergent structure-preserving finite-volume scheme for the Shigesada-Kawasaki-Teramoto population system

Abstract

An implicit Euler finite-volume scheme for an n-species population crossdiffusion system of Shigesada-Kawasaki-Teramoto-type in a bounded domain with no-flux boundary conditions is proposed and analyzed. The scheme preserves the formal gradientflow or entropy structure and preserves the nonnegativity of the population densities. The key idea is to consider a suitable mean of the mobilities in such a way that a discrete chain rule is fulfilled and a discrete analog of the entropy inequality holds. The existence of finite-volume solutions, the convergence of the scheme, and the large-time asymptotics to the constant steady state are proven. Furthermore, numerical experiments in one and two space dimensiona for two and three species are presented. The results are valid for a more general class of cross-diffusion systems satisfying some structural conditions.