張亮， 張安龍， 曲平平， 荊宇燕

(燕山大學(xué)車輛與能源學(xué)院， 秦皇島 066004)

如今能源緊缺問題日趨嚴(yán)重，如何合理地利用能源已經(jīng)成為當(dāng)今各國工業(yè)發(fā)展的核心問題，為了緩解能源緊張的問題，世界上各個(gè)國家都在尋找新能源或者節(jié)約能源的方法，而換熱器在能量傳遞方面發(fā)揮著重要作用。換熱器中直接參與換熱的部分為換熱管，它的換熱效率往往能夠直接影響到換熱器的效率。波壁管作為一種新型換熱管，與傳統(tǒng)光管相比，具有容易產(chǎn)生流動(dòng)分離，能獨(dú)立支撐等優(yōu)勢，引起中外學(xué)者的廣泛研究。曾敏等[1]指出，波壁管內(nèi)流體在流動(dòng)換熱過程中，存在一個(gè)臨界雷諾數(shù)，當(dāng)雷諾數(shù)Re大于此臨界值時(shí)，在相同條件下，波壁管的換熱性能比光滑圓管更好。Nishimura等[2-3]指出在相同功率下，波壁管內(nèi)流體的熱量與質(zhì)量傳遞效果最好發(fā)生在中等雷諾數(shù)范圍內(nèi)。隨著科學(xué)技術(shù)的不斷發(fā)展，熱交換設(shè)備的熱負(fù)荷與傳熱強(qiáng)度不斷增大，傳統(tǒng)的換熱工質(zhì)不足以滿足換熱設(shè)備的傳熱與冷卻要求，而通過向純液體中加入納米顆粒從而形成的納米流體能夠很好地代替?zhèn)鹘y(tǒng)工質(zhì)。Choi[4]第一次提出了納米流體的概念后，大量學(xué)者對此進(jìn)行了研究。研究[5-7]表明，納米流體的導(dǎo)熱系數(shù)比傳統(tǒng)流體的導(dǎo)熱系數(shù)更大，且隨體積濃度的增加而增大。丁雅勤等[8]指出，Ag-乙二醇納米流體比純乙二醇的導(dǎo)熱系數(shù)更大，當(dāng)納米流體中銀顆粒的體積分?jǐn)?shù)為0.06%時(shí)，Ag-乙二醇納米流體的導(dǎo)熱系數(shù)與純乙二醇相比提高了36%。丁潔等[9]對混合基納米流體在汽車散熱器中的穩(wěn)定性及傳熱特性進(jìn)行了研究，研究結(jié)果表明，與基液相比，納米流體的傳熱速率明顯提高。張勝寒等[10]指出，向基礎(chǔ)流體中添加納米顆粒會(huì)導(dǎo)致其黏度增加，且納米流體的黏度隨納米顆粒的濃度增大而增加，同時(shí)發(fā)現(xiàn)親水型納米顆粒比疏水型納米顆粒與水形成的納米流體的黏度高。李強(qiáng)等[11]提出，向水中加入體積分?jǐn)?shù)為2%的Cu納米顆粒，與基液相比，Cu-水納米流體的h提高了60%左右。Chandraprabu等[12]研究了Al2O3-水納米流體與CuO-水納米流體的傳熱性能，發(fā)現(xiàn)兩者的傳熱速率均比基礎(chǔ)流體更高。Srinivas等[13]研究了三種水基納米流體(Al2O3、CuO和TiO2/水)在螺旋管殼式換熱器中的自然對流換熱，發(fā)現(xiàn)這三種水基納米流體均能提高換熱器的對流換熱系數(shù)。

基于此，現(xiàn)主要研究波壁管中Cu-乙醇(Ethanol，ET)/乙二醇(ethylene glycol，EG)/丙二醇(propylene glycol，PG)三種納米流體的強(qiáng)化傳熱，分析三種納米流體在體積分?jǐn)?shù)φ為1.5%與3%情況下波壁管內(nèi)流體的對流換熱系數(shù)與摩擦阻力系數(shù)以及其綜合換熱性能?？梢詾镃u-醇基納米流體在換熱管內(nèi)的應(yīng)用提供重要的理論基礎(chǔ)。

