4. 翼型气动特性随攻角变化的二维仿真
问题定义与仿真策略
NACA0012 是经典对称翼型，其气动特性随攻角变化具有典型规律：小攻角下升力系数 Cl 随攻角近似线性增长，阻力系数 Cd 缓慢上升；当攻角超过临界值后，上表面流动分离加剧，升力骤降、阻力急升，即发生失速。本例通过一组定常 RANS 仿真，扫描攻角 α = 0°、4°、8°、12°、14°、16°、18°，提取 Cl、Cd 曲线并识别失速攻角。
仿真采用二维不可压缩定常求解，湍流模型选用 k-ω SST，因其对逆压梯度下的边界层分离具有较好预测能力。来流条件：弦长 c = 1 m，来流速度 U∞ = 30 m/s，基于弦长的雷诺数 Re ≈ 2.0×10⁶。计算域采用 C 型拓扑，翼型表面布置边界层网格，第一层高度满足 y⁺ ≈ 1。
Cl、Cd 监测值稳定。
攻角扫描与系数提取
每个攻角独立初始化或从上个攻角结果续算，可显著加快收敛。升力系数与阻力系数由翼型表面压力与黏性力积分得到：
Cl = L / (0.5 * ρ * U∞² * c)
Cd = D / (0.5 * ρ * U∞² * c)其中 L、D 分别为垂直于来流和平行于来流的合力分量。下表为示例计算结果（教学示例数据，用于说明趋势与判读方法）：
| 攻角 α (°) | Cl | Cd | 升阻比 Cl/Cd | 状态判读 |
|---|---|---|---|---|
| 0 | 0.000 | 0.0082 | — | 对称，无升力 |
| 4 | 0.438 | 0.0091 | 48.1 | 附着流 |
| 8 | 0.872 | 0.0125 | 69.8 | 附着流，线性段 |
| 12 | 1.245 | 0.0210 | 59.3 | 接近临界 |
| 14 | 1.318 | 0.0295 | 44.7 | 轻微分离 |
| 16 | 1.102 | 0.0480 | 23.0 | 失速后 |
| 18 | 0.865 | 0.0720 | 12.0 | 深度失速 |
由表可见，Cl 在 α ≈ 14° 达到峰值约 1.32，随后在 16° 明显下降，说明失速攻角约为 14°~15°。升阻比在 8° 附近最高，是巡航效率最优点。
压力分布特征与失速机理
通过提取不同攻角下翼型表面压力系数 Cp 分布，可以直观看到失速过程。小攻角时，上表面从前缘到后缘压力逐渐恢复，逆压梯度平缓；接近失速攻角时，上表面后段出现明显压力平台，表明边界层分离；失速后，分离点前移，上表面大面积低压区消失，升力下降。
建议在结果后处理中绘制以下曲线进行对比：
- α = 8°：上表面
Cp呈典型吸力峰，后缘压力恢复良好； - α = 14°：上表面后段
Cp曲线出现拐点，分离泡形成； - α = 16°：上表面
Cp平台区扩大，前缘吸力峰降低，升力显著下降。
Cl 在 14° 后仍继续上升，常见原因是网格分辨率不足或湍流模型过度耗散，未能捕捉分离。可加密翼型后段边界层网格，或改用 γ-Reθ 转捩模型验证。
结果验证与工程判断
将计算得到的 Cl、Cd 与 NACA0012 实验数据（如 Abbott 数据）对比，小攻角线性段误差通常可控制在 5% 以内；失速攻角附近误差可能达 10%~15%，属于 RANS 定常方法的固有局限。工程上可据此判断：
- 巡航攻角应选在升阻比峰值附近（本例约 8°）；
- 失速裕度需保留 2°~3°，即最大使用攻角不超过 11°~12°；
- 若需精确预测失速后行为，应改用 URANS 或 DES 进行瞬态仿真。