4. 二维翼型升阻力系数随攻角变化的参数化计算

案例任务与物理背景

NACA0012 是经典对称翼型，弦长 c 取 1 m，来流速度 U∞ = 30 m/s，空气密度 ρ = 1.225 kg/m³，动力粘度 μ = 1.789×10⁻⁵ Pa·s，对应弦长雷诺数 Re ≈ 2.05×10⁶。本案例在 Fluent 中对攻角 α = 0°、2°、4°、6°、8°、10°、12°、14°、16° 共 9 个工况做定常绕流计算，提取升力系数 C_L 与阻力系数 C_D，绘制曲线并识别失速攻角，同时用粗/中/细三套网格评估网格分辨率对结果的影响。

攻角通过旋转来流速度矢量实现，而不是旋转翼型几何：u = U∞·cosα，v = U∞·sinα。这样几何、网格与边界命名保持不变，便于批量参数化。升阻力方向约定：升力垂直于来流方向，阻力沿来流方向，Fluent 中通过 Force Report 指定方向矢量 (sinα, -cosα) 与 (cosα, sinα) 分别得到升力与阻力。

计算域、网格与边界条件

计算域采用 C 型拓扑：翼型前缘上游 12c、下游 20c、上下各 12c，避免远场边界对翼型附近流动的干扰。边界条件设置如下：

• 入口（半圆外弧）：velocity-inlet，给定 u、v 分量；
• 出口（下游直边）：pressure-outlet，表压 0 Pa；
• 翼型表面：wall，无滑移；
• 上下远场：symmetry 或 velocity-inlet，本案例用 symmetry 减少反射。

网格采用三套分辨率：粗网格约 1.2 万单元、中网格约 4.8 万单元、细网格约 18 万单元。翼型表面第一层网格高度按 y⁺ ≈ 1 估算，采用 SST k-ω 湍流模型，压力-速度耦合用 Coupled 格式，动量与湍流量均用二阶迎风，残差收敛标准 1×10⁻⁵，同时监测 C_L、C_D 的迭代历史，要求最后 200 步波动小于 0.5%。

参数化计算流程

每个攻角工况的完整流程为：更新入口速度分量 → 初始化（用上一次攻角结果做初值可加速收敛）→ 迭代至收敛 → 读取 C_L、C_D → 写入结果表。9 个工况按攻角从小到大依次计算，小攻角结果作为大攻角初值，可显著减少迭代步数。下表给出中网格下的示例结果（教学示例数据，用于说明曲线形态与失速判断方法，非实测值）：

从表中可见：C_L 在 α ≈ 12° 达到峰值约 1.18，随后在 14° 下降，说明失速攻角约在 12°~13° 之间；C_D 随攻角单调增大，失速后增长更快；升阻比在 α ≈ 8° 达到最大约 64，是巡航效率最优点。

网格分辨率影响与结果判读

用三套网格重复 α = 8° 与 α = 14° 两个工况，观察网格对结果的影响：

判读要点：附着流工况（α=8°）三套网格 C_L 差异小于 6%，中网格已接近网格无关；失速工况（α=14°）粗网格 C_L 偏低约 8%、C_D 偏高约 20%，因为粗网格无法分辨吸力面分离剪切层。因此识别失速攻角必须使用中网格以上分辨率，且失速点附近建议加密翼型后缘与吸力面网格。

结果验证与工程判断

将计算得到的 C_L-α 曲线与 NACA0012 在 Re ≈ 2×10⁶ 下的经典风洞数据对比：小攻角段（α ≤ 8°）斜率约 0.11 /°，与薄翼理论 2π 弧度斜率（约 0.11 /°）一致；失速攻角实验值约 12°~15°，本案例中网格结果落在该区间内，说明 SST k-ω 模型对中等雷诺数翼型失速的预测能力可接受。若需要更精确的失速后行为，应改用非定常求解或 DES/LES 方法。

最终交付物包括：9 个攻角工况的 C_L、C_D 数据表、C_L-α 与 C_D-α 曲线图、失速攻角识别结论（约 12°~13°）、以及网格无关性验证报告（中网格与细网格差异小于 3%）。