, The use of element free Galerkin method based on moving Kriging and radial point interpolation techniques for solving some types of Turing models

نویسندگانمهدی دهقان-مصطفی عباس زاده-اکبر محبی
نشریهENG ANAL BOUND ELEM
تاریخ انتشار۲۰۱۵-۱۰-۰۱
نوع نشریهالکترونیکی
نمایه نشریهISI

چکیده مقاله

In this paper two numerical procedures are presented for solving a class of Turing system. Firstly, we obtain a time discrete scheme by approximating time derivative via fi nite difference technique. Then we introduce the moving Kriging interpolation and radial point interpolation and also obtain their shape functions. We use the element free Galerkin method for approximating the spatial derivatives. This method uses a weak form of the considered equation that is similar to the fi nite element method with the difference that in the classical element free Galerkin method test and trial functions are moving least squares (ML S) approximation shape functions. Since the shape functions of moving least squares (ML S) approximation do not have Kronecker delta property, we cannot implement the essential boundary condition, directly. Thus we employ the shape functions of moving Kriging interpolation and radial point interpolation technique which have the mentioned property. Also, in the element free Galerkin method, we do not use any triangular, quadrangular or other type of meshes. The element free Galerkin method is a global method while fi nite elements method is a local one. This technique employs a background mesh for integration which makes it different from the truly mesh procedures. The coeffi cient matrix of the element free Galerkin is symmetric. Also, using numerical algorithms, we can conclude that the eigen- values of the coeffi cient matrix are positive. Thus, for solving the obtained linear system of equations from the discretization, we use the conjugant gradient method. To keep away from solving a nonlinear algebraic system of equations and obtaining the acceptable numerical results, we use a predictor –cor- rector algorithm. Several test problems are solved and numerical simulations are reported which confi rm the ef fi ciency of the proposed schemes.