完整word版,FLUENT知識點(diǎn)

1、一、基本設(shè)置1. Double Precision 的選擇啟動設(shè)置如圖，這里著重說說Double Precision （雙精度）復(fù)選框，對于大多數(shù)情況，單精度求解器已能很好的滿足精度要求，且計算量小，這里我們選擇單精度。然而對于以下一些特定的問題，使用雙精度求解器可能更有利1。a. 幾何特征包含某些極端的尺度（如非常長且窄的管道），單精度求解器可能不能足夠精確地表達(dá)各尺度方向的節(jié)點(diǎn)信息。b. 如果幾何模型包含多個通過小直徑管道相互連接的體，而某一個區(qū)域的壓力特別大（因?yàn)橛脩糁荒茉O(shè)定一個總體的參考壓力位置），此時，雙精度求解器可能更能體現(xiàn)壓差帶來的流動（如 漸縮漸擴(kuò)管的無粘與可壓縮流動模擬）。c

2、. 對于某些高導(dǎo)熱系數(shù)比或高寬縱比的網(wǎng)格，使用單精度求解器可能會遇到收斂性不佳或精確度不足不足的問題，此時，使用雙精度求解器可能會有所幫助。FLUENT LauncherOptasDouble P redsianPracsssrg OipHiaris Serial -ParallelANSYSDimersioniod2D"r- II-"If I 3DDisplay Options叼 DispIdA MeshAfier Pleading 団 Embed Graphics Windows0 V/orkbench Color S chenrie圈 Show More Opt io

3、 ns0K1DeiaJtCancelf Help 生2. 網(wǎng)格光順化用光滑和交換的方式改善網(wǎng)格：通過Mesh下的Smooth/Swap來實(shí)現(xiàn)，可用來提高網(wǎng)格質(zhì)量，一般用于三角形或四邊形網(wǎng)格，不過質(zhì)量提高的效果一般般，影響較小，網(wǎng)格質(zhì)量的 提高主要還是在網(wǎng)格生成軟件里面實(shí)現(xiàn)，所以這里不再用光滑和交換的方式改善網(wǎng)格，其原理可參考FLUENT全攻略（已下載）。3. Pressure-based 與 Density-based求解器設(shè)置如圖。下面說一說Pressure-based和Density-based的區(qū)別：P ressure-Based Solver是Flue nt的優(yōu)勢，它是基于壓力法的求解

4、器，使用的是壓力修正 算法，求解的控制方程是標(biāo)量形式的，擅長求解不可壓縮流動，對于可壓流動也可以求解；Flue nt 6.3以前的版本求解器，只有Segregated Solver和 Cou pled Solver，其實(shí)也是P ressure-Based Solver的兩種處理方法；Den sity-Based Solver是Flue nt 6.3新發(fā)展出來的，它是基于密度法的求解器，求解的控 制方程是矢量形式的，主要離散格式有Roe, AUSM+，該方法的初衷是讓Flue nt具有比較1 李鵬飛，徐敏義，王飛飛精通CFD工程仿真與案例實(shí)戰(zhàn)：FLUENT GAMBIT ICEM CFD Tec

5、plotM.北京,人民郵電出版社，2011:114-116好的求解可壓縮流動能力，但目前格式?jīng)]有添加任何限制器，因此還不太完善；它只有Coupled的算法；對于低速問題，他們是使用Preconditioning方法來處理，使之也能夠計算低速問題。Density-Based Solver下肯定是沒有 SIMPLEC，PISO這些選項(xiàng)的，因?yàn)檫@些都是 壓力修正算法，不會在這種類型的求解器中出現(xiàn)的；一般還是使用P ressure-Based Solver解決問題?；趬毫Φ那蠼馄鬟m用于求解不可壓縮和中等程度的可壓縮流體的流動問題。而基于密度的求解器最初用于高速可壓縮流動問題的求解。雖然目前兩種求解器

6、都適用于各類流動問題的求解（從不可壓縮流動到高度可壓縮流動），但對于高速可壓縮流動而言，使用基于密度的求解器通常能獲得比基于壓力的求解器更為精確的結(jié)果。GeneralMesh Check Report Quality Display" SolverTypeI* Ftessure-eased 廣 Density-Based呻odtv Formulation 唾 Absolute廣 RelatjveTirne11 Steady 'C Transient四5p電e"4 PlrtarC Axisyrnmetrk匚 Axismmetric 5wirtunits. Gravit

7、y4. axisymmetric 禾R axisymmetric swirl從字面的意思很好理解axisymmetric和axisymmetric swirl的差別：axisymmetric :是軸對稱的意思，也就是關(guān)于一個坐標(biāo)軸對稱，2D的axisymmetric問題仍為2D問題。而axisymmetric swirl :是軸對稱旋轉(zhuǎn)的意思，就是一個區(qū)域關(guān)于一條坐標(biāo)軸 回轉(zhuǎn)所產(chǎn)生的區(qū)域，這產(chǎn)生的將是一個回轉(zhuǎn)體，是3D的問題。在Flue nt中使用這個，是將一個3D的問題簡化為2D問題，減少計算量，需要注意的是，在Flue nt中，回轉(zhuǎn)軸必須是x軸。5.操作工況參數(shù)（Op erati ng C

