找到vof中interface的位置

[linhan.ge]

@东岳 最终还是通过简化力的模型来解决这个问题。正如你之前说的，通过找到界面进行几何分析来计算相应的力需要网格分别率很高， 但是再小的网格也比颗粒要大，这就无法提供足够精确的几何关系。

[litong189456]

@队长别开枪 您是做带有几何重构的VOF算法的，我想请问一下，现有的算法都是在不可压缩流体中使用，对于可压缩流体，这些几何重构算法（比如PLIC）还能适用吗，为什么现在没有相应的文章呢，是没有人做，还是说在算法上有难以解决的问题呢，这个问题困扰我好久了。

[队长别开枪]

@litong189456 目前我的工作还没考虑气相的可压缩性，考虑可压缩性后，控制体内的密度就不再是VOF的函数了，一些控制方程也要改，应该不容易搞

[litong189456]

@队长别开枪 谢谢

[linhan.ge]

@队长别开枪 最近想钻研下isoadvector，想请教大神一些基本问题。
相方程中的时间项在isoadvector代码的哪个部分求解的呢？找了一圈只找到advection那项。
@队长别开枪 在 找到vof中interface的位置 中说：

求解过程中的0.5等值面提取没搞过啊，只有后处理的时候提取过iso-surface。你说的这个需要在求解过程中确定interface的位置，使用OpenFoam最新的iso-Advector或者PLIC-VOF吧，这两个都属于几何重构方法，求解过程中可以得到interface的位置信息。

这里，您说的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]

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