We study a diffusive Lotka-Volterra competition system with advection under Neumann boundary conditions. Our system models a competition relationship that one species escape from the region of high population density of their competitors in order to avoid competition. We establish the global existence of bounded classical solutions to the system over one-dimensional finite domains. For multi-dimensional domains, globally bounded classical solutions are obtained for a parabolic-elliptic system under proper assumptions on the system parameters. These global existence results make it possible to study bounded steady states in order to model species segregation phenomenon. We then investigate the one-dimensional stationary problem. Through bifurcation theory, we obtain the existence of nonconstant positive steady states, which are small perturbations from the positive equilibrium; we also rigourously study the stability of these bifurcating solutions when diffusion coefficients of the escaper and its competitor are large and small respectively. In the limit of large advection rate, we show that the reaction-advection-diffusion system converges to a shadow system involving the competitor population density and an unknown positive constant. Existence and stability of positive nonconstant solutions to the shadow system have also been obtained through bifurcation theories. Finally, we construct infinitely many single interior transition layers to the shadow system when crowding rate of the escapers and diffusion rate of their interspecific competitors are sufficiently small. The transition-layer solutions can be used to model the interspecific segregation phenomenon.
1 January 2017 in Discrete and Continuous Dynamical Systems
This paper concerns pattern formation in a class of reaction-advection-diffusion systems modeling the population dynamics of two predators and one prey. We consider the biological situation that both predators forage along the population density gradient of the preys which can defend themselves as a group. We prove the global existence and uniform boundedness of positive classical solutions for the fully parabolic system over a bounded domain with space dimension $ N=1,2 $ and for the parabolic-parabolic-elliptic system over higher space dimensions. Linearized stability analysis shows that prey-taxis stabilizes the positive constant equilibrium if there is no group defense while it destabilizes the equilibrium otherwise. Then we obtain stationary and time-periodic nontrivial solutions of the system that bifurcate from the positive constant equilibrium. Moreover, the stability of these solutions is also analyzed in detail which provides a wave mode selection mechanism of nontrivial patterns for this strongly coupled system. Finally, we perform numerical simulations to illustrate and support our theoretical results.
1 January 2017 in Discrete and Continuous Dynamical Systems-Series B
This paper investigates the formation of time-periodic and stable patterns of a two-competing-species Keller-Segel chemotaxis model with a focus on the effect of cellular growth. We carry out rigorous Hopf bifurcation analysis to obtain the bifurcation values, spatial profiles and time period associated with these oscillating patterns. Moreover, the stability of the periodic solutions is investigated and it provides a selection mechanism of stable time-periodic mode which suggests that only large domains support the formation of these periodic patterns. Another main result of this paper reveals that cellular growth is responsible for the emergence and stabilization of the oscillating patterns observed in the 3×3 system, while the system admits a Lyapunov functional in the absence of cellular growth. Global existence and boundedness of the system in 2D are proved thanks to this Lyapunov functional. Finally, we provide some numerical simulations to illustrate and support our theoretical findings.
We study a one--dimensional quasilinear system proposed by J. Tello and M. Winkler [27] which models the population dynamics of two competing species attracted by the same chemical. The kinetic terms of the interacting species are chosen to be of the Lotka--Volterra type and the boundary conditions are of homogeneous Neumann type which represent an enclosed domain. We prove the global existence and boundedness of classical solutions to the fully parabolic system. Then we establish the existence of nonconstant positive steady states through bifurcation theory. The stability or instability of the bifurcating solutions is investigated rigorously. Our results indicate that small intervals support stable monotone positive steady states and large intervals support nonmonotone steady states. Finally, we perform extensive numerical studies to demonstrate and verify our theoretical results. Our numerical simulations also illustrate the formation of stable steady states and time--periodic solutions with various interesting spatial structures.
26 May 2016 in European Journal of Applied Mathematics
In this paper, we study two PDEs that generalize the urban crime model proposed by Short et al. (2008 Math. Models Methods Appl. Sci.18, 1249–1267). Our modifications are made under assumption of the spatial heterogeneity of both the near-repeat victimization effect and the dispersal strategy of criminal agents. We investigate pattern formations in the reaction–advection–diffusion systems with non-linear diffusion over multi-dimensional bounded domains subject to homogeneous Neumann boundary conditions. It is shown that the positive homogeneous steady state loses its stability as the intrinsic near-repeat victimization rate ε decreases and spatially inhomogeneous steady states emerge through bifurcation. Moreover, we find the wavemode selection mechanism through rigorous stability analysis of these non-trivial spatial patterns, which shows that the only stable pattern must have a wavenumber that maximizes the bifurcation value. Based on this wavemode selection mechanism, we will be able to predict the formation of stable aggregates of the house attractiveness and criminal population density, at least when the diffusion rate ε is around the principal bifurcation value. Our theoretical results also suggest that large domains support more stable aggregates than small domains. Finally, we perform extensive numerical simulations over 1D intervals and 2D squares to illustrate and verify our theoretical findings. Our numerics also demonstrate the formation of other interesting patterns in these models such as the merging of two interior spikes and the emerging of new spikes, etc. These non-trivial solutions can model the well-observed phenomenon of aggregation in urban criminal activities.