【空调系统】新能源车用 CO2 空调系统泄漏特性仿真研究

导读

CO2作为一种安全､高效且环境友好的自然工质,是未来汽车空调制冷剂替代方案的重要选择｡考虑到CO2系统制冷剂泄漏导致乘员舱气体体积分数超标的舒适及安全问题,本文使用三维仿真软件STAR-CCM+搭建乘员舱CO2制冷剂泄漏仿真模型,仿真分析两种不同泄漏速率和泄漏圆孔孔径下乘员舱乘客面部测点CO2体积分数的动态变化情况｡在50g/s的高CO2泄漏速率下,面部CO2体积分数在泄漏完成时能够达到9%以上,需要切换送风模式提供新风来降低乘员舱内CO2体积分数;在0.1g/s的低CO2泄漏速率下,泄漏完成时面部CO2体积分数不超过3%,处于呼吸安全区范围｡

关键词：CO2;热泵空调;泄漏特性;制冷剂

作者：宗硕 殷翔 黄龙飞 宋昱龙 曹锋

西安交通大学能源与动力工程学院西安

传统的汽车空调系统大部分以R134a[1]作为制冷剂,R134a因其良好的制冷性能而被广泛应用于各个领域,但高达1430的GWP(全球变暖潜值,GlobalWarmingPotential)[2]使其在环保效益方面有所欠缺｡随着各类法律法规的出台,R134a等传统制冷剂将逐步淘汰使用｡新能源汽车热管理系统目前仍广泛使用HFC类强温室效应气体(如R134a和R407C)作为制冷剂,而CO2作为自然工质无毒没有可燃性[3],其ODP(臭氧损耗潜值,ozonedepletionpotential)为0,GWP为1[4],是目前作为替代传统制冷剂的理想选择[5]｡CO2热泵技术在车辆热管理领域发展前景广阔,制热能力明显优于采用R134a､R407C､R1234yf等制冷剂[6]的方案｡对于跨临界CO2汽车热泵空调系统,由于制冷剂CO2可能通过HVAC系统泄漏至乘员舱中,导致乘员舱内CO2气体体积分数过高影响乘员舱舒适性[7],因此研究CO2泄漏后乘员舱驾驶员和乘客呼吸处CO2体积分数是有必要性的｡

若CO2体积分数超过人体可接受的阈值,则会产生威胁到人体健康的后果[8],且乘员舱内人体不断通过呼吸产生补充CO2｡根据美国疾病控制和预防中心[9]的数据,人体在体积分数为5%CO2下暴露30min或在7%~10%CO2体积分数下暴露几分钟后将出现头晕的症状｡EPA(美国环境保护署,Environ-mentalProtectionAgency)将体积分数3%定义为乘客自由空间内平均超过15min的短期呼吸水平,并将4%定义为乘客呼吸区的上限｡因此,在制冷剂泄漏特别是由蒸发器严重泄漏引起的情况下,评估汽车乘员舱CO2体积分数对于制定应对措施以确保乘客的安全至关重要,泄漏的开始可能是由于汽车碰撞事件､HVAC系统的故障或长期使用导致的腐蚀｡G.D.Mathur等[10]通过实车测量乘客呼吸引起的CO2体积分数变化,乘员舱在3名乘客的呼吸作用下平均CO2体积分数在50min内能够达到0.25%,空调在再循环模式下运行｡H.Jung[11]开发了一个数学模型,用于预测空调再循环模式下乘员舱中的CO2体积分数变化,作为乘员舱容积､车身泄漏和乘客数量的函数,通过拟合实验数据确定相关系数｡徐廷喜等[12]提出一种基于支持向量数据描述算法的变频空调系统制冷剂泄漏故障检测和诊断方法｡

S.Yadav等[13]研究了R290和R1234yf两种潜在替代制冷剂在蒸发器､冷凝器中的泄漏特性,并对其由于泄漏引起的爆炸特性进行了分析｡M.E.Ko-ban等[14]采用CFD模型模拟HFO-1234yf在不同空调管路破裂工况下的扩散特性,结果显示,由于空气被夹带到制冷剂射流中,制冷剂的泄漏过程得到快速稀释｡刘全义等[15]基于Fluent软件模拟了水平和垂直方向不同位置抽气口作用对建筑物内泄漏丙烷后的混合气体体积分数场分布的影响规律,结果显示,制冷剂泄漏后,垂直方向抽气口位置越低,对于降低房间内丙烷气体体积分数的效果越好｡杨清泉等[16]采用实验的方法对R32制冷剂在运行过程中发生的泄漏和爆炸的问题进行研究,结果表明,不可燃时间值随着泄漏速度的增加呈先下降后上升的趋势｡

