翼型与机翼气动分析
理论与物理
概述
老师，飞机机翼周围的空气动力学分析，用CFD具体是怎么做的呢？
这是预测翼型（Airfoil）的升力/阻力特性，评估失速行为和高升力装置效果的解析。是支撑飞机设计根基的技术。
从NACA翼型系列开始，到现代的超临界翼型和自然层流翼型设计，CFD都不可或缺。CFD用于拓展仅靠风洞试验无法完全探索的设计空间。
只用风洞不行吗？
风洞试验每个工况的成本高达数百万日元级别。先用CFD筛选设计候选方案，再进入风洞测试，这是现代标准的工作流程。
控制方程
请告诉我描述机翼周围流动的方程。
可压缩Navier-Stokes方程是基础。用连续方程、动量方程、能量方程这三组方程来描述。
升力系数和阻力系数是这样定义的。
其中 $L$ 是升力，$D$ 是阻力，$\rho_\infty$ 是自由流密度，$V_\infty$ 是自由流速度，$S$ 是机翼面积。
雷诺数大概在什么范围呢？
基于翼弦长的雷诺数，在客机巡航时大约是 $Re_c \approx 2 \times 10^7$ 量级。马赫数主要在跨音速区域的 $M \approx 0.78$--$0.85$。
在跨音速区域，机翼上表面会产生局部超音速区域并形成激波。这种激波与边界层的相互作用会引发抖振现象，成为失速的诱因。
原来如此。准确捕捉激波的位置很重要呢。
湍流模型的选择
机翼周围的解析通常使用哪种湍流模型呢？
根据用途区分使用。预测转捩的情况和假设完全湍流的情况，选择会不同。
| 模型 | 特点 | 对机翼解析的适用性 |
|---|---|---|
| Spalart-Allmaras (SA) | 单方程模型。航空领域广泛使用 | 巡航条件下表现良好。接近失速时稍弱 |
| SST k-omega | 近壁k-omega与远场k-epsilon的混合 | 对逆压梯度、分离预测能力强 |
| gamma-Re_theta 转捩模型 | 与SA/SST结合预测自然转捩 | 自然层流翼设计必备 |
| DDES/IDDES | RANS+LES的混合 | 大规模分离、抖振解析 |
Spalart-Allmaras模型在航空领域非常主流呢。
是的。SA模型原本就是NASA为翼型解析开发的模型。波音、空客两家公司都广泛使用。不过，对于接近失速时的大规模分离，就需要SST k-omega或DDES了。
翼型的气动特性
具体会得到什么样的数值结果呢？
我们来整理一下典型翼型的气动参数。
| 翼型 | 用途 | 设计$C_L$ | $C_{L,max}$ | 失速迎角 |
|---|---|---|---|---|
| NACA 0012 | 基准测试 | 0 (对称) | 约1.5 | 约16度 |
| NACA 23012 | 通用 | 0.3 | 约1.8 | 约18度 |
| RAE 2822 | 跨音速基准测试 | 0.74 (M=0.73) | -- | -- |
| SC(2)-0710 | 超临界 | 0.7 (M=0.78) | -- | -- |
RAE 2822因为风洞数据是公开的，所以经常用于CFD的验证，对吧。
没错。Case 9（$M=0.73$, $\alpha=2.79°$, $Re=6.5 \times 10^6$）是业界标准的基准测试。用CFD计算翼面上表面的压力分布和激波位置，与风洞数据进行比较。
实务注意事项
机翼解析中需要特别注意哪些要点？
最重要的是边界层分辨率。需要将壁面第一层的$y^+$值设定在1以下，并确保边界层内有足够的棱柱层。
机翼解析的关键在于边界层的处理啊。我明白了。
是的。阻力的预测精度直接取决于边界层网格的质量。有时会要求达到 $\Delta C_D = 0.0001$（1 count）的精度，所以对网格需要格外小心。
NACA翼型的“编号”有意义
NACA的4位数翼型，比如NACA 2412，这些数字可不是随机的。第一个“2”表示最大弯度是弦长的2%，“4”表示该位置距离前缘40%，“12”表示最大厚度是12%。这是1930年代NACA系统地对数百种翼型进行风洞测量后设定的命名规则。多亏了这个规则，看到翼型编号就能在脑海中浮现出“哦，这是较薄、适合后掠翼的形状”。在CAE中建模时，了解这个体系，可以在导入坐标数据前检查形状的合理性，会默默地帮上大忙。
各项的物理意义
- 时间项 $\partial(\rho\phi)/\partial t$：想象一下拧开水龙头的瞬间。一开始水流会不稳定地喷溅，过一会儿才变成稳定的水流，对吧？描述这个“变化过程中”的就是时间项。心脏搏动导致血流脉动，发动机阀门每次开闭引起流动变化，这些都是非定常现象。那么定常解析是什么？就是只看“经过足够时间、流动稳定之后”——也就是将此项设为零。计算成本会大幅下降，因此先用定常求解是CFD的基本策略。
- 对流项 $\nabla \cdot (\rho \mathbf{u} \phi)$：把落叶丢进河里会怎样？会被水流带着往下游漂去。这就是“对流”——流体的运动搬运物体的效果。暖气的热风能到达房间另一端，也是因为空气这个“搬运工”通过对流输送热量。这里有趣的是——这项包含“速度×速度”，因此是非线性的。也就是说，流速变快，这项会急剧增强，变得难以控制。这就是湍流的根本原因。常见的误解：“对流和传导差不多吧？”→ 完全不一样！对流是流动搬运，传导是分子传递。效率有天壤之别。
- 扩散项 $\nabla \cdot (\Gamma \nabla \phi)$：有过在咖啡里倒入牛奶后放置不管的经历吗？即使不搅拌，过一会儿也会自然混合。那就是分子扩散。那么下一个问题——蜂蜜和水，哪个更容易流动？当然是水。因为蜂蜜的粘度（$\mu$）高，所以难流动。粘度越大，扩散项越强，流体的运动就变得“粘稠”。雷诺数小的流动（缓慢、粘稠）中扩散占主导。相反，Re数大的流动中，对流占压倒性优势，扩散则成为配角。
- 压力项 $-\nabla p$：推注射器的活塞，液体就会从针头有力地射出，对吧？为什么？因为活塞侧压力高，针头侧压力低——这个压力差就是推动流体的力。大坝放水也是同样原理。天气图上等压线密集的地方会怎样？没错，会刮强风。“有压力差的地方就会产生流动”——这就是纳维-斯托克斯方程压力项的物理意义。这里容易误解的点是：CFD中的“压力”通常指表压而非绝对压力。切换到可压缩解析时结果突然出错，原因可能就是混淆了绝对压力/表压。
- 源项 $S_\phi$：受热的空气会上升——为什么？因为比周围空气轻（密度低），被浮力推上去了。这个浮力作为源项添加到方程中。此外，燃气灶的火焰产生化学反应热，工厂的电磁泵对金属熔液施加洛伦兹力……这些都是“从外部向流体注入能量或力”的作用，都用源项来表示。忘记源项会怎样？自然对流的解析中如果忘记加入浮力，流体就完全不动——冬天房间里开了暖气，暖空气却不往上走，得到这种物理上不可能的结果。
假设条件与适用范围
- 连续介质假设：克努森数 Kn < 0.01（分子平均自由程 ≪ 特征长度）时成立
- 牛顿流体假设：剪切应力与应变速率呈线性关系（非牛顿流体需要粘度模型）
- 不可压缩假设（Ma < 0.3 时）：将密度视为常数处理。马赫数0.3以上需考虑压缩性效应
- Boussinesq近似（自然对流）：仅在浮力项中考虑密度变化，其他项使用恒定密度
- 不适用的情形：稀薄气体（Kn > 0.1）、超音速/高超音速流动（需要激波捕捉）、自由表面流动（需要VOF/Level Set等方法）
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 速度 $u$ | m/s | 入口条件从体积流量换算时，注意截面积单位 |
| 压力 $p$ | Pa | 区分表压与绝对压力。可压缩解析使用绝对压力 |
| 密度 $\rho$ | kg/m³ | 空气: 约1.225 kg/m³@20°C，水: 约998 kg/m³@20°C |
| 粘性系数 $\mu$ | Pa·s | 注意与运动粘性系数 $\nu = \mu/\rho$ [m²/s] 混淆 |
| 雷诺数 $Re$ | 无量纲 | $Re = \rho u L / \mu$。层流/湍流转捩的判定指标 |
| CFL数 | 无量纲 | $CFL = u \Delta t / \Delta x$。直接关系到时间步长的稳定性 |
数值解法与实现
空间离散化
用CFD求解机翼周围的流动时，具体使用什么数值方法呢？
主流是有限体积法（FVM）。采用单元中心型格式，将控制方程在每个单元的体积上进行积分来离散化。
对流项的离散化是精度的关键。必须使用二阶精度以上的格式，具体有以下选择。
| 格式 | 精度 | 特点 | 适用场景 |
|---|---|---|---|
| 二阶中心差分 | 二阶 | 数值耗散少 | LES/DES |
| 二阶迎风差分 | 二阶 | 稳定性高 | RANS定常解析 |
| MUSCL (van Leer) | 二阶TVD | 适合激波捕捉 | 跨音速/超音速 |
| Roe近似黎曼求解器 | 二阶 | 激波高分辨率 | 跨音速翼型 |
跨音速翼型因为有激波，所以会用到Roe格式之类的呢。
相关主题
なった
詳しく
報告