Numerical Approach to a Nonlocal Advection-Reaction-Diffusion Model of Cartilage Pattern Formation

Abstract

We propose a numerical approach that combines a radial basis function (RBF) meshless approximation with a finite difference discretization to solve a nonlinear system of integro-differential equations. The equations are of advection-reaction-diffusion type modeling the formation of pre-cartilage condensations in embryonic chicken limbs. The computational domain is four dimensional in the sense that the cell density depends continuously on two spatial variables as well as two structure variables, namely membrane-bound counterreceptor densities. The biologically proper Dirichlet boundary conditions imposed in the semi-infinite structure variable region is in favor of a meshless method with Gaussian basis functions. Coupled with WENO5 finite difference spatial discretization and the method of integrating factors, the time integration via method of lines achieves optimal complexity. In addition, the proposed scheme can be extended to similar models with more general boundary conditions. Numerical results are provided to showcase the validity of the scheme.

Keywords

• advection-reaction-diffusion systems
• radial basis functions
• cell-cell adhesion
• pattern formation
• models of chondrogenesis