Prevention of numerical fracture in smoothed particle hydrodynamics method

To solve the problem of numerical fracture in the process of numerical simulation of smoothed particle hydrodynamics method, based on Delaunay triangulation and Voronoi diagram, a new technique of particle generation and particle elimination was presented to redistribute the scattered particles of the numerical fracture region. And meanwhile, the conservation of mass,momentum and energy were discussed. Two typical numerical examples were simulated and the results show that the problem of numerical fracture is solved effectively. When investigate the change of relative runtime along with the particle number of the model, if only use particle generation technique, due to the increase of particle number, the runtime go up sharply. But if both of the particle generation and particle elimination techniques are used, because of the invariability of particle number, the runtime increases by only about 10%. This percentage is acceptable obviously. When investigate the change of total energy all through the process, results show that the total energy is almost the same from beginning to the end of the whole process.