Volume 19, Issue 3
Efficient Implementation of Smoothed Particle Hydrodynamics (SPH) with Plane Sweep Algorithm

Dong Wang, Yisong Zhou & Sihong Shao

Commun. Comput. Phys., 19 (2016), pp. 770-800.

Published online: 2018-04

Preview Purchase PDF 127 1097
Export citation
  • Abstract

Neighbour search (NS) is the core of any implementations of smoothed particle hydrodynamics (SPH). In this paper, we present an efficient O(NlogN) neighbour search method based on the plane sweep (PW) algorithm with N being the number of SPH particles. The resulting method, dubbed the PWNS method, is totally independent of grids (i.e., purely meshfree) and capable of treating variable smoothing length, arbitrary particle distribution and heterogenous kernels. Several state-of-the-art data structures and algorithms, e.g., the segment tree and the Morton code, are optimized and implemented. By simply allowing multiple lines to sweep the SPH particles simultaneously from different initial positions, a parallelization of the PWNS method with satisfactory speedup and load-balancing can be easily achieved. That is, the PWNS SPH solver has a great potential for large scale fluid dynamics simulations.

  • Keywords

  • AMS Subject Headings

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address
  • BibTex
  • RIS
  • TXT
@Article{CiCP-19-770, author = {}, title = {Efficient Implementation of Smoothed Particle Hydrodynamics (SPH) with Plane Sweep Algorithm}, journal = {Communications in Computational Physics}, year = {2018}, volume = {19}, number = {3}, pages = {770--800}, abstract = {

Neighbour search (NS) is the core of any implementations of smoothed particle hydrodynamics (SPH). In this paper, we present an efficient O(NlogN) neighbour search method based on the plane sweep (PW) algorithm with N being the number of SPH particles. The resulting method, dubbed the PWNS method, is totally independent of grids (i.e., purely meshfree) and capable of treating variable smoothing length, arbitrary particle distribution and heterogenous kernels. Several state-of-the-art data structures and algorithms, e.g., the segment tree and the Morton code, are optimized and implemented. By simply allowing multiple lines to sweep the SPH particles simultaneously from different initial positions, a parallelization of the PWNS method with satisfactory speedup and load-balancing can be easily achieved. That is, the PWNS SPH solver has a great potential for large scale fluid dynamics simulations.

}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.010415.110915a}, url = {http://global-sci.org/intro/article_detail/cicp/11109.html} }
TY - JOUR T1 - Efficient Implementation of Smoothed Particle Hydrodynamics (SPH) with Plane Sweep Algorithm JO - Communications in Computational Physics VL - 3 SP - 770 EP - 800 PY - 2018 DA - 2018/04 SN - 19 DO - http://dor.org/10.4208/cicp.010415.110915a UR - https://global-sci.org/intro/article_detail/cicp/11109.html KW - AB -

Neighbour search (NS) is the core of any implementations of smoothed particle hydrodynamics (SPH). In this paper, we present an efficient O(NlogN) neighbour search method based on the plane sweep (PW) algorithm with N being the number of SPH particles. The resulting method, dubbed the PWNS method, is totally independent of grids (i.e., purely meshfree) and capable of treating variable smoothing length, arbitrary particle distribution and heterogenous kernels. Several state-of-the-art data structures and algorithms, e.g., the segment tree and the Morton code, are optimized and implemented. By simply allowing multiple lines to sweep the SPH particles simultaneously from different initial positions, a parallelization of the PWNS method with satisfactory speedup and load-balancing can be easily achieved. That is, the PWNS SPH solver has a great potential for large scale fluid dynamics simulations.

Dong Wang, Yisong Zhou & Sihong Shao. (2020). Efficient Implementation of Smoothed Particle Hydrodynamics (SPH) with Plane Sweep Algorithm. Communications in Computational Physics. 19 (3). 770-800. doi:10.4208/cicp.010415.110915a
Copy to clipboard
The citation has been copied to your clipboard