电动汽车空调室外换热器空气侧流动换热模拟研究

刘妮,贺剑锋,崔强,等.电动汽车空调室外换热器空气侧流动换热模拟研究[J].低温物理学报,2022,44(02):124-130.

摘要

采用数值模拟方法对室外微通道换热器翅片侧空气流动换热性能进行仿真计算，探讨了在制冷工况下，不同百叶窗结构对微通道换热器空气侧传热及流动特性的影响。结果表明ｊ 因子的模拟结果与实验关联式之间的平均偏差在７．８％以内，ｆ 因子的平均误差在７．３５％以内，符合工程应用要求．雷诺数较低时，传热因子ｊ 和阻力因子ｆ 都随Ｆｐ的增大而减小，雷诺数较高时，Ｆｐ对两者的影响不明显；随着开窗角度增加换热器换热系数会呈现先增加后减小的趋势，同时压降会随开窗角度的增大而有所升高。

１ 引言

近年来，电动汽车由于不需要排放二氧化碳，不依赖化石燃料，逐渐成为汽车行业未来的发展方向。电动汽车在冬季制热时多采用 ＰＴＣ电加热器，严重影响汽车的续航里程。在低温－１０ ℃时，采用 ＰＴＣ电加热器采暖能 使 电 动 汽 车 的 续 航 里 程 下 降５０％以上，而利用热泵空调系统采暖可以使续航里程提升３５％以上［１］，因此车用热泵空调系统便应运而生。而车外换热器作为热泵空调系统的重要部分，不同于传统汽车空调室外换热器只做冷凝器，在热泵系统中，室外换热器兼顾冷凝器与蒸发器两种功能，因此提高室外微通道换热器性能的研究具有实际意义。

目前对于室外换热器研究主要集中在换热器的结霜特性［２～５］与百叶窗结构对室外换热器传热与流动特性的影响上。谷波等［６］研究了百叶窗结构变 化对微通道换热器传热和流动性能的影响，发现采用较大的百叶窗 间 距 和 翅 片 间 距 会 使 传 热 性 能 下 降．赵松田［７］分析了换热器进口空气湿度、风速和布置倾角对整体性能的影响，实验表明，迎面风速对换热器空气侧压降影响很大。张克鹏［８］研究发现翅片 开窗角度为４５°且开窗数为１４个时，微通道换热器空气侧换热量最大，压降相对较小。刘欣欣等［９］提出一种新型变截面百叶窗翅片，与传统的矩形翅片相比，综合性能因子ＪＦ 提高了７．６５％．杨凤叶等［１０］提出一种换热性能更佳的双梯形百叶窗结构，得出当双梯形百叶窗间距ＬＰ ＝１．３ｍｍ、翅 片 间 距 ＦＰ ＝１．４ｍｍ 时，综合性能最佳。徐博 等［１１］采 用ε－ＮＴＵ 法 建立了微通道换热器数值模型，计算翅片参数对传热性能的影响，结果表明换热系数与翅片间距呈负相关，与开窗角度呈正相关．ＡｒｓｌａｎＳａｌｅｅｍ 等［１２］研 究了不同翅片结构微通道换热器的空气侧热工性能研究，结果表明百 叶 窗 间 距 为１．０ ｍｍ 的 翅 片 在 换 热和压降特性方面具有最佳效果。张剑飞等［１３］研究发现微通道换热器空气阻力与换热量和迎风面积相同的平翅片圆管换热器空气阻力相当，认为独特的扁管结构和较小的换热器厚度是其减小空气阻力的有效手段。薛庆峰等［１４］研究了流道布局对微通道平行流车外换热器换热和压降性能的影响，结果表明，制冷工况下车外换热器采用１∶３的流道布局时，换热量优于其他 两 种 布 局，但 其 压 降 更 大．丁 鎏 俊 等［１５］研究了流程数、各流程扁管布置方式、扁管宽度和内孔高对室外换热器换热及压降的影响，在传热和压降的限制下应采用低流程数，适当增加扁管宽度和孔高。ＪｕｎＹｕ等［１６］采用实验和数值模拟相结合的方法对换热器外平行流动微通道进行了研究，实验结果表明，随着进气压力和进气流量的增加，传热速率增加。Ａｍｉｒｎｏｒｄｉｎ等［１７］模 拟 分 析，得出增加翅片间距时，传热系数增加的同时压降减小；而增加百叶窗间距时，传热系数和压降均减小。Ｋｉｍ 和 Ｂｕｌｌａｒｄ［１８］对比了开窗角度、翅间距对空气侧换热的影响，发现当开窗角大于 ２３°时增大翅间距换热系数会增大，且翅片间距对换热的影响会伴随雷诺数的差异而减弱。目前对于百叶窗结构对微通道换热器传热和流动性能的影响已有一定的研究，但依然不够完善。

