太陽(yáng)能煙囪發(fā)電系統(tǒng)內(nèi)傳熱問(wèn)題的數(shù)值分析

郭朋華,王元,李景銀

(西安交通大學(xué)流體機(jī)械及工程系, 710049, 西安)

太陽(yáng)能煙囪發(fā)電系統(tǒng)內(nèi)傳熱問(wèn)題的數(shù)值分析

郭朋華,王元,李景銀

(西安交通大學(xué)流體機(jī)械及工程系, 710049, 西安)

為了揭示輻射換熱、太陽(yáng)輻照度以及渦輪壓降對(duì)太陽(yáng)能煙囪發(fā)電系統(tǒng)性能的影響規(guī)律,基于系統(tǒng)內(nèi)部能量傳遞機(jī)制的理論分析,通過(guò)合理選擇邊界條件,在充分考慮各種氣象因素影響下的系統(tǒng)內(nèi)部多種熱效應(yīng)耦合的數(shù)值計(jì)算中引入了離散坐標(biāo)輻射模型、太陽(yáng)載荷模型及渦輪模型,以解決目前數(shù)值研究中物理模型及邊界條件選擇時(shí)需要明確的問(wèn)題。以西班牙示范電站為原型,對(duì)太陽(yáng)能煙囪發(fā)電系統(tǒng)內(nèi)部的流動(dòng)換熱行為進(jìn)行了詳細(xì)的三維數(shù)值模擬,并將數(shù)值模擬數(shù)據(jù)與西班牙示范電站的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了對(duì)比,二者具有良好的一致性。模擬結(jié)果表明:集熱棚內(nèi)的輻射換熱對(duì)棚內(nèi)溫度分布有重要影響,在數(shù)值模擬中不能忽略,否則會(huì)導(dǎo)致系統(tǒng)內(nèi)部換熱行為完全背離真實(shí)的換熱過(guò)程;集熱棚覆蓋層與外界環(huán)境對(duì)流、輻射換熱造成的熱量損失是系統(tǒng)內(nèi)的主要熱損失;渦輪壓降對(duì)系統(tǒng)性能也會(huì)產(chǎn)生顯著影響。

太陽(yáng)能煙囪;數(shù)值分析;輻射換熱;集熱棚效率;渦輪壓降

自1982年在西班牙Manzanares建成世界上第一座太陽(yáng)能煙囪示范電站以來(lái)[1],太陽(yáng)能煙囪發(fā)電技術(shù)已成為太陽(yáng)能大規(guī)模利用的研究熱點(diǎn)之一,它以空氣為工質(zhì),無(wú)需水資源,在太陽(yáng)能資源豐富、水資源匱乏的干旱、半干旱地區(qū)具有極其廣闊的應(yīng)用前景[2-4]。

太陽(yáng)能煙囪發(fā)電系統(tǒng)內(nèi)傳熱主要通過(guò)建立集熱棚一維熱平衡方程進(jìn)行熱力分析和計(jì)算[5-9],這種方法操作簡(jiǎn)便,但需要對(duì)系統(tǒng)進(jìn)行大量的假設(shè),也無(wú)法得到詳細(xì)的流場(chǎng)和溫度場(chǎng)信息。隨著計(jì)算機(jī)性能及數(shù)值計(jì)算方法的發(fā)展,數(shù)值模擬廣泛應(yīng)用于太陽(yáng)能煙囪發(fā)電系統(tǒng)的研究[10-12],目前的研究主要集中在二維數(shù)值模擬。由于太陽(yáng)能煙囪發(fā)電系統(tǒng)內(nèi)流動(dòng)和傳熱模型較為復(fù)雜,尤其是涉及到多種地理、氣象因素(太陽(yáng)輻射、風(fēng)速、環(huán)境溫度)影響下多效應(yīng)的耦合(導(dǎo)熱、對(duì)流、輻射)時(shí),各種物理模型的選擇及邊界條件確定中有許多問(wèn)題需要明確。

