ANSYS Fluent 在熱分析中的使用介紹 宋步云
1 b0 ^5 t( i2 [) ^, v 1. 基本概念:
% G0 K5 S' O z9 s7 c2 R4 i5 |熱能的傳遞有三種基本的方式:熱傳導,熱對流,熱輻射 , v- @6 b) x Q0 t2 Q P1 _
1.1 熱傳導% M% Q F% s# U( u; G4 ]) l1 y
物體各部分之間不發生相對位移時,依靠分子、原子及自由電子等微觀粒子的熱運動而產生的熱能傳遞稱為熱傳導。導熱的基本定率被總結為傅立葉定率: ; [# @# s' P" K
其中, 為熱流量,單位為 , 為導熱系數,單位為 , 為面積, 為溫度。
( p" `6 ?2 {* R; k4 |2 ^" V8 v一般而言,氣體的導熱系數值約在0.006~0.6 之間,其值隨著溫度的升高而增大。液體的導熱系數約在0.07~0.7 之間,除了水和某些水溶液及甘油外,絕大多數液體的導熱系數會隨著溫度的升高而減小。
1 [, a7 e- J1 Z: k* }1.2 熱對流
5 _3 }; B2 @5 q8 J由于流體的宏觀運動而引起民的流體各部分之間發生相對位移,冷熱流體相互摻混所導致的熱量傳遞過程稱為熱對流。需要說明的是熱對流只能發生在流體當中,而且由于流體中的分子同時在進行著不規則的熱運動,因而熱對流必然伴隨著熱傳導。工程中感興趣的是流體流對一個物體表面時流體與物體表面之間的熱量傳遞過程,我們稱之為對流傳熱,以區別于一般意義上的熱對流。實際上,我們平時所說的熱對流也指這種情況。根據引起流動的原因來劃分,對流傳熱可以區分為自然對流和強制對流兩大類。對流傳熱的基本計算公式為牛頓冷卻公式: 6 W8 Y k. y( v* F% p1 [
其中, 為表面傳熱系數,也被稱為對流換熱系數,單位為 。
) m$ a7 ^: v% K" B2 Z1.3 熱輻射
$ z% q" p5 r+ l" n物體由于熱的原因而發出輻射能的現象我們稱之為熱輻射。理論上講,只要物體的溫度高于絕對零度(0 ),物體就會不斷的把熱能變為輻射能,向外發出熱輻射。熱輻射的基本計算公式為斯忒藩-玻耳茲曼定律,又稱為四次方定律: ![]()
" ?( ?1 P1 L* u7 @/ R其中, 為物體的發射率,也稱為黑度,其值總小于1, 為斯忒藩-玻耳茲曼常量,它是個自然常數,其值為5.67e-08 , 為熱力學溫度,單位 。
+ ]+ K5 Q T1 ? m# t1 T" y以上為三種基本傳熱方式的介紹,在實際問題中,這些方式往往不是單獨出現的,很可能是多種傳熱方式的組合形式。
" ~& v! c. j! j2. 導熱問題的三大類邊界條件7 F% [8 k& D7 A/ k- ^
1) 規定了邊界上的溫度值,稱為第一類邊界條件,也稱為Dirichlet條件。此類條件最簡單的例子就是規定邊界的溫度為常數。
) `0 @9 \/ d$ M* G% U; W* z2) 規定了邊界上的熱流密度值,稱為第二類邊界條件,也稱為Neumann條件。此類條件最簡單的例子就是規定邊界上的熱流密度為常數。8 i6 g s5 K% w' r7 q* g
3) 規定了邊界上物體與周圍流體間的表面傳熱系數及周圍流體的溫度,稱為第三類邊界條件,也稱為Robin條件。
6 B k# [' @ P3 ~; F- j8 K1 B此外,在處理復雜的實際工程問題時,我們還會遇到輻射邊界條件,即物體表面與外界環境之間只發生輻射換熱,如航天器上的發熱元件向太空的散熱。$ U) S+ }. |3 T) L: l
2.1 Fluent中熱邊界的設置 }. w% F# H& j
在Wall設置中的Thermal選項卡中,前三種Heat Flux,Temperature,Convection分別對應前面所說的第二類,第一類及第三類邊界條件。Radiation為輻射邊界條件,Mixed為混合邊界條件。默認情況下壁面為絕熱,即通過壁面的熱流量為0。值得一提的是,當計算當中存在共軛傳熱問題時,導入網格時,Fluent會自動為共軛傳熱交界面生成shadow面,如圖 2所示,一般情況下,此類壁面我們不需要進行額外設置。 2 Q: K2 ], j) x* r9 V
![]() 圖 2 耦合壁面 $ m9 q( _) g$ P0 A
在Heat Flux選項中,需要設置通過壁面的熱流密度及壁面發熱功率(默認壁面不發熱,即發熱功率設為0),如圖 3所示。 3 o% U1 s1 i8 w( x0 D( g; K
![]() 圖 3 第二類邊界條件設置面板 . ?) p' \6 [+ b5 [
在Temperature選項中需要設置壁面的溫度,其它設置與Heat Flux選項一樣,如圖 4所示。 2 s. C5 x* u) S8 U& Z; }! Y
圖 4第一類邊界條件設置面板
# Q! D V) I; d S在Convection選項中,需要設置對流換熱系數及外界溫度,其它設置與Heat Flux選項一樣,如圖 5所示
, g0 O9 A' F4 v8 [' I ![]() 圖 5第三類邊界條件設置面板 . r E6 s" V2 ?) B. H' F
在Radiation選項中需要設置壁面的發射率及外界溫度,其它設置與Heat Flux選項一樣,如圖 6所示。0 ?6 f6 e' p0 f; `/ J: U
* m) X O$ f( ~- E3 n! B, B
![]() 圖 6 輻射邊界條件的設置 8 v! |- d5 }& E2 F7 T% M
在Mixed選項中需要同時指定壁面的表面換熱系數、發射率、外界對流換熱溫度及外界輻射溫度,如圖 7所示
4 E+ x, i$ J0 }6 p" g2 X ![]() 圖 7 混合邊界條件的設置
0 i9 U% q! {+ K3. 不同傳熱方式的計算設置
2 x# Y3 _; J2 U" H" T l; r/ S3.1 熱傳導問題的設置
2 A3 S J4 F2 H/ i; E' H6 x! j在工程計算中,傳熱導的問題通常會以熱阻設置的情況呈現,熱阻的定義為:
# D" f/ Z( p; O' O/ @6 M 為壁面厚度, 為導熱系數, 為面積。, I) j5 }0 O2 `; C' h% ]
在Fluent中,處理這種問題有三種方法:7 n2 Q. J, L: K# d j3 F! B; b4 v+ G
1) 設置一個有厚度的薄壁,并為其劃分網格,設置材料通過求解器來計算熱阻的值。這種方法可以考慮到各個方向上的熱量傳遞過程,但這種做法往往會伴隨的大量的網格數量增加。
1 W; o) Q- F6 E7 t) q. ] ![]() 圖 8 薄壁網格示意圖 / d# ~# m6 |% y1 S/ @! L
2) 為壁面設置一個虛擬厚度。幾何模型中不體現其厚度,因而劃分網格時也僅以0厚度壁面存在,導入Fluent中后,在壁面設置中為其設置一個虛擬的厚度。這種做法可以在考慮熱阻的同時大幅降低網格數量,但這種做法只能考慮到垂直于壁面方向的熱量傳遞過程。在Fluent中,我們可以通過指定材料屬性,壁面厚度來考慮薄壁熱阻對傳熱的影響,如圖 10所示。
, Q$ A+ K6 t5 T/ ]0 R6 {圖 9 虛擬壁面厚度網格示意圖 圖 10 Fluent中設置以設置壁面厚度的方式考慮熱阻
# v- E! @) L! r" ?/ l+ @( Y& a# x& Z3) 設置Shell Conduction,這種方法類似于方法2,但可以考慮到各個方向上的熱量傳遞。同時可以指定多層不同材料的薄壁,如圖 12所示。這種做法在一些工程應用上能夠在保證精度的前提下,大大降低網格劃分的難度及數量。
8 E( S0 o6 \# e# p ![]() 圖 11 Shell conduction網格示意圖
4 ]5 c) ^9 u8 ?" P' @- [ z3 P圖 12 Fluent 中Shell Conduction設置
+ Y* r1 u6 H1 G- V由于方法3的優越性,在很多計算中都會采用這種方式來處理一些薄壁結構。在Fluent17.0及以后的版本當中,我們可以通過Shell Conduction Manager來批量的管理和設置Shell Conduction 。對于大量Shell Conduction的設置,我們還可以通過讀寫csv文件來實現。 ~/ Z# c3 {7 g) s/ g* i; _
![]() 3.2 熱對流問題的設置) }, i+ }1 A- A4 t
3.2.1 強制對流設置
# m( n/ I1 d$ Z6 b/ y+ U在強制對流計算中,一般需要打開湍流模型面板并選取合適的湍流模型。在一般的計算當中,推薦使用Realizable 或SST 模型。
# M# I4 v) `! |' \3.2.2 自然對流設置% M. p3 \2 Z |: K1 D
我們知道,一般情況下,在流體計算中是以雷諾數大小來判斷流動是否為湍流。但在自然對流中,我們不再以雷諾數的大小為判斷依據,取而代之的是瑞利數: V0 m2 _; b, ?, I! a" t6 a7 @
0 c8 b7 v1 J C, a1 x% O0 X# h1 ?![]() 其中, 為膨脹系數, 為重力加速度, 為特征長度, 為溫度差, 為運動粘度, 為熱擴散率。一般認為,當瑞利數大于10e9時,流動為湍流,此時需要打開相應的湍流模型。
4 J [5 |! V' u與一般計算不同,自然對流計算中有一些特別的設置。' D. |% [7 {0 ?9 _3 f, A
1) 由于自然對流是由于重力場下密度的變化所引起的,因此在計算中需要打開重力項,并設置其大小和方向。 " s/ D3 d; S4 h; ~7 v& `
圖 13 打開重力項
: u1 w5 S2 \* t W6 K( M& f2) 數值離散格式。自對流中壓力離散格式需要使用Body Force Weighted或PRESTO!。采用默認的二階格式會出現非物理現象的錯誤結果。
4 A- a" M% E7 @; {# q3 s+ X ![]() 圖 14 選擇離散格式 : L p; f% X& Q& G0 h
3) 參考密度的設置。在自然對流的計算中需要打開重力項并設置重力加速度的方向及大小。勾上Operating Density選項,可以增加計算的穩定性。
9 I$ B% \* ?8 d+ f( C) A% i圖 15 設置參考密度
: Z- c; l4 l4 ~* L3 ]/ ]2 I9 s4) 密度的設置。在自然對流中,由于流體的流動是由于密度變化引起的,因此在材料屬性中需要對相關屬性進行設置。對于氣體而言,密度設置可以選擇Boussinesq假設或不可壓理想氣體模型。對于液體而言,只能選擇Boussinesq假設。對于封閉區域的自然對流計算需要使用Boussinesq假設。Boussinesq模型假設在動量方程中,除了體積力項之外,其它各項的密度為常數。需要注意的是Boussinesq假設只能用在密度變化小于20%的情況下。打開Boussinesq需要在密度設置中選擇boussinesq,并設定參考密度,同時需要設置流體膨脹系數,一般而言,氣體的膨脹系數為其熱力學溫度的倒數。5 B% r2 `# G# @- Z; G
![]() 圖 16 材料屬性設置面板 ; U$ [+ {+ }' J- v* T, o/ T
3.3 熱輻射問題的設置
7 P1 _- E8 B. M; N o在介紹熱輻射計算之前,我們需要了解一下光學厚度(Optical thickness)的概念。光學厚度是介質吸收輻射能力的量度。在Fluent中,光學厚度 ' a5 O! {6 d' }+ s, E; f# ?$ l' v
其中, 為吸收系數,即由于介質吸收而導致的輻射強度在經過每單位長度的介質后改變的量。由于空氣一般不吸收輻射,因此,流體介質為空氣時,該系數可近似設為0。 為散射系數,即由于介質散射而導致的輻射強度在經過每單位長度的介質后改變的量。同樣的,流體介質為空氣時,該系數可近似設置為0。 為特征長度。
$ s n9 N+ V+ S5 {7 n6 G![]() 表格 1 各輻射模型的適用范圍及介紹 模型 光學厚度 計算量8 e' R2 K9 |1 `5 i) I" a- {7 _1 @. i
Surface to surface model (S2S) 0 當光學厚度等于0時,S2S模型的計算精度與DO模型相當,但計算量要小于DO模型
; g5 k5 R" u" F b% _ A9 P! ZRosseland > 5 計算量小,但在實際應用中不多
( W8 \! K9 U5 y, V8 ~9 iP-1 > 1 計算量較小,在光學厚度較大的問題中計算效果較好; y! f0 X. s( P9 q
Discrete ordinates model (DO) All 使用范圍最廣,計算量最大,計算精度最高的輻射模型
9 U( A- {' z' B7 n, gDiscrete Transfer Method (DTRM) All 計算量小,但由于無法用于并行計算,因此很少使用
# N' c4 f# R5 r5 L7 S一般而言,熱輻射模型用在高溫工況及僅有自然對流存在的工況中。若要在計算中考慮熱輻射的影響,需要打開輻射模型面板,在其中選取相應的熱輻射模型。如圖 17所示。
. b; X) o1 s. e8 K4 X. y ![]() 圖 17 輻射模型面板 : B; ^; s2 Q5 W
4. Fluent熱分析關鍵步驟總結* a# O- L S! n' m" O" _: v+ _" }
1) 根據計算問題類型來確定是否需要打開重力項。
. ]1 ]$ ^- ~- a: v/ G3 X, ~0 c ! z/ O; a+ \) Y) k
![]() 圖 18 打開重力項
2 y; X. o& |6 ?& I8 R2) 在Fluent中激活能量方程。
4 f: V" h$ I- l9 F ![]() 圖 19 激活能量方程
: s3 U3 T) ]2 S5 h3) 根據情況選取合適的湍流模型。
; ?6 x$ d7 }5 g* S" Y7 G: z
圖 20 選取湍流模型8 _8 D2 Y# F) _6 I
4) 根據情況選取合適的輻射模型。 ( A$ j/ W5 q; d: I& ?3 C( v0 N
圖 21 選取輻射模型9 A" K+ b) ?# A3 T9 W+ ?& g8 R
5) 設置相關材料屬性,若求解問題為輻射問題需要特別注意設置相關的輻射參數,若為自然對流問題則需要注意密度項的處理方式。
3 n# E" Y5 G0 d! g# F8 q; U6) 設置相關的邊界條件,如流動進出口,壁面的熱邊界等。
. i, e, V+ ^; y' ~$ Q7) 設置合適的離散方式,對于自然對流而言需要特別注意壓力項的離散方式。
6 Q9 N( e2 S1 t% p( b1 X ![]() 圖 22 選取壓力速度耦合方式及各離散項格式
2 k: D3 K) q2 s& A+ a4 n# [8) 初始化 |