Smoothed Particle Hydrodynamics with Stress Points and Centroid Voronoi Tessellation (CVT) Topology Optimization

Abstract

Various formulations of smoothed particle hydrodynamics (SPH) have been presented by scientists to overcome inherent numerical difficulties including instabilities and inconsistencies. Low approximation accuracy could cause a result of particle inconsistency in SPH and other meshfree methods. In this study, centroid Voronoi tessellation (CVT) topology optimization is used for rearrangement of particles so that the inconsistency due to irregular particle arrangement can be corrected. Using CVT topology optimization method, the SPH particles, which are generated randomly inside a predetermined domain, are moved to the centroids, i.e., the center of mass of the corresponding Voronoi cells based on Lloyd’s algorithm. The volume associated with each particle is determined by its Voronoi cell. On the other hand, it has been shown that particle methods with stress point integration are more stable than the ones using nodal integration. Conventional SPH approximations only use SPH particles, and it results in the so-called tensile instability. In this paper, a new approach of using stress points is introduced to assist SPH approximations and stabilize the SPH methods.