Survival, extinction, and interface stability in a two--phase moving boundary model of biological invasion

We consider a moving boundary mathematical model of biological invasion. The model describes the spatiotemporal evolution of two adjacent populations: each population undergoes linear diffusion and logistic growth, and the boundary between the two populations evolves according to a two--phase Stefan condition. This mathematical model describes situations where one population invades into regions occupied by the other population, such as the spreading of a malignant tumour into surrounding tissues. Full time--dependent numerical solutions are obtained using a level--set numerical method. We use these numerical solutions to explore several properties of the model including: (i) survival and extinction of one population initially surrounded by the other; and (ii) linear stability of the moving front boundary in the context of a travelling wave solution subjected to transverse perturbations. Overall, we show that many features of the well--studied one--phase single population analogue of this model can be very different in the more realistic two--phase setting. These results are important because realistic examples of biological invasion involve interactions between multiple populations and so great care should be taken when extrapolating predictions from a one--phase single population model to cases for which multiple populations are present. Open source Julia--based software is available on GitHub to replicate all results in this study.