Free boundary problem for a nonlocal time-periodic diffusive competition model

In this paper we consider a free boundary problem for a nonlocal time-periodic competition model. One species is assumed to adopt nonlocal dispersal, and the other one adopts mixed dispersal, which is a combination of both random dispersal and nonlocal dispersal. We first prove the global well-posedness of solutions to the free boundary problem with more general growth functions, and then discuss the spreading and vanishing phenomena. Moreover, under the weak competition condition, we study the long-time behaviors of solutions for the spreading case.

Analysis of Solutions to a Free Boundary Problem with a Nonlinear Gradient Absorption

In this paper, we investigate the blow-up rate and global existence of solutions to a parabolic system with absorption and the free boundary. By using the comparison principle and super-sub solution method, we obtain some sufficient conditions on the global existence, blow-up in finite time of solutions, and blow-up sets when blow-up phenomenon occurs. Furthermore, the global solution is bounded and uniformly tends to zero, and it is either a global fast solution or a global slow solution. Finally, we obtain a trichotomy conclusion by considering the size of parameter σ.

A Numerical Method for Multispecies Populations in a Moving Domain Using Combined Masses

This paper concerns the numerical evolution of two interacting species satisfying coupled reaction–diffusion equations in one dimension which inhabit the same part of a moving domain. The domain has both moving external boundaries and moving interior interfaces where species may arise, overlap, or disappear. Numerically, a moving finite volume method is used in which node movement is generated by local mass preservation, which includes a general combined mass strategy for species occupying overlapping domains. The method is illustrated by a test case in which a range of parameters is explored.

A Fisher-KPP Model with a Nonlocal Weighted Free Boundary: Analysis of How Habitat Boundaries Expand, Balance or Shrink

In this paper, we propose a novel free boundary problem to model the movement of single species with a range boundary. The spatial movement and birth/death processes of the species found within the range boundary are assumed to be governed by the classic Fisher-KPP reaction-diffusion equation, while the movement of a free boundary describing the range limit is assumed to be influenced by the weighted total population inside the range boundary and is described by an integro-differential equation. Our free boundary equation is a generalization of the classical Stefan problem that allows for nonlocal influences on the boundary movement so that range expansion and shrinkage are both possible. In this paper, we prove that the new model is well-posed and possesses steady state. We show that the spreading speed of the range boundary is smaller than that for the equivalent problem with a Stefan condition. This implies that the nonlocal effect of the weighted total population on the boundary movement slows down the spreading speed of the population. While the classical Stefan condition categorizes asymptotic behavior via a spreading-vanishing dichotomy, the new model extends this dichotomy to a spreading-balancing-vanishing trichotomy. We specifically analyze how habitat boundaries expand, balance or shrink. When the model is extended to have two free boundaries, we observe the steady state scenario, asymmetric shifts, or even boundaries moving synchronously in the same direction. These are newly discovered phenomena in the free boundary problems for animal movement.

Invasive behaviour under competition via a free boundary model: a numerical approach

What will happen when two invasive species are competing and invading the environment at the same time? In this paper, we try to find all the possible scenarios in such a situation based on the diffusive Lotka-Volterra competition system with free boundaries. In a recent work, Du and Wu (Calc Var Partial Differ Equ, 57(2):52, 2018) considered a weak-strong competition case of this model (with spherical symmetry) and theoretically proved the existence of a "chase-and-run coexistence" phenomenon, for certain parameter ranges when the initial functions are chosen properly. Here we use a numerical approach to extend the theoretical research of Du and Wu (Calc Var Partial Differ Equ, 57(2):52, 2018) in several directions. Firstly, we examine how the longtime dynamics of the model changes as the initial functions are varied, and the simulation results suggest that there are four possible longtime profiles of the dynamics, with the chase-and-run coexistence the only possible profile when both species invade successfully. Secondly, we show through numerical experiments that the basic features of the model appear to be retained when the environment is perturbed by periodic variation in time. Thirdly, our numerical analysis suggests that in two space dimensions the population range and the spatial population distribution of the successful invader tend to become more and more circular as time increases no matter what geometrical shape the initial population range possesses. Our numerical simulations cover the one space dimension case, and two space dimension case with or without spherical symmetry. The numerical methods here are based on that of Liu et al. (Mathematics, 6(5):72, 2018, Int J Comput Math, 97(5): 959-979, 2020). In the two space dimension case without radial symmetry, the level set method is used, while the front tracking method is used for the remaining cases. We hope the numerical observations in this paper can provide further insights to the biological invasion problem, and also to future theoretical investigations. More importantly, we hope the numerical analysis may reach more biologically oriented experts and inspire applications of some refined versions of the model tailored to specific real world biological invasion problems.

The formation of spreading front: the singular limit of three-component reaction-diffusion models

