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

[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

[CFDngu]

OLAFlow?那我假设你是做Ocean/Coastal Engineering的，那这个应该是关于波浪的问题？波浪-结构物相互作用的问题？
个人感觉OLAFlow对interFOAM的改动不算大，一般来说，你的这个可能问题可能是由于MULES算法，自由面出现的wiggles