三角形腔體內(nèi)納米流體自然對流換熱數(shù)值研究

馬文強,馬兵善,王　剛

(蘭州理工大學(xué) 土木工程學(xué)院,甘肅 蘭州　730050)

?

三角形腔體內(nèi)納米流體自然對流換熱數(shù)值研究

馬文強,馬兵善,王剛

(蘭州理工大學(xué) 土木工程學(xué)院,甘肅 蘭州730050)

摘要采用有限容積法對二維三角形腔體內(nèi)CuO-機油納米流體自然對流換熱進(jìn)行了數(shù)值模擬研究,重點分析研究了液體納米層厚度對CuO-機油納米流體自然對流換熱的影響,同時也研究了瑞利數(shù)Ra、納米顆粒體積分?jǐn)?shù)φ以及腔體高寬比AR對CuO-機油納米流體自然對流換熱的影響。研究結(jié)果表明,液體納米層的存在使納米流體的自然對流增強,換熱量增大；在相同納米顆粒體積分?jǐn)?shù)下,隨著Ra數(shù)的增大,自然對流換熱強度顯著增強,且Ra數(shù)較小時,換熱量隨著腔體高寬比AR的增大而減小。

關(guān)鍵詞納米流體；自然對流換熱；三角形腔體；數(shù)值模擬

20世紀(jì)90年代以來,研究人員開始探索將納米材料技術(shù)應(yīng)用于強化換熱領(lǐng)域,研究新一代高效換熱冷卻技術(shù)。美國Argonne國家實驗室的Choi[1]在1995年首次提出了納米流體的概念,之后許多學(xué)者在納米流體的熱物性、納米流體強制對流換熱、自然對流換熱和沸騰換熱等方面開展了大量的研究工作。納米流體與傳統(tǒng)傳熱介質(zhì)相比,不僅具有較高的導(dǎo)熱系數(shù),而且相對于毫米和微米級顆粒,懸浮液穩(wěn)定性好、不會堵塞或磨損管道、也不會引起系統(tǒng)壓降的過大增加。但是,作為一種新型的強化傳熱介質(zhì),要將納米流體成功應(yīng)用于工業(yè),研究其換熱性能是非常必要的。

在實際工程應(yīng)用領(lǐng)域中,如在研究溫室型太陽房[2]、屋面太陽能蒸餾器、三角形形狀的太陽能集熱器及微型電子元件冷卻等問題時,經(jīng)常要對三角形封閉腔體內(nèi)的自然對流換熱現(xiàn)象進(jìn)行分析研究。許多學(xué)者已對此類問題進(jìn)行了廣泛深入地研究。Lei等[3]采用實驗方法和數(shù)值模擬方法研究了等腰三角形腔體內(nèi)水的自然對流換熱,研究發(fā)現(xiàn)等腰三角形腔體傾斜壁面和底面的平均Nu數(shù)近似正比于Gr0.2。Karyakin等[4]以及Fuad Kent等[5]采用數(shù)值模擬方法研究了腔體高寬比AR及Ra數(shù)對直角三角形腔體內(nèi)空氣自然對流換熱的影響。Ghasemi等[6]采用數(shù)值模擬方法研究了布朗運動、瑞利數(shù)Ra、納米顆粒體積分?jǐn)?shù)φ、熱源位置d以及腔體高寬比AR對直角三角形腔體內(nèi)CuO-水納米流體的自然對流換熱的影響。Yu等[7]同時采用實驗方法和數(shù)值模擬方法研究了等腰三角形內(nèi)CuO-水納米流體的自然對流換熱。

Yu等[8]提出了新的納米流體導(dǎo)熱系數(shù)模型,指出納米流體由流體分子、納米顆粒及包裹在納米顆粒表面,類似于固體的液體納米層組成。在納米顆粒與基液間的傳熱過程中,液體納米層扮演著熱橋的角色,其厚度也會影響傳熱。以下將重點分析研究液體納米層厚度對三角形腔體內(nèi)CuO-機油納米流體自然對流換熱的影響,同時也研究瑞利數(shù)Ra、納米顆粒體積分?jǐn)?shù)φ以及腔體高寬比AR對CuO-機油納米流體自然對流換熱的影響。

1物理模型與控制方程

1.1 物理模型

