槽式太陽能集熱器傳熱模型及性能分析

陳玉英

摘 要：

槽式太陽能集熱器一維和二維傳熱數(shù)學(xué)模型是一組非線性代數(shù)方程，為改進(jìn)求解的穩(wěn)定性和計(jì)算精度，將槽式太陽能集熱器一維和二維傳熱模型的求解看作有約束優(yōu)化問題，建立了集熱器傳熱過程求解的有約束優(yōu)化數(shù)學(xué)模型，應(yīng)用MATLAB軟件優(yōu)化函數(shù)fmincon進(jìn)行求解。分析了傳熱流體入口溫度及太陽能輻射熱流密度變化對集熱器性能的影響。采用fmincon函數(shù)求解集熱器傳熱過程，計(jì)算速度快，計(jì)算過程穩(wěn)定。分析表明，傳熱流體溫度變化對集熱器效率的影響大于太陽能輻射熱流密度對集熱器效率的影響。

關(guān)鍵詞：

槽式集熱器；太陽能；傳熱；MATLAB

中圖分類號：TK512

文獻(xiàn)標(biāo)志碼：A 文章編號：16744764（2016）04005306

槽式聚焦太陽能集熱器是一種高效的太陽能熱利用集熱器，廣泛用于太陽能熱力發(fā)電、物料干燥，海水淡化等領(lǐng)域[15]。槽式太陽能集熱器性能涉及聚焦器、集熱器的光學(xué)性質(zhì)、集熱器復(fù)雜的傳熱過程、集熱器結(jié)構(gòu)等著諸多方面，就其傳熱模型，傳熱性能來說，學(xué)者們對其進(jìn)行了不同程度的研究。文獻(xiàn)[6] 基于熱平衡方法，建立了對吸熱器直接加熱的實(shí)驗(yàn)裝置，對集熱器熱損失進(jìn)行了研究。由于對拋物面聚焦集熱器整體進(jìn)行實(shí)驗(yàn)研究受到多方面的限制，而基于熱平衡法的實(shí)驗(yàn)研究缺乏對包括聚光器、傳熱流體在內(nèi)的集熱器整體性能的認(rèn)識，基于集熱器傳熱模型建立與求解的數(shù)值模擬分析方法受到人們的關(guān)注[713]。文獻(xiàn)[78]對槽式太陽能集熱器建立了一維、二維穩(wěn)態(tài)傳熱的數(shù)學(xué)模型，分別應(yīng)用EES工程方程求解器及MATLAB軟件符號函數(shù)solve進(jìn)行了求解，文獻(xiàn)[910]應(yīng)用MATLAB軟件非線性方程組求解函數(shù)fsolve對集熱器一維穩(wěn)態(tài)傳熱模型進(jìn)行了求解。

槽式太陽能集熱器一維和二維穩(wěn)態(tài)傳熱數(shù)學(xué)模型是一組非線性代數(shù)方程，在集熱器實(shí)際工作中，集熱管和玻璃護(hù)管壁面溫度存在一定依賴關(guān)系，這里將集熱器傳熱數(shù)學(xué)模型的求解看做優(yōu)化問題，通過集熱器熱平衡方程和傳熱方程，建立求解集熱器傳熱過程的有約束優(yōu)化數(shù)學(xué)模型，應(yīng)用MATLAB軟件優(yōu)化工具箱函數(shù)fmincon進(jìn)行求解。fmincon函數(shù)內(nèi)嵌多種優(yōu)化方法，可使求解過程更穩(wěn)定，同時(shí)計(jì)入集熱器各部分溫度的依賴關(guān)系，使計(jì)算結(jié)果更合理。本文在求解傳熱方程的基礎(chǔ)上，分析了部分因素對集熱器性能的影響。

1 集熱管傳熱模型

1.1 集熱管一維傳熱模型

集熱管傳熱模型及溫度分布分別如圖1和圖2所示。集熱管主要由外層玻璃護(hù)管、金屬吸熱管組成。玻璃管和吸熱管之間為真空，金屬管內(nèi)為傳熱介質(zhì)，這里以水為傳熱介質(zhì)。設(shè)金屬管內(nèi)流體編號為1；金屬管內(nèi)外壁編號分別為2和3；玻璃管內(nèi)外壁編號為4和5；玻璃管外空氣編號為6；天空編號為7。單位管長的傳熱量用′（W/m）表示。玻璃管吸收的太陽輻射熱量為q′5SolAbs，吸熱管吸收的太陽輻射熱量為q′3SolAbs。集熱管正常工作時(shí)，金屬管外表面溫度最高[67]，如圖2所示。因此，對于集熱管來說，以金屬吸熱管外表面為界，熱量分別向吸熱管內(nèi)和吸熱管外傳遞，如圖1所示。

