王 燁，王瑞君，魯紅鈺，宋榮飛，藺虎相

太陽能蓄熱水箱內(nèi)置隔板開孔方式對(duì)蓄熱性能的影響

王 燁1,2，王瑞君1，魯紅鈺1，宋榮飛1，藺虎相1

（1. 蘭州交通大學(xué)環(huán)境與市政工程學(xué)院，蘭州 730070；2. 蘭州交通大學(xué)鐵道車輛熱工教育部重點(diǎn)實(shí)驗(yàn)室，蘭州 730070）

為研究太陽能蓄熱水箱內(nèi)具有相同開孔面積、不同開孔數(shù)量的隔板對(duì)熱分層及蓄熱性能的影響，該文對(duì)安裝有不同開孔方式隔板的不同流動(dòng)參數(shù)條件下水箱內(nèi)換熱過程進(jìn)行數(shù)值分析。結(jié)果表明：熱水和冷水入口溫度分別為333和303 K、熱水和冷水入口流速分別為0.05和0.5 m/s情況下，開孔面積為0.007 9 m2的1個(gè)開孔隔板水箱熱分層效率最高，理查森數(shù)為2.341；與相同開孔面積的1個(gè)開孔隔板相比，開孔面積為0.031 4 m2的5個(gè)開孔隔板使得水箱內(nèi)高溫水的容積得到顯、著擴(kuò)展，開孔數(shù)量對(duì)水箱冷熱水出口溫差的影響很微弱，開孔面積越大，開孔數(shù)量對(duì)理查森數(shù)的影響越明顯。當(dāng)用戶端用水量在一定范圍內(nèi)變化時(shí)，理查森數(shù)隨用戶用水量的增大而減小，在冷水入口流速低于0.3 m/s范圍內(nèi)增大該流速破壞熱分層的可能性更大。熱水入口溫度分別為333 和343 K對(duì)應(yīng)的熱水出口溫度最大相對(duì)偏差為0.22%，但這2種工況對(duì)冷水出口溫度幾乎沒有影響。研究結(jié)果可為液體分離設(shè)備中的孔板結(jié)構(gòu)優(yōu)化設(shè)計(jì)提供一定的理論參考。

太陽能；溫度；流場(chǎng)；蓄熱水箱；隔板開孔方式；熱分層效率

0 引 言

在太陽能熱水系統(tǒng)中，熱能收集和使用時(shí)間的不同步性使得蓄熱水箱在整個(gè)系統(tǒng)中起著至關(guān)重要的作用。因此，如何改善蓄熱水箱的蓄熱性能成為了行業(yè)研究熱點(diǎn)之一。其中，在水箱內(nèi)設(shè)置隔板是一種行之有效的方法，其原理主要是通過隔板改變流場(chǎng)結(jié)構(gòu)來優(yōu)化水箱內(nèi)的溫度場(chǎng)結(jié)構(gòu)。研究?jī)?nèi)容主要包括隔板高度、數(shù)量、安裝角度[1-4]，隔板形狀[5-7]等。另外，還可以通過改變水箱的形狀[8-10]、流體流動(dòng)參數(shù)[11-13]、在水箱內(nèi)設(shè)置流動(dòng)管或?qū)Я鞑踇14-16]、在水箱內(nèi)部或外部加設(shè)保溫材料[17-18]以及對(duì)水箱的放置方式和進(jìn)出口結(jié)構(gòu)進(jìn)行改進(jìn)[19-22]等措施來提高整個(gè)水箱的蓄熱能力。但對(duì)于相同隔板開孔面積、不同開孔方式對(duì)太陽能蓄熱水箱蓄熱性能影響的研究還未見報(bào)道。本文通過計(jì)算和對(duì)比隔板相同過流斷面面積條件下太陽能蓄熱水箱內(nèi)換熱過程的多個(gè)工況，對(duì)比分析了不同流動(dòng)參數(shù)以及不同隔板開孔方式對(duì)太陽能蓄熱水箱熱分層及蓄熱性能的影響。目的在于根據(jù)當(dāng)?shù)靥柲軓?qiáng)度、用戶對(duì)熱水的需求特點(diǎn)以及回水參數(shù)等因素設(shè)計(jì)合適的隔板結(jié)構(gòu)，調(diào)整流體參數(shù)，最大限度地利用太陽能資源并降低運(yùn)行能耗。

1 物理模型和數(shù)學(xué)模型

1.1 物理模型

