On meshless methods : a novel interpolatory method and a GPU-accelerated implementation

摘要

Meshless methods have been developed to avoid the numerical burden imposed by meshing in the Finite Element Method. Such methods are especially attrac- tive in problems that require repeated updates to the mesh, such as problems with discontinuities or large geometrical deformations. Although meshing is not required for solving problems with meshless methods, the use of meshless methods gives rise to different challenges. One of the main challenges associated with meshless methods is imposition of essential boundary conditions. If exact interpolants are used as shape functions in a meshless method, imposing essen- tial boundary conditions can be done in the same way as the Finite Element Method. Another attractive feature of meshless methods is that their use involves compu- tations that are largely independent from one another. This makes them suitable for implementation to run on highly parallel computing systems. Highly par- allel computing has become widely available with the introduction of software development tools that enable developing general-purpose programs that run on Graphics Processing Units. In the current work, the Moving Regularized Interpolation method has been de- veloped, which is a novel method of constructing meshless shape functions that achieve exact interpolation. The method is demonstrated in data interpolation and in partial differential equations. In addition, an implementation of the Element-Free Galerkin method has been written to run on a Graphics Processing Unit. The implementation is described and its performance is compared to that of a similar implementation that does not make use of the Graphics Processing Unit.