空气动力学 — CAE 术语解说
空气动力学是什么
简单地说,空气动力学是什么?Navier-Stokes 方程的求解有什么区别?
空气动力学理论基础
基本概念和控制方程
空气动力学分析中经常出现的"阻力系数"和"升力系数"具体表示什么?不仅仅是力的大小,对吧?
完全正确,它们是"无量纲化的力大小"。具体来说,它们是一个指标,表示物体的形状和迎角在产生阻力和升力方面的效率有多高。例如,在相同的速度、空气密度和翼面积下,阻力系数
控制方程是 Navier-Stokes 方程,但在实际 CAE 中,是否完全求解?还是要简化?
取决于应用。最常见的是时间平均 RANS(Reynolds-Averaged Navier-Stokes)方程。这通过建立湍流的详细时间变动模型,有效地求解平均流场。用方程表示,连续方程和动量方程如下:
"压力系数"
这是物体表面各点压力的无量纲指标。定义为
空气动力学数值计算方法
离散化和求解器设置
在 CFD 中选择"单元"或"元素"的类型(四面体、六面体等)时,空气动力学分析应该以什么为基准决定?
边界层的分辨率是最大的要点。在物体表面,形成"边界层",即流速梯度非常大的区域。为了准确捕捉这一点,通常在表面沿着细长的棱柱(六面体的一种)或多面体单元积层。具体来说,第一层单元的高度应设置为使得无量纲壁面距离
求解器的"压力基"和"密度基"的区别是什么?应该选哪个?
从历史上讲,不可压缩流(Mach 数 Ma < 0.3)适合压力基,可压缩流(Ma > 0.3)适合密度基。但现在的软件(Fluent、STAR-CCM+ 等)模糊了这一界限。实际选择标准如下:压力基对于低速流、自然对流等汽车和建筑物周围的广泛问题具有鲁棒性。另一方面,密度基在超音速流和高度可压缩燃烧分析中,对激波分辨率有更好的倾向。然而,Ansys Fluent 已改进了密度基求解器,可用于所有速度范围。如果不确定,最好从压力基的"Coupled"算法开始。
"稳态分析"和"非稳态分析"何时切换?非稳态需要更多计算时间,所以我想避免。
取决于流场是否有"本质上随时间变化的结构"。例如,汽车尾部后方的涡(弱尾流)即使在稳态分析中也能获得平均模式。但桥梁的"跳振"、航空器的"颤振"、风力发电叶片的"失速现象"是流结构随时间大幅变化的,因此需要非稳态分析。Karman 涡街(圆柱后流等)也是非稳态的。判断的一个帮助是,运行稳态分析后,如果残差继续振荡,或力系数无法收敛到固定值,则可能是非稳态现象。在非稳态分析中,时间步长应根据现象的特征频率(例如从 Strouhal 数估计)设置为 20 倍以上。
空气动力学的实际应用
工作流和检查表
开始空气动力学分析前的"计算域(domain)"大小应该如何确定?太大会增加计算成本,太小结果可能不对。
这是个很好的问题。标准和最佳实践可以参考。在汽车行业,有 SAE J1252 和 J2966 等标准的指导。一般目标是,对于物体的特征长度 L(汽车为全长),
流入面:物体前端上游 3L ~ 5L
流出面:物体后端下游 7L ~ 10L
侧面/顶面/底面:物体中心向外 3L ~ 5L
例如,全长 4.5 m 的轿车,计算域至少需要宽度、高度约 20 m,长度约 50 m。对于建筑物的风荷载分析,可以参考 AIJ(日本建筑学会)的"建筑物荷载指针"。如果域太小,边界的影响会扭曲流,特别是后流的压力分布会变得不准确。
边界条件中"速度入口"和"压力入口"如何区分?"墙面"设置时要注意什么?
"速度入口"用于流入速度已知的情况(模拟风洞实验等)。另一方面,"压力入口"用于流入速度作为结果出现的情况,或模拟外部大气开放条件。实际上,汽车行驶分析通常使用"速度入口",建筑物周围的自然风分析使用指定湍流强度和特征长度的"压力入口"。
"墙面"最重要的是先前提到的边界层网格与壁面函数选择的一致性。使用
如何判断结果是否收敛?残差下降就可以了吗?
残差下降是必要条件,但不充分。最重要的是求解的物理量(阻力系数
1. 主要残差(连续性、动量)至少下降 3~4 个数量级,并在低值处保持平坦。
2. 监测的力或系数值在最后 500~1000 步变化幅度小于 1%。
3. 计算域的进出口质量流量和动量流量平衡良好(误差 1% 以内)。
4. 物理上有意义的值(例如汽车的
如果收敛困难,怀疑网格质量或时间步长(非稳态情况)、求解器的松弛系数。
空气动力学软件比较
Ansys Fluent vs STAR-CCM+ vs OpenFOAM
听说汽车制造商经常使用 STAR-CCM+。与 Fluent 相比,在空气动力学分析中具有决定性优势的是什么?
最大的特点是"一体化"集成环境和强大的自动网格划分功能。STAR-CCM+ 的前处理(CAD 处理、网格生成)、求解、后处理在单一 GUI 中完成,工作流自动化(Java 宏)功能非常强大。在汽车外气流分析中,使用"Surface Remesher"吸收 CAD 的微小缺陷,结合"Polyhedral Mesher"和"Prism Layer Mesher",只需几次点击就能生成复杂车体周围的高质量网格。此外,车辆运动(俯仰、偏航)的耦合分析、电池热流体等多物理场耦合分析也很擅长。Fluent 也有类似功能,但 STAR-CCM+ 在工作流集成和自动化易用性上更受青睐。
那么,Ansys Fluent 在什么场景下有优势?
Fluent 的强项是丰富的物理模型和长期历史积淀的鲁棒性,以及与 Ansys 生态系统的亲和力。特别是在以下领域:
复杂化学反应的燃烧分析(航空发动机等)
多相流(喷雾、气泡)
高级湍流模型(DES、LES 各种衍生模型)
与 Ansys Mechanical 的紧密耦合(流固耦合:FSI)
在这些高级分析中积累了大量实绩和知识。用户定义函数(UDF)的定制性也很高。在航空航天领域和研究机构,Fluent 的市场份额仍然很高。
免费的 OpenFOAM 在实际工作中是否能达到有偿软件的水平?
作为计算引擎,其核心性能非常高,在研究和特定产业领域(特别是风力发电、船舶)被广泛使用。已实现与有偿软件相当或更高的物理模型。但最大的障碍是前处理和后处理。GUI 薄弱(或标准版本不包含),网格生成和设置需要命令行或词典文件,需要时间学习。因此,在实务中引入需要公司内有专业知识的工程师,或购买有偿支持和 GUI 的发行版(ESI-OpenCFD、foam-extend 等)。对于预算有限但需要深度定制的研究开发项目是最佳选择。
空气动力学故障排除
常见错误和对策
计算中出现"负体积"或"负雅可比"错误并停止。这是什么原因,如何解决?
原因是网格质量极差。具体来说:
1. 斜度角过大的单元:特别是六面体网格中,单元过度扭曲。
2. 宽高比过大的单元:单元过于细长。虽然边界层网格特殊,但流域内单元通常超过 1000 时很危险。
3. 网格大小变化过急:相邻单元体积比变化过大会导致梯度计算破裂。
对策是检查网格仪(Fluent Meshing、STAR-CCM+ 的网格生成器等)的质量报告,确保"斜度角"小于 0.9(理想小于 0.8),"宽高比"在流域内部小于 100,通过重新网格划分或网格光滑功能进行调整。
在稳态分析中,阻力系数无法收敛到固定值,而是周期性振荡。这是非稳态现象吗?
很可能是这样。例如,钝头物体(卡车驾驶室后部)的涡脱落本质上是非稳态的。稳态求解器试图找到"时间平均解",但由于物理是非稳态的,解会振荡而无法收敛。首先,估算振荡的频率(如 FFT),判断这个现象对设计有多重要。如果重要,切换到非稳态分析。如果只需要平均值而不重要,可以进一步加强稳态分析的松弛系数(从默认 0.2 降至 0.1 左右),或将 RANS 模型从 SST k-ω 改为更扩散的 k-ε Realizable 模型,以获得伪"平均化"解,但这不是根本解决办法。
风洞实验数据与 CFD 结果,特别是在高迎角处差异很大。可能的原因有什么?
在高迎角,流动分离很大,分离点的再现性是关键。主要原因有:
1. 湍流模型的限制:标准 RANS 模型对分离点和大规模分离域的预测能力较弱。考虑切换到 DES(分离涡模拟)或 LES。
2. 网格分辨率不足:分离和再附近区域的网格太粗。局部细化分离预期边缘周围的网格。
3. 流入湍流条件不匹配:风洞有其固有的湍流强度和尺度。不能随意设置 CFD 的流入条件为"湍流强度 1%、特征长度 0.01 m",需要与风洞实测值(如可能)相符。
4. 模型保真度:CFD 模型是否忽略了支撑杆或地面边界层的影响?风洞实验包含这些干扰,所以 CFD 也需要重现相似条件。首先在分离未发生的低迎角处使实验与 CFD 一致,然后逐步验证是基本原则。