为了提升室外换热器的流动换热性能，本文建立了百叶窗翅片换热器三维模型，采用数值模拟研究了翅片间距和开窗角度对空气侧传热和流动特性的影响，得到了两者对传热流动的影响规律，并选取最佳间距和开窗角度。

２ 数值模拟方法

２．１ 物理模型

微通道换热器翅片由铝制成，根据翅片几何结构，将所研究翅片外侧空气域进行建模．将翅片边缘处温度视为与扁管外壁面相同温度。如图１，设置计算域上下两表面为周期面，水平方向两表面为对称面，高度为一个管间距；左侧为空气速度进口，为保证来流均匀，将计算区域向上游延伸２倍翅片间距；右侧为空气压 力 出 口，相 对 压 力 为０Ｐａ，将 计 算 区域向下游延伸以防出口区域回流产生［１９］；翅 片 表面温度分布由翅片内部热传导与空气对流换热耦合求解．将室外换热器翅片结构尺寸，列于表 １ 中。

２．２ 网格划分及独立性验证

计算域网格由ＩＣＥＭ 生成，由于百叶窗翅片结构较为复杂，整体划分精度较低，所以对计算区域进行分块处理。计算区域网格划分如图２～３，固体域翅片及倒流边界部分采用结构网格，对流体域中与翅片接触部分采用四面体非结构网格，其它区域应用六面体非结构网格，最后对接触区域流固耦合处理。由于翅片表面附近流体温度速度变化比较剧烈，因此对翅片接触 面 进 行 边 界 层 处 理．边 界 层 网 格 分 三层，设置增长率为１．２，最大厚度０．０５ｍｍ，进口区域网格尺寸设置为１ｍｍ，出口区域网格尺寸为２ｍｍ，翅片与空气接触区域最大网格尺寸为０．１ｍｍ。

在进行数值模拟之前，需对网格进行独立性分析，确定模拟结果与网格数量无关．网格独立性分析方法为比较 相 同 模 型 在 不 同 网 格 数 量 下 所 求 得 的Ｃｏｌｂｕｒｎｊ因子与阻力因子ｆ 的相对差值．本文对三套不同数量网 格（１１０万、２３０万、５６０万）进 行 独 立性分析，计算结果如图４所示。

网格数量１１０万与２３０万之间 Ｃｏｌｂｕｒｎｊ 因子的最大偏差为４．８６％，平均偏差为３．９５％；阻力因子ｆ 最大偏差为１２．２３％，平均偏差为１０．９９％．网格数量２３０万与５６０万模型之间 Ｃｏｌｂｕｒｎｊ 因子的最大偏差１．５６％，平均偏差为０．９１％；阻 力 因 子ｆ 最 大偏差为１．２３％，平 均 偏 差 为 ０．８％．２３０万 网 格 数 值结果可认为独立，本文对空气侧流动传热模拟采用网格数为２３０万。

２．３ 模型可靠性验证

对比本文 模 拟 结 果 与 采 用 Ｃｈａｎｇ ＆ Ｗａｎｇ［２０］实验关联式计算结果，得出图５，由图可知，５０≤Ｒｅ≤９５０范围 内，数 值 模 拟 结 果 得 出 的 Ｃｏｌｂｕｒｎｊ 因子 与 Ｃｈａｎｇ ＆ Ｗａｎｇ 实 验 关 联 式 最 大 偏 差 为 ９．１％，平均偏差为 ７．８％；阻力因子ｆ 数值与实验关联式最大偏差为１１．５２，％，平均偏差为７．３５％．数值计算结果与实验关联式计算结果之间的偏差符合一般工程应用要求，且变化规律相同，因而本文使用模型的正确性得到验证。

３ 结果及分析

通 过 以 下 几 个 参 数 对 空 气 侧 流 动 传 热 进 行 研究：空气侧 换 热 系 数 ｈ、压 降 ΔＰ、Ｃｏｌｂｕｒｎｊ 因 子、阻力因子 ｆ 以及综合能效因子ＪＦ［２１］，其中 ＪＦ 是一个无量纲数，用以衡量相同输入功率下所能达到的换热性能，可用来评价空气侧流动换热的综合性能，其参数定义如下：

