劉睿盈,吳靜怡，孫曉琳

(上海交通大學(xué)制冷與低溫工程研究所，上海 200240)

直膨式太陽(yáng)能熱泵熱水器集熱/蒸發(fā)器流道結(jié)構(gòu)分析與實(shí)驗(yàn)

劉睿盈*,吳靜怡，孫曉琳

(上海交通大學(xué)制冷與低溫工程研究所，上海 200240)

隨著我國(guó)經(jīng)濟(jì)和社會(huì)的發(fā)展,建筑能耗正在飛速上升。國(guó)家正大力研究開(kāi)發(fā)綠色建筑，充分利用可再生能源，因此太陽(yáng)能得到了廣泛的重視。太陽(yáng)能熱泵熱水器將熱泵與太陽(yáng)能相結(jié)合，生產(chǎn)生活熱水以滿足人們的需求。太陽(yáng)能熱泵熱水器中的集熱/蒸發(fā)器是重要的核心部件，只有優(yōu)化其結(jié)構(gòu)，才能更高效的利用太陽(yáng)能。本文采用Fluent軟件對(duì)太陽(yáng)能熱泵集熱/蒸發(fā)器三種不同流道形式的結(jié)構(gòu)單元進(jìn)行了數(shù)值分析，獲得了流道表面溫度場(chǎng)分布，通過(guò)對(duì)流道單元換熱特性的分析，對(duì)集熱/蒸發(fā)器的優(yōu)化方法進(jìn)行理論分析與實(shí)際應(yīng)用。

集熱/蒸發(fā)器；不同流道單元；Fluent軟件；溫度分布

0 前言

為滿足社會(huì)可持續(xù)發(fā)展的需求，需要盡可能地減少建筑能耗，一方面要推廣節(jié)能技術(shù)，另一方面也要充分利用可再生資源。太陽(yáng)能熱泵熱水器是近年發(fā)展起來(lái)的一個(gè)重要的節(jié)能產(chǎn)品，它將熱泵技術(shù)與太陽(yáng)能相結(jié)合，加熱熱水滿足人們的生活需要。太陽(yáng)能熱泵熱水器主要分為直膨式和非直膨式兩種。直膨式太陽(yáng)能熱泵熱水器將太陽(yáng)能集熱器與熱泵蒸發(fā)器結(jié)合，相比于非直膨式不僅簡(jiǎn)化了結(jié)構(gòu)，而且彌補(bǔ)了太陽(yáng)能熱水器陰雨天無(wú)法使用的缺點(diǎn)。太陽(yáng)能集熱/蒸發(fā)器可以安裝于建筑外立面，實(shí)現(xiàn)與建筑的一體化結(jié)合。作為太陽(yáng)能熱泵熱水器的核心部件，集熱/蒸發(fā)器的性能對(duì)太陽(yáng)能熱泵熱水器系統(tǒng)的工作效率有很大影響，因此對(duì)于集熱/蒸發(fā)器的結(jié)構(gòu)研究和優(yōu)化十分重要。在保證太陽(yáng)能集熱/蒸發(fā)器吸收效率的同時(shí)，減少材料的使用和占地空間是對(duì)其結(jié)構(gòu)改進(jìn)的主要方向[1]。

集熱/蒸發(fā)器的形式主要是平板式集熱/蒸發(fā)器。文獻(xiàn)[2]中針對(duì)有效集熱面積為 l m2全鋁板式結(jié)構(gòu)的集熱/蒸發(fā)器進(jìn)行了研究，證明其良好的工作性能。徐國(guó)英等[3]對(duì)一復(fù)合熱源熱泵系統(tǒng)中的螺旋翅片蒸發(fā)盤(pán)管平板型集熱/蒸發(fā)器進(jìn)行熱性能分析,證明其能在四季全天候高效工作，劉艷玲等和孫振華等[4-5]也對(duì)采用不同結(jié)構(gòu)集熱/蒸發(fā)器的熱水器系統(tǒng)進(jìn)行過(guò)研究。

本文主要就是針對(duì)平板吹脹式集熱/蒸發(fā)器的流道結(jié)構(gòu)單元進(jìn)行研究。吹脹式蒸發(fā)器就是將一張印有回路的板與另一張沒(méi)有印刷回路的板壓合，然后進(jìn)行吹脹管路而形成的，這種蒸發(fā)器具有熱傳導(dǎo)效率高、結(jié)構(gòu)簡(jiǎn)單、安裝方便等特點(diǎn)。本文通過(guò)對(duì)平板吹脹式集熱/蒸發(fā)器不同流道結(jié)構(gòu)的流動(dòng)換熱特性的模擬分析，從理論上分析了不同的流道結(jié)構(gòu)單元對(duì)集熱/蒸發(fā)器換熱效果的影響，為太陽(yáng)能集熱/蒸發(fā)器的設(shè)計(jì)和太陽(yáng)能熱泵熱水器系統(tǒng)的優(yōu)化提供理論依據(jù)和指導(dǎo)。

