interPlicFoam

[myjimp]

@队长别开枪 是的，

[litong189456]

@队长别开枪 您好，最近在用您的PLIC工具，发现了一个问题，目前我还没有找到原因，不知您是否遇见过。
当我计算的时候，发现在纯水或者纯气体区域会出现另一相，二维和轴对称都出现了这个问题，用您的输出界面的writeplicfaces功能也发现，在本应该是纯水的区域产生了气体。黑色的是writeplicfaces输出的界面，粉色的是体积分数1e-8到9.9999999的单元，本来应该只在界面上。

[nanxuan]

@队长别开枪 前辈 ，您文章中这个后处理用什么做的，很清晰呀

[nanxuan]

前辈，这个plic几何重构对动网格里面的拓扑关系改变也是支持的吗，谢谢前辈，麻烦了

[队长别开枪]

@litong189456 这些都是极小值或者非常接近1的单元造成的，可以在控制字典里将snapTol数值调大，默认为0，可以试试1e-5之类的极小值把这些单元强制置零或者置1。

[队长别开枪]

@nanxuan 用的Paraview，不过PV有好多局限，最近买了Tecplot的许可证，以后会全部转为Tecplot。

[nanxuan]

@队长别开枪 好的谢谢前辈，我用的ensight，感觉这个也不错😊 ，前辈那个网格拓扑这个问题，像几何重构多段线或者等值线 应该会出现越界嘛？

[队长别开枪]

@nanxuan 正在开发的版本会支持AMR（附加VOF场专用插值算法）、overset以及动边界（造波器）网格。

[队长别开枪]

@nanxuan 不会的，在Foam::plicVofSolving::applyBruteForceBounding()函数里加了数值限制。

[nanxuan]

@队长别开枪 好的，多谢前辈，期待这个造波器，我用三维的造波使用cyclic周期性边界感觉老是怪怪的

[nanxuan]

@队长别开枪 前辈，如果使用滑移网格的话，interface前后的拓扑不一样，或者动网格remesh生成也不影响么？如果这个网格单元里面既没有液体也没有气体，他这想mules给他限制一个值嘛

[队长别开枪]

@nanxuan 需要说明的是我们不做造波器的，只是将PLIC-VOF方法适配已有的造波器，会有两个版本，一个是OpenFOAM plus本身的造波器，另一个是waves2Foam库里的造波器。如果你现有的算例使用MULES算法结果不理想，可能PLIC也不能改善，有个心理准备。

[队长别开枪]

@nanxuan 不支持滑移网格

[nanxuan]

@队长别开枪 好的前辈

[nanxuan]

@队长别开枪 前辈，我之前用iso-advect. 计算发现如果一个网格单元内如果没有气相和液相就会报错越界，mules却可以计算，不知道plic这可不可以，前辈这个版本我用的of6还不支持😁 ，也没测试效果，但是感觉理论几何重构这个算法，不像mules可以给他一个限制值适应性好(✪▽✪)。

[队长别开枪]

@nanxuan PLIC应该也不行。本来有计划把这个求解器适配所有OF版本的，但是精力实在有限，进度一直不咋滴。

[nanxuan]

@队长别开枪 是呀几何重构这个真是有了气液才能画多段线，如果没有气液只能用mules给他限制了，谢谢前辈哈

[Hungryandfool]

@队长别开枪
问题如下：

1. waves2Foam里面的这个适配做好没？
2. waves2Foam里面已经有了isoAdvector集成了，waveIsoFoam求解器。冒昧的问一句，您这个算法强在哪里？还没有认真读您的论文，抱歉。
3. 您适配之后，会有算例教程和设置参数吗？

ps: 您就是大名鼎鼎的daidezhi?

[队长别开枪]

@Hungryandfool

1. 代码搞完了，还在测试，最近很忙，每天下班只想睡觉。
2. 几何类VOF方法：
   a. IsoAdvector第一次提出是在2016年，我在OpenFOAM基础上开发PLIC-VOF始于14年下半年，不过开发进度很慢，那时刚开始学习OpenFOAM。有一个完成度很高的2D版本，大家有兴趣的话我会把这个代码整理整理，供大佬们批评指正，几何运算是基于CGAL库的；
   b. PLIC和IsoAdvector是两类不同的几何重构方法，具体区别在于如何表征混合单元内的interface。PLIC使用平面，IsoAdvector使用等值面。确定interface位置时PLIC求解平面方程$\overrightarrow{n}\overrightarrow{X}+D=0$里的距离$D$，IsoAdvector求解等值面阈值。界面重构精度应该差不多，求解效率还没比较，现在很难说谁高谁低；
   c. 我们求解器的公开版本advection部分直接借用了isoAdvection的，做了一些改动适配PLIC方法。主要是因为懒，不想单独开发advection部分了，不过这也符合开源精神吧，只做自己需要做和感兴趣的部分；
   d. 目前PLIC方法的一个研究热点是多面体任意单元里interface的快速定位技术。这方面研究主要集中在卡塔赫纳理工大学的Dr. Lopez [López, Joaquín, et al. "A new volume conservation enforcement method for PLIC reconstruction in general convex grids." Journal of Computational Physics 316 (2016): 338-359] (他们去年开发了适配凹多面体的快速算法，不过我感觉这个学术意义更大一些，基本上网格生成都是避免凹单元的出现的)，LANL的Dr. Diot [Diot, Steven, and Marianne M. François. "An interface reconstruction method based on an analytical formula for 3D arbitrary convex cells." Journal of Computational Physics 305 (2016): 63-74] 和Dr. Shashkov [Dyadechko, Vadim, and Mikhail Shashkov. "Moment-of-fluid interface reconstruction." Los Alamos Report LA-UR-05-7571 (2005)] (MOF方法本质还是PLIC方法，不过在$\overrightarrow{n}\overrightarrow{X}+D=0$的求解上都引入了全新概念)。
3. 会有算例教程，但是事情多，我尽量快一些，但是这些事只有我一个人在弄，所以会比较慢。

PS：我不大名鼎鼎哦，只是小卒，只想靠着CFD混口饭吃，学界和业界都太他妈卷了

[Hungryandfool]

@队长别开枪
好的。因为isoAdvector会成为我工作中的一部分，想深入学习了解下。谢谢您的热心解答！非常感谢！