多相流VOF

[红豆沙]

@东岳 十分感谢哇，，嗯嗯可以啊，，就是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的速度，就算到了泰勒流啊。

[红豆沙]

@qhwj1990 所以 我也不清楚哇 哈哈哈 研究ing

[红豆沙]

@qhwj1990 我觉得 跟初始化时compute from的选择 有关，，我以前是从all zone 进行初始化或者 干脆不选
all zone初始化计算 还会报错，，因为主相是水，所以从空气入口初始化开始计算，算出的结果和你的一样

[红豆沙]

@红豆沙 这是从空气入口进行初始化计算后的结果 和你的一样，，刚好反过来了，液体充当了气泡形状

[红豆沙]

@红豆沙 是 从水相入口初始化的 说错了

[qhwj1990]

@红豆沙 是的，我就是从水入口初始化的，我也觉得这相当于气液两相反了过来，你可以从接触角上来看。实在想不明白，从气相入口初始化也是这样，关键是怎样初始化，内部流场速度都是均一的啊，为0.