8、on ditio ns ） 操作壓力的介紹關(guān)于參考壓力的設(shè)定，首先需了解有關(guān)壓力的一些定義。ANSYS FLUENT中有以下幾個壓力，即 Static Pressure（靜壓）、Dynamic Pressure （動壓）與 Total Pressure（總壓）；Absolute P ressure （絕對壓力）、Relative P ressure （參考壓力）與 Op erat ing P ressure （操作壓力）。這些壓力間的關(guān)系為，Total Pressure （總壓）=Static Pressure （靜壓）+Dynamic Pressure（動壓）；Absolute Pressu

9、re （絕對壓力）=Operating Pressure （操作壓力）+Gauge Pressure （表 壓）。靜壓是測量到的壓力，動壓是其中，靜壓、動壓和總壓是流體力學(xué)中關(guān)于壓力的概念。 有關(guān)速度動能的壓力，是流動速度能量的體現(xiàn)。而絕對壓力、操作壓力和表壓是FLUENT引入的壓力參考量，在 ANSYS FLUENT 中,所有設(shè)定的壓力都默認(rèn)為表壓。這是考慮到計算精度的問題。 操作壓力的設(shè)定設(shè)定操作壓力時需要注意的事項(xiàng)如下：對于不可壓縮理想氣體的流動，操作壓力的設(shè)定直接影響流體密度的計算，因?yàn)閷τ诶硐霘怏w而言，流動的密度由理想氣體方程獲得，理想氣體方程中的壓力為操作壓力。對于低馬赫數(shù)的可壓縮

10、流動而言，相比絕對靜壓，總壓降是很小的，因此其計算精度很容易受到數(shù)值截斷誤差的影響。需要采取措施來避免此誤差的形成，ANSYSFLUENT通過采用表壓（由絕對壓力減去操作壓力）的形式來避免截斷誤差的形 成，操作壓力一般等于流場中的平均總壓。對于高馬赫數(shù)可壓縮流動的求解而言，因?yàn)榇藭r的壓力比低馬赫可壓縮流動的大得多，所以求解過程中的截斷誤差的影響不大，可以不設(shè)定表壓。由于ANSYSFLUENT中所有需輸入的壓力都為表壓，因此此時可以將操作壓力設(shè)定為0（這樣可以最小化由于壓力脈動而引起的誤差），使表壓與絕對壓力相等。如果密度設(shè)定為常數(shù)或者其值由通過溫度變化的函數(shù)獲得，操作壓力并沒有在計算密度的過程

11、中被使用。默認(rèn)的操作壓力為101325Pa。操作壓力的設(shè)定主要基于兩點(diǎn)考慮，一是流動馬赫數(shù)的大小，二是密度計算方法。表格1操作壓力的推薦設(shè)置密度關(guān)系式馬赫數(shù)操作壓力理想氣體定律大于0.10或約等于流場的平均壓力理想氣體定律小于0.1約等于流場的平均壓力關(guān)于溫度的函數(shù)不可壓縮不使用常數(shù)不可壓縮不使用不可壓縮的理想氣體不可壓縮約等于流場的平均壓力關(guān)于參考壓力位置的設(shè)定對于不涉及任何壓力邊界條件的不可壓縮流動，ANSYS FLUENT在每次迭代后要調(diào)整表壓值。這個過程通過使用參考壓力位置處（或該位置附近）節(jié)點(diǎn)的壓力完成。因此，參考 壓力位置處的表壓應(yīng)一直為0。如果使用了壓力邊界條件，則不會使用到上述

12、關(guān)系，因此參考壓力位置不被使用。參考壓力位置默認(rèn)為等于或接近（0，0，0）的節(jié)點(diǎn)中心位置。實(shí)際計算中可能需要設(shè)置參考壓力位置到絕對靜壓已知的位置處。在Operating Conditions對話框中的 RefereneePressure Location選項(xiàng)組中設(shè)置新的參考壓力位置的x， y， z的坐標(biāo)即可。如果要考慮某一方向的加速度，如重力，可以勾選Gravity復(fù)選框。對于 VOF 計算，應(yīng)當(dāng)選擇 Sp eeified Op erat ing Den sity，并且在 Op erat ing Den sity 下為 最輕相設(shè)置密度。這樣做排除了水力靜壓的積累，提高了round-off精度為

13、動量平衡。同樣需要打開Imp licit Body Force，部分平衡壓力梯度和動量方程中體積力，提高解的收斂性。（空Referenee Pressure Location （參考壓強(qiáng)位置）應(yīng)是位于流體永遠(yuǎn)是100%的某一相氣）的區(qū)域，光滑和快速收斂是其基本條件。Op erating ConditionsGravityOperdtng Pressure (jDascal)1013250 GravityGravitational Accel era ton Crr/s2)-9,S|P|Reference Pressure Locatjori回PIV Cm /s2)BoLissiiesq Par

14、ametersOperdtrig Temperature Ic)Pressure2B3J6Variable-Dena ty Par ametersI Specified Operating DenatyOKCancel Help二、求解模型的設(shè)定1. 流動模型的設(shè)置無粘模型即真正所以理想流體是一種設(shè)想的沒有粘性的流體， 在流動時各層之間沒有相互作用的切應(yīng)力， 沒有內(nèi)摩擦力。十分明顯，理想流體對于切向變形沒有任何抗拒能力。應(yīng)該強(qiáng)調(diào)指出，的理想流體在客觀實(shí)際中是不存在的，它只是實(shí)際流體在某些條件下的一種近似模型。在In viscid流動模型應(yīng)用方面，無粘流動忽略了粘性對流動的影響，這對高雷諾數(shù)的流

