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

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

2、. 對(duì)于某些高導(dǎo)熱系數(shù)比或高寬縱比的網(wǎng)格，使用單精度求解器可能會(huì)遇到收斂性不佳或精確度不足不足的問(wèn)題，此時(shí)，使用雙精度求解器可能會(huì)有所幫助。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)格：通過(guò)Mesh下的Smooth/Swap來(lái)實(shí)現(xiàn)，可用來(lái)提高網(wǎng)格質(zhì)量，一般用于三角形或四邊形網(wǎng)格，不過(guò)質(zhì)量提高的效果一般般，影響較小，網(wǎng)格質(zhì)量的 提高主要還是在網(wǎng)格生成軟件里面實(shí)現(xiàn)，所以這里不再用光滑和交換的方式改善網(wǎng)格，其原理可參考FLUENT全攻略（已下載）。3. Pressure-based 與 Density-based求解器設(shè)置如圖。下面說(shuō)一說(shuō)Pressure-based和Density-based的區(qū)別：P ressure-Based Solver是Flue nt的優(yōu)勢(shì)，它是基于壓力法的求解

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

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

6、都適用于各類(lèi)流動(dòng)問(wèn)題的求解（從不可壓縮流動(dòng)到高度可壓縮流動(dòng)），但對(duì)于高速可壓縮流動(dòng)而言，使用基于密度的求解器通常能獲得比基于壓力的求解器更為精確的結(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 :是軸對(duì)稱(chēng)的意思，也就是關(guān)于一個(gè)坐標(biāo)軸對(duì)稱(chēng)，2D的axisymmetric問(wèn)題仍為2D問(wèn)題。而axisymmetric swirl :是軸對(duì)稱(chēng)旋轉(zhuǎn)的意思，就是一個(gè)區(qū)域關(guān)于一條坐標(biāo)軸 回轉(zhuǎn)所產(chǎn)生的區(qū)域，這產(chǎn)生的將是一個(gè)回轉(zhuǎn)體，是3D的問(wèn)題。在Flue nt中使用這個(gè)，是將一個(gè)3D的問(wèn)題簡(jiǎn)化為2D問(wèn)題，減少計(jì)算量，需要注意的是，在Flue nt中，回轉(zhuǎn)軸必須是x軸。5.操作工況參數(shù)（Op erati ng C

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

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

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

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

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

13、動(dòng)量平衡。同樣需要打開(kāi)Imp licit Body Force，部分平衡壓力梯度和動(dòng)量方程中體積力，提高解的收斂性。（空Referenee Pressure Location （參考?jí)簭?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. 流動(dòng)模型的設(shè)置無(wú)粘模型即真正所以理想流體是一種設(shè)想的沒(méi)有粘性的流體， 在流動(dòng)時(shí)各層之間沒(méi)有相互作用的切應(yīng)力， 沒(méi)有內(nèi)摩擦力。十分明顯，理想流體對(duì)于切向變形沒(méi)有任何抗拒能力。應(yīng)該強(qiáng)調(diào)指出，的理想流體在客觀實(shí)際中是不存在的，它只是實(shí)際流體在某些條件下的一種近似模型。在In viscid流動(dòng)模型應(yīng)用方面，無(wú)粘流動(dòng)忽略了粘性對(duì)流動(dòng)的影響，這對(duì)高雷諾數(shù)的流

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

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

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

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

19、通CFD工程仿真與案例實(shí)戰(zhàn)：FLUENT GAMBIT ICEM CFD TecplotM. 北京, 人民郵電出版社，2011:122其中，U是流體的時(shí)均速度，U是壁面摩擦速度，UW是壁面切應(yīng)力，y是壁面的垂直距離。在605時(shí)，區(qū)域?yàn)檎承缘讓?，此時(shí)速度沿壁面法線方向呈線性分布，即y 300時(shí)，流動(dòng)處于對(duì)數(shù)率層，此時(shí)速度沿壁面法線方向呈對(duì)數(shù)率分布，即u 2.5lny 5.5。使用標(biāo)準(zhǔn)的或者非平衡的壁面函數(shù)（Re 106）。O標(biāo)準(zhǔn)壁面函對(duì)于大多數(shù)高雷諾數(shù)情況壁面函數(shù)法的本質(zhì)是，對(duì)于湍流核心區(qū)的流動(dòng)使用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ù)法利用對(duì)數(shù)校正法提供了必需的壁面邊界條件（對(duì)于平衡湍流邊界層） 而非平衡壁面函數(shù)法用來(lái)改善高壓力梯度、分離、再附和滯止等情況下的結(jié)果。數(shù)法和非平衡壁面函數(shù)法都允許在近避免區(qū)域上使用較粗的網(wǎng)格。優(yōu)點(diǎn)缺點(diǎn)標(biāo)準(zhǔn)壁面函數(shù) 法應(yīng)用較多，計(jì)算量小，有較 高的精度適合高雷諾數(shù)流動(dòng)， 對(duì)低雷諾數(shù)流動(dòng)問(wèn)題， 有 壓力梯度、高度蒸騰和大的體積力、低雷諾數(shù) 和高速三維流動(dòng)問(wèn)題不適合非平衡壁面函數(shù)法考慮了壓力梯度，可以計(jì)算 分離，在附著以及撞擊問(wèn)題對(duì)低雷諾數(shù)流動(dòng)問(wèn)題，有較強(qiáng)壓力梯度、強(qiáng)體 積力及強(qiáng)三維性問(wèn)題不適合增強(qiáng)壁面處理不依賴(lài)壁