采用的物理模型如圖1所示,該模型為二維直角三角形腔體,腔體豎直方向(高)長度為H(m),水平方向長度(底邊)為L(m),加熱部位表面維持相對較高溫度Th,其中心距坐標(biāo)系原點的距離為d(m),長度為w(m),而傾斜壁面維持相對較低的溫度Tc(Th>Tc),其余壁面均為絕熱壁面。腔體中間充滿CuO-機油納米流體,納米顆粒與機油的熱物理參數(shù)如表1所列。

圖1　物理模型與坐標(biāo)系統(tǒng)Fig.1　Physical model and coordinate system

物理性質(zhì)ρ/(kg·m-3)Cp/(J·kg-1·K-1)k/(W·m-1·K-1)β/×106K-1α/×107(m2·s-1)機油884.119090.14570085.9CuO6500536205157.5

1.2控制方程

假設(shè)納米顆粒的大小形狀均一,腔體內(nèi)的納米流體為不可壓縮、各向同性、做層流運動的牛頓流體,且基液與納米顆粒處于熱平衡狀態(tài)并以同樣的速度流動,同時引入Boussinesq假設(shè)來考慮由于浮升力作用導(dǎo)致的密度變化。描述此穩(wěn)態(tài)、不可壓縮、二維自然對流換熱問題的控制方程無量綱形式為

(1)

(2)

(3)

(4)

式(1)~式(4)中涉及到的無量綱參數(shù)定義為

(5)

納米流體熱物性參數(shù)的計算公式列于表2。式(1)~式(5)及表2各式中,物性參數(shù)符號的下標(biāo)f、p、nf分別代表水、納米顆粒和納米流體；L為特征長度；ρ、α、β、ν、μ、k分別為密度、熱擴散系數(shù)、熱膨脹系數(shù)、運動粘度、動力粘度和導(dǎo)熱系數(shù)。表2中φ為納米顆粒體積分?jǐn)?shù)；η為納米顆粒表面形成的液體納米層厚度與納米顆粒半徑之比[8]。

表2　納米流體熱物性參數(shù)的計算公式

描述上述問題的無量綱邊界條件為

(6)

其中:D為熱源距水平壁面的無量綱特征長度,D=d/L；W為熱源的無量綱特征長度,W=w/L。

熱壁面表面的幾何平均努塞爾數(shù)Num定義為

(7)

2數(shù)值計算方法驗證

對無量綱控制方程(1)~(4)的離散采用有限體積法,采用SIMPLEC算法[9]處理速度場與壓力場的耦合。為降低假擴散的影響,對流項采用具有較高精度的QUICK格式進(jìn)行離散。

2.1網(wǎng)格獨立性驗證

網(wǎng)格的劃分也會影響計算結(jié)果的精度,網(wǎng)格越密,迭代過程中就會占用越多的計算機資源,使計算機的計算速度減慢。同時網(wǎng)格數(shù)增加到一定程度時還會引起計算誤差的累計,產(chǎn)生偏離。選用20×20、40×40、60×60、80×80和100×100五套網(wǎng)格進(jìn)行獨立性考察,腔體內(nèi)的介質(zhì)為CuO-機油納米流體,Ra取106,見表3,從表3中可以看出,選用80×80的網(wǎng)格可以消除網(wǎng)格數(shù)對模擬結(jié)果的影響,網(wǎng)格劃分如圖2所示。

表3　網(wǎng)格獨立性考察

圖2　均分網(wǎng)格示意圖Fig.2　Sketch map of uniform mesh

2.2計算程序的驗證

為使計算結(jié)果準(zhǔn)確、可靠,采用體積份額φ=3%的CuO-水納米流體驗證了計算程序的可靠性,并與文獻(xiàn)[4]中所得結(jié)果做了對比(見表4),通過表4中與模擬值的對比結(jié)果可知,熱壁面的平均Num與文獻(xiàn)[4]中所給出的數(shù)據(jù)吻合較好。

3模擬結(jié)果分析

數(shù)值計算針對Ra分別為103、104、105、106時進(jìn)行求解,φ取值分別為0.05、0.10、0.15及0.20,熱壁面溫度Th=310 K,冷壁面溫度Tc=290 K。熱源無量綱特征長度W=0.3。在研究腔體高寬比AR(AR=H/L)對自然對流換熱的影響時,AR的取值為0.5、1.0、1.5和2.0。

表4　模擬值與文獻(xiàn)值的對比

3.1Ra數(shù)對自然對流換熱的影響

