rhoCentralFoam求解shockTube算例出现数值震荡

[李玉棠]

大家好。我发现用rhoCentralFoam求解shockTube算例时x方向速度会有一些数值震荡。我调了网格数和maxCo数还是消不掉它。计算结果如下面两个图所示：

根据这个网站[Shocktube – rhoCentralFoam TVD Schemes Test – curiosityFluids](Shocktube – rhoCentralFoam TVD Schemes Test – curiosityFluids)我将限制器改为vanAlbada可是仍有数值震荡。
我将shockTube算例也上传了shockTube_test.rar 。我觉的应该哪里的参数配置需要调整一下，可以消除这个数值震荡。希望大佬们指点迷津。

[李东岳]

这个求解器有一个自己的算例，自带的算例没有震荡，你尝试一下

[李玉棠]

感谢李老师。这是我原始算例计算的结果。绿线是100网格./Allrun计算结果，蓝线是将网格改为800然后./Allrun的计算结果。其它参数都没有动。 这是算例shockTube_org.rar 我用的OF-4.1

[李东岳]

你时间步长减小没

[尚善若水]

时间步长太大了，0.02试试？还不行网格数改为1000

[李玉棠]

两位老师这是设置不同maxCo的计算结果。800网格

调小时间步长后，我觉得振荡被抑制的很好，但没有完全消去。似乎rhoCentralFoam的鲁棒性不是很好？我印象中有些求解器的CFL数设置在0.98都没有数值振荡的。不知道我的理解对不对

[尚善若水]

rCF是这样的，库朗数一定要足够小，单纯加密网格没改善。这是我的对比计算结果

一些可能的分析（一家之言）：

注：bF为blastFoam，rCF为rhoCentralFoam，sMF是我自己的求解器

[李东岳]

这个求解器有一个自己的算例，自带的算例没有震荡，你尝试一下

@李玉棠 我看原始算例，默认是1000个网格，也没有震荡。你尝试原始的算例没有

[李玉棠]

@李东岳 李老师好。对于shockTube算例，OF-4.1和OF-6是100mesh的;OF-8是1000mesh。我用OF-4.1算了OF-8的算例，只添加了turbulenceProperties字典，计算结果如图：
原始算例maxCo=0.2时Ux参数仍有震荡，调小maxCo后震荡被抑制。那位老师能给一下OF-8的计算结果吗？可能是版本的问题吗？

[李玉棠]

@尚善若水 谢谢老师。我觉得rCF这点不太好，，有没有更新的算法解决这个问题？现在有些既是高阶精度稳定性又好的格式，有能移植到OF的相关文献吗？

[李东岳]

我用OF-4.1算了OF-8的算例，只添加了turbulenceProperties字典

用of8添加湍流后，也有震荡么

[尚善若水]

@李玉棠 blastFoam好像就可以

[李玉棠]

@李东岳 老师我OF-4.1版本算例constant文件夹下得有turbulenceProperties字典，里面写的是laminar，不涉及湍流。好像OF-8以后都不需要turbulenceProperties字典了。