使用平滑粒子流体动力学进行流体模拟

Question

实际上，我正在使用 C++ 和 openGL 开发 SPH 模拟器。我现在遇到了几个问题，

我在盒子的一侧制作引发流体粒子（我将盒子作为边界体积），将它们塑造成盒子的形状，并赋予它们初始速度为零 。然后我开始主循环，中提琴，液体开始移动，这里出现了一个奇怪的现象。液体开始向各个方向扩散。

请看图片：

http://i278.photobucket.com/albums/kk86/ anggytrisnawan/Screenshot-UntitledWindow-2.png

从模拟开始几秒钟后拍摄的图片。对我来说这似乎很奇怪。这是我用于模拟的参数：

#define H                   0.040   // Smoothing Length
#define Rho0            1000    // (kg/m^3) water particle rest density
#define Mass                    0.012       // (kg)
#define DT                  0.001       // time step
#define TotalParticles      5000        // total number of particle

注意：目前我还没有计算表面张力。

已解决：我的错在这里..它们一开始就相互重叠..所以压力使它们散开..

Answer 1

从您提供的图像来看，模拟似乎是从所有粒子都在一个角落开始的。由于粒子可能必须具有一定的排斥势才能相互重叠，因此一旦模拟开始，这种排斥力将迫使粒子分开。

再说一次，您还没有显示任何代码，因此以上只是一个合格的猜测。

Answer 2

边界条件的值是多少，即哪些数字坐标是框限制？我假设，由于系统舍入误差，您的模拟会朝某个方向发展。

如果我看图片，我会觉得我正在沿着 -Z 方向向下看，并且框的左下角位于 (0, 0, 0)。如果是这种情况，那么您的粒子模拟并不保守，即倾向于将粒子推向数值 0，就像存在力场一样。