本文以常用且文獻(xiàn)報(bào)道蓄熱性能良好、具有半球形頂?shù)膱A柱體水箱為例[8]，其外形及流體流向如圖1a所示。水箱直徑為1.0 m，高為1.5 m，冷熱水進(jìn)出口管徑均為0.02 m。計(jì)算中，水箱內(nèi)置隔板的開孔布局如圖1b、1c所示。開孔尺寸與數(shù)量如表1所示，其中內(nèi)置隔板距底面高度為0.2 m；a、b、c表示在隔板幾何中心開1個(gè)孔，孔直徑分別為0.1、0.2和0.3 m，d，e，f表示隔板幾何中心1個(gè)開孔、直徑0.5 m的圓周上均布4個(gè)相同直徑的孔，這5個(gè)孔的面積之和分別與a、b、c中孔面積相等，即隔板開孔面積分別為：1=a=d=0.007 9 m2，2=b=e= 0.031 4 m2，3=c=f= 0.070 7m2。

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

對(duì)水箱內(nèi)流動(dòng)與傳熱現(xiàn)象的研究是基于浮力模型，而浮升力項(xiàng)中的密度隨溫度的變化符合Boussinesq假設(shè)。求解流動(dòng)與傳熱的湍流模型非穩(wěn)態(tài)控制方程如下：

連續(xù)性方程：

動(dòng)量方程：

能量方程：

注：T1為來自集熱器的熱水入口溫度，K；T2為來自用戶端進(jìn)入蓄熱水箱的回水溫度，K；T3為供熱熱水出口溫度，K；T4為回到集熱器的冷水出口溫度，K。尺寸單位為mm.

表1 隔板開孔尺寸

其中

湍流動(dòng)能方程：

湍流動(dòng)能耗散率方程：

式中為溫度，K；為動(dòng)力黏度，(N·s)/m2；為密度，kg/m3；為體積膨脹系數(shù)，m3/K；p為定壓比熱，J/(kg·K）；k為湍流普朗特?cái)?shù)，其值在1.0左右；m為湍流長(zhǎng)度標(biāo)尺，m；1、2為經(jīng)驗(yàn)系數(shù)，標(biāo)準(zhǔn)模型中一般取1.44及1.92。其余符號(hào)及含義見文獻(xiàn)[23]。流體為水，密度變化采用Boussinesq假設(shè)。

1.3 邊界條件和初始條件

根據(jù)文獻(xiàn)[9]，設(shè)計(jì)不同冷熱水進(jìn)出口溫度、速度的組合計(jì)算工況，如表2所示。水箱內(nèi)初始溫度設(shè)為313 K，水箱中的壓力為1個(gè)大氣壓，冷、熱水出口設(shè)為自由出流邊界。水箱內(nèi)壁面設(shè)為絕熱邊界條件，考慮到隔板厚度0.003 m與水箱內(nèi)空間高度1.5 m相比尺度相差很大，可以忽略隔板本身導(dǎo)熱對(duì)傳熱過程的影響，即不考慮隔板上下表面間的溫度梯度，只將其看作對(duì)流動(dòng)有阻礙效應(yīng)的局部障礙物，因此，計(jì)算中將隔板表面設(shè)為絕熱邊界條件。所有液固交界面均為速度無滑移邊界條件。

表2 流體參數(shù)

2 數(shù)值求解方法

2.1 數(shù)學(xué)模型驗(yàn)證

為了驗(yàn)證1.2節(jié)數(shù)學(xué)模型的可靠性，采用該模型在相同條件下對(duì)文獻(xiàn)[24]的流動(dòng)與傳熱過程進(jìn)行數(shù)值試驗(yàn)，并與文獻(xiàn)[24]的實(shí)測(cè)結(jié)果進(jìn)行比對(duì)，結(jié)果如圖2所示。圖2中的溫度及時(shí)間均按照文獻(xiàn)[24]中的方法進(jìn)行無量綱化處理。可以看出，測(cè)試開始階段模擬結(jié)果與試驗(yàn)結(jié)果具有很好的一致性，兩者最大相對(duì)偏差為0.58%，這是因?yàn)殚_始階段熱質(zhì)交換過程還不夠充分，溫度場(chǎng)較均勻，兩者偏差不明顯。隨著測(cè)試的進(jìn)行，實(shí)測(cè)值與模擬值之間出現(xiàn)了一定的偏差，這是因?yàn)閷?shí)際中存在邊界絕熱不徹底以及流場(chǎng)與溫度場(chǎng)間的耦合作用逐漸凸顯的情況。但模擬結(jié)果與試驗(yàn)結(jié)果間的最大相對(duì)偏差為1.28%，證明所采用數(shù)學(xué)模型的正確性和可靠性。因此，可以采用該數(shù)學(xué)模型進(jìn)行后續(xù)計(jì)算[25]。

圖2 本文模擬結(jié)果與文獻(xiàn)[24]試驗(yàn)結(jié)果比對(duì)

2.2 網(wǎng)格劃分及獨(dú)立性驗(yàn)證

