@李东岳 好的，谢谢李老师 ：-）

@youhaoyu 1719302824186.jpg 看了本比较通用的书，里面把我第一张图里（1.15）直接写在方程里面的（我感觉这样写比较合适）。那么这个R0和Pg0我还是不理解什么意思，有没有大佬解释解释

opm_set_num_threads应该是设置openmp并行线程数量的，具体错误原因还需要看详细代码

@AppleKiller 感谢，我再看看，以前也接触过这些资料，感觉用欧拉法求解器为基础修改的拉格朗日求解器不大好，没怎么深入去学习。我再看看这些思路，感谢

@李东岳 东岳老师您好，我最近发现了一个现象可能造成了我上述提出的问题哈哈哈哈哈，我最开始使用的是DPMFoam默认的的Ergunwenyuforce曳力，我后来换成了sphereforce曳力就能稳定了，而且对于颗粒运动和流场有挺大影响，看来这和我的设置也有一些关系。（最开始图方便就直接使用的默认的模板了，看来还是要根据实际情况来改变）

哪位老师知道目前哪一款CFD软件能够精确模拟这种现象？感激不尽！

@RQLI_CFD 感谢，我研究研究

一般泰勒格林涡的文章都会指定计算域长度与涡的大小啊。比如这个：https://www.jianguoyun.com/p/DblYnzoQ9s3ZBhil274FIAA 方程6.2

我感觉有必要擦除流体与固体的试算结果，再进入下一个迭代步。固体好说，只有一个D，我令D=D.oldTime();就能把D的试算结果消去，但流体就很麻烦，因为有很多个变量，我试了半天也不知道该怎么把所有变量的值都擦除掉，恢复迭代前的值

没太看明白

@Qing 感谢老师！您说的特别清楚，这对我特别有帮助。再次谢谢您！

也不知算不算高粘度，1-几百mpas都有，感觉速度对的还行，液体表面泡沫层是真的没有

拉格朗日一般要从粒子碰撞的角度来考虑破碎聚丙啥的

属实没看明白。有一些教材，sci，写的就很晦涩。

很有帮助。

@WeakForm 这是我微信liuzongku

@李东岳 李老师，请问如果求解非守恒型的动量方程，那么对于其中一相为0的情况，是不是就直接求解就好了，不用再做别的处理了，如让相为0的相的速度等于按照模型去求出的速度。

我自己在用TDMA算法，用Matlab

伴随优化我一直想渗入 就是一直没时间详细研究

@thegame 感谢老铁，将在下一版本更新

多谢两位老师