有关low Mach的疑问

[李东岳]

RPF为什么多换一项rho就变成了transonic？

这个目前从理论来讲暂时不知道原因。我也在找关于transonic算法的文献。我目前觉得是要修正可压缩性，因为密度存在变化。

RPF为什么只能跑transonic，哪里限制了它sonic，supersonic和hypersonic？

rhoPimpleFoam属于压力基求解器，类似原始的PISO算法是可用于可压缩流和不可压缩流的。但是目前很少有人将PISO用于超音速流。超音速那面用密度基的比较多。这个在sci论文里面讨论的比较多。参考rhoPimpleFoam的介绍Transient solver for laminar or turbulent flow of compressible fluids for HVAC and similar applications.。主要用于可压缩流。

rho_0等求解器自动储存了，其他代码问题有空我再看看。

[程迪]

没有theory guide的弊端。。。

[noringname]

@李东岳 谢谢前辈！

[noringname]

@李东岳 顺便再问一下，到底什么是pressure-based和density-based，网上搜了不少，还是不很清楚。
我现在理解是解泊松方程的比如BPF和RPF这俩都算是pressure based。density-based完全没概念。

[李东岳]

@noringname
最浅显的，压力基求解压力，密度基求解密度这个是知道的吧？

[noringname]

@李东岳 嗯嗯，暂时只知道pressure-based是: momentum-->u, 泊松方程-->p，能量方程-->T, 然后p=rhoRT求得rho。
density-based网上搜的说是所有方程一起coupled解，具体怎么执行的不知道，还在搜索资料。

[李东岳]

密度基也可以分离求解，比如rhoCentralFoam，所以压力基分离求解器和密度基分离求解器，主要就是求解变量的不同。
密度基耦合求解就是密度、速度、能量一起解，构成一个大方程组，就像速度（每个速度具有3个分量）。目前OpenFOAM还没植入这个算法，好像extend里面有coupled算法，暂不清楚。

目前还没遇到有关基于OpenFOAM做coupled求解器的文章，都是分离式求解，要是看到可以共享下。

[noringname]

@李东岳 Thx!

[程迪]

@李东岳

其实密度基/压力基，耦合/分离，显式/隐式 是几个正交的形容词。

比如foam_extend的coupledPU就是耦合的压力基
常规OF里的pimpleFoam就是分离的压力基

DensityBasedTurbo就是显式的耦合的密度基
OF的rhoCentralFoam就是显式的分离的密度基

密度基和压力基求解压力时用的压力泊松方程/压力修正方程 还是 状态方程，压力泊松/修正方程是微分方程，状态方程式代数方程，所以差别比较大。

耦合/分离看你的主变量在每个迭代/时间步里是同时更新还是分别更新的。

隐式/显式看你的空间离散用没用未知时间步的量。具体隐式还分点隐/全隐之类的玩意儿。甚至可以半隐式。

[李东岳]

@noringname
参考楼上论述更详细。
@程迪
我不觉得 压力泊松方程/压力修正方程 还是 状态方程，压力泊松/修正方程是微分方程，状态方程式代数方程，所以差别比较大。有什么差别。

[李东岳]

Ah okay, 你说的是压力方程和状态方程，那是的。我以为你说的是压力修正方程和压力泊松方程 :confused:

[noringname]

@李东岳 @程迪 恶补了一下密度基，追问一个问题，为什么密度基解不可压的时候不准，为什么压力基解可压的时候不准，有没有什么相关的资料可以看的，只知道这两个不准是大家纷纷改进这两个方法的motivation。