采用四面體網(wǎng)格對(duì)計(jì)算區(qū)域進(jìn)行離散。以裝有隔板a的水箱為例，采用3套網(wǎng)格進(jìn)行計(jì)算，取時(shí)間步長(zhǎng)為0.2 s，計(jì)算時(shí)長(zhǎng)為1 800 s，得到不同網(wǎng)格數(shù)下沿軸正向的溫度分布，如圖3所示?？梢钥闯?，3套網(wǎng)格所得計(jì)算結(jié)果吻合很好，其最大相對(duì)偏差為0.76%。綜合考慮計(jì)算的經(jīng)濟(jì)性與計(jì)算精度，選取394 710為后續(xù)計(jì)算的網(wǎng)格數(shù)。采用同樣方法得出裝有隔板b、c、d、e和f的水箱網(wǎng)格數(shù)分別為394 710、369 686、413 963、416 915和416 945。

圖3 網(wǎng)格獨(dú)立性驗(yàn)證

2.3 時(shí)間步長(zhǎng)驗(yàn)證

分別取時(shí)間步長(zhǎng)為0.15、0.20和0.25 s，以裝設(shè)有隔板a的水箱為例，進(jìn)行迭代計(jì)算。得到水箱內(nèi)沿軸正向的溫度分布，如圖4所示。由圖可見，選取的3個(gè)時(shí)間步長(zhǎng)計(jì)算結(jié)果吻合很好，其相對(duì)偏差不超過0.1%。因此，為了節(jié)省計(jì)算時(shí)間，后續(xù)計(jì)算取時(shí)間步長(zhǎng)為0.25 s。

2.4 求解計(jì)算

使用Fluent求解器計(jì)算，采用非耦合隱式求解方法。離散格式與文獻(xiàn)[12]相同，采用SIMPLE算法求解速度和壓力耦合問題[23]。松弛因子設(shè)置及求解控制方程時(shí)的收斂條件與文獻(xiàn)[26]相同。

圖4 時(shí)間步長(zhǎng)驗(yàn)證

3 結(jié)果與分析

分別通過改變內(nèi)置隔板的開孔方式、用戶端回流至蓄熱水箱的流速以及來自太陽能集熱器流入蓄熱水箱的熱水溫度3個(gè)因素來分析太陽能蓄熱水箱的蓄熱能力和熱分層效果。目前已經(jīng)有很多方法可以用來分析對(duì)比熱分層效率[27-30]。本文采用蓄熱水箱的溫度分布和理查森數(shù)來評(píng)判其熱分層效率。

理查森數(shù)是分析熱儲(chǔ)存問題中控制流動(dòng)和分層的重要參數(shù)[2]，是同時(shí)考慮浮升力、慣性力、黏滯力對(duì)流體質(zhì)點(diǎn)運(yùn)動(dòng)狀態(tài)的綜合作用效果的無量綱數(shù)。本研究中的垂直式蓄熱水箱中流體熱交換主要以自然對(duì)流換熱為主，因此可以使用理查森數(shù)來描述蓄熱水箱的熱分層效率，其表達(dá)式為：

式中為重力加速度，m/s2；為水箱高度，m；D為水箱熱水與冷水出口溫差，K；為冷水入口速度，m/s。

3.1 內(nèi)置隔板開孔方式對(duì)溫度場(chǎng)的影響

以表2中的工況3為例，分析隔板的不同開孔方式對(duì)水箱內(nèi)溫度場(chǎng)的影響。結(jié)果如圖5所示，當(dāng)隔板只有1個(gè)孔時(shí)，隨著孔面積增大，水箱中高溫區(qū)域逐漸減小，裝有隔板b和c的水箱內(nèi)溫度場(chǎng)分布非常接近，這說明隔板開1個(gè)孔時(shí)，在流動(dòng)參數(shù)一定的情況下增大開孔面積會(huì)減少高溫水容積，再適當(dāng)增大開孔面積對(duì)高溫水容積的減少越來越不明顯，但隔板上方附近的低溫水影響區(qū)隨著開孔面積增大而增大。這是因?yàn)殚_孔面積與冷水經(jīng)圓孔沖入隔板上方的阻力有關(guān)，局部阻力的差異決定了水箱內(nèi)流場(chǎng)結(jié)構(gòu)的不同，從而以耦合的方式改變著對(duì)流換熱及導(dǎo)熱過程對(duì)整個(gè)溫度場(chǎng)的影響機(jī)制。由圖5d~5f可知，當(dāng)隔板有5個(gè)開孔時(shí)，高溫水區(qū)域隨開孔總面積增大而減小，但減小的幅度差異較大。隔板上方低溫水區(qū)域覆蓋范圍也隨著開孔總面積增大而增大，導(dǎo)致這一現(xiàn)象的主要原因與1個(gè)開孔的情況類似。與前者不同的是，5個(gè)孔時(shí)，在開孔面積較大時(shí)水箱內(nèi)高溫水容量差異明顯，而1個(gè)孔時(shí)的情況剛好相反。另外，在相同開孔面積情況下，5個(gè)孔對(duì)開孔面積1和3結(jié)構(gòu)水箱中高溫水區(qū)域影響微弱，但使得開孔面積為2的水箱內(nèi)高溫水容積得到了較大幅度的擴(kuò)展。尋找使得水箱內(nèi)高溫水容積最大、出水溫度最高的開孔方式是隔板結(jié)構(gòu)優(yōu)化的主要探索方向。