1.2 集熱管熱平衡方程

根據(jù)能量守恒原理，傳熱流體、吸熱管和玻璃管的熱平衡關(guān)系可用式（1）～（6）表示。

對于傳熱流體和吸熱管，傳熱流體與吸熱管內(nèi)表面的對流換熱量等于通過吸熱管的導(dǎo)熱量。

對于吸熱管外表面，吸熱管吸收的太陽輻射熱量等于通過吸熱管壁的導(dǎo)熱量、吸熱管外表面與玻璃管內(nèi)表面的對流換熱量及輻射換熱量之和。

對于玻璃管內(nèi)表面，玻璃管內(nèi)表面與吸熱管外表面通過對流和輻射交換的熱量等于通過玻璃管壁的導(dǎo)熱量。

對于玻璃管外表面，玻璃管外表面吸收的太陽輻射熱量與通過玻璃管壁的導(dǎo)熱量之和等于玻璃管外表面與周圍環(huán)境的對流換熱量與玻璃管外壁與天空的輻射換熱量之和。

1.3 集熱管二維傳熱模型

在一維傳熱模型中，換熱流體溫度可以按兩種方式考慮，一是取流體平均溫度作為流體的整體溫度，不考慮溫度沿管長的變化；二是僅考慮流體進(jìn)出口溫度的變化，認(rèn)為沿管長流體溫度呈線性變化。在二維模型中，將集熱管分成n段，如圖3所示，在每一段的邊界上溫度是連續(xù)的。與一維傳熱模型一樣，假設(shè)輻射傳熱熱流沿管長均勻分布，且垂直于管壁表面，在每一小段上，縱向溫度假設(shè)是線性的。

1.4 集熱管傳熱計(jì)算

1.4.1 傳熱流體與吸熱管之間的對流傳熱

由牛頓冷卻公式，從吸收管內(nèi)表面到傳熱流體的對流傳熱為

1.4.2 吸熱管與玻璃護(hù)管的傳熱

在吸熱管與玻璃護(hù)管間通過對流和輻射進(jìn)行換熱。對于對流換熱，當(dāng)環(huán)形腔內(nèi)壓力較低（壓力<133.32 Pa） 時(shí)，傳熱機(jī)理是分子導(dǎo)熱，在較高壓力下（壓力>133.32 Pa），傳熱過程是自然對流。當(dāng)環(huán)形區(qū)域內(nèi)壓力<133.32 Pa時(shí)，對流換熱系數(shù)按式（13）計(jì)算[78）。

1.4.3 玻璃護(hù)管與環(huán)境的對流換熱 玻璃護(hù)管與外界空氣的對流換熱，在無風(fēng)時(shí)為自然對流換熱，有風(fēng)時(shí)為強(qiáng)迫對流換熱。自然對流換熱時(shí)努塞爾數(shù)通過式（22）～（24）計(jì)算[78，14]。

有風(fēng)時(shí)，玻璃護(hù)管與空氣的對流換熱按流體橫掠單管時(shí)的經(jīng)驗(yàn)公式計(jì)算。通過吸熱管及玻璃管導(dǎo)熱計(jì)算按通過圓筒壁導(dǎo)熱公式計(jì)算，吸熱管外壁面與玻璃管內(nèi)壁面的輻射換熱按兩無限長環(huán)形空腔輻射傳熱公式計(jì)算，玻璃護(hù)管外表面與天空輻射換熱按無限大封閉面與內(nèi)包凸表面物體輻射換熱計(jì)算[13，15]。

2 集熱管傳熱過程求解

集熱管傳熱計(jì)算中，通常給定太陽總輻射熱流密度Isolar，傳熱介質(zhì)入口溫度t1，in，環(huán)境溫度t6，天空溫度t7，待定參數(shù)包括傳熱介質(zhì)出口溫度t1，out，吸熱管內(nèi)外壁表面平均溫度t2及t3，玻璃護(hù)管內(nèi)外壁表面平均溫度t4及t5，共5個(gè)未知量，可通過式（1）～（4）、式（8）以及相應(yīng)的傳熱公式來進(jìn)行求解。由于輻射熱流量與溫度的4次方成正比，因此，引入傳熱計(jì)算式后，熱平衡方程是關(guān)于溫度的非線性代數(shù)方程。文獻(xiàn)[7]應(yīng)用工程方程求解器（EES）來求解這一組代數(shù)方程。MATLAB軟件有豐富的數(shù)學(xué)庫函數(shù)及優(yōu)化工具箱函數(shù)，給非線性代數(shù)方程組的求解帶來很大方便，文獻(xiàn)[910]分別應(yīng)用MATLAB軟件提供的代數(shù)方程組符號求解函數(shù)solve及fsolve來求解這一組代數(shù)方程。這里采用MATLAB軟件提供的求解約束優(yōu)化問題的函數(shù)fmincon來求解這一組非線性代數(shù)方程[16]。與solve及fsolve函數(shù)相比，采用約束優(yōu)化函數(shù)fmincon進(jìn)行求解時(shí)程序設(shè)計(jì)簡單，計(jì)算速度更快，穩(wěn)定性更好。

3 計(jì)算結(jié)果及分析

3.1 計(jì)算參數(shù)

算例采用文獻(xiàn)[8]給出的集熱管參數(shù)，各部分結(jié)構(gòu)參數(shù)如表1所示。

3.2 計(jì)算結(jié)果與實(shí)測結(jié)果比較

傳熱計(jì)算分別按一維和二維模型進(jìn)行，對于二維傳熱模型，將集熱管分為12段。在太陽輻射熱流密度為均勻的情況下，一維和二維傳熱模型結(jié)果相同。模擬計(jì)算結(jié)果與文獻(xiàn)[8]實(shí)驗(yàn)結(jié)果如表2所示。

3.3 傳熱流體入口溫度對集熱器性能的影響

傳熱流體體積流量V=0.226 8 m3/h，環(huán)境溫度t6=26.6 ℃，外界風(fēng)速u5=3.8 m/s，到達(dá)集熱槽開口面單位長度的太陽輻射熱量I′=758 W/m時(shí)，傳熱流體入口溫度與集熱器熱效率關(guān)系曲線如圖4所示。

從圖4可以看出，隨著吸熱管傳熱流體入口溫度的增加，流體從太陽輻射吸收的熱量Φgain呈下降趨勢，且流體入口溫度越高，下降的程度越明顯，集熱器效率隨傳熱流體溫度的增加也呈下降趨勢。這是因?yàn)?，隨著傳熱流體入口溫度的增加，吸熱管管壁溫度增加，在外界氣溫和風(fēng)速不變的情況下，玻璃套管壁面溫度增加有限，使得吸熱管與玻璃護(hù)管間輻射換熱量Φ34rad增加，熱損失增加，因此，集熱器效率下降。在吸熱管與玻璃護(hù)管環(huán)形空腔內(nèi)真空度較高的情況下，吸熱管與玻璃護(hù)管間對量換熱量較少，熱損失主要由輻射散熱產(chǎn)生。吸熱管與玻璃護(hù)管間輻射換熱量通過玻璃護(hù)管的導(dǎo)熱及護(hù)管外表面與環(huán)境的對流換熱和輻射換熱散出。

