0 引言
20世纪90年代,LORENTZEN等[3]率先提出将跨临界CO2空调系统运用于汽车空调,随后大量学者针对CO2制冷剂在汽车空调系统中的应用展开了广泛研究。刘洪胜等[4]研究了CO2充注量对汽车空调性能的影响,发现存在使系统性能系数（coefficient of performance, COP）最大的CO2充注量。KAUF[5]和LIAO等[6]通过数值模拟研究发现存在最优排气压力使得系统COP最大,并分析了影响最优排气压力的因素。KIM等[7]对跨临界CO2汽车空调进行了实验研究,发现在最佳气体冷却器压力下,系统在变工况下均具有较好的制冷性能。对于跨临界CO2循环的优化,CHO等[8]和BAEK等[9]研究了带中间补气的跨临界CO2循环,发现带补气的跨临界CO2循环拥有更高的制冷、制热性能。董丽玮等[10]研究了带引射器的跨临界CO2机械过冷热泵系统,发现引射器可以减少系统能量损失,提升系统COP。方健珉等[11]研究了回热器对汽车跨临界CO2系统的影响,结果表明回热器的使用可以有效提高系统性能,同时会引起压缩机压比增大和排气温度的升高。
汽车空调不同于家用空调系统,由于其空间可移动性决定了系统运行工况处于不断变化中,跨临界CO2空调系统的制冷性能受外界环境影响较大,当制冷性能不佳时将导致电动汽车续航里程下降。本文根据电动汽车空调系统构建跨临界CO2循环模型,通过仿真研究变工况下系统的运行性能。利用仿真数据对车用跨临界CO2空调系统及其系统部件进行㶲分析,探讨变压缩机转速、变气体冷却器进风风速、变气体冷却器进风温度、变蒸发器进风量和变蒸发器进风温度对系统㶲损和㶲效率及各部件相对㶲损的影响,以期为电动汽车跨临界CO2空调系统的优化和设计提供参考。
1 模型建立
Fig. 1 Diagram of trans-critical CO2 cycle heat pump system图1 跨临界CO2循环热泵系统简图 |
Fig. 2 Diagram of trans-critical CO2 cycle图2 跨临界CO2循环图 |
1.1 压缩机模型
压缩机模型由等熵效率（ηis）、容积效率（ηv）、制冷剂流量（${{\dot{m}}_{\text{r}}}$）和压缩机耗功（Wcomp）所构成：
${{\eta }_{\text{is}}}=\text{0}\text{.9611}-\text{0}\text{.052}\frac{{{p}_{\text{dis}}}}{{{p}_{\text{suc}}}}$ (1)
${{\eta }_{\text{v}}}=\text{0}\text{.8882}-\text{0}\text{.04185}\frac{{{p}_{\text{dis}}}}{{{p}_{\text{suc}}}}$ (2)
${{\dot{m}}_{\text{r}}}={{\rho }_{1}}{{\eta }_{\text{v}}}V{{n}_{\text{com}}}$ (3)
${{W}_{\text{comp}}}=\frac{{{{\dot{m}}}_{\text{r}}}\left( {{h}_{2}}-{{h}_{1}} \right)}{{{\eta }_{\text{is}}}{{\eta }_{\text{me}}}}$ (4)
式中：pdis、psuc分别为压缩机的排气压力和吸气压力,MPa;${{\dot{m}}_{\text{r}}}$为制冷剂质量流量,kg/s;${{\rho }_{1}}$为压缩机入口制冷剂密度,kg/m3;V为压缩机的排量,cm3;ncom为压缩机的转速,r/min;Wcomp为压缩机耗功,kW;h1、h2分别为压缩机入口和出口的制冷剂焓值,kJ/kg;ηme为压缩机机械效率,取0.9。
1.2 换热器模型
对换热器中制冷剂的流动和换热进行了简化,采用离散微元（空气/壁面、壁面和翅片以及制冷剂/壁面）进行建模。
对于换热器内部的制冷剂侧：
${{Q}_{\text{int,}i}}={{h}_{\text{int,}i}}{{A}_{\text{r,}i}}\left( {{T}_{\text{r,}i}}-{{T}_{\text{w,}i}} \right)+{{\dot{m}}_{\text{r}}}{{L}_{\text{vap}}}$ (5)
${{h}_{\text{int}}}=N{{u}_{\text{r}}}\frac{{{\lambda }_{\text{r}}}}{{{D}_{\text{r}}}}$ (6)
$N{{u}_{\text{r}}}=\frac{\left( f/\text{8} \right)\left( Re-\text{1000} \right)Pr}{\text{1}+\text{12}\text{.7}\sqrt{f/\text{8}}\left( P{{r}^{2/3}}-\text{1} \right)}$ (7)
$f=8{{\left\{ {{\left( \frac{8}{Re} \right)}^{12}}+{{\left[ {{\left\{ 2.45\ln \left[ \frac{1}{{{\left( {7}/{Re}\; \right)}^{0.9}}+0.27{{\varepsilon }_{0}}} \right] \right\}}^{16}}+{{\left( \frac{37530}{Re} \right)}^{16}} \right]}^{-{3}/{2}\;}} \right\}}^{{1}/{12}\;}}$ (8)
式中：Qint,i为换热器内部第i个微元段制冷剂的换热量,kW;hint,i为换热器第i个微元段内制冷剂与管壁的对流换热系数,kJ/(m2∙K);Ar,i为换热器内部第i个微元段制冷剂的换热面积,m2;Tr,i为第i个微元段内制冷剂的温度,℃;Tw,i第i个微元段内换热器的壁温,℃;Lvap为微元段内的制冷剂蒸发潜热,kJ/kg;hint为换热器内部制冷剂与管壁的对流换热系数,kJ/(m2∙K);Nur为制冷剂的努塞尔数;λr为制冷剂的导热系数,kJ/(m∙K);Dr为换热器制冷剂侧特征长度,m;Re为雷诺数;Pr为普朗特数;f为阻力系数;ε0为管壁的相对粗糙度,取5 μm。
对于换热器外部的空气侧：
${{Q}_{\text{ext,}i}}={{\eta }_{\text{fin,}i}}{{h}_{\text{ext,}i}}{{A}_{\text{a,}i}}\left( {{T}_{\text{w,}i}}-{{T}_{\text{a,}i}} \right)$ (9)
${{h}_{\text{ext}}}={{\lambda }_{\text{a}}}N{{u}_{\text{a}}}/{{D}_{\text{a}}}$ (10)
$N{{u}_{\text{a}}}={{\eta }_{\text{fin}}}R{{e}^{\text{0}\text{.8}}}P{{r}^{\text{0}\text{.3}}}$ (11)
${{\eta }_{\text{fin}}}={{\lambda }_{\text{w}}}/{{e}^{\text{*}}}$ (12)
式中：Qext,i为换热器外部第i个微元段空气的换热量,kW;ηfin,i为翅片第i个微元段的肋效率;hext,i为换热器第i个微元段内空气与管壁的对流换热系数,kJ/(m2∙K);Aa,i为换热器外部第i个微元段空气的换热面积,m2;Ta,i为第i个微元段内空气的温度,℃;hext换热器外部空气与管壁的对流换热系数,kJ/(m2∙K);λa为空气的导热系数,kJ/(m∙K);Da为换热器空气侧特征长度,m;Nua为空气的努塞尔数;ηfin为翅片的肋效率;λw为管壁的导热系数,kJ/(m∙K);e* 是几何参数。
换热器壁温：
$\frac{\text{d}{{T}_{\text{w,}i}}}{\text{d}t}=\frac{{{Q}_{\text{int,}i}}+{{Q}_{\text{ext,}i}}}{{{m}_{i}}{{C}_{p\text{,}i}}}$ (13)
式中：mi为换热器第i个微元段质量,kg;Cp,i为换热器第i个微元段比热容,kJ/(kg∙K)。
1.3 膨胀阀模型
制冷剂经过膨胀阀是等焓节流过程,制冷剂质量流量计算如下：
${{\dot{m}}_{\text{r}}}={{C}_{\text{D}}}{{A}_{\text{C}}}\sqrt{\text{2}{{\rho }_{\text{exp}}}\text{ }\!\!\Delta\!\!\text{ }p}$ (14)
式中：CD为制冷剂通过膨胀阀的流量系数;AC为膨胀阀的流通截面积,m2;ρexp为膨胀阀入口处制冷剂的密度,kg/m3;∆p为膨胀阀进出口的压差,MPa。
1.4 系统仿真平台
Table 1 Parameters of trans-critical CO2 air conditioning system表1 跨临界CO2空调系统参数 |
| 系统部件 | 参数 | 数值 |
|---|---|---|
| 压缩机 | 排量/(cm3/r) | 8 |
| 转速/(r/min) | 2 000 ~ 8 800 | |
| 气冷器 | 尺寸/mm3 | 500 × 364 × 25 |
| 流程 | 10, 8, 6, 4 | |
| 空气迎风面积/mm2 | 137 200 | |
| 制冷剂对流换热面积/mm2 | 420 000 | |
| 蒸发器 | 尺寸/mm3 | 250 × 232 × 30 |
| 流程 | 8, 6, 6 | |
| 空气迎风面积/mm2 | 37 025 | |
| 制冷剂对流换热面积/mm2 | 300 000 |
Table 2 Simulation conditions of trans-critical CO2 air conditioning system表2 跨临界CO2空调系统仿真工况 |
| 工况 | 压缩机转速/(r/min) | 室外气冷器迎面风速/(m/s) | 室外气冷器进风温度/℃ | 蒸发器进风量/(m3/h) | 蒸发器进风温度/℃ |
|---|---|---|---|---|---|
| 1 | 变量 | 5 | 35 | 700 | 25 |
| 2 | 3 000 | 变量 | 35 | 700 | 25 |
| 3 | 3 000 | 5 | 变量 | 700 | 25 |
| 4 | 3 000 | 5 | 35 | 变量 | 25 |
| 5 | 3 000 | 5 | 35 | 700 | 变量 |
2 㶲分析
在没有核效应、磁、电和表面张力的情况下,物质流的㶲可表达为：
$\dot{E}={{\dot{E}}_{\text{k}}}+{{\dot{E}}_{\text{p}}}+{{\dot{E}}_{\text{ph}}}+{{\dot{E}}_{\text{ch}}}$ (15)
式中：${{\dot{E}}_{\text{k}}}$为动能㶲;${{\dot{E}}_{\text{p}}}$为势能㶲;${{\dot{E}}_{\text{ph}}}$为物理㶲;${{\dot{E}}_{\text{ch}}}$为化学㶲。由于动能和势能的变化在制冷空调系统内造成的影响较小,因此在热泵系统内可忽略动能㶲和势能㶲。化学㶲与系统的化学成分偏离其化学平衡有关,制冷过程不涉及化学物质的混合和化学反应,化学㶲可忽略不计。因此在整个制冷系统中只需要考虑其物理㶲。
对于单一工质的物理㶲$\psi $可表达为：
$\psi =\left( h-{{h}_{\text{0}}} \right)-{{T}_{\text{0}}}\left( s-{{s}_{\text{0}}} \right)$ (16)
式中：T0为环境温度;h和s为该工况点下工质的焓和熵;h0和s0为环境工况下工质的焓和熵。
系统内的进出口流动㶲的差异为系统的不可逆损失,即㶲损,可由㶲平衡方程表示：
$\text{ }\!\!\Delta\!\!\text{ }\dot{E}=\sum {{\dot{E}}_{\text{in}}}-\sum {{\dot{E}}_{\text{out}}}$ (17)
式中：${{\dot{E}}_{\text{in}}}$和${{\dot{E}}_{\text{out}}}$为通过系统进出口的质量和功所传递的总流动㶲,将公式（17）展开可得：
$\begin{align} & \Delta \dot{E}=\sum\nolimits_{\text{in}}{{{{\dot{m}}}_{\text{r}}}\psi }-\sum\nolimits_{\text{out}}{{{{\dot{m}}}_{\text{r}}}\psi }+\sum\nolimits_{\text{in}}{{{{\dot{m}}}_{\text{a}}}\psi }- \\ & \ \ \ \ \ \ \ \ \sum\nolimits_{\text{out}}{{{{\dot{m}}}_{\text{a}}}\psi }+\sum\nolimits_{\text{in}}{{\dot{W}}}-\sum\nolimits_{\text{out}}{{\dot{W}}} \\ \end{align}$ (18)
式中：$\psi $为对应工质的物理㶲,kJ/kg;${{\dot{m}}_{\text{a}}}$为空气质量流量,kg/m3;$\dot{W}$为流体的流动功,kW。
根据表达式（18）和热力学第一定律,可得跨临界CO2系统各部件的㶲损失。
蒸发器：
$\begin{aligned} &\Delta\dot{E}_{\mathrm{e}}&& =\dot{m}_\mathrm{r}\left(\psi_4-\psi_1\right)+\dot{m}_\mathrm{ae}\left(\psi_\mathrm{in}-\psi_\mathrm{out}\right) \\ &&&=\begin{cases}\dot{m}_\mathrm{r}\left[\left(h_4-h_1\right)-T_0\left(s_4-s_1\right)\right]+\\\dot{m}_\mathrm{ae}\left[\left(h_\mathrm{aei}-h_\mathrm{aeo}\right)-T_0\left(s_\mathrm{aei}-s_\mathrm{aeo}\right)\right]\end{cases} \end{aligned}$ (19)
式中：${{\dot{m}}_{\text{ae}}}$为蒸发器中空气质量流量,kg/m3;h4为蒸发器入口处制冷剂焓值,kJ/kg;h1为蒸发器出口处制冷剂焓值,kJ/kg;s4为蒸发器入口处制冷剂熵值,kJ/(kg∙K);s1为蒸发器出口处制冷剂熵值,kJ/(kg∙K);haei为蒸发器入口处空气焓值,kJ/kg;haeo为蒸发器出口处空气焓值,kJ/kg;saei为蒸发器入口处空气熵值,kJ/(kg∙K);saeo为蒸发器出口处空气熵值,kJ/(kg∙K)。
压缩机：
$\text{ }\!\!\Delta\!\!\text{ }{{\dot{E}}_{\text{comp}}}={{\dot{m}}_{\text{r}}}\left( {{\psi }_{\text{1}}}-{{\psi }_{\text{2}}} \right)+{{W}_{\text{comp}}}$ (20)
压缩机耗功：
${{W}_{\text{comp}}}={{\dot{m}}_{\text{r}}}\left( {{h}_{\text{2}}}-{{h}_{\text{1}}} \right)$ (21)
公式（20）可表达为：
$\begin{align} & \text{ }\!\!\Delta\!\!\text{ }{{{\dot{E}}}_{\text{comp}}}={{{\dot{m}}}_{\text{r}}}\left[ \left( {{h}_{\text{1}}}-{{h}_{\text{2}}} \right)-{{T}_{\text{0}}}\left( {{s}_{\text{1}}}-{{s}_{\text{2}}} \right) \right]+{{{\dot{m}}}_{\text{r}}}\left( {{h}_{\text{2}}}-{{h}_{\text{1}}} \right) \\ & \ \ \ \ \ \ \ \ \ \ ={{T}_{\text{0}}}\left( {{s}_{\text{2}}}-{{s}_{\text{1}}} \right) \\ \end{align}$ (22)
气冷器：
$\begin{aligned} &\Delta\dot{E}_{\mathrm{g}}&& =\dot{m}_\mathrm{r}\left(\psi_2-\psi_3\right)+\dot{m}_\mathrm{ag}\left(\psi_\mathrm{in}-\psi_\mathrm{out}\right) \\ &&&=\begin{cases}\dot{m}_\mathrm{r}\left[\left(h_2-h_3\right)-T_0\left(s_2-s_3\right)\right]+\\\dot{m}_\mathrm{ag}\left[\left(h_\mathrm{agi}-h_\mathrm{ago}\right)-T_0\left(s_\mathrm{agi}-s_\mathrm{ago}\right)\right]\end{cases} \end{aligned}$ (23)
式中：${{\dot{m}}_{\text{ag}}}$为气冷器中空气质量流量,kg/m3;h2为气冷器入口处制冷剂焓值,kJ/kg;h3为气冷器出口处制冷剂焓值,kJ/kg;s2为气冷器入口处制冷剂熵值,kJ/(kg∙K);s3为气冷器出口处制冷剂熵值,kJ/(kg·K);hagi为气冷器入口处空气焓值,kJ/kg;hago为气冷器出口处空气焓值,kJ/kg;sagi为气冷器入口处空气熵值,kJ/(kg∙K);sago为气冷器出口处空气熵值,kJ/(kg∙K)。
膨胀阀：
$\text{ }\!\!\Delta\!\!\text{ }{{\dot{E}}_{\text{v}}}={{\dot{m}}_{\text{r}}}\left[ \left( {{h}_{\text{3}}}-{{h}_{\text{4}}} \right)-{{T}_{\text{0}}}\left( {{s}_{\text{3}}}-{{s}_{\text{4}}} \right) \right]$ (24)
膨胀阀为等焓节流过程,其中：
${{h}_{\text{3}}}={{h}_{\text{4}}}$ (25)
公式（24）可表达为：
$\text{ }\!\!\Delta\!\!\text{ }{{\dot{E}}_{\text{v}}}={{T}_{\text{0}}}\left( {{s}_{\text{4}}}-{{s}_{\text{3}}} \right)$ (26)
系统总㶲损为：
$\text{ }\!\!\Delta\!\!\text{ }{{\dot{E}}_{\text{tot}}}=\text{ }\!\!\Delta\!\!\text{ }{{\dot{E}}_{\text{e}}}+\text{ }\!\!\Delta\!\!\text{ }{{\dot{E}}_{\text{comp}}}+\text{ }\!\!\Delta\!\!\text{ }{{\dot{E}}_{\text{g}}}+\text{ }\!\!\Delta\!\!\text{ }{{\dot{E}}_{\text{v}}}$ (27)
系统各部件的相对㶲损为：
${{\delta }_{i}}=\frac{\text{ }\!\!\Delta\!\!\text{ }{{{\dot{E}}}_{i}}}{\text{ }\!\!\Delta\!\!\text{ }{{{\dot{E}}}_{\text{tot}}}}\times 100\%$ (28)
跨临界CO2循环的制冷㶲效率表示为：
$\begin{aligned} &\eta_{\mathrm{ex}}&& =\frac{\dot{E}_{\mathrm{gain}}}{\dot{E}_{\mathrm{pay}}}\times100\%=\frac{\left|\dot{m}_{\mathrm{ae}}\left(\psi_{\mathrm{out}}-\psi_{\mathrm{in}}\right)\right|}{\dot{m}_{\mathrm{r}}\left(h_{2}-h_{1}\right)}\times100\% \\ &&&=\frac{\left|\dot{m}_{\mathrm{ae}}\left[\left(h_{\mathrm{aeo}}-h_{\mathrm{aei}}\right)-T_{0}\left(s_{\mathrm{aeo}}-s_{\mathrm{aei}}\right)\right]\right]}{\dot{m}_{\mathrm{r}}\left(h_{2}-h_{1}\right)}\times100\% \end{aligned}$ (29)
式中：Egain、Epay分别为跨临界CO2循环制冷的收益㶲和投入㶲。
3 㶲分析结果与讨论
3.1 压缩机转速的影响
Fig. 3 The effects of compressor speed on the exergy loss of system components图3 压缩机转速对系统各部件㶲损值的影响 |
Fig. 4 The effects of compressor speed on the relative exergy loss of system components图4 压缩机转速对系统各部件相对㶲损的影响 |
由图4可知,随着压缩机转速的增加,气冷器和压缩机的相对㶲损增大,膨胀阀和蒸发器的相对㶲损减小。压缩机转速较低时,膨胀阀的相对㶲损最大,压缩机的相对㶲损最小,当压缩机转速为3 000 r/min时,膨胀阀、气冷器、蒸发器和压缩机的相对㶲损分别是37.24%、23.4%、20.37%和18.97%。随着压缩机转速的增加,气冷器的相对㶲损最大,蒸发器的相对㶲损最小,当压缩机转速为8 000 r/min时,气冷器、压缩机、膨胀阀和蒸发器的相对㶲损分别为29.71%、27.76%、26.11%和16.42%。表3展现了不同的压缩机转速下各部件的平均相对㶲损。由表可知,膨胀阀的平均相对㶲损最大,平均占比31.10%,其次是气冷器和压缩机,平均占比分别为26.93%和23.47%,蒸发器的相对㶲损最小,平均占比18.49%。
Table 3 Average relative exergy loss of system components at different compressor speeds表3 不同的压缩机转速下系统各部件平均相对㶲损 |
| 系统部件 | 平均相对㶲损/% |
|---|---|
| 蒸发器 | 18.49 |
| 气冷器 | 26.93 |
| 压缩机 | 23.47 |
| 膨胀阀 | 31.10 |
图5展示了压缩机转速对系统总㶲损和㶲效率的影响。由图可知,随着压缩机转速增加,系统的总㶲损逐渐大幅度增大,㶲效率逐渐减小。当压缩机转速从3 000 r/min增加至8 000 r/min时,系统总㶲损从1.22 kW增大至4.25 kW,㶲效率从15.85%减小至11.13%。这是由于随着压缩机转速的增加,压缩机的排气温度逐渐升高,压缩机和气冷器内由于与环境之间的换热造成的不可逆损失逐渐增加,㶲损逐渐增加,其相对㶲损逐渐增大,而膨胀阀和蒸发器的温度更接近环境温度,压缩机转速变化导致的温度变化相对于压缩机和气冷器较小,因此膨胀阀和蒸发器的㶲损增加幅度相对较小,其相对㶲损逐渐减小。由于压缩机转速的提升,系统的制冷剂流量大幅增加,系统总㶲损增加。同时,系统获得的冷量㶲增加,但压缩机功耗提升更显著,因此系统㶲效率随压缩机转速增加而降低。
Fig. 5 The effects of compressor speed on the exergy efficiency and total exergy loss of system图5 压缩机转速对㶲效率和总㶲损的影响 |
3.2 气体冷却器进风风速的影响
Fig. 6 The effects of air inlet speed of gas cooler on the exergy loss of system components图6 气冷器进风风速对系统各部件㶲损值的影响 |
Fig. 7 The effects of air inlet speed of gas cooler on the relative exergy loss of system components图7 气冷器进风风速对系统各部件相对㶲损的影响 |
由图7可知,随着气冷器进风风速的增加,气冷器、压缩机和膨胀阀的相对㶲损先减小,后趋于平稳,蒸发器的相对㶲损先增加后趋于平稳。气冷器进风风速较低时,膨胀阀的相对㶲损最大,蒸发器的相对㶲损最小,当气冷器进风风速为1 m/s时,膨胀阀、气冷器、压缩机和蒸发器的相对㶲损分别为41.36%、25.86%、22.22%和10.57%。随着气冷器进风风速的增加,膨胀阀的相对㶲损仍最大,压缩机的相对㶲损最小。当气冷器进风风速为5 m/s时,膨胀阀、气冷器、蒸发器和压缩机的相对㶲损分别为37.24%、23.39%、20.42%和18.96%。表4展现了不同的气冷器进风风速下各部件的平均相对㶲损,由表可知,膨胀阀的相对㶲损最大,平均占比38.39%,其次是气冷器和压缩机,平均占比分别为24.10%和20.10%,蒸发器的相对㶲损最小,平均占比17.42%。
Table 4 Average relative exergy loss of system components at different air inlet speed of gas cooler表4 不同的气冷器进风风速下系统各部件平均相对㶲损 |
| 系统部件 | 平均相对㶲损/% |
|---|---|
| 蒸发器 | 17.42 |
| 气冷器 | 24.10 |
| 压缩机 | 20.10 |
| 膨胀阀 | 38.39 |
图8展示了气冷器进风风速对系统总㶲损和㶲效率的影响。由图可知,随着气冷器进风风速的增加,总㶲损先快速减小后趋于平缓,㶲效率先快速增大后趋于平缓。当气冷器风速从1 m/s增加到5 m/s时,总㶲损从1.71 kW减小到1.22 kW,㶲效率从9.14%增加到15.85%。这是由于随着气冷器进风风速的增加,气冷器对处于超临界状态的CO2制冷剂冷却效果越显著,气冷器出口的制冷剂温度降低,为维持系统设定的5 ℃过热度,膨胀阀开度增大,压缩机的压比减小。压缩机和膨胀阀内制冷剂的平均温度减小,与环境换热造成的不可逆损失减小,相应的㶲损减小。气冷器内的制冷剂平均温度减小,而气冷器外的空气温度略有降低,但由于空气温度降幅较小且与环境温度相差较小,综合两者因素,气冷器由于与环境换热的不可逆损失减小,相应㶲损减小。蒸发器内制冷剂温度减小,增大了与空气的换热温差,因此蒸发器与空气换热的不可逆损失增大,相应㶲损增大。虽然蒸发器的㶲损随着气冷器出风温度的升高而增大,但压缩机、膨胀阀和气冷器的㶲损减小,在各部件的综合影响下,系统的总㶲损减小。随着蒸发器空气温度的降低,系统获得的㶲增加,系统付出的㶲即压缩机功耗减小,系统的㶲效率增大。由于气冷器体积的确定,其最大换热能力受限,当接近最大换热能力时,继续增加气冷器进风风速不能使气冷器内的制冷剂降温,系统的运行状态变化较小,因此当气冷器进风风速继续增加时,总㶲损和㶲效率变化逐渐缓慢,直至相对平缓。
Fig. 8 The effects of air inlet speed of gas cooler on the exergy efficiency and total exergy loss of system图8 气冷器进风风速对总㶲损和㶲效率的影响 |
3.3 气体冷却器进风温度的影响
Fig. 9 The effects of air inlet temperature of gas cooler on the exergy loss of system components图9 气冷器进风温度对系统各部件㶲损值的影响 |
Fig. 10 The effects of air inlet temperature of gas cooler on the relative exergy loss of system components图10 气冷器进风温度对系统各部件相对㶲损的影响 |
由图10可知,随着气冷器进风温度的升高,膨胀阀和压缩机的相对㶲损增加,气冷器和蒸发器的相对㶲损减小。气冷器进风温度较低时,膨胀阀的相对㶲损最大,压缩机的相对㶲损最小,当气冷器进风温度为30 ℃时,膨胀阀、气冷器、蒸发器和压缩机的相对㶲损分别为33.16%、25.22%、24.75%和17.86%。随着气冷器进风温度的升高,膨胀阀的相对㶲损仍最大,蒸发器的相对㶲损较小,当气冷器进风温度为46 ℃时,膨胀阀、压缩机、气冷器和蒸发器的相对㶲损分别为46.32%、20.61%、20.23%和12.84%。表5展现了不同的气冷器进风温度下各部件的平均相对㶲损,由表可知,膨胀阀的相对㶲损最大,平均占比39.73%,其次是气冷器和压缩机,平均占比分别为22.56%和19.33%,蒸发器的相对㶲损最小,平均占比18.38%。
Table 5 Average relative exergy loss of system components at different air inlet temperature of gas cooler表5 不同的气冷器进风温度下系统各部件平均相对㶲损 |
| 系统部件 | 平均相对㶲损/% |
|---|---|
| 蒸发器 | 18.38 |
| 气冷器 | 22.56 |
| 压缩机 | 19.33 |
| 膨胀阀 | 39.73 |
图11展示了气冷器进风温度对系统总㶲损和㶲效率的影响。由图可知,随着气冷器进风温度的升高,总㶲损和㶲效率逐渐增大。当气冷器进风温度从30 ℃升高至46 ℃时,总㶲损从1.12 kW增加到1.46 kW,㶲效率从14.02%增加到17.52%。这是由于随着气冷器进风温度的升高,气冷器出口的制冷剂温度也随之升高,为维持系统设定的5 ℃过热度,膨胀阀开度减小,压缩机的压比增大。压缩机和膨胀阀内制冷剂的平均温度升高,与环境温度温差增大,不可逆损失增加,由于膨胀阀开度减小,膨胀阀内的损耗更大,相应的㶲损增幅更显著。气冷器内制冷剂的平均温度升高,但是由于环境温度（即气冷器进风温度）同时升高,因此气冷器内制冷剂和外侧空气的平均温度与环境温度的温差上升较低,㶲损增加较小,相较于膨胀阀和压缩机的增长,其相对㶲损减小。蒸发器内制冷剂与外侧空气的平均温度升高,换热温差减小,不可逆损失减小,相应的㶲损降低。在各部件的综合影响下,系统的总㶲损随气冷器进风温度的升高而增大。随着气冷器进风温度的升高,系统的制冷性能减弱,压缩机耗功增加,但气冷器温度的升高意味着环境温度也升高,环境温度的升高会导致系统获得的冷量㶲增加,综合影响下,系统㶲效率随着气冷器进风温度的升高而增大。
Fig. 11 The effects of air inlet temperature of gas cooler on the exergy efficiency and total exergy loss of system图11 气冷器进风温度对总㶲损和㶲效率的影响 |
3.4 蒸发器进风量的影响
Fig. 12 The effects of air inlet volume of evaporator on the exergy loss of system components图12 蒸发器进风量对系统各部件㶲损值的影响 |
Fig. 13 The effects of air inlet volume of evaporator on the relative exergy loss of system components图13 蒸发器进风量对系统各部件相对㶲损的影响 |
由图13可知,随着蒸发器进风量的增加,膨胀阀的相对㶲损先增大后趋于平稳,气冷器和压缩机的相对㶲损减小,蒸发器的相对㶲损先减小后增大。蒸发器进风量较低时,膨胀阀的相对㶲损最大,蒸发器的相对㶲损最小,当蒸发器进风量为300 m3/h时,膨胀阀、气冷器、压缩机和蒸发器的相对㶲损分别为33.46%、24.68%、21.79%和20.07%。随着蒸发器进风量的增加,膨胀阀的相对㶲损仍最大,压缩机的相对㶲损最小,当蒸发器进风量为700 m3/h时,膨胀阀、气冷器、蒸发器和压缩机的相对㶲损分别为37.24%、23.39%、20.42%和18.96%。表6展示了不同的蒸发器进风量下各部件的平均相对㶲损失,由表可知,膨胀阀的相对㶲损最大,平均占比36.33%,其次是气冷器和压缩机,平均占比分别为23.93%和20.08%,蒸发器的相对㶲损最小,平均占比19.67%。
Table 6 Average relative exergy loss of system components at different air inlet volume of evaporator表6 不同的蒸发器进风量下系统各部件平均相对㶲损 |
| 系统部件 | 平均相对㶲损/% |
|---|---|
| 蒸发器 | 19.67 |
| 气冷器 | 23.93 |
| 压缩机 | 20.08 |
| 膨胀阀 | 36.33 |
图14展现了蒸发器进风量对总㶲损和㶲效率的影响。由图可知,随着蒸发器进风量的增加,总㶲损先减小后增大,㶲效率先增大后减小。当蒸发器进风量从300 m3/h增加至700 m3/h时,总㶲损先从1.268 kW减小到1.206 kW,再增加到1.221 kW,㶲效率先从16.32%增加到17.53%,再减小到15.85%。这是由于随着蒸发器进风量的增加,蒸发器的换热能力逐渐增强,为维持系统设定的5 ℃过热度,膨胀阀开度增大,压比减小,蒸发温度上升,制冷剂流量增加。压缩机、气冷器内制冷剂的平均温度降低,换热造成的不可逆损失减小,相应的㶲损减小。虽然膨胀阀进出口制冷剂的温度逐渐接近环境温度,但是膨胀阀由换热温差造成的不可逆损失影响较小,由于制冷剂质量流量的增加,膨胀阀的㶲损在综合影响下增加。蒸发器内制冷剂的温度逐渐升高,外侧的空气温度逐渐升高,换热温差逐渐缩小,蒸发器的㶲损减小,但蒸发器的换热能力受限,蒸发器进风量增加到一定量时,蒸发温度变化较小,外侧空气的出风温度逐渐升高,制冷剂与空气的温差增大,且制冷剂的质量流量增加,综合影响下蒸发器的㶲损当进风量到达一定值后逐渐增大。由于蒸发器的换热能力受限,系统的总㶲损整体呈现先减小后增大。蒸发器外侧的空气温度逐渐升高,单位质量下系统获得的冷量㶲减少,压缩机压比减小,压缩机耗功减小,综合影响下,蒸发器的㶲效率随着蒸发器进风温度的升高呈现先增加后减小。
Fig. 14 The effects of air inlet volume of evaporator on the exergy efficiency and total exergy loss of system图14 蒸发器进风量对总㶲损和㶲效率的影响 |
3.5 蒸发器进风温度的影响
Fig. 15 The effects of air inlet temperature of evaporator on the exergy loss of system components图15 蒸发器进风温度对系统各部件㶲损值的影响 |
由图16可知,随着蒸发器进风温度的升高,膨胀阀和蒸发器的相对㶲损增加,压缩机和气冷器的相对㶲损减少。蒸发器进风温度较低时,膨胀阀的相对㶲损最大,蒸发器的相对㶲损最小。当蒸发器进风温度为20 ℃时,膨胀阀、气冷器、压缩机和蒸发器的相对㶲损分别为35.67%、24.77%、20.88%和18.69%。随着蒸发器进风温度的升高,膨胀阀的相对㶲损仍最大,压缩机的相对㶲损最小。当蒸发器进风温度为28 ℃时,膨胀阀、气冷器、蒸发器和压缩机的相对㶲损分别为38.50%、22.41%、21.27%和17.82%。表7展现了不同的蒸发器进风温度下各部件的平均相对㶲损。由表可知,膨胀阀的相对㶲损最大,平均占比37.00%,其次是气冷器和蒸发器,平均占比分别为23.62%和20.04%,压缩机的相对㶲损最小,平均占比19.34%。
Fig. 16 The effects of air inlet temperature of evaporator on the relative exergy loss of system components图16 蒸发器进风温度对系统各部件相对㶲损的影响 |
Table 7 Average relative exergy loss of system components at different air inlet temperature of evaporator表7 不同的蒸发器进风温度下系统各部件平均相对㶲损 |
| 系统部件 | 平均相对㶲损/% |
|---|---|
| 蒸发器 | 20.04 |
| 气冷器 | 23.62 |
| 压缩机 | 19.34 |
| 膨胀阀 | 37.00 |
图17展示了蒸发器进风温度对总㶲损和㶲效率的影响。由图可知,随着蒸发器进风温度的升高,㶲损和㶲效率逐渐降低。当蒸发器进风温度从20 ℃升高至28 ℃时,总㶲损从1.230 kW减小到1.209 kW,㶲效率从17.90%减小到14.34%。这是由于随着蒸发器进风温度的升高,进风温度与制冷剂温差增大,换热器换热能力增加,为维持系统设定的5 ℃过热度,膨胀阀开度增大,压比减小。压缩机和气冷器内制冷剂的平均温度降低,不可逆损失减小,相应的㶲损降低。由于压比的减小,压缩机的容积效率增加,系统的制冷剂质量流量增加,虽然膨胀阀内的制冷剂平均温度略有降低,但在制冷剂流量增加的影响下膨胀阀的总㶲损微弱增加。蒸发器内制冷剂与外侧空气温度均升高,换热温差增大,同时在制冷剂流量增加的影响下㶲损略微增加。在各部件的综合影响下,系统总㶲损随着蒸发器进风温度的升高而略微减小。由于蒸发器进风温度的升高,虽然压缩机耗功减小,但系统获得的冷量㶲降低,导致㶲效率随蒸发器出风温度升高而降低。
Fig. 17 The effects of air inlet temperature of evaporator on the exergy efficiency and total exergy loss of system图17 蒸发器进风温度对总㶲损和㶲效率的影响 |
4 总结
构建了车用跨临界CO2空调系统,通过仿真对系统的变工况性能进行了研究,在仿真数据基础上对系统进行了㶲分析,得到以下结论：
（1）随着压缩机转速的提升,系统各部件㶲损失均增大,压缩机和气体冷却器的相对㶲损增加,膨胀阀和蒸发器的相对㶲损减小,系统的总㶲损增大,㶲效率降低。膨胀阀平均㶲损占比最大,平均相对㶲损为31.10%,蒸发器平均相对㶲损最小,为18.49%。
（2）随着气体冷却器进风风速的增加,气冷器、压缩机、膨胀阀㶲损失和相对㶲损减小,蒸发器㶲损失和相对㶲损增大,系统的总㶲损减小,㶲效率升高,增加和降低趋势逐渐缓慢。膨胀阀平均㶲损占比最大,平均相对㶲损为38.39%,蒸发器平均相对㶲损最小,为17.42%。
（3）随着气体冷却器进风温度的升高,气冷器、压缩机、膨胀阀㶲损失增大,蒸发器㶲损失减小,膨胀阀和压缩机的相对㶲损增加,气体冷却器和蒸发器的相对㶲损减小,系统的总㶲损增大,㶲效率升高。膨胀阀平均㶲损占比最大,平均相对㶲损为39.79%,蒸发器平均相对㶲损最小,为18.42%。
（4）随着蒸发器进风量的增加,气冷器和压缩机㶲损失和相对㶲损减小,膨胀阀㶲损失和相对㶲损增大,蒸发器㶲损失和相对㶲损先减后增,系统总㶲损先减后增,㶲效率先增后减。膨胀阀平均㶲损占比最大,平均相对㶲损为36.37%,蒸发器平均相对㶲损最小,为19.67%。
（5）随着蒸发器进风温度的升高,气冷器和压缩机㶲损失和相对㶲损减小,膨胀阀和蒸发器㶲损失和相对㶲损增大,系统总㶲损减小,㶲效率降低。膨胀阀平均㶲损占比最大,平均相对㶲损为37.05%,压缩机平均相对㶲损最小,为19.34%。
综上,在车用跨临界CO2空调系统中,膨胀阀和气体冷却器内的㶲损占比最大,同时,气体冷却器进风风速和进风温度对㶲损影响较大。因此在优化车用跨临界CO2空调系统时需要重点关注节流过程和气体冷却器内的冷却过程,搭建跨临界CO2空调系统时需要重点关注气体冷却器进风风速和进风温度与气体冷却器大小的匹配。