1 幾何模型與納米流體物性參數(shù)

1.1 幾何模型

圖1為波壁管的二維結(jié)構(gòu)圖，波壁管的最大直徑Dmax=10 mm，最小直徑Dmin=3 mm，波長λ=14 mm，波幅A=3.5 mm，兩端各留有一段長為14 mm的直管壁，總長度為L=224 mm，圖2為波壁管的三維模型圖。

圖1 波壁管二維結(jié)構(gòu)圖Fig.1 Two dimensional structure diagram of wave wall tube

圖2 波壁管三維模型圖Fig.2 Three dimensional model of wave wall tube

1.2 納米流體物性參數(shù)

納米流體與常規(guī)流體相比具有更大的導(dǎo)熱系數(shù)，現(xiàn)運(yùn)用分子動(dòng)力學(xué)的方法對Cu-ET/EG/PG 三種納米流體的導(dǎo)熱系數(shù)κ進(jìn)行了計(jì)算。首先計(jì)算了ET、EG、PG三種常規(guī)流體的導(dǎo)熱系數(shù)κ并與實(shí)驗(yàn)進(jìn)行了對比，兩者的對比如表1所示。發(fā)現(xiàn)計(jì)算值與實(shí)驗(yàn)值的誤差均在5%以內(nèi)，故使用該分子動(dòng)力學(xué)的方法來對納米流體的導(dǎo)熱系數(shù)的計(jì)算是可行的。

對φ為1.5%與3%的Cu-ET/EG/PG納米流體的導(dǎo)熱系數(shù)計(jì)算結(jié)果如表2所示，其余的物性參數(shù)可由式(1)～式(3)計(jì)算出。由表2可知，納米流體的導(dǎo)熱系數(shù)比基液更大，且隨著φ的增大導(dǎo)熱系數(shù)逐漸增大。在φ為1.5%與3%的情況下，與基液相比，Cu-ET納米流體的導(dǎo)熱系數(shù)κ增加了4.8%與7.84%，Cu-EG納米流體的κ增加了0.6%和2.3%，而Cu-PG納米流體的κ增加了5.82%與14.4%。

表1 計(jì)算值與實(shí)驗(yàn)值的對比Table 1 Comparison between calculated and experimental values

表2 Cu-ET/EG/PG 納米流體的導(dǎo)熱系數(shù)Table 2 Thermal conductivity of Cu-ET/EG/PG nanofluids

納米流體的密度可表示為

ρnf=φρp+(1-φ)ρf

(1)

納米流體的黏度可表示為

(2)

納米流體的比熱容可表示為

ρnfcnf=(1-φ)ρfcf+φρpcp

(3)

式中：φ為納米流體的體積分?jǐn)?shù)；ρnf、μnf、cnf分別為納米流體的密度、黏度、比熱容；ρf、μf、cf分別為基液的密度、黏度、比熱容；ρp、cp分別為Cu納米粒子的密度與比熱容。

2 邊界條件與網(wǎng)格無關(guān)性驗(yàn)證

2.1 邊界條件

由于本文研究的Cu納米顆粒的粒徑比較小，粒徑為0.7 nm，且φ低于5%，故可將納米流體視為單相流[16]。壓力耦合求解采用SIMPLE算法，邊界條件設(shè)為：波壁管入口采用速度入口，入口溫度設(shè)為300 K，出口采用自由出口，壁面為固定無滑移壁面，溫度采用恒溫，設(shè)為350 K。

2.2 網(wǎng)格無關(guān)性驗(yàn)證

