A Special Sorting Method for Neighbor Search Procedure in Smoothed Particle Hydrodynamics on GPUs

Author(s):  
Hongyu Sun ◽  
Yanshan Tian ◽  
Yulong Zhang ◽  
Jiong Wu ◽  
Sen Wang ◽  
...  
2019 ◽  
Vol 36 (3) ◽  
pp. 899-929 ◽  
Author(s):  
Lin Fu ◽  
Zhe Ji ◽  
Xiangyu Y. Hu ◽  
Nikolaus A. Adams

Purpose This paper aims to develop a parallel fast neighbor search method and communication strategy for particle-based methods with adaptive smoothing-length on distributed-memory computing systems. Design/methodology/approach With a multi-resolution-based hierarchical data structure, the parallel neighbor search method is developed to detect and construct ghost buffer particles, i.e. neighboring particles on remote processor nodes. To migrate ghost buffer particles among processor nodes, an undirected graph is established to characterize the sparse data communication relation and is dynamically recomposed. By the introduction of an edge coloring algorithm from graph theory, the complex sparse data exchange can be accomplished within optimized frequency. For each communication substep, only efficient nonblocking point-to-point communication is involved. Findings Two demonstration scenarios are considered: fluid dynamics based on smoothed-particle hydrodynamics with adaptive smoothing-length and a recently proposed physics-motivated partitioning method [Fu et al., JCP 341 (2017): 447-473]. Several new concepts are introduced to recast the partitioning method into a parallel version. A set of numerical experiments is conducted to demonstrate the performance and potential of the proposed parallel algorithms. Originality/value The proposed methods are simple to implement in large-scale parallel environment and can handle particle simulations with arbitrarily varying smoothing-lengths. The implemented smoothed-particle hydrodynamics solver has good parallel performance, suggesting the potential for other scientific applications.


2008 ◽  
Vol 96 (6) ◽  
pp. 263-268 ◽  
Author(s):  
E. Mounif ◽  
V. Bellenger ◽  
A. Ammar ◽  
R. Ata ◽  
P. Mazabraud ◽  
...  

2020 ◽  
Vol 59 (40) ◽  
pp. 18236-18246
Author(s):  
Tianwen Dong ◽  
Yadong He ◽  
Jianchun Wu ◽  
Shiyu Jiang ◽  
Xingyuan Huang ◽  
...  

Author(s):  
Steven J. Lind ◽  
Benedict D. Rogers ◽  
Peter K. Stansby

This paper presents a review of the progress of smoothed particle hydrodynamics (SPH) towards high-order converged simulations. As a mesh-free Lagrangian method suitable for complex flows with interfaces and multiple phases, SPH has developed considerably in the past decade. While original applications were in astrophysics, early engineering applications showed the versatility and robustness of the method without emphasis on accuracy and convergence. The early method was of weakly compressible form resulting in noisy pressures due to spurious pressure waves. This was effectively removed in the incompressible (divergence-free) form which followed; since then the weakly compressible form has been advanced, reducing pressure noise. Now numerical convergence studies are standard. While the method is computationally demanding on conventional processors, it is well suited to parallel processing on massively parallel computing and graphics processing units. Applications are diverse and encompass wave–structure interaction, geophysical flows due to landslides, nuclear sludge flows, welding, gearbox flows and many others. In the state of the art, convergence is typically between the first- and second-order theoretical limits. Recent advances are improving convergence to fourth order (and higher) and these will also be outlined. This can be necessary to resolve multi-scale aspects of turbulent flow.


2013 ◽  
Vol 61 (1) ◽  
pp. 111-121 ◽  
Author(s):  
T. Jankowiak ◽  
T. Łodygowski

Abstract The paper considers the failure study of concrete structures loaded by the pressure wave due to detonation of an explosive material. In the paper two numerical methods are used and their efficiency and accuracy are compared. There are the Smoothed Particle Hydrodynamics (SPH) and the Finite Element Method (FEM). The numerical examples take into account the dynamic behaviour of concrete slab or a structure composed of two concrete slabs subjected to the blast impact coming from one side. The influence of reinforcement in the slab (1, 2 or 3 layers) is also presented and compared with a pure concrete one. The influence of mesh density for FEM and the influence of important parameters in SPH like a smoothing length or a particle distance on the quality of the results are discussed in the paper


Sign in / Sign up

Export Citation Format

Share Document