圖3為η=0.1、AR=1、D=0.5及φ=0.2時不同Ra數(shù)下的流函數(shù)圖與等溫線圖。從圖3中可以看出,當(dāng)Ra數(shù)較小時,腔體內(nèi)順時針旋轉(zhuǎn)的漩渦中心位于腔體中部,最大無量綱流函數(shù)值(Ψmax,nf)較小,等溫線近似與豎直壁面平行。隨著Ra數(shù)的不斷增大,漩渦中心逐漸上移,最大無量綱流函數(shù)值迅速增大；同時,等溫線形狀也發(fā)生變化,等溫線逐漸由豎直變得平坦,當(dāng)Ra達(dá)到106時,等溫線幾乎與豎直壁面垂直。這主要是因為在較低Ra數(shù)時,腔體內(nèi)浮升力較小,換熱方式以導(dǎo)熱為主；隨著Ra數(shù)的增加,腔體內(nèi)浮升力加大,是自然對流換熱加強所致。然而,在熱壁面附近等溫線始終與豎直壁面平行,這說明在邊界層內(nèi)壁面與納米流體間的換熱方式主要以導(dǎo)熱為主。

圖3　不同Ra數(shù)下CuO-機油納米流體的流函數(shù)圖(上)與等溫線圖(下)(η=0.1,AR=1,D=0.5,φ=0.2)Fig.3　Flow function diagram and isothermal diagram of CuO-engine oil nano-fluid in presence of different Ra

3.2納米流體層厚度對自然對流換熱的影響

圖4為AR=1、D=0.5時η=0和η=0.1兩種情況下平均Num數(shù)的變化趨勢。從圖4中能夠很直觀地看出,Num數(shù)隨著納米顆粒體積分?jǐn)?shù)的增加而增大,幾乎呈線性變化。不難發(fā)現(xiàn),液體納米層存在時的平均Num數(shù)均高于不考慮液體納米層時的平均Num數(shù)。納米顆粒體積分?jǐn)?shù)較低時平均Num數(shù)提高的幅度較小,而當(dāng)納米顆粒體積分?jǐn)?shù)增大時平均Num數(shù)提高的幅度加大,當(dāng)Ra分別為104、106,納米顆粒體積分?jǐn)?shù)φ=0.05時,熱壁面的平均Num數(shù)分別提高了4.07%和3.35%;納米顆粒體積分?jǐn)?shù)φ=0.2時,平均Num數(shù)分別提高了18.1%和11.7%。這是因為納米顆粒表面形成的液體納米層構(gòu)建了納米顆粒與液體之間的傳熱橋梁,液體納米層的存在使得顆粒與流體間的換熱增強,納米流體導(dǎo)熱系數(shù)增大,自然對流換熱加強。

3.3腔體高寬比AR對自然對流換熱的影響

圖5為納米顆粒體積分?jǐn)?shù)φ=0.2、D=0.5、η=0,Ra分別為103、104、105及106時,三角形腔體熱壁面平均Num數(shù)隨腔體高寬比AR的變化。從圖5中可以看出,當(dāng)Ra分別為103、104時,隨著AR的增大,Num數(shù)呈現(xiàn)減小的趨勢,當(dāng)AR=2時,Num數(shù)分別減小了27.9%、21.7%。而當(dāng)Ra分別為105、106時,Num數(shù)隨著AR的增大而增大,說明腔體內(nèi)換熱方式由導(dǎo)熱為主變?yōu)閷α鳛橹?自然對流換熱加強,對流換熱量增加,Num數(shù)增大。

圖4　η=0(虛線)和η=0.1(實線)兩種情況下平均Num隨納米顆粒體積分?jǐn)?shù)的變化Fig.4　Change of average Num along nano-particle volumefraction when η=0(dotted line) and η=0.1(full line)

圖5　不同Ra數(shù)下平均Num數(shù)隨AR的變化Fig.5　Change of average Num along AR in presence of different Ra

4結(jié)論

通過對三角形腔體內(nèi)CuO-機油納米流體自然對流換熱進(jìn)行數(shù)值研究,重點分析了液體納米層的存在對CuO-機油納米流體自然對流換熱的影響,同時也研究了瑞利數(shù)Ra、納米顆粒體積分?jǐn)?shù)φ以及腔體高寬比AR對CuO-機油納米流體自然對流換熱的影響。研究結(jié)果表明,液體納米層的存在使自然對流換熱過程中的換熱量增大。當(dāng)Ra數(shù)較小時,自然對流相對較弱,隨著Ra數(shù)的增加,腔體內(nèi)浮升力加大,致使自然對流強度逐漸加大。當(dāng)Ra分別為103、104時,隨著AR的增大,Num數(shù)呈現(xiàn)減小的趨勢,而當(dāng)Ra分別為105、106時,Num數(shù)隨著AR的增大而增大,說明腔體內(nèi)換熱方式由導(dǎo)熱為主變?yōu)閷α鳛橹?自然對流換熱加強,對流換熱量增加,Num數(shù)增大。