圖5 裝有不同開孔隔板的水箱豎向中心截面溫度分布

圖6為隔板不同開孔方式時(shí)水箱冷熱水出口溫差比較?？梢钥闯觯瑢?duì)于隔板1個(gè)開孔的情況，當(dāng)開孔面積由1=0.007 9 m2增大到3=0.070 7 m2，冷熱水出口溫差由12.74 K降低到11.12 K；對(duì)于隔板5個(gè)開孔的情況，當(dāng)開孔面積由1=0.007 9 m2增大到3=0.070 7 m2，冷熱水出口溫差由12.68 K降低到11.01 K。無論隔板1個(gè)開孔還是5個(gè)開孔，冷熱水出口溫差均隨著隔板上開孔面積的增大而減小。這是由于隔板開孔面積越小，隔板下方低溫水經(jīng)圓孔沖入隔板上方的阻力越大，水箱中心區(qū)穩(wěn)定的熱層起到了熱阻的作用，從而抑制了高低溫水之間的導(dǎo)熱換熱過程。但比較開孔面積相同的2種隔板結(jié)構(gòu)對(duì)應(yīng)的冷熱水出口溫差，發(fā)現(xiàn)開孔數(shù)目對(duì)于冷熱水出口溫差幾乎沒有影響。這是因?yàn)楫?dāng)流動(dòng)穩(wěn)定時(shí)換熱過程也達(dá)到了動(dòng)態(tài)的熱平衡狀態(tài)。

圖7為不同開孔方式下開孔面積對(duì)理查森數(shù)的影響對(duì)比?？梢钥闯?，無論隔板1個(gè)開孔還是5個(gè)開孔，理查森數(shù)均隨著隔板上開孔直徑的增大而減小。當(dāng)開孔面積1=0.007 9 m2時(shí)，1個(gè)開孔的理查森數(shù)為2.341，比5個(gè)開孔高出0.012；當(dāng)開孔面積S=0.031 4 m2時(shí)，1個(gè)開孔的理查森數(shù)比5個(gè)開孔高出0.013；當(dāng)開孔面積為S=0.070 7 m2時(shí)，1個(gè)開孔的理查森數(shù)比5個(gè)開孔高出0.022，即開孔面積越大，開孔數(shù)對(duì)理查森數(shù)的影響越明顯，對(duì)熱分層越不利。因此，表1中隔板開孔面積最小的結(jié)構(gòu)a型水箱內(nèi)的熱分層效率最高。

圖6 裝有不同開孔隔板的水箱冷熱水出口溫差

圖7 開孔方式對(duì)理查森數(shù)的影響

綜合圖5中、圖6中及圖7中可知，開孔方式a的隔板結(jié)構(gòu)性能最優(yōu)。

3.2 用戶用水量對(duì)溫度場(chǎng)的影響

在工程實(shí)際中，由于各種因素導(dǎo)致用戶用水量一般呈動(dòng)態(tài)變化，本研究中以冷水入口流速變化來體現(xiàn)這一過程。后續(xù)計(jì)算中分別選取5個(gè)冷水入口流速，作為工程實(shí)際中不同時(shí)段的流速平均值，即表2中的0.1，0.3，0.5，0.7和0.9 m/s。以上文所獲得的性能最優(yōu)的開孔結(jié)構(gòu)a型隔板水箱為例，研究不同用戶用水量情況下的蓄熱性能變化規(guī)律。熱水入口溫度1取333 K。

圖8為其他條件相同、不同冷水入口流速時(shí)水箱內(nèi)的溫度場(chǎng)計(jì)算結(jié)果。可以看出，冷水入口流速由0.1 m/s增大到0.9 m/s的過程中，隔板上方320 K溫度區(qū)域逐漸減小，溫度梯度越來越大，溫度整體水平降低，隔板開孔位置處水溫由318 K降低到312 K，隔板下方低溫區(qū)域不斷擴(kuò)展，溫度整體水平降低，最低水溫由316 K降低到304 K。

