Comparison of two finite element schemes for a chemo-repulsion system with quadratic production

TL;DR

This paper compares two finite element schemes for a chemo-repulsion model with quadratic production, demonstrating improved stability and properties of the schemes, along with analysis of their long-term behavior and numerical performance.

Contribution

It introduces and analyzes two fully discrete FE schemes, proving their energy stability and positivity properties, and compares their effectiveness through numerical simulations.

Findings

Scheme UV is energy-stable under a compatibility condition.

Scheme US$_ ext{ε}$ is energy-stable with respect to a modified energy.

Both schemes exhibit exponential convergence to constant states.

Abstract

In this paper we propose two fully discrete Finite Elements (FE) schemes for a repulsive chemotaxis model with quadratic production term. The first one (called scheme UV) corresponds to the backward Euler in time with FE in space approximation; while the second one (called scheme US$_\varepsilon$) is obtained as a modification of the scheme US proposed by [Guill\'en-Gonz\'alez et al.], by applying a regularization procedure. We prove that the schemes UV and US$_\varepsilon$ have better properties than the FE scheme US. Specifically, we prove that, unlike the scheme US, the scheme UV is energy-stable in the primitive variables of the model, under a "compatibility" condition on the FE spaces. On the other hand, the scheme US$_\varepsilon$ is energy-stable with respect to the same modified energy of the scheme US, and an "approximated positivity" property holds (which is not possible to prove for the schemes US and UV). Additionally, we study the well-posedness of the schemes and the long time behaviour obtaining exponential convergence to constant states. Finally, we compare the numerical schemes throughout several numerical simulations.