膜冷却 — 常见CFD陷阱与对策

这是RANS的典型症状。首先检查湍流模型。标准k-ε会过度估计冷却喷流的横向扩散，应改用SST k-ω或Realizable k-ε。如果仍然过大，尝试将湍流Prandtl数 $Pr_t$ 从0.85提高到1.0~1.2。

网格可能过粗，导致无法分辨喷流对壁面的重新附着。如果孔下游区域 $\Delta x / D > 0.5$，应细化网格。此外，如果省略了孔的充压室（供应室）并使用均匀流入条件，孔内速度分布会与实际不同，导致 $\eta$ 被低估。

这在吹气比较小（$M < 0.3$）或主流压力波动较大时可能发生。首先检查边界条件的设置。如果冷却孔入口使用了mass flow inlet，尝试改为pressure inlet。逆流在某些情况下在物理上是合理的，需要结合孔的L/D比和供应充压室的压力分布来判断。

DES和LES的数值方案选择很关键。对流项使用First Order会导致数值扩散而LES模式无法正常工作；使用Central Difference则不稳定。推荐用Bounded Central Difference或Blended Second Order/Central。Fluent中选择「Bounded Second Order Implicit」，STAR-CCM+中选择「Hybrid 2nd/Central」。

基本要求是Courant数不超过1。用孔直径 $D$ 和主流速度 $u_{\infty}$ 估算 $\Delta t < D / (20 \cdot u_{\infty})$ 作为初值，然后根据Courant数分布进行调整。可用STAR-CCM+的Convective Courant Number场或Fluent的CFL等高线来可视化。

以冷却喷流主流方向的流通时间 $\tau = L_{domain} / u_{\infty}$ 为基准，至少舍弃初始过渡的10~20 $\tau$，之后用30~50 $\tau$ 进行时间平均是标准做法。