两相流VOF自由面处总会出现锯齿状

[Mr.Shan]

接触OpenFOAM不久，现在终于会了一些简单操作，目前用开源求解器进行两相流的计算，在自由面处，总是会出现锯齿状，而且只有这个区域会出现，该处网格是均匀网格，0.05cm0.05cm，我把网格放大，变为0.2cm0.1cm，锯齿不会出现了，但是60秒后便不收敛了，想问问各位老师这是怎么回事

[李东岳]

Quite interesting. 在你这个图的基础上，显示一下网格？interFoam? 水合空气？

[队长别开枪]

把算例发来看看

[Mr.Shan]

@东岳 东岳老师好，求解器是一个基于interFoam的求解器，叫做olaFlow，模拟的是波浪-结构物-海床之间的相互作用，图片里展示的就是波浪的传播，第一幅图是波浪从不均匀网格处进入均匀网格，第二幅图是在均匀网格中，第三幅图是均匀网格的后半程，可以看到在不均匀网格区域，并没这个锯齿，进入均匀网格区域后，锯齿出现，然后又消失了。

[Mr.Shan]

@队长别开枪 不好意思呀，我们课题组属于涉密的，所以我的图片都只是其中一小部分，对不起了

[队长别开枪]

@mr-shan 你们动量方程对流格式div(rhoPhi,U)使用的哪个？还有PIMPLE控制里面momentumPredictor开了吗？

[队长别开枪]

@mr-shan 还有啊，你把第一次回复我的内容删了吧，以后可以说不方便，直接说出来不好。

[Mr.Shan]

@队长别开枪 我现在还处于摸索的阶段，对于fvSchemes和fvSolution下的配置还不是很清楚，我给你截个图

[李东岳]

olaFlow相比interFoam做了哪些改动？

[anubis]

我在fluent里面用CLSVOF的时候也碰到类似的情况， 原因我感觉是网格太小或者时间步长太大，你可以试试缩小时间步长。

[Mr.Shan]

@东岳 olaFLow是模拟水在孔隙介质中的流动，添加了一个孔隙介质对水的拖曳力项

[Mr.Shan]

@anubis 我也考虑到这方面了，因为我把网格放大，就不会出现锯齿了，但是这和我的认知有点不符，我对OpenOFAM的浅显认识，不应该是网格越小精度越高吗？怎么还会出现这种问题

[李东岳]

你的模拟里面有多孔介质么

[Mr.Shan]

@东岳 有的，除了孔隙介质，还有一个用snappyHexMesh切割出来的结构物，它们都在锯齿区域下方

[Mr.Shan]

@东岳 如果我不添加孔隙介质以及结构物，锯齿就不会出现

[anubis]

@mr-shan 跟库朗特数有关， 越小的网格精度越高，但是要保证一样的库朗特数就要用越小的时间步长。

[yhdthu]

看样子是压缩过度了，如果界面太sharp，CSF模型计算曲率会出问题，代数重构方法的通病

[李东岳]

@Mr-Shan 问题解决了么

[zhouhoucun]

试着把system/fvSolution 下的alpha里面的 cAlpha设置成0

[SnoWolf2018]

@anubis 支持+1