Fluent輻射傳熱模型理論以及相關(guān)設(shè)置

/Fluent輻射傳熱模型理論以與相關(guān)設(shè)置目錄1 概述 22 基礎(chǔ)理論 22.1 專業(yè)術(shù)語(yǔ)解釋： 22.2 FLUENT輻射模型介紹： 22.3 輻射模型適用范圍總結(jié) 23 Fluent實(shí)際案例操作 23.1 Case1-測(cè)試externalemissivity使用DO模型計(jì)算-2D模型 23.2 Case2-測(cè)試internalemissivity-使用DO模型計(jì)算-2D模型 23.3 仿真結(jié)論 2概述在傳熱的仿真中，有時(shí)候會(huì)不可避免的涉與到輻射傳熱，而我們對(duì)Fluent中輻射模型的了解甚少，很難得到可靠的計(jì)算結(jié)果。因此，一直以來(lái)，F(xiàn)luent中的帶輻射的傳熱仿真是我們的一個(gè)難點(diǎn)，本專題重點(diǎn)來(lái)學(xué)習(xí)輻射模型的理論，讓我們對(duì)輻射計(jì)算模型有一個(gè)深入的了解，以幫助我們攻克這個(gè)仿真難點(diǎn)?；A(chǔ)理論專業(yè)術(shù)語(yǔ)解釋：在Fluent中開(kāi)啟輻射模型時(shí)，流體介質(zhì)以與固體壁面會(huì)出現(xiàn)一些專業(yè)的參數(shù)需要用戶來(lái)設(shè)置。在Fluenthelp中介紹輻射模型時(shí)會(huì)經(jīng)常提到一些專業(yè)術(shù)語(yǔ)。對(duì)這些專業(yè)參數(shù)以與術(shù)語(yǔ)，我們來(lái)一一解釋：1、Opticalthickness（光學(xué)深度，無(wú)量綱量）：介質(zhì)層不透明性的量度。即介質(zhì)吸收輻射的能力的量度，等于入射輻射強(qiáng)度與出射輻射強(qiáng)度之比。設(shè)入射到吸收物質(zhì)層的入射輻射強(qiáng)度為I，透射的輻射強(qiáng)度為e，則T=I/e，其中T為光學(xué)深度。按照此定義，那介質(zhì)完全透明，對(duì)輻射不吸收、也不散射，透射的輻射強(qiáng)度e=入射輻射強(qiáng)度I，即光學(xué)深度為T(mén)=1，介質(zhì)不參與輻射?！园俣劝倏贫鳩LUENT中T=αL，其中L為介質(zhì)的特征長(zhǎng)度，α為輻射削弱系數(shù)（可理解為介質(zhì)因吸收和散射引起的光強(qiáng)削弱系數(shù)）。如果T=0，說(shuō)明介質(zhì)不參與輻射，和百度百科中的定義有出入。但是所表達(dá)的意思是接近的，一個(gè)是前后輻射量的比值；一個(gè)是變化量和入射輻射量的比值（根據(jù)Fluenthelp里的解釋，經(jīng)過(guò)介質(zhì)的輻射損失量=I*T，個(gè)人理解，按照此定義，T不可能大于1啊，矛盾。//TheoryGuide::0//5.HeatTransfer//5.3.ModelingRadiation//5.3.2.RadiativeTransferEquation）。該問(wèn)題的解釋為：其實(shí)一點(diǎn)也不矛盾，如果Opticalthickness=1，就說(shuō)明輻射在經(jīng)過(guò)一定特征長(zhǎng)度L的介質(zhì)后被完全吸收。如果>1，就說(shuō)明輻射根本穿透不了特征長(zhǎng)度L的介質(zhì)，而被早早吸收完了。打個(gè)比方，Opticalthickness=10，說(shuō)明輻射在經(jīng)過(guò)L/10距離后已經(jīng)被吸收（或散射）完。其中α=αA+αS；2、AbsorptionCoefficient（αA吸收系數(shù)，單位1/m，見(jiàn)圖2-1）：因?yàn)榻橘|(zhì)吸收而導(dǎo)致的輻射強(qiáng)度在經(jīng)過(guò)每單位長(zhǎng)度介質(zhì)后改變的量。空氣作為流體介質(zhì)時(shí)，一般不吸收熱輻射，該系數(shù)可近視設(shè)為0。而當(dāng)氣體中水蒸氣和CO2含量較高時(shí)，那對(duì)輻射的系數(shù)就不能忽略了。3、ScatteringCoefficient（αS散射系數(shù)，單位1/m）：因?yàn)榻橘|(zhì)散射而導(dǎo)致的輻射強(qiáng)度在經(jīng)過(guò)每單位長(zhǎng)度介質(zhì)后改變的量?？諝庾鳛榱黧w介質(zhì)時(shí)，一般情況下，該系數(shù)可近視設(shè)為0。對(duì)于含顆粒物的流體，散射作用不容忽視。4、RefractiveIndex（折射系數(shù)，無(wú)量綱量）：介質(zhì)中的光速和真空中的光速之比。如是空氣，可近視設(shè)為1（默認(rèn)值）。一般對(duì)于具有方向性的輻射源問(wèn)題，比如LED發(fā)光或激光等光學(xué)傳熱問(wèn)題，輻射在經(jīng)過(guò)水以與玻璃等透明介質(zhì)時(shí)，需要設(shè)置該參數(shù)。一般情況，熱輻射在計(jì)算域中是往各個(gè)方向輻射的，各項(xiàng)同性，沒(méi)有方向性，該參數(shù)設(shè)為1即可。圖2-1介質(zhì)的輻射相關(guān)參數(shù)設(shè)置5、DiffuseReflection（漫反射）：輻射到不透明固體表面的能量，一部分被固體吸收，另一部分被反射，其中反射分為鏡面反射和漫反射。6、SpecularReflection（鏡面反射）：7、InternalEmissivity（內(nèi)部發(fā)射率）：處于計(jì)算域中的couplewall，solid和fluidzone或者solid和solidzone或者fluid和fluidzone之間的輻射率。8、ExternalEmissivity（外部發(fā)射率）：處于計(jì)算邊界上wall，外部環(huán)境和wall之間的輻射率。對(duì)于基于灰體輻射假設(shè)的計(jì)算，灰體輻射率不隨波長(zhǎng)變化，灰體輻射率=吸收率；9、ThetaDivisionandPhiDivision：設(shè)置為2，可作為初步估算；為了得到更為準(zhǔn)確的結(jié)果，最少設(shè)置成3，甚至為5，F(xiàn)luent13.0默認(rèn)值為4。10、ThetaPixelsandPhiPixels：對(duì)于灰體輻射，默認(rèn)值1*1足夠了；但是對(duì)于涉與到對(duì)稱面、周期性邊界、鏡面反射、半透明邊界時(shí)，需設(shè)置為3*3；FLUENT輻射模型介紹：Fluent中有五種輻射計(jì)算模型，各個(gè)模型的使用范圍以與其優(yōu)缺點(diǎn)分別為：1、DTRM模型：優(yōu)勢(shì)：模型相對(duì)簡(jiǎn)單，可以通過(guò)增加射線數(shù)量來(lái)提高計(jì)算精度，適用于光學(xué)深度范圍非常廣的各種輻射問(wèn)題。限制：1）模型假設(shè)所有面都是漫反射，意味著輻射的反射相對(duì)于入射角是各項(xiàng)同性的，無(wú)鏡面反射。2）忽略散射作用。3）灰體輻射假設(shè)。4）使用大數(shù)目射線求解問(wèn)題，非常耗費(fèi)CPU資源。5）和非一致網(wǎng)格（non-conformalinterface）、滑移網(wǎng)格（slidingmesh）不能一起使用，不能用并行計(jì)算。2、P1模型；：優(yōu)勢(shì)：相比DTRM模型，P1模型耗費(fèi)自己資源更少，并且考慮了散射作用；對(duì)于光學(xué)深度較大的燃燒模型，P1模型更穩(wěn)定。P1模型使用曲線中uobiao比較容易處理復(fù)雜幾何的輻射問(wèn)題。限制：1）假設(shè)所有面都是漫反射，和DTRM相同。2）使用與灰體和非灰體輻射問(wèn)題。3）如果光學(xué)深度很小時(shí)，模型計(jì)算精度取決于幾何的復(fù)雜性。4）對(duì)于局部熱源以與散熱片問(wèn)題，該模型會(huì)夸大輻射傳熱量。3、Rossland模型：優(yōu)勢(shì)：相對(duì)P1模型。它不求解額外的關(guān)于入射輻射的傳輸方程，因此比P1模型耗資源要少。限制：只能用于光學(xué)深度比較大的情況，推薦用于光學(xué)深度大于3的情況；不能用于密度求解器，只能用于壓力求解器。4、Surface-to-Surface（S2S）輻射模型；優(yōu)勢(shì)：非常適用于封閉空間中沒(méi)有介質(zhì)的輻射問(wèn)題，（如航天器、太陽(yáng)能搜集系統(tǒng)、輻射供熱裝置等）；限制：1）所有面都是漫反射。2）灰體輻射假設(shè)。3）在表面增加時(shí)，耗費(fèi)計(jì)算資源大幅增加。4）不能用于介質(zhì)參與的輻射問(wèn)題（participatingradiation）。5）不能和周期性邊界、對(duì)稱邊界、非一致網(wǎng)格交界面、網(wǎng)格自適應(yīng)一起使用。5、DO模型優(yōu)勢(shì)：適用于所有光學(xué)深度范圍的輻射問(wèn)題；既能求解S2S的無(wú)介質(zhì)封閉區(qū)域問(wèn)題，也能求解介質(zhì)參與的輻射問(wèn)題。適用于灰體、非灰體、漫反射、鏡面反射以與半透明介質(zhì)的輻射。輻射模型適用范圍總結(jié)DTRM和DO模型幾乎可適用于所有光學(xué)深度問(wèn)題，相比之下，DO模型的范圍更廣。光學(xué)深度>1，可用P1和Rossland模型；而>3時(shí)，Rossland模型比較合適。對(duì)于光學(xué)深度<1的問(wèn)題，只能用DTRM和DO模型。S2S適用于光學(xué)深度為0的問(wèn)題，即流體介質(zhì)不參與輻射的問(wèn)題?？偨Y(jié)：一般關(guān)于空氣對(duì)流輻射的問(wèn)題，屬于光學(xué)深度=0的問(wèn)題，因此可使用DTRM、S2S、DO模型，在ICEPAK解決輻射問(wèn)題就有這三個(gè)模型的選項(xiàng)（在13.0版本中才加入DTRM和DO模型）。Fluent實(shí)際案例操作從簡(jiǎn)單的2Dcase入手，在實(shí)際操作中真正搞清楚emissivity和absorptioncoefficient的含義，以與Fluent中solid和fluidzone之間的輻射傳熱機(jī)理。Case1-測(cè)試externalemissivity使用DO模型計(jì)算-2D模型2D模型，直徑2m，externalradiationtemperature400K，圓形為solid，恒溫300K圖3-1溫度場(chǎng)分布圖圖3-2輻射換熱設(shè)置設(shè)置externalemissivity1，計(jì)算出外界對(duì)wall輻射傳熱功率為6230.3188W，根據(jù)理論公式計(jì)算：Pra=5.67e-8*1*3.14*2*（400^4-300^4）=6231W。仿真結(jié)果和理論計(jì)算非常接近。將externalemissivity設(shè)成0.5，計(jì)算出輻射傳熱功率為3114.6W。改變internalemissivity的值，計(jì)算值不變。從以上仿真結(jié)果可知：1、2.1小結(jié)的第八點(diǎn)externalemissivity的解釋是正確的，輻射傳熱基于灰體假設(shè)，輻射系數(shù)等于吸收系數(shù)。Case2-測(cè)試internalemissivity-使用DO模型計(jì)算-2D模型1）Solid（Al）-solid（Steel）-solid（Al）-caseExternalemissivityInternalemissivityExternalemissivityInternalemissivity圖3-3從里到外Solid（Al）-solid（Steel）-solid（Al）i）InternalsolidFixtemperature=400℃，externalradiationtemperature=300，externalemissivtiy=1；internalemissivtiy=1：圖3-4溫度分布以與換熱量ii）internalemissivtiy=0：圖3-5溫度分布以與換熱量從圖4、5可知，上下兩張圖的溫度分布非常相近，上圖中溫度稍高，而zone之間的換熱量存在差異，將internalemissivity改為0，代表兩個(gè)不同材料的zone之間輻射傳熱量為零，因此總傳熱量從5555W降低至5055W。可知，F(xiàn)luent中認(rèn)為緊密相連的兩個(gè)solidzone（存在couplewall）之間是存在輻射傳熱的（也可設(shè)置為無(wú)輻射傳熱），相當(dāng)于實(shí)際情況中的兩個(gè)物體的接觸面，只不過(guò)在Fluent中未設(shè)置接觸熱阻?？偨Y(jié)：實(shí)際情況中有接觸熱阻，有輻射傳熱；Fluent中無(wú)接觸熱阻，有輻射傳熱。用Fluent一般不進(jìn)行涉與接觸熱阻細(xì)節(jié)的仿真。2）Solid（Al）-fluid（air）-solid（Al），nogravity-casei），externalemissivtiy=1；internalemissivtiy=1，fluid的absorptioncoefficient=0；圖3-6溫度分布以與換熱量中心400K的solid往externalsolid的輻射傳熱功率為：Pra=5.67e-8*1*3.14*1*（400^4-335^）=2315W，和fluentreprot值2333W(包含了空氣熱傳導(dǎo)的功率)比較接近；ii）internalemissivtiy=0，fluid的absorptioncoefficient=0；圖3-7溫度分布以與換熱量將internalemissivtiy=0后，傳熱功率下降為21W，說(shuō)明無(wú)輻射換熱時(shí)，僅靠空氣導(dǎo)熱的傳熱功率非常小。iii）Fluid和externalsolid之間的internalemissivtiy=1，fluid的absorptioncoefficient=1；圖3-8溫度分布以與換熱量iv）Fluid和externalsolid之間的internalemissivtiy=0，fluid的absorptioncoefficient=1；圖3-9溫度分布以與換熱量v）Fluid和externalsolid之間的internalemissivity=1，fluidabsorptioncoefficient=0，externalsolidabsorptioncoefficient=1，conductivity=0.02；vi）emissivity=1，fluidabsorptioncoefficient=0，externalsolidabsorptioncoefficient=10，conduct