1 系統(tǒng)內(nèi)能量的傳遞過(guò)程

如圖1所示,太陽(yáng)能煙囪發(fā)電系統(tǒng)主要由集熱棚、風(fēng)力渦輪、煙囪和蓄熱層組成。集熱棚是一個(gè)利用溫室效應(yīng)制成的巨大的空氣加熱器,主體受力結(jié)構(gòu)一般為鋼支架,在其上覆蓋著玻璃或者塑料薄膜等半透明材料,這種材料對(duì)太陽(yáng)光的短波輻射有較高的穿透比,可以阻隔地面的長(zhǎng)波輻射,從而使棚內(nèi)空氣加熱、棚內(nèi)空氣密度下降得以實(shí)現(xiàn)。集熱棚中央是一個(gè)高達(dá)數(shù)百甚至上千米的煙囪,利用煙囪內(nèi)、外空氣的密度差形成強(qiáng)大的抽力,從而驅(qū)動(dòng)安裝在煙囪底部的風(fēng)力渦輪運(yùn)轉(zhuǎn)發(fā)電。在晚上及陰雨天,地面蓄熱層釋放熱量,繼續(xù)為集熱棚內(nèi)空氣加熱,保證電站持續(xù)運(yùn)行。

圖1 太陽(yáng)能煙囪發(fā)電系統(tǒng)模型

圖2為集熱棚內(nèi)能量的傳遞過(guò)程示意圖。集熱棚上表面吸收、反射和透射太陽(yáng)輻射,大部分太陽(yáng)輻射穿透集熱棚后被地面吸收。地面吸收的能量S2通過(guò)對(duì)流換熱傳遞給棚內(nèi)空氣,通過(guò)輻射換熱傳遞給集熱棚覆蓋層,通過(guò)導(dǎo)熱向地下傳遞,因此地面蓄熱層的熱平衡方程為

(1)

圖2 集熱棚內(nèi)能量的傳遞過(guò)程示意圖

集熱棚材料對(duì)太陽(yáng)短波輻射的穿透比很大,但是對(duì)地面的長(zhǎng)波輻射穿透比很小,所以集熱棚覆蓋層除了吸收小部分太陽(yáng)輻射外,其主要能量來(lái)自與地面間的輻射換熱。集熱棚覆蓋層吸收的能量一部分通過(guò)對(duì)流換熱傳遞給棚內(nèi)空氣,另一部分通過(guò)與環(huán)境空氣的對(duì)流換熱及對(duì)天空的長(zhǎng)波輻射損失掉,因此集熱棚覆蓋層的熱平衡方程為

S1+hr,pc(Tp-Tc)=

hc,f(Tc-Tf)+hr,cs(Tc-Ts)+hc,a(Tc-Ta)

(2)

集熱棚覆蓋層和地面蓄熱層與棚內(nèi)空氣的對(duì)流換熱量即為棚內(nèi)空氣的吸熱量,其熱平衡方程為

(3)

定義Qcoll為集熱棚內(nèi)空氣吸收的熱量,Qsolar為系統(tǒng)內(nèi)總的太陽(yáng)入射能量,Pt為渦輪的輸出功率,則集熱棚效率和系統(tǒng)效率可分別定義為

(4)

(5)

2 數(shù)值模型

2.1 物理模型

西班牙示范電站是迄今為止建造成功的規(guī)模最大的太陽(yáng)能熱氣流發(fā)電系統(tǒng),其實(shí)驗(yàn)數(shù)據(jù)具有一定的代表性和權(quán)威性,是進(jìn)行計(jì)算驗(yàn)證的理想對(duì)象,該電站的主要參數(shù)如表1所示。

