汽车空气动力学仿真
理论与物理
概述
老师，汽车空气动力学仿真的目的是什么？
汽车空气动力学有三大目标。(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显示“这样做会更好”，但受法规或量产成本限制而无法实现的情况，在实际工作中是家常便饭。后来法规修订允许摄像头后视镜系统时，工程师们会是怎样的心情呢？
各项的物理意义
- 时间项 $\partial(\rho\phi)/\partial t$：想象一下拧开水龙头的瞬间。最初水流会不稳定地喷溅，过一会儿才会变成稳定的水流，对吧？描述这个“变化过程中”的就是时间项。心脏搏动导致血流脉动，发动机气门每次开闭导致流动变化，这些都是非定常现象。那么定常分析是什么？就是“只观察经过足够长时间、流动稳定之后的状态”——也就是将此项设为零。计算成本会大幅下降，因此先用定常求解是CFD的基本策略。
- 对流项 $\nabla \cdot (\rho \mathbf{u} \phi)$：把落叶扔进河里会怎样？它会随水流被带到下游，对吧？这就是“对流”——流体运动搬运物体的效果。暖风的暖空气能到达房间另一端，也是因为空气这个“搬运工”通过对流输送热量。这里有趣的是——这项包含“速度×速度”，因此是非线性的。也就是说，流速加快时此项会急剧增强，变得难以控制。这就是湍流的根本原因。常见的误解：“对流和传导差不多”→ 完全不一样！对流是流动搬运，传导是分子传递。效率有天壤之别。
- 扩散项 $\nabla \cdot (\Gamma \nabla \phi)$：有过把牛奶倒入咖啡后放置的经历吗？即使不搅拌，过一会儿也会自然混合。那就是分子扩散。那么下一个问题——蜂蜜和水，哪个更容易流动？当然是水，对吧？因为蜂蜜的粘度（$\mu$）高，所以不易流动。粘度越大，扩散项越强，流体的运动就越“粘稠”。雷诺数小的流动（缓慢、粘稠）中扩散占主导。相反，雷诺数大的流动中对流占压倒性优势，扩散则成为配角。
- 压力项 $-\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使用哪些方法？
我们来整理一下方法的选择和区别。
| 方法 | 网格数 | 用途 | OEM使用情况 |
|---|---|---|---|
| 定常RANS | 3000万--1亿 | $C_D$/$C_L$的设计评估 | 所有OEM |
| 非定常RANS (URANS) | 5000万--2亿 | 后视镜周围的波动 | 多数OEM |
| DES/DDES | 1亿--5亿 | 尾流、A/C柱涡 | 顶级OEM |
| LBM (PowerFLOW等) | 数亿体素 | 整车非定常分析 | BMW, Ford等 |
| LES | 5亿--10亿+ | 研究用途 | 大学·研究机构 |
BMW使用PowerFLOW很有名呢。
BMW将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%的重要要素。
| 建模要素 | 效果 | 备注 |
|---|---|---|
| 车轮旋转 | $\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。CFD对于评估这种细微的$\Delta C_D$是不可或缺的。
收敛判定
Ahmed Body成为世界标准基准模型的原因
汽车空气动力学CFD验证中常用的“Ahmed Body”，是Ahmed等人在1984年公开风洞实验数据的简单箱型模型。其后端倾斜角设为25°时会产生强烈的纵向涡，35°时Cd会突然下降。CFD能否再现这种“倾斜角敏感性”成为了衡量模型能力的试金石。Fluent, OpenFOAM, SUPERFLOW等各种工具都用这个案例进行验证，先用Ahmed Body确认自己的CFD设置已成为汽车空气动力学工程师的习惯。