3.4 太陽輻射熱流密度對集熱器性的能的影響

傳熱流體體積流量V=0.226 8 m3/h，傳熱流體入口溫度tin=200 ℃，環(huán)境溫度t6=26.6 ℃，外界風(fēng)速u5=3.8 m/s時(shí)，太陽輻射熱流密度與集熱器熱效率關(guān)系曲線如圖5示。

隨著太陽輻射熱流密度的增加，集熱器效率呈增加的趨勢，由于散熱損失（如玻璃護(hù)管外壁與環(huán)境的對流換熱量Φ56conv）也呈增加的趨勢，因此，集熱器效率隨太陽輻射熱流密度的增加其增加的幅度較平緩。傳熱流體溫度變化對集熱器效率的影響大于太陽輻射對集熱器效率的影響。

4 結(jié) 論

在建立槽式集熱器一維和二維穩(wěn)態(tài)傳熱模型的基礎(chǔ)上，應(yīng)用MATLAB軟件優(yōu)化工具箱函數(shù)fmincon對熱平衡方程和傳熱方程進(jìn)行了求解，并對部分因素對集熱管性能的影響進(jìn)行了分析，通過模擬計(jì)算和分析可得出以下結(jié)論。

1）將集熱器傳熱模型按優(yōu)化問題進(jìn)行求解可以加快求解速度，提高求解的穩(wěn)定性。

2）隨傳熱流體入口溫度的增加，集熱器效率顯著下降；隨太陽輻射熱流密度的增加，集熱器效率呈平穩(wěn)增加的趨勢。

3）在吸熱管與玻璃護(hù)管環(huán)形空腔內(nèi)真空度較高的情況下，吸熱管的散熱主要是輻射換熱；集熱管熱效率對傳熱流體溫度變化的影響要高于太陽輻射熱流密度變化對集熱器效率的影響。

參考文獻(xiàn)：

[1]

黃素逸， 黃樹紅， 許國良. 太陽能熱發(fā)電原理及技術(shù)[M]. 北京： 中國電力出版社，2012.

