Numerical Methods for a Two-Species Competition-Diffusion Model with Free Boundaries

The systems of reaction-diffusion equations coupled with moving boundaries defined by Stefan condition have been widely used to describe the dynamics of spreading population and with competition of two species. To solve these systems numerically, new numerical challenges arise from the competition of two species due to the interaction of their free boundaries. On the one hand, extremely small time steps are usually needed due to the stiffness of the system. On the other hand, it is always difficult to efficiently and accurately handle the moving boundaries especially with competition of two species. To overcome these numerical difficulties, we introduce a front tracking method coupled with an implicit solver for the 1D model. For the general 2D model, we use a level set approach to handle the moving boundaries to efficiently treat complicated topological changes. Several numerical examples are examined to illustrate the efficiency, accuracy and consistency for different approaches.

[1]  Shuang Liu.  Numerical Methods for a Class of Reaction-Diffusion Equations With Free Boundaries , 2019 .

[2]  Jong-Shenq Guo,et al.  On a Free Boundary Problem for a Two-Species Weak Competition System , 2012 .

[3]  Kelei Wang,et al.  Regularity and Asymptotic Behavior of Nonlinear Stefan Problems , 2014 .

[4]  Ning Zhao,et al.  Conservative front tracking and level set algorithms , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[5]  J. Sethian,et al.  FRONTS PROPAGATING WITH CURVATURE DEPENDENT SPEED: ALGORITHMS BASED ON HAMILTON-JACOB1 FORMULATIONS , 2003 .

[6]  Phillip Colella,et al.  A Front Tracking Method for Compressible Flames in One Dimension , 1995, SIAM J. Sci. Comput..

[7]  C. Peskin.  The immersed boundary method , 2002, Acta Numerica.

[8]  Bingtuan Li,et al.  Cooperative , 1952, The Fairchild Books Dictionary of Fashion.

[9]  Bingtuan Li,et al.  Spreading speed and linear determinacy for two-species competition models , 2002, Journal of mathematical biology.

[10]  J A Sethian,et al.  A fast marching level set method for monotonically advancing fronts. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[11]  Robert A. Gardner.  Existence and stability of travelling wave solutions of competition models: A degree theoretic approach , 1982 .

[12]  Yihong Du,et al.  The Stefan problem for the Fisher–KPP equation with unbounded initial range , 2012, Calculus of Variations and Partial Differential Equations.

[13]  Ping Lin,et al.  Numerical simulation of 3D bubbles rising in viscous liquids using a front tracking method , 2008, J. Comput. Phys..

[14]  Yihong Du,et al.  Nonlinear Diffusion Problems with Free Boundaries: Convergence, Transition Speed, and Zero Number Arguments , 2015, SIAM J. Math. Anal..

[15]  D. Hilhorst,et al.  A competition-diffusion system approximation to the classical two-phase Stefan problem , 2001 .

[16]  Note on a two-species competition-diffusion model with two free boundaries , 2015, 1508.07545.

[17]  Randall J. LeVeque,et al.  Two-Dimensional Front Tracking Based on High Resolution Wave Propagation Methods , 1996 .

[18]  G. Tryggvason,et al.  A front-tracking method for viscous, incompressible, multi-fluid flows , 1992 .

[19]  C. Peskin,et al.  Simulation of a Flapping Flexible Filament in a Flowing Soap Film by the Immersed Boundary Method , 2002 .

[20]  Y. Hosono,et al.  The minimal speed of traveling fronts for a diffusive Lotka-Volterra competition model , 1998 .

[21]  C. Conley,et al.  An application of the generalized Morse index to travelling wave solutions of a competitive reaction-diffusion model , 1984 .

[22]  Chang-Hong Wu,et al.  Dynamics for a two-species competition–diffusion model with two free boundaries , 2014 .

[23]  H. G. Landau,et al.  Heat conduction in a melting solid , 1950 .

[24]  Zhilin Li,et al.  A level-set method for interfacial flows with surfactant , 2006, J. Comput. Phys..

[25]  K. Bube,et al.  The Immersed Interface Method for Nonlinear Differential Equations with Discontinuous Coefficients and Singular Sources , 1998 .

[26]  Chang-Hong Wu.  The minimal habitat size for spreading in a weak competition system with two free boundaries , 2015 .

[27]  R. LeVeque,et al.  A comparison of the extended finite element method with the immersed interface method for elliptic equations with discontinuous coefficients and singular sources , 2006 .

[28]  Yihong Du,et al.  Spreading and vanishing in nonlinear diffusion problems with free boundaries , 2013, 1301.5373.

[29]  Ronald Fedkiw,et al.  Level set methods and dynamic implicit surfaces , 2002, Applied mathematical sciences.

[30]  T. R. Peters,et al.  AP-3 Directs the Intracellular Trafficking of HIV-1 Gag and Plays a Key Role in Particle Assembly , 2005, Cell.

[31]  Benjamin Pfaff,et al.  Free And Moving Boundary Problems , 2016 .

[32]  T. Chan,et al.  A Variational Level Set Approach to Multiphase Motion , 1996 .

[33]  Bingtuan Li,et al.  Analysis of linear determinacy for spread in cooperative models , 2002, Journal of mathematical biology.

[34]  Oliver A. McBryan,et al.  Front Tracking for Gas Dynamics , 1984 .

[35]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

[36]  Mingxin Wang,et al.  Free Boundary Problems for a Lotka–Volterra Competition System , 2014, 1401.0806.

[37]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods , 1999 .