由于波壁管的曲率較大，所以采用適應(yīng)性更強(qiáng)的非結(jié)構(gòu)性網(wǎng)格，為了消除網(wǎng)格數(shù)量對計(jì)算結(jié)果的影響，故需要進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證。劃分了6組不同數(shù)量的網(wǎng)格，在相同條件下對6組不同網(wǎng)格數(shù)量進(jìn)行數(shù)值計(jì)算，計(jì)算結(jié)果對流換熱系數(shù)隨網(wǎng)格數(shù)量的變化如圖3所示。發(fā)現(xiàn)網(wǎng)格數(shù)量在達(dá)到20萬左右時(shí)，對流換熱系數(shù)的變化不明顯，故本文數(shù)值計(jì)算的網(wǎng)格數(shù)量為20萬左右。

圖3 網(wǎng)格無關(guān)性驗(yàn)證Fig.3 Grid independence verification

3 結(jié)果分析

3.1 溫度分析

圖4為Cu-PG納米流體在波壁管內(nèi)的溫度場，其中圖4(a)～圖4(c)分別為丙二醇基液、體積分?jǐn)?shù)為1.5%、體積分?jǐn)?shù)為3%Cu-PG納米流體的溫度云圖。圖5為波壁管通道中心的Cu-PG納米流體溫度隨流體流動(dòng)方向的變化。綜合分析圖4和圖5發(fā)現(xiàn)，波壁管內(nèi)Cu-PG納米流體的平均溫度比基液丙二醇的平均溫度更高，且隨著φ的增加，納米流體的平均溫度逐漸增大，基液與φ為1.5%和3%的納米流體在波壁管出口處平均溫度分別為310.8 K、311.3 K和312 K。這是因?yàn)榧{米流體與基液相比有更大的導(dǎo)熱系數(shù)，所以在傳熱過程中能夠更快地通過恒定的壁面溫度將納米流體加熱到更高的溫度。

圖4 Cu-PG納米流體的溫度云圖Fig.4 Temperature nephogram of Cu-PG nanofluids

圖5 Cu-PG納米流體的溫度變化Fig.5 Temperature variation of Cu-PG nanofluids

3.2 對流換熱系數(shù)分析

圖6為Cu-ET/EG/PG納米流體的對流換熱系數(shù)h隨Re與φ的變化，從圖6可以看出，隨著Re的增大，波壁管內(nèi)的流體的h逐漸增大；即Re越大，流體的流速越大，波壁管內(nèi)流體的擾動(dòng)更加劇烈，從而能夠強(qiáng)化換熱。Cu-ET/EG/PG納米流體的對流換熱系數(shù)均大于基液的對流換熱系數(shù)，Cu-ET/EG納米流體的h略有增加，而Cu-PG納米流體的h則增加明顯；Cu-ET/EG/PG納米流體的h隨φ的增加而逐漸增大。這表明，向基液中添加納米顆粒，可以強(qiáng)化流體的換熱，適當(dāng)?shù)脑黾蛹{米顆粒的φ也可以使流體的換熱加強(qiáng)。

圖6 h隨Re與φ的變化Fig.6 Variation of h with Re and φ

3.3 摩擦阻力系數(shù)分析

在傳熱強(qiáng)化理念中，強(qiáng)化傳熱不應(yīng)該以摩擦阻力大幅度增加為代價(jià)，故需要對摩擦阻力系數(shù)進(jìn)行分析，摩擦阻力系數(shù)f可表示為

(4)

式(4)中：Dmax為波壁管的最大直徑，mm；ΔP為波壁管進(jìn)出口的平均壓降，Pa；ρ為流體密度，kg/m3；u為流體的平均速度，m/s；L為波壁管的管長，m。

圖7 f隨Re與φ的變化Fig.7 Variation of f with Re and φ

波壁管中Cu-ET/EG/PG納米流體的摩擦阻力系數(shù)f隨Re與φ的變化關(guān)系如圖7所示，由圖7可知，隨著Re的增大，波壁管內(nèi)的f在逐漸減小，且減小的幅度在逐漸變小。在相同Re下，Cu-ET/EG/PG納米流體的f比基液大，且隨著φ的增大，納米流體的f逐漸增大。這說明，向基液中添加納米顆粒改善流體換熱的同時(shí)，也會(huì)導(dǎo)致阻力損失的增大，且該損失隨著φ的增加而增大。

3.4 綜合性能評(píng)價(jià)分析

通過對波壁管內(nèi)流體的h與f的分析發(fā)現(xiàn)，向基液中添加較小體積分?jǐn)?shù)φ的納米顆?？梢蕴岣吡黧w的h，但也會(huì)相對應(yīng)地增加f，增大流體流動(dòng)過程中的流動(dòng)損失，相對應(yīng)地會(huì)增大功率的消耗。而強(qiáng)化傳熱的理念是在強(qiáng)化傳熱的同時(shí)不能消耗太多的功率，故引用性能評(píng)價(jià)因子PEC[17]來對波壁管內(nèi)流體的h與f進(jìn)行綜合分析。