圖9為過底面圓心沿軸正向鉛垂線上的水溫變化與冷水入口流速間的關(guān)系。冷水入口流速由0.1 m/s增大到0.9 m/s的過程中，在<0.25 m范圍，水箱底部的最大溫差達(dá)到了11 K，隨冷水入口流速的增大水箱底部水溫的降低幅度逐漸減??；在0.25 m<<1.2 m范圍，各水箱溫差逐漸減小，冷水入口流速最小值對(duì)應(yīng)的溫度水平整體偏高。冷水入口流速對(duì)熱水出口附近溫度影響很微弱，溫度最大差值為0.78 K。

以上分析表明：當(dāng)用戶端用水量變大時(shí)，水箱內(nèi)的熱水量可能滿足不了用戶的需求，此時(shí)可考慮不同冷水入口流速的2個(gè)水箱聯(lián)合運(yùn)行，以提高供水系統(tǒng)的可靠性。

注：熱水入口溫度333 K

圖9 冷水入口流速對(duì)z軸正向溫度影響

圖10為理查森數(shù)隨冷水入口流速的變化關(guān)系?？梢钥闯觯聿樯瓟?shù)隨著冷水入口流速的增大而減小。冷水入口速度從0.1 m/s增大到0.9 m/s過程中，理查森數(shù)相應(yīng)地從24減小到0.9。其中，當(dāng)冷水入口流速大于0.3 m/s時(shí)，理查森數(shù)的下降趨勢(shì)較為平緩，在冷水入口流速從0.1 m/s上升到0.3 m/s 的過程中理查森數(shù)下降迅速。這說明在冷水入口流速較低時(shí)增大該流速破壞熱分層的可能性更大。

3.3 熱水入口溫度對(duì)溫度場(chǎng)的影響

在太陽能集熱系統(tǒng)中，多種因素的影響會(huì)使從集熱器流出的熱水溫度也就是蓄熱水箱熱水入口的水溫變化，為了分析熱水入口的溫度變化對(duì)水箱蓄熱性能的影響，本文在邊界條件設(shè)置時(shí)分別選取333和343 K 共2種不同的熱水入口溫度[9]。下面以表2中工況3和工況8為例進(jìn)行分析。

圖10 理查森數(shù)隨冷水入口流速變化

圖11為其他條件相同、熱水入口溫度不同時(shí)水箱內(nèi)縱斷面上的溫度分布。可以看出，熱水入口溫度從333 K提高到343 K時(shí)，隔板開孔位置處水溫均為316和320 K熱水容量幾乎相等，隔板下方最低水溫均為308 K。

圖12為其他條件相同、不同熱水入口溫度時(shí)過底面圓心沿軸正向鉛垂線上的溫度比較?？梢钥闯?，熱水入口溫度從333 K提高到343 K時(shí)，在<1.1 m范圍內(nèi)，2種情況的溫度相對(duì)平均偏差僅為0.04%，到了>1.1 m范圍才體現(xiàn)出熱水入口溫度對(duì)熱水出口溫度的影響，熱水出口溫度差最大值出現(xiàn)在水箱出口處，最大相對(duì)偏差為0.22%。但這2種工況對(duì)冷水出口溫度幾乎沒有影響。因此，對(duì)于有內(nèi)置隔板的蓄熱水箱，考慮太陽輻射強(qiáng)度變化的隨機(jī)性引起熱水入口水溫的變化進(jìn)而影響熱水出口溫度的變化但對(duì)冷水出口水溫?zé)o明顯影響這一特點(diǎn)，可以根據(jù)不同水溫的用戶需求來設(shè)計(jì)流向可調(diào)的熱水供應(yīng)系統(tǒng)，即對(duì)于熱水出口溫度能滿足用戶需求的情況直接供水；對(duì)于熱水出口溫度高于用戶需求的情況，可以考慮通過調(diào)節(jié)來自用戶回水管路上的閥門開度以提高水箱的冷水入口流速，從而降低熱水出口溫度；對(duì)于熱水出口溫度低于用戶需求溫度的情況，就要考慮與其他輔助熱源的聯(lián)合使用。

圖11 熱水入口溫度對(duì)水箱中心豎向截面溫度場(chǎng)影響

圖12 熱水入口溫度對(duì)z軸正向溫度影響

4 結(jié) 論

本文以內(nèi)置不同開孔形式隔板的蓄熱水箱為研究對(duì)象，采用標(biāo)準(zhǔn)模型對(duì)其內(nèi)部流體區(qū)域的換熱過程進(jìn)行了數(shù)值分析，得到如下主要結(jié)論：

