The post-processing method of SPH on non-convex configuration

A post-processing method for SPH was presented based on Delaunay triangulation. However for non-convex domain occupied by SPH particles(nodes), Delaunay triangulation results in some empty elements which contains no mass. Thus an "element mass weighing" method should be proposed to identify and delete those empty elements from the element collection. The reasonability of the algorithm is that each empty element has a very much less SPH-weighted-approximation than the average quality of the particles on it. With the SPH data and the kept elements being exported into the universal post-processing software Tecplot, the values at any point in the computational domain could be obtained. The feasibility of the algorithm was tested in a series of numerical cases including a case with heterogeneous particle distribution and material splashing.