多相流VOF

[红豆沙]

为何我的VOF模拟段塞流（或泰勒气泡）均不能成功，从文献中选取的相同的模型，参数，边界条件，但是 出来的相体积云图一直是 分层流，
我要做的是在高速气流下大管径中的段塞流
以下是一个微通道泰勒气泡的模拟，但是也未能实现!

下图是 文献中的达到的效果，我的气泡总是不能聚集

采用的显示算法，隐式体积力，全局库朗数小于2，自适应时间步长，

，

PISO算法，压力用PRESTO,动量二姐迎风，体积分数，几何重构，

求多相流 大神指教，，我刚入门多相流，要做的是 气液固多相流下冲蚀，，以前做的是涡激振动

[qhwj1990]

气相流量大了吧。

[红豆沙]

没有哇，，如果做泰勒气泡，要气流量小，都是速度入口，空气0.3m/s，水0.5m/s，，按照文献中信息，应该会有气泡。而且，接触角设置也为0度。主相是 空气和水 都试过了，一样结果。

[爱学习的LDH]

请问您，接触角的定义方法是通过udf定义，还是直接在物理模型中打开表面张力，直接输入水和空气的接触角即可呀？

[李东岳]

网格尺寸，重力方向都对么

[红豆沙]

网格尺寸正确，重力方向也正确，（我试过 加重力和不加重力都不行）

[红豆沙]

@爱学习的LDH 接触角在 表面张力设置面板中 有个 wall adhesion 勾选了后 在边界条件 wall哪里会有接触角设置

[李东岳]

网格数量多少？2D？

[爱学习的LDH]

@红豆沙 这个接触角使用之后，你有没有应用到动网格或者是重叠网格，在fluent中，接触角模型和重叠网格兼容么？

[红豆沙]

@东岳 网格数量有（6000,4w,8w）都是这样子，，是2D,想着先把2d调试出来，再做3D

[红豆沙]

@红豆沙 文献中，，网格在5k的时候，也是有泰勒气泡的，2d

[李东岳]

文献发上来看看？
那你也可以用5K网格试试。网格数量太大测试太慢

[红豆沙]

@东岳 嗯嗯，，好的

[红豆沙]

@红豆沙 手机不知道怎么传文献，我发给你

[爱学习的LDH]

@红豆沙 这个接触角使用之后，你有没有应用到动网格或者是重叠网格，在fluent中，接触角模型或者说表面张力模型和重叠网格兼容么？

[红豆沙]

@爱学习的LDH 我没有应用到 重叠网格或者动网格，但是在fluent 19.0中，重叠网格这个功能是加入到了多相流计算中的（粒状流动除外）如VOF、Mixture、空化、蒸发。

[红豆沙]

@红豆沙 那接触角模型和表面张力模型 应该是兼容的，，你可以试试 哈哈哈

[红豆沙]

@东岳 （气液两相段塞流数值模拟及其在结晶过程中的应用-李亚娟）文件超过500kb了，，

没注意到你的签名，冒昧了(⊙o⊙)…

[李东岳]

我找了一下，是cnki的格式，我这面看不了，你要模拟对标的是一个硕士论文？相关结果有没有发表的英文SCI？我可以参考英文的pdf版本

[红豆沙]

@东岳 无相关SCI，只有一篇是生物类的实验研究
因为硕博士论文较为详细，我是参考她文章中的段塞流部分。

[李东岳]

文献中，，网格在5k的时候，也是有泰勒气泡的，2d

如果可以的话，把你的网格处理到5千，然后把边界条件告诉我，我可以用openfoam给你做一下，

[红豆沙]

@东岳 enen 好的呢 谢谢啦

[wwzhao]

@东岳 硕博论文可以到中国知网的海外版下载pdf，网址：http://gb.oversea.cnki.net/

[红豆沙]

@东岳 文件为 画好的.msh网格，

[红豆沙]

@东岳

文献中的网格尺寸模型，我把50d这个长度缩短了，层流和湍流都试过。

[李东岳]

@红豆沙 网格上传失败?

https://www.sciencedirect.com/science/article/pii/S1385894717312433
我发现一篇泰勒泡的文章

