CFD中文网

@李东岳 感谢李老师的回复。刚刚我又试着修改了一下（代码在第一张图），但是这次修改过后运行结果仍然存在有问题，程序没有运行起来，请您指点。
![3699c109-7ea4-4bd2-a3d8-8143185f9510-_157$9_AFAWJ1{)8W35YL3.png
709379c4-f503-43a1-99b2-40762da25237-8L9F0U3P()586BOLA@{2NAP.png

@东岳 我最近用reactingMultiphaseEulerFoam做算例，也是设置的时候某一相为0，湍流模型一直出现divide by zero，找了半天没发现错误，更改之后搞定了，谢谢东岳老师

参考这个标题Ten years of industrial experience with the SST turbulence model，我想我可以写一个Ten years of academic experience with OpenFOAM CFD code :mihu: :mihu:

画周期网格不行吗，取出1/10画好就行了啊

@流动魅力 我比较好奇你的相间作用力的直径怎么给的

@FluidGao 我不是啥前辈:chouchou: 高校，博士后。

直接Filters-Search搜索就能搜到

经验之谈，如果是复杂翼型这种网格的话，用RANS算的问题应该不大，我算过挺多的了，开非正交修正，能算。
不过还是想想法，尽量提高网格质量吧。

@东岳 对于$\nabla{\mathbf{U}}$的矩阵形式应该是什么这个问题，从张量的角度来看我站第一种写法，更规范。:huahua:
但是从CFD角度来看，公式这东西本就是规则，制定了规则就是交流用的， 可读性和易理解性还是更为重要的。看个人理解，我觉得并无好坏。:mianmo:

$sdsfa$

sdsfa

@李东岳 原文作者是用CFX做的:136:

如果不考虑时间项
\begin{equation}
\int^{f+1}_ {f} \frac{\partial}{\partial x}\left(\alpha \frac{\partial \phi}{\partial x}\right)\mathrm{d}x=\left(\alpha \frac{\partial \phi}{\partial x}\right)_ {f+1}-\left(\alpha \frac{\partial \phi}{\partial x}\right)_ {f}
\end{equation}
把后面两项$\left(\alpha \frac{\partial \phi}{\partial x}\right)_ {f+1}-\left(\alpha \frac{\partial \phi}{\partial x}\right)_{f}$把所有的面连起来，看起来是守恒的，剩下的只有边界条件

@bestucan 多谢多谢，非常感谢！！！！:xiexie: :xiexie: :xiexie: 受益良多！

@随风飘姚 您好 我正遇到这个问题，我设置了outflow出口，但是进出口质量就是达不到平衡。请问可以和您交流一下么？

感谢分享！

@东岳 原来是连续方程的引用啊，谢谢李老师啦

@alex 看到你的网格了，感谢分享，感觉划分的非常不错，好厉害！学习了！

@G8S7

究竟能不能用欧拉-欧拉模型来说明这个体系下固体颗粒的运动情况，

如果有实验图片的话更好判断。取决于固体颗粒粒径，相分数，粘度等。如果用你所说的VOF+欧拉-欧拉才能描述这一个过程，也是可以的。另一种方式是三相欧拉，在皮层熔体中对界面进行尖锐挤压处理。

另外，用什么模型取决你你关注的重点，如果关注内部的运动，皮层溶体不是重点用欧拉无妨，VOF也是有耗散。

At last, 几何很简单，是验证算法的利器，不用OpenFOAM可惜了。