Fluent 飞机仿真：跨声速运输机翼型阻力发散与超临界翼型压力系数分布

项目过程

本项目为某通用航空研究所的运输机机翼设计优化项目，需求是评估常规翼型与超临界翼型在跨声速段（M∞=0.65~0.90）的阻力发散特性，并定量给出设计巡航点（M∞=0.78）的升阻比与上翼面压力分布。客户给定：巡航马赫数 M∞=0.78、设计升力系数 Cl=0.55、雷诺数 Re=3×10⁷（基于平均气动弦长 1.5 m）、巡航高度 11000 m（大气参数按 ISA）。

本项目采用 ANSYS Fluent 2024 R2 求解 RANS + Spalart-Allmaras 湍流模型。翼型网格约 12 万单元，前缘 0.5% 弦长加密（y+=1），远场边界距离 50 倍弦长。求解设定：密度基求解器、AUSM+ 通量格式、隐式时间推进，CFL=200。巡航工况稳态求解收敛判据：残差下降 4 个量级 + Cl 变化 <0.001/MaxCl。扫描马赫数从 0.4 到 0.9，每 0.05 步长记录 CD。

案例分析

阻力系数 CD 随马赫数曲线显示典型”跨声速鼓包”：常规翼型（灰色虚线）在 M∞=0.72 后 CD 急剧上升，M∞=0.78 时 CD≈0.078，M∞=0.85 时 CD≈0.105，已突破阻力发散临界（CD=0.08）；超临界翼型（蓝色实线）在 M∞<0.65 段与常规翼型重合，M∞>0.65 后阻力增长明显放缓，M∞=0.78 时 CD≈0.022（仅为常规的 28%），M∞=0.85 时 CD≈0.05。M∞=0.78 设计点的”阻力发散裕度”从常规翼型的 0.06（已经发散）提升到超临界翼型的 0.13，巡航升阻比从 7.0 提升到 25。

M∞=0.78 设计点上翼面压力系数 Cp 分布对比：常规翼型吸力峰达 Cp≈−1.2（位于 x/c≈0.15），随后压力缓慢恢复至尾缘附近出现明显激波前兆；超临界翼型则呈现”前缘高吸力 + 中段平台”特征——前缘 x/c=0.05 处 Cp≈−0.85，x/c=0.20~0.60 维持平台约 Cp≈−0.85，到 x/c≈0.65 才出现弱激波（压力恢复坡度 -0.6）。弱激波位置后移 + 强度减弱是超临界翼型跨声速减阻的核心机制。

总结教训

本项目最容易踩的坑是 SA 湍流模型对激波的捕捉。SA 模型对激波位置预测偏后约 5% 弦长，且强度低估 10%。如果客户需要更精确激波位置，应升级到 SST k-ω 或带转捩的 γ-Reθ 模型。本项目客户接受 5% 偏差，SA 已足够。

另一坑是远场边界位置。50 倍弦长已经够用，但如果客户后续加装襟翼，需扩到 80 倍弦长避免襟翼涡干扰远场压力。

项目经验

经验一：超临界翼型减阻效益对巡航 M∞ 高度敏感。M∞=0.78 时减阻 72%，M∞=0.85 时减阻 52%；如果客户巡航点降到 M∞=0.70，超临界翼型效益仅 18%，性价比骤降。

经验二：上翼面平台 Cp≈−0.85 是判断超临界翼型设计成功的关键指标。如果平台 Cp 仍 >−0.6 或平台长度 <30% 弦长，说明超临界设计不足，仍可能跨声速发散。

经验三：弱激波位置 x/c=0.65 是襟翼铰链点设计的极限参考。襟翼铰链点必须位于激波下游，否则激波-襟翼干扰会引发俯仰力矩非线性，影响纵向配平。

有限元仿真