3.2 风扇模型处理

风扇子程序对模型要求不高，只需要大致的轮廓尺寸即可，即风扇的大致直径和厚度，但风扇域的网格却要经过一定的处理。风扇域网格必须使用六面体网格，且不能是O 型六面体网格，因为此子程序在计算O 型网格的通量时会发生方向错乱。风扇域的六面体网格可以在ICEM-CFD 里做出来， 也可以在STAR-CD 里基于四边形面网格直接拉伸。 整个风扇域需要分解成三个域，即入口域、旋转域和出口域。入口域和出口域至少需要有一层六面体体网格，主要用于计算速度和压力。在整个分析模型中，风扇域是六面体网格，其余部分为 四面体网格，因此风扇入口面和出口面的网格需要与外围四面体网格做couple 处理，以保证体网格的连通性。做好后的风扇六面体网格及前后相接的四面体网格模型如图2 所示。

图2 风扇模型网格处理

3.3 数学模型

汽车流场空气流动是低马赫、高雷诺数、三维不可压缩粘性流动，因此采用标准高雷诺数k-e 两方程湍流模型建立数学模型。此湍流模型包含的非线性构成关系能够较好地处理流场中湍流的各向异性等困难问题。流动控制方程用有限体积法进行离散，用中心差分法求解。方程的对流项采用中心差分格式与迎风差分格式的混合格式，二阶中心差分格式处理扩散项，以提高计算精度。近壁面采用非平衡壁面函数处理。

3.4 边界条件

计算工况为40km/h、坡度12%，冷凝器和散热器由于结构复杂均采用多孔介质模型模拟，其主要参数粘性阻抗和惯性阻抗可由试验数据拟合获得。风扇采用风扇子程序模拟。

4 风扇子程序编译

整个风扇子程序分成两大部分，即用于加载动力源项的“sormom”和用于计算及后处理的“posdat”。在“sormom”程序段里的主要工作就是定义风扇旋转域体网格，然后赋上初始速度参考值。这个初始值可根据给定的风扇P-Q 曲线及风扇的实际运行环境进行调整。根据实际的模型参数调整好后，接下来就可以在STAR-CD 的源项里启动该子程序。

“posdat”是整个风扇子程序的核心，在这个程序里有一些重要的参数需要设定。需要单独定义的有风扇出口域、旋转域和出口域体网格号、空气密度、风扇出口面积、通量积分方向以及风扇P-Q 曲线拟合公式等。dppq 方程是风扇子程序的关键部分之一，其方程拟合需要风扇P-Q 曲线以及风扇出口面积等数据。横坐标为速度（m/s），纵坐标为压降（Pa）。方程拟合方式为二次、3 阶多项式拟合。如图3 所示。

图3 dppq 方程二次多项式拟合

5 计算结果

5.1 原风扇及布置形式分析结果

对原车稳态冷却流场进行详细分析，结果显示：从进气格栅流入的冷却空气有相当大的一部分从发动机舱盖与水箱横梁之间、水箱横梁与散热器之间等的各处缝隙及侧面流入发动机后舱。冷却系统的布置形式不佳和风扇能力不足导致各热交换器的通风量很低，无法满足设计目标值。

5.2 新风扇及布置优化后分析结果

此冷却系统包括中冷器、冷凝器、散热器和冷却风扇等。由于舱内拥挤，布置空间有限，前期采用了热交换器同轴布置形式，造成热交换器的通风阻力很大。经过首轮冷却性能分析后发现问题相当严重。

为了增加冷却系统的进风量，提高冷却性能，更换了能力更强的风扇，并且冷却系统的布置形式也进行了优化（见图4）。

图4 冷却模块布置更改示意图

如图5所示，计算结果显示新风扇的抽风能力和优化后的布置形式效果非常好，性能明显优于之前的设计方式。各热交换器的迎风风速有明显提高，整个冷却系统的通风量显著增加，通风能力得到较大幅度的增强。

图5 优化后发动机舱流速及热交换器迎风风

表1 各关键部件进气量对比（㎡／s）

6 结论

1) STAR-CD风扇子程序能够在设计前期且缺少风扇数模的情况下进行整车冷却性能分析。结果表明，通过拟合输入不同性能的风扇P-Q曲线，能够优选出合适型号的风扇，能够对冷却系统的匹配设计提供有力地支持。

2) 利用CFD技术，能够在汽车设计的初期预测整车性能，并且能够快速便捷地对不同方案进行优化分析，从中选出最优的方案进行产品设计，避免开发后期出现难以解决的问题，从而加快开发速度，节省开发费用。