Adhesion and volume filling in one-dimensional population dynamics under no-flux boundary condition

Hyung Jun Choi,Seonghak Kim,Youngwoo Koh
2024-09-07
Abstract:We study the (generalized) one-dimensional population model developed by Anguige \& Schmeiser [1], which reflects cell-cell adhesion and volume filling under no-flux boundary condition. In this generalized model, depending on the adhesion and volume filling parameters $\alpha,\beta\in[0,1],$ the resulting equation is classified into six types. Among these, we focus on the type exhibiting strong effects of both adhesion and volume filling, which results in a class of advection-diffusion equations of the forward-backward-forward type. For five distinct cases of initial maximum, minimum and average population densities, we derive the corresponding patterns for the global behavior of weak solutions to the initial and no-flux boundary value problem. Due to the presence of a negative diffusion regime, we indeed prove that the problem is ill-posed and admits infinitely many global-in-time weak solutions, with the exception of one specific case of the initial datum. This nonuniqueness is inherent in the method of convex integration that we use to solve the Dirichlet problem of a partial differential inclusion arising from the ill-posed problem.
Analysis of PDEs
What problem does this paper attempt to address?
### Problems the paper attempts to solve This paper aims to study the behaviors of cell - to - cell adhesion and volume - filling effects in a one - dimensional population dynamics model under no - flux boundary conditions. Specifically, the paper focuses on the fact that when both cell - to - cell adhesion and volume - filling effects are strong, the model equation will exhibit a forward - backward - forward type of convection - diffusion equation. Such an equation may be ill - posed under certain initial conditions, that is, there exist multiple global weak solutions. The main objectives of the paper are: 1. **Analyze the global behaviors under different initial maximum, minimum and average population densities**: By analyzing five different initial population density situations, derive the global behavior patterns of weak solutions in these cases. 2. **Prove the ill - posedness of the problem**: Due to the existence of a negative diffusion region, the paper proves that this problem is ill - posed under most initial conditions, that is, there exist infinitely many global weak solutions. But under a specific initial condition, the problem can be well - posed. 3. **Solve the ill - posed problem using the convex integration method**: By the convex integration method, the paper solves the partial differential inclusion problem caused by the ill - posed problem, thus obtaining the existence and non - uniqueness of weak solutions. ### Specific problem description The model equation discussed in the paper is: \[ \begin{cases} u_t = (\rho(u))_{xx} & \text{in } \Omega\times(0, \infty), \\ u = u_0 & \text{on } \Omega\times\{t = 0\}, \\ \sigma(u)u_x = 0 & \text{on } \partial\Omega\times(0, \infty), \end{cases} \] where \(u = u(x, t)\) represents the population density at position \(x\) and time \(t\), \(\Omega=(0, 1)\) is a unit - sized suitable habitat, \(u_0 = u_0(x)\) is the initial population density, and the flux function \(\rho(s)\) is defined as: \[ \rho(s)=\alpha\beta s^3 - 2\alpha s^2 + s\quad(s\in\mathbb{R}), \] The parameters \(\alpha, \beta\in[0, 1]\) represent the adhesion and volume - filling parameters respectively, and the diffusion coefficient \(\sigma(s)\) is: \[ \sigma(s)=3\alpha\beta s^2 - 4\alpha s + 1\quad(s\in\mathbb{R}). \] ### Main results The main results of the paper are as follows: 1. **Smooth stabilization** (Case (i)): If the initial maximum density \(M_0 < s_0^-\) or the initial minimum density \(m_0 > s_0^+\), then there exists a global classical solution \(u\), which lies in the forward region \([0, s_0^-)\cup(s_0^+, 1]\) at all times and finally stabilizes to the initial average density \(\bar{u}_0\). 2. **Density mixing in finite time** (Case (ii - 1) and Case (iii)): If the initial maximum density \(M_0\geq s_0^-\) but the initial average density \(\bar{u}_0 < s_0^-\) or \(\bar{u}_0 > s_0^+\), then fine - scale density mixing will occur in finite time, but finally these weak solutions will become smooth and stabilize to \(\bar{u}_0\). 3. **Persistent density mixing** (Case (iv)): If the initial average density \(\bar{u}_0\) lies in the backward region \([s_0^-, s_0^+]\), then the weak solution will always experience fine - scale density mixing, and this mixing will persist throughout the habitat. Through these results, the paper reveals that under strong adhesion and