Finite element approximation and very weak solution existence in a two-dimensional, degenerate Keller-Segel model

This paper is devoted to the design and analysis of a numerical algorithm for approximating solutions of a degenerate cross-diffusion system, which models particular instances of taxis-type migration processes under local sensing mechanisms. The degeneracy leads to solutions that are very weak due to the low regularity themselves. Specifically, the solutions satisfy pointwise bounds (such as positivity and the maximum principle), integrability (such as mass conservation), and dual a priori estimates. The proposed numerical scheme combines a finite element spatial discretization with Euler time stepping. The discrete solutions preserve the above-mentioned properties at the discrete level, enabling the derivation of compactness arguments and the convergence (up to a subsequence) of the numerical solutions to a very weak solution of the continuous problem on two-dimensional polygonal domains.