CFD中文网

@两月三年 你现在用的编译器是VC98？ 不知道其他版本会不会解决这个问题。

欧拉多相流的话不是应该是网格尺寸大于颗粒尺寸吗，和气泡也会有具体的关系吗？

@promisegwj 啊不对 是用拉梅系数···

@wangfei9088 非常感谢:xinxin3:

按网格类型分类的软件列表：
Mesh Generation & Grid Generation on the Web

经Gmsh论坛解释，Gmsh不支持该功能，建议用meshlab等软件先进行几何前处理。
弃坑，走人

学OF时间长了就越明白OF编程的这种独特，刚开始入门知道用的是C++我还找了以前学C++的书出来重刷，但是后来发现OF虽然也是基于C++但是自己很有自己的特点，所以如果C++不好的入门者也不用担心，耐心点慢慢看OF的码有不懂的就网上搜一搜别人的解析还是能理解的。

https://askubuntu.com/questions/1143840/downgrade-openmpi-v2-1-1-to-v2-0-2

https://github.com/open-mpi/ompi/issues/6300

有空测试下

@白礼耕 好的，谢谢您

@chengan-wang 没做过分段式的，单个翼型使用blockMesh生成结构化网格是可以的，网上有可以找到python或matlab用来生成blockMeshDict字典文件的脚本，你可以试试多段翼型的网格设计。

https://www.cfd-china.com/topic/5198

参考这个帖子，有人可以做个STAR-CCM、COMSOL、CFX的算例么？

@tidedrinker 谢谢大佬！

最后一张图是X方向的受力，也就是船的阻力。
谢谢回复:146: ，如您所言，setFields确实有问题，我设置的值正好等于流体域边界值，当时没有检查其实边界处的相分数没有完全赋值。将setFields中的box范围改成稍大于流体域就可以正常求解了。我本来还以为是数值振荡这么严重呢，都想离散格式精度调低看看了。现在看来，还是物理模型导致的振荡才更严重。:136:

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

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

想开个贴讨论下多相流界面追踪类方法目前的前沿水平，例如能算到的最高参数，各类方法都有什么限制等等。目前多相流算法里，计算界面追踪的多相动力学，比如droplet collision, droplet splashing, droplet impact on solid wall，这几个典型案例，目前有用VOF，level-set，phase field，LBM, SPH, 欧拉拉格朗日耦合其它模型的。这几个案例在不同Re和We下会产生许多种演化形态，实验都做出来了，但感觉计算还没跟上。

调研过不下百篇的文献，基本上关键词相关的都翻了一遍，在文献中看到droplet collision大概最好的也就是Re=6500,We=1500，密度比DR=1000左右的，大致是这个级别，像droplet impact有篇08年的用TransAT算到了RE=11000，而最近16年左右的用levelset这些也算到Re=8000，We=1500左右。大多数都是采用了一些改进的数值技巧来稳定模拟，而且网格也都比较密，计算资源消耗也都比较大，估计都得上超算，当然和实验也都对的上，基本都是到了prompt splattering这样的形态。

想讨论下大家都做过或者见过的这类模拟最高参数是多少? 或者往上算各类方法的限制在哪？ 因为多相流在算到高RE和We的时候会数值发散，稳定性应该是主要限制，还有例如VOF和level set会在薄液层处由于数值格式产生人工破碎的非物理情况，需要单独使用人工识别薄液层再稳定的方式，对不规则碰撞并不实用。

@bestucan 感谢回复

@李东岳 嗯，两个一起下下来就可以的，之前理解错了。

如果不考虑时间项
\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}$把所有的面连起来，看起来是守恒的，剩下的只有边界条件

哦，明白了。