R.Monforte等[17]使用CFD仿真模拟以评估R1234yf制冷剂泄漏至乘员舱内的风险,该CFD模型通过两种实验方法进行了验证,可在简化的乘员舱中进行测量和车辆测试｡汪琳琳等[18]通过数值模拟对R1234yf制冷剂在蒸发器破损泄漏随送风进入乘员舱后的体积分数变化过程和最高体积分数进行了动态监测,研究发现汽车后排较前排更易发生制冷剂聚集｡袁小勇[19]通过实验和数值模拟的方法,分析以R32空调系统发生制冷剂泄漏的安全性,得到泄漏速度､窗开启面积､障碍物与泄漏口的距离等因素对泄漏扩散规律､体积分数分布及房间内安全性的影响｡

本文通过CFD仿真软件STAR-CCM+建模,模拟制冷剂泄漏时乘员舱内乘客面部测点的体积分数动态变化,泄漏的场景参考SAEJ2772泄漏标准的测试场景[20],通过该模拟对两种CO2制冷剂泄漏速率的场景做相应的预估和分析,进而提出相应的应对方法,CFD模型能够准确地仿真乘员舱内气体CO2体积分数随时间的变化趋势｡

01  乘员舱CO2气体泄漏仿真模型建立

几何模型和计算模型

图1所示为乘员舱CO2气体泄漏仿真的几何模型,该几何模型包含HVAC风道､乘客､座椅以及车门缝｡HVAC风道模型有一个总入口､回风口和连接到乘员舱内部的出风口,本节的仿真模型只开前排4个吹面的出风口,分别位于中间位置的两个出风口和两侧的两个出风口｡汽车内部容积为2.8m3,乘员舱内前后排座椅和5名乘客,在每个乘客的面部上设立一个出口用于考虑人体呼吸作用产生的CO2｡车门缝参照实际汽车的车门缝隙尺寸,用于考虑外部环境通过车门缝隙泄漏进入到乘员舱内部的新风｡

美国汽车工程学会标准SAEJ2772中选取了腐蚀泄漏和碰撞泄漏两种工况作为R744制冷剂泄漏典型工况[20]｡在腐蚀泄漏工况下,泄漏孔直径为0.1mm,对应R744泄漏速率为0.1g/s;在碰撞泄漏工况下,泄漏孔直径为6.35mm,对应R744泄漏速率为50g/s｡因此,本文仿真研究中泄漏孔的直径设置为6.35mm和0.1mm,分别对应高流量(50g/s)和低流量(0.1g/s)两种泄漏速率｡CO2泄漏孔的位置位于蒸发器中央位置的上1/3高度处,如图2所示｡考虑到实际高速泄漏工况存在节流和蒸发等复杂的物性变化过程,较难以三维仿真形式实现｡同时HVAC风道的存在有利于均衡客舱内部出风口的CO2体积分数｡因此本文对于CFD模型中的泄漏孔进行合理简化,仅考虑泄漏孔径尺寸和泄漏质量流量边界条件｡

仿真求解中将乘员舱内部的空气当作不可压缩流体,湍流模型选择realizable模型｡CO2制冷剂泄漏过程中不与乘员舱内的空气发生化学反应,选择非反应模型｡气体模型选用理想气体模型,满足理想气体状态方程｡CO2生热量交换,制冷剂泄漏过程中温度恒定,不产生热量交换,忽略重力的影响｡

边界条件

CFD模型边界条件设置为:

1)泄漏孔设置为质量流量进口,流体为CO2,泄漏速率为50g/s和0.1g/s,其中50g/s泄漏持续时间为8s,0.1g/s泄漏持续时间为4000s,总泄漏量均为400g｡

2)乘员舱内有5名乘客,每名乘客面部设置CO2质量流量进口,速率为20L/h｡

3)空调循环模式为纯回风模式,总风量设置为480m3/h,该模式下CO2泄漏所能达到的体积分数为所有模式下最差的情况,回风口设置为出口｡

4)车身的车门缝与周围环境存在空气交换,空气交换率设为2,其含义为单位时间内进入乘员舱空气量与乘员舱净车内容积的比率｡

CFD模型初始条件设置为:

1)乘员舱内外压力为环境压力101325Pa,乘员舱内外温度为298K｡

2)乘员舱内气体组分为100%空气｡

02  仿真结果与分析

泄漏过程CO2体积分数分布

高流量泄漏为CO2在泄漏孔处泄漏速率为50g/s,泄漏圆孔直径为6.35mm,泄漏时长为8s｡图3所示为高流量泄漏过程中面部测点CO2体积分数变化｡

泄漏过程中乘员舱内CO2体积分数迅速上升,面部测点体积分数在8s内达到8%以上,其中驾驶员测点体积分数达到9%以上,大大超过了乘客呼吸安全区的上限值｡因此在泄漏过程停止后,需要采取相应的措施降低乘员舱内CO2体积分数｡