PEC公式可表示為

(5)

式(5)中：hn與h0分別為納米流體的對流換熱系數(shù)與相對應(yīng)的基液的對流換熱系數(shù)；fn與f0分別為Cu-ET/EG/PG納米流體的摩擦阻力系數(shù)與相對應(yīng)的基液的摩擦阻力系數(shù)。

圖8 PEC隨Re與φ的變化Fig.8 Variation of PEC with Re and φ

PEC隨Re與φ的變化規(guī)律如圖8所示，從圖8可以看出，Cu-ET/EG納米流體的PEC均小于1，且PEC隨著φ的增加而下降，這表明與基液相比，Cu-ET/EG納米流體的綜合換熱性能沒有得到加強(qiáng)；這是因?yàn)橄蚧褐刑砑蛹{米顆粒后，Cu-ET/EG兩種納米流體的導(dǎo)熱系數(shù)的增加不大，進(jìn)而納米流體的換熱系數(shù)h增加較小，而摩擦阻力系數(shù)f增加較大，所以導(dǎo)致其PEC小于1。Cu-PG納米流體的PEC大于1，且隨著納米顆粒的φ的增大，PEC逐漸增大。這表明向基液中添加納米顆粒后，波壁管內(nèi)納米流體的綜合換熱性能得到加強(qiáng)，且隨著φ的增大Cu-PG納米流體的綜合換熱效果更強(qiáng)。在Cu-ET/EG/PG納米流體中，Cu-PG納米流體的綜合換熱效果最好。

4 結(jié)論

(1)通過分子動(dòng)力學(xué)對Cu-ET/EG/PG納米流體的導(dǎo)熱系數(shù)的數(shù)值計(jì)算，發(fā)現(xiàn)納米流體與基液相比，其導(dǎo)熱系數(shù)更大，且隨著φ的增大，導(dǎo)熱系數(shù)逐漸增大；在Cu-ET/EG/PG納米流體中，Cu-PG納米流體的導(dǎo)熱系數(shù)增幅最大。

(2)通過對波壁管內(nèi)流體的溫度分析，發(fā)現(xiàn)波壁管內(nèi)納米流體的平均溫度比基液的平均溫度更高，且隨著φ的增大平均溫度逐漸增加。

(3)在相同Re下，與基液相比，Cu-ET/EG納米流體的h略有增加，而Cu-PG納米流體的h則明顯增加；隨著φ的增大，Cu-ET/EG/PG納米流體的對流換熱系數(shù)逐漸增加；納米流體在強(qiáng)化傳熱的同時(shí)，也會(huì)增加摩擦阻力損失，相同條件下，Cu-ET/EG/PG納米流體的摩擦阻力系數(shù)f比基液更大，且隨著φ的增加，納米流體的f也在不斷增加。

(4)Cu-ET/EG納米流體的綜合性能評(píng)價(jià)因子PEC小于1，且隨著φ的增加而減小，而Cu-PG納米流體的PEC大于1，隨著φ的增加而增大。在Cu-ET/EG/PG三種納米流體中，Cu-PG納米流體的綜合換熱性能更好。