改进的接触算法及其在光滑粒子流体动力学中的应用

摘要

采用以黎曼解描述粒子间相互作用的接触算法,模拟一维情况下包含间断的流场.在弱间断的条件下,采用弱波近似的黎曼解来描述粒子间的相互作用;在强间断的条件下,则引入非迭代黎曼解来描述粒子间的相互作用.并引入Taylor展开思想,提高了接触算法在自由边界的计算精度.该算法无需引入人为粘性,程序结构简洁.数值计算表明,改进后的接触算法能很好地描述强间断,并有效改善了原始接触算法在自由面的计算缺陷.