為減少阻力損失,集熱棚出口與煙囪的連接處采用光滑過(guò)渡,風(fēng)力渦輪在過(guò)渡段出口、距離地面9 m處。根據(jù)西班牙示范電站的實(shí)驗(yàn)數(shù)據(jù)[13],地下0.5m處的土壤溫度不存在實(shí)際性的波動(dòng),Pretorius等研究也得到類似的結(jié)論[14],因此本文蓄熱層厚度為2 m。模型的網(wǎng)格劃分在Fluent的前處理軟件Gambit中完成,為了便于劃分,將整個(gè)系統(tǒng)分為蓄熱層、集熱棚、煙囪及煙囪和集熱棚之間的過(guò)渡轉(zhuǎn)換段,且均采用結(jié)構(gòu)化網(wǎng)格。轉(zhuǎn)換段氣流從水平方向變成豎直方向時(shí),通流面積變化較大,該段氣流的速度、壓力及溫度的梯度較大。渦輪安裝在轉(zhuǎn)換段,該區(qū)域?yàn)檎麄€(gè)系統(tǒng)的核心,因此轉(zhuǎn)換段的網(wǎng)格需進(jìn)行加密處理,整個(gè)計(jì)算區(qū)域的網(wǎng)格劃分見圖3,網(wǎng)格總數(shù)達(dá)150萬(wàn),滿足網(wǎng)格無(wú)關(guān)性要求。

表1 西班牙太陽(yáng)能煙囪示范電站主要的參數(shù)

圖3 計(jì)算區(qū)域網(wǎng)格劃分

2.2 計(jì)算模型和邊界條件

由于集熱棚頂通常使用半透明材料,而在Fluent提供的幾種輻射模型中只有離散坐標(biāo)輻射模型(DO)可以處理半透明壁面,可以提供灰?guī)Ｐ蛠?lái)計(jì)算非灰輻射,因此本文選擇DO輻射模型來(lái)模擬集熱棚內(nèi)的輻射換熱。

對(duì)于太陽(yáng)輻射,以往的文獻(xiàn)中通常將透射的太陽(yáng)輻射以給定地面熱流密度或者土壤內(nèi)熱源的方式來(lái)處理[10,12,15]。本文采用Fluent提供的太陽(yáng)能加載模型來(lái)模擬太陽(yáng)輻射。該模型可根據(jù)太陽(yáng)位置向量和照射參數(shù)將太陽(yáng)輻射處理成射線,并應(yīng)用于用戶指定的邊界上,然后通過(guò)進(jìn)行面與面之間的遮蔽分析,來(lái)確定計(jì)算區(qū)域內(nèi)被遮擋的部分,從而計(jì)算出入射太陽(yáng)輻射在邊界上產(chǎn)生的熱流,該熱流將加載到能量方程的源項(xiàng)。由于渦輪內(nèi)的流場(chǎng)不是本文的研究對(duì)象,所以采用Fluent提供的壓降模型來(lái)模擬渦輪,該模型將渦輪視為一個(gè)無(wú)限薄的圓盤,通過(guò)渦輪的壓降可指定為常數(shù)或者是通過(guò)渦輪氣流速度的函數(shù)。

集熱棚進(jìn)口和煙囪出口采用壓力邊界條件,為了消除進(jìn)、出口壓差在系統(tǒng)內(nèi)產(chǎn)生的自然對(duì)流,進(jìn)、出口處的壓力均設(shè)置為0 Pa。在蓄熱層的底面,溫度不會(huì)隨環(huán)境溫度的日變化出現(xiàn)波動(dòng),此處設(shè)置為恒溫邊界條件;集熱棚頂為對(duì)流和輻射換熱邊界條件,對(duì)天空的長(zhǎng)波輻射可通過(guò)天空當(dāng)量溫度[16]獲得,即

(6)

太陽(yáng)能煙囪發(fā)電系統(tǒng)的邊界條件設(shè)置如表2所示。

表2 太陽(yáng)能煙囪發(fā)電系統(tǒng)的邊界條件設(shè)置

3 結(jié)果與討論

3.1 數(shù)值模型驗(yàn)證