1 太陽(yáng)能熱泵熱水器系統(tǒng)及工作原理

直膨式太陽(yáng)能熱泵熱水器的蒸發(fā)器同時(shí)作為太陽(yáng)能集熱器使用，制冷劑在太陽(yáng)能集熱/蒸發(fā)器中直接獲取太陽(yáng)能吸熱蒸發(fā)。系統(tǒng)主要構(gòu)成部件包括：太陽(yáng)能集熱/蒸發(fā)器、壓縮機(jī)、冷凝器、熱力膨脹閥和儲(chǔ)熱水箱等,如圖1所示。

圖1 直膨式熱泵熱水器流程圖

這種熱水器的工作原理是：制冷劑作為循環(huán)工質(zhì)，流經(jīng)壓縮機(jī)后形成高溫高壓的氣體,高溫的制冷劑氣體流入蓄熱水箱內(nèi)部的冷凝器中,與冷凝器外部的水進(jìn)行熱交換,此時(shí)水箱內(nèi)的水吸熱升溫,同時(shí)制冷劑由于放熱冷凝變?yōu)橐后w，液體制冷劑流經(jīng)熱力膨脹閥節(jié)流降溫降壓后流入太陽(yáng)能集熱/蒸發(fā)器中,吸收太陽(yáng)輻射熱蒸發(fā)并再次進(jìn)入壓縮機(jī)，從而完成一次循環(huán)。當(dāng)水箱內(nèi)的水溫達(dá)到設(shè)定值時(shí)，壓縮機(jī)通過(guò)恒溫器控制自動(dòng)停機(jī)。在太陽(yáng)輻射值較高時(shí)，太陽(yáng)能熱泵熱水器以太陽(yáng)能為主要熱源；當(dāng)輻射值較小時(shí)，則主要從空氣中獲取熱量。除太陽(yáng)能集熱/蒸發(fā)器外，直膨式太陽(yáng)能熱泵的其余部件與常規(guī)熱泵相同。與常規(guī)熱泵相比，直膨式太陽(yáng)能熱泵的優(yōu)勢(shì)在于能夠達(dá)到更高的蒸發(fā)溫度，從而獲得較高的COP[6]。

2 流道單元模擬模型

為了研究不同結(jié)構(gòu)的流道對(duì)集熱/蒸發(fā)器性能的影響，本文選取三種不同的流道結(jié)構(gòu)單元，進(jìn)行了CFD模擬。

2.1 研究對(duì)象

本文主要研究三種不同流道結(jié)構(gòu)的吹脹式集熱/蒸發(fā)器。其結(jié)構(gòu)的幾何尺寸如表 1所示，三種流道結(jié)構(gòu)單元如圖2所示。

表1 鋁板幾何參數(shù)

圖2 三種流道結(jié)構(gòu)單元

2.2 數(shù)學(xué)模型

2.2.1 模型假設(shè)

在實(shí)際生產(chǎn)中，所研究的結(jié)構(gòu)單元按一定規(guī)律排布，形成集熱/蒸發(fā)器制冷劑流道。然而數(shù)值分析中只研究尺寸較小的單一結(jié)構(gòu)單元，故而假設(shè)制冷劑在三種流道結(jié)構(gòu)單元中沒(méi)有溫度變化。數(shù)學(xué)模型建立中有如下假設(shè)[7]：

1）流動(dòng)為定常流動(dòng)。

2）流體為不可壓縮的牛頓流體。

3）忽略重力和由于密度差異引起的浮力。

2.2.2 流道單元的控制方程[8]

1）流體質(zhì)量守恒方程：

式中：

u，v，w——速度的三個(gè)分量。

2）流體動(dòng)量方程：

式中：

式中：

a——流體的熱擴(kuò)散率，m2/s。

4）湍流模型

在本文中采用k-ε雙方程模型，是標(biāo)準(zhǔn)k-ε兩方程湍流模型。標(biāo)準(zhǔn)k-ε模型通過(guò)求解湍流動(dòng)能k方程和湍流耗散率ε方程，得到k和ε的解，再用k和ε的值計(jì)算得出湍流粘度，進(jìn)而求出雷諾應(yīng)力的解。

以下為湍流動(dòng)能k方程和湍流耗散率ε方程。k方程：

