Graphics Processing Unit-Accelerated Boundary Element Method and Vortex Particle Method

2011 ◽  
Vol 8 (7) ◽  
pp. 224-236 ◽  
Author(s):  
Mark J. Stock ◽  
Adrin Gharakhani
Author(s):  
Yingjun Wang ◽  
Qifu Wang ◽  
Gang Wang ◽  
Yunbao Huang ◽  
Yixiong Wei

Finite Element Method (FEM) is pervasively used in most of 3D elastostatic numerical simulations, in which Computer Aided Design (CAD) models need to be converted into mesh models first and then enriched with semantic data (e.g. material parameters, boundary conditions). The interaction between CAD models and FEM models stated above is very intensive. Boundary Element Method (BEM) has been used gradually instead of FEM in recent years because of its advantage in meshing. BEM can reduce the dimensionality of the problem by one so that the complexity in mesh generation can be decreased greatly. In this paper, we present a Boundary Element parallel computation method for 3D elastostatics. The parallel computation runs on Graphics Processing Unit (GPU) using Computing Unified Device Architecture (CUDA). Three major components are included in such method: (1) BEM theory in 3D elastostatics and the boundary element coefficient integral methods, (2) the parallel BEM algorithm using CUDA, and (3) comparison the parallel BEM using CUDA with conventional BEM and FEM respectively by examples. The dimension reduction characteristics of BEM can dispose the 3D elastostatic problem by 2D meshes, therefore we develop a new faceting function to make the ACIS facet meshes suitable for Boundary Element Analysis (BEA). The examples show that the GPU parallel algorithm in this paper can accelerate BEM computation about 40 times.


2020 ◽  
Vol 7 (2) ◽  
pp. 211-227
Author(s):  
Ahmed A Torky ◽  
Youssef F Rashed

Abstract This study develops a high-performance computing method using OpenACC (Open Accelerator) for the stiffness matrix and load vector generation of shear-deformable plates in bending using the boundary element method on parallel processors. The boundary element formulation for plates in bending is used to derive fully populated displacement-based stiffness matrices and load vectors at degrees of freedom of interest. The computed stiffness matrix of the plate is defined as a single superfloor element and can be solved using stiffness analysis, $Ku = F$, instead of the conventional boundary element method, $Hu = Gt$. Fortran OpenACC code implementations are proposed for the computation of the superfloor element’s stiffness, which includes one serial computing code for the CPU (central processing unit) and two parallel computing codes for the GPU (graphics processing unit) and multicore CPU. As industrial level practical floors are full of supports and geometrical information, the computation time of superfloor elements is reduced dramatically when computing on parallel processors. It is demonstrated that the OpenACC implementation does not affect numerical accuracy. The feasibility and accuracy are confirmed by numerical examples that include real buildings with industrial level structural floors. Engineering computations for massive floors with immense geometrical detail and a multitude of load cases can be modeled as is without the need for simplification.


2013 ◽  
Vol 49 (5) ◽  
pp. 1737-1740 ◽  
Author(s):  
Gaku Yoshikawa ◽  
Fumikazu Miyasaka ◽  
Katsuhiro Hirata ◽  
Shuhei Matsuzawa

2014 ◽  
Vol 12 (2) ◽  
Author(s):  
Kirana Kumara P

In this work, possibility of simulating biological organs in realtime using the Boundary Element Method (BEM) is investigated. Biological organs are assumed to follow linear elastostatic material behavior, and constant boundary element is the element type used.  First, a Graphics Processing Unit (GPU) is used to speed up the BEM computations to achieve the realtime performance. Next, instead of the GPU, a computer cluster is used.  Results indicate that BEM is fast enough to provide for realtime graphics if biological organs are assumed to follow linear elastostatic material behavior. Although the present work does not conduct any simulation using nonlinear material models, results from using the linear elastostatic material model imply that it would be difficult to obtain realtime performance if highly nonlinear material models that properly characterize biological organs are used. Although the use of BEM for the simulation of biological organs is not new, the results presented in the present study are not found elsewhere in the literature.


1995 ◽  
Vol 117 (3) ◽  
pp. 468-472
Author(s):  
D. E. Papantonis ◽  
K. P. Pothou

In the present work, the numerical prediction of the steady as well as of the unsteady flowfield through the cascade of guidevanes of a reaction turbine is examined. The unsteadiness of the flow results from the rotation of the guidevanes around their pivot. The applied numerical procedure is based on the coupling of the boundary element and the vortex particle method and is supplemented by two different procedures for the evaluation of the viscous-inviscid interaction and the boundary layers developed on the vanes. The cinematic conditions at every point of the flowfield are calculated first; however, special attention is paid to the prediction of the torque acting on the cascade and, namely, the variation of the torque during the rotation of the guidevanes, as compared with the torque corresponding to the steady state.


Sign in / Sign up

Export Citation Format

Share Document