A parallel numerical acoustic simulation on a GPU using an edge-based smoothed finite element method

Abstract In this paper, a parallel computing scheme for performing implicit finite element calculations on acoustic problems running on a graphics processing unit (GPU) is proposed. This scheme combines the element-by-element (EBE) strategy and the edge-based smoothed finite element method (ES-FEM). A mapping strategy based on the ES-FEM is developed to parallel calculate various formulas for smoothing domains. A semi-parallel construction method for the smoothing domain is developed to address the fact that the smoothing domain is difficult to construct in parallel computing. This construction method is further developed into a pre-indexing method of the degrees of freedom (DOFs), which can be used with the parallel preconditioned conjugate gradient (PCG) method of solving linear equations to improve its computational efficiency. After the program flow diagram has been designed, Nvidia's CUDA is used to implement GPU-based parallel computing. Several numerical experiments are used to demonstrate the accuracy and efficiency of the proposed scheme.

[1]  S. Gyimothy,et al.  Parallel Realization of the Element-by-Element FEM Technique by CUDA , 2012, IEEE Transactions on Magnetics.

[2]  Martin Ochmann,et al.  Boundary Element Acoustics Fundamentals and Computer Codes , 2002 .

[3]  K. Y. Dai,et al.  A Smoothed Finite Element Method for Mechanics Problems , 2007 .

[4]  E. Li,et al.  Stability and accuracy improvement for explicit formulation of time domain acoustic problems , 2017 .

[5]  P. Bouillard,et al.  Element-Free Galerkin solutions for Helmholtz problems: fomulation and numerical assessment of the pollution effect , 1998 .

[6]  Bing Li,et al.  3D mass-redistributed finite element method in structural–acoustic interaction problems , 2016 .

[7]  Qunli Sun,et al.  Acoustic-structural coupled finite element analysis for sound transmission in human ear--pressure distributions. , 2006, Medical engineering & physics.

[8]  K. Law.  A parallel finite element solution method , 1986 .

[9]  Xiangyang Cui,et al.  Numerical treatment of acoustic problems with the smoothed finite element method , 2010 .

[10]  Genki Yagawa,et al.  Parallel computing of high‐speed compressible flows using a node‐based finite‐element method , 2003 .

[11]  Jérémie Allard,et al.  Implicit FEM Solver on GPU for Interactive Deformation Simulation , 2011 .

[12]  Genki Yagawa,et al.  Node‐by‐node parallel finite elements: a virtually meshless method , 2004 .

[13]  K. Y. Dai,et al.  Theoretical aspects of the smoothed finite element method (SFEM) , 2007 .

[14]  F. Ihlenburg.  Finite Element Analysis of Acoustic Scattering , 1998 .

[15]  José Escolano,et al.  Parallelization of the finite-difference time-domain method for room acoustics modelling based on CUDA , 2013, Math. Comput. Model..

[16]  T. Hughes,et al.  An element-by-element solution algorithm for problems of structural and solid mechanics , 1983 .

[17]  Guirong Liu,et al.  Hybrid smoothed finite element method for acoustic problems , 2015 .

[18]  Guiyong Zhang,et al.  Analysis of elastic-plastic problems using edge-based smoothed finite element method , 2009 .

[19]  Guirong Liu,et al.  An edge-based smoothed finite element method (ES-FEM) for analyzing three-dimensional acoustic problems , 2009 .

[20]  E. Li,et al.  Optimal balance between mass and smoothed stiffness in simulation of acoustic problems , 2019, Applied Mathematical Modelling.

[21]  Arne S. Gullerud,et al.  MPI-based implementation of a PCG solver using an EBE architecture and preconditioner for implicit, 3-D finite element analysis , 2001 .

[22]  Wu-chun Feng,et al.  Performance Characterization and Optimization of Atomic Operations on AMD GPUs , 2011, 2011 IEEE International Conference on Cluster Computing.

[23]  Guangyao Li,et al.  A parallel finite element procedure for contact-impact problems using edge-based smooth triangular element and GPU , 2018, Comput. Phys. Commun..

[24]  Guangyao Li,et al.  Development of parallel explicit finite element sheet forming simulation system based on GPU architecture , 2012, Adv. Eng. Softw..

[25]  Renyuan Tang,et al.  GPU acceleration of EBE method for 3-D linear steady eddy current field , 2015, 2015 18th International Conference on Electrical Machines and Systems (ICEMS).

[26]  J. Shewchuk.  An Introduction to the Conjugate Gradient Method Without the Agonizing Pain , 1994 .

[27]  William J. Dally,et al.  Scatter-add in data parallel architectures , 2005, 11th International Symposium on High-Performance Computer Architecture.

[28]  T. Belytschko,et al.  Element‐free Galerkin methods , 1994 .

[29]  Jesús Martínez-Frutos,et al.  Efficient matrix-free GPU implementation of Fixed Grid Finite Element Analysis , 2015 .