复现iisph遇到问题

运行环境

taichi 0.8.9

问题描述

我在课程提供的sph_base和particle_system基础上,复现iisph原论文时,发现无法收敛压强,具体表现为a_ii的值很小,迭代过程中压强p_i,l+1 = (density_0 - desity_i,adv - (…))/a_ii非常大,超过正常的压强应有的数值
我读公式没发现问题,求大佬们教教

代码链接

链接
论文pdf

1 Like

Hi @emt, 你好,非常欢迎你来到太极论坛。

你提供的代码,只有一个class,能够提供一个能跑的代码么?还有就是SPH base 类也可以放在你的repo里面。

不好意思,我之前疏忽忘记放了 ……
已经更新在repo里了,运行ii_demo.py会调用iisph.py和其他文件,麻烦大佬看一下谢谢

1 Like