The numerical scheme for which we have proven the analytical results below, and which we used for our numerical simulations is given by
We proved that such a scheme is consistent with the Parabolic-Elliptic Keller-Segel equations. The statement of the theorem is as follows
We proved two major results for the numerical scheme given above which together give convergence of the scheme to a weak solution of the Keller-Segel equations. The first involves proving that the scheme has a unique solution at each time step. For this we proved the following theorem, which also gives us conservation of mass, positivity of the solution and a uniform bound on the discrete norm of the approximation to the chemical concentration . The proof of existence involves a linearisation of the scheme together with an application of the Brouwer fixed point theorem and was adapted from the proof given by Filbet in the finite volume setting . In the following we have made use of the trapezium rule, so for notational convenience we define the ''trapezium sum" as follows
The statement of the first theorem is
The second result was split into two parts, the first of which was to prove convergence in some sense of the unique solution to the scheme to a couple as the time step and grid size approach zero, and finally we proved that this is a weak solution to the Keller-Segel equations. For this we first defined discrete functions and associated to the solution of the scheme as piecewise constant functions taking corresponding values at grid points and discrete time steps. We also define piecewise constant approximations and to the derivatives with on . We then needed some a-priori estimates on discrete norms of and and due to the form of the non-linearity in the parabolic equation, in order to pass to the limit we required strong convergence of to but only weak convergence of to and of the approximations to the derivatives, and to and . This result can be summarised in the following theorem, the proof of which was mostly adapted from .