由數(shù)值模擬得到的系統(tǒng)輸出功率與西班牙示范電站的實(shí)驗(yàn)數(shù)據(jù)[17]對(duì)比如圖4所示。由圖4可以看出,數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果吻合較好。數(shù)值模擬得到的輸出功率在總體上略高于實(shí)驗(yàn)結(jié)果,主要由兩方面原因:其一,本文采用的是穩(wěn)態(tài)數(shù)值模擬,在給定的太陽(yáng)輻射條件下,集熱棚覆蓋層、棚內(nèi)空氣及地面蓄熱層均已達(dá)到了該工況下的熱平衡,而實(shí)際土壤存在熱慣性,因此數(shù)值模擬的預(yù)測(cè)結(jié)果偏高;其二,實(shí)驗(yàn)電站的集熱棚內(nèi)存在大量的鋼架支撐結(jié)構(gòu),其會(huì)產(chǎn)生一定的阻力損失,為了便于建模和計(jì)算,數(shù)值模擬中忽略了這些鋼架結(jié)構(gòu)的影響。綜上分析,本文提出的數(shù)值模型完全適用于太陽(yáng)能熱氣流發(fā)電系統(tǒng)的模擬計(jì)算。

圖4 系統(tǒng)輸出功率的數(shù)值模擬數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)對(duì)比

3.2 輻射換熱對(duì)集熱棚內(nèi)溫度分布的影響

在以往的太陽(yáng)能煙囪發(fā)電系統(tǒng)的數(shù)值模擬研究中,集熱棚被視為普通的日光溫室,即忽略了棚內(nèi)的輻射換熱。實(shí)際的集熱棚工作過(guò)程并不等同于一般的日光溫室,集熱棚中部的高大煙囪會(huì)使得集熱棚不再是一個(gè)密閉空間。在煙囪抽力的作用下,空氣不斷地從集熱棚四周進(jìn)入棚內(nèi),加熱后的空氣最終通過(guò)煙囪排放到大氣中。因此,在太陽(yáng)能煙囪發(fā)電系統(tǒng)中,集熱棚覆蓋層對(duì)地面長(zhǎng)波輻射的阻隔作用成為集熱棚實(shí)現(xiàn)加熱功能的重要組成部分。根據(jù)傳熱學(xué)的基本理論,對(duì)于表面溫度為幾十?dāng)z氏度的表面散熱,對(duì)流換熱量與輻射換熱量具有相同的數(shù)量級(jí),必須同時(shí)予以考慮[18]。

(a)忽略棚內(nèi)輻射換熱

(b)考慮棚內(nèi)輻射換熱

為了研究輻射換熱對(duì)模擬結(jié)果的影響,本文將考慮輻射換熱前后的模擬結(jié)果進(jìn)行了對(duì)比。以太陽(yáng)輻射強(qiáng)度為600W/m2為例,考慮輻射換熱前后的集熱棚內(nèi)溫度分布如圖5所示。由圖5可以看出,未考慮輻射換熱時(shí),集熱棚內(nèi)的氣流溫度明顯高于考慮輻射換熱后的結(jié)果,尤其是地面溫度,未考慮輻射換熱時(shí)最高溫度達(dá)到了400K,在實(shí)際中顯然是不可能的。Pastohr等的數(shù)值研究認(rèn)為,這種不自然的高溫是由定常計(jì)算引起的[10],而從圖5的對(duì)比中可以看到這是數(shù)值模擬中地面輻射換熱被忽略的結(jié)果。

太陽(yáng)輻射強(qiáng)度為600W/m2時(shí),考慮輻射換熱前后集熱棚頂?shù)臏囟确植既鐖D6所示。由圖6可以看出,考慮輻射換熱后,集熱棚頂?shù)臏囟让黠@高于未考慮輻射換熱的情況。通過(guò)與圖5對(duì)比可知,當(dāng)忽略棚內(nèi)輻射換熱時(shí),集熱棚頂溫度普遍低于棚內(nèi)的空氣溫度。這是因?yàn)榧療崤镯數(shù)臒崃恐饕獊?lái)自棚內(nèi)空氣的對(duì)流換熱和一小部分的太陽(yáng)輻射?？紤]輻射換熱后,受地面輻射的影響,集熱棚頂溫度顯然會(huì)高于棚內(nèi)空氣溫度,此時(shí)集熱棚頂不僅不會(huì)從棚內(nèi)空氣吸收熱量,反而會(huì)將一部分熱量傳遞給空氣。

