求解:interFOAM网格收敛做不到



  • 最近用interFOAM和我自己的写的一部分代码计算了几个波浪方面的问题,对于波陡比较大的工况,网格收敛性做不到,波陡比较小的还可以,但是也不算太理想。我计算了逐渐加密的四套网格,每次加密后的网格尺度是上一级网格的一半,时间步长也变为之前的一半(固定时间步长)。我发现对于不是那么稳定的波浪,时间和空间分辨率过于精细时,计算的结果反而变差了?这是怎么回事?按照我的理解,即使对于一个不是那么稳定的问题,这种情况也不应该出现吧?提升了时间和空间的分辨率,结果反而更差了?一般的商业软件是不是不会出现这种情况?(如fluent等)

    此外,我知道interFOAM一直有些问题,如交界面出现褶皱,交界面相对速度不正常的大。但是这个应该和我的算例没关系吧。。。。:xinlei:

    有遇到类似问题的老师或同学吗?求交流指导



  • 对了,计算用的层流模型,且控制方程里面有源项,源项是显式的,在每个时间步的值是能算出来的。



  • mark一下,最近做气泡模拟也遇到了验证网格无关性的问题,加密了网格以后结果和原来的符合的不是那么完美,总有点差别



  • @CFDngu求解:interFOAM网格收敛做不到 中说:

    interFOAM和我自己的写的一部分代码

    自己写的代码其他人很难debug



  • @东岳 是的,目前我是这样想的: 可能越精细的时空分辨率,对细微扰动就越敏感,本来在上一个尺度下体现不出来的现象,在这个尺度下就能体现出来了。



  • 楼主您好,请问您提到的“interFOAM一直有些问题,如交界面出现褶皱,交界面相对速度不正常的大”的信息是在哪里获取到的?我最近在自己的计算中也发现了这个问题,但是一直不知道如何解释以及如何避免。。。。谢谢~



  • 表面张力计算一直是老大难的问题,以前算两相流的时候也出现过网格越密,求解精度越低的情况,后来发现是表面张力计算的问题,open foam的表面张力用的是很普通的csf模型,也没有针对湍流模型优化,至少les框架下那个表达不合理,忽略了部分亚网格表面张力应力,楼主如果要通过精确解析界面的方法做数值模拟,需要关注一下表面张力的计算


 

Forest
Mountains