２０２０年第２期
车辆与动力技术Ｖｅｈｉｃｌｅ＆ＰｏｗｅｒＴｅｃｈｎｏｌｏｇｙ
总第１５８期
文章编号：１００９－４６８７（２０２０）０２－００３１－０５
收稿日期：２０２０－０３－３０
作者简介：刘继伟（１９９３－），男，硕士，研究方向为新能源汽车热管理．
基于ＣＦＤ的新能源汽车冷却风扇气动性能仿真分析
刘继伟，　龙海峰，　席忠民，　许俊海，　何凯欣，　何重光，　梁佳佳
（广汽新能源汽车有限公司，广州５１１４３４）
摘　要：以某新能源汽车的７叶片的冷却风扇为研究模型，通过ＳＴＡＲＣＣＭ＋软件中Ｒｅａｌｉｚａｂｌｅｋ ε湍流模型对
其进行定常三维数值计算．首先进行了网格数量的无关性验证；然后通过试验验证了数值计算模型的准确性，并对冷却风扇内部流场压力与速度分布进行了分析；最后分析了叶片个数参数对冷却风扇气动性能的影响．结果表明：相同转速的工况下，当冷却风扇静压相同时，随着叶片个数增多，其产生的流量越大．在冷却风扇的静压效率方面，在风扇静压１７０ ２００Ｐａ左右时，９叶片风扇静压效率最高．在其他静压区间，当叶片数为７、８时，风扇静压效率要高于９叶片风扇．研究可以为新能源汽车冷却风扇气动性能优化提供依据．关键词：冷却风扇；气动性能；数值仿真；叶片个数中图分类号：ＴＨ４２；Ｕ４６２    　文献标识码：Ａ
ＳｉｍｕｌａｔｉｏｎＡｎａｌｙｓｉｓｏｆＮｅｗＥｎｅｒｇｙＶｅｈｉｃｌｅＣｏｏｌｉｎｇＦａｎｓ
ＡｅｒｏｄｙｎａｍｉｃＰｅｒｆｏｒｍａｎｃｅＢａｓｅｄｏｎＣＦＤ
ＬＩＵＪｉｗｅｉ，　ＬＯＮＧＨａｉｆｅｎｇ，　ＸＩＺｈｏｎｇｍｉｎ，　ＸＵＪｕｎｈａｉ，　ＨＥＫａｉｘｉｎ，
ＨＥＣｈｏｎｇｇｕａｎｇ，　ＬＩＡＮＧＪｉａｊｉａ
（ＴｅｃｈｎｏｌｏｇｙＣｅｎｔｅｒ，ＧＡＣＮｅｗＥｎｅｒｇｙＡｕｔｏｍｏｂｉｌｅＣｏ ，Ｌｔ
ｄ，Ｇｕａｎｇｚｈｏｕ５１１４３４，Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：Ｉｎｔｈｉｓｐａｐｅｒ，ｔｈｅ７ ｂｌａｄｅｃｏｏｌｉｎｇｆａｎｏｆａｎｅｗｅｎｅｒｇｙｖｅｈｉｃｌｅｉｓｔａｋｅｎａｓｔｈｅｒｅｓｅａｒｃｈｍｏｄｅｌ
Ｔｈｅｓｔｅａｄｙｔｈｒｅｅ ｄｉｍｅｎｓｉｏｎａｌｎｕｍｅｒｉｃａｌｃａｌｃｕｌａｔｉｏｎｉｓｃａｒｒｉｅｄｏｕｔｂｙｕｓｉｎｇｔｈｅＲｅａｌｉｚａｂｌｅｋ εｔ
ｕｒｂｕｌｅｎｃｅｍｏｄｅｌｉｎＳＴＡＲＣＣＭ＋ｓｏｆｔｗａｒｅ Ｆｉｒｓｔｌｙ，ｔｈｅｎｕｍｂｅｒｏｆｇｒｉｄ ｉｎｄｅｐｅｎｄｅｎｔｖｅｒｉｆｉｃａｔｉｏｎｉｓｃａｒｒｉｅｄｏｕｔ Ｔｈｅｎｔｈｅａｃｃｕｒａｃｙｏｆｔｈｅｎｕｍｅｒｉｃａｌｍｏｄｅｌｉｓｖｅｒｉｆｉｅｄｂｙｅｘｐｅｒｉｍｅｎｔｓ Ｔｈｅｄｉｓｔｒｉｂｕｔｉｏｎｏｆｐｒｅｓｓｕｒｅａｎｄｖｅｌｏｃｉｔｙｉｎｔｈｅｃｏｏｌｉｎｇｆａｎｉｓａｎａｌｙｚｅｄ Ｆｉｎａｌｌｙ，ｔｈｅｉｎｆｌｕｅｎｃｅｏｆｔｈｅｎｕｍｂｅｒｏｆｂｌａｄｅｓｏｎｔｈｅａｅｒｏｄｙｎａｍｉｃｐｅｒｆｏｒｍａｎｃｅｏｆｔｈｅｃｏｏｌｉｎｇｆａｎｉｓａｎａｌｙｚｅｄ Ｔｈｅｒｅｓｕｌｔｓｓｈｏｗｔｈａｔｕｎｄｅｒｔｈｅｓａｍｅｒｏｔａｔｉｎｇｓｐｅｅｄ，ｗｈｅｎｔｈｅｓｔａｔｉｃｐｒｅｓｓｕｒｅｏｆｔｈｅｃｏｏｌｉｎｇｆａｎｉｓｔｈｅｓａｍｅ，ｗｉｔｈｔｈｅｉｎｃｒｅａｓｅｏｆｔｈｅｎｕｍｂｅｒｏｆｂｌａｄｅｓ，ｔｈｅｆｌｏｗｒａｔｅｉｓｌａｒｇｅｒ Ｉｎｔｈｅａｓｐｅｃｔｏｆｓｔａｔｉｃｐｒｅｓｓｕｒｅｅｆｆｉｃｉｅｎｃｙｏｆｃｏｏｌｉｎｇｆａｎ，ｗｈ
ｅｎｔｈｅｓｔａｔｉｃｐｒｅｓｓｕｒｅｏｆｆａｎｉｓａｂｏｕｔ１７０ ２００Ｐａ，ｔｈｅｓｔａｔｉｃｐｒｅｓｓｕｒｅｅｆｆｉｃｉｅｎｃｙｏｆ９ ｂｌａｄｅｆａｎｉｓｔｈｅｈｉｇｈｅｓｔ Ｉｎｏｔｈｅｒｓｔａｔｉｃｐｒｅｓｓｕｒｅｒｅｇｉｏｎｓ，ｗｈｅｎｔｈｅｎｕｍｂｅｒｏｆｂｌａｄｅｓｉｓ７ａｎｄ８，ｔｈｅｓｔａｔｉｃｐｒｅｓｓｕｒｅｅｆｆｉｃｉｅｎｃｙｏｆｔｈｅｆａｎｉｓｈｉｇｈｅｒｔｈａｎｔｈａｔｏｆｔｈｅ９ ｂｌａｄｅｆａｎ Ｔｈｅｒｅｓｅａｒｃｈｉｎｔｈｉｓｐａｐｅｒｃａｎｐｒｏｖｉｄｅａｂａｓｉｓｆｏｒｔｈｅｏｐｔｉｍｉｚａｔｉｏｎｏｆａｅｒｏｄｙｎａｍｉｃｐｅｒｆｏｒｍａｎｃｅｏｆｃｏｏｌｉｎｇｆａｎｓｆｏｒｎｅｗｅｎｅｒｇｙｖｅｈｉｃｌｅｓ Ｋｅｙｗｏｒｄｓ：ｃｏｏｌｉｎｇｆａｎｓ；ａｅｒｏｄｙｎａｍｉｃｐｅｒｆｏｒｍａｎｃｅ；ｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎ；ｎｕｍｂｅｒｏｆｂｌａｄｅｓ  近些年新能源汽车在中国发展迅速，新能源汽车的电子冷却风扇是整车热管理重要组成部分，电子冷却风扇的设计要满足电驱系统、电池系统与空调系统的冷却需求；同时，电子冷却风扇也会对新
能源汽车的ＮＶＨ性能影响很大［１］
．因此，设计出
冷却性能好与低噪音的电子冷却风扇是至关重要的．ＣＦＤ仿真分析技术的出现可以缩短产品的开发周期，同时降低开发成本，更可以从机理上研究
车辆与动力技术２０２０年
冷却风扇的流动细节，目前已经广泛应用到冷却风扇的开发中．当前对冷却风扇的研究主要集中在轮
毂比［２］、叶片个数［３］、叶顶间隙［４］、叶片安装角［５］与叶片形状［６］等方面对冷却风扇性能的影响．
本文使用ＣＡＴＩＡ三维绘图软件创建了风扇的模型，用Ｈｙｐｅｒｍｅｓｈ几何处理，最后使用ＳＴＡＲＣＣＭ＋计算软件求解，介绍了冷却风扇的气动性能数值计算方法、求解方法以及模型计算参数的选取．计算并分析了冷却风扇的静压与流量的关系，数值结果与试验值进行对比，证明了计算模型和计算方法的准确性；同时，对风扇内部压力场、速度场进行分析，所得结果为进一步研究新能源汽车冷却风扇提供了一定的理论基础与参考意义．
１　几何模型
冷却风扇的主要结构参数如下，叶片数：７；风扇半径：２
３０ｍｍ；轮毂比：０ ４；叶片均匀等距分布在轮毂上．由于冷却风扇的原始几何模型较复杂，在网格划分的过程中，将对风扇流场影响小的区域（圆角和孔）用Ｈｙｐｅｒｍｅｓｈ软件简化前处理，风扇计算域网格如图１所示，对叶片附近网格加密处理．按照冷却风扇的试验条件，将整个计算域划分为入
口区、过渡区、旋转区和出口区，进出口区计算域均为半径是２倍风扇半径的圆柱．一般数值计算要求满足入口区长度大于１０倍风扇半径，这里入口区长度３０
００ｍｍ；出口区长度大于１２倍风扇半径，出口长度４０００ｍｍ
．
图１　计算域网格整体划分
２　控制方程与边界条件
２ １　控制方程
不可压缩质量守恒方程为
Ｕｉ／ ｘｉ＝０
，（１）
动量守恒为
ρ Ｕｉ ｔ＋ρＵｊ Ｕｉ
ｘｊ＝－ Ｐ ｘｉ＋μ  ｘｊ  ｘｊ
Ｕｉ＋  ｘｊμτ
Ｕｉ ｘｊ＋ Ｕｊ
ｘ()[]
ｉ
，（２）式中：Ｐ为流体静压；Ｕｉ，Ｕｊ为平均速度分量；ｘｉ，ｘｊ为坐标分量；μ为动力学黏性系数；μτ为湍流黏性系数．
ｋ－ε湍流模型为
［７］
ρ Ｋ ｔ＋ρＵｊ
Ｕｉ ｘｊ＝－  ｘｊμ＋μτσ()ｋ
Ｋ ｘ[]
ｊ＋μτ
Ｕｊ ｘｊ Ｕｊ ｘｉ＋ Ｕｉ
ｘ()
ｊ
－ρε，（３）ρ ε ｔ＋ρＵｉ
ε ｘｉ＝  ｘｊμ＋μτσ()ｓ
ε ｘ[]
ｊ＋ｃ１εｋμτ Ｕｊ ｘｊ Ｕｊ ｘｉ＋ Ｕｉ
ｘ()
ｊ
－ｃ２
ρε２
Ｋ，（４）式中：Ｋ为湍动能；ε为湍流耗散率；στ，σｓ为湍动能和湍流耗散率的普朗特常数；ｃ１，ｃ２为常数．２ ２　边界条件
文中采用在汽车领域应用广泛的ＳＴＡＲＣＣＭ＋商业软件，研究的风扇流场属于低马赫数流动，介质空气属性近似为不可压缩介质，密度为
１ ２０４ｋｇ／ｍ３
，动力粘度为１ ８５５Ｅ ５Ｐａ·ｓ，文中
采用的旋转流动中应用广泛的Ｒｅａｌｉｚａｂｌｅｋ－ε湍流模型，壁面处理为Ｔｗｏ ｌａｙｅｒＡｌｌｙ＋ＷａｌｌＴｒｅａｔｍｅｎｔ，选用多重参考系模型可将风扇旋转的瞬态问题用稳态方法数值求解，压力、动量、湍流耗散率都采用二阶迎风格式．
试验过程中将冷却风扇安装在小型风洞出口处，风洞内在距离冷却风扇入口某一位置处，限定不同静压值，输入１３Ｖ电压，冷却风扇旋转，进而得到不同静压条件下冷却风扇的风量转速、电流、轴功率和效率．为与实验结果相对比，文中入口边界条件为质量流量入口，出口边界条件为压力出口，相对大气压力的静压为０，风扇表面、轮毂
表面、护风罩表面为壁面边界条件．
３　计算结果及分析
３ １　网格无关性验证
计算区域的网格数量对数值求解的数值精度与模拟结果影响很大，数值求解时一般在关键流动区域进行网格细化，在对流动影响不大的区域适当调整网格大小，采用合理的网格参数控制策略既能提
·
２３·
第２期刘继伟等：基于ＣＦＤ的新能源汽车冷却风扇气动性能仿真分析
高计算精度又能节约时间成本．文中计算模型在扇叶周围划分边界层网格，在旋转区域、进出入口区域过度区均采用ｐ
ｏｌｙ多面体网格．文中选取相同流量１ ２７９６ｍ３
／ｓ条件下进行网
格无关性分析，表１为选取的５种不同网格数量条件进行计算，得到冷却风扇的静压值与试验结果进行对比．从表１中可以看出网格数量达到１６０万左右，风扇静压几乎没有变化，并与试验结果很接近，最大误差为０ ５％左右．为保证计算资源和计算时间的限制，后续分析也在此数量网格条件下进行计算．
表１　网格无关性验证
总网格数转速／（ｒ·ｍｉｎ－１）流量／（ｍ３·ｓ－１
）
广州汽车展
静压／Ｐａ５３３５３２２０６７１ ２７９６９５ ４８３３８２６２０６７１ ２７９６９８ ０１６７５４９１２０６７１ ２７９６１００ ３２６３６６６５２０６７１ ２７９６１００ ２６５１２９５０２０６７１ ２７９６１００ １试验值
２０６７
１ ２７９０
９９ ８
３ ２　计算结果验证
图２为试验结果与数值计算结果的静压与冷却风扇流量的关系对比，图３为试验结果与数值计算结果的静压与静压效率的关系对比，可以发现仿真结果与试验测试结果趋势大体相同，静压试验值和仿真结果最大误差为３ ６６％，最小误差为０ ０３％，静压效率最大误差为２
１％，最小误差为０，说明文中采用的网格精度与计算模型可以较精确的仿真冷却风扇的流动状态，可以为后面研究冷却风扇的气动性能研究提供理论支持
．
图２　
风扇静压与流量的关系
图３　风扇静压与静压效率的关系
３ ３　内部流动特性分析
研究冷却风扇工作机理对冷却风扇的开发有着重要作用，分析静压云图可以了解冷却风扇的受力情况，也可以发现流经风扇气体的流动状态．图４和图５为转速为２０６７ｒ／ｍｉｎ，入口流量为
２２２６ｍ３
／ｈ，冷却风扇压力面和吸力面上的静压分
布，可以发现风扇压力面存在很大的正向压力区，且压力呈阶梯状分布，在扇叶的前端压力变化最
大，此区域是冷却风扇做功的主要区域［８］；风扇
吸力面存在很大的负压区域，在叶片的前缘与外圈交接处正、负压相差较大，出现回流区域，此处会造成功率损失，风扇气动噪声也主要在此区域发
生［９］．冷却风扇吸力面压力最高达１００Ｐａ，压力
面正压最大的位置与吸力面负压最小在扇叶上的位置相同，分别位于扇叶的两面．冷却风扇的轮毂部分为封闭实体，无气流通过，压力分布也比较均匀接近为零
．
图４　冷却风扇压力面静压分布图
·
３３·
车辆与动力技术２０２０
年
图５　冷却风扇吸力面静压分布图
如图６和图７所示，轮毂附近与靠近轮毂的气体速度很小，气体流动速度最大位置出现在风扇旋转方向叶片前缘处．同时，在扇叶相互作用的影响下，压力面后缘具有比压力面前缘更大的速度
．
图６　
冷却风扇压力面速度分布图
图７　冷却风扇吸力面速度分布图
３ ４　叶片个数对风扇气动性能的影响
研究在相同转速和相同结构情况下，扇叶个数与风扇性能的关系．图８为不同叶片个数条件下风扇流量与静压的关系．图９为不同叶片个数条件下风扇静压与效率的关系．在同一静压条件下，叶片个数越多，冷却风扇的流量越大，整体情况也符合叶片个数越多，流量越大的事实．９叶片风扇静压在１７０ ２００Ｐａ左右时，风扇静压效率最高，但是在冷却风扇工作的整个区间静压效率并不是随着叶片个数增加而变高．从图９中可以看出，在冷却风扇工作的其他静压区间，７、８叶片时高效率表现则
优于９叶片风扇
．
图８　
不同叶片个数条件下风扇静压与流量的关系
图９　不同叶片个数条件下风扇静压与效率的关系
由图１０可知，气体流动到冷却风扇轮毂处，空气不能通过，会出现低速流动区；气体均匀从入口处流经冷却风扇扇叶，气体流过冷却风扇后速度明显增加，由于风扇旋转的离心力的影响，气体流动方向会变化，变成螺旋流动；随着流动距离的增
·
４３·
第２期刘继伟等：基于ＣＦＤ的新能源汽车冷却风扇气动性能仿真分析
加，气体流动速度逐渐降低
．
图１０　流场内速度流线图
４　结　语
文中介绍了使用ＳＴＡＲＣＣＭ＋软件，采用多重参考系模型求解某新能源汽车冷却风扇的工作特性，介绍了仿真计算的方法与相关参数的选取，通过与试验结果的比较，证明了本文所用方法数值模拟冷却风扇流场特性的正确性；得到了冷却风扇流动区域的压力场与速度场等内部流场特征，分析了冷却风扇的流场特征；揭示了风扇的叶片个数对风扇流量的影响，并分析了风扇叶数对风扇效率的影响，本文的研究可为优化新能源汽车冷却风扇的性能提供依据．
参考文献：
［１］　周　杨 纯电动汽车冷却系统数值分析［Ｄ］．武
汉：华中科技大学，２
０１５．［２］　王振宁，王　红 基于计算流体力学汽车冷却风扇
优化设计［Ｊ］．机械设计与制造，２０１６（１０）：１８２ １８６．
［３］　赵要珍 轿车发动机冷却风扇的ＣＦＤ分析与低噪声
优化设计［
Ｄ］．长春：吉林大学，２００７．［４］　赵立杰，王新玲，署恒涛，等 轴流式发动机冷却
风扇气动性能与气动噪声数值分析［Ｊ］．沈阳航空航天大学学报，２０１７，３４（０１）：５０ ５６．
［５］　王银姣，卢剑伟，江　斌，等 利用ＣＦＤ技术研究
叶片斜度对贯流风机性能的影响［Ｊ］．合肥工业大学学报：自然科学版，２０１２，３５（７）：８８２ ８８７．
［６］　李盛福，王欣欣 汽车冷却风扇叶片参数优化设计
分析［Ｊ］．机械设计与制造，２０１９（０７）：４８ ５２．
［７］　谭礼斌，袁越锦，黄　灿，等 旋风分离器流场分
析与结构优化的数值模拟［Ｊ］．陕西科技大学学报，２
０１８，３６（０３）：１５２ １５９．［８］　庞　磊，李孝宽，李　嵩，等 对旋轴流通风机气动
性能的数值预估［Ｊ］．风机技术，２００８（０４）：２０ ２２．
［９］　ＮＡＳＨＩＭＯＴＯＡ，ＦＵＪＩＳＡＷＡＮ，ＡＫＵＴＯＴ，ｅｔａｌ
Ｍｅａｓｕｒｅｍｅｎｔｓｏｆａｅｒｏｄｙｎａｍｉｃｎｏｉｓｅａｎｄｗａｋｅｆｌｏｗｆｉｅｌｄｉｎａｃｏｏｌｉｎｇｆａｎｗｉｔｈｗｉｎｇｌｅｔｓ［Ｊ］．ＪｏｕｒｎａｌｏｆＶｉｓｕａｌｉｚａｔｉｏｎ，２００４，７（１）：檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲檲８５ ９２．
（上接第３０页）
４　结　论
１）针对某插电式混合动力汽车设计了一套整车电机冷却热管理系统，来保证动力系统、电池系统、空调系统等在各模式／工况下的安全可靠运行．
２）基于三维ＣＦＤ仿真分析与一维系统仿真分析相结合的方法，计算了电机冷却系统在纯电动模式，低速蠕行工况和６０ｋｍ／ｈ爬坡工况下系统的温度和流量，评估系统设计可行．
３）考虑到热管理系统的能耗，对电机冷却系统中电子水泵及其控制策略进行优化．计算了优化前后，春秋季、夏季两种环境下，城市循环工况（
３０次）电机冷却系统电子水泵总能耗．计算结果显示，优化后，春秋季、夏季环境，城市循环工况下，电子水泵能耗分别降低５４％和８５％，能耗降
低明显．
参考文献：
［１］　国家质量监督检验检疫度总局和国家标准化管理委
员会．ＧＢ １９５７８—２００４乘用车燃料消耗量限值［
Ｓ］．北京：中国标准出版社，２００４．［２］　国家质量监督检验检疫总局和国家环境保护部．
ＧＢ１８３５２ ６ ２０１６，轻型汽车污染物排放限值及测量方法（中国第六阶段）［Ｓ］．北京：中国标准出版社，２０１６．
［３］　周能辉，赵春明，辛明华，等．插电式混合动力轿
车整车控制策略的研究［Ｊ］．汽车工程学报，２０１３，３５（２）：９９ １０４．
［４］　卢　山，卢桂萍，李　馨，等．基于Ｖ模式开发插
电式混合动力汽车整车热管理控制策略研究［Ｊ］．计算机测量与测试，２０１８，２６（４）：８８ ９１．
［５］　李　峰．插电式混合动力汽车热管理系统开发及其
控制算法研究［Ｄ］．长春：吉林大学，２０１６．
·
５３·