找到vof中interface的位置

[队长别开枪]

@linhan-ge 绝大部分几何重构方法都是使用一阶欧拉显式方法处理VOF方程中的时间项的，具体到代码就是

    // Initialising dVf with upwind values
    // i.e. phi[facei]*alpha1[upwindCell[facei]]*dt
    dVf_ = upwind<scalar>(mesh_, phi_).flux(alpha1_)*mesh_.time().deltaT();

    // Do the isoAdvection on surface cells
    timeIntegratedFlux();

    // Synchronize processor patches
    syncProcPatches(dVf_, phi_);

    // Adjust dVf for unbounded cells
    limitFluxes();

    // Advect the free surface
    alpha1_ -= fvc::surfaceIntegrate(dVf_);

fvc::surfaceIntegrate(dVf_)就是用来计算当前时间步流过控制体表面的主相体积占控制体体积的比值。

isoAdvector获取重构好的interface比较困难。在PLIC方法中则会定义交界面为$\vec{n} \cdot \vec{X} + D_0 = 0$，通过求解$\vec{n}$和$D_0$得到interface的近似平面。我们组最近会开源一个基于PLIC方法的二维求解器，在里面可以很方便地得到这个interface的空间位置，三维的预计年底开源，文章正在审。参考 https://doi.org/10.1002/fld.4664

[李东岳]

这个蛮好的嘛

[队长别开枪]

@东岳 在前人的基础上做的研究，只求老板给毕业

[linhan.ge]

@队长别开枪 非常有用的信息，期待你们求解器的开源。

[李东岳]

@队长别开枪
你主要做界面这面的哪个模型？PLIC？

[队长别开枪]

@东岳 做的比较杂，主要是PLIC，提高多边形多面体网格中PLIC的求解效率，PLIC + overset mesh，使用PLIC研究气蚀，就这些

[fireztw]

@队长别开枪 队教授 我目前在使用VOF代数重构的方法在做相关研究，对于几何重构的优点不是很明白，PLIC算法在of中的应用能否提供简单的例子学习一下啊，非常感谢

[队长别开枪]

@fireztw 给个图吧
左边是isoAdvector，属于几何重构方法里的一种，右边代数方法。纠正一下说法，代数方法没有重构这一步的，优点在于求解效率和网格单元形状的适应性。

[nbyjn]

@队长别开枪 大神您好 刚读了上面您的PLIC-VOF的论文 请问这套代码开源了吗 谢谢

[队长别开枪]

@nbyjn 留个邮箱。

[nbyjn]

@队长别开枪 992906307@qq.com 万分感谢！ 另外想请教您一个问题 plic界面重构方法和动网格（没有拓扑变形的动网格）会有冲突吗 因为上文提到的isoadvector方法只能用于一部分动网格 不知道plic方法会不会有类似问题

[队长别开枪]

@nbyjn PLIC可以用于重叠网格，其他动网格技术还没试过，不好评价。

[wallong]

@队长别开枪 @nbyjn 最近在计算液滴撞击固体壁面的飞溅现象，使用Fluent可以产生飞溅，interFoam没有该现象（修改了壁面接触角模型），而且isoAdvector表面张力和压力平衡又不太好。考虑应该是VOF界面处理引起的，可以求一份PLIC的代码吗？1611689516@qq.com 不胜感激。

[队长别开枪]

@wallong 在 找到vof中interface的位置 中说：

1611689516@qq.com

已发。不过advection部分我们暂时使用的iso-advector的时间积分方案，可能对你的问题帮助有限。

[wallong]

@队长别开枪 好的，我研究一下，十分感谢！

[linhan.ge]

@wallong 在 找到vof中interface的位置 中说：

而且isoAdvector表面张力和压力平衡又不太好。考虑应该是VOF界面处理引起的

这个结论是您测试的结果还是有相关文献呢？

[wallong]

@linhan-ge 测试过，另外可以参看一下Johan Roenby的说法："For you surface tension people, I believe you will also need an improved estimate of the curvature, right? This information is somehow contained in the difference in surface normal direction between adjacent surface cells. Hmm... still thinking about how to cook this up in a consistent and robust way for general meshes..." 来自https://www.cfd-online.com/Forums/openfoam-news-announcements-other/180439-isoadvector-release.html

[mohui]

@队长别开枪 为什么基于openfoam做的，还有三维，二维之分？

[队长别开枪]

@mohui Splitting 2D and 3D yields two papers.

[李东岳]

@队长别开枪 屌炸...