projection method 真的比piso快吗？

[wwzhao]

@mohui 感谢分享！

1. 试了一下cavity算例，可以正常跑完。

2. 将system/fvSolution里面的压力求解器的relTol设为0就可以使残差降下来，不过相应的矩阵求解迭代步要增加。

[mohui]

@wwzhao ，恩，我试试看。另外，请问你对这个有什么见解呢？噱头大于实际效果？

[wwzhao]

@mohui 正如你所说的，piso一般循环两次即可收敛，也就是要求解两次压力泊松方程，而rk4则需要求解4次。实际上对于不可压缩流体，求解压力泊松方程所需的时间占总计算时间的70%-80%左右，所以我觉得这方法对于一般问题没太大意义。

但那篇文章作者声称是这方法是low-dissipative的，我想对于一些特定问题应该还是有用的:cheeky:

[mohui]

@wwzhao， 恩，确实是。看来这个问题可以不用去深究，我本来打算移植到多相流求解器看下是否节省时间，从目前结果来看，貌似也没省多少时间。

[Foamer24]

[Zephyr]

我没看明白啊，他这不还是解了poisson方程吗，那不可能快的起来啊

[李东岳]

@mohui 在 projection method 真的比piso快吗？ 中说：

这两天看了一篇文献关于使用四阶龙格库塔时间推进代替PISO算法，他的结果显示要比piso快50%,而且精度还高。

只是一篇SCI，后续没准还有一篇SCI出现相反的断论呢，都不好说

[mohui]

@东岳 东岳老师挖了一手好坟啊

[Foamer24]

这个已经有开源的的openfoam代码了，但是是openfoam2.1.1的，不知道有没有比较会的能做版本适应

[mohui]

@Foamer24 这个你所说的开源是哪篇文献的？

[Foamer24]

@mohui 看到rhoEnergyFoam中有三阶四步的龙格库塔迭代。但是OpenFOAM大网格量并行速度相比fluent，不是很快

[mohui]

@Foamer24 有文献吗？另外这个求解器是官方的吗？我很久没关注新的版本的of了

[Foamer24]

@mohui 不是，OpenFOAM大网格量计算效率真的比不过fluent，好奇怪

[Foamer24]

@mohui 终于帮你找到啦！A LOW-DISSIPATIVE SOLVER FOR TURBULENT COMPRESSIBLE FLOWS ON UNSTRUCTURED MESHES, WITH OPENFOAM IMPLEMENTATION

[veen]

最近也基于pisoFoam植入了一下RKprojection foam，单个计算步不限制残差，projection method的确快很多，但是限制残差之后速度就变得特别慢，并且该算法内存占用量是pisoFoam的20倍。
不知道有没有大佬看到过把这个算法植入到多相流求解器的呢？想尝试一下把这个方法放到多相流计算里边去。