21、面法則，對(duì)于復(fù)雜 流動(dòng)，特別是低雷諾數(shù)流動(dòng) 很適合要求網(wǎng)格密，因而要求計(jì)算機(jī)處理時(shí)間長(zhǎng)，內(nèi)存大增強(qiáng)壁面處理選項(xiàng)把混合邊界模型和兩層邊界模型結(jié)合起來(lái), 近壁面現(xiàn)象很適合，湍流模型在內(nèi)層上得到了修正。表格2幾種壁面處理方法的比較對(duì)低雷諾數(shù)流動(dòng)或者復(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模型該模型通過(guò)求解單獨(dú)的動(dòng)量方程和

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

25、用于瞬變流的VOF問(wèn)題中，但必須注意的要使用該格式VOF模型必須使用顯示離散格式（在 VOF模型設(shè)置選項(xiàng)設(shè)置）。Body Force Formulati on為提高解的收斂性，對(duì)于涉及到表面張力的計(jì)算，建議在Body Force Formulation中勾選Imp licit Body Force。這樣做由于壓力梯度和動(dòng)量方程中表面張力的部分平衡，從而提高 解的收斂性。3李進(jìn)良，李承曦，胡仁喜.精通FLUENT63流場(chǎng)分析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模型這是一種簡(jiǎn)化的多相流模型，用于模擬各種有不

27、同速度的多相流，但是假定了在短空間尺度上局部的平衡。相之間的耦合應(yīng)當(dāng)是很強(qiáng)的。它也用于模擬有強(qiáng)烈耦合的各向同性多相 流和各向以相同速度運(yùn)動(dòng)的多相流。典型的應(yīng)用包括沉降(sedimentation)、氣旋分離器、低載荷作業(yè)下的多粒子流動(dòng)、氣相容積率很低的泡狀流。Mixture P arameters一般需要勾選 Mixture Parameters中的Slip Velocity復(fù)選框，以此來(lái)求解滑移速度模型, 因?yàn)樵诙嘞嗔髦懈鞣N組分的速度有很大不同。對(duì)于求解一個(gè)均勻的多相流問(wèn)題可以選擇不做滑移速度的計(jì)算，可以在 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 ）過(guò)程。在流體的固化和熔化問(wèn)題中，流場(chǎng)可以分成流體區(qū)域、固體 區(qū)域和兩者之間的糊狀區(qū)域?！办室欢嗫锥取奔夹g(shù)采用的計(jì)算策略是將流體在網(wǎng)格單元內(nèi)占有的體積百分比定義為多孔度（P orosity），并將流體和固體并存的糊狀區(qū)域看作多孔介質(zhì)區(qū) 進(jìn)行處理。在流體的固化過(guò)程中，多孔度從1降低到0;反之，在熔化過(guò)程中，多孔度則從0升至1?！办室欢嗫锥取奔夹g(shù)通過(guò)在動(dòng)量方程中添加匯項(xiàng)（即負(fù)的源項(xiàng)）模擬因固體材料 存在而出現(xiàn)的壓強(qiáng)降?！办识嗫锥取奔夹g(shù)可以模

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

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

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

33、i cstidTeb ngP 申宙 mehes陽(yáng)加心Sr”以r咖訕100000Bd(± DiffusionOlndtide Pull VeJodtiesOKCancel Help三、相設(shè)置相設(shè)置一般用于多相流的設(shè)置，對(duì)于相設(shè)置，這里主要講一下In teraction的設(shè)置，如圖:PbasesPhaseswater - P nrnay P hase air - Seccmdary Phasecdif-.rritftractan.14Interaction 設(shè)置Interaction設(shè)置用來(lái)定義兩相的相互作用，其有多個(gè)選項(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)卡針對(duì)每對(duì)物相，在下拉菜單中選擇阻力函數(shù)。其中包括schiller-naumann模型、morsi-alexander模型、symmetric （對(duì)稱(chēng)）模型等用于流體與流體之間阻力計(jì)算的模型，也 包括wen-yu模型、gidaspow模型、syamlal-obrien模型等用于液體與固體之間阻力計(jì)算的 模型，還包括 syamlal-obrien-

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

36、rame Motion選項(xiàng)確定坐標(biāo)運(yùn)動(dòng)方式(如離心泵內(nèi)部流體的旋轉(zhuǎn) 使用運(yùn)動(dòng)參考系模型)，如圖：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)口處流場(chǎng)速度及相關(guān)流動(dòng)變量作為邊界條件，在速度進(jìn)口邊界條件中，流場(chǎng)進(jìn)口邊界的駐點(diǎn)參數(shù)是不固定的。 為了滿(mǎn)足進(jìn)口處的速度條件， 駐點(diǎn)參數(shù)將在一 定范圍內(nèi)波動(dòng)。需要注意的是，因?yàn)檫@種條件中允許駐點(diǎn)參數(shù)浮動(dòng)，所以速度進(jìn)口邊界條件僅適用于不不要讓速度進(jìn)口可壓流，如果用于可壓流，則可能導(dǎo)致出現(xiàn)非物理解。同時(shí)還要注意的是， 條件過(guò)于靠近進(jìn)口內(nèi)側(cè)的固體障礙物，這樣會(huì)使駐點(diǎn)參數(shù)的不均勻程度大大增加。湍流參數(shù)的設(shè)置對(duì)于一般的流動(dòng)邊界條件，均涉及到湍流參數(shù)的定義，在Turbule nee Sp ecificatio n Method(湍流定義方法)下拉列表中，可以簡(jiǎn)單地用一個(gè)常數(shù)來(lái)定義湍流參數(shù)，即通過(guò)給

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

