The element-free Galerkin method based on moving least squares and moving Kriging approximations for solving two-dimensional tumor-induced angiogenesis model

The numerical simulation of the tumor-induced angiogenesis process is an useful tool for the prediction of this mechanism and drug targeting using anti-angiogenesis strategy. In the current paper, we study numerically on the continuous mathematical model of tumor-induced angiogenesis in two-dimensional spaces. The studied model is a system of nonlinear time-dependent partial differential equations, which describes the interactions between endothelial cell, tumor angiogenesis factor and fibronectin. We first derive the global weak form of the model and discretize the time variable via a semi-implicit backward Euler method. To approximate the spatial variables of the studied model, we use a meshless technique, namely element-free Galerkin. Also, the shape functions of moving least square and moving Kriging approximations are used in this method. The main difference between two meshless methods proposed here is that the shape functions of moving least squares approximation do not satisfy Kroncker's delta property, while moving Kriging technique satisfies this property. Also, both techniques do not require the generation of a mesh for approximation, but a background mesh is needed to compute the numerical integrations, which are appeared in the derived global weak form. The full-discrete scheme obtained here gives the linear system of algebraic equations that is solved via an iterative method, namely biconjugate gradient stabilized with zero-fill incomplete lower upper (ILU) preconditioner. Some numerical simulations are provided to illustrate the ability of the presented numerical methods, which show the endothelial cell migration in response to the tumor angiogenesis factors during angiogenesis process as well.