基于黎曼解的一种SPH改进算法

A Modified SPH Method Based on Riemann Solution

将基于黎曼解的粒子间接触算法与基于Taylor展开的插值方法相结合,构造出一种改进的光滑粒子算法,给出了两个算例.数值计算表明,改进算法可以提高自由面计算精度、克服拉应力不稳定性,避免了CSPM在核函数的影响域包含较多粒子且流场梯度很大时出现的压缩失稳现象.

An interparticle contact algorithm based on Riemann solvers and an interpolation method in which the Taylor series incorporates with a kernel estimation are combined to construct a modified smoothed particle hydrodynamic algorithm. Numerical simulations show that the modified algorithm ensures numerical solution accuracy on the boundary and eliminates tensile instability occurred at contact interfaces in impact problemes. It removes numerical compression instability which appears in corrective smoothed particle method （CSPM） as the domain of influence of a kernel function includes many particles or the gradient of the flow is sharp.