ε方程：

以上兩式中，C1ε、C2ε、σk、σε四個(gè)經(jīng)驗(yàn)常數(shù)一般通過(guò)對(duì)某些特定的湍流過(guò)程的分析和測(cè)量得到，其標(biāo)準(zhǔn)取值為：

3 網(wǎng)格模型及數(shù)值計(jì)算

3.1 網(wǎng)格模型

為了減小計(jì)算量從而縮小模型規(guī)模，故采用非結(jié)構(gòu)化的四面體網(wǎng)格。對(duì)計(jì)算區(qū)域內(nèi)的邊均采取均勻網(wǎng)格劃分，如圖3所示；并且為了驗(yàn)證網(wǎng)格的數(shù)量對(duì)模擬結(jié)果沒(méi)有影響，采用不同網(wǎng)格數(shù)的劃分方法進(jìn)行了驗(yàn)證。

圖3 三種結(jié)構(gòu)單元網(wǎng)格劃分圖

3.2 求解器求解

將網(wǎng)格文件導(dǎo)入數(shù)值求解器中，在設(shè)置粘性模型前需要先計(jì)算雷諾數(shù)：

式中：

ρ——制冷劑密度，kg/m3；

v——制冷劑流速，m/s；

η——?jiǎng)恿φ扯?，Pa·s；

d——特征長(zhǎng)度，m；由于制冷劑流道為矩形流道，故d= 4A/L，其中A為流道面積，L為濕周。

由于Re＞2500，流動(dòng)為湍流，故選擇k-ε雙方程模型。

定義R134a在0℃下飽和液體物性參數(shù)：密度為1294.78 kg/m3，定壓比熱容為1341 J/(kg·K)，導(dǎo)熱系數(shù)0.092 W/(m·K)，粘度2.665×10-4Pa·s。

并且數(shù)值分析時(shí)邊界條件設(shè)定見(jiàn)表 2，選取制冷劑入口流速為初始條件進(jìn)行迭代，研究穩(wěn)定后流道單元的溫度場(chǎng)。

表2 邊界條件設(shè)置

3.3 數(shù)值計(jì)算結(jié)果

選取輻射面為研究基準(zhǔn)面，分別得到三種不同結(jié)構(gòu)流道的表面溫度分布圖，如圖4所示。部分?jǐn)?shù)值計(jì)算數(shù)據(jù)見(jiàn)表3。

圖4 三種不同結(jié)構(gòu)流道的表面溫度分布圖

表3 仿真數(shù)據(jù)

3.4 結(jié)果分析

根據(jù)模擬結(jié)果，當(dāng)制冷劑蒸發(fā)溫度相等均為273 K時(shí)，采用矩形及六邊形流道的蒸發(fā)器表面溫度更低，蒸發(fā)器單位面積得熱量更大。在蒸發(fā)器面積相等、冷凝溫度與環(huán)境工況相同的條件下，相較于蛇形流道的蒸發(fā)器，矩形及六邊形流道的蒸發(fā)器能為制冷劑提供更多的熱量，此時(shí)熱力膨脹閥將增大制冷劑流量（同時(shí)蒸發(fā)溫度提高）以達(dá)到設(shè)定的過(guò)熱度。因此與蛇形流道相比，采用矩形或六邊形流道能夠提高蒸發(fā)溫度及蒸發(fā)器得熱量，從而提高系統(tǒng)COP及加熱功率。

4 實(shí)驗(yàn)結(jié)果及分析

仿真結(jié)果顯示，采用結(jié)構(gòu)單元2和3時(shí)，能夠有效增加集熱效率，改善蒸發(fā)器換熱性能，從而提高系統(tǒng)COP。為驗(yàn)證這一結(jié)論，在相同環(huán)境條件下，對(duì)采用三種結(jié)構(gòu)流道單元的集熱板進(jìn)行系統(tǒng)實(shí)驗(yàn)研究，流道網(wǎng)絡(luò)見(jiàn)圖 5，并且得到相關(guān)運(yùn)行參數(shù)，從實(shí)際運(yùn)行效果上，分析集熱器結(jié)構(gòu)形式對(duì)系統(tǒng)性能的影響，為集熱蒸發(fā)器結(jié)構(gòu)優(yōu)化提供依據(jù)。