(a)忽略棚內(nèi)輻射換熱

(b)考慮棚內(nèi)輻射換熱

輻射換熱對(duì)集熱棚集熱效率的影響如圖7所示。由圖7可以看出,當(dāng)考慮棚內(nèi)輻射換熱時(shí),在不同輻照強(qiáng)度下得到的集熱棚集熱效率均遠(yuǎn)遠(yuǎn)高于考慮輻射換熱時(shí)的情況。西班牙示范電站的實(shí)驗(yàn)數(shù)據(jù)中,集熱棚的集熱效率均在30%左右,顯然考慮輻射換熱后的數(shù)值預(yù)測(cè)結(jié)果與實(shí)驗(yàn)結(jié)果更加吻合。這是因?yàn)楹雎耘飪?nèi)輻射換熱背離了實(shí)際傳熱過(guò)程而導(dǎo)致不正確的溫度分布所致。由圖5和圖6還可以看出,當(dāng)未考慮地面與集熱棚覆蓋層之間的輻射換熱時(shí),地面溫度過(guò)高,集熱棚頂溫度過(guò)低。在太陽(yáng)能煙囪發(fā)電系統(tǒng)中,由于通過(guò)集熱棚頂?shù)臒崃繐p失遠(yuǎn)遠(yuǎn)大于通過(guò)地面的熱量損失而忽略輻射換熱,使得集熱棚頂?shù)臒崃繐p失會(huì)被大大低估,從而導(dǎo)致集熱棚效率的預(yù)測(cè)值偏高。

圖7 輻射換熱對(duì)集熱棚的集熱效率的影響

3.3 渦輪壓降對(duì)系統(tǒng)性能的影響

利用Fluent軟件的壓降模型可以為渦輪指定一系列的壓降,以此研究壓降對(duì)系統(tǒng)性能的影響。本文模擬中給定的渦輪壓降間隔為20Pa,太陽(yáng)輻射強(qiáng)度分別設(shè)定為200、400、600、800W/m2。

(a)對(duì)速度的影響

(b)對(duì)溫升的影響

不同太陽(yáng)輻射下渦輪壓降對(duì)煙囪入口處氣流速度和溫升的影響如圖8所示。由圖8可以看出:當(dāng)太陽(yáng)輻射強(qiáng)度不變時(shí),隨著渦輪壓降的增加,渦輪對(duì)氣流的阻塞效應(yīng)增強(qiáng),系統(tǒng)需要更大的抽力來(lái)推動(dòng)渦輪做功,因此氣流速度逐漸降低;隨著氣流速度的降低,氣流在集熱棚內(nèi)停留的時(shí)間延長(zhǎng),吸收的熱量更多,因此溫升逐漸增大。由圖8還可看出,當(dāng)渦輪壓降恒定時(shí),隨著太陽(yáng)輻射強(qiáng)度的增加,系統(tǒng)的入射能量增加,氣流速度和溫升均顯著增大。

當(dāng)太陽(yáng)輻射強(qiáng)度為600W/m2時(shí),改變渦輪壓降可以得到集熱棚內(nèi)的熱損失,如圖9所示。由圖9可以看到,不論何種壓降,通過(guò)集熱棚頂?shù)臒釗p失總是遠(yuǎn)遠(yuǎn)大于通過(guò)地面的熱損失。當(dāng)渦輪壓降逐漸增大時(shí),通過(guò)集熱棚頂和通過(guò)地面的熱損失均逐漸增加。這是因?yàn)楫?dāng)渦輪壓降增大時(shí),渦輪對(duì)氣流的阻礙效應(yīng)增強(qiáng),系統(tǒng)流量減小,集熱棚內(nèi)空氣溫度升高,集熱棚內(nèi)、外溫差增大,從而導(dǎo)致熱損失增加。因此,集熱棚覆蓋采用雙層玻璃會(huì)大大降低集熱棚頂?shù)臒釗p失,提高集熱棚的集熱效率。

