A quasilinear Keller-Segel model with saturated discontinuous advection

Maria Gualdani,Mikel Ispizua,Nicola Zamponi

2024-07-08

Abstract:We consider the singular limit of a chemotaxis model of bacterial collective motion recently introduced in <a class="link-https" data-arxiv-id="2009.11048" href="https://arxiv.org/abs/2009.11048">arXiv:2009.11048</a> [math.AP]. The equation models aggregation-diffusion phenomena with advection that is discontinuous and depends sharply on the gradient of the density itself. The quasi-linearity of the problem poses major challenges in the construction of the solution and complications arise in the proof of regularity. Our method overcomes these obstacle by relying solely on entropy inequalities and the theory of monotone operators. We provide existence, uniqueness and smoothing estimates in any dimensional space.

Analysis of PDEs,Mathematical Physics

What problem does this paper attempt to address?

### Problems the paper attempts to solve

This paper aims to study the singular limit problem of a quasilinear Keller - Segel model with a saturated discontinuous convection term. Specifically, the author considers the aggregation - diffusion phenomenon in bacterial collective motion and introduces a new model to describe this phenomenon. The feature of this model is that the convection term is discontinuous and depends on the sharp change of the density gradient. This quasilinear property poses significant challenges to the existence and regularity of solutions.

### Model background

The Keller - Segel model is widely used to explain the pulse formation and propagation phenomena in bacterial colonies driven by chemotaxis. The classical Keller - Segel model describes how cell density moves towards nutrient \( N \) and chemical molecule \( S \) (chemotactic factor), and its speed is proportional to the gradients of these quantities. However, new models proposed in recent years deviate mathematically from the classical Keller - Segel framework. The speed (or flux) of cells is saturated and is proportional to \( \nabla S / |\nabla S| \) and \( \nabla N / |\nabla N| \).

### Form of the new model

The general form of the new model is:
\[
\begin{cases}
\partial_t \rho = D_\rho \Delta \rho - \text{div}(\rho u_S + \rho u_N), \\
\partial_t S = D_S \Delta S - \alpha S + \beta \rho, \\
\partial_t N = D_N \Delta N - \gamma \rho N,
\end{cases}
\]
where
\[
u_S = \chi_S J(\partial_t S, |\nabla S|) \frac{\nabla S}{|\nabla S|}, \quad u_N = \chi_N J(\partial_t N, |\nabla N|) \frac{\nabla N}{|\nabla N|}.
\]

### Main problems studied

In this paper, the author studies the quasilinear parabolic problem obtained from the above model through the singular limit \( \alpha \to \infty \):
\[
\partial_t \rho - \Delta \rho + \chi \text{div}\left( \frac{\rho \nabla \rho}{|\nabla \rho|} \right) = 0.
\]
Unlike previous models, the convection term here depends on the gradient of density, which makes it a competing term and interacts with the diffusion term. The difficulty of this model lies in controlling the diffusion coefficient, because the diffusion coefficient may change during the flow process, destroying the parabolic property of the equation.

### Main results

The author overcomes these difficulties through entropy inequality and monotone operator theory, and proves that the model has a unique solution in any - dimensional space and the solution has smoothness estimates. Specifically, the main results include:

1. **Mass conservation**: The \( L^1 \) norm of the solution remains unchanged.
2. ** \( L^p \) norm bound**: The \( L^p \) norm of the solution is bounded.
3. **Moment evolution**: The moments of the solution decay exponentially in time.
4. **Entropy bound**: The entropy of the solution is bounded.
5. ** \( L^\infty \) smoothing**: The solution is bounded in the \( L^\infty \) space.
6. ** \( L^p \) bound of the gradient**: The gradient of the solution is bounded in the \( L^p \) space.

### Method overview

To prove these results, the author uses the following methods:

1. **Monotone operator theory**: Using monotone operator theory and fixed - point arguments, the existence of the approximate problem with an additional viscosity term is proved.
2. **Relative entropy method**: Through the relative entropy method, the convergence of the nonlinear term in the vanishing viscosity limit is proved.
3. **Compactness lemma**: Using