實(shí)驗(yàn)樣機(jī)采用直膨式熱泵系統(tǒng)，太陽(yáng)能集熱/蒸發(fā)器均為裸板式，采用兩塊板并聯(lián)的方式連接，系統(tǒng)部件主要參數(shù)見(jiàn)表 4。壓縮機(jī)啟動(dòng)前，首先打開(kāi)充水閥將蓄熱水箱充滿，然后啟動(dòng)電功率表以及計(jì)算機(jī)自動(dòng)數(shù)據(jù)采集系統(tǒng)，每隔 55秒采集一次熱泵系統(tǒng)的各運(yùn)行參數(shù)并存儲(chǔ)至EXCEL中，直至蓄熱水箱內(nèi)的熱水溫度達(dá)到55℃，此時(shí)溫控器會(huì)控制壓縮機(jī)停止運(yùn)轉(zhuǎn)。

圖5 流道網(wǎng)絡(luò)圖

表4 系統(tǒng)的主要部件參數(shù)

4.1 結(jié)構(gòu)單元1和2對(duì)比實(shí)驗(yàn)

為比較結(jié)構(gòu)單元1和2，分別采用面積相同，流道結(jié)構(gòu)分別為蛇形和矩形的兩塊集熱板進(jìn)行實(shí)驗(yàn)，系統(tǒng)其它部件均相同，其中壓縮機(jī)排量為12.77 cm3/rev。部分實(shí)驗(yàn)結(jié)果見(jiàn)表5，可以看到在相同的大氣溫度和太陽(yáng)輻射值下，采用結(jié)構(gòu)單元2的系統(tǒng)累計(jì)COP為4.7，集熱因子為1.36，優(yōu)化集熱板結(jié)構(gòu)后系統(tǒng)COP明顯提高。

表5 實(shí)驗(yàn)結(jié)果對(duì)比

為了對(duì)比結(jié)構(gòu)優(yōu)化后的蒸發(fā)器的換熱性能，將采集到的數(shù)據(jù)進(jìn)行整理，比較結(jié)構(gòu)單元1和結(jié)構(gòu)單元2的集熱/蒸發(fā)器的蒸發(fā)溫度和系統(tǒng)累計(jì)COP隨水溫的變化，見(jiàn)圖6。

由圖6可見(jiàn)，當(dāng)水溫上升到55℃時(shí)，采用結(jié)構(gòu)單元2的蒸發(fā)溫度高于結(jié)構(gòu)單元1約為7℃，并且采用結(jié)構(gòu)單元2的系統(tǒng)COP提高15%。采用結(jié)構(gòu)單元2后改善了集熱板換熱性能，從而提高了系統(tǒng)蒸發(fā)溫度及COP。

圖6 結(jié)構(gòu)單元1和2蒸發(fā)器蒸發(fā)溫度及系統(tǒng)COP隨水溫變化

4.2 結(jié)構(gòu)單元2和3對(duì)比實(shí)驗(yàn)

為比較結(jié)構(gòu)單元2和3，分別采用面積相同、流道結(jié)構(gòu)分別為矩形和六邊形的兩塊集熱板進(jìn)行實(shí)驗(yàn)，系統(tǒng)其它部件均相同，其中壓縮機(jī)排量為排8.77 cm3/rev。部分實(shí)驗(yàn)結(jié)果見(jiàn)表6，可以看到在相同的大氣溫度和太陽(yáng)輻射值下，采用結(jié)構(gòu)單元3的系統(tǒng)累計(jì)COP為5.56，集熱因子為1.37，優(yōu)化集熱板結(jié)構(gòu)后系統(tǒng)性能略有提高。

表6 實(shí)驗(yàn)結(jié)果對(duì)比

為了對(duì)比結(jié)構(gòu)優(yōu)化后，蒸發(fā)器的換熱性能，將采集到的數(shù)據(jù)進(jìn)行整理，比較結(jié)構(gòu)單元2和3的集熱/蒸發(fā)器蒸發(fā)溫度及系統(tǒng)累計(jì)COP隨水溫的變化見(jiàn)圖7。

圖7 結(jié)構(gòu)單元2和3蒸發(fā)器蒸發(fā)溫度及系統(tǒng)COP隨水溫變化

由圖7可見(jiàn)，當(dāng)水溫上升到55℃時(shí)，采用結(jié)構(gòu)單元3的蒸發(fā)溫度略高于結(jié)構(gòu)單元2，并且采用結(jié)構(gòu)單元3的系統(tǒng)COP略有提升。采用結(jié)構(gòu)單元3后改善了集熱板換熱性能，從而提高了系統(tǒng)蒸發(fā)溫度及COP。

5 結(jié)論