1）對(duì)于隔板上1個(gè)開孔、流動(dòng)參數(shù)相同的情況下，隨著隔板開孔面積的增大，水箱中高溫區(qū)域范圍會(huì)變小，開孔面積增大到一定程度時(shí)高溫水區(qū)域縮減的幅度會(huì)降低，但熱分層效率隨著開孔面積的增大而減小，開孔面積為0.007 9 m2的1個(gè)開孔隔板水箱熱分層效率最高，理查森數(shù)為2.341。

2）當(dāng)隔板有5個(gè)開孔時(shí)，高溫水區(qū)域隨開孔總面積增大而減小，但減小的幅度差異較大。即開孔面積越大，開孔數(shù)對(duì)理查森數(shù)的影響越明顯，對(duì)熱分層越不利。開孔面積為0.070 7 m2時(shí)，1個(gè)開孔的理查森數(shù)與5個(gè)開孔的差值最大，為0.022；開孔面積為0.007 9 m2時(shí)，1個(gè)開孔的理查森數(shù)與5個(gè)開孔的差值最小，為0.012。

3）與相同開孔面積的1個(gè)開孔隔板相比，隔板上開5個(gè)孔對(duì)于最小開孔面積0.007 9 m2和最大開孔面積0.070 7 m2對(duì)應(yīng)的水箱內(nèi)高溫水容積影響很小，使得中等開孔面積0.031 4 m2水箱內(nèi)高溫水的容積得到了顯著擴(kuò)展，但5個(gè)開孔對(duì)水箱冷熱水出口溫差的影響很微弱。

4）當(dāng)用戶端的用水量動(dòng)態(tài)變化時(shí)，理查森數(shù)隨用戶的用水量增大而減小。在冷水入口流速低于0.3 m/s范圍內(nèi)增大該流速值破壞熱分層的可能性更大。

5）在實(shí)際工程中，考慮太陽輻射強(qiáng)度的隨機(jī)性，對(duì)于給定的流體參數(shù)情況下，可以考慮將太陽能蓄熱水箱與輔助熱源聯(lián)合運(yùn)行的方式來滿足不同用戶的需求。

[1] Bouhal T, Fertahi S, Agrouaz Y, et al. Numerical modeling and optimization of thermal strati?cation in solar hot water storage tanks for domestic applications: CFD study [J]. Solar Energy, 2017, 157: 441－455.

[2] Erdemir D, Altuntop N. Improved thermal stratification with obstacles placed inside the vertical mantled hot water tanks [J]. Applied Thermal Engineering, 2016, 100: 20－29.

[3] Mevlut A, Necdet A Veysel O, et al. Optimum dimensions of an obstacle placed in a hot water tank for thermal stratifi-cation [J]. Energy, Exploration & Exploitation, 2009, 23(5): 357－373.

[4] Tian H, Zheng S Z. Lorentz estimates for the gradient of weak solutions to elliptic obstacle problems with partially BMO coefficients [J]. Boundary Value Problems, 2017, 2017(1): 1－27.

[5] Altuntop N, Arslan M, Ozceyhan V, et al. Effect of obstacles on thermal strati?cation in hot water storage tanks[J]. Applied Thermal Engineering，2005, 25: 2285－2298.

[6] 張忠斌，黃虎，杜塏，等. 隔板裝置對(duì)恒溫水箱的性能影響分析(1)數(shù)值設(shè)計(jì)[J]，太陽能學(xué)報(bào)，2016，37(2)：469－475. Zhang Zhongbin, Huang Hu, Du Kai, et al. Analysis of the effect of diaphragm on the thermostatic water tank (1) Numerical design [J], Acta Energine Solars Sinica, 2016, 37(2): 469－475. (in Chinese with English abstract)

[7] Zachár A, Farkbs I, Szlivka F. Numerical analyses of the impact of plates for thermal stratification inside a storage tank with upper and lower inlet flows [J]. Solar Energy, 2003, 74: 287－302.

[8] Yang Z, Chen H S, Wang L, et al. Comparative study of the influences of different water tank shapes on thermal energy storage capacity and thermal stratification[J]. Renewable Energy, 2016, 85: 31－44.

