OF两相流中是否能够实现只在气相求解某一个方程

[李东岳]

双流体模型？

[WangQi 0]

@李东岳 应该不是李老师。具体问题是这样，扩散项的数值扩散问题太大，我实在是无法解决（我看到了Cfd界发的那个文章，但是对于我并不适用）。 但是通过我和外导的数值推演之后，我们想到了另一种解决办法。 就是希望气液界面及界面以上的数值都是0，只解开界面下方，气相的扩散方程，这样就可以避免气液两侧扩散系数差距太大导致的数值扩散。
但是在of中我目前还没找到合适的办法实现它

[李东岳]

就是希望气液界面及界面以上的数值都是0，只解开界面下方，气相的扩散方程

这个看起来涉及到一个移动边界条件的问题?

[WangQi 0]

@李东岳 对，就是想把界面的以及界面以上的液相电势全部设置为0

[李东岳]

你算的是什么。

就是希望气液界面及界面以上的数值都是0，只解开界面下方，气相的扩散方程

这个可以处理。但是移动边界不好处理，要么动网格，要么多计算域求解动边界，要么浸没边界，要么就是一些新算法了，要找找sci

[李东岳]

我刚才找学生要我们写的一篇文章，之前我们做过这个事情。但是学生说帖子里面“扩散性引起的问题”。我一看真是。你这个是气液界面的扩散项导致的，还不确定行不行。扩散项本身就是要扩散的啊，你想让他不扩散么？那只能多计算域了。

[WangQi 0]

@李东岳 多谢李老师，请问openFoam中，您提到的这个多计算域如何实现？有什么关键词吗

[李东岳]

@WangQi-0 openfoam里面存在一些多计算域模拟，比如chtMultiRegionFoam。

就是想把界面的以及界面以上的液相电势全部设置为0

你具体这个是算什么？求解液相电势方程么。上面是空气还是其他的？

[WangQi 0]

@李东岳 我想做沸腾过程，下面是空气，上面是液体，其实是需要求解整个场的电势方程。奈何气液电导率差距过大，所以一直解不准

[李东岳]

看起来属于这部分内容，我之前手算过类似的东西（附录A.3 扩散系数的影响）。是，那是要多计算域求解。但是你这个还涉及到一个动边界的问题。其他sci里面怎么处理的？

[Antutu]

遇到了类似的问题，就是水平集模拟两相流时，怎么样利用水平集计算的相界面把一个对流扩散方程隔离开，只在一个相内求解对流扩散方程。