圖9 不同渦輪壓降下集熱棚的熱損失

圖10 不同渦輪壓降下的系統(tǒng)性能

太陽(yáng)輻射強(qiáng)度為600W/m2時(shí),改變渦輪壓降后的系統(tǒng)性能如圖10所示。隨著渦輪壓降的增加,系統(tǒng)各部分的熱損失增大,集熱棚的集熱效率逐漸下降。但是,由于渦輪壓降增加會(huì)引起流量減小,使得通過(guò)煙囪出口的動(dòng)能損失大大降低,因此系統(tǒng)總效率逐漸上升,在壓降為120Pa時(shí)達(dá)到最大值。繼續(xù)增大渦輪壓降,煙囪產(chǎn)生的抽力部分用來(lái)驅(qū)動(dòng)渦輪,但剩余部分并不足以克服系統(tǒng)內(nèi)的流動(dòng)損失。

4 結(jié) 論

本文建立了太陽(yáng)能煙囪發(fā)電系統(tǒng)的數(shù)值模型,并對(duì)其內(nèi)部的傳熱問(wèn)題進(jìn)行了數(shù)值模擬。結(jié)果表明,集熱棚內(nèi)的輻射換熱是其內(nèi)部傳熱過(guò)程的重要組成部分,不能忽略,否則會(huì)導(dǎo)致集熱棚內(nèi)溫度分布背離實(shí)際的溫度分布。渦輪壓降對(duì)系統(tǒng)性能也有重要影響,隨著渦輪壓降的增大,煙囪內(nèi)的氣流溫升逐漸增大,集熱棚內(nèi)的熱損失逐漸增大,集熱棚效率下降。由于煙囪內(nèi)氣流速度減小,出口的余速損失減小,系統(tǒng)效率逐漸提高。此外,對(duì)集熱棚內(nèi)的熱損失分析發(fā)現(xiàn),通過(guò)集熱棚頂表面的對(duì)流換熱及長(zhǎng)波輻射的損失是系統(tǒng)內(nèi)的主要熱損失。

[1] HAAF W, FRIEDRICH K, MAYR G, et al.Solar chimneys: part Ⅰ Principle and construction of the pilot plant in Manzanares [J].International Journal of Sustainable Energy, 1983, 2(1): 3-20.

[2] 代彥軍, 黃海賓, 王如竹.太陽(yáng)能熱風(fēng)發(fā)電技術(shù)應(yīng)用于寧夏地區(qū)的研究 [J].太陽(yáng)能學(xué)報(bào), 2003, 24(3): 408-412.DAI Yanjun, HUANG Haibin, WANG Ruzu.Utilization of solar chimney power plant into Ningxia Hui Autonomous Region [J].Acta Energiae Solaris Sinica, 2003, 24(3): 408-412.

[3] SALAH L, AMOR B, TOUFIK C.Performance analysis of a solar chimney power plant in the southwestern region of Algeria [J].Renewable and Sustainable Energy Reviews, 2010, 14(1): 470-477.

[4] ZHOU X P, WANG F, FAN J, et al.Performance of solar chimney power plant in Qinghai-Tibet Plateau [J].Renewable and Sustainable Energy Reviews, 2010, 14(8): 2249-2255.

[5] 張楚華.大型太陽(yáng)能煙囪發(fā)電站熱力計(jì)算與分析 [J].可再生能源雜志, 2007, 25(2): 3-6.ZHANG Chuhua.Thermodynamic analysis and calculation of large-scale solar chimney electricity generation plant [J].Renewable Energy Resources, 2007, 25(2): 3-6.

