Fluent輻射傳熱模型理論以及相關(guān)設(shè)置
2016-09-04 by:CAE仿真在線 來(lái)源:互聯(lián)網(wǎng)
1、 概述Fluent高級(jí)工程實(shí)例視頻
在傳熱的仿真中,有時(shí)候會(huì)不可避免的涉及到輻射傳熱,而我們對(duì)Fluent中輻射模型的了解甚少,很難得到可靠的計(jì)算結(jié)果。因此,一直以來(lái),Fluent中的帶輻射的傳熱仿真是我們的一個(gè)難點(diǎn),本專題重點(diǎn)來(lái)學(xué)習(xí)輻射模型的理論,讓我們對(duì)輻射計(jì)算模型有一個(gè)深入的了解,以幫助我們攻克這個(gè)仿真難點(diǎn)。
2、基礎(chǔ)理論
2.1 專業(yè)術(shù)語(yǔ)解釋:
在Fluent中開(kāi)啟輻射模型時(shí),流體介質(zhì)以及固體壁面會(huì)出現(xiàn)一些專業(yè)的參數(shù)需要用戶來(lái)設(shè)置。
在Fluent help中介紹輻射模型時(shí)會(huì)經(jīng)常提到一些專業(yè)術(shù)語(yǔ)。
對(duì)這些專業(yè)參數(shù)以及術(shù)語(yǔ),我們來(lái)一一解釋:
1、Optical thickness(光學(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=1,介質(zhì)不參與輻射?!园俣劝倏?/span>
而FLUENT中T=αL,其中L為介質(zhì)的特征長(zhǎng)度,α為輻射削弱系數(shù)(可理解為介質(zhì)因吸收和散射引起的光強(qiáng)削弱系數(shù))。如果T=0,說(shuō)明介質(zhì)不參與輻射,和百度百科中的定義有出入。但是所表達(dá)的意思是接近的,一個(gè)是前后輻射量的比值;一個(gè)是變化量和入射輻射量的比值(根據(jù)Fluent help里的解釋,經(jīng)過(guò)介質(zhì)的輻射損失量 =I*T,個(gè)人理解,按照此定義,T不可能大于1啊,矛盾。// Theory Guide :: 0 // 5. Heat Transfer // 5.3. Modeling Radiation // 5.3.2. Radiative Transfer Equation)。該問(wèn)題的解釋為:其實(shí)一點(diǎn)也不矛盾,如果Optical thickness =1,就說(shuō)明輻射在經(jīng)過(guò)一定特征長(zhǎng)度L的介質(zhì)后被完全吸收。如果>1,就說(shuō)明輻射根本穿透不了特征長(zhǎng)度L的介質(zhì),而被早早吸收完了。打個(gè)比方,Optical thickness=10,說(shuō)明輻射在經(jīng)過(guò)L/10距離后已經(jīng)被吸收(或散射)完。
其中α=αA+αS;
2、Absorption Coefficient(αA吸收系數(shù),單位1/m,見(jiàn)圖2-1):因?yàn)榻橘|(zhì)吸收而導(dǎo)致的輻射強(qiáng)度在經(jīng)過(guò)每單位長(zhǎng)度介質(zhì)后改變的量??諝庾鳛榱黧w介質(zhì)時(shí),一般不吸收熱輻射,該系數(shù)可近視設(shè)為0。而當(dāng)氣體中水蒸氣和CO2含量較高時(shí),那對(duì)輻射的系數(shù)就不能忽略了。
3、Scattering Coefficient(α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、Refractive Index(折射系數(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、Diffuse Reflection(漫反射):輻射到不透明固體表面的能量,一部分被固體吸收,另一部分被反射,其中反射分為鏡面反射和漫反射。
6、Specular Reflection(鏡面反射):
7、Internal Emissivity(內(nèi)部發(fā)射率):處于計(jì)算域中的couple wall,solid和 fluid zone或者solid和solid zone 或者 fluid和fluid zone 之間的輻射率。
8、External Emissivity(外部發(fā)射率):處于計(jì)算邊界上wall,外部環(huán)境和wall之間的輻射率。對(duì)于基于灰體輻射假設(shè)的計(jì)算,灰體輻射率不隨波長(zhǎng)變化,灰體輻射率=吸收率;
9、Theta Division and Phi Division:設(shè)置為2,可作為初步估算;為了得到更為準(zhǔn)確的結(jié)果,最少設(shè)置成3,甚至為5,Fluent13.0默認(rèn)值為4。
10、Theta Pixels and Phi Pixels:對(duì)于灰體輻射,默認(rèn)值1*1足夠了;但是對(duì)于涉及到對(duì)稱面、周期性邊界、鏡面反射、半透明邊界時(shí),需設(shè)置為3*3;
2.2 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-conformal interface)、滑移網(wǎng)格(sliding mesh)不能一起使用,不能用并行計(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)題(participating radiation)。
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ì)的輻射。
2.3 輻射模型適用范圍總結(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)題。
總結(jié):一般關(guān)于空氣對(duì)流輻射的問(wèn)題,屬于光學(xué)深度=0的問(wèn)題,因此可使用DTRM、S2S、DO模型,在ICEPAK解決輻射問(wèn)題就有這三個(gè)模型的選項(xiàng)(在13.0版本中才加入DTRM和DO模型)。
3、Fluent實(shí)際案例操作
從簡(jiǎn)單的2D case入手,在實(shí)際操作中真正搞清楚 emissivity 和 absorption coefficient的含義,以及Fluent中 solid和fluid zone之間的輻射傳熱機(jī)理。
3.1 Case1-測(cè)試external emissivity使用DO模型計(jì)算-2D模型
2D模型,直徑2m,external radiation temperature 400K,圓形為solid,恒溫300K
圖3-1 溫度場(chǎng)分布圖
圖3-2 輻射換熱設(shè)置
設(shè)置external emissivity 1,計(jì)算出外界對(duì)wall輻射傳熱功率為6230.3188W,根據(jù)理論公式計(jì)算:
Pra=5.67e-8*1*3.14*2*(400^4-300^4)=6231W。
仿真結(jié)果和理論計(jì)算非常接近。
將external emissivity 設(shè)成0.5,計(jì)算出輻射傳熱功率為3114.6W。改變internal emissivity的值,計(jì)算值不變。
從以上仿真結(jié)果可知:
1、2.1小結(jié)的第八點(diǎn)external emissivity的解釋是正確的,輻射傳熱基于灰體假設(shè),輻射系數(shù)等于吸收系數(shù)。
3.2 Case2-測(cè)試internal emissivity-使用DO模型計(jì)算-2D模型
1)Solid(Al)-solid(Steel)-solid(Al)-case
i)Internal solid Fix temperature=400℃, external radiation temperature=300,external emissivtiy=1;internalemissivtiy=1:
圖3-4溫度分布以及換熱量
ii)internal emissivtiy=0:
圖3-5溫度分布以及換熱量
從圖4、5可知,上下兩張圖的溫度分布非常相近,上圖中溫度稍高,而zone之間的換熱量存在差異,將internal emissivity改為0,代表兩個(gè)不同材料的zone之間輻射傳熱量為零,因此總傳熱量從5555W降低至5055W。可知,Fluent中認(rèn)為緊密相連的兩個(gè)solid zone(存在couple wall)之間是存在輻射傳熱的(也可設(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),no gravity-case
i),external emissivtiy=1;internal emissivtiy=1,fluid的absorption coefficient=0;
圖3-6 溫度分布以及換熱量
中心400K的solid往external solid的輻射傳熱功率為:
Pra=5.67e-8*1*3.14*1*(400^4-335^)=2315W,和fluentreprot值2333W(包含了空氣熱傳導(dǎo)的功率)比較接近;
ii)internalemissivtiy=0,fluid的absorption coefficient=0;
圖3-7 溫度分布以及換熱量
將internal emissivtiy=0后,傳熱功率下降為21W,說(shuō)明無(wú)輻射換熱時(shí),僅靠空氣導(dǎo)熱的傳熱功率非常小。
iii)Fluid 和external solid之間的internal emissivtiy=1,fluid的absorption coefficient=1;
圖3-8 溫度分布以及換熱量
iv)Fluid 和external solid之間的internal emissivtiy=0,fluid的absorption coefficient=1;
圖3-9 溫度分布以及換熱量
v)Fluid 和external solid之間的internal emissivity=1,fluid absorption coefficient=0,external solid absorption coefficient=1,conductivity=0.02;
vi)emissivity=1,fluid absorption coefficient=0,external solid absorption coefficient=10,conductivity=0.02
圖3-10 溫度分布以及換熱量
圖3-11 溫度分布以及換熱量
3.3 仿真結(jié)論
從以上仿真結(jié)果,可以得出以下結(jié)論:
1、2.1小結(jié)的internal emissivity以及external emissivity的解釋是正確的。
2、air 的absorption coefficient的默認(rèn)值=0,代表air 不吸收輻射,即不參與輻射。
3、solid的absorption coefficient的默認(rèn)值=0,代表solid吸收輻射,并且absorption coefficient為無(wú)窮大,輻射被固體表面直接吸收。輻射系數(shù)可設(shè)置。
相關(guān)標(biāo)簽搜索:Fluent輻射傳熱模型理論以及相關(guān)設(shè)置 Fluent培訓(xùn) Fluent流體培訓(xùn) Fluent軟件培訓(xùn) fluent技術(shù)教程 fluent在線視頻教程 fluent資料下載 fluent分析理論 fluent化學(xué)反應(yīng) fluent軟件下載 UDF編程代做 Fluent、CFX流體分析 HFSS電磁分析