【正文】 FLOTRAN計算流體力學全文目 錄 第一章 FLOTRAN計算流體動力學(CFD)分析概述 1第二章 FLOTRAN分析基礎 3第三章 FLOTRAN設置命令 14第四章 FLOTRAN邊界條件 59第五章 FLOTRAN層流和湍流分析算例 65第一章 FLOTRAN 計算流體動力學(CFD)分析概述FLOTRAN CFD 分析的概念ANSYS程序中的FLOTRAN CFD分析功能是一個用于分析二維及三維流體流動場的先進的工具，使用ANSYS中用于FLOTRAN CFD分析的FLUID 141和FLUID 142 單元，可解決如下問題： 作用于氣動翼(葉)型上的升力和阻力 超音速噴管中的流場 彎管中流體的復雜的三維流動同時，F(xiàn)LOTRAN還具有如下功能： 計算發(fā)動機排氣系統(tǒng)中氣體的壓力及溫度分布 研究管路系統(tǒng)中熱的層化及分離 使用混合流研究來估計熱沖擊的可能性 用自然對流分析來估計電子封裝芯片的熱性能 對含有多種流體的(由固體隔開)熱交換器進行研究FLOTRAN 分析的種類FLOTRAN可執(zhí)行如下分析： 層流或紊流 傳熱或絕熱 可壓縮或不可壓縮 牛頓流或非牛頓流 多組份傳輸這些分析類型并不相互排斥，例如，一個層流分析可以是傳熱的或者是絕熱的，一個紊流分析可以是可壓縮的或者是不可壓縮的。層流分析層流中的速度場都是平滑而有序的，高粘性流體（如石油等）的低速流動就通常是層流。紊流分析紊流分析用于處理那些由于流速足夠高和粘性足夠低從而引起紊流波動的流體流動情況，ANSYS中的二方程紊流模型可計及在平均流動下的紊流速度波動的影響。 如果流體的密度在流動過程中保持不變或者當流體壓縮時只消耗很少的能量，該流體 就可認為是不可壓縮的，不可壓縮流的溫度方程將忽略流體動能的變化和粘性耗散。熱分析流體分析中通常還會求解流場中的溫度分布情況。如果流體性質不隨溫度而變，就可不解溫度方程。在共軛傳熱問題中，要在同時包含流體區(qū)域和非流體區(qū)域（即固體區(qū)域）的整個區(qū)域上求解溫度方程。在自然對流傳熱問題中，流體由于溫度分布的不均勻性而導致流體密度分布的不均勻性，從而引起流體的流動，與強迫對流問題不同的是，自然對流通常都沒有外部的流動源?？蓧嚎s流分析對于高速氣流，由很強的壓力梯度引起的流體密度的變化將顯著地影響流場的性質，ANSYS對于這種流動情況會使用不同的解算方法。非牛頓流分析應力與應變率之間成線性關系的這種理論并不能足以解釋很多流體的流動，對于這種非牛頓流體，ANSYS程序提供了三中粘性模式和一個用戶自定義子程序。多組份傳輸分析這種分析通常是用于研究有毒流體物質的稀釋或大氣中污染氣體的傳播情況，同時，它也可用于研究有多種流體同時存在（但被固體相互隔開）的熱交換分析。 第二章 FLOTRAN分析基礎FLOTRAN單元的特點ANSYS中的FLOTRAN單元，即FLUID141和FLUID142，用于解算單相粘性流體的二維和三維流動、壓力和溫度分布。對于這些單元，ANSYS通過質量、動量和能量三個守恒性質來計算流體的速度分量、壓力、以及溫度。FLUID141單元FLUID141單元具有下列特征：維數(shù)：二維 形狀：四節(jié)點四邊形或三節(jié)點三角形 自由度：速度、壓力、溫度、紊流動能、紊流能量耗散、多達六種流體的各自質量所占的份額FLUID142單元FLUID142單元具有下列特征：維數(shù)：三維形狀：四節(jié)點四面體或八節(jié)點六面體自由度：速度、壓力、溫度、紊流動能、紊流能量耗散、多達六種流體的各自質量所占的份額FLUID141單元 FLUID142單元FLUID單元的其他特征FLUID單元的其他特征包括： 用于模擬紊流的二方程紊流模式 有很多推導結果，諸如：流場分析中的馬赫數(shù)、壓力系數(shù)、總壓、剪應力、 壁面處的yplus、以及流線函數(shù)；熱分析中的熱流、熱交換（膜）系數(shù)等。 流體邊界條件，包括：速度、壓力、紊流動能以及紊流能量耗散率。用戶無 需提供流場進口處紊流項的邊界條件，因FLOTRAN對此提供的缺省值適用 于絕大多數(shù)分析。 熱邊界條件，包括：溫度、熱流、體積熱源、熱交換（膜）系數(shù)。用戶可使用的坐標系有：的卡爾坐標系、柱坐標系、極坐標系和軸對稱坐標系。如果所計算的問題是軸對稱的，激活旋轉（swirl）選項即可算出垂直于對稱平面的速度分量。使用FLOTRAN單元的一些限制及注意事項FLOTRAN單元的一些局限性： 在同一次分析中不能改變求解的區(qū)域 單元不支持自由流面邊界條件 ANSYS程序的某些特征不能同F(xiàn)LOTRAN單元一起使用 使用FLOTRAN單元時不能使用某些命令或菜單 當使用ANSYS的圖形用戶界面時，程序將只能顯示那些在菜單和對話框中的FLOTRAN SetUp部分要求了的特征和選項。FLOTRAN單元使用中的一些限制當使用FLOTRAN單元時，要避免使用ANSYS的某些特征和命令，至少，要注意到在使用FLOTRAN單元時與別的分析稍微有些不同，當使用了無效的命令時，程序會給出相應的警告或錯誤信息。使用FLOTRAN單元要注意如下幾點（下面所列命令相應的菜單路徑請參見ANSYS命令手冊或聯(lián)機幫助中的“Commands and Their Location in the GUI” FLOTRAN單元不能和其他單元聯(lián)合使用。 節(jié)點坐標系必須與總體坐標系一致。 /CLEAR命令并不破壞業(yè)已存在的FLOTRAN結果文件（），這 有助于防止用戶不小心破壞那些花了很多時間和精力才求得的結果，用戶必 須在*作系統(tǒng)里才能刪除那些無用的結果文件。 CP命令通過對自由度進行耦合來形成周期邊界條件，ANSYS命令手冊對CP 命令的描述是可以只對某些自由度進行耦合，但作FLOTRAN分析時，周期 邊界的所有自由度都將被耦合。用戶不能對同一個單元中的節(jié)點進行耦合，相鄰單元間節(jié)點耦合也很困難。 ADAPT命令不適用于FLOTRAN分析。 不能用ANTYPE命令來引入FLOTRAN的瞬態(tài)分析。 FLOTRAN分析不支持自動時間步長功能，詳見“FLOTRAN瞬態(tài)分析”。 如果用戶通過BFCUM、BFDELE或BFUNIF來定義節(jié)點熱源，則ANSYS會 在內部用BFE命令來代替。 不能使用LDREAD,FORC命令來電磁載荷轉換到FLOTRAN分析中，而必須 使用相應的宏來進行轉換。 FLOTRAN分析不能使用CE、CECMOD、CEDELE和DEINTF命令。 FLOTRAN分析不能使用CNVTOL命令來設置收斂容差。 不能用DSYM命令來定義FLOTRAN的對稱和反對稱邊界條件。 FLOTRAN不支持旋轉坐標系中的角加速度向量。 在FLOTRAN分析中，使用FLDATA4,TIME命令而不是DELTIM命令來定義 一個載荷步的時間。 對于FLOTRAN分析，不能使用DESOL命令或PRESOL命令來修改節(jié)點的熱 （HEAT）、流（FLOW）、或流密（FLUX）結果。 FLOTRAN分析不允許將積分點結果外推到節(jié)點上（ERESX命令）。 FLOTRAN分析不允許通過鏡象*作來形成單元。 FLOTRAN分析不能用KBC命令來施加漸變載荷，而必須用多個載荷步來逐 漸改變載荷。 FLOTRAN分析不允許用LCCALC、LCDEF、LCFA、LCFI等命令來作不同載荷狀況之間的運算*作。 NCNV命令中的收斂工具不能用于FLOTRAN中相互獨立的求解器。 FLOTRAN分析不允許用NEQUIT命令來定義非線性分析的平衡迭代數(shù)。 FLOTRAN分析要求節(jié)點坐標系必須是總體的卡爾坐標系，故不能用N命令、NMODIF命令、和NROTAT命令的旋轉（rotational）域。 FLOTRAN分析中，F(xiàn)LDATA2,ITER命令和FLDATA4,TIME命令用來控制一 個載荷步中的總體迭代數(shù)。 FLOTRAN分析不允許用戶自定義單元。 FLOTRAN分析中，F(xiàn)LDATA2,ITER命令和FLDATA4,TIME命令用來對寫入 數(shù)據(jù)庫中的結果進行控制。 PRNLD命令不能用于FLOTRAN分析中，因其不能將邊界條件作為可打印的 單元節(jié)點載荷來保存。 FLOTRAN分析中沒有節(jié)點反力解。 部分和預定義求解選項（PSOLVE命令定義）不適用于FLOTRAN各自獨立 的求解器。 TIME命令不能用在FLOTRAN分析中。 FLOTRAN用FLDATA1,SOLU命令而不是TIMINT命令來定義瞬態(tài)載荷 步。 FLOTRAN用FLDATA4,TIME命令而非TRNOPT命令來定義瞬態(tài)分析選項。FLOTRAN分析的主要步驟一個典型的FLOTRAN分析有如下七個主要步驟：1. 確定問題的區(qū)域。2. 確定流體的狀態(tài)。3. 生成有限元網格。4. 施加邊界條件。5. 設置FLOTRAN分析參數(shù)。6. 求解。7. 檢查結果。第一步：確定問題的區(qū)域用戶必須確定所分析問題的明確的范圍，將問題的邊界設置在條件已知的地方，如果并不知道精確的邊界條件而必須作假定時，就不要將分析的邊界設在靠近感興趣區(qū)域的地方，也不要將邊界設在求解變量變化梯度大的地方。有時，也許用戶并不知道自己的問題中哪個地方梯度變化最大，這就要先作一個試探性的分析，然后再根據(jù)結果來修改分析區(qū)域。這些在后面章節(jié)中都有詳述。第二步：確定流體的狀態(tài)用戶在此需要估計流體的特征，流體的特征是流體性質、幾何邊界以及流場的速度幅值的函數(shù)。FLOTRAN能求解的流體包括氣流和液流，其性質可隨溫度而發(fā)生顯著變化，F(xiàn)LOTRAN中的氣流只能是理想氣體。用戶須自己確定溫度對流體的密度、粘性、和熱傳導系數(shù)的影響是否是很重要，在大多數(shù)情況下，近似認為流體性質是常數(shù)，即不隨溫度而變化，都可以得到足夠精確的解。通常用雷諾數(shù)來判別流體是層流或紊流，雷諾數(shù)反映了慣性力和粘性力的相對強度，詳見第四章。通常用馬赫數(shù)來判別流體是否可壓縮，詳見第七章。流場中任意一點的馬赫數(shù)是該點流體速度與該點音速之比值，就應考慮用可壓縮算法來進 行求解；，可壓縮算法與不可壓縮算法之間就會有極其明顯的差 異。第三步： 生成有限元網格用戶必須事先確定流場中哪個地方流體的梯度變化較大，在這些地方，網格必須作適當?shù)恼{整。例如：如果用了紊流模型，靠近壁面的區(qū)域的網格密度必須比層流模型密得多，如果太粗，該網格就不能在求解中捕捉到由于巨大的變化梯度對流動造成的顯著影響，相反，那些長邊與低梯度方向一致的單元可以有很大的長寬比。為了得到精確的結果，應使用映射網格劃分，因其能在邊界上更好地保持恒定的網格特性，映射網格劃分可由命令MSHKEY,1或其相應的菜單Main MenuPreproce ssor MeshingMeshentityMapped來實現(xiàn)。第四步：施加邊界條件可在劃分網格之前或之后對模型施加邊界條件，此時要將模型所有的邊界條件都考慮進去，如果與某個相關變量的條件沒有加上去，則該變量沿邊界的法向值的梯度將被假定為零。求解中，可在重啟動之間改變邊界條件的值，如果需改變邊界條件的值或不小心忽略了某邊界條件，可無須作重啟動，除非該改變引起了分析的不穩(wěn)定。第五步：設置FLOTRAN分析參數(shù)為了使用諸如紊流模型或求解溫度方程等選項，用戶必須激活它們。諸如流體性質等特定項目的設置，是與所求解的流體問題的類型相關的，該手冊的其他部分詳細描述了各種流體類型的所建議的參數(shù)設置。第六步：求解通過在觀察求解過程中相關變量的改變率，可以監(jiān)視求解的收斂性及穩(wěn)定性。這些變量包括速度、壓力、溫度、動能 (ENKE自由度) 和動能耗散率 (ENDS自由度) 等 紊流量以及有效粘性（EVIS）。一個分析通常需要多次重啟動。第七步：檢查結果可對輸出結果進行后處理，也可在打印輸出文件里對結果進行檢查，此時用戶應使用自己的工程經驗來估計所用的求解手段、所定義的流體性質、以及所加的邊界條件的可信程度。FLOTRAN分析中產生的一些文件在ANSYS中進行的大多數(shù)流體分析都是通過多次中斷和重啟動來完成的，通常,分析人員需要在各個重啟動之間改變諸如松弛系數(shù)等參數(shù)或開關某些項（如求解溫度方程的開關）。每當用戶繼續(xù)一個分析時，ANSYS程序會自動將數(shù)據(jù)附加在所有的 由FLOTRAN單元產生的文件中。下面將對FLOTRAN單元產生的所有文件進行說明: 結果文件， Jobname. RFL，包含節(jié)點結果。 打印文件， ，包含各量的收斂記錄及進/出口狀態(tài)(如流量等)。 壁面文件， ，包含壁面剪切應力以及YPlus信息。 殘差文件， ，包含節(jié)點殘差。 調試文件， ，包含數(shù)學求解器的有關信息。 結果備份文件， ，包含結果文件數(shù)據(jù)的一個拷貝。 重啟動文件， ，包含F(xiàn)LOTRAN的數(shù)據(jù)結構。結果文件FLOTRAN分析的結果并不自動保存在ANSYS的數(shù)據(jù)庫中，在每次求解之后。用戶可對結果文件的內容及程序 對結果文件的更新頻率進行控制，ANSYS命令手冊中對FLDATA5,OUTP命令的介紹就詳細說明了結果文件會基于用戶的選擇而保存些什么內容。在一個穩(wěn)態(tài)FLOTRAN分析中，結果文件能保存多少個結果集是沒有限制的，在求解的初期多保存幾個結果有很多好處：可以比較各結果集之間的變化、可以使用不同的選項或松弛系數(shù)來從一個分析