Understanding the invasion processes of biological species is a fundamental issue in ecology. Several mathematical models have been proposed to estimate the spreading speed of species. In recent decades, it was reported that some mathematical models of population dynamics have an explicit form of the evolution equations for the spreading front, which are represented by free boundary problems such as the Stefan-like problem (e.g., Mimura et al., Jpn J Appl Math 2:151-186, 1985; Du and Lin, SIAM J Math Anal 42:377-405, 2010). To understand the formation of the spreading front, in this paper, we will consider the singular limit of three-component reaction-diffusion models and give some interpretations for spreading front from the viewpoint of modeling. As an application, we revisit the issue of the spread of the grey squirrel in the UK and estimate the spreading speed of the grey squirrel based on our result. Also, we discuss the relation between some free boundary problems related to population dynamics and mathematical models describing Controlling Invasive Alien Species. Lastly, we numerically consider the traveling wave solutions, which give information on the spreading behavior of invasive species.

A Wolbachia infection model with free boundary

Scientists have been seeking ways to use Wolbachia to eliminate the mosquitoes that spread human diseases. Could Wolbachia be the determining factor in controlling the mosquito-borne infectious diseases? To answer this question mathematically, we develop a reaction-diffusion model with free boundary in a one-dimensional environment. We divide the female mosquito population into two groups: one is the uninfected mosquito population that grows in the whole region while the other is the mosquito population infected with Wolbachia that occupies a finite small region. The mosquito population infected with Wolbachia invades the environment with a spreading front governed by a free boundary satisfying the well-known one-phase Stefan condition. For the resulting free boundary problem, we establish criteria under which spreading and vanishing occur. Our results provide useful insights on designing a feasible mosquito releasing strategy that infects the whole mosquito population with Wolbachia and eradicates the mosquito-borne diseases eventually.

Pushing the Boundaries: Models for the Spatial Spread of Ecosystem Engineers

Ecosystems engineers are species that can substantially alter their abiotic environment and thereby enhance their population growth. The net growth rate of obligate engineers is even negative unless they modify the environment. We derive and analyze a model for the spread and invasion of such species. Prior to engineering, the landscape consists of unsuitable habitat; after engineering, the habitat is suitable. The boundary between the two types of habitat is moved by the species through their engineering activity. Our model is a novel type of a reaction-diffusion free boundary problem. We prove the existence of traveling waves and give upper and lower bounds for their speeds. We illustrate how the speed depends on individual movement and engineering behavior near the boundary.

Revisiting the Fisher-Kolmogorov-Petrovsky-Piskunov equation to interpret the spreading-extinction dichotomy

The Fisher-Kolmogorov-Petrovsky-Piskunov model, also known as the Fisher-KPP model, supports travelling wave solutions that are successfully used to model numerous invasive phenomena with applications in biology, ecology and combustion theory. However, there are certain phenomena that the Fisher-KPP model cannot replicate, such as the extinction of invasive populations. The Fisher-Stefan model is an adaptation of the Fisher-KPP model to include a moving boundary whose evolution is governed by a Stefan condition. The Fisher-Stefan model also supports travelling wave solutions; however, a key additional feature of the Fisher-Stefan model is that it is able to simulate population extinction, giving rise to a spreading-extinction dichotomy. In this work, we revisit travelling wave solutions of the Fisher-KPP model and show that these results provide new insight into travelling wave solutions of the Fisher-Stefan model and the spreading-extinction dichotomy. Using a combination of phase plane analysis, perturbation analysis and linearization, we establish a concrete relationship between travelling wave solutions of the Fisher-Stefan model and often-neglected travelling wave solutions of the Fisher-KPP model. Furthermore, we give closed-form approximate expressions for the shape of the travelling wave solutions of the Fisher-Stefan model in the limit of slow travelling wave speeds, c≪1.

A diffusive Holling–Tanner prey–predator model with free boundary

This paper is concerned with a diffusive Holling–Tanner prey–predator model in a bounded domain with Dirichlet boundary condition and a free boundary. The global existence of the unique solution is proved. Moreover, the criteria governing spreading–vanishing are derived by mainly using the comparison principle. The results show that if the length of the occupying line is bigger than a threshold value (spreading barrier), then the spreading of predators will make an achievement, and, if the length of the occupying line is smaller than this spreading barrier and the spreading coefficient is relatively small depending on initial size of predators, then the predators will fail in establishing themselves and eventually die out.

One-dimensional collective migration of a proliferating cell monolayer

The importance of collective cellular migration during embryogenesis and tissue repair asks for a sound understanding of underlying principles and mechanisms. Here, we address recent in vitro experiments on cell monolayers, which show that the advancement of the leading edge relies on cell proliferation and protrusive activity at the tissue margin. Within a simple viscoelastic mechanical model amenable to detailed analysis, we identify a key parameter responsible for tissue expansion, and we determine the dependence of the monolayer velocity as a function of measurable rheological parameters. Our results allow us to discuss the effects of pharmacological perturbations on the observed tissue dynamics.