无人机空力设计
无人机空力设计的理论基础
概述
老师，无人机的空力设计与有人机有什么区别？
最大的区别是雷诺数。有人机在$Re \sim 10^7$下飞行，而小型无人机在$Re \sim 10^4$--$10^6$的低雷诺数区域飞行。在这个区域，层流剥离气泡和过渡现象主导性能。
无人机空力设计需要考虑的要点:
- 低Re翼型选择（Eppler、Selig/Donovan翼型等）
- 螺旋桨-机体干涉
- 多旋翼相互干涉
- 阵风响应（机体小，受阵风影响大）
低雷诺数空力
低Re区域的翼型特性与高Re区域质量上不同。
| Re范围 | 流动特征 | 相应无人机 |
|---|---|---|
| $10^4$--$10^5$ | 层流剥离气泡主导，过渡不稳定 | 微型无人机、昆虫型 |
| $10^5$--$10^6$ | 过渡位置决定性能 | 小型固定翼无人机 |
| $10^6$--$10^7$ | 与有人机相似 | 大型MALE/HALE无人机 |
什么是层流剥离气泡？
在低Re下，边界层以层流形式分离，分离的自由剪切层过渡到湍流并重新附着。这个分离-过渡-重新附着的区域称为"层流剥离气泡"。气泡的大小和位置会大大影响升力和阻力。
低Re下最大升力系数会下降。
螺旋桨空力
无人机螺旋桨空力也是CFD的重要研究对象。
螺旋桨的推力系数和效率:
其中$T$是推力，$n$是转速[rps]，$D$是螺旋桨直径，$J = V_\infty/(nD)$是前进比，$C_P$是功率系数。
对于多旋翼，螺旋桨之间的干涉影响有多大？
相邻螺旋桨的尾流（下洗流）干涉会导致悬停效率下降5--15%。当螺旋桨间距小于直径的1.5倍时，干涉变得明显。用CFD评估这种干涉效应对设计高效的机体至关重要。
超低雷诺数世界——与昆虫相同条件下飞行
小型无人机的螺旋桨直径为10~20厘米，工作雷诺数为Re=10,000~100,000，即"超低Re域"。在这个区域，翼型性能会急剧变化，层流剥离气泡频繁出现。实际上昆虫飞行也在同一区域。蜜蜂和蝴蝶的飞行机制研究正被直接应用于无人机翼型设计。从生物体学到的仿生设计已悄然进入现代商用无人机。
无人机空力设计的数值计算方法
低Re翼型的数值解法
要用CFD求解低雷诺数翼型，什么最重要？
过渡模型的选择最为关键。使用完全湍流假设的RANS模型根本无法重现低Re翼型的特性。
| 模型 | 特点 | 对低Re翼型的适用性 |
|---|---|---|
| SST k-omega（完全湍流） | 无过渡 | 不适用。$C_D$过大，$C_{L,max}$不准确 |
| $\gamma$-$Re_\theta$ 过渡模型 | RANS过渡预测 | 优秀。可重现层流剥离气泡 |
| k-kl-omega | 三方程过渡模型 | 优秀。支持低Re |
| LES | 直接计算大尺度涡 | 最高精度但成本高 |
| XFOIL（面板法+BL） | 2D专用。高速 | 最适合初期设计 |
XFOIL现在还经常使用吗？
Mark Drela开发的XFOIL是低Re翼型设计的标准工具。它将面板法和边界层耦合法结合，包括过渡和层流剥离气泡的分析，在几秒内完成。初期筛选用XFOIL比CFD更高效。
螺旋桨CFD
螺旋桨怎样分析？
有三种方法。
| 方法 | 模型化 | 精度 | 成本 |
|---|---|---|---|
| BEM (叶素动量法) | 1D理论 | 中 | 极低 |
| 虚拟圆盘（Actuator Disk） | 用体积力表示螺旋桨 | 中 | 低 |
| 完整叶片分析 | 3D CFD直接求解叶片形状 | 高 | 高 |
虚拟圆盘模型的原理是什么？
在螺旋桨位置设置薄盘区域，根据BEM理论计算的推力和扭矩对应的体积力。无需用网格解析叶片形状，因此对螺旋桨-机体干涉评估很高效。
多旋翼的CFD
多旋翼（四轴飞行器等）的CFD策略:
- 悬停：每个转子用虚拟圆盘或MRF定常分析
- 前进飞行：需非定常分析。捕捉转子周期性变动
- 转子间干涉：上方转子的下洗流对下方有影响（同轴转子的情况）
- 网格规模：四转子完整叶片LES有1-3亿个单元
在STAR-CCM+中用刚体运动设置每个转子时，转子方向设置很关键。为了抵消反扭矩，相邻转子要反向旋转。转子方向错误会产生偏航力矩。
火星直升机Ingenuity的CFD验证秘闻
NASA的火星直升机Ingenuity在密度约为地球1/100的超低密火星大气（0.02 kg/m³）中飞行。地面风洞实验很难重现"火星大气压"，所以CFD成为设计的主力工具。最大的困难是低Re×高马赫数（叶端速度超音速的70%）的组合，远超常规空力CFD的范围。通过包含可压缩性效果的精细CFD和部分真空室实验的结合，完成了设计验证。这个CFD应用案例极具参考价值。
无人机空力设计的实际应用
分析流程
请教一下无人机空力设计的典型CFD工作流程。
以固定翼无人机翼设计为例:
1. 需求规范定义：巡航速度、航续距离、载荷
2. 翼型选择：用XFOIL筛选适合Re的翼型（Selig S1223、Eppler E387等）
3. 翼平面形状设计：纵横比、锥角比、扭转下缘
4. 3D RANS分析：用$\gamma$-$Re_\theta$过渡模型获得$C_L$-$\alpha$、$C_D$-$C_L$
5. 螺旋桨-机体耦合：虚拟圆盘模型评估推力影响
6. 稳定性分析：计算纵向、横向静稳定微系数
7. 风洞/飞行试验验证：将$C_L$/$C_D$与实测对比
低Re翼型的网格
低Re翼型的网格怎样构建？
层流剥离气泡的解析是关键。
- $y^+ < 1$：必须。壁面剪切应力精度直接影响过渡预测
- 棱柱层：30-40层（解析层流→过渡→湍流发展）
- 翼弦方向：300-500点（气泡区聚集划分）
- 翼上表气泡预期位置：单元尺寸≤翼弦的0.1%
- 增长率：1.05-1.1（低Re建议增长率缓慢）
增长率1.1以下，相当细啊。
低Re区，过渡点位置变化1%就使空力特性大变。网格粗细直接影响过渡位置误差，比常规高Re分析对网格品质的要求更严格。
常见失败和对策
| 症状 | 原因 | 对策 |
|---|---|---|
| $C_D$是实验的2倍以上 | 完全湍流假设 | 启用过渡模型 |
| 看不到层流剥离气泡 | 网格不足 | 细化气泡预期位置 |
| 螺旋桨尾流太弱 | 虚拟圆盘参数不对 | 与BEM数据对标 |
| 纵向力矩不符 | 翼型后缘处理 | 使用有限厚后缘（0.5%c） |
| 横风稳定性不准 | 机体完整度不足 | 建模垂直尾翼、鸭翼细部 |
XFOIL应用要点
用XFOIL的技巧？
XFOIL和CFD的分工很重要呢。
XFOIL最适合快速评估2D翼型，但对3D效应（翼端涡、螺旋桨干涉）和超低Re域（<$5 \times 10^4$）需用CFD。按设计阶段合理分工最高效。
配送无人机螺旋桨设计的"风向"效应
物流无人机水平飞行时迎风。悬停优化的螺旋桨前进飞行效率低，前进飞行优化的在悬停时很吃力。CFD需要覆盖从悬停到巡航速度的多工况参数研究。亚马逊的无人机团队投入数千CPU小时研究这个权衡，考虑了固定桨和变桨的组合方案。"找到最优折衷点"是无人机实际设计的核心。
无人机空力设计的软件对比
主要工具
| 工具 | 特点 | 无人机设计的优势 |
|---|---|---|
| XFOIL | 2D翼型。面板法+BL | 低Re翼型高速筛选 |
| Ansys Fluent | 过渡模型全面 | $\gamma$-$Re_\theta$低Re分析 |
| STAR-CCM+ | 重叠网格、6DOF | 多旋翼、螺旋桨耦合 |
| OpenFOAM | 免费、可定制 | 高校、创业公司 |
| AVL（雅典涡格法） | 涡格法。3D升力面 | 翼平面初期设计 |
| VSPAero | OpenVSP的空力求解器 | NASA开发。概念设计 |
AVL和VSPAero是新接触的。
AVL是Mark Drela开发的基于涡格法（VLM）的分析代码。能高速计算翼的3D升力分布和静稳定微系数。OpenVSP是NASA开发的参数化机体形状生成工具，附带的VSPAero能概算空力特性。
Fluent中的低Re翼型分析
推荐设置:
- 求解器：压力基、定常（初期）→非定常（失速附近）
- 湍流模型：SST k-omega+过渡模型（$\gamma$-$Re_\theta$）
- $y^+$：<1。低Re壁面处理
- 入口湍流：$Tu = 0.1$--$1.0%$（匹配风洞湍流强度）
- 空间离散：二阶迎风
- 收敛：残差$10^{-6}$，$C_L$/$C_D$监控稳定
STAR-CCM+中的多旋翼分析
工具选择指南
| 设计阶段 | 推荐工具 | 理由 |
|---|---|---|
| 概念设计 | XFOIL + AVL | 高速、免费 |
| 翼型优化 | Fluent（$\gamma$-$Re_\theta$） | 过渡预测精度 |
| 螺旋桨设计 | STAR-CCM+（重叠） | 转子处理容易 |
| 完整飞行器 | Fluent/STAR-CCM+ | 3D RANS+过渡模型 |
| 飞行动力学耦合 | STAR-CCM+（DFBI） | 6DOF+CFD耦合 |
| 教学、科研 | OpenFOAM | 免费、源码公开 |
按阶段选择工具这样很高效。
无人机开发资源往往有限。先用XFOIL和AVL缩小设计空间，再进入3D CFD，成本和进度都最优。
农业无人机的喷嘴喷雾CFD防止农药漂移
农业无人机喷洒农药时，CFD预测螺旋桨气流如何运送农药液滴。螺旋桨复杂涡的液滴漂移问题"漂移问题"在实际中很严重，导致农药进入相邻农田。通过CFD+粒子追踪（拉格朗日追踪）可预测不同风向和飞行高度的漂移量。这份数据已成为农药喷洒航空法规的依据，展示了CFD在农业行政中的影响。
无人机空力设计的前沿研究
eVTOL/空飞行汽车的空力
空飞行汽车（eVTOL）的空力设计有什么难题？
eVTOL要兼顾垂直起降和巡航，空力设计极其复杂。
主要课题:
- 过渡飞行：从悬停到巡航过程空力特性急变
- 多转子干涉：8-18个转子的相互干涉需CFD评估
- 倾斜转子：转子倾角变化的空力特性制图
- 城市环境：建筑物风中飞行稳定性评估
- 噪声管制：都市区噪声是最大的社会课题
多转子的CFD计算成本很大吧。
用虚拟圆盘模型化所有转子，关键的几个用完整叶片分析切换——多保真度方法更实用。
振翅翼MAV
仿生昆虫和鸟类的摆动型微型无人机（MAV）的空力：
- $Re \sim 10^3$--$10^4$的超低Re区
- 翼大幅变形的流固耦合
- 前缘涡（LEV）：非定常升力的主要机制
- 拍打&摆动机制：翼端接触产生升力增大
定常翼的$C_{L,max} \approx 1.0$大幅超过，非定常效应带来升力增大。
机器学习的空力建模
无人机飞行控制需要CFD精度，但不能承受计算成本。代理模型ML发挥作用。
- 空力表生成：CFD数据训练神经网络，实时预测$C_L$/$C_D$/$C_M$
- 阵风响应预测：RNN或LSTM学习时序空力响应
- 形状优化：VAE（变分自编码器）探索翼型形状的潜在空间
- 强化学习：CFD环境中学习飞行控制策略
无人机空力设计的故障排除
1. 低Re翼型的$C_D$与实验偏离大
症状：CFD的$C_D$是风洞实验的2-3倍
原因：未使用过渡模型（完全湍流假设）
对策:
- 必须启用$\gamma$-$Re_\theta$过渡模型
- 入口湍流强度$Tu$调整到风洞实测值（低乱流风洞：$Tu < 0.1%$）
- OpenFOAM使用kOmegaSSTLM（Langtry-Menter过渡模型）
2. 层流剥离气泡预测不准
层流剥离气泡位置和长度与XFOIL不一致。
对策:
- 网格：在气泡预期分离位置集中细化（翼弦的0.05%以下单元）
- $y^+ < 0.5$为目标（比标准$y^+ < 1$更严格）
- 对标XFOIL的Ncrit设置和CFD入口湍流强度
- $Tu = 0.07%$（XFOIL Ncrit=9相当）是典型低乱流风洞
3. 螺旋桨尾流衰减太快
症状：螺旋桨下游几个直径后尾流消失
原因：RANS的数值耗散使涡衰减
对策:
- 细化尾流区网格（单元尺寸≤叶片弦长的1/10）
- 转换到DDES直接计算尾流涡
- 使用涡约束法补偿涡衰减（部分求解器可用）
4. 多旋翼推力过大
四轴飞行器计算中，转子间干涉没有正确出现。
对策:
- 虚拟圆盘：检查转子间距离的网格解析度
- 完整叶片：检查重叠网格的插值品质
- 转子上下确保充足空间（转子直径的3倍以上）
- 悬停时明确地面效应的有无（近地面vs自由空间）
验证数据
低Re翼型的CFD验证用公开数据:
| 翼型 | Re | 数据源 |
|---|---|---|
| Eppler E387 | $6 \times 10^4$--$5 \times 10^5$ | UIUC低速翼型试验 |
| Selig S1223 | $2 \times 10^5$ | UIUC数据库 |
| SD7003 | $6 \times 10^4$ | AFIT基准（LES/DNS） |
| NACA 0012 | $10^4$--$10^6$ | 众多文献数据 |
UIUC的低速翼型试验数据很有用。
Michael Selig团队在伊利诺伊大学积累的低Re翼型风洞数据是国际标准基准。翼型坐标和风洞数据都免费公开，CFD验证最合适。
多旋翼"单轴失效"问题的CFD分析
多旋翼的安全关键是"单轴失效"应对。一个转子突然停止，机体会剧烈旋转下降，这个动态与剩余转子的复杂气流干涉密切相关。单纯计算无法精确预测。CFD模拟失效直后的瞬时流场，为控制算法提供输入数据的工作在推进。"CFD与控制工程融合"的领域中，流体工程师和控制工程师的协作增加。航空认证机关也开始把CFD-控制耦合数据作为安全证明的根据。
相关主题
了
细
报告