參考文獻(xiàn):

[1]Choi S U S.Enhancing Thermal Conductivity of Fluids with Nanoparticles[C]//Developments and Applications of Non-Newtonian Flows,New York,1995:99-105.

[2]王德芳.陽光間和房屋之間的自然對流換熱計算方法[J].甘肅科學(xué)學(xué)報,1991,3(2):1-5.

[3]Lei C W,Armfield S W,Patterson J C.Unsteady Natural Convection in a Water-Filled Isosceles Triangular Enclosure Heated from Below[J].Heat Mass Transfer,2008,51(11-12):2 637-2 650.

[4]Karyakin Y E,Sokovishin Y A,Martynenko O G.Transient Natural Convection in Triangular Enclosures[J].International Journal of Heat and Mass Transfer,1988,31(9):1 759-1 766.

[5]Fuad Kent E,Asmaz E,Ozerbay S.Laminar Natural Convection in Right Triangular Enclosures[J].Heat and Mass Transfer,2007,44(2):187-200.

[6]Ghasemi B,Aminossadati S M.Brownian Motion of Nanoparticles in a Triangular Enclosure with Natural Convection[J].Thermal Sciences,2010,49(6):931-940.

[7]Yu Z T,Xu X,Hu Y C,etal.Numerical Study of Transient Buoyancy-Driven Convective Heat Transfer of Water-Based Nanofluids in a Bottom-Heated Isosceles Triangular Enclosure[J].Heat Mass Transfer,2011,54(1-3):526-532.

[8]Yu W,Choi S U S.The Role of Interfacial Layers in the Enhanced Thermal of Nanfluids:a Renovated Maxwell Model[J].Nanoparticle Research,2003,5(1-2):167-171.

[9]陶文銓.數(shù)值換熱學(xué)[M].第2版.西安:西安交通大學(xué)出版社,2001.

Study on Natural Convection Heat Exchange Value of Nano-Fluid in Triangular Cavity

Ma Wenqiang,Ma Bingshan,Wang Gang

(SchoolofCivilEngineering,LanzhouUniversityofTechnology,Lanzhou730050,China)

AbstractNumerical simulation research was implemented for the natural convection heat exchange of the CuO-engine oil nano-fluid in the two-dimensional triangular cavity by using the finite volume method,and the influence of the thickness of the fluid nanolayer to the natural convection heat exchange of the CuO-engine oil nano-fluid was analyzed and researched intensively,as well as the influence of Rayleigh numberRa,nano-particle volume fractionφand ratioARof height to width of the cavity to the natural convection heat exchange of the CuO-engine oil nano-fluid.The research result showed that the natural convection of the nano-fluid should be strengthened and the amount of heat exchange should be increased in presence of fluid nanolayer;the natural convection heat exchange strength might be strengthened obviously along increasingRaon the premise of identical nano-particle volume fraction;and the amount of heat exchange might be reduced along increasing ratioARof height to width of the cavity when theRawas small.

Key wordsNano-fluid;Natural convection heat exchange;Triangular cavity;Numerical simulation

中圖分類號:TK124

文獻(xiàn)標(biāo)志碼:A

文章編號:1004-0366(2016)01-0001-05

作者簡介:馬文強(1991-)，男，甘肅臨夏人，碩士，研究方向為納米流體自然對流換熱.E-mail:lut_mwq@sina.com.

基金項目:國家自然科學(xué)基金項目(51266006)；國家科技支撐計劃項目(2011BAJ03B08)；蘭州理工大學(xué)博士基金項目(BS04200903).

收稿日期:2015-01-20;修回日期:2015-05-03.

doi:10.16468/j.cnki.issn1004-0366.2016.01.001.

引用格式:Ma Wenqiang,Ma Bingshan,Wang Gang.Study on Natural Convection Heat Exchange Value of Nano-Fluid in Triangular Cavity[J].Journal of Gansu Sciences,2016,28(1):1-5.[馬文強,馬兵善,王剛.三角形腔體內(nèi)納米流體自然對流換熱數(shù)值研究[J].甘肅科學(xué)學(xué)報,2016,28(1):1-5.]