@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 感谢老铁，将在下一版本更新

多谢两位老师

@李东岳 感谢东岳老师的回复！

有个方向，叫混沌之类的。详细可以了解下蝴蝶效应和非线性的一些内容。

请问各位前辈，有没有界面蒸发源项的模型啊，找了好多都是壁面成核的，因为要做的工况是在气液界面处产生相变，而不是壁面加热导致相变的，看了一些文章好像找到没有one-group下界面蒸发的模型:136:

就是下面这个
$$
\frac{\p T_{ij}}{\p x_i\p x_j}
$$
其实$T_{ij}=u_iu_j$，所以演变为
$$
\frac{\p u_iu_j}{\p x_i\p x_j}
$$
这个公式是可以写成$\nabla\cdot(\nabla\cdot(\bfU\bfU))$这种形式，$\frac{\p T_{ij}}{\p x_i\p x_j}$好像写不出来 :134:

大佬们看看？