1. 引言
目前,电动汽车已经越来越受群众的喜爱,但其内部乘坐热舒适性也越来为人们所关注。现在汽车为了保证驾驶乘客驾驶安全及审美疲劳,而增大汽车车窗来开阔驾乘人员的视野,但却增大了太阳辐射,影响乘坐人员的热舒适性,甚至影响驾驶安全,所以对乘客舱舒适性研究极为重要 [1]。汽车室内的舒适性主要与空气的温度、速度、质量及太阳辐射有关 [2]。热舒适性在炎热和寒冷的气候表现的最为明显,车内温度过高、过低都会影响人的舒适性 [3]。近年,有关汽车乘客舱及空调的研究比较多,大多在建立三维汽车模型的基础上对送风结构即参数进行研究,如文献 [4] [5] [6] [7] 通过建立模型,研究不同送风条件对车内的热舒适性的影响。文献 [8] 通过CFD软件建立单个通风口的地上、地下两个相同有限空间内的气流运动状态,采用不同网格长度对比分析有限空间内的气流运动,得出外界的风量和风速越大,内部空气循环的越彻底。车内热舒适不能仅考虑送风结构变化对其的影响,还应综合考虑太阳辐射和人体热生理情况,文献 [9] [10] 考虑了人体生理热,引入PMV-PPD等多种热舒适性指标对有限空间的热环境进行评估。文献 [11] [12] [13] 综合考虑了太阳辐射和人体生理热对室内热环境的影响,通过建立仿真模型以及引入评价指标对有限空间内的空气质量和热舒适性进行评估,得出室内环境的舒适性情况。
综上所述,这些研究大多侧重于风道的改造,对于综合考虑太阳辐射和人体自身因素对热环境影响的研究很少,不能精确地反映出热环境的舒适性。因此,本文对SUV 6座汽车内部的热舒适性进行模拟研究,并综合引入各种评价指标对车内的舒适性进行评估,具有一定的实用意义。
2. 模型的建立
2.1. 几何模型及网格划分
本文使用CATIA软件对SUV汽车和假人模型进行简化,简化过程中保留制冷吹面的风道,且风道密封良好,如图1为简化后的SUV汽车模型。然后对整车和假人模型网格划分,STAR-CCM+中,整个计算域内体网格数为850万。图2为假人及汽车网格模型。
2.2. 数学模型
STAR-CCM+中有自带的数值求解器,可以精确模拟车内的热舒适情况。车内空气流速较小,使用k-ε湍流模型,内部流体可看作三维不可压缩流 [14];乘客舱内的空气为辐射透明介质,符合灰体辐射。求解方程如下:
连续性方程:
(1)
Figure 1. Curve: system result of standard experiment
图1. SUV汽车三维几何模型
Figure 2. Curve: system result of standard experiment
图2. 假人及车舱网格划分示意图
动量方程:
(2)
能量方程:
(3)
式中,
是气流密度;
是动力粘度;P为压力,P;T为温度,℃;t为时间,s;
为x、y、z方向的张量;k为热传导;
是比热容。
3. 边界参数及评价指标
3.1. 边界参数
本文将进口设为质量流量入口,出口设为分散流出口,采用瞬时计算3900 s。在STAR-CCM+仿真软件中有自带的参数设置框。
3.2. 评价指标
3.2.1. PMV-PPD
PMV表示预期平均热感觉,其将舒适状态从冷到热划分为7个等级,见表1。由于环境的舒适性因人体的生理情况不同而异。所以不适人员比例PPD被提出来表达人们对热环境的不满意程度,当PMV在−0.5~0.5时,且PPD < 15%时,人体处于舒适状态。PMV和PPD的计算如式(4)和(5)。
(4)
(5)
式中,M——人能量代谢率,W/m2;
W——人做的机械功,W/m2;
——表皮及出汗的散热量,W/m2;
——呼吸潜热损失量,W/m2;
——呼吸显热损失量,W/m2;
H——体表对流和辐射热损失,W/m2。
Table 1. Comfort level of PMV evaluation index
表1. PMV评价指标舒适等级
3.2.2. 空气龄(MAA)
空气龄为新鲜空气在有限空间内停留的时间,可通过示踪气体实验估算其值 [15]。有限空间内某点的平均空气龄计算公式,单位:秒:
(6)
式中,
——有限空间内示踪气体均匀分布时的初始浓度;
——当以某种特定送风方式送入新鲜空气时测得的示踪气体浓度随时间的变化。
4. 模拟结果分析
4.1. 乘客舱空气流速分布情况
图3为模拟过程SUV车内空气流动状态。可以看出,风向从风口出来后直接吹响乘客,离风口较近的部位,气流密度强,对于乘客来说吹风感较猛。图中,乘客从头部到小腿出现涡旋状回流,表现为头暖脚冷,且风口直吹的地方风速较大达到3 m/s,人体表面最大风速达到2.5 m/s。总体而言,乘客头部的风速相对腿脚处较高,但因热胀冷缩原理,车内温差会尽快降到合理范围;图中可以看出人体感到舒适的风速为0.5~2 m/s,前两排乘客体表风速基本能够满足,后排乘客与C柱出风口的距离太近,使得乘客表面及肩部风速稍高。
4.2. 乘客舱成员体表温度分布情况
图4为模拟过程乘客表面的温度分布状态。图中,前排乘客腿部温度较其他乘客高,基本在38℃左右。一方面是该部分新鲜空气流速较低,换热量较少;另一方面是因为前排挡风玻璃较大,腿部吸收的辐射热较多。二排分布乘客体表温度相对均匀。一方面是因为侧窗玻璃的透射率较低,吸收的辐射热较少;另一方面是因为空气流速较均匀,能够迅速降低温度。第三排乘客除后侧头部区域温度较高,其他区域的温度均匀分布。这是因为三排乘客表面风速虽然分布较均匀,但后脑勺部位吸收的辐射热较大。可见乘客体表温度不仅与出风口温度、吹风位置有关,还与乘客舱周围玻璃透射率有关。总体而言,乘客舱处于人体舒适的温度区间内。
Figure 3. Schematic diagram of air flow line in vehicle
图3. 车内空气流线示意图
4.3. 乘客舱平均温度分布
图5为车舱截面温度分布,可以看出乘客舱温度总体处在26℃附近,但车舱轮廓处温度要比车舱内部温度稍高。这是因为边界处直接吸收大量的太阳辐射,且又与外界高温环境进行对流换热,冷空气到达车舱边界时的流速很低,对流换热系数较小,不能及时给边界降温,故使轮廓温度高于内部温度。图中前挡风玻璃与仪表盘的夹角处温度过高,这是由于三维仿真中的网格尺寸大于该处的实际几何尺寸,因而产生计算死角,所得的温度结果失真。
Figure 5. Cross section temperature inside vehicle
图5. 车内横截面温度
4.4. PMV-PPD评价指标分析
图6中,PMV值基本分布在−0.6~0.5之间,乘客处于热舒适状态。因为后两排车窗玻璃的透光率比前排玻璃低,乘客所吸收的太阳辐射热少;且对流换热较好,冷空气均匀分布在体表上,所以后两排乘客基本处于热舒适状态。前排乘客除大腿和腹部位置的PMV值较高,其余处的PMV值处在[−0.5, 0.5]的舒适区间内,这是因为前挡风玻璃的透光率较高,大腿及腹部吸收的辐射热最多,并且该处的冷空气流速较低,换热量较少,热量难被带走。
Figure 6. PMV distribution on passenger surface
图6. 乘客表面PMV分布
Figure 7. PPD distribution on passenger surface
图7. 乘客表面PPD分布
图7,乘客表面PPD分布基本低于15%,可认为乘客的满意程度比较高。前排乘客除大腿及腹部的PPD值高于15%,其余位置处于满意的状态,这是因为乘客的大腿和腹部吸收大量的太阳辐射,且该部位的空气流速低,换热量较少;因后两排乘客的肩颈直对出风口,风速感较强,且从车窗处吸收了部分的辐射热,使人感到极不舒适。
4.5. 空气龄(MAA)
有限空间内,空气龄在40 s内被视为可接受的舒适范围。如图8,乘客表面的MAA分布在25 s内,说明有限空间内的新鲜空气循环的比较顺畅,乘客感觉比较舒服。这里是乘客表面的平均MAA已经达到了舒适的状态,但是车内局部位置MAA较高,如图9,车舱前部的MAA较高,这是因为车舱前部吸收过多的太阳辐射,且风速循环较慢,换热量低,导致温度太高,空气干燥,MAA较高。总的来说,后两排乘客的MAA要比前排乘客的稍低。这是因为后排出乘客接收的新鲜冷空气较多、空气流速较均匀,并且吸收的辐射热较少,因此MAA要比前排的稍低,感觉比较舒适。
Figure 8. RH distribution of passenger surface relative humidity
图8. 乘客表面MAA分布
Figure 9. RH distribution of passenger surface relative humidity
图9. 乘客表面MAA分布
5. 总语
本文通过STAR-CCM+建立了三维SUV 6座车舱模型,引入人体模型和太阳辐射模型模拟车内热环境在高温夏季静置后分别以不同的运行方式启动时车内的热舒适性情况。总得来说,车内及乘客表面温湿度分布较均匀且舒适,稳定时平均温度在26℃附近,但前排乘客的大腿及腹部位置吸收较多的辐射热,使得平均辐射温度值偏高、空气比较干燥,这是因为前挡风玻璃的透射率较大,使PMV-PPD的值超出了舒适的范围。同时,后排出风口的风速偏高,后排乘客的肩颈部吹风感较强,该区域的乘客PPD值偏高。因此,出风口温度过低以及空气流速过快都会影响乘客的乘坐热舒适性,为了满足热舒适性要求,汽车空调系统不仅需要合理控制出风口风速,还需要控制阳光进入乘员舱的辐射量。
基金项目
国家自然科学基金重点项目,51707122;上海理工大学科技发展项目,2018KJFZ179。