OF可压流求解器

[程迪]

https://github.com/chengdi123000/compressibleFoam/tree/of41
把 http://pavanakumar.github.io/compressibleFoam 移植到了openfoam 4.1下。里面有roe格式。没啥大区别。

[youmengtian]

@程迪 其实也并不只是收收敛快慢的问题，还有更重要的是耦合关系的处理问题。

[nuaawqf]

@程迪 终于找到和我研究方向相同的FOAMer了！这个代码我也研究过，虽然思路很清晰明了，与我之前Fortran编写可压缩流求解器相同。但似乎没那么FOAM了？OF里面很多东西都没有用起来，似乎只是当成了一个编程语言在用。还有一个AeroFOAM也挺久没更新了，似乎面向对象的思想更多一些？不知是否研究过？　希望以后多多交流。:cheeky:

[程迪]

@youmengtian

耦合关系处理好了，不久快了么？

[程迪]

@nuaawqf
刚把AeroFoam v1-2008从1.5移植到了201612+，感觉基本上也是另起炉灶搞的。
AeroFoam v2-2013 最新支持到2.0，更是另起炉灶，fvMesh都给加功能改了，还没有移植成功。

[gyzhangqm]

时隔四年了，还是没有有基于Godunov的近似通量和隐式推进方法的开源OpenFOAM application。
已有的替代解决方案是：
1、使用rhoCentralFoam的KT和KNP格式，本质上还是压力基的求解器。可能存在的问题是耗散大以及收敛性差。
2、使用foam-extend的dbnsFoam，里边有气动人常用的HLLC格式和Roe格式，但推进方法只有龙格库塔显式推进。并且Roe格式熵修正处理不好的话会有Cabuncle现象。
3、使用论文Implementation of density-based solver for all speeds in the framework of OpenFOAM作者公开的allSpeedUnsteadFoam,空间格式使用的是带低速预处理的AUSM系列格式。高超声速热化学非平衡计算中很多人喜欢使用AUSM系列的格式。但allSpeedUnsteadFoam依然只有Runge-Kutta推进。

可以按照论文Implementation of density-based implicit LU-SGS solver in the framework of OpenFOAM 的方法实现LUSGS推进方法，工作量其实很小。如果原作者能开源其lusgsFoam最好了。

不过如果想进一步实现其他隐式推进方法，还是需要搞明白blockLduMatrix，并且组装Jacobian矩阵。我想对于foam-extend的大佬来讲，他们已经实现了blockLduMatrix类，以及后续的代数方程预处理及求解器，如果实现密度基的隐式推进方法应该是易如反掌的。不过等了四年依然没有看到有进展。

空闲时间我会自己着手做这件事情。有人同样在做的话可以一起交流合作~

[dzw05]

@gyzhangqm 看看这个代码，https://hisa.gitlab.io/index.html

[gyzhangqm]

@dzw05 十分感谢！
这个似乎可以完美应用于空气动力学亚跨超和高超问题的计算。自带AUSM+up全速域通量格式和HLLC格式，有LUSGS预处理的GMRES推进方法。貌似不是Jacobian Free的。还有作者的论文。非常好的资源。
我会慢慢测试几个定常的高超算例，确定比较合适的梯度计算方法、限制器等参数。大佬多多指教~

[heray]

还有一个myFoam，github可以搜到，里面有个lusgsFoam求解器