15、 動是合適的，因?yàn)楦呃字Z數(shù)流動慣性力的作用遠(yuǎn)大于粘性力的作用，粘性力可以忽略，可以將其考慮成無粘流動。無粘流動的求解更快，其激波在某些值上預(yù)測的偏高。無粘流動能對流動狀態(tài)和激波位置進(jìn)行快速預(yù)測。馬赫數(shù)與激波 馬赫數(shù)的定義是它表示流體的流動速度與當(dāng)?shù)芈曀僦龋?和M 1這三種情況的流動分別稱為亞聲速流、 的相對變化只能引起很小的密度相對變化， 化，這也說明了馬赫數(shù)是流體壓縮性的一個表征。a對應(yīng)于M是一個無量綱的參量。對應(yīng)于M 1，M 1 聲速流和超聲速流。 當(dāng)馬赫數(shù)很小時，速度 但當(dāng)馬赫數(shù)很大時，則將引起較大的密度相對變當(dāng)飛機(jī)、炮彈和火箭以超音速飛行時，或者發(fā)生強(qiáng)爆炸、強(qiáng)爆震時，氣流受到急劇的壓

16、縮，壓強(qiáng)和密度突然顯著增加，這時所產(chǎn)生的壓強(qiáng)擾度將比聲速大得多的速度傳播，波陣面所到之處氣流的各種參數(shù)都將發(fā)生顯著變化，參數(shù)突躍。這樣一個強(qiáng)間斷面叫做激波陣面。漸縮漸擴(kuò)管的流動是計算流體力學(xué)模擬的經(jīng)典問題之一。在這類流動中，激波的出現(xiàn)是流動中可壓縮效應(yīng)的體現(xiàn)。精確的激波模擬是CFD研究的熱點(diǎn)之一。為了更好捕捉壓力梯度，需要采用較細(xì)的網(wǎng)格并結(jié)合合適的數(shù)值模擬和格式。很多實(shí)際模擬中，局部網(wǎng)格的自適應(yīng)會很有幫助。層流模型流動有層流和湍流之分，判斷湍流的標(biāo)準(zhǔn)可以參考2，這里寫出內(nèi)流的判斷標(biāo)準(zhǔn):Re2300對于內(nèi)流而言，一般大多數(shù)流動都是湍流， 一般不使用湍流模型。而對一些外流而言（如 外掠平板或是外掠

17、障礙物），則很有可能是層流運(yùn)動。 湍流模型的評價與選擇a. k 湍流模型這里我們使用的湍流模型是 Standard k模型，這種模型應(yīng)用較多，計算量適中，有較多數(shù)據(jù)積累和比較高的精度，對于曲率較大和壓力梯度較強(qiáng)等復(fù)雜流動模擬效果欠佳。一般工程計算都使用該模型，其收斂性和計算精度能滿足一般的工程計算要求，但模擬旋流和繞流時有缺陷。壁面函數(shù)的選擇對于有壁面的流動，當(dāng)主流為充分發(fā)展湍流時，根據(jù)離壁面法線距離不同，可將流動劃 分為壁面區(qū)（或稱內(nèi)區(qū)、近壁區(qū)）和核心區(qū)（或稱外區(qū)）核心區(qū)是完全湍流區(qū)，為充分發(fā)展的湍流。在動量、熱量和質(zhì)量的交換過程中粘性力起主要作 因此流動幾乎可以看成層流流動，且在平行于壁面

18、方向上的速在壁面區(qū)，由于有壁面的影響， 流動與核心區(qū)不同。 壁面區(qū)可分為3個子層：粘性底層、 過渡層和對數(shù)率層。粘性底層是一個緊貼壁面的極薄層， 用，而湍流切應(yīng)力可以忽略， 度呈線性分布。在此層中，粘性力和湍流切應(yīng)力的作用相當(dāng)， 流動狀況較為 實(shí)際工程計算中由于過渡層厚度極小， 可不考慮此層，直接過渡層處于粘性底層之外， 復(fù)雜，很難用公式或定律表述。 以對數(shù)率層的方法處理。對數(shù)率層處于近壁區(qū)的最外層，粘性力的影響不明顯， 湍流切應(yīng)力占主要地位， 流動處于充分發(fā)展的湍流狀態(tài)，流速分布接近對數(shù)律。壁面區(qū)內(nèi)不同子層的高度和速度可以沿壁面法向的無量綱高度和無量綱速度表達(dá)。2 李鵬飛，徐敏義，王飛飛.精

19、通CFD工程仿真與案例實(shí)戰(zhàn)：FLUENT GAMBIT ICEM CFD TecplotM. 北京, 人民郵電出版社，2011:122其中，U是流體的時均速度，U是壁面摩擦速度，UW是壁面切應(yīng)力，y是壁面的垂直距離。在605時，區(qū)域?yàn)檎承缘讓?，此時速度沿壁面法線方向呈線性分布，即y 300時，流動處于對數(shù)率層，此時速度沿壁面法線方向呈對數(shù)率分布，即u 2.5lny 5.5。使用標(biāo)準(zhǔn)的或者非平衡的壁面函數(shù)（Re 106）。O標(biāo)準(zhǔn)壁面函對于大多數(shù)高雷諾數(shù)情況壁面函數(shù)法的本質(zhì)是，對于湍流核心區(qū)的流動使用k模型求解，而在壁面區(qū)并不進(jìn)行求解，直接使用半經(jīng)驗(yàn)公式得出該區(qū)域的速度等物理量。FLUENT提供了

