Pbf2d关于h的取值问题

流体模拟小白玩家请教各位:
我跑了一下example里的 pbf2d,然后想增加一些粒子,并且减小粒子半径,我改了两行代码
h = 1.1 / 3
num_particles = num_particles_x * 200
然后粒子就满屏乱飞了…
请教诸位,h的选取有什么限制吗? :face_with_monocle:

我觉得
poly6_value函数和spiky_gradient里,h都是分母的位置
如果分母小于1,计算出来的结果就会非常的大
粒子与粒子之间的作用力也会非常大
然后仿真就爆炸了