OpenFOAM不同类型网格的计算速度

[tidedrinker]

有人比较过同样网格数量下，多面体网格、六面体网格和四面体网格在OpenFOAM（Fluent、CFX或其他CFD软件）中的计算精度和计算速度吗？

[李东岳]

木，一直想对比，2014年就找人算了个150万网格的算例，不过放弃了，我不懂fluent，不知道怎么比。

[程迪]

精度就不好比较，大家用的公式就不一样，残差不知道公式的话没法简单比较，OF类似L2残差/归一化因子，这个归一化因子我是没看太懂；然后Fluent也是归一化过得，归一化因子是采用的前5步残差的值算的。CFX不了解，SU2据说是真的L2残差。

最近想做一做这些软件的精度和速度问题，精度可以用MMS，有个人做过OF的，貌似有效精度只有1阶，有些边界条件精度也不高。

速度可以参考高阶精度CFD workshop的搞法，用DLR的taubench作为计算量单位来衡量。

reference: https://www.grc.nasa.gov/hiocfd/guidelines/

[李东岳]

是的，不同的软件之间对比起来难度非常大。无法做到单因素对比。不过也可以比着玩一玩。

[程迪]

@程迪

OF用的是归一化的L1模...每次归一化因子还不一样。

solverPerf.initialResidual() =
        gSumMag(rA, matrix().mesh().comm())
       /normFactor;

[程迪]

@李东岳
工程上有意义的是达到指定误差的计算代价，一般理论上认为误差可以从某种定义的残差推算出的不确定来衡量，时间以taubench用25000网格迭代10步作为1个时间单位来衡量。

这样的假设下做横向比较是可能的。但是，具体情况嘛，呵呵。湍流模型的不确定度就说不清楚。