40、游是沒(méi)有充分發(fā)展的未受擾流 則進(jìn)口處湍流強(qiáng)度可以達(dá)I ( turbulenee intensity)I ± 0.16(ReDH) 1/8 u其中：ReH為按等效水力直徑Dh計(jì)算得到的雷諾數(shù)。湍流強(qiáng)度小于1%時(shí)，可以認(rèn)為湍流強(qiáng)度是比較低的，而在湍流強(qiáng)度大于10%時(shí)，則可以認(rèn)為湍流強(qiáng)度是比較高的。比如，當(dāng)雷諾數(shù)為50000時(shí)，代入上述公式可得到湍流強(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.湍流的長(zhǎng)度尺度

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

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

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

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

46、于0.09。在沒(méi)有直接輸入湍流動(dòng)能k和湍流耗散率的情況下，可以用Intensity and HydraulicDiameter （湍流強(qiáng)度與水力直徑）或 Intensity and Length Scale （湍流強(qiáng)度與長(zhǎng)度尺度）等辦 法，利用上述公式確定湍流動(dòng)能k和湍流耗散率。2.壓強(qiáng)進(jìn)口邊界條件（ P ressure In let ）壓強(qiáng)進(jìn)口邊界條件用于定義進(jìn)口流體的壓強(qiáng)，可用于不可壓縮和可壓縮流動(dòng)。當(dāng)進(jìn)口壓強(qiáng)已知，而流動(dòng)速度或流量未知時(shí)，可使用壓強(qiáng)進(jìn)口邊界條件。壓強(qiáng)進(jìn)口邊界條件也可用于定義外部或非受限流動(dòng)的“自由邊界”。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（動(dòng)量）選項(xiàng)卡內(nèi)，Refere nee Frame為參考值，有絕對(duì)值（ Absolute ）與Relative to Adjaeent Zone ）兩個(gè)選項(xiàng)供選擇，一般

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

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

50、UENT在壓強(qiáng)出口邊界條件上可以使用徑向平衡條件，同時(shí)可以給定預(yù)期的流量。有關(guān)回流的湍流參數(shù)的設(shè)置如下圖所示，需要說(shuō)明的是這些湍流參數(shù)只有當(dāng)回流發(fā)生時(shí)才被使用。但是即使在計(jì)算結(jié)果中沒(méi)有回流出現(xiàn)，也應(yīng)該將出口條件用真實(shí)流場(chǎng)的值設(shè)定， 這樣可以在計(jì)算過(guò)程中出現(xiàn)回流時(shí)加速收斂。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、出口邊界。在該邊界上，用戶(hù)不需要定義任何內(nèi)容（除非模擬輻射傳熱、粒子的離散相及多口出流）。該邊界條件適用于出口處的流動(dòng)是完全發(fā)展的情況。 所謂完全發(fā)展，意味著出流面上的流動(dòng)情況由區(qū)域內(nèi)部外推得到， 且對(duì)上游流動(dòng)沒(méi)有影響。出游邊界條件不能用于可壓流動(dòng)，也不能與壓力進(jìn)口邊界一起使用（壓力進(jìn)口邊界條件可與壓力出口邊界條件一起使用）（占總出口邊界條件的設(shè)置比較簡(jiǎn)單，只需給定所指定的出流邊界上流體的流出量權(quán)重流出量的百分比）。如果系統(tǒng)只有一個(gè)出口，則直接輸入“1”即可。注意：在使用出流邊界條件時(shí)，如果在計(jì)算過(guò)程中，在出流邊界上的任何一點(diǎn)有回流， 計(jì)算的收斂性都會(huì)受到影響，尤其在進(jìn)行湍流計(jì)算時(shí)，這種現(xiàn)象比

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

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