interPlicFoam

[myjimp]

@mohui 兄弟发一份interPlicFoam的论文，，谢谢

[队长别开枪]

@myjimp 我记得给你发过了啊

[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?