20、多種壁面函數(shù)處理方式，如標(biāo)準(zhǔn)壁面函數(shù)法、非平衡壁面函數(shù)法和增 強(qiáng)壁面處理。標(biāo)準(zhǔn)壁面函數(shù)法利用對數(shù)校正法提供了必需的壁面邊界條件（對于平衡湍流邊界層） 而非平衡壁面函數(shù)法用來改善高壓力梯度、分離、再附和滯止等情況下的結(jié)果。數(shù)法和非平衡壁面函數(shù)法都允許在近避免區(qū)域上使用較粗的網(wǎng)格。優(yōu)點(diǎn)缺點(diǎn)標(biāo)準(zhǔn)壁面函數(shù) 法應(yīng)用較多，計算量小，有較 高的精度適合高雷諾數(shù)流動， 對低雷諾數(shù)流動問題， 有 壓力梯度、高度蒸騰和大的體積力、低雷諾數(shù) 和高速三維流動問題不適合非平衡壁面函數(shù)法考慮了壓力梯度，可以計算 分離，在附著以及撞擊問題對低雷諾數(shù)流動問題，有較強(qiáng)壓力梯度、強(qiáng)體 積力及強(qiáng)三維性問題不適合增強(qiáng)壁面處理不依賴壁

21、面法則，對于復(fù)雜 流動，特別是低雷諾數(shù)流動 很適合要求網(wǎng)格密，因而要求計算機(jī)處理時間長，內(nèi)存大增強(qiáng)壁面處理選項(xiàng)把混合邊界模型和兩層邊界模型結(jié)合起來, 近壁面現(xiàn)象很適合，湍流模型在內(nèi)層上得到了修正。表格2幾種壁面處理方法的比較對低雷諾數(shù)流動或者復(fù)雜Q VbcoiJ5 ModdMmJsIModel C«nstaritsI _ fl riVKod(j LaminarII Sp alart Allmar as. (1 eqr 'O«k*ep£i1on (2 eqn)匚«t-o*necia (3eqn)I _ * TransfLon k-kJ -omega

22、 (3 eqnj»_ * TransfLori SST 斗 eqn'JI Revriolds Stress (5 ew)S仙也pB低 Simd目ton (SAS)k<epsilon Model 佔(zhàn)"standard (_-I FtNQI * RealizableNeerWall Tretmert崎n Standard Wall FiincliOriS(_ “ ScaleWe Wall FunctionsIJ Non-Equilibrium Wai Fijnctioris 廠 * Enhanced Wall Treatjrwnt 廠 N User-Defined

23、 Wall FupictonsOptionsZU WiscoLB HeatingZlCurvature Correctionc0.09Ci-£pdon1.44CZ-EpsibnTKE FraudtJ NumberInoneTPtardtJ NumbersTXE Prandtl Nuntjernone>TDR. Prand-H Numbernor>e<&ierqy Prandtl NumbernoneUser 'Pefined FunctionsTurbulent ViscoartyHfclp 2. 多相流模型 VOF模型該模型通過求解單獨(dú)的動量方程和

24、處理穿過區(qū)域的每一流體的容積比來模擬兩種或三 種不能混合的流體。典型的應(yīng)用包括流體噴射、流體中大泡運(yùn)動、流體在大壩壩口的流動、 氣液界面的穩(wěn)態(tài)和瞬態(tài)處理等。一般而言VOF主要適用于非穩(wěn)態(tài)的多相流模型，僅對某些特定問題的多相流模型的穩(wěn)態(tài)問題能夠適用。VOF方法適用于計算空氣和水這樣不能互相摻混的流體流動，對于分層流和活塞流，最方便的就是選擇 VOF模型。需要注意的是，對于湍流模型的設(shè)置，VOF不能用于無粘流， 也不能用大渦模擬3。Geo-Rec on struct 格式Geo-Reconstruct格式（在Solution Methods中設(shè)置）是一種較為精確的追蹤自由表面的 計算格式，廣泛地應(yīng)

25、用于瞬變流的VOF問題中，但必須注意的要使用該格式VOF模型必須使用顯示離散格式（在 VOF模型設(shè)置選項(xiàng)設(shè)置）。Body Force Formulati on為提高解的收斂性，對于涉及到表面張力的計算，建議在Body Force Formulation中勾選Imp licit Body Force。這樣做由于壓力梯度和動量方程中表面張力的部分平衡，從而提高 解的收斂性。3李進(jìn)良，李承曦，胡仁喜.精通FLUENT63流場分析M.北京，化學(xué)工業(yè)出版社,2009:231-236ModelIVIukiphagehkmbef of Etilefian PhasesO Off tTJ Volurrw of

26、 HurdO MixtureO Ejlerian'.VetStea；v Level SetSctierre.fi'iiEiphdtJlmpbot open diaonel Flow ll Open Channel Wave BC n Zonal DiscretijstQriVftlunw Fration Paramt也肱OptiopieVolume Fracton Cutoff |le-06匚ffljrant hhjmbefBody Force Formulatiori onpiicit Bod* ForceHelp Mixture模型這是一種簡化的多相流模型，用于模擬各種有不

27、同速度的多相流，但是假定了在短空間尺度上局部的平衡。相之間的耦合應(yīng)當(dāng)是很強(qiáng)的。它也用于模擬有強(qiáng)烈耦合的各向同性多相 流和各向以相同速度運(yùn)動的多相流。典型的應(yīng)用包括沉降(sedimentation)、氣旋分離器、低載荷作業(yè)下的多粒子流動、氣相容積率很低的泡狀流。Mixture P arameters一般需要勾選 Mixture Parameters中的Slip Velocity復(fù)選框，以此來求解滑移速度模型, 因?yàn)樵诙嘞嗔髦懈鞣N組分的速度有很大不同。對于求解一個均勻的多相流問題可以選擇不做滑移速度的計算，可以在 mixture parameters選項(xiàng)下將slip velocity關(guān)掉。Mfid

28、eJMulti ph 曰旨 e ModelNumber of Eulerian Phasesg OffQ Votume of Flud (i/MiitueCl Eulerian內(nèi)討4妒刁mMixture Psrarreters0 Slifj VeltHzitvBody Faroe Porrmjiation0 Impliot Body ForceOKCarioel Help Eulerian模型該模型可以模擬多相分離流及相互作用的相，相可以是液體、 氣體、固體。與在離散相模型中Eulerian-Lagrangian方案只用于離散相不同，在多相流模型中Eulerian方案用于模型中的每一項(xiàng)。3.

29、固化與熔化模型FLUENT 采用“焓多孔度（enthalpy-porosity ）”技術(shù)模擬流體的固化和熔化（Solidificatio n/Melt ing ）過程。在流體的固化和熔化問題中，流場可以分成流體區(qū)域、固體 區(qū)域和兩者之間的糊狀區(qū)域?！办室欢嗫锥取奔夹g(shù)采用的計算策略是將流體在網(wǎng)格單元內(nèi)占有的體積百分比定義為多孔度（P orosity），并將流體和固體并存的糊狀區(qū)域看作多孔介質(zhì)區(qū) 進(jìn)行處理。在流體的固化過程中，多孔度從1降低到0;反之，在熔化過程中，多孔度則從0升至1。“焓一多孔度”技術(shù)通過在動量方程中添加匯項(xiàng)（即負(fù)的源項(xiàng)）模擬因固體材料 存在而出現(xiàn)的壓強(qiáng)降。“焓多孔度”技術(shù)可以模

30、擬的問題包括純金屬或二元合金中的固化、熔化問題、連續(xù)鑄造加工過程等。計算中可以計算固體材料與壁面之間因空氣的存在而產(chǎn)生的熱阻，固化、 熔化過程中組元的輸運(yùn)等等。需要注意的是，在求解固化、熔化問題的過程中，只能采用分離算法，只能與 VOF模型配合使用，不能計算可壓縮流，不能單獨(dú)設(shè)定固體材料和流體材 料的性質(zhì)，同時在模擬帶反應(yīng)的組元輸運(yùn)過程時，無法將反應(yīng)區(qū)限制在流體區(qū)域，而是在全流場進(jìn)行反應(yīng)計算。 Parameters 定義在Parameters下面定義 Mushy Zone Constant （糊狀區(qū)域常數(shù)）。這個常數(shù)的取值范圍一般在104到107之間，取值越大沉降曲線就越陡峭，固化過程的計算速

31、度就越快，但是取 值過大容易引起計算振蕩，因此需要在計算中通過試算獲得最佳數(shù)值。 Materials設(shè)置在 Materials （材料）面板上，定義 Melting Heat （熔化熱）、Solidus Temperature （固相 點(diǎn)溫度）和Liquidus Temperature （液相點(diǎn)溫度）。如果計算中涉及組元輸運(yùn)過程，則必須同 時定義溶劑的融解溫度（Melting Temperature ）,同時需要定義熔化物的液相線相對于濃度的斜率(Slope of Liquidus Line )、分配系數(shù)(Partition Coefficient )和固體中的擴(kuò)散速率 (Diffusion

32、in Solid、等參數(shù)。設(shè)置邊界條件除了常規(guī)的邊界條件設(shè)置，對于固化和熔化問題還有一些特殊設(shè)置，其中包括：在計算壁面接觸熱阻時設(shè)置接觸熱阻(Con tact Resista nee)。這個參數(shù)在Wall (壁面、面板中的Thermal Conditions (熱力學(xué)條件、下給定。Shear Condition (剪切條件)下選擇如果需要定義壁面上表面張力對溫度的梯度，則在 Marangoni Stress ( Marangoni 應(yīng)力)選項(xiàng)。如果計算拉出速度，則在邊界條件中的速度邊界條件將被用于拉出速度的計算。Iidiflegation and MeltingMcxlel0 Eclid if

33、i cstidTeb ngP 申宙 mehes陽加心Sr”以r咖訕100000Bd(± DiffusionOlndtide Pull VeJodtiesOKCancel Help三、相設(shè)置相設(shè)置一般用于多相流的設(shè)置，對于相設(shè)置，這里主要講一下In teraction的設(shè)置，如圖:PbasesPhaseswater - P nrnay P hase air - Seccmdary Phasecdif-.rritftractan.14Interaction 設(shè)置Interaction設(shè)置用來定義兩相的相互作用，其有多個選項(xiàng)卡，如圖。Pfidie Interaction"bM I

34、 MML I嚇h Sj岳vh那 I閥 Gjg 曲*茸 II jH I 冷出站也口v| Ti/tkETlD町t4r 1*jbU«rlrL«rKtDr.| C口附3讓 2dJay"論LHso*Cnori Hff-pDrag選項(xiàng)卡針對每對物相，在下拉菜單中選擇阻力函數(shù)。其中包括schiller-naumann模型、morsi-alexander模型、symmetric （對稱）模型等用于流體與流體之間阻力計算的模型，也 包括wen-yu模型、gidaspow模型、syamlal-obrien模型等用于液體與固體之間阻力計算的 模型，還包括 syamlal-obrien-

35、symmetric模型用于固體與固體之間的阻力計算。除此之外， 還可以將阻力函數(shù)定義為constant （常數(shù)），或者選擇user-defined （用戶定義）由用戶自己定義阻力函數(shù)。如果計算中不需要設(shè)定阻力，還可以選擇none （不計阻力）選項(xiàng)。阻力設(shè)置的相關(guān)原理比較復(fù)雜，可參考幫助，一般保持默認(rèn)的schiller-naumann設(shè)置不變。Surface Tension 選項(xiàng)卡Surface Tension選項(xiàng)卡用來定義表面張力，如果相包含壁面粘附， 可勾選“ Wall Adhesion ”復(fù)選框。四、Cell Zone ConditionFrame Motion 選項(xiàng)對于流體，可以通過 F

36、rame Motion選項(xiàng)確定坐標(biāo)運(yùn)動方式(如離心泵內(nèi)部流體的旋轉(zhuǎn) 使用運(yùn)動參考系模型)，如圖：Flluiid7one NameHjid. 7Mstsrial NameairEdit 巨 rraiTkC Mo bon Source T cftra匚 Mesh Motion Fixed VaJurtI IPorajs ZoneLDFReference Fra |iMesh Motior | Paroue Zone | Embedded LES | Rtarton | SoifCft Terms | Fixed '/alues | Mulhiihast |Relative Speofiot

37、iorRelflthe To Cd Tore absokjtsZone Motior FunctionnoneRolatiorAxE 0廣pinX fni任EonctintoonstantRolatiqral Vcloci：¥Sfcd (ra/s)orstanlCopy 10 Mesh Motwnirransktionfll VeodtvconstantconstantOK I I Lanoel HelpRotation-Axis Origin 指轉(zhuǎn)軸，Rotational Velocity 值旋轉(zhuǎn)速度。五、邊界條件設(shè)置1.速度進(jìn)口邊界條件(Velocity Inlet )速度進(jìn)口邊

38、界條件用進(jìn)口處流場速度及相關(guān)流動變量作為邊界條件，在速度進(jìn)口邊界條件中，流場進(jìn)口邊界的駐點(diǎn)參數(shù)是不固定的。 為了滿足進(jìn)口處的速度條件， 駐點(diǎn)參數(shù)將在一 定范圍內(nèi)波動。需要注意的是，因?yàn)檫@種條件中允許駐點(diǎn)參數(shù)浮動，所以速度進(jìn)口邊界條件僅適用于不不要讓速度進(jìn)口可壓流，如果用于可壓流，則可能導(dǎo)致出現(xiàn)非物理解。同時還要注意的是， 條件過于靠近進(jìn)口內(nèi)側(cè)的固體障礙物，這樣會使駐點(diǎn)參數(shù)的不均勻程度大大增加。湍流參數(shù)的設(shè)置對于一般的流動邊界條件，均涉及到湍流參數(shù)的定義，在Turbule nee Sp ecificatio n Method(湍流定義方法)下拉列表中，可以簡單地用一個常數(shù)來定義湍流參數(shù)，即通過給

39、定湍流強(qiáng)度、湍流粘度比、水力直徑或湍流特征長在邊界上的值來定義流場邊界上的湍流。a. 湍流強(qiáng)度(Turbulenee Intensity )湍流強(qiáng)度定義如下：|Ju2 V2 w2上式中U , v，w是速度脈動量，u是平均速度。比較常用的是Intensity and Hydraulic Diameter，湍流強(qiáng)度與水力直徑的確定有相應(yīng)的計 算方法，這里只是采用估算來加以確定。內(nèi)流問題進(jìn)口處的湍流強(qiáng)度取決于上游流動狀態(tài)。動，則進(jìn)口處可以使用低湍流強(qiáng)度。如果上游是充分發(fā)展的湍流，到幾個百分點(diǎn)。如果管道中的流動是充分發(fā)展的湍流，則湍流強(qiáng)度 可以用下面公式計算得到，這個公式是從管流經(jīng)驗(yàn)公式得到的：如果上

40、游是沒有充分發(fā)展的未受擾流 則進(jìn)口處湍流強(qiáng)度可以達(dá)I ( turbulenee intensity)I ± 0.16(ReDH) 1/8 u其中：ReH為按等效水力直徑Dh計算得到的雷諾數(shù)。湍流強(qiáng)度小于1%時，可以認(rèn)為湍流強(qiáng)度是比較低的，而在湍流強(qiáng)度大于10%時，則可以認(rèn)為湍流強(qiáng)度是比較高的。比如，當(dāng)雷諾數(shù)為50000時，代入上述公式可得到湍流強(qiáng)度為4%，默認(rèn)的湍流強(qiáng)度設(shè)置為 5% (中等強(qiáng)度)。勺In let7ore Marnevelocity-inlet-5Velocity Specifiedtion MetfiodMagntude, Moriral to BoundaryRef

41、ierence FrarneVelodtv Magnitude On/s)0.2constantni tial Gauge Pressure pascal)lOconstantSxafiaton MethodSupersonic/lMoirertuni | iherTal Radiation | Speoes | DPM | P1ultiDhase| UDSTurbulenceIntiensitv jnd Hydraulic DianneberTurfciJent Ihhensity (鳴)Hydraulic Diameter &n)32OK Cmcgl Help Ib.湍流的長度尺度

42、與水力直徑湍流的長度尺度（Turbulenee Length Scale）與水力直徑（Hydraulic Diameter ）是設(shè)置湍 流的重要參數(shù)。l則是與大渦結(jié)構(gòu)相關(guān)的物理量。在 所以l是受到管道尺寸制約的幾何。湍流能量主要集中在大渦結(jié)構(gòu)中，而湍流長度尺度 充分發(fā)展的管流中，因?yàn)殇鰷u尺度不可能大于管道直徑， 湍流長度尺度l與管道物理尺寸L關(guān)系可以表示為：l 0.07L式中的比例因子0.07是充分發(fā)展管流中混合長的最大值，而L則是管道直徑。在管道截面不是圓形時，L可以取為管道的水力直徑。水力直徑是在管內(nèi)流動（internal pipe flow ）中引入的，其目的是為了給非圓管流動取一個合適

43、的特征長度來計算其雷諾數(shù)。四倍的濕橫截面面積與濕圓周長度之商。管道直比在這湍流的特征長取決于對湍流發(fā)展具有決定性影響的幾何尺度。在上面的討論中，徑是決定湍流發(fā)展過程的唯一長度量。如果在流動中還存在其他對流動影響更大的物體，如在管道中存在一個障礙物，而障礙物對湍流的發(fā)生和發(fā)展過程起著重要的干擾作用。種情況下，湍流特征長就應(yīng)該取為障礙物的特征長度。從上面的分析可知，雖然上式對于大多數(shù)管道流動是適用的，但并不是普遍適用的，某些情況下可以進(jìn)行調(diào)整。在FLUENT中選擇特征長L或湍流長度尺度l的方法如下：1）對于充分發(fā)展的內(nèi)流，可以用Inten sity and Hydraulic Diameter （

44、湍流強(qiáng)度與水力直徑）方法定義湍流，其中湍流特征長度就是Hydraulic Diameter （水力直徑）DH。2） 對于導(dǎo)向葉片或分流板下游的流場，可以用Intensity and Hydraulic Diameter （湍流強(qiáng)度與水力直徑）定義湍流，并在Hydrauli Diameter （水力直徑）中將導(dǎo)向葉片或分流板的開口部分的長度L定義為特征長度。Intensity and3）如果進(jìn)口處的流動為受到壁面限制且?guī)в型牧鬟吔鐚拥牧鲃?，可以在Length Scale面板中用邊界層厚度99通過公式l 0.4 99計算得到湍流長度尺度l。最后在Turbulenee Length Scale （湍

45、流長度尺度）中輸入 l的值。c.湍流動能和湍流耗散率在使用各種k模型對湍流進(jìn)行計算時，需要給定進(jìn)口邊界上的湍流動能k（ TurbulentKinetic Energy ）和湍流耗散率（Turbulent Dissipation Rate ）的估算值。目前沒有理論上的精確計算這兩個參數(shù)的公式，只能通過試驗(yàn)得到。但不可能對各種各樣的流動都去做試驗(yàn)，因此，我們必須借助文獻(xiàn)中已有的近似公式來估算（許多商用CFD軟件也是這樣處理的）。對于沒有已知條件的情況， 可根據(jù)湍流強(qiáng)度I和特征長度L,由下式粗略估計k和 的分布：32k 2（Uavgl）I 3/203/4 l式中C為湍流模型中的一個經(jīng)驗(yàn)常數(shù)，其值約等

46、于0.09。在沒有直接輸入湍流動能k和湍流耗散率的情況下，可以用Intensity and HydraulicDiameter （湍流強(qiáng)度與水力直徑）或 Intensity and Length Scale （湍流強(qiáng)度與長度尺度）等辦 法，利用上述公式確定湍流動能k和湍流耗散率。2.壓強(qiáng)進(jìn)口邊界條件（ P ressure In let ）壓強(qiáng)進(jìn)口邊界條件用于定義進(jìn)口流體的壓強(qiáng)，可用于不可壓縮和可壓縮流動。當(dāng)進(jìn)口壓強(qiáng)已知，而流動速度或流量未知時，可使用壓強(qiáng)進(jìn)口邊界條件。壓強(qiáng)進(jìn)口邊界條件也可用于定義外部或非受限流動的“自由邊界”。pres sure Inlet£ori& Name

47、inletM&mertufn | Tjia-mal Radiation | Spews | DPM | MultipHase | UDS |Specification, Method intensity 占皿 Hydraulic DiameterTurbuJent Inbensaty (%)Hydraulic Diametfir (m)OKCancelH«lp '定義總壓與靜壓首先道總壓（1PPo 72（動量）選項(xiàng)卡內(nèi)，Refere nee Frame為參考值，有絕對值（ Absolute ）與Relative to Adjaeent Zone ）兩個選項(xiàng)供選擇，一般

48、保持默認(rèn)的絕對值，V2P ）與靜壓（po ）的關(guān)系如下（根據(jù)伯努力積分）在 Momentum相對于臨近區(qū)域值（Gauge Total Pressure （表總壓）文本框中輸入總壓的值。靜壓在FLUENT中被稱為 Supersonic/lnitial Gauge Pressure （超音速/初始表壓），如果進(jìn) 口流動是超音速的或者是準(zhǔn)備壓強(qiáng)進(jìn)口邊界條件進(jìn)行計算的初始化工作，則必須定義靜壓。在流場為亞音速時，F(xiàn)LUENT將忽略Supersonic/lnitial Gauge Pressure （超音速/初始表壓）的輸入數(shù)據(jù)，而用駐點(diǎn)參數(shù)求出靜壓。（例如，對于空氣，當(dāng)速度為 100m/S時，可先求出

49、2 1 1.225 1002 6125，靜壓等于總壓減去動壓得到。）21其動壓為12對于高雷諾數(shù)的計算， 由于采用壓強(qiáng)入口邊界條件，湍流強(qiáng)度并不易求出， 這里給出推薦范圍1%5% （當(dāng)然也可以采用預(yù)估的方法估計一下）。3.壓強(qiáng)出口邊界條件（ Pressure Outlet ）壓強(qiáng)出口邊界條件在流場出口邊界上定義靜壓，而靜壓的值僅在流場為亞音速時使用。如果在出口邊界上流場達(dá)到超音速，則邊界上的壓強(qiáng)將從流場內(nèi)部通過插值得到。其他流場變量均從流場內(nèi)部通過插值獲得。在壓強(qiáng)出口邊界上還需要定義“回流（backflow ）”條件?；亓鳁l件是在壓強(qiáng)出口邊界上出現(xiàn)回流使用的邊界條件，這樣計算將更容易收斂。FL

50、UENT在壓強(qiáng)出口邊界條件上可以使用徑向平衡條件，同時可以給定預(yù)期的流量。有關(guān)回流的湍流參數(shù)的設(shè)置如下圖所示，需要說明的是這些湍流參數(shù)只有當(dāng)回流發(fā)生時才被使用。但是即使在計算結(jié)果中沒有回流出現(xiàn)，也應(yīng)該將出口條件用真實(shí)流場的值設(shè)定， 這樣可以在計算過程中出現(xiàn)回流時加速收斂。Pressure OutletZone NaneoutletMomentum | Ihermal Radiation | Speaes | DPM | MJliphase | UDSoorstantGauge FV ess LIT已(pas call) |qBackflow Direchon Specificabon Meth

51、od formal to BoirdaryI |Ra Jal EquiBarium Prsssure Di£tjibutnnQAverage Pressure SpeciationI llarael Flow FtaterurbulenceSpecifcation Method intensity and Hydraulic DiametervBackflow Turbulent Intensitv)Backflow Hydraulfc Dianne ter (m)4,BtJ|P1P0K Cancel Hep4出流邊界條件（Outflow）出流邊界條件用于模擬在求解前流速和壓力未知的

