汽车发动机舱流热耦合问题是进行发动机舱热管理优化的关键问题。随着汽车工业技术的不断发展,发动机的功率不断提高,散热量大幅增加
[1],同时,舱内各部件的布置更加紧凑,留给冷却系统的空间越来越小,使得汽车发动机舱的热管理需求越来越迫切
[2-4]。汽车发动机舱内传热与流动息息相关,二者互相影响,研究汽车发动机舱内流动与传热的关联性可以指导优化发动机舱内的热设计,提高汽车的安全性、经济性、环保性等性能。
目前,研究者们对汽车发动机舱内热管理问题开展了大量研究,主要研究手段是数值模拟,通过模拟发动机舱内的流动与传热过程,并以此提出改进机舱散热性能的方法,以改变空气对流形式,将发动机舱内的热量尽快转移到车身外部,主要改进方法包括:①调整舱内结构布置和工作参数,如进气格栅
[5-10]、散热器
[11-14]、风扇
[5,7, 15,16]、排气管
[5]等;②增加局部结构部件以改变流动状态,如增加导流板
[17-22]、导流罩
[18]、密封挡板
[17,20]等。相对而言,实验研究较少,主要研究手段为风洞实验
[8,23-28]与道路实验
[18,24,29],当前大多数实验研究仅关注发动机舱的宏观参数,如进气量、冷却效率、散热器散热率等。相比较而言,数值模拟方法虽然可提供高分辨率流场,但计算精度依赖湍流模型的准确性,需要实验数据来进行标定。实验研究虽能复现真实工况,但由于发动机舱内结构布局复杂,流入发动机舱内的气流具有很强的三维特性,并伴随流动分离、回流、涡旋等特征
[18],使得实验测量十分具有挑战性。当前缺乏发动机舱内详细、精准的流场及温度场实验数据,导致汽车发动机舱内流热耦合机理不清晰,相应地,数值模拟研究也缺乏用于模型标定实验数据。
针对上述问题,本文使用自动PIV技术与热电偶对汽车发动机舱等效模型的流场与温度场进行测量,研究了流动与传热间的关联机理,并提出了流动控制方法。
1 实验方法
1.1 汽车发动机舱等效模型
本文关注发动机舱宏观上的流热耦合机制,综合考虑实验测量的可行性以及宏观流场的真实性,选用某型燃油汽车发动机舱等效模型作为研究对象,其结构示意图如
图1所示,等效模型由发动机舱与排气管舱两部分组成,去除了真实汽车发动机舱内特征尺寸较小的管道、线缆等部件,保留了简化的进气格栅、冷却模块、发动机、空气滤清器、变速箱、排气管部件,这些大尺寸部件及流动路径对宏观上的流动与传热起主导作用。虽然小尺寸部件会影响细观上的流动与传热,但对发动机舱宏观上的流场与温度场分布影响较小。本文的实验结果表明(详见第2.1节),发动机迎风高温区位置与文献[
18]中真实发动机舱的典型热害区域基本一致,说明等效模型具有合理性。发动机舱冷却模块包括冷凝器、散热器、PTC加热板、风扇,PTC加热板装有自动控温系统,控温稳定性在± 3 ℃范围内,通过加热PTC板至恒定温度,可模拟真实情况下散热器、冷凝器所处的高温状态。
舱内气流由风扇独立驱动产生,气流从进气格栅进入,从排气管舱的出口流出。风扇转速Ω=1 830 r/min,转动方向为x轴正向,叶片半径R=0.2 m,设置PTC加热板温度T0=80 ℃,模拟汽车怠速情况下冷却模块所产生的热气流在舱内的流动情况,实验室温Ta 位于14~19 ℃范围。
本文关注发动机部件迎风区温度场与流场分布特征,该区域直接承受风扇所射出的热气流,所形成的温度场分布特征与发动机外壳的散热效率密切相关,具有重要的研究价值。下文相关物理量的上标“*”表示以风扇叶尖线速度、风扇叶片半径R为基本量作无量纲,温度变量以加热板温度T0及室温Ta为参考量作无量纲,旨在消除风扇转速、散热器温度等输入条件对结果的影响,便于将实验结果推广至不同输入条件的发动机舱模型。
1.2 温度场测量方法
使用T型热电偶感温线,测量发动机舱模型内部的时均温度。根据发动机舱模型的内部结构,制作不同形状的热电偶耙件,以完成流向截面温度场测量,热电偶耙件使用紫铜管作为外部框架,以直径0.38 mm高强度尼龙线编制热电偶悬挂节点,每个节点的横向及纵向间距为70 mm,在耙件节点上布置热电偶阵列。
本文使用的两种热电偶耙件形状如
图2所示,对
图3所示发动机迎风区的4个流向截面(
x/L = 0.15, 0.19, 0.24, 0.31)的时均温度场进行测量。其中
x/L = 0.15、0.19、0.24截面使用
图2(a)所示热电偶耙测量,
x/L = 0.31截面使用
图2(b)所示热电偶测量。每个热电偶的采样频率均设置为2 Hz,采样时间150 s。本文对同一截面的温度场进行了3次重复性实验,得到时均温度的测量不确定度约为1.2%(95%置信区间)。
1.3 流场测量方法
使用2D2C TR-PIV测速系统对汽车发动机舱内流场进行测量,系统包括高速相机(分辨率1 280 × 800 pixels,配备180 mm f/3.5微距镜头)、Nd-YAG激光器(波长532 nm,最大能量40 mJ,频率0.1~10 kHz)、BNC 575同步控制器。选用烟饼燃烧后产生的烟雾颗粒作为示踪粒子。
由于发动机舱内流场具有很强的三维特性,采用研发的自动PIV技术
[30, 31]对发动机舱内的流场进行测量,如
图4与
图5所示,自动PIV测试系统由PIV系统、精密滑台、控制软件所组成,PIV相机与激光器分别搭载在两根性能参数一致、平行安装的滑台上,测试过程中通过控制滑台使PIV相机与片光源同步移动,以改变PIV的测量位置,移动过程中保证PIV相机和激光片光每次运动相同的距离,使PIV相机的焦平面始终与激光面重合,因此只需要在实验前对相机进行一次标定,即可在控制软件控制下,快速完成一系列平行截面流场的扫描测量。使用自动PIV技术分别从两个相互垂直的方向进行扫描测量后,可获得时均速度的3个分量,进而重构测量空间内的三维时均速度场。据估计,使用自动PIV技术相较于传统的PIV技术,测试效率提高了约50倍,且大大降低人为调整测量位置可能引入的未知不确定性。
本文关注发动机迎风区的两处高温区(分别称为区域1、区域2,详见2.1节)内的流场情况,使用自动PIV技术分别对区域1、区域2内流场进行扫描测量(见
图4、
图5),并重构三维时均流场。其中,区域1测量区域为:0.09 <
x/L < 0.39、0.48 <
y/B < 0.93、0.08 <
z/W < 0.33,区域2测量区域为:0.09 <
x/L < 0.28,0.1 <
y/B < 0.37,0.33 <
z/W < 0.6。
由于发动机舱模型尺寸较大,为保证PIV拍摄质量,将每个测试截面分为若干小区域分别进行拍摄,在数据处理后进行拼接获得最终的统计平均结果,PIV采样频率设置为300 Hz,每个测量小区域采样图像对数为3 000对。本文对同一截面的流场进行了3次重复性实验,得到流场的时均合速度的测量不确定度约为2.3%(95%置信区间)。
2 结果与讨论
2.1 温度场分布特征
图6(彩图见电子版,下同)给出了发动机迎风区不同流向截面的平均温度场分布,其中
为无量纲平均温度,
为时均温度,图中黑色虚线圆框表示风扇叶尖所在圆周。可以发现,冷却模块所产生的热量经风扇排出后,逐渐向发动机迎风区的两个区域聚集,形成高温区,分别位于:①发动机上方与侧壁交界处,如
图6(b)~(d)红框所示,定义为区域1,其范围为0.19<
x/L<0.31、0.45<
y/B<0.9、
z/W<0.35;②发动机下方与变速箱交界处,如
图6(b)(c)蓝框所示,定义为区域2,其范围为0.19<
x/L<0.24、0.15<
y/B<0.25、0.3<
z/W<0.45。
2.2 流场分布特征
为弄清发动机迎风区域高温区产生的原因,利用自动PIV技术分别对区域1、区域2的流场进行测量。
2.2.1 区域1
区域1二维切面的时均涡量及速度矢量分布如
图7所示,图中箭头方向代表时均速度矢量方向,箭头的长度代表时均速度的大小,其中
为展向(
z向)无量纲平均涡量,
为法向(
y向)无量纲平均涡量。由
图7 (b)
y/
B=0.48, 0.57截面的流场结果可以看出,由风扇射出的高速热气流直接冲击至区域1发动机的迎风面,并与周围气流发生强剪切作用,由于机舱侧壁的阻挡,热气流在向
z轴负向的运动过程中呈现回旋形运动(
图7(b)黑框所示区域),结合
图7(a)气流的流动方向可知,热气流在区域1内呈现螺旋形运动方式,沿发动机迎风面向上方运动。根据Antonia与Zhou等
[32, 33]对受热圆柱绕流的实验研究,圆柱尾流区的温度分布与尾流结构密切相关,大尺度涡旋结构能够裹挟加热圆柱周围的热量向下游运动,并且这些涡结构输运热量的能力要比动量交换更强,Seuntiëns等
[34]的LIF实验结果显示,受热圆柱尾流大尺度涡结构所在区域的温度明显高于其他区域,这些结果说明涡旋结构在一定程度上具有储存热量的能力。因此,区域1螺旋形气流的存在将使得热量在区域1内长时间滞留,导致温度升高。此外,
图7(a)左上方的气流的运动速度极小,即当热气流运动至机舱上方位置时会承受逆压梯度的作用,热气流无法从发动机与机舱上壁板的间隙流向机舱后方,而是在靠近风扇上方的区域形成回旋气流(
图7(a)蓝框所示),沿展向方向(
z轴正向)流出区域1(
图7(b)
y/B=0.93截面),即热气流流出区域1时亦呈现螺旋形的运动方式,加剧了热量在区域1内的滞留,导致区域1内温度升高。
为更加直观地反映区域1内气流的运动形式,根据自动PIV的测量结果,对区域1的三维时均流场进行重构。
图8给出了区域1内三维时均流线示意图,图中
为无量纲时均合速度,蓝色截面表示PIV拍摄截面。从
图8中可以直观地看出,热气流以螺旋形运动方式流入、流出区域1。上述结果表明,螺旋形运动的热气流是导致区域1温度升高的原因。
2.2.2 区域2
区域2二维切面的时均涡量及速度矢量分布如
图9所示,与区域1结果类似,区域2同样受到风扇射出的高速热气流的冲击(见
图9(a)),由于受到机舱底板的阻碍作用,热气流在沿发动机迎风面向下运动的过程中发生偏转,在
x-y平面内呈回旋形运动(
图9(a)红色方框所示),结合
图9(a)所示
x-z平面内的流动状态可知,区域2内热气流以螺旋形运动方式向
z轴负向方向运动,这同样会使热气流在区域2内发生滞留,导致区域2温度的升高。值得注意的是,在变速箱与发动机之间的交界处(0.19 <
x/L < 0.24、0.15 <
y/B < 0.25、0.45 <
z/W < 0.52),存在低压回流区(
图9 (b)黑框所示),回流区会诱导机舱底部的低温回流上卷,对该区域进行冷却,因此图6(b)(c)黑框所示区域相对于区域2而言温度较低。
同样,根据
图9各截面的流场测量结果,使用三维重构的方法,给出了区域2三维时均流线示意图(见
图10),从图中可以直观地看出流动呈螺旋形运动方式由变速箱一侧向
z轴负向运动。上述结果表明,区域2内的高温亦是由螺旋形热气流的运动所导致。
2.2.3 真实情况下发动机外壳对流换热分析
根据Guo等
[35]提出的强化对流散热场协同理论,对于二维边界层,固壁任意位置的努塞尔数可表示为:
式中:Re为特征雷诺数;Pr为流体的普朗特数;为场协同数;、、分别为无量纲速度、无量纲温度、无量纲壁面法向坐标;U∞、T∞分别为自由来流速度和温度,Tw为壁面温度,δt为温度边界层厚度。公式(1)表明,固壁的对流换热强度不仅仅取决于Re与Pr,还取决于场协同数的大小,即无量纲速度矢量与无量纲温度梯度的点积在热边界层内的积分,该积分的大小取决于速度矢量与温度梯度之间的夹角。真实情况下,汽车发动机为高温部件,周围气流对发动机外壳产生冷却作用,即固壁附近温度梯度指向固壁方向,该情况下,应使气流速度矢量与温度梯度夹角尽量接近0°,以提高发动机外壳的对流换热强度。
对于区域1,观察
图14(a)中
z/W = 0.08截面,可以发现固壁附近的速度矢量方向接近壁面切向,即速度矢量夹角与温度梯度夹角远离0°,这会降低
z/W = 0.08截面附近发动机外壳的对流换热强度,对散热产生不利影响。同理,对于区域2,
图10(b)中
y/B = 0.19 截面处 0.25 <
z/W < 0.4 区域内固壁附近的气流速度矢量的方向接近于固壁切向方向,说明
y/B = 0.19 截面附近的发动机外壳对流换热强度较低。区域1与区域2的局部高温也会降低发动机外壳附近的温度梯度,使对流换热强度进一步下降。因此,在实际工程应用中,应重点关注区域1与区域2位置处发动机外壳的冷却问题。
2.3 流动控制方法
当前燃油车发动机舱发动机功率大、散热需求高,而发动机舱紧凑化的布局也使得热害问题进一步加剧,使热管理面临挑战。最佳的流动控制策略是避免大幅改变舱体结构,进行非侵入式、精准的局部降温,以适用于现有车型的快速优化。
基于上述考虑,本文基于区域1和区域2局部高温的产生机理,提出流动控制方法(见
图11),进行精准的流动控制,实现局部冷却。以下分别对区域1和区域2的流动控制方法进行阐述。
2.3.1 区域1
针对区域1螺旋形热气流的滞留问题,在发动机部件上方靠近侧壁板位置(0.05<
z/
W<0.26),加装3个沿展向并排放置的小风扇(见
图12),旨在驱动区域1处的热气流向机舱后方区域释放,减弱热气流的滞留。小风扇长宽尺寸为60 mm × 60 mm,额定转速为4 600 r/min。
图13给出了加装小风扇后区域1二维切面的时均涡量及速度矢量分布,比较加装小风扇前后的流场分布(
图7与
图13),加装小风扇后,在小风扇所产生的顺压梯度作用下,更多的热气流直接自发动机上方(
图13(a)左上方)向机舱后方运动,机舱上方回流减弱,并且由回流引起的回旋气流涡量降低(
图13(a)蓝框所示),同时,流动在
x-z平面内的回旋强度(
图13(b)黑框所示)减弱,使得热气流的螺旋形上升形式被抑制,区域1内热量的滞留时间减少。从
图14中的三维流线,也可以直观地看出,加装风扇后区域1内热气流的螺旋形上升形式被抑制。在机舱上方区域,有少数热气流仍从机舱上方沿展向方向流出,但这部分气流的流量相比于原机舱已大幅度降低,更多的气流在小风扇的驱动下直接从发动机上方向机舱后方流出区域1。
图15给出了进行流动控制后
x/L=0.31截面的时均温度场分布,可以发现,加装小风扇后区域1的平均温度有了显著下降,比较该截面0.11<
z/
W<0.27, 0.55<
y/
B<0.91范围内热电偶测点的平均温度,加装风扇后平均温度相比于原机舱下降了4.5%。
上述结果表明,加装小风扇可以改善热气流在区域1内的滞留情况,抑制热气流的螺旋形运动形式,降低区域1内的平均温度。
2.3.2 区域2
由
图6与
图9可知,在变速箱与发动机的交界位置存在由流动分离引起的低压区,在低压区的引导下,发动机与机舱下壁间的低温回流会上卷,降低该区域的温度。受到这一启发,针对区域2的局部高温问题,在变速箱侧壁加装导流板(见
图16),旨在延缓气流在变速箱表面分离后向发动机表面的再附过程,产生更大范围的低压区,引导发动机与机舱下壁间更多的低温回流上卷,对区域2进行冷却。
图17给出了加装导流板后区域2二维切面的时均涡量及速度矢量分布。与无导流板情况相比,加装导流板后,展向流动的流动分离被延后至导流板尖端处(见
图17(b)),使得导流板与发动机之间产生更大范围的低压区(
图17黑框所示),在低压区的诱导下,发动机下方的低温回流发生上卷(
图17(a)红框所示),相比于无导流板情况(见
图9),上卷气流的展向覆盖范围更大、上卷速度更高。从
图17(a)以及
图18的三维时均流线中也可以看出,呈螺旋形运动的热气流也被导流板阻隔在机舱上游一侧,减少了热气流与发动机表面的直接接触。
图19给出了加装导流板后流向截面
x/L=0.24的平均温度场分布情况,可以发现区域2的平均温度有了明显下降,在该截面0.27<
z/
W<0.43, 0.11<
y/
B<0.28范围内,加装导流板后平均温度相比于原机舱下降了8.3%。
上述结果表明,加装导流板可以实现局部冷却,导流板能够引导低温气流上卷,阻隔螺旋形热气流,有效地降低区域2内的平均温度。
3 结 论
(1)在汽车发动机舱发动机迎风区域存在两处由螺旋形气流聚集引起的局部高温区,分别位于发动机上方与侧壁交界处、发动机下方与变速箱交界处。
(2)通过加装小风扇、导流板等流动控制措施,削弱了螺旋形气流运动,引导低温流体进入高温区,实现局部冷却,两处高温区平均温度分别下降了4.5%和8.3%。
国家重点研发计划项目(2022YFE0208000)
广汽研究院横向项目(GACRND-2022-JS0328(A))