汽车空气动力学模拟
汽车空气动力学理论基础
概述
老师，汽车的空气动力学模拟是为了什么目的呢？
汽车的空气动力学有3个目标。(1)空气阻力$C_D$降低带来的燃油效率改善、(2)升力$C_L$降低带来的高速稳定性保证、(3)风噪降低。
空气阻力与速度的平方成正比。高速行驶时的燃油效率与此直接相关，$C_D$每改进0.01燃油效率就能改善约0.3--0.5%。对电动汽车的续航里程也有重大影响。
支配方程
汽车空气阻力：
其中$A$是前面投影面积（乘用车约2.0--2.5 m^2）。
代表车型的$C_D$值：
| 车型 | $C_D$ | 备注 |
|---|---|---|
| 普通轿车 | 0.28--0.35 | 标准乘用车 |
| Tesla Model S | 0.208 | 2024年量产车最低水平 |
| Mercedes EQS | 0.20 | 量产车全球最低 |
| SUV | 0.35--0.45 | 车高较大不利 |
| 卡车 | 0.6--0.9 | 形状棱角分明 |
$C_D = 0.20$确实相当低啊。
理想的流线型（水滴型）$C_D \approx 0.04$。由于实际车辆设计受到居住空间和法规限制，0.20在量产车中属于极其优秀的数值。
雷诺数与流动特性
乘用车的雷诺数基于车长约为$Re \approx 3 \times 10^6$--$10^7$。处于完全湍流区，边界层转移的影响相对较小。
汽车周围的流动特征：
- 驻点：前格栅附近
- 加速域：引擎盖上表面、车顶
- 分离点：A柱、后车窗末端
- 尾流：巨大的涡结构（阻力的主要来源）
- 底盘：地面效应、轮胎周围的复杂流动
后部形状对阻力影响很大吧。
Ahmed body（汽车空气动力学标准基准）研究表明，后部倾斜角度25°和35°时尾流结构发生巨大变化。25°时形成C柱涡结构，35°时完全分离，$C_D$会不连续地变化。
行驶阻力与燃油效率
行驶阻力的组成：
| 速度 | $F_{aero}$的贡献 |
|---|---|
| 60 km/h | 约30% |
| 100 km/h | 约60% |
| 130 km/h | 约75% |
高速公路速度范围内空气动力学是占主导地位的。
普锐斯的Cd=0.25和"无镜"争议
初代普锐斯的Cd值为0.29，但第三代实现了0.25，属当时量产车顶级水平。开发团队特别讨论的是侧镜。将侧镜改为摄像头能让Cd进一步降低0.004～0.006。但当时日本道路交通法的限制使其无法实现。后来法律修改允许摄像头镜像系统时，工程师们的感受可想而知。CFD虽然能指出"应该这样改进"，但在实务中由于法规和成本限制往往无法实现。
汽车空气动力学数值计算方法
分析方法
汽车空气动力学CFD中使用什么方法呢？
让我们整理一下方法的选择和应用。
| 方法 | 单元数 | 用途 | OEM使用现状 |
|---|---|---|---|
| 定常RANS | 3000万--1亿 | $C_D$/$C_L$的设计评估 | 所有OEM |
| 非定常RANS (URANS) | 5000万--2亿 | 侧镜周围的变动 | 多数OEM |
| DES/DDES | 1亿--5亿 | 尾流、A/C柱涡 | 顶级OEM |
| LBM (PowerFLOW等) | 数亿体素 | 整车非定常分析 | 宝马、福特等 |
| LES | 5亿--10亿+ | 科研用途 | 大学、研究机构 |
宝马使用PowerFLOW是众所周知的。
宝马已使用PowerFLOW（格子玻尔兹曼法）作为量产车开发的主要工具超过20年。与传统N-S求解器相比，它的优势在于网格生成更容易、非定常尾流再现效果更好。
网格策略
整车网格划分：
- 表面网格：车身表面3--5mm的三角形网格
- 棱柱层：$y^+ \approx 30$--100（壁面函数）或$y^+ < 1$（低雷诺数壁处理）
- 轮毂旋转：MRF / 滑移网格
- 地面移动：与车速相同的速度
- 散热器：多孔体模型（压力损失系数从实测数据获得）
- 发动机舱：内部流路建模（冷却系统压力损失）
- 远方边界：车长的5倍以上
有时使用壁面函数，有时不使用。
量产车开发受计算时间限制，通常使用壁面函数（$y^+ \approx 30$--100）。虽然$C_D$的绝对精度不如$y^+ < 1$，但对设计变更的差分评估（$\Delta C_D$）已经足够实用。
旋转轮毂与接地面
轮毂占总阻力的约25--30%，是重要元素。
| 模型化要素 | 对$C_D$的影响 | 备注 |
|---|---|---|
| 轮毂旋转 | $\Delta C_D \approx +0.015$ | 旋转与否影响很大 |
| 轮胎变形 | $\Delta C_D \approx +0.005$ | 接地面形状的影响 |
| 制动冷却导管 | $\Delta C_D \approx +0.003$ | 内部流动的影响 |
| 轮圈设计 | $\Delta C_D = -0.005$--$+0.010$ | 取决于开口率 |
仅轮毂对$C_D$就有0.02以上的影响。
最近的电动汽车趋势是装配气动轮毂盖来降低$C_D$。Tesla Model 3的气动盖能降低$C_D$0.008。这种微小的$\Delta C_D$评估正是CFD的用武之地。
收敛判定
Ahmed Body为什么成为全球标准基准
"Ahmed Body"是1984年Ahmed等人公开风洞实验数据的一个简单箱型模型，已成为汽车空气动力学CFD验证的常用基准。后端倾斜角为25°时产生强纵向涡，35°时Cd值突然下降。这种"倾斜角敏感性"是否能被CFD再现成为衡量模型能力的试金石。Fluent、OpenFOAM、SUPERFLOW等各种工具都用这个案例验证，汽车空气动力学工程师的惯例就是先用Ahmed Body验证CFD设置，然后再应用到实车。
汽车空气动力学实务应用
汽车OEM的开发流程
请介绍一下汽车厂商的CFD开发流程。
1. 概念设计：从设计师草图创建粗糙的CFD模型
2. 外形优化：车顶高度、后悬长度、底盘形状的优化
3. 细节设计：格栅、侧镜、轮毂、扰流板的研究
4. 风洞验证：1/4缩放或全尺寸风洞中$C_D$的实测
5. CFD-风洞相关性：确认相关性，必要时调整CFD设置
6. 量产规格确认：最终形状的$C_D$/$C_L$确认
一台车运行多少个CFD案例呢？
量产车一个车型的开发会运行数百到数千个CFD案例。有些OEM每周运行50--100个案例的速度。自动化是必需的。
底盘设计
底盘占$C_D$的20--25%。仅平坦化就能达到$\Delta C_D = -0.01$--$-0.03$的效果。
| 底盘要素 | 对$C_D$的影响 |
|---|---|
| 平整底盘罩 | -0.010--0.030 |
| 后部扩散器 | -0.005--0.015 |
| 前部扰流板（挡风板） | -0.010--0.020 |
| 侧面裙板 | -0.003--0.008 |
电动汽车没有发动机，更容易使底盘平坦吧。
完全同意。电动汽车的电池组装在底部，底盘原本就接近平坦。这正是电动汽车$C_D$较低的重要原因之一。
常见失败与对策
| 症状 | 原因 | 对策 |
|---|---|---|
| $C_D$与风洞偏差0.02以上 | 轮毂旋转条件不完备 | 检查MRF/滑移网格 |
| $C_L$前后平衡不符 | 散热器模型不完备 | 从实测取得多孔体参数 |
| A柱涡无法再现 | 定常RANS的局限 | 转移到DDES/DDES |
| 风切噪音预测不准确 | 网格分辨率不足 | 侧镜周围提升至5000万单元以上 |
| 横风稳定性不准确 | 偏航角解析网格问题 | 为各偏航角创建专用网格 |
风洞试验的相关性
CFD和风洞的$C_D$一般能精确到什么程度？
一般的相关精度：
- RANS：$\Delta C_D \approx \pm 0.010$--$0.020$（绝对值）
- DES/DDES：$\Delta C_D \approx \pm 0.005$--$0.010$
- LBM (PowerFLOW)：$\Delta C_D \approx \pm 0.003$--$0.008$
不过绝对值的精度不如设计变更差分$\Delta C_D$的精度在实务中重要。差分评估通常能达到$\Delta(\Delta C_D) \approx \pm 0.002$--$0.005$的精度。
如何在CFD中再现实车行驶时的"地面效应"
实际车辆在路面上行驶，地面靠近时下表面流速增加，升力减少，这叫"地面效应"。风洞试验用皮带式地面（Moving Ground）再现路面运动，CFD中也要将底面边界设置为动态壁面才能再现。若不这样设置，升力系数会比实测值高10～20%。"为什么CFD结果与风洞不符"的困惑，检查地面边界条件往往能解决。
汽车空气动力学软件对比
主要工具
| 工具 | 使用OEM | 特点 |
|---|---|---|
| Ansys Fluent | 丰田、通用、大众 | 通用CFD。RANS/DES支持 |
| STAR-CCM+ | 宝马 (RANS)、现代 | 自动网格、重叠网格 |
| PowerFLOW (达索) | 宝马、福特、PSA | LBM法。非定常性强 |
| OpenFOAM | 沃尔沃、大学 | 免费。定制性强 |
| Exa/Simulia XFlow | 雷诺 | LBM。无网格 |
| iconCFD | 部分OEM | 基于OpenFOAM。汽车专用 |
丰田使用Fluent。
丰田长年使用Fluent，构建了独自的网格自动化管线。最近DES/DDES的活用也在推进。
PowerFLOW vs N-S求解器
格子玻尔兹曼法的PowerFLOW与传统的N-S求解器有什么区别？
宝马选择PowerFLOW的理由是网格生成容易实现自动化，非定常尾流（特别是后形状评估）精度高。而N-S求解器在物理模型自由度方面更高（散热器多孔体、风扇旋转等）。
工具选择指南
| 用途 | 推荐工具 | 原因 |
|---|---|---|
| 量产车$C_D$开发 | PowerFLOW / Fluent | 实績、精度 |
| 电动汽车续航里程优化 | STAR-CCM+ | 集成优化工作流 |
| 空气动力学噪声 | PowerFLOW | LBM+FW-H噪声解析 |
| 研究、教育 | OpenFOAM | 免费、DrivAer基准支持 |
| 供应商（零部件） | Fluent / STAR-CCM+ | 通用性 |
电动汽车时代到来，空气动力学工具选型也在变化吗？
电动汽车时代续航里程是最大的卖点，$C_D$的0.001单位改善很重要。PowerFLOW、XFlow等LBM型工具的需求正在上升，这反映了对非定常精度要求的提高。
格子玻尔兹曼法为何在汽车业开始流行
传统的RANS-CFD在定常计算方面很强，但对汽车的A柱附近分离涡和行驶中非定常尾流的处理就不够了。格子玻尔兹曼法（LBM）本质上是非定常的，这种非定常特性被自然处理，而且大规模并行计算效率高。BMW、大众、通用等纷纷采纳正是这个2010年代的变化，打破了之前Fluent一家独大的局面。"工具选择涉及供应商锁定风险"这个话题在汽车业是真实发生的历史。
汽车空气动力学前沿研究
空气动力学噪声（航空声学）
汽车的风切音能用CFD预测吗？
可以用Ffowcs Williams-Hawkings（FW-H）公式从车体表面压力脉动推算远场音压。
主要风切音源：
- 侧镜：最大音源。涡流放出频率$f \approx St \cdot V/D$（$St \approx 0.2$）
- A柱-车窗隙间：隙间风引起的高频音
- 轮毂舱：湍流噪声
- 雨刷/天线：圆柱周围的埃奥罗斯音
改成摄像头镜像（数字镜像）风切音就会减少吧。
完全同意。用小型摄像头取代传统侧镜，$C_D$能降低0.01--0.02，风切音也大幅降低。奥迪e-tron和雷克萨斯ES已采用数字镜像。
形状优化
汽车空气动力学形状优化方法：
- 参数化优化：以设计参数（车顶高度、后扰流板角度等）为变量进行DOE+响应面法
- 随伴法：单次追加计算获得全表面灵敏度。可得$\partial C_D/\partial x_i$的全节点值
- 拓扑优化：自动确定计算域内最优物体配置（研究阶段）
随伴法效率看起来很高啊。
Fluent/STAR-CCM+都配备了随伴求解器。即使有数百个设计变量，追加计算也只需1次，能高效地获得表面形状的灵敏度图。
电动汽车与自动驾驶的空气动力学
电动/自动驾驶时代的空气动力学趋势：
- 无前格栅：无需冷却进气，前端设计更流畅
- 隐藏式车门把手：消除突起，$C_D$降低0.003
- 主动空气动力学：根据速度自动控制格栅快门、扰流板
- 激光雷达空气动力学影响：车顶激光雷达单元$\Delta C_D \approx +0.005$
- 编队行驶：跟随车辆阻力降低30--40%
DrivAer基准
慕尼黑工业大学（TUM）公开的DrivAer模型是汽车空气动力学CFD的标准基准。
- Fastback/Notchback/Estateback三种形态
- 开/闭底盘
- 有/无轮毂
- 公开风洞实验数据
- 可作为OpenFOAM教程案例使用
用DrivAer验证CFD设置后再移到实车，是最佳实践吧。
电动汽车时代汽车空气动力学"优先级"改变
汽油车时代"发动机冷却所需进气"与"空气阻力降低"相互制约。为了冷却，必须开大格栅，这会加重空气阻力。而电动汽车冷却风量需求陡降，格栅可以大胆关闭，空气动力学优先级可以重新定义。特斯拉Model 3的Cd=0.23正是因为这一点，设计自由度与汽油车根本不同。电动化对空气动力学约束条件的重新定义，是汽车空气动力学工程师们讨论的重点。
汽车空气动力学故障排除
1. $C_D$与风洞不符
症状：CFD的$C_D$与风洞偏差0.02以上
检查清单：
- 轮毂是否在旋转（MRF或滑移网格）
- 地面是否为移动墙
- 散热器多孔体模型是否正确（使用实测压力损失系数）
- 发动机舱/冷却系统流路是否建模
- 风洞的支架/支持装置是否考虑
- 轮胎接地面是否再现
2. 尾流非对称
对称车体为什么尾流偏向一侧？
原因：钝头物体尾流可出现双稳定（bi-stable）状态。定常RANS会锁定在一侧。
对策：
- 用对称面条件进行半模计算（无法捕捉非对称尾流但保证$C_D$精度）
- 转换到URANS或DDES进行时间平均
- 严格确保网格的对称性
3. 轮毂舱附近发散
症状：旋转轮毂与挡板间计算发散
对策：
- 轮毂与挡板间隙至少保留5个网格单元
- 轮胎接地面附近网格质量检查（非正交度 < 70°）
- 确认MRF区域边界与挡板内表面无干涉
- 初期条件从低速开始，逐步加速
4. 冷却系统压力损失不准确
症状：散热器通过后温度/压力与实车偏差
对策：
- 从实测数据输入散热器压力损失特性作为风速函数
- 冷凝器、中冷器也分别设置为多孔体
- 风扇旋转：MRF（简化）或SRF（旋转坐标系定常）
- 发动机舱细节（管路、配管）适当简化
Ahmed Body验证要点
Ahmed Body如何验证？
Ahmed Body（倾斜角25°和35°）是必需的验证案例。
检查要点：
- 25°：C柱涡结构是否被再现（通过压力场可视化确认）
- 35°：完全分离导致尾流扩大的模式是否被再现
- 25°→35°是否能捕捉$C_D$的不连续变化
- 底面（后面）$C_p$分布是否与实验一致
- DrivAer模型的$C_D$是否在公开数据的$\pm 0.005$以内
先用基准验证手法，再应用到实车是很重要的呢。
完全同意。用Ahmed Body和DrivAer验证CFD设置的合理性是业界最佳实践。跳过这一步，实车出现问题时排查会很困难。
"与风洞完全不符"原因的判断步骤
汽车CFD中"与风洞不符"的咨询在实务中很常见。首先要检查的是轮毂旋转。轮毂静止建模会导致轮胎周围流动形态与实测差很大，影响整个尾流。其次检查车内冷却流有无，格栅进气的冷却风对$C_D$的贡献最多0.01。"同一形状却出现差异"时，要一项一项地检查这样的边界条件细节，这是CFD故障排除的铁律。
更详细
错误