We use cookies to distinguish you from other users and to provide you with a better experience on our websites. Close this message to accept cookies or find out how to manage your cookie settings.
To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure no-reply@cambridge.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
I examine some analytical properties of a nonlinear reaction-diffusion system that has been used to model the propagation of a wildfire. I establish global-in-time existence and uniqueness of bounded mild solutions to the Cauchy problem for this system given bounded initial data. In particular, this shows that the model does not allow for thermal blow-up. If the initial temperature and fuel density also satisfy certain integrability conditions, the spatial $L^2$-norms of these global solutions are uniformly bounded in time. Additionally, I use a bootstrap argument to show that small initial temperatures give rise to solutions that decay to zero as time goes to infinity, proving the existence of initial states that do not develop into traveling combustion waves.
In this article, we investigate a free boundary problem for the Lotka–Volterra model consisting of an invasive species with density u and a native species with density v in one dimension. We assume that v undergoes diffusion and growth in $[0,+\infty )$, and u invades into the environment with spreading front $x=h(t)$ satisfying free boundary condition $h'(t)=-u_x(t,h(t))-\alpha $ for some decay rate $\alpha>0$, this is caused by the bad environment at the boundary. When u is an inferior competitor, $u(t,x)$ and $h(t)$ tend to 0 within a finite time, while another specie $v(t,x)$ tends to a stationary $\Lambda (x)$ defined on the half-line. When u is a superior competitor, we have a trichotomy result: spreading of u and vanishing of v (i.e., as $t \to +\infty $, $h(t)$ goes to $+\infty $ and $(u,v)\to (\Lambda ,0)$); the transition case (i.e., as $t \to +\infty $, $(u,v)\to (w_\alpha , \eta _\alpha )$, $h(t)$ tends to a finite point); vanishing of u and spreading of v (i.e., $u(t,x)$ and $h(t)$ tends to 0 within a finite time, $v(t,x)$ converges to $\Lambda (x)$). Additionally, we show that this trichotomy result depends on the initial data $u(0,x)$.
In the second part of this series of papers, we address the same evolution problem that was considered in part 1 (see [16]), namely the nonlocal Fisher-KPP equation in one spatial dimension,
\begin{equation*} u_t = D u_{xx} + u(1-\phi *u), \end{equation*}
where $\phi *u$ is a spatial convolution with the top hat kernel, $\phi (y) \equiv H\left (\frac {1}{4}-y^2\right )$, except that now we modify this to an associated initial-boundary value problem on the finite spatial interval $[0,a]$ rather than the whole real line. Boundary conditions are required at the end points of the interval, and we address the situations when these are of either Dirichlet or Neumann type. This model is a natural extension of the classical Fisher-KPP model, with the introduction of the simplest possible nonlocal effect into the saturation term. Nonlocal reaction-diffusion models arise naturally in a variety of (frequently biological or ecological) contexts, and as such it is of fundamental interest to examine their properties in detail, and to compare and contrast these with the well known properties of the classical Fisher-KPP model.
For multi-scale differential equations (or fast–slow equations), one often encounters problems in which a key system parameter slowly passes through a bifurcation. In this article, we show that a pair of prototypical reaction–diffusion equations in two space dimensions can exhibit delayed Hopf bifurcations. Solutions that approach attracting/stable states before the instantaneous Hopf point stay near these states for long, spatially dependent times after these states have become repelling/unstable. We use the complex Ginzburg–Landau equation and the Brusselator models as prototypes. We show that there exist two-dimensional spatio-temporal buffer surfaces and memory surfaces in the three-dimensional space-time. We derive asymptotic formulas for them for the complex Ginzburg–Landau equation and show numerically that they exist also for the Brusselator model. At each point in the domain, these surfaces determine how long the delay in the loss of stability lasts, that is, to leading order when the spatially dependent onset of the post-Hopf oscillations occurs. Also, the onset of the oscillations in these partial differential equations is a hard onset.
This paper is focused on the existence and uniqueness of nonconstant steady states in a reaction–diffusion–ODE system, which models the predator–prey interaction with Holling-II functional response. Firstly, we aim to study the occurrence of regular stationary solutions through the application of bifurcation theory. Subsequently, by a generalized mountain pass lemma, we successfully demonstrate the existence of steady states with jump discontinuity. Furthermore, the structure of stationary solutions within a one-dimensional domain is investigated and a variety of steady-state solutions are built, which may exhibit monotonicity or symmetry. In the end, we create heterogeneous equilibrium states close to a constant equilibrium state using bifurcation theory and examine their stability.
This article is concerned with the spreading speed and traveling waves of a lattice prey–predator system with non-local diffusion in a periodic habitat. With the help of an associated scalar lattice equation, we derive the invasion speed for the predator. More specifically, when the dispersal kernel of the predator is exponentially bounded, the invasion speed is finite and can be characterized in terms of principal eigenvalues; while the dispersal kernel is algebraically decaying, the invasion speed is infinite and the accelerated spreading rate is obtained. Furthermore, the existence and non-existence of traveling waves connecting the semi-equilibrium point to a uniformly persistent state are established.
In this paper, we report the spatiotemporal dynamics of an intraguild predation (IGP)-type predator–prey model incorporating harvesting and prey-taxis. We first discuss the local and global existence of the classical solutions in N-dimensional space. It is found that the model has a global classical solution when controlling the prey-taxis coefficient in a certain range. Thereafter, we focus on the existence of the steady-state bifurcation. Moreover, we theoretically investigate the properties of the bifurcating solution near the steady-state bifurcation critical threshold. As a consequence, the spatial pattern formation of this model can be theoretically confirmed. Importantly, by means of rigorous theoretical derivation, we provide discriminant criteria on the stability of the bifurcating solution. Finally, the complicated patterns are numerically displayed. It is demonstrated that the harvesting and prey-taxis significantly affect the pattern formation of this IGP-type predator–prey model. Our main results of this paper reveal that (i) The repulsive prey-taxis could destabilize the spatial homogeneity, while the attractive prey-taxis effect and self-diffusion will stabilize the spatial homogeneity of this model. (ii) Numerical results suggest that over-harvesting for prey or predators is not advisable, it can lead to an ecological imbalance due to a significant reduction in population numbers. However, harvesting within a certain range is a feasible approach.
are obtained, in the range of exponents $p\gt 1$, $\sigma \ge -2$. More precisely, we establish conditions fulfilled by the initial data in order for the solutions to either blow-up in finite time or decay to zero as $t\to \infty$ and, in the latter case, we also deduce decay rates and large time behavior. In the limiting case $\sigma =-2$, we prove the existence of non-trivial, non-negative solutions, in stark contrast to the homogeneous case. A transformation to a generalized Fisher–KPP equation is derived and employed in order to deduce these properties.
This paper is concerned with a predator–prey system with hunting cooperation and prey-taxis under homogeneous Neumann boundary conditions. We establish the existence of globally bounded solutions in two dimensions. In three or higher dimensions, the global boundedness of solutions is obtained for the small prey-tactic coefficient. By using hunting cooperation and prey species diffusion as bifurcation parameters, we conduct linear stability analysis and find that both hunting cooperation and prey species diffusion can drive the instability to induce Hopf, Turing and Turing–Hopf bifurcations in appropriate parameter regimes. It is also found that prey-taxis is a factor stabilizing the positive constant steady state. We use numerical simulations to illustrate various spatiotemporal patterns arising from the abovementioned bifurcations including spatially homogeneous and inhomogeneous time-periodic patterns, stationary spatial patterns and chaotic fluctuations.
We investigate a recent model proposed in the literature elucidating patterns driven by chemotaxis, similar to viscous fingering phenomena. Notably, this model incorporates a singular advection term arising from a modified formulation of Darcy’s law. It is noteworthy that this type of advection can also be well interpreted as a description of a radial fluid flow source surrounding an aggregation of cells. For the two-dimensional scenario, we establish a precise threshold delineating between blow-up and global solution existence. This threshold is contingent upon the pressure magnitude and the initial total mass of the aggregating cells.
We study the global well-posedness and uniform boundedness of a two-dimensional reaction–advection–diffusion system with nonlinear advection. This strongly coupled system of nonlinear partial differential equations represents the continuum of a 2D lattice model designed to describe residential burglary, where each location is characterised by a tractability value that varies in both space and time. We show that the model with sublinear advection enhancement is globally well-posed, with a unique solution that is classical and uniformly bounded in time. Our results provide valuable insights into the development of urban crime models with nonlinear advection enhancements, making them suitable for broader applications, including nonlocal or heterogeneous near-repeat victimisation effects.
In this paper, we prove the global exstence of weak solutions for a porous medium dynamics of m species moving between two domains separated by a zero-thickness membrane. On this membrane, Kedem–Katchalsky conditions are considered, and the study is characterized by natural structural conditions applied to the nonlinear reactive terms. The global existence is established under the assumption that these reactive terms are bounded in $L^1$. This problem has already been analyzed in the linear diffusion case by Ciavolella and Perthame in Ciavolella and Perthame (2021, Journal of Evolution Equations 21, 1513–1540). The present work constitutes an extension for nonlinear diffusion, particularly of the porous medium type, in the form $\partial _t v_i - \Delta v_i^{r_i} = R_i$, for an exponent $r_i < 2$. The case $r_i \geq 2$ remains an open problem. This paper is an adaptation of the ideas from Ciavolella and Perthame (2021, Journal of Evolution Equations 21, 1513–1540), with new strategies to overcome the appearance of nonlinearity and degeneracy in the diffusion term.
Coffee berry diseases (CBD) pose significant threats to coffee production worldwide, affecting the livelihoods of millions of farmers and the global coffee market. Fractional calculus provides a powerful framework for describing non-local and memory-dependent phenomena, making it suitable for modelling the long-range interactions inherent in CBD spread. This study aims to formulate and analyse fractional order model for CBD transmission dynamics in the sense of Atangana–Baleanu–Caputo. Fixed point theorems were utilised to test the existence and uniqueness of the model’s solutions using fractional order. The basic reproduction number was calculated utilising the next-generation matrix. The model has locally asymptotically stable equilibrium positions (disease-free and endemic). Furthermore, the Lyapunov function was used to conduct a global stability analysis of the equilibrium locations. A numerical simulation of the CBD model was created using the fractional Adam–Bashforth–Moulton approach to validate the analytical findings. Our findings contribute to the development of more accurate predictive models and inform the design of targeted interventions to mitigate the impact of CBD on coffee production systems.
We study the Cauchy problem on the real line for the nonlocal Fisher-KPP equation in one spatial dimension,
\begin{equation*} u_t = D u_{xx} + u(1-\phi *u), \end{equation*}
where $\phi *u$ is a spatial convolution with the top hat kernel, $\phi (y) \equiv H\left (\frac{1}{4}-y^2\right )$. After observing that the problem is globally well-posed, we demonstrate that positive, spatially periodic solutions bifurcate from the spatially uniform steady state solution $u=1$ as the diffusivity, $D$, decreases through $\Delta _1 \approx 0.00297$ (the exact value is determined in Section 3). We explicitly construct these spatially periodic solutions as uniformly valid asymptotic approximations for $D \ll 1$, over one wavelength, via the method of matched asymptotic expansions. These consist, at leading order, of regularly spaced, compactly supported regions with width of $O(1)$ where $u=O(1)$, separated by regions where $u$ is exponentially small at leading order as $D \to 0^+$. From numerical solutions, we find that for $D \geq \Delta _1$, permanent form travelling waves, with minimum wavespeed, $2 \sqrt{D}$, are generated, whilst for $0 \lt D \lt \Delta _1$, the wavefronts generated separate the regions where $u=0$ from a region where a steady periodic solution is created via a distinct periodic shedding mechanism acting immediately to the rear of the advancing front, with this mechanism becoming more pronounced with decreasing $D$. The structure of these transitional travelling wave forms is examined in some detail.
We provide well-posedness results for nonlinear parabolic partial differential equations (PDEs) given by reaction–diffusion equations describing the concentration of oxygen in encapsulated cells. The cells are described in terms of a core and a shell, which introduces a discontinuous diffusion coefficient as the material properties of the core and shell differ. In addition, the cells are subject to general nonlinear consumption of oxygen. As no monotonicity condition is imposed on the consumption, monotone operator theory cannot be used. Moreover, the discontinuity in the diffusion coefficient bars us from applying classical results on strong solutions. However, by directly applying a Galerkin method, we obtain uniqueness and existence of the strong form solution. These results provide the basis to study the dynamics of cells in critical states.
This article offers an advanced and novel investigation into the intricate propagation dynamics of the Belousov–Zhabotinsky system with non-local delayed interaction, which exhibits dynamical transition structure from bistable to monostable. We first solved the enduring open problem concerning the existence, uniqueness and the speed sign of the bistable travelling waves. In the monostable case, we developed and derived new results for the minimal wave speed selection, which, as an application, further improved the existing investigations on pushed and pulled wavefronts. Our results can provide new estimate to the minimal speed as well as to the determinacy of the transition parameters. Moreover, these results can be directly applied to standard localised models and delayed reaction diffusion models by choosing appropriate kernel functions.
In a smoothly bounded domain $\Omega \subset \mathbb{R}^n$, $n\ge 1$, this manuscript considers the homogeneous Neumann boundary problem for the chemotaxis system
\begin{eqnarray*} \left \{ \begin{array}{l} u_t = \Delta u - \nabla \cdot (u\nabla v), \\[5pt] v_t = \Delta v + u - \alpha uv, \end{array} \right . \end{eqnarray*}
with parameter $\alpha \gt 0$ and with coincident production and uptake of attractants, as recently emphasized by Dallaston et al. as relevant for the understanding of T-cell dynamics.
It is shown that there exists $\delta _\star =\delta _\star (n)\gt 0$ such that for any given $\alpha \ge \frac{1}{\delta _\star }$ and for any suitably regular initial data satisfying $v(\cdot, 0)\le \delta _\star$, this problem admits a unique classical solution that stabilizes to the constant equilibrium $(\frac{1}{|\Omega |}\int _\Omega u(\cdot, 0), \, \frac{1}{\alpha })$ in the large time limit.
Infection mechanism plays a significant role in epidemic models. To investigate the influence of saturation effect, a nonlocal (convolution) dispersal susceptible-infected-susceptible epidemic model with saturated incidence is considered. We first study the impact of dispersal rates and total population size on the basic reproduction number. Yang, Li and Ruan (J. Differ. Equ. 267 (2019) 2011–2051) obtained the limit of basic reproduction number as the dispersal rate tends to zero or infinity under the condition that a corresponding weighted eigenvalue problem has a unique positive principal eigenvalue. We remove this additional condition by a different method, which enables us to reduce the problem on the limiting profile of the basic reproduction number into that of the spectral bound of the corresponding operator. Then we establish the existence and uniqueness of endemic steady states by a equivalent equation and finally investigate the asymptotic profiles of the endemic steady states for small and large diffusion rates to provide reference for disease prevention and control, in which the lack of regularity of the endemic steady state and Harnack inequality makes the limit function of the sequence of the endemic steady state hard to get. Finally, we find whether lowing the movements of susceptible individuals can eradicate the disease or not depends on not only the sign of the difference between the transmission rate and the recovery rate but also the total population size, which is different from that of the model with standard or bilinear incidence.
Flowering plants depend on some animals for pollination and contribute to nourish the animals in natural environments. We call these animals pollinators and build a plants-pollinators cooperative model with impulsive effect on a periodically evolving domain. Next, we define the ecological reproduction index for single plant model and plants-pollinators system, respectively, whose threshold dynamics, including the extinction, persistence and coexistence, is established by the method of upper and lower solutions. Theoretical analysis shows that a large domain evolution rate has a positive influence on the survival of pollinators whether or not the impulsive effect occurs, and the pulse eliminates the pollinators even when the evolution rate is high. Moreover, some selective numerical simulations are still performed to explain our theoretical results.
In this paper, we analyse Turing instability and bifurcations in a host–parasitoid model with nonlocal effect. For a ordinary differential equation model, we provide some preliminary analysis on Hopf bifurcation. For a reaction–diffusion model with local intraspecific prey competition, we first explore the Turing instability of spatially homogeneous steady states. Next, we show that the model can undergo Hopf bifurcation and Turing–Hopf bifurcation, and find that a pair of spatially nonhomogeneous periodic solutions is stable for a (8,0)-mode Turing–Hopf bifurcation and unstable for a (3,0)-mode Turing–Hopf bifurcation. For a reaction–diffusion model with nonlocal intraspecific prey competition, we study the existence of the Hopf bifurcation, double-Hopf bifurcation, Turing bifurcation, and Turing–Hopf bifurcation successively, and find that a spatially nonhomogeneous quasi-periodic solution is unstable for a (0,1)-mode double-Hopf bifurcation. Our results indicate that the model exhibits complex pattern formations, including transient states, monostability, bistability, and tristability. Finally, numerical simulations are provided to illustrate complex dynamics and verify our theoretical results.