HUANG S Y， HHUANG S H， XU G L. Principle and technology of solar thermal power generation[M]. Beijing： China Electric Power Press，2012.（in Chinese）

[2] 王金平，王軍，張耀明，等. 槽式太陽能聚光集熱器傳熱特性分析[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2015，31（7）： 185192.

WANG J P， WANG J， ZHANG Y M，et al. analysis of heat transfer characteristics for parabolic trough solar collector[J]. Transactions of the Chinese Society of Agricultural Engineering， 2015，31（7）： 185192. （in Chinese）

[3] 徐偉， 楊帆， 劉靜靜，等. 槽式拋物面太陽能熱發(fā)電系統(tǒng)發(fā)展現(xiàn)狀與探討[J]. 發(fā)電設(shè)備，2015，21（1）： 7478.

XU W， YANG F， LIU J J， et al. Review and discussion on parabolic trough solar power systems[J]， Power Equipment， 2015，21（1） ： 7478. （in Chinese）

[4] 楊敏林， 楊曉西， 林汝謀， 等. 太陽能熱發(fā)電技術(shù)與系統(tǒng)[J]. 熱能動力工程， 2008， 23（3）： 221228.

YANG M L， YANG X X， LIN R M. Solar energybased thermal power generation technologies and their systems[J]. Journal of Engineering for Thermal Energy & Power， 2008， 23（3）： 221228.（in Chinese）

[5] FERNANDEZGARC A， ZARZA E， VALENZUELA L，et al. Parabolictrough solar collectors and their applications[J]. Renewable and Sustainable Energy Reviews， 2010，14： 16951721.

[6] 熊亞選， 吳玉庭， 馬重芳，等. 槽式太陽能聚光集熱器光熱性能初步實(shí)驗(yàn)研究[J]. 太陽能學(xué)報(bào)， 2012， 33（12） ： 20872092.

XIONG Y X，WU Y T ，MA C F，et al. Preliminary experiment study on the thermal performance of a parabolic trough collector[J]. Acta Energiae Solaris Sinica， 2012， 33（12） ： 20872092.（in Chinese）

[7] FORRISTAL R. Heat transfer analysis and modeling of a parabolic trough solar receiver implemented in engineering equation solver[R]. Colorado： National Renewable Energy Laboratory， 2003.

[8] ALFELLAG， MOHANAD ABDULAZEEZ ABDULRAHEEM. Modeling and experimental investigation of parabolic trough solar collector[D]. Florida： EmbryRiddle Aeronautical University， 2014 .

[9] 楊潤， 李斌. 基于MATLAB 的拋物槽式太陽能集熱器集熱管一維穩(wěn)態(tài)模型[J]. 可再生能源，2008，26（4）：69.

YANG R， LI B. The calculation of the 1D steadystate model of the parabolic trough solar absorber pipe based on MATLAB [J]. Renewable Energy Resources， 2008，26（4）：69.（in Chinese）

[10] 韋彪， 朱天宇， 劉德有. 槽式DSG太陽能集熱系統(tǒng)模擬分析[J]. 工程熱物理學(xué)報(bào)，2012，33（3）： 473476.

WEI B， ZHU T Y， LIU D Y. Simulation of heat transfer for parabolic trough DSG solar collectors [J]. Journal of Engineering Thermophysics， 2012， 33（3）：473476 .（in Chinese）

[11] WANG K， HE Y L， CHENG Z D. A design method and numerical study for a new type parabolic trough solar collector with uniform solar flux distribution [J]. Science China Technological Sciences， 2014， 57（3）：531540.

[12] MC. EL JAI，F(xiàn)Z. CHALQI. A modified model for parabolic trough solar receiver [J]. American Journal of Engineering Research， 2013，2（5）：200211.

[13] 周傳文. 槽式太陽集熱器熱性能研究[D]. 北京：華北電力大學(xué)，2013.

ZHOU C W. Thermal performance analysis of trough solar collector[D]. Beijing： North China Electric Power University， 2013.（in Chinese）

[14] WU Z Y， LI S D， YUAN G F，et al. Threedimensional numerical study of heat transfer characteristics of parabolic trough receiver [J]. Applied Energy， 2014，113： 902911.

[15] 楊世銘，陶文銓. 傳熱學(xué)[M].4版. 北京：高等教育出版社，2006 .

YANG S M， TAO W Q. Heat transfer [M].4th Edition Beijing： Higher Education Press， 2006.（in Chinese）

[16] 張永恒. 工程優(yōu)化設(shè)計(jì)與MATLAB實(shí)現(xiàn)[M]. 北京：清華大學(xué)出版社，2011.

ZHANG Y H. Engineering optimization design with MATLAB implementation[M]. Beijing： Tsinghua University Press， 2011.（in Chinese）

（編輯 胡玲）