激波-边界层干涉 — 故障排除
更丰富的内容请查看 shock-boundary-layer.html。
收敛失败·发散
SBLI分析中计算发散时,首先应该检查什么?
为您介绍可压缩流发散原因的前三名。
1. CFL数过高:密度基础求解器的显式时间推进中CFL < 1.0为稳定条件。即使在隐式求解器中,对于SBLI这样的强不连续也应从初始CFL=0.5开始逐步增加以确保安全。
2. 网格质量不佳:激波通过的单元宽高比过大会放大数值振荡。特别是在激波与边界层外边界交叉处,棱柱层顶部到四面体层的过渡部分网格质量至关重要。
3. 初始条件不当:从均匀流开始会导致强激波引起初期的非物理过渡状态。建议使用较低马赫数的求解结果作为初始条件,或使用坡函数逐步提高马赫数。
激波模糊
激波被模糊化,扩展到多个网格单元中,无法被清晰捕捉...
让我整理一下可能的原因和对应措施。
| 现象 | 原因 | 对策 |
|---|---|---|
| 激波跨越5个或更多单元 | 使用一阶精度格式 | 切换为二阶精度迎风格式(Roe、AUSM+) |
| 激波位置依赖网格 | 网格太粗糙 | 将激波附近的网格密度增加2倍 |
| 激波前后出现振荡(Gibbs现象) | 未使用限制器的高阶格式 | 应用TVD限制器(van Leer、Superbee) |
| 斜激波呈现阶梯状 | 网格与激波未对齐 | 将网格线与激波方向对齐 |
有没有限制器选择的窍门?
van Leer限制器很安全但散逸较大。Superbee限制器很锐利但容易振荡。Minmod限制器散逸最大。对于SBLI,压力场通常使用van Leer限制器,如果激波仍然模糊,可以切换到Superbee。这是标准做法。
剥离气泡未被再现
实验中有明显的剥离,但CFD中没有出现剥离...
这是RANS分析中常见的问题。检查要点如下。
1. 壁面网格是否充分:当 $y^+ > 5$ 时,壁函数起作用会过度估计壁面剪应力,抑制剥离。必须严格遵守 $y^+ < 1$
2. 湍流模型的选择:SA(Spalart-Allmaras)模型倾向于低估剥离。切换到SST $k$-$\omega$ 往往能显著改善
3. 二维分析的局限性:SBLI剥离本质上是三维现象,二维RANS分析的剥离长度可能与实验不符。使用三维网格配合周期边界条件可能会有所改善
4. 考虑转捩:在假设完全湍流的情况下解边界层,反而会导致剥离过小。建议考虑使用Fluent的 $\gamma$-$Re_{\theta}$ 转捩模型或STAR-CCM+的Gamma Transition Model
即使将壁面 $y^+$ 设置为1以下也没有改善,是否只能转向LES/DES?
在RANS范围内还可以尝试一个方法,即启用QCR(二次本构关系)。标准Boussinesq假设无法表现的二阶应力项可以通过添加来改进角流和强逆压力梯度下的预测精度。在Fluent中可以为SST $k$-$\omega$ 附加QCR选项。
激波-边界层干涉分析中"剥离气泡"消失和出现的谜团
在激波-边界层干涉CFD分析中,初学者经常遇到的一个困惑是:每次网格细化时,剥离区域的大小都会改变。粗网格显示"无剥离",但精密网格会显示出小剥离气泡。这不是故障,而是剥离尺度非常小,当壁面第一层的Y+较大时无法捕捉到。激波后面的逆压力梯度极其陡峭,需要Y+≤1左右的网格来解决。实务中"首先检查Y+分布图,确认边界层内有足够分辨率"是故障排除的第一步。
相关主题
细节
错误