52、出口邊界。在該邊界上，用戶不需要定義任何內(nèi)容（除非模擬輻射傳熱、粒子的離散相及多口出流）。該邊界條件適用于出口處的流動是完全發(fā)展的情況。 所謂完全發(fā)展，意味著出流面上的流動情況由區(qū)域內(nèi)部外推得到， 且對上游流動沒有影響。出游邊界條件不能用于可壓流動，也不能與壓力進(jìn)口邊界一起使用（壓力進(jìn)口邊界條件可與壓力出口邊界條件一起使用）（占總出口邊界條件的設(shè)置比較簡單，只需給定所指定的出流邊界上流體的流出量權(quán)重流出量的百分比）。如果系統(tǒng)只有一個出口，則直接輸入“1”即可。注意：在使用出流邊界條件時，如果在計算過程中，在出流邊界上的任何一點(diǎn)有回流， 計算的收斂性都會受到影響，尤其在進(jìn)行湍流計算時，這種現(xiàn)象比

53、較明顯。這里，可嘗試使用壓力出口邊界條件代替出流邊界條件。5.內(nèi)部界面（in terior ）與交界面（in terface ）這兩種面用于兩個區(qū)域的交界處，在此界面上不需要用戶輸入任何內(nèi)容，只需要指定其位置。一般內(nèi)部節(jié)點(diǎn)單元的表面都默認(rèn)是內(nèi)部界面。圍的非旋轉(zhuǎn)流體區(qū)域）的界面處，將兩個區(qū)域“隔開” 內(nèi)容，只需要指定其位置。我們注意到，相接壁面在導(dǎo)入 面生成另外一個相應(yīng)的shadow wall，將相接壁面由wall 與 shadow wall 合并為 interior 類型。內(nèi)部界面實(shí)際是兩個區(qū)域公用一個界面（只有一個面）內(nèi)部界面（interior）邊界條件用在兩個區(qū)域（如水泵中同葉輪一起旋轉(zhuǎn)的