[9] 王燁，鮑成柯，付銀安，等. 圓錐形頂太陽能蓄熱水箱錐頂結(jié)構(gòu)及運(yùn)行參數(shù)優(yōu)化[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2017，33(22): 255－261. Wang Ye, Bao Chengke, Fu Yin’an, et al. Structure and operating parameter optimization of solar hot water storage tank with conical top[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(22): 255－261. (in Chinese with English abstract)

[10] 胡延鐸，李亞奇，宋鴻杰，等. 殼管式與圓柱式梯級(jí)相變蓄熱裝置的數(shù)值模擬與比較[J]. 兵器裝備工程學(xué)報(bào)，2015，36(3)：102－106. Hu Yanduo, Li Yaqi, Song Hongjie, et al. Numerical simulation and comparison of shell-and-cylinder and cylindrical Step change thermal storage devices[J]. Journal of Ordnance Equipment Engineering, 2015, 36(3): 102－106. (in Chinese with English abstract)

[11] E Azad. Design installation and operation of a solar thermal public bath in eastern Iran[J]. Energy Sustain, 2012, 16(1): 68－73.

[12] 王登甲，劉艷峰. 太陽能熱水采暖蓄熱水箱溫度分層分析[J]. 建筑熱能通風(fēng)空調(diào)，2010，29(1)：16—19. Wang Dengjia, Liu Yanfeng. Temperature stratification studying of heat storage tank[J]. Building Energy & Environment, 2010, 29(1): 16－19. (in Chinese with English abstract)

[13] 李琳. 太陽能熱水器貯熱水箱內(nèi)混水特性的數(shù)值模擬分析[D]. 昆明：云南師范大學(xué)，2011.

Li Ling. Numerical Simulation and Analysis on Characteristic of Water Mixing to the Solar Heater’s Heat sStorage Water Tank[D]. Kun Ming: Yunnan Normal University, 2011. (in Chinese with English abstract)

[14] Zachár A. Investigation of a new tube-in-tube helical flow distributor design to improve temperature stratification inside hot water storage tanks operated with coiled-tube heat exchangers[J]. International Journal of Heat and Mass Transfer, 2013, 63: 150－161.

[15] 王曉瑩，王琪，于淼. 溫度分層式太陽能蓄熱水箱模擬分析與試驗(yàn)研究[J]. 中國(guó)農(nóng)機(jī)化學(xué)報(bào)，2014，35(5)：192－194. Wang Xiaoyin, Wang Qi, Yu Miao. Simulative analyses and experimental research of temperature stratification in solar hot water tank[J]. Journal of Chinese Agricultural Mechanization, 2014, 35(5): 192－194. (in Chinese with English abstract)

[16] 楊征，陳海生，王亮，等. 豎直圓柱形水箱保溫過程熱分層現(xiàn)象與機(jī)理研究[J]. 中國(guó)電機(jī)工程學(xué)報(bào)，2015，35(6)：1420－1428.Yang Zheng, Chen Haisheng, Wang Liang, et al. Study on behavior and mechanism of thermal stratification of vertical cylindrical heat storage tank in insulation process [J]. Proceedings of the CSEE, 2015, 35(6): 1420－1428. (in Chinese with English abstract)

[17] Gasque M, Gonzaz-Altozano P, Maurer D, et al. Study of the in?uence of inner lining material on thermal strati?cation in a hot water storage tank[J]. Applied Thermal Engineering, 2015, 75: 344－356.

[18] Shah L J, Furbo S. Entrance effects in solar storage tanks[J]. Solar Energy, 2003, 75: 337－348.

[19] Moncho-Esteve I J, Gasque M, González-Altozano P. Simple inlet devices and their influence on thermal stratification in a hot water storage tank [J]. Energy and Buildings, 2017, 150: 625－638..

[20] Burak K, Korhan O. Effect of rectangular hot water tank position and aspect ratio on thermal stratification enhancem-ent[J]. Renewable Energy, 2018, 116: 639－646.

[21] Mohammad R A, Hassan B T, Mohammad J M. Experimental study on destruction of thermal stratification tank in solar collector performance[J]. Journal of Energy Storage, 2018, 15: 124－132.

[22] Wang Zilong, Zhang Hua, Dou Binlin, et al. Experimental and numerical research of thermal stratification with a novel inlet in a dynamic hot water storage tank[J]. Renewable Energy, 2017, 111: 353－371.

[23] 陶文銓. 數(shù)值傳熱學(xué)(第二版)[M]. 西安：西安交通大學(xué)出版社. 2013.

[24] Aviv A, Blyakhman Y, Beeri O, et al. Experimental and numerical study of mixing in a hot-water storage tank[J]. Journal of Solar Energy Engineering, 2009, 131: 1－8.

[25] 陽祥，陶文銓. 高瑞利數(shù)下封閉腔內(nèi)自然對(duì)流的數(shù)值模擬[J]. 西安交通大學(xué)學(xué)報(bào)，2014，48(5)：27－31. Yang Xiang, Tao Wenquan. Numerical simulations for natural convection with high Rayleigh number in a tall rectangular cavity[J]. Journal of Xi’an Jiaotong University, 2014, 48(5): 27－31. (in Chinese with English abstract)

[26] 王燁，徐燃. 青藏線空調(diào)列車室內(nèi)氣流組織的沿線變化數(shù)值分析[J]. 重慶大學(xué)學(xué)報(bào)，2014，37(9)：59－67. Wang Ye, Xu Ran. Numerical study of airflow distribution in air-conditioned train in Qinghai-Tibet Railway[J]. Journal of Chongqing University, 2014, 37(9): 59－67. (in Chinese with English abstract)

[27] Haller M Y, Cruickshank C A, Streeicher W, et al. Methods to determine stratification efficiency of thermal energy storage processes-review and theoretical comparison[J]. Solar Energy, 2009, 83 (10): 1847－1860.

[28] Arefeh H, Armin H, Sture H. Low-temperature heat emission combined with seasonal thermal storage and heat pump[J]. Solar Energy, 2015, 119: 122－133.

[29] 朱寧，王占恒，王大水，等. 太陽能蓄熱水箱的溫度分層研究[J]. 太陽能，2013(18)：16－18.Zhu Ning, Wang Zhanheng, Wang Dashui, et al. Study of thermal stratification inside a solar storage tank[J]. Solar Energy, 2013(18): 16－18. (in Chinese with English abstract)

[30] Liu Zhan, Li Yanzhong, Zhou Guoqing. Study on thermal stratification in liquid hydrogen tank under different gravity levels[J]. International Journal of Hydrogen Energy, 2018, 43: 9369－9378.

Influence of arrangement of holes on barrier of hot water tank of solar energy on heat storage performance

Wang Ye1,2, Wang Ruijun1, Lu Hongyu1, Song Rongfei1, Lin Huxiang1

(1.,,, 730070,2.,,, 730070,)

In solar water heating system, the discordance of heat energy collection and use time makes the hot water storage tank play a vital role in the whole system. Therefore, how to improve the thermal storage performance of the hot water storage tank has become one of the hot spots in industry research. Among many measures, as an effective method, setting barriers in the water storage tank has been adopted by more and more people. Its principle is to optimize the structure of temperature field in water tank by setting barriers to change the flow field structure.In this paper, in order to study the influence of the barrier with the same hole area and different numbers of hole in the solar water storage tank on thermal stratification and thermal storage performance, the process of heat transfer in a water tank with different openings is analyzed numerically by standardmodel under different flow parameters. The results show that: when the inlet temperatures of hot water and cold water is 333 and 303 K and the inlet velocity of hot water and cold water is 0.05 and 0.5 m/s, respectively, the water storage tank with a hole area of 0.007 9 m2has the highest thermal stratification efficiency, 2.341. When the hole area is the same as 0.0314 m2, installing the barrier with five holes will significantly expand the volume of high-temperature water in the water tank compared with installing the barrier with one hole. The number of holes has a weak influence on the outlet temperature difference between the hot and cold water of the tank, but the larger the hole area, the more obvious the influence of the numbers of hole on the Richardson number. If the user's water consumption changes within a certain range, the Richardson number decreases with the increase of the user's water consumption, and the heat stratification is more likely broken when the inlet velocity of cold water is lower than 0.3 m/s. The maximum relative deviation of the outlet temperature of hot water corresponding to the inlet temperature of hot water of 333 and 343 K is 0.22%, but these two working conditions have little effect on the outlet temperature of cold water. In this study, by calculating and comparing the heat transfer process in the solar water storage tank under the condition of the same hole area on barrier, the effects of different flow parameters and different arrangement of holes on barrier on the thermal stratification and thermal storage performance of the solar water storage tank are compared and analyzed, which has certain theoretical reference for the optimization design of barrier structure with openings in liquid separation equipment in different researches of various industries. In addition, due to the randomness of solar radiation intensity in practical engineering, if the fluid parameters are given, the way of combining solar water storage tank and auxiliary heat source to meet the needs of different users should be considered.

solar energy; temperature; flow field; hot water storage tank; arrangement of holes on barrier; thermal stratification efficiency

10.11975/j.issn.1002-6819.2019.10.024

TU822

A

1002-6819(2019)-10-0191-08

2019-01-20

2019-02-19

國(guó)家自然科學(xué)基金（51476073, 51266004）；甘肅省建設(shè)科技攻關(guān)項(xiàng)目（KJ2016-2）

王 燁，教授, 博士生導(dǎo)師，從事自然對(duì)流及暖通空調(diào)方面研究。Email：wangye@mail.lzjtu.cn

王 燁，王瑞君，魯紅鈺，宋榮飛，藺虎相. 太陽能蓄熱水箱內(nèi)置隔板開孔方式對(duì)蓄熱性能的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2019，35(10)：191－198. doi：10.11975/j.issn.1002-6819.2019.10.024 http://www.tcsae.org

Wang Ye, Wang Ruijun, Lu Hongyu, Song Rongfei, Lin Huxiang. Influence of arrangement of holes on barrier of hot water tank of solar energy on heat storage performance [J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(10): 191－198. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2019.10.024 http://www.tcsae.org