我正在尝试实现用于在 JavaScript 中模拟流体的Clavet 方法，所以调试是一场噩梦，这就是我在这里问的原因，希望经历过同样事情的人能告诉我我做错了什么。
到目前为止，我的工作正常：
但是我有两个问题：
1）由于在这种方法中，所有东西都有点“偏移”半步，我不确定如何正确地将粒子从墙上反弹。现在，我获取粒子的位置和之前的位置，并将它们围绕交叉的墙壁翻转，然后通过反弹因子围绕交叉点进行缩放。
我的逻辑告诉我这应该有效。算法的下一步是更新粒子速度，所以我也反映了之前的位置。但在实践中，这给了我一个我不明白的结果：
这显示了粒子上的“力”。墙壁反射了太多的力量，这让一切都保持在永恒的运动中。
本文中的公式 4.58 显然显示了一种防止这种情况的方法，但我无法让它发挥作用。
论文中还有一些我没有得到的东西，比如“我们只想反映碰撞中省略的速度”的意思。为什么？有人可以把这些东西给我吗？
2）即使不涉及墙壁，模拟也会周期性地“爆炸”。这种情况在更高的压力下发生得更多：
这是JavaScript，所以有，但我已经检查了代码，没有除以零或我可以想象发生 NaN 的情况。
我在报纸上看到过一些关于 sim 不稳定的讨论，我想知道是不是这样。这部文学作品中的大部分内容都超出了我的理解范围。
据我所知（我认为），消除不稳定性的方法之一是粘度，但我添加了它，但它对爆炸没有帮助：
我可以发布代码，但在开始工作的第一阶段，现在有点难以阅读。
最后一个问题：我如何弄清楚如何从这种方法中的伪常数转换为物理单位？
编辑：我发现 sim 偶尔会冻结，似乎它确实在某处产生了 NaN，但 Chrome 捕捉到它为时已晚。