54、流體區(qū)域與周 。在該邊界上，不需要用戶輸入任何wall。交界面（in terface）也是兩個區(qū)需要在flue ntFLUENT 時，F(xiàn)LUENT 會為該壁 更改為interior時，F(xiàn)LUENT會將域的交界，不同的是有兩個面成對出現(xiàn)。 中定義。交界面網(wǎng)格劃分可以不一樣， 這個功能使得劃分網(wǎng)格變得容易許多，對于交界面，必須是重合或部分重合， flue nt會自動在重合的部分進(jìn)行變量的插值和傳遞。 對于復(fù)雜的模型可以采取分塊劃分。如果可以熟練的將兩個區(qū)域的網(wǎng)格劃分出來，并使得交界處的節(jié)點(diǎn)統(tǒng)一（合并節(jié)點(diǎn)），那么可以不用采用in terface,直接采用in terior是最好的選擇，因?yàn)椴捎眠@種方式即提高計算速度又提高計算精 度。相反，采用in terface由于插值過程的存在，必然會降低計算速度和精度。注意：in terior界面處的網(wǎng)格必須節(jié)點(diǎn)保持一致，in terface可以不用一致，但要保證in terface是成對出現(xiàn)的。6.壁面邊界條件（wall ）一般來說，壁面邊界條件指定后不用修改，但有些情況會做一些相應(yīng)的設(shè)置，下面是一些特殊例子。離心泵旋轉(zhuǎn)的葉輪，因葉輪要隨內(nèi)部流體旋轉(zhuǎn)，故其需設(shè)置moving wa