Stabilization and pattern formation in chemotaxis models with acceleration and logistic source

Abstract

We consider the following chemotaxis-growth system with an acceleration assumption, {u(t) = Delta u - del. (uw) + gamma (u - u(alpha)) , x is an element of Omega, t > 0, v(t) = Delta v - v + u, x is an element of Omega, t > 0, w(t) = Delta w - w + chi del v, x is an element of Omega, t > 0, under the homogeneous Neumann boundary condition for u, v and the homogeneous Dirichlet boundary condition for w in a smooth bounded domain Omega subset of R-n (n >= 1) with given parameters chi > 0, gamma >= 0 and alpha > 1. It is proved that for reasonable initial data with either n <= 3, gamma >= 0, alpha > 1 or n >= 4, gamma > 0, alpha > 1/2 + n/4, the system admits global bounded solutions, which significantly differs from the classical chemotaxis model that may have blow-up solutions in two and three dimensions. For given gamma and alpha, the obtained global bounded solutions are shown to convergence exponentially to the spatially homogeneous steady state (m, m, 0) in the large time limit for appropriately small chi, where m = 1/vertical bar Omega vertical bar integral(Omega) u(0)(x) if gamma = 0 and m = 1 if gamma > 0. Outside the stable parameter regime, we conduct linear analysis to specify possible patterning regimes. In weakly nonlinear parameter regimes, with a standard perturbation expansion approach, we show that the above asymmetric model can generate pitch-fork bifurcations which occur generically in symmetric systems. Moreover, our numerical simulations demonstrate that the model can generate rich aggregation patterns, including stationary, single merging aggregation, merging and emerging chaotic, and spatially inhomogeneous time-periodic. Some open questions for further research are discussed.