SPH实践之WCSPH两问题

背景

本人入SPH已有一段时间，近日遇到两个问题（emmm，往日问题无数）。今日心血来潮，写下第一篇博客，既是保存一点工作进展，也供有需要的读者参考。若有错漏，还请指正。

近日实验所使用2D框架是一个改自三维模拟系统，该系统在三维方面表现良好，后强行将z轴的位置和速度设置为0，从而得到二维模拟框架，但二维模拟效果存在一些问题。个人看来，由于模拟系统的复杂性和迭代性，调试难度高。现将近日所遇两问题和解决方法总结如下。

问题一：边界接触处流体粒子躁动

这个问题目前原因未知，由于是边界接触处发生问题，所以自然想到会不会是边界处理引发。

尝试一：背景提到过该框架由三维改造而来，边界用固定粒子模拟，流体部分已经是z值为零的一层粒子，但边界依旧是一个三维盒子，所以本人怀疑是边界接触处的三维边界粒子太多导致的此问题，遂改成二维粒子边界，问题依旧存在。该尝试目前无果。

尝试二：由于本人看过其他几个二维框架，边界均是直接用反弹处理，所以就尝试通过这种方法解决。首先设置反弹界限为x[-2.0, 2.0]，y[0.0, 4.0]。模拟异常，有些粒子会突然获得极大速度导致模拟异常。由于之前调试其他框架时也遇到过相似问现象（当时在其他框架调节其边界后模拟异常），所以将边界调整为x[-1.9, 1.9]，y[0.1, 3.9]后，问题解决。

小结：这个问题目前算是浅度解决，就是模拟异常解决了，但是其中原因没搞清楚。另外说一下，直接反弹边界对比粒子边界