Khái niệm cốt lõi
The fast-convolving reproducing kernel particle method (FC-RKPM) is introduced, which is hundreds to millions of times faster than the traditional RKPM for 3D meshfree simulations by expressing the meshfree discretizations with RK approximation in terms of convolution sums and using fast Fourier transform (FFT) to efficiently compute the convolutions.
Tóm tắt
The paper introduces the fast-convolving reproducing kernel particle method (FC-RKPM), which is significantly more efficient than the traditional RKPM for 3D meshfree simulations. The key ideas are:
- The meshfree discretizations with RK approximation are expressed in terms of convolution sums.
- Fast Fourier transform (FFT) is then used to efficiently compute the convolutions.
- Certain modifications to the domain and shape functions are considered to maintain generality for complex geometries and arbitrary boundary conditions.
- The new method does not need to identify, store, and loop over the neighbors, which is a major bottleneck of traditional meshfree methods. As a result, the run-times and memory allocations are independent of the number of neighbors and the shape function's support size.
- The method is verified for a Galerkin weak form of the Poisson problem with the RK approximation in 1D, 2D, and 3D. Tables with run-times and allocated memory are presented to compare the performance of FC-RKPM with the traditional method in 3D.
- The performance is studied for various node numbers, support size, and approximation degree.
- Implementation details and a roadmap for software development are provided.
- Application of the new method to nonlinear and explicit problems are briefly discussed.
Thống kê
Tables with run-times and allocated memory are presented to compare the performance of FC-RKPM with the traditional method in 3D.
Trích dẫn
"The new method does not need to identify, store, and loop over the neighbors which is one of the bottleneck of the traditional meshfree methods."
"As a result, the run-times and memory allocations are independent of the number of neighbors and the shape function's support size."