A cross-diffusion system modelling rivaling gangs: global existence of bounded solutions and FCT stabilization for numerical simulation
Mario Fuest, Shahin Heydari
TL;DR
This work analyzes a two-species cross-diffusion model for rival gang territoriality with graffiti-based signals, proving global existence and boundedness of classical solutions under saturating production, and showing convergence to homogeneous steady states for small data. The authors transform the system to a regularized form, derive extensive a priori estimates (including L∞ and gradient bounds), and establish global existence for the original variables. They also prove that smooth heterogeneous steady states cannot exist under analytic production functions, while providing a robust FEM-FCT numerical framework to simulate the model with positivity preservation. Numerical experiments reveal diffusion-dominated regimes that homogenize and convection-dominated regimes that exhibit various segregation patterns, including complete separation under certain parameters, validating the theoretical results and highlighting the method's ability to capture complex territorial dynamics. Overall, the study advances understanding of when rival groups mix or segregate and offers a practical, stable computational approach for simulating cross-diffusion territorial models.
Abstract
For the gang territoriality model \begin{align*} \begin{cases} u_t = D_u Δu + χ_u \nabla \cdot (u \nabla w), \\ v_t = D_v Δv + χ_v \nabla \cdot (v \nabla z), \\ w_t = -w + \frac{v}{1+v}, \\ z_t = -z + \frac{u}{1+u}, \end{cases} \end{align*} where $u$ and $v$ denote the densities of two rivaling gangs which spray graffiti (with densities $z$ and $w$, respectively) and partially move away from the other gang's graffiti, we construct global, bounded classical solutions. By making use of quantitative global estimates, we prove that these solutions converge to homogeneous steady states if $\|u_0\|_{L^\infty(Ω)}$ and $\|v_0\|_{L^\infty(Ω)}$ are sufficiently small. Moreover, we perform numerical experiments which show that for different choices of parameters, the system may become diffusion- or convection-dominated, where in the former case the solutions converge toward constant steady states while in the later case nontrivial asymptotic behavior such as segregation is observed. In order to perform these experiments, we apply a nonlinear finite element flux-corrected transport method (FEM-FCT) which is positivity-preserving. Then, we treat the nonlinearities in both the system and the proposed nonlinear scheme simultaneously using fixed-point iteration.
