CFD中文网

@李东岳 @张晓宇 我最近也遇到了这个问题，但可以通过ICEM网格修复功能“check Mesh”，将其解决。

@李东岳 跟湍流项有关的参数，在发散的时候，在进口是正常的。这个case我也特意关注yplus+， 优化到<10左右。

附件：发散的各参数截图
https://drive.google.com/file/d/1U7dGKk-jv6OYMdOf8Ukiy1Gak648iJEU/view?usp=sharing

你好，
宝丁@宝丁 ，东岳老师@李东岳

很高兴能联系上你。

你提供了一个很nice的case。我利用该案例做了一些测试。

solver：rhoPimpleFoam
turbulence:SpalartAllmarasIDDES

以下记录一下我的测试过程以及结果：

结果：

初尝试：

一开始，我基本参考了https://github.com/Tushargh29/transonic_shock_buffet
算外流的LES的case，这个case 有非常详细的step-by-step的操作过程。
那我的操作就是，将这些设置原封不动的搬到你的naca65的case上面。
对于湍流的刻画，整体还是非常好的。但是，算到大概0.0017s(跟网格有关，我后来加密了网格，会算的长一些)。但无一例外，最终都会出现如下图的inlet发散问题。一开始，我还以为可能是inlet 网格加密存在不均匀导致的。后来，又重新将inlet网格均匀化，之后仍无果。

寻找原因：

经过相关文献的调研，我找到18年的一篇文章。题为：A characteristic inlet boundary condition for compressible, turbulent, multispecies turbomachinery flows。

这篇文章，摘要就提到了LES高精度仿真叶栅，需要考虑“non-reflecting boundary conditions”。也就是“声波” （还有涡波，熵波的影响）。

在这篇文章的引导下，我开启了二次尝试。参考你之前曾经在Oulet边界做过wavetransimission的边界条件。这相当于NSCBC无反射边界条件的简化。
这块工作，需要修改底层的求解器。目前，我还在为此努力。

新的突破

这个问题一直就卡在inlet发散的问题，在我查阅资料过程中，发现激波管的case 。这里面的介绍，其中有一段话：
“TVD schemes are highly recommended when you are dealing with strong discontinuities (such as shock waves)”

根据该提示，我将 格式改成了 div(phi,U) Gauss MinmodV;
结果，之前无法收敛的求解器，可以继续计算下去了（大概计算了好0.2秒了，后面由于资源问题，就自己将其停掉了）。

对比添加MinmodV前后

turbulence的影响，在添加MinmodV后感觉趋向于laminar。inlet无发散现象。

我是做气动声学的，若不考虑声学边界条件，其结果还是不准确的。
目前仍非常嫉于NSCBC无反射边界的工作，还在推进。。

以下是测试算例以及结果：
https://drive.google.com/file/d/1vCfITb8qeTXG5_yqTn_va7vhlUTw95mz/view?usp=sharing

祝收敛。