[6] PASUMARTHI N, SHERIF S A.Experimental and theoretical performance of a demonstration solar chimney model: part I Mathematical model development [J].International Journal of Energy Research, 1998, 22(3): 277-288.

[7] BERNARDES M A D, VOβ A, WEINREBE G.Thermal and technical analyses of solar chimneys [J].Solar Energy, 2003, 75(6): 511-524.

[8] ZHOU X P, YANG J K, XIAO B, et al.Simulation of a pilot solar chimney thermal power generating equipment [J].Renewable Energy, 2007, 32(10): 1637-1644.

[9] LI J Y, GUO P H, WANG Y.Effects of collector radius and chimney height on power output of a solar chimney power plant with turbines [J].Renewable Energy, 2012, 47: 21-28

[10]PASTOHR H, KORNADT O, GURLEBECK K.Numerical and analytical calculations of the temperature and flow field in the upwind power plant [J].International Journal of Energy Research, 2004, 28(6): 495-510.

[11]MING T Z, LIU W, XU G L, et al.Numerical simulation of the solar chimney power plant systems coupled with turbine [J].Renewable Energy, 2008, 33(5): 897-905.

[12]XU G L, MING T Z, PAN Y, et al.Numerical analysis on the performance of solar chimney power plant system [J].Energy Conversion and Management, 2011, 52(2): 876-883.

[13]HAAF W.Solar chimneys: part II Preliminary test results from the Manzanares pilot plant [J].International Journal of Sustainable Energy, 1984, 2(2): 141-161.

[14]PRETORIUS J P, KR?GER D G.Solar chimney power plant performance [J].ASME Journal of Solar Energy Engineering, 2006, 128(3): 302-311.

[15]LI JY, GUO PH, WANG Y.Preliminary investigation of a novel solar and wind energy extraction system [J].Proceedings of The Institution of Mechanical Engineers: Part A Journal of Power and Energy, 2012;226(A1): 73-85.

[16]SWINBANK W C.Long-wave radiation from clear skies [J].Quarterly Journal of the Royal Meteorological Society, 1963, 89(381): 339-348.

[17]SCHLAICH J, BERGERMANN R, SCHIEL W, et al.Design of commercial solar updraft tower systems: utilization of solar induced convective flows for power generation [J].ASME Journal of Solar Energy Engineering, 2005, 127(1): 117-124.

[18]楊世銘, 陶文銓.傳熱學(xué) [M].北京: 高等教育出版社, 1998.

(編輯 苗凌)

NumericalAnalysisofHeatTransferinSolarChimneyPowerPlant

GUO Penghua,WANG Yuan,LI Jingyin

(Department of Fluid Machinery and Engineering, Xi’an Jiaotong University, Xi’an 710049, China)

The research was conducted to investigate the influences of the radiation heat transfer, solar radiation and turbine pressure drop on the performance of a solar chimney power plant.The comprehensive 3D numerical simulations were carried out to examine the heat transfer process in a Spanish prototype.The heat losses in the system were considered in the numerical simulations incorporating the radiation model, solar load model and turbine model and the experimental data for the Spanish prototype were used to verify the simulations.The results show that the radiation heat transfer plays an important role in the heat transfer process inside the collector and the heat loss from the collector roof to the ambient is the major heat loss in the system.In addition, the collector efficiency decreases while the system efficiency gradually increases when the turbine pressure drop increases under the constant solar radiation condition.

solar chimney; numerical analysis; radiation heat transfer; collector efficiency; turbine pressure drop

10.7652/xjtuxb201403019

2013-05-16。

郭朋華(1986—),男,博士生;王元(通訊作者),男,教授,博士生導(dǎo)師。

國(guó)家自然科學(xué)基金資助項(xiàng)目(11272252,51276137);中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金資助項(xiàng)目。

時(shí)間: 2013-12-19

O121.8;G558

:A

:0253-987X(2014)03-0102-06

網(wǎng)絡(luò)出版地址: http:∥www.cnki.net/kcms/detail/61.1069.T.20131219.1121.002.html