图4所示为高流量泄漏过程中4个时间点中间特征截面面部测点CO2体积分数云图,泄漏过程中制冷剂CO2从前排吹面出风口吹入乘员舱,使得乘员舱内CO2体积分数迅速上升｡泄漏的CO2气体从乘员舱上方空间流动,先到达乘客的头部再到达面部,因此头部和面部的CO2体积分数提升较为迅速,在极短时间内达到5%以上,泄漏的CO2气体循环后流向回风口的方向｡

泄漏过程动态特性

高流量泄漏过程中,乘员舱内CO2体积分数超过了安全值,应对措施最直接的方法是改变汽车空调的送风模式,改为全新风快速降低乘员舱内CO2体积分数,图5所示为全回风和全新风模式下CO2面部测点体积分数变化｡

如图5(a)所示,在8s的泄漏时间结束后,如果不改变送风模式,在门缝泄漏进来的新风和人呼吸产生的CO2综合影响下,随着时间的推进CO2会逐渐均匀散布在乘员舱内,各面部测点的CO2体积分数逐渐趋于稳定,在100s内面部测点的体积分数仍在

4%以上,超过了呼吸安全区的上限值｡因此,不改变送风方式仅依靠门缝新风和CO2均匀化来降低呼吸体积分数是不合理的｡如图5(b)所示,在8s的泄漏时间结束后改变送风模式变成全新风,在新风的作用下乘员舱内CO2体积分数快速降低,在20s内面部测点体积分数已降至安全区范围内,在50s时已降至1%以下｡

图6所示为低流量泄漏全过程中面部测点CO2体积分数变化,在整个低流量泄漏过程中由于制冷剂CO2泄漏的速率较低,乘员舱内各面部测点CO2体积分数上升的速度较为缓慢,使得汽车乘员舱内CO2气体体积分数较均匀｡由于低流量泄漏下乘员舱内CO2气体体积分数较均匀,各面部测点CO2体积分数的差值较小｡在第4000s时泄漏停止,此时乘员舱内各面部测点CO2体积分数均在3%以下,没有达到CO2气体体积分数舒适性安全区的上限值｡泄漏停止后在门缝漏进乘员舱的新风和人呼吸产生的CO2的综合影响下,各面部测点CO2体积分数逐渐下降趋于平稳｡不同于高流量泄漏情况,低流量泄漏在整个过程中面部测点CO2体积分数均在舒适性安全范围之内,不需要空调系统切换送风模式提供额外的新风来降低乘员舱环境内部的CO2体积分数,因此乘员舱在CO2制冷剂低流量泄漏下基本不会对人体呼吸造成安全影响｡

03  结论

基于CO2汽车空调系统制冷剂泄漏存在的乘员舱气体舒适度隐患,本文使用三维仿真软件STAR-CCM+搭建乘员舱制冷剂泄漏仿真模型,仿真分析两种不同泄漏速率和泄漏圆孔孔径下乘员舱乘客面部测点CO2体积分数的动态变化情况,为CO2系统制冷剂泄漏采取相应措施提供参考｡得到如下结论:

1)汽车空调系统在50g/s的高CO2泄漏速率下需要切换送风模式提供新风来使乘员舱内CO2体积分数达到乘客呼吸安全区,在0.01g/s的低CO2泄漏速率下,CO2体积分数不会达到乘客呼吸安全区的上限值｡

2)CO2制冷剂泄漏速率为50g/s的情况下,面部测点体积分数在8s内达到9%以上,在泄漏结束后切换为全新风模式能够在20s内使乘客面部CO2体积分数降至安全范围｡

3)CO2制冷剂泄漏速率为0.1g/s的情况下,在乘客呼吸作用和门缝漏风的共同影响下,泄漏完成时面部CO2体积分数不会超过3%｡

由上述结论可知,在完善极端泄漏工况送风模式切换功能的前提下,CO2汽车空调系统可以较好的规避乘员舱舒适度安全隐患｡本文研究和分析主要基于仿真方法,限于实验条件暂无相应整车实验验证和校核仿真模型的准确性和可靠性,在未来工作中有必要进行对比实验验证｡

【参考文献】

[1]陈科峰,龚海涛.汽车空调HFC-134a替代品研究进展[J].有机氟工业,2008(4):35-37,62.(CHENKefeng,GONGHaitao.ResearchprogressofHFC-134asubstituteforautomobileairconditioning[J].Organo-FluorineIndustry,2008(4):35-37,62.)

[2]COLEMANJW,GARIMELLAS.Two-phaseflowregimesinround,squareandrectangulartubesduringcondensationofrefrigerantR134a[J].InternationalJournalofRefrigeration,2003,26(1):117-128.