本文以三種結(jié)構(gòu)不同的流道單元為研究對(duì)象，對(duì)太陽(yáng)能集熱/蒸發(fā)器三維流場(chǎng)進(jìn)行了流動(dòng)與傳熱數(shù)值模擬。得到了集熱/蒸發(fā)器表面溫度場(chǎng)分布圖，并對(duì)結(jié)果進(jìn)行了討論。仿真結(jié)果顯示，采用結(jié)構(gòu)單元2和3流道形式的集熱/蒸發(fā)器可以有效改善集熱板的傳熱性能，并且采用這兩種流道形式的集熱板溫度分布更加均勻，集熱板表面平均溫度更低，在蒸發(fā)溫度相同時(shí)，結(jié)構(gòu)單元2和3換熱量更大，大大提高太陽(yáng)能集熱/蒸發(fā)器的集熱效率。為了進(jìn)一步驗(yàn)證仿真結(jié)果，采用三種流道結(jié)構(gòu)的集熱/蒸發(fā)器帶入到系統(tǒng)實(shí)驗(yàn)當(dāng)中，并且對(duì)獲得的實(shí)驗(yàn)參數(shù)進(jìn)行處理和計(jì)算，從實(shí)驗(yàn)結(jié)果可以看出，采用結(jié)構(gòu)單元 2和 3的集熱/蒸發(fā)器在水箱內(nèi)初始水溫與終止水溫均相同情況下，系統(tǒng)COP和集熱因子有顯著提高。

[1]梁振南, 秦紅, 卜其輝. 太陽(yáng)能-空氣能集熱蒸發(fā)器現(xiàn)狀與發(fā)展方向分析[J]. 太陽(yáng)能, 2009 (12): 22-25.

[2]曠玉輝, 王如竹. 太陽(yáng)能熱泵熱水器的研制和性能測(cè)試[C]// 21世紀(jì)太陽(yáng)能新技術(shù)——2003 年中國(guó)太陽(yáng)能學(xué)會(huì)學(xué)術(shù)年會(huì)論文集. 2003: 567-571.

[3]徐國(guó)英, 張小松. 復(fù)合熱源熱泵系統(tǒng)集熱/蒸發(fā)器的模型[J]. 工程熱物理學(xué)報(bào), 2006, 27(1): 61-64.

[4]劉艷玲, 王如竹, 夏再忠. 一種新型的太陽(yáng)能吸附式制冷系統(tǒng)的設(shè)計(jì)及性能模擬[C]// 上海市制冷學(xué)會(huì)二○○五年學(xué)術(shù)年會(huì)論文集. 2005: 215-218.

[5]孫振華, 王如竹, 翟曉強(qiáng). 直膨式太陽(yáng)能熱泵熱水器實(shí)驗(yàn)研究及經(jīng)濟(jì)性分析[J]. 上海電力, 2008, 20(5): 475-479.

[6]矯洪濤, 王學(xué)生, 魏國(guó). 太陽(yáng)能熱泵技術(shù)研究進(jìn)展[J].能源技術(shù), 2007, 28(5): 270-274.

[7]劉楊. 運(yùn)用 fluent軟件對(duì)板式換熱器性能的數(shù)值分析[J]. 廣州化工, 2013, 41(13): 221-223.

[8]邵光杰, 王銳, 董紅星, 等. 材料物理化學(xué)[M]. 哈爾濱: 哈爾濱工業(yè)大學(xué)出版社, 2003.

Analysis of Collector/Evaporator for Direct Expansion Solar Assisted Heat Pump Water Heater

LIU Rui-ying*, WU Jing-yi, SUN Xiao-lin
(Institute of Refrigeration and Cryogenics, Shanghai Jiao Tong University, Shanghai 200240, China)

With the development of economic and social in our country, building energy consumption is rising. Government is vigorously developing green building to make full use of renewable energy. As a result, solar energy has been paid widely attention. Solar assisted heat pump water heater combines heat pump and solar energy to satisfy the domestic hot water demand. A solar assisted heat pump system was introduced in this paper. In a direct solar assisted heat pump system, the solar collector/evaporator is a critical component, and the performance of the collector/evaporator is affected a lot by the flow channel structure. Through the theoretical 3D numerical simulation of three different channel styles with Fluent software, the temperature distribution on the surface of the collector was obtained. Based on the heat transfer in the channels, a methodology for the design optimization of the collector/evaporator was introduced and applied.

Collector/evaporator; Different channels; Fluent software; Temperature distribution

10.3969/j.issn.2095-4468.2014.02.101

*劉睿盈（1988-），女，工程碩士。研究方向：太陽(yáng)能熱泵。聯(lián)系地址：上海市東川路800號(hào)上海交通大學(xué)制冷與低溫工程研究所，郵編：200240。聯(lián)系電話：13671779684。E-mail：ruiyingliu@163.com。