[红豆沙]

@东岳 嗯嗯呢，，我看看去，，谢谢啦，FKG2.zip

[李东岳]

这是我用OpenFOAM算的，网格很糙，你觉得这个看起来更合理么？

• 从你第一个图来看，你检查检查你的表面张力设置？

• 从我模拟的经验来看，50d不能缩短，气泡需要时间发展才能出最终的形状，你看我下面这个图大体就是这个过程

[红豆沙]

@东岳 嗯嗯 看起来合理哇，有气泡形成的说，，，，我为了算的快点 就缩短了 ，，，但是您的 都有气泡噻，为什么我的没有哇，，我的表面张力设置的是0.072N/m，，您可以将您的设置说一下吗，，

[红豆沙]

@东岳 或者我用我 笔记本算一下，，笔记本装的是OpenFoam滴，，但是OpenFoam我还不熟练。

[李东岳]

就是很普通的设置，细节很多一一明说太费事了。你把你的Fluent的cas文件上传一下。你主楼那个云图看起来就是表面张力没起作用。

[红豆沙]

@东岳 嗯嗯 好的 我先自己再搞搞看，，实在不行，，就麻烦你们了

[红豆沙]

@东岳 FKG2.zip 确实 没算出气泡来，，不知是不是 入口条件设置问题，，我的气速是0.45m/s，液速0.7m/s

[红豆沙]

@红豆沙 包含了DAT的文件FKG2 (2).zip

[李东岳]

好的。看看有没有人对你的case感兴趣帮你调试调试

OpenFOAM算例我这几天处理一下发给你。要不你用OpenFOAM算吧。哈哈

[红豆沙]

@东岳 十分感谢哇，，嗯嗯可以啊，，就是OpenFoam 得捣鼓段时间 估摸才能算起来 哈哈哈，，

[李东岳]

主楼的图中，为什么你的气体没有上浮？你的重力方向对么。

我昨天那个重力方向设置反了，今天重新模拟了一下，下图：

[红豆沙]

@东岳 我没有设重力，，（昨天就感觉你的重力是不是设反了，哈哈哈，），，但是我做的加了重力的，也没有气泡，，我一会重新再试试，，再重新模拟一遍看看，检查检查

[李东岳]

@红豆沙 http://dyfluid.com/interFoam.html#testcase OpenFOAM的验证算例可在此下载

[红豆沙]

@东岳 哈哈哈 妈耶 我一阵激动 差点在实验室 叫出来了 哈哈哈哈， 实验室是全用fluent或自己编程（matlab），，，只能自己琢磨了 （赶紧打开笔记本 琢磨去 哈哈）

[qhwj1990]

@东岳 您的接触角采用的是多少？

[李东岳]

0度

接触角对气泡的形状有着重要的影响，用户可研究不同的接触角对液体润湿、段塞长度的影响；

http://dyfluid.com/interfoam.html#testcase

[maisui1993]

@东岳 在fluent中二维模型是默认为处理成z方向为1m的厚度。请问在openfoam中二维模型在z方向怎么处理的？

[qhwj1990]

同样的边界，我用Fluent计算的结果是这样的，气体怎么会冲到上部水的入口处？很费解。

[李东岳]

@maisui1993 OpenFOAM第三个维度方向可以任意处理成任意长度 你的结果怎么样？

[红豆沙]

@qhwj1990 你这个 可能是气体速度太大了，液体入口速度太小吧，，你还没算几步吧

[qhwj1990]

@红豆沙 我这是0.685s时刻的云图，红色为气相，速度均为0.1m/s。重力方向向下

[qhwj1990]

@qhwj1990 这和东岳老师的边界条件一致

[红豆沙]

@qhwj1990 我最多算到，0.几秒，气速比液速小，最后出现的是分层流，，时间步长会达到1e-5级别，（全局库朗数我控制在了小于0.25）你的长管中的云图很像我 刚开始计算时的云图，，，
原因 我也说不清额，，你可以试试减小气速或者加大液速呗 哈哈哈

[qhwj1990]

@红豆沙 关键我看东岳老师在OF中计算的两相均为0.1的速度，就算到了泰勒流啊。