Exploring numerical blow-up phenomena for the Keller–Segel–Navier–Stokes equations

Abstract

AbstractThe Keller-Segel-Navier-Stokes system governs chemotaxis in liquid environments. This system is to be solved for the organism and chemoattractant densities and for the fluid velocity and pressure. It is known that if the total initial organism density mass is below 2πthere exist globally defined generalised solutions, but what is less understood is whether there are blow-up solutions beyond such a threshold and its optimality.Motivated by this issue, a numerical blow-up scenario is investigated. Approximate solutions computed via a stabilised finite element method founded on a shock capturing technique are such that they satisfya prioribounds as well as lower andL1(Ω) bounds for the organism and chemoattractant densities. In particular, these latter properties are essential in detecting numerical blow-up configurations, since the non-satisfaction of these two requirements might trigger numerical oscillations leading to non-realistic finite-time collapses into persistent Dirac-type measures.Our findings show that the existence threshold value 2πencountered for the organism density mass may not be optimal and hence it is conjectured that the critical threshold value 4πmay be inherited from the fluid-free Keller-Segel equations. Additionally it is observed that the formation of singular points can be neglected if the fluid flow is intensified.