３．１ 翅片间距对流动传热的影响

对翅片深度１６ｍｍ，翅片高度５ｍｍ，开窗角度２７°，翅片间距１．０ｍｍ，１．４ｍｍ，１．８ｍｍ 进 行 数 值计算，分析翅片间距对微通道室外换热器空气侧传热和流动的影响。空气侧换热系数随进风风速Ｖｉｎ变化规律如图６所示，进风风速Ｖｉｎ一定时，空气侧换热系数随翅片间距的减小而相应增大；翅片间距一定时，空气侧换热系数随进风风速的增大而增加。

随着进风风速 的 增 大，翅 片 间 距１．４ ｍｍ 与 翅片间距１．０ｍｍ 对应的换热系数之间的差值逐渐减小，结合图７不同翅片间距下温度云图可以看出，百叶窗翅片在进风风速４ｍ／ｓ时，翅片间距的增大导致整个空气域的温度有所下降，并且在中心过渡区和出口段温度变化较为明显，其原因是百叶窗翅片间距较小时，空气主要从百叶窗通道流过，流动过程中流道长度增加，空气与翅片换热更加充分，换热效率高，因而中心平面与出口处温度较高。当翅片间距较大时，空气与百叶窗的扰动减弱，百叶窗通道流作用性减弱，管导向流的作用性增强，更多的空气流经管导流通道，因此换热效果减弱．由于百叶窗翅片与空气扰动减弱，所以Ｆｐ增大后流动阻力也应呈现下降趋势，从图８也可佐证这一观点．另外从图８还可以看 出，随 着 翅 片 间 距 的 增 加，空气侧压降逐渐减小。

如图９ 所示，雷诺数较低 时 Ｃｏｌｂｕｒｎｊ 因 子 随Ｆｐ的增大而减小，然而这种趋势随着雷诺数的增加而 减 弱，并 且 高 雷 诺 数 下 Ｆｐ ＝１．０ ｍｍ 时 的Ｃｏｌｂｕｒｎｊ 因子基本与另两组重合。这表明在低雷诺数下，较小的翅片间距更有助于强化空气与百叶窗间的换热，而当雷诺数增大后这种促进效果有所降低，此时应从其他指标因素去衡量．在空气侧阻力方面，从图９中的 阻 力 因 子ｆ 曲 线 可 以 看 出，增 大 翅片间距有助于降低空气侧的流动阻力，但当雷诺数升高后，想要降低流动阻力，单单提升 Ｆｐ效果并不明显。

不同翅片间距下ＪＦ 因子变化如图１０ 所示，在Ｒｅ＜３００时增大Ｆｐ会导致ＪＦ 因子小幅下降；当Ｒｅ＞３００ 时ＪＦ 因子随Ｆｐ的增大而升高。邓敏锋［２２］通过数值模拟，在不同迎面风速下研究得出空气侧的换热系数和压降均与翅片间距呈负相关。董军启［２３］给出了不同翅片间距 ＦＰ 时，ｊ 因 子 和ｆ 因 子 随Ｒｅ的变化曲线，发现随着翅片间距 ＦＰ 的增加 ，传热ｊ因子和摩擦ｆ 因子而降低，两者与本文研究发现相契合．综合 ＪＦ 因子与上述各参量变化趋势后，本文室外换热器翅片参数中的翅片间距取１．４ｍｍ更为适合。

３．２ 翅片开窗角度对流动传热的影响

百叶窗翅片通过对空气的扰动形成边界层，来引导空气流向，开窗角度θ 作为百叶窗翅片关键结构参数，对空气侧流动换热有着重要影响。对翅片宽度１６ｍｍ，翅片间距１．４ｍｍ 的百叶窗翅片，空气进风风速Ｖｉｎ＝１．０～６．０ｍ／ｓ，开窗角度２３°、２７°、３０°进行数值计算，分析了翅片开窗角度对微通道室外换热器传热和流动 的 影 响．不 同 开 窗 角 度θ 下 空 气 侧的换热系数和压降的变化趋势如图 １１所 示，当 开窗角度一定时，换热系数与压降均随进风风速的增大而升高．对比相同进风风速下不同开窗角度曲线，发现增加开窗角度后换热系数ｈ 会 呈 现 先 增 加 后减小的趋势，可以判断出存在一个最佳开窗角度值，在该值下百叶窗翅片空气侧换热最优。

