使用compressibleinterfoam可以开启transonic吗？

[mohui]

最近打算用compressibleinterfoam求解器算喷雾过程的激波现象，但是不知道需不要开启transonic。

[李东岳]

因为要分析算法并且和实验对比，我找了很久多相超音速/亚音速的文章（双流体模型），没找到，你那有多相+可压缩的文章么？VOF也可以。

transonic算法非常简单，OpenFOAM里面只是植入了进去，目前Henry也没验证过准确不准确。迫切需要验证，这就回归到原始了，你那坐的是多相高速可压缩流动？

[mohui]

@李东岳 http://www.sciencedirect.com/science/article/pii/S0307904X17300410 这是我所知道的做多相可压缩的文章，另外我也和这个作者一直有过联系，他说是可以算激波的。这是他给我看的图片。

[mohui]

@李东岳 东岳老师，请问您有什么资料关于transonic算法的吗？

[李东岳]

那个文章我还没细看。
transonic我在Peric那本书里面可压缩流章节看到类似的处理方法，但是还不一样。目前还没看到有关OpenFoam里面transonic算法的理论来源书籍

[mohui]

@李东岳 知道了，今天上午google了一下也没发现什么对应的资料。对了东岳老师，如果说算喷雾的激波，给入口边界是个压力边界还是速度边界呢？假如给的是压力边界的话，速度增加的慢些需要算较长时间才会出现激波，如果给速度入口的话，一开始射流的前端会出现激波面，但是后面容易发散，所以不知道如何做了。

[yuan_neu]

@李东岳
请问您知道这个作者有公开代码吗？

另外我还碰到一个问题，对于可压缩气相的空化流，涉及到负压问题，MULES的算法好像没有隐性的fvm表现好。我很想看看这个作者的求解器里面的相方程的处理。
他的论文里面也提过，这个三相可压求解器是基于compressibleinterfoam开发的。而compressibleinterfoam里面只用了mules explicit的方法。

[yuan_neu]

@mohui 貌似很多diesel nozzle方面的计算 都是用的压力入口压力出口吧

[李东岳]

最近有关于transonic那面一些较新的文献么

[mohui]

@东岳 nope，最近没看文章了

[李东岳]

才发现这个帖子里面原来俩个人。你俩头像一个颜色。

https://www.sciencedirect.com/science/article/abs/pii/S0307904X17300410 这个文章里面也没有这个图啊 ，还是scihub的accept版本里面没有？

[李东岳]

东岳老师，请问您有什么资料关于transonic算法的吗？

更新一下近期对OpenFOAM中transonic算法的理解。经过梳理几个OpenFOAM中附带if (transonic)算法的求解器分析，目前从代码中获得的猜测是：

由密度方程
$$\begin{array}{}\text{(1)}& \frac{\mathrm{\partial }\rho }{\mathrm{\partial }t}+\mathrm{\nabla }\cdot \left(\rho \mathbf{U}\right)=0\end{array}$$
组建压力方程的时候，是否体现压力对流项。也就是密度方程中的对流项中的密度是否有通过状态方程变为压力项，比如
$$\begin{array}{}\text{(2)}& \mathrm{\nabla }\cdot \left(\rho \mathbf{U}\right)=\mathrm{\nabla }\cdot \left(\psi p\mathbf{U}\right)\end{array}$$
是否进行方程$\text{(2)}$的变换取决于压力传输效应。

在transonic与sonic中由于压力对流太明显了，因此需要进行方程$\text{(2)}$的变化，进而可以对其进行关于压力的隐性离散。一方面隐形离散时间布更加精准，再一方面耦合更加细致；在transonic以下的流动区域，压力对流弱爆了，没有理由进行方程$\text{(2)}$的转换，更何况这种操作会改变压力矩阵的对称性。

以上是个人从方程上的理解，没有见过书籍系统的论述过。抛砖引🐟吧！

[mohui]

@东岳 这篇文章是没有，这个图片是他个人私聊中分享的

[李东岳]

还有这种神操作