[3]王从飞,曹锋,李明佳,等.碳中和背景下新能源汽车热管理系统研究现状及发展趋势[J].科学通报,2021,66(32):4112-4128.(WANGCongfei,CAOFeng,LIMingjia,etal.Researchstatusandfuturedevelopmentofthermalmanagementsystemfornewenergyvehiclesunderthebackgroundofcarbonneutrality[J].ChineseScienceBulletin,2021,66(32):4112-4128.)

[4]陈凯胜.电动汽车二氧化碳热泵空调制热性能及控制策略模拟研究[D].北京:北京理工大学,2018.(CHENKaisheng.Researchonsimulationofcarbondioxideheatpumpairconditionerheatingperformanceandcontrolstrategyforelectricvehicles[D].Beijing:BeijingInstituteofTechnology,2018.)

[5]刘业凤,王雨晴,唐丹萍.新能源汽车CO2热泵空调系统仿真研究[J].农业装备与车辆工程,2021,59(12):45-51.(LIUYefeng,WANGYuqing,TANGDanping.SimulationresearchonCO2heatpumpairconditioningsystemofnewenergyvehicle[J].AgriculturalEquipment&VehicleEngineering,2021,59(12):45-51.)

[6]宋昱龙,王海丹,殷翔,等.跨临界CO2蒸气压缩式制冷与热泵技术综述[J].制冷学报,2021,42(2):1-24.(SONGYulong,WANGHaidan,YINXiang,etal.ReviewoftranscriticalCO2vaporcompressiontechnologyinrefrigerationandheatpump[J].JournalofRefrigeration,2021,42(2):1-24.)

[7]GARCÍA-LÓPEZE,HEARDC.Astudyofthesocialacceptabilityofaproposaltoimprovethethermalcomfortofatraditionaldwelling[J].AppliedThermalEngineering,2015,75:1287-1295.

[8]AMELINEA,DUMESTRE-TOULETV,RAULJS,etal.Determinationofathresholdfatal3-MMCconcentrationinhuman:missionimpossible[J].Psychopharmacology,2019,236(3):865-867.

[9]CHAKRABARTIM,PEREZMONTIELA,CORRILOI,etal.CO2concentrationinthecabinintheeventofaleak:CFDsimulationandtesting[J].SAETechnicalPaper,2017-01-1444.

[10]MATHURGD.Experimentalinvestigationtodetermineinfluenceofbuild-upofcabincarbondioxideconcentrationsforoccupantsfatigue[J].SAETechnicalPaper,2016-01-0254.

[11]JUNGH.ModelingCO2concentrationsinvehiclecabin[J].SAETechnicalPaper,2013-01-1497.

[12]徐廷喜,杜志敏,吴斌,等.基于支持向量数据描述算法的变频空调系统制冷剂泄漏故障诊断研究[J].制冷技术,2019,39(4):25-31.(XUTingxi,DUZhimin,WUBin,etal.Faultdiagnosisresearchonrefrigerantleakageofinverterair-conditioningsystembasedonsupportvectordatadescriptionalgorithm[J].ChineseJournalofRefrigerationTechnology,2019,39(4):25-31.)

[13]YADAVS,LIUJ,KIMSC.Acomprehensivestudyon21st-centuryrefrigerants-R290andR1234yf:areview[J].InternationalJournalofHeatandMassTransfer,2022,182:121947.

[14]KOBANME,HERRMANNDD.DispersionmodelingofleaksoflowglobalwarmingpotentialrefrigerantHFO-1234yfinanautomobilegarage[J].ProcessSafetyProgress,2011,30(1):27-34.

[15]刘全义,张辉,刘奕,等.建筑物空调系统可燃制冷剂泄漏浓度场分布[J].清华大学学报(自然科学版),2013,53(9):1326-1333.(LIUQuanyi,ZHANGHui,LIUYi,etal.Concentrationdistributionofleakingflammablerefrigerantsfrombuildingair-conditioningsystems[J].JournalofTsinghuaUniversity(ScienceandTechnology),2013,53(9):1326-1333.)

[16]杨清泉,袁小勇,金梧凤.R32制冷剂泄漏速度对室内可燃风险的影响研究[J].日用电器,2022(4):30-34,42.(YANGQingquan,YUANXiaoyong,JINWufeng.ExperimentalstudyontheeffectofR32refrigerantleakagerateonindoorconcentrationdistributionandsafety[J].ElectricalAppliances,2022(4):30-34,42.)

[17]MONFORTER,MATTIELLOF,PEROSINOA,etal.ValidationofaCFDmodeltopredictR-1234yfconcentrationsinavehiclecabincompartment[J].SAETechnicalPaper,2016-01-0221.

[18]汪琳琳,杨昭,王丹.R1234yf整车汽车空调泄漏扩散研究[J].制冷学报,2022,43(1):59-67.(WANG