该现象符合 Ｋａｊｉｎｏ［２４］提出的边界层冲击理论，该理论表示百叶窗类翅片强化换热机理在于每一个百叶窗翅片前端形成的薄边界层，而这类边界层通常在气流的影响下是不平衡的，具体表现在百叶窗上下表面的换热能力上，百叶窗总换热能力是上下表面换热能力之和，在某个环境条件下存在阈值。结合图１２温度场进行分析，发现开窗角度θ的升高整个空气域温度略微下降．设定来流方向为正方向，在进口迎风面位置（图示 Ａ 区域），随着开窗角度θ 的升高，百叶窗正面温度有所下降，百叶窗背面温度呈现上升趋势，这是由于翅片角度增大后，空气中百叶窗导向流所占比例随之增大，随后达到临界值，到达临界值之后，空气侧流动方式管导向流所占比例增加，增强了百叶窗翅片背面换热，因而导致了百叶窗翅片正背面温度差异。

由图１１与图１３ 可 以 看 出，空 气 侧 压 降 Δｐ 会随开窗角度θ的增大而有所升高．开窗角度θ 由２３°增大到２７°时，空 气 侧 压 降 变 化 相 较 于２７°到３０°比较明显，这 主 要 是 由 于２７°开 窗 角 更 接 近 百 叶 窗 导向流与管导向流的临界平衡值，压降变化比较明显．在开窗角度达 到 ３０°后，空 气 的 流 动 方 式 主 要 为 管导向流，θ对 Δｐ 的影响有所下降．如图 １３ 所示，在较低雷诺数时，随着开窗角度的增加，Ｃｏｌｂｕｒｎｊ 因子变化不是很明显，在较高雷诺数时，随着开窗角度的增加，Ｃｏｌｂｕｒｎｊ 因 子 有 明 显 的 上 升，并 且２３°到２７°的增幅明显高于２７°到３０°范围，这与开窗角度对换热系数影响趋势相同。由 图１２与 图 １３分 析 开 窗角对空气侧流动的影响，可以看出开窗角度的增加，阻力因子ｆ 呈现上升趋势，即增加开窗角度不利于空气侧的流动 换 热，随 后 根 据 Ｃｏｌｂｕｒｎｊ 因 子 和 阻力因子ｆ 数值计算出ＪＦ 因子，结果如图１４所示。

从图１４可以看出，随着雷诺数的增加，综合评价ＪＦ 因子呈现下降趋势；对于相同雷诺数情况，开窗角度θ为２３°与２７°时的综合评价因子高于３０°．然而对于较小开窗角下 的ＪＦ 因 子 评 价，则 要 根 据 具体雷诺数进行讨论，当Ｒｅ＜３００时，ＪＦ 因子随着开窗角度的升高有所下降，说明增大开窗角度对不利于百叶窗翅片空气侧综合性能的提高，即车辆日常车速较低，室外换热器迎面风速较低时，百叶窗翅片开窗建议选取较低开窗角度．当 Ｒｅ＞３００ 时，θ＝２７°较优，并且３０°下ＪＦ 因子降幅减缓，说明较高雷诺数时，增大开窗角度有利于提升空气侧流动换热的综合性能，较大的翅片开窗角更适合迎面风速更大的情形。

４ 结　论

（１）当进风风速一定时，空气侧换热系数与随翅片间距的减小而相应增加；翅片间距一定时，空气侧换热系数随进 风 风 速 的 增 大 而 增 加．随 着 进 风 风 速的增大，翅片间距１．４ｍｍ 与翅片间距１．０ｍｍ 对应的换热系数之间的差值逐渐减小。

（２）在Ｒｅ＜３００时增大Ｆｐ会导致ＪＦ 因子小幅下降；当Ｒｅ＞３００时ＪＦ 因子随Ｆｐ的增大而升高，进风风速越高时，增大翅片间距对综合性能更有利，因此百叶窗翅片间距取１．４ｍｍ。

（３）当Ｒｅ＜３００时，ＪＦ 因子随着开窗角度的升高有所下降，说明增大开窗角度对不利于百叶窗翅片空气侧综合 性 能 的 提 高；当 Ｒｅ＞３００ 时，θ＝２７°较优，并且３０°下ＪＦ 因子降幅减缓，说明较高雷诺数时，增大开窗角度有利于提升空气侧流动换热的综合性能，因此翅片开窗角度取２７°更优。