OpenFOAM计算斜激波，波后压力振荡

[bit_hypersonic]

来流条件如图：

流场大致如下图，使用层流模型，Y+小于1，压缩拐角处存在回流区：

最开始想计算Cp和实验对比，结果发现了振荡：

进一步加密网格，振荡更严重了：

除了壁面附近有这种诡异的波动，激波波后也存在波动，想问一下各位大佬，这种情况应该如何调整？为何网格加密后，振荡更严重了？

[先先贝]

你用的rhoCentralFoam求解器么

[bit_hypersonic]

@先先贝 课题组的求解器，rhoCentral耦合化学反应的，，但是我这个算例没开化学反应，相当于原生的rhoCentral

[李东岳]

https://cfd-china.com/topic/8018

能做的话可以上传到这里我给你看看

[bit_hypersonic]

@李东岳 非常抱歉，隔了一个月才回复您
这个问题已经解决了，是库朗数的问题，我原本库朗数设置为0.2，后来改成了0.05，振荡就没了
有个问题想请教一下您，库朗数0.2已经很小了吧？课题组里的老师说OpenFOAM里算超声速，库朗数要设置的很小很小，这个是OpenFOAM的特性吗？Fluent好像没有这种问题吧？

[李东岳]

fluent我不清楚。不过你这个看看加pimple循环试一下。你可以看看这个图 http://dyfluid.com/class.html#id6

[尚善若水]

rhoCentralFoam有的CFL数都设置到0.02，成本太高了。李老师说的pimple循环在高版本才支持，低版本设置了也不起作用

[李东岳]

@尚善若水 老铁见过隐性的密度基超音速算法么 来个sci我瞅瞅

[李东岳]

我刚试了一下，这个文章里面的最简单的eSSPRK两步两阶方法，也没有震荡。多步高阶也可以试试。

STRONG STABILITY PRESERVING INTEGRATING FACTOR RUNGE KUTTA METHODS

[李东岳]

rhoCentralFoam有的CFL数都设置到0.02

一说到这个库朗数，我看了一下，这面是用最大波速来计算的CFL。我用compressibleCourant对比了下。二者一个0.17，一个0.14。还有点差距。波速是求解之前预估的，phi是求解之后。我估计是差距的来源。不过没有细看代码，自己想的。

[尚善若水]

@李东岳 算法研究的比较少，偏应用~

[bit_hypersonic]

@李东岳 好的，谢谢李老师！
（好想去您的OF课，可惜老板比较抠门，哭死）

[李东岳]

一说到这个库朗数，我看了一下，这面是用最大波速来计算的CFL。我用compressibleCourant对比了下。二者一个0.17，一个0.14。还有点差距。波速是求解之前预估的，phi是求解之后。我估计是差距的来源。不过没有细看代码，自己想的。

我今天在准备这个算例的时候发现，http://dyfluid.com/class.html#hyperbolic 特征速度和传输速度有2倍的差异。针对这种偏hyperbolic的方程可以用俩种方法计算库朗数：

1. 特征速度
2. 传输速度

所以，rhoCentralFoam应该就是用的特征速度计算的库朗数，而不是传输速度。如果想看传输度的库朗数，需要用compressibleCourant.H来计算。