CFD常見問題(轉(zhuǎn)載）

2017-02-21  by:CAE仿真在線  來源:互聯(lián)網(wǎng)

1、流場數(shù)值計算的目的是什么?主要方法有哪些?其基本思路是什么?各自的適用范圍是什么? 答:這個問題的范疇好大啊。簡要的說一下個人的理解吧:流場數(shù)值求解的目的就是為了得到某個流動狀態(tài)下的相關參數(shù),這樣可以節(jié)省實驗經(jīng)費,節(jié)約實驗時間,并且可以模擬一些不可能做實驗的流動狀態(tài)。主要方法有有限差分,有限元和有限體積法,好像最近還有無網(wǎng)格法和波爾茲曼 法(格子法)。基本思路都是將復雜的非線性差分/積分方程簡化成簡單的代數(shù)方程。相對來說,有限差分法對網(wǎng)格的要求較高,而其他的方法就要靈活的多。 2、可壓縮流動和不可壓縮流動,在數(shù)值解法上各有何特點?為何不可壓縮流動在求解時反而比可壓縮流動有更多的困難? 答:注:這個問題不是一句兩句話就能說清楚的,大家還是看下面的兩篇小文章吧,摘自《計算流體力學應用》,讀完之后自有體會。 3、可壓縮Euler及Navier-Stokes方程數(shù)值解 描述無粘流動的基本方程組是Euler方程組,描述粘性流動的基本方程組是Navier-Stokes方程組。用數(shù)值方法通過求解Euler方程和Navier-Stokes方程模擬流場是計算流體動力學的重要內(nèi)容之一。由于飛行器設計實際問題中的絕大多數(shù)流態(tài)都具有較高的雷諾數(shù),這些流動粘性區(qū)域很小,由 對流作用主控,因此針對Euler方程發(fā)展的計算方法,在大多數(shù)情況下對Navier-Stokes方程也是有效的,只需針對粘性項用中心差分離散。 用數(shù)值方法求解無粘Euler方程組的歷史可追溯到20世紀50年代,具有代表性的方法是1952年Courant等人以及1954年Lax和Friedrichs提出的一階方法。從那時開始,人們發(fā)展了大量的差分格式。Lax和Wendroff的開創(chuàng)性工作是非定常Euler(可壓縮Navier-Stokes)方程組數(shù)值求解方法發(fā)展的 里程碑。二階精度Lax-Wendroff格式應用于非線性方程組派生出了一類格式,其共同特點是格式空間對稱,即在空間上對一維問題是三點中心格式,在時間上是顯式格式,并且該類格式是從時間空間混合離散中導出的。該類格式中最流行的是MacCormack格式。 采用時空混合離散方法,其數(shù)值解趨近于定常時依賴于計算中采用的時間步長。盡管由時間步長項引起的誤差與截斷誤差在數(shù)量級上相同,但這卻體現(xiàn)了一個概念上的缺陷,因為在計算得到的定常解中引進了一個數(shù)值參數(shù)。將時間積分從空間離散中分離出來就避免了上述缺陷。常用的時空 分別離散格式有中心型格式和迎風型格式?？臻g二階精度的中心型格式(一維問題是三點格式)就屬于上述范疇。該類格式最具代表性的是Beam-Warming隱式格式和Jameson等人采用的Runge-Kutta時間積分方法發(fā)展的顯式格式。迎風型差分格式共同特點是所建立起的特征傳播特性與差分空間 離散方向選擇的關系是與無粘流動的物理特性一致的。第一個顯式迎風差分格式是由Courant等人構(gòu)造的,并推廣為二階精度和隱式時間積分方法?；谕糠较蛐噪x散的Steger-Warming和Van Leer矢通量分裂方法可以認為是這類格式的一種。該類格式的第二個分支是Godunov方法,該方法在每個網(wǎng)格步求解描述相鄰間斷(Riemann問題)的當?shù)匾痪SEuler方程。根據(jù)這一方法Engquist、Osher和Roe等人構(gòu)造了一系列引入近似Riemann算子的格式,這就是著名的通量差分方法。 對于沒有大梯度的定常光滑流動,所有求解Euler方程格式的計算結(jié)果都是令人滿意的,但當出現(xiàn)諸如激波這樣的間斷時,其表現(xiàn)確有很大差異。絕大多數(shù)最初發(fā)展起來的格式,如Lax-Wendroff格式中心型格式,在激波附近會產(chǎn)生波動。人們通過引入人工粘性構(gòu)造了各種方法來控制和限制?廡┎ǘT諞桓鍪逼誒錚飫喔袷皆詬叢恿鞒〖撲闃械玫攪擻τ謾H歡?，由釉嶑式譃a兇雜剎問?，秷A煌侍庖械髡?，矄君l褂蒙洗戳酥疃嗖槐?，而且格式都sげǚ直媛適艿接跋歟蚨湓詬叢恿鞫撲闃械撓τ檬艿攪艘歡ㄏ拗啤? 另外一種方法是力圖阻止數(shù)值波動的產(chǎn)生,而不是在其產(chǎn)生后再進行抑制。這種方法是建立在非線性限制器的概念上,這一概念最初由Boris和Book及Van Leer提出,并且通過Harten發(fā)展的總變差減小(TVD, Total Variation Diminishing)的重要概念得以實現(xiàn)。通過這一途徑,數(shù)值解的變化以非線性的方式得以控制。這一類格式的研究和應用,在20世紀80年代形成了一股發(fā)展浪潮。1988年,張涵信和莊逢甘利用熱力學熵增原理,通過對差分格式修正方程式的分析,構(gòu)造了滿足熵增條件能夠捕捉激波的無波動、 無自由參數(shù)的耗散格式(NND格式)。該類格式在航空航天飛行器氣動數(shù)值模擬方面得到了廣泛應用。 1987年,Harten和Osher指出,TVD格式最多能達到二階精度。為了突破這一精度上的限制引入了實質(zhì)上無波動(ENO)格式的概念。該類格式“幾乎是TVD”的,Harten因此推斷這些格式產(chǎn)生的數(shù)值解是一致有界的。繼Harten和Osher之后,Shu和Osher將ENO格式從一維推廣到多維。J.Y.Yang在 三階精度ENO差分格式上也做了不少工作。1992年,張涵信另辟蹊徑,在NND格式的基礎上,發(fā)展了一種能捕捉激波的實質(zhì)上無波動、無自由參數(shù)的三階精度差分格式(簡稱ENN格式)。1994年,Liu、Osher和Chan發(fā)展了WENO(Weighted Essentially Non-Oscillatory)格式。WENO格式是基于ENO格式構(gòu)造的高階混合格式,它在保持了ENO格式優(yōu)點的同時,計算流場中虛假波動明顯減少。此后,Jiang提出了一種新的網(wǎng)格模板光滑程度的度量方法。目前高階精度格式的研究與應用是計算流體力學的熱點問題之一。 不可壓縮Navier-Stokes方程求解 不可壓縮流體力學數(shù)值解法有非常廣泛的需求。從求解低速空氣動力學問題,推進器內(nèi)部流動,到水動力相關的液體流動以及生物流體力學等。滿足這么廣泛問題的研究,要求有與之相應的較好的物理問題的數(shù)學模型以及魯棒的數(shù)值算法。 相對于可壓縮流動,不可壓縮流動的數(shù)值求解困難在于,不可壓縮流體介質(zhì)的密度保持常數(shù),而狀態(tài)方程不再成立,連續(xù)方程退化為速度的散度為零的方程。由此,在可壓縮流動的計算中可用于求解密度和壓力的連續(xù)方程在不可壓縮流動求解中僅是動量方程的一個約束條件,由此求解不可 壓縮流動的壓力稱為一個困難。求解不可壓縮流動的各種方法主要在于求解不同的壓力過程。 目前,主要有兩類求解不可壓縮流體力學的方法,原始變量方法和非原始變量方法。求解不可壓縮流動的原始變量方法是將Navier-Stokes方程寫成壓力和速度的形式,進行直接求解,這種形式已被廣為應用。非原始變量方法主要有Fasel提出的流函數(shù)-渦函數(shù)法、Aziz和Hellums提出的勢函 數(shù)-渦函數(shù)方法。在求解三維流動問題時,上述每一個方法都需要反復求解三個Possion方程,非常耗時。原始變量方法可以分為三類:第一種方法是Harlow和Welch首先提出的壓力Possion方程方法。該方法首先將動量方程推進求得速度場,然后利用Possion方程求解壓力,這一種方法由于?懇皇奔洳繳閑枰蠼釶ossion方程,求解非常耗時。第二種方法是Patanker和Spalding的SIMPLE(Semi-Implicit Method for Pressure-Linked Equation)法,它是通過動量方程求得壓力修正項對速度的影響,使其滿足速度散度等于零的條件作為壓力控制方程。第三種方法是虛擬壓縮方法,這一方法是Chorin于1967年提出的。該方法的核心就是通過在連續(xù)方程中引入一個虛擬壓縮因子,再附加一項壓力的虛擬時間導數(shù),使壓力顯?降賾胨俁攘燈鵠?，同湿N匠桃脖涑閃慫頭匠獺U庋?，方程的袉栻溚与求斤兩杨欟流动的方程相似，因此，写夃求斤兩杨欟流动的硥勳方法都可用诱E豢裳顧趿鞫那蠼狻? 目前,由于基于求解壓力Possion方程的方法非常復雜及耗時,難于求解具體的工程實際問題,因此用此方法解決工程問題的工作越來越少。工程上常用的主要是SIMPLE方法和虛擬壓縮方法。 4、什么叫邊界條件?有何物理意義?它與初始條件有什么關系? 邊界條件與初始條件是控制方程有確定解的前提。 邊界條件是在求解區(qū)域的邊界上所求解的變量或其導數(shù)隨時間和地點的變化規(guī)律。對于任何問題,都需要給定邊界條件。 初始條件是所研究對象在過程開始時刻各個求解變量的空間分布情況,對于瞬態(tài)問題,必須給定初始條件,穩(wěn)態(tài)問題,則不用給定。對于邊界條件與初始條件的處理,直接影響計算結(jié)果的精度。 在瞬態(tài)問題中,給定初始條件時要注意的是:要針對所有計算變量,給定整個計算域內(nèi)各單元的初始條件;初始條件一定是物理上合理的,要靠經(jīng)驗或?qū)崪y結(jié)果。 5、在數(shù)值計算中,偏微分方程的雙曲型方程、橢圓型方程、拋物型方程有什么區(qū)別? 我們知道很多描述物理問題的控制方程最終就可以歸結(jié)為偏微分方程,描述流動的控制方程也不例外。 從數(shù)學角度,一般將偏微分方程分為橢圓型(影響域是橢圓的,與時間無關,且是空間內(nèi)的閉區(qū)域,故又稱為邊值問題),雙曲型(步進問題,但依賴域僅在兩條特征區(qū)域之間),拋物型(影響域以特征線為分界線,與主流方向垂直;具體來說,解的分布與瞬時以前的情況和邊界條件相關 ,下游的變化僅與上游的變化相關;也稱為初邊值問題); 從物理角度,一般將方程分為平衡問題(或穩(wěn)態(tài)問題),時間步進問題。 兩種角度,有這樣的關系:橢圓型方程描述的一般是平衡問題(或穩(wěn)態(tài)問題),雙曲型和拋物型方程描述的一般是步進問題。 至于具體的分類方法,大家可以參考一般的偏微分方程專著,里面都有介紹。關于各種不同近似水平的流體控制方程的分類,大家可以參考張涵信院士編寫《計算流體力學—差分方法的原理與應用》里面講的相當詳細。 三種類型偏微分方程的基本差別如下: 1)三種類型偏微分方程解的適定性(即解存在且唯一,并且解穩(wěn)定)要求對定解條件有不同的提法; 2)三種類型偏微分方程解的光滑性不同,對定解條件的光滑性要求也不同; 橢圓型和拋物型方程的解是充分光滑的,因此對定解條件的光滑性要求不高。而雙曲型方程允許有所謂的弱解存在(如流場中的激波),即解的一階導數(shù)可以不連續(xù),所以對定解條件的光滑性要求很高,這也正是采用有限元法求解雙曲型方程困難較多的原因之一。 3)三種類型偏微分方程的影響區(qū)域和依賴區(qū)域不一樣。 在雙曲型和拋物型方程所控制的流場中,某一點的影響區(qū)域是有界的,可采用步進求解。如對雙曲型方程求解時,為了與影響區(qū)域的特征一致,采用上風格式比較適宜。而橢圓型方程的影響范圍遍及全場,必須全場求解,所采用的差分格式也要采用相應的中心格式。 6、在網(wǎng)格生成技術中,什么叫貼體坐標系?什么叫網(wǎng)格獨立解? 數(shù)值計算的與實驗值之間的誤差來源只要有這幾個:物理模型近似誤差(無粘或有粘,定常與非定常,二維或三維等等),差分方程的截斷誤差及求解區(qū)域的離散誤差(這兩種誤差通常統(tǒng)稱為離散誤差),迭代誤差(離散后的代數(shù)方程組的求解方法以及迭代次數(shù)所產(chǎn)生的誤差),舍入誤差 (計算機只能用有限位存儲計算物理量所產(chǎn)生的誤差)等等。在通常的計算中,離散誤差隨網(wǎng)格變細而減小,但由于網(wǎng)格變細時,離散點數(shù)增多,舍入誤差也隨之加大。 由此可見,網(wǎng)格數(shù)量并不是越多越好的。 再說說網(wǎng)格無關性的問題,由上面的介紹,我們知道網(wǎng)格數(shù)太密或者太疏都可能產(chǎn)生誤差過大的計算結(jié)果,網(wǎng)格數(shù)在一定的范圍內(nèi)的結(jié)果才與實驗值比較接近,這樣在劃分網(wǎng)格時就要求我們首先依據(jù)已有的經(jīng)驗大致劃分一個網(wǎng)格進行計算,將計算結(jié)果(當然這個計算結(jié)果必須是收斂的)與 實驗值進行比較(如果沒有實驗值,則不需要比較,后面的比較與此類型相同),再酌情加密或減少網(wǎng)格,再進行計算,再與實驗值進行比較,并與前一次計算結(jié)果比較,如果兩次的計算結(jié)果相差較小(例如在2%),說明這一范圍的網(wǎng)格的計算結(jié)果是可信的,說明計算結(jié)果是網(wǎng)格無關的。 再加密網(wǎng)格已經(jīng)沒有什么意義(除非你要求的計算精度較高)。但是,如果你用粗網(wǎng)格也能得到相差很小的計算結(jié)果,從計算效率上講,你就可以完全使用粗網(wǎng)格去完成你的計算。加密或者減少網(wǎng)格數(shù)量,你可以以一倍的量級進行。 7、在GAMBIT中顯示的“check”主要通過哪幾種來判斷其網(wǎng)格的質(zhì)量?及其在做網(wǎng)格時大致注意到哪些細節(jié)? 判斷網(wǎng)格質(zhì)量的方面有: Area單元面積,適用于2D單元,較為基本的單元質(zhì)量特征。 Aspect Ratio長寬比,不同的網(wǎng)格單元有不同的計算方法,等于1是最好的單元,如正三角形,正四邊形,正四面體,正六面體等;一般情況下不要超過5:1. Diagonal Ratio對角線之比,僅適用于四邊形和六面體單元,默認是大于或等于1的,該值越高,說明單元越不規(guī)則,最好等于1,也就是正四邊形或正六面體。 Edge Ratio長邊與最短邊長度之比,大于或等于1,最好等于1,解釋同上。 EquiAngle Skew通過單元夾角計算的歪斜度,在0到1之間,0為質(zhì)量最好,1為質(zhì)量最差。最好是要控制在0到0.4之間。 EquiSize Skew通過單元大小計算的歪斜度,在0到1之間,0為質(zhì)量最好,1為質(zhì)量最差。2D質(zhì)量好的單元該值最好在0.1以內(nèi),3D單元在0.4以內(nèi)。 MidAngle Skew通過單元邊中點連線夾角計算的歪斜度,僅適用于四邊形和六面體單元,在0到1之間,0為質(zhì)量最好,1為質(zhì)量最差。 Size Change相鄰單元大小之比,僅適用于3D單元,最好控制在2以內(nèi)。 Stretch伸展度。通過單元的對角線長度與邊長計算出來的,僅適用于四邊形和六面體單元,在0到1之間,0為質(zhì)量最好,1為質(zhì)量最差。 Taper錐度。僅適用于四邊形和六面體單元,在0到1之間,0為質(zhì)量最好,1為質(zhì)量最差。 Volume單元體積,僅適用于3D單元,劃分網(wǎng)格時應避免出現(xiàn)負體積。 Warpage翹曲。僅適用于四邊形和六面體單元,在0到1之間,0為質(zhì)量最好,1為質(zhì)量最差。 以上只是針對Gambit幫助文件的簡單歸納,不同的軟件有不同的評價單元質(zhì)量的指標,使用時最好仔細閱讀幫助文件。 另外,在Fluent中的窗口鍵入:grid quality 然后回車,Fluent能檢查網(wǎng)格的質(zhì)量,主要有以下三個指標: 1.Maxium cell squish: 如果該值等于1,表示得到了很壞的單元; 2.Maxium cell skewness: 該值在0到1之間,0表示最好,1表示最壞; 3.Maxium 'aspect-ratio': 1表示最好。 8、在兩個面的交界線上如果出現(xiàn)網(wǎng)格間距不同的情況時,即兩塊網(wǎng)格不連續(xù)時,怎么樣克服這種情況呢? 這個問題就是非連續(xù)性網(wǎng)格的設置,一般來說就是把兩個交接面設置為一對interface。 另外,作此操作可能出現(xiàn)的問題及可供參考的解決方法為: 問題:把兩個面(其中一個實際是由若干小面組成,將若干小面定義為了group了)拼接在一起,也就是說兩者之間有流體通過,兩個面?zhèn)€屬不同的體,網(wǎng)格導入到fluent時,使用interface時出現(xiàn)網(wǎng)格check的錯誤,將interface的邊界條件刪除,就不會發(fā)生網(wǎng)格檢查的錯誤,如何將兩個面的 網(wǎng)格相連? 原因:interface后的兩個體的交接面,fluent以將其作為內(nèi)部流體處理(非重疊部分默認為wall,合并后網(wǎng)格會在某些地方發(fā)生畸變,導致合并失敗,也可能準備合并的兩個面幾何位置有誤差,應該準確的在同一幾何位置(合并的面大小相等時),在合并之前要合理分塊。 解決方法:為了避免網(wǎng)格發(fā)生畸變(可能一個面上的網(wǎng)格跑到另外的面上了),可以一面網(wǎng)格粗,一面網(wǎng)格細避免; 再者就是通過將一個面的網(wǎng)格直接映射到另一面上的,兩個面默認為interior.也可以將網(wǎng)格拼接一起. 9、在設置GAMBIT邊界層類型時需要注意的幾個問題:a、沒有定義的邊界線如何處理?b、計算域內(nèi)的內(nèi)部邊界如何處理(2D)? 答:gambit默認為wall,一般情況下可以到fluent再修改邊界類型。 內(nèi)部邊界如果是split產(chǎn)生的,那么就不需再設定了,如果不是,那么就需要設定為interface或者是internal 10、為何在劃分網(wǎng)格后,還要指定邊界類型和區(qū)域類型?常用的邊界類型和區(qū)域類型有哪些? 答:要得到一個問題的定解就需要有定解條件,而邊界條件就屬于定解條件。也就是說邊界條件確定了結(jié)果。不同的流體介質(zhì)有不同的物理屬性,也就會得到不同的結(jié)果,所以必須指定區(qū)域類型。對于gambit來說,默認的區(qū)域類型是fluid,所以一般情況下不需要再指定 20 何為流體區(qū)域(fluid zone)和固體區(qū)域(solid zone)?為什么要使用區(qū)域的概念?FLUENT是怎樣使用區(qū)域的? Fluid Zone是一個單元組,是求解域內(nèi)所有流體單元的綜合。所激活的方程都要在這些單元上進行求解。向流體區(qū)域輸入的信息只是流體介質(zhì)(材料)的類型。對于當前材料列表中沒有的材料,需要用戶自行定義。注意,多孔介質(zhì)也當作流體區(qū)域?qū)Υ? Solid Zone也是一個單元組,只不過這組單元僅用來進行傳熱計算,不進行任何的流動計算。作為固體處理的材料可能事實上是流體,但是假定其中沒有對流發(fā)生,固體區(qū)域僅需要輸入材料類型。 Fluent中使用Zone的概念,主要是為了區(qū)分分塊網(wǎng)格生成,邊界條件的定義等等; 21 如何監(jiān)視FLUENT的計算結(jié)果?如何判斷計算是否收斂?在FLUENT中收斂準則是如何定義的?分析計算收斂性的各控制參數(shù),并說明如何選擇和設置這些參數(shù)?解決不收斂問題通常的幾個解決方法是什么? 可以采用殘差控制面板來顯示;或者采用通過某面的流量控制;如監(jiān)控出口上流量的變化;采用某點或者面上受力的監(jiān)視;渦街中計算達到收斂時,繞流體的面上受的升力為周期交變,而阻力為平緩的直線。 怎樣判斷計算結(jié)果是否收斂? 1、觀察點處的值不再隨計算步驟的增加而變化; 2、各個參數(shù)的殘差隨計算步數(shù)的增加而降低,最后趨于平緩; 3、要滿足質(zhì)量守恒(計算中不牽涉到能量)或者是質(zhì)量與能量守恒(計算中牽涉到能量)。 特別要指出的是,即使前兩個判據(jù)都已經(jīng)滿足了,也并不表示已經(jīng)得到合理的收斂解了,因為,如果松弛因子設置得太緊,各參數(shù)在每步計算的變化都不是太大,也會使前兩個判據(jù)得到滿足。此時就要再看第三個判據(jù)了。 還需要說明的就是,一般我們都希望在收斂的情況下,殘差越小越好,但是殘差曲線是全場求平均的結(jié)果,有時其大小并不一定代表計算結(jié)果的好壞,有時即使計算的殘差很大,但結(jié)果也許是好的,關鍵是要看計算結(jié)果是否符合物理事實,即殘差的大小與模擬的物理現(xiàn)象本身的復雜性有關?匭氪郵導飾錮硐窒笊峽醇撲憬峁1熱縊狄桓鋈P停詿蠊ソ喬榭魷?解震蕩得非常厲害,而且殘差的量級也總下不去,但這仍然是正確的,為什么呢,因為大攻角下實際流動情形就是這樣的,不斷有渦的周期性脫落,流場本身就是非定常的,所以解也是波動的,處理的時候取平均?涂梢閱? 22 什么叫松弛因子?松弛因子對計算結(jié)果有什么樣的影響?它對計算的收斂情況又有什么樣的影響? 1、亞松馳(Under Relaxation):所謂亞松馳就是將本層次計算結(jié)果與上一層次結(jié)果的差值作適當縮減,以避免由于差值過大而引起非線性迭代過程的發(fā)散。用通用變量來寫出時,為松馳因子(Relaxation Factors)?！稊?shù)值傳熱學-214》 2、FLUENT中的亞松馳:由于FLUENT所解方程組的非線性,我們有必要控制的變化。一般用亞松馳方法來實現(xiàn)控制,該方法在每一部迭代中減少了的變化量。亞松馳最簡單的形式為:單元內(nèi)變量等于原來的值 加上亞松馳因子a與 變化的積, 分離解算器使用亞松馳來控制每一步迭代中的計算變量的更新。這就意味著使用分離解算器解的方程,包括耦合解算器所解的非耦合方程(湍流和其他標量)都會有一個相關的亞松馳因子。在FLUENT中,所有變量的默認亞松馳因子都是對大多數(shù)問題的最優(yōu)值。這個值適合于很多問題,但是 對于一些特殊的非線性問題(如:某些湍流或者高Rayleigh數(shù)自然對流問題),在計算開始時要慎重減小亞松馳因子。使用默認的亞松馳因子開始計算是很好的習慣。如果經(jīng)過4到5步的迭代殘差仍然增長,你就需要減小亞松馳因子。有時候,如果發(fā)現(xiàn)殘差開始增加,你可以改變亞松馳因子 重新計算。在亞松馳因子過大時通常會出現(xiàn)這種情況。最為安全的方法就是在對亞松馳因子做任何修改之前先保存數(shù)據(jù)文件,并對解的算法做幾步迭代以調(diào)節(jié)到新的參數(shù)。最典型的情況是,亞松馳因子的增加會使殘差有少量的增加,但是隨著解的進行殘差的增加又消失了。如果殘差變化有 幾個量級你就需要考慮停止計算并回到最后保存的較好的數(shù)據(jù)文件。注意:粘性和密度的亞松馳是在每一次迭代之間的。而且,如果直接解焓方程而不是溫度方程(即:對PDF計算),基于焓的溫度的更新是要進行亞松馳的。要查看默認的亞松弛因子的值,你可以在解控制面板點擊默認按?ァ6雜詿蠖嗍鞫?，不需要修改默认亚松痴涷子?牽綣魷植晃榷ɑ蛘叻⑸⒛憔托枰跣∧系難撬沙諞蜃恿?，其中压力、动量、k和e的亞松弛因子默認值分別為0.2,0.5,0.5和0.5。對于SIMPLEC格式一般不需要減小壓力的亞松弛因子。在密度和溫度強烈耦合的問題中,如相當高 的Rayleigh數(shù)的自然或混合對流流動,應該對溫度和/或密度(所用的亞松弛因子小于1.0)進行亞松弛。相反,當溫度和動量方程沒有耦合或者耦合較弱時,流動密度是常數(shù),溫度的亞松弛因子可以設為1.0。對于其它的標量方程,如漩渦,組分,PDF變量,對于某些問題默認的亞松弛可能 過大,尤其是對于初始計算。你可以將松弛因子設為0.8以使得收斂更容易。 SIMPLE與SIMPLEC比較 在FLUENT中,可以使用標準SIMPLE算法和SIMPLEC(SIMPLE-Consistent)算法,默認是SIMPLE算法,但是對于許多問題如果使用SIMPLEC可能會得到更好的結(jié)果,尤其是可以應用增加的亞松馳迭代時,具體介紹如下: 對于相對簡單的問題(如:沒有附加模型激活的層流流動),其收斂性已經(jīng)被壓力速度耦合所限制,你通常可以用SIMPLEC算法很快得到收斂解。在SIMPLEC中,壓力校正亞松馳因子通常設為1.0,它有助于收斂。但是,在有些問題中,將壓力校正松弛因子增加到1.0可能會導致不穩(wěn)定。對于 所有的過渡流動計算,強烈推薦使用PISO算法鄰近校正。它允許你使用大的時間步,而且對于動量和壓力都可以使用亞松馳因子1.0。對于定常狀態(tài)問題,具有鄰近校正的PISO并不會比具有較好的亞松馳因子的SIMPLE或SIMPLEC好。對于具有較大扭曲網(wǎng)格上的定常狀態(tài)和過渡計算推薦使用PI SO傾斜校正。當你使用PISO鄰近校正時,對所有方程都推薦使用亞松馳因子為1.0或者接近1.0。如果你只對高度扭曲的網(wǎng)格使用PISO傾斜校正,請設定動量和壓力的亞松馳因子之和為1.0比如:壓力亞松馳因子0.3,動量亞松馳因子0.7)。如果你同時使用PISO的兩種校正方法,推薦參閱PIS O鄰近校正中所用的方法 23 在FLUENT運行過程中,經(jīng)常會出現(xiàn)“turbulence viscous rate”超過了極限值,此時如何解決?而這里的極限值指的是什么值?修正后它對計算結(jié)果有何影響 Let's take care of the warning "turbulent viscosity limited to viscosity ratio****" which is not physical. This problem is mainly due to one of the following: 1)Poor mesh quality(i.e.,skewness > 0.85 for Quad/Hex, or skewness > 0.9 for Tri/Tetra elements). {what values do you have?} 2)Use of improper turbulent boudary conditions. 3)Not supplying good initial values for turbulent quantities. 出現(xiàn)這個警告,一般來講,最可能的就是網(wǎng)格質(zhì)量的問題,尤其是Y+值的問題;在劃分網(wǎng)格的時候要注意,第一層網(wǎng)格高度非常重要,可以使用NASA的Viscous Grid Space Calculator來計算第一層網(wǎng)格高度;如果這方面已經(jīng)注意了,那就可能是邊界條件中有關湍流量的設置問題, 24 在FLUENT運行計算時,為什么有時候總是出現(xiàn)“reversed flow”?其具體意義是什么?有沒有辦法避免?如果一直這樣顯示,它對最終的計算結(jié)果有什么樣的影響? 這個問題的意思是出現(xiàn)了回流,這個問題相對于湍流粘性比的警告要寬松一些,有些case可能只在計算的開始階段出現(xiàn)這個警告,隨著迭代的計算,可能會消失,如果計算一段時間之后,警告消失了,那么對計算結(jié)果是沒有什么影響的,如果這個警告一直存在,可能需要作以下處理: 1.如果是模擬外部繞流,出現(xiàn)這個警告的原因可能是邊界條件取得距離物體不夠遠,如果邊界條件取的足夠遠,該處可能在計算的過程中的確存在回流現(xiàn)象;對于可壓縮流動,邊界最好取在10倍的物體特征長度之處;對于不可壓縮流動,邊界最好取在4倍的物體特征長度之處。 2.如果出現(xiàn)了這個警告,不論對于外部繞流還是內(nèi)部流動,可以使用pressure-outlet邊界條件代替outflow邊界條件改善這個問題。 26 什么叫問題的初始化?在FLUENT中初始化的方法對計算結(jié)果有什么樣的影響?初始化中的“patch”怎么理解? 問題的初始化就是在做計算時,給流場一個初始值,包括壓力、速度、溫度和湍流系數(shù)等。理論上,給的初始場對最終結(jié)果不會產(chǎn)生影響,因為隨著跌倒步數(shù)的增加,計算得到的流場會向真實的流場無限逼近,但是,由于Fluent等計算軟件存在像離散格式精度(會產(chǎn)生離散誤差)和截斷誤 差等問題的限制,如果初始場給的過于偏離實際物理場,就會出現(xiàn)計算很難收斂,甚至是剛開始計算就發(fā)散的問題。因此,在初始化時,初值還是應該給的盡量符合實際物理現(xiàn)象。這就要求我們對要計算的物理場,有一個比較清楚的理解。 初始化中的patch就是對初始化的一種補充,比如當遇到多相流問題時,需要對各相的參數(shù)進行更細的限制,以最大限度接近現(xiàn)實物理場。這些就可以通過patch來實現(xiàn),patch可以對流場分區(qū)進行初始化,還可以通過編寫簡單的函數(shù)來對特定區(qū)域初始化。 27 什么叫PDF方法?FLUENT中模擬煤粉燃燒的方法有哪些? 概率密度函數(shù)輸運輸運方程方法(PDF方法)是近年來逐步建立起來的描述湍流兩相流動的新模型方法。所謂的概率密度函數(shù)(Probability Density Function,簡稱PDF)方法是基于湍流場隨機性和概率統(tǒng)計描述,將流場的速度、溫度和組分濃度等特征量作為隨機變量,研究其概率密度函數(shù)在相空間的傳遞行為的研究方法。PDF模型介于統(tǒng)觀模擬和細觀模擬之間,是從隨機運動的分子動力論和兩相湍流的基本守恒定律出發(fā),探討兩相湍流 的規(guī)律,因此可作為發(fā)展雙流體模型框架內(nèi)兩相湍流模型的理論基礎。它實質(zhì)上是溝通E-L模型和E-E模型的橋梁,可以用顆粒運動的拉氏分析通過統(tǒng)計理論,即PDF方程的積分建立封閉的E-E兩相湍流模型 非預混湍流燃燒過程的正確模擬要求同時模擬混合和化學反應過程。FLUENT 提供了四種反應模擬方法:即有限率反應法、混合分數(shù)PDF 法、不平衡(火焰微元)法和預混燃燒法?；鹧嫖⒃ㄊ腔旌戏謹?shù)PDF 方法的一種特例。該方法是基于不平衡反應的,混合分數(shù)PDF 法不能模擬的不平衡現(xiàn)象如火焰的懸舉和熄滅,NOx 的形成等都可用該方法模擬。但由于該方法還未完善,在FLUENT 只能適用于絕熱模型。 對許多燃燒系統(tǒng),輻射式主要的能量傳輸方式,因此在模擬燃燒系統(tǒng)時,對輻射能量的傳輸?shù)哪M也是非常重要的。在FLUENT 中,對于模擬該過程的模型也是非常全面的。包括DTRM、P-1、Rosseland、DO 輻射模型,還有用WSGG 模型來模擬吸收系數(shù)。 30 FLUENT運行過程中,出現(xiàn)殘差曲線震蕩是怎么回事?如何解決殘差震蕩的問題?殘差震蕩對計算收斂性和計算結(jié)果有什么影響? 31 數(shù)值模擬過程中,什么情況下出現(xiàn)偽擴散的情況?以及對于偽擴散在數(shù)值模擬過程中如何避免? 假擴散(false diffusion)的含義: 基本含義:由于對流—擴散方程中一階導數(shù)項的離散格式的截斷誤差小于二階而引起較大數(shù)值計算誤差的現(xiàn)象。有的文獻中將人工粘性(artificial viscosity)或數(shù)值粘性(numerical viscosity)視為它的同義詞。 拓寬含義:現(xiàn)在通常把以下三種原因引起的數(shù)值計算誤差都歸在假擴散的名稱下 1.非穩(wěn)態(tài)項或?qū)α黜棽捎靡浑A截差的格式; 2.流動方向與網(wǎng)格線呈傾斜交叉(多維問題); 3.建立差分格式時沒有考慮到非常數(shù)的源項的影響。 克服或減輕假擴散的格式或方法, 為克服或減輕數(shù)值計算中的假擴散(包括流向擴散及交叉擴散)誤差,應當: 1. 采用截差階數(shù)較高的格式; 2. 減輕流線與網(wǎng)格線之間的傾斜交叉現(xiàn)象或在構(gòu)造格式時考慮到來流方向的影響。 3. 至于非常數(shù)源項的問題,目前文獻中,還沒有為克服這種影響而專門構(gòu)造的格式,但是高階格式顯然對減輕其影響是有利的。 32 FLUENT輪廓(contour)顯示過程中,有時候標準輪廓線顯示通常不能精確地顯示其細節(jié),特別是對于封閉的3D物體(如柱體),其原因是什么?如何解決? FLUENT等高線(contour)顯示過程中,可以通過調(diào)節(jié)顯示的水平等級來調(diào)節(jié)其顯示細節(jié),Levels...最大值允許設置為100.對于封閉的3D物體,可以通過建立Surface,監(jiān)視Surface上的量來顯示計算結(jié)果?；蛘哂嬎阒髮⒔Y(jié)果導入到Tecplot中,作切片圖顯示。 33 如果采用非穩(wěn)態(tài)計算完畢后,如何才能更形象地顯示出動態(tài)的效果圖? 對于非定常計算,可以通過創(chuàng)建動畫來形象地顯示出動態(tài)的效果圖。 Solve->Animate->Define...,具體操作請參考Fluent用戶手冊。 34 在FLUENT的學習過程中,通常會涉及幾個壓力的概念,比如壓力是相對值還是絕對值?參考壓力有何作用?如何設置和利用它? GAUGE PRESSURE 就是靜壓。 GAUGE total PRESSURE 是總壓。 這里需要強調(diào)一下 Gauge為名義值, 什么意思呢?如果, INITIAL Gauge PRESSURE =0 那么 GAUGE PRESSURE 就是實際的靜壓Pinf。 GAUGE total PRESSURE 是實際的總壓Pt。 如果INITIAL Gauge PRESSURE 不等于零 GAUGE PRESSURE = Pinf - INITIAL Gauge PRESSURE GAUGE total PRESSURE = Pt - INITIAL Gauge PRESSURE 35 在FLUENT結(jié)果的后處理過程中,如何將美觀漂亮的定性分析的效果圖和定量分析示意圖插入到論文中來說明問題? 三種方法來得到用于插入到論文的圖片: 1.在Fluent中顯示你想得到的效果圖的窗口,可以直接在任務欄中右鍵該窗口將其復制到剪貼板,保存;或者打印到文件,保存。 2.在Fluent中,在你想要保存相關窗口的效果圖時,首先激活效果圖監(jiān)視窗口,就是用鼠標左鍵監(jiān)視窗口,然后在Fluent中操作,Fluent->File->Hardcopy...,選擇好你想要的圖片格式,然后就可以保存了。 3.將計算結(jié)果或者相關數(shù)據(jù)導入到Tecplot中,然后作出你想要的效果圖,這種方法得出的圖片,個人感覺比Fluent得到的圖片美觀簡潔大方 36 在DPM模型中,粒子軌跡能表示粒子在計算域內(nèi)的行程,如何顯示單一粒徑粒子的軌道(如20微米的粒子)? 首先選擇DMP模型,在set injection properties 面板中,選擇injection type的類型為single, 然后設置初始條件,如位置(x,y,z),速度,直徑(如20微米的粒子),溫度,質(zhì)量流率等! 設定完成后,你就可以行迭代了。等氣相和離散相收斂以后,你就可以追蹤粒子軌跡。在display中打開particle tracks面板進行操作! 37 在FLUENT定義速度入口時,速度入口的適用范圍是什么?湍流參數(shù)的定義方法有哪些?各自有什么不同? 速度入口的邊界條件適用于不可壓流動,需要給定進口速度以及需要計算的所有標量值。速度入口邊界條件不適合可壓縮流動,否則入口邊界條件會使入口處的總溫或總壓有一定的波動。 關于湍流參數(shù)的定義方法,根據(jù)所選擇的湍流模型的不同有不同的湍流參數(shù)組合,具體可以參考Fluent用戶手冊的相關章節(jié),也可以參考王福軍的書《計算流體動力學分析—CFD軟件原理與應用》的第214-216頁, 38 在計算完成后,如何顯示某一斷面上的溫度值?如何得到速度矢量圖?如何得到流線? 這些都可以用tecplot來處理 將fluent計算的date和case文件倒入到tecplot中 斷面可以做切片 速度矢量圖流線圖 直接就可以選擇相應選項來查看 39 分離式求解器和耦合式求解器的適用場合是什么?分析兩種求解器在計算效率與精度方面的區(qū)別 分離式求解器以前主要用于不可壓縮流動和微可壓流動,而耦合式求解器用于高速可壓流動?，F(xiàn)在,兩種求解器都適用于從不可壓到高速可壓的很大范圍的流動,但總的來講,當計算高速可壓流動時,耦合式求解器比分離式求解器更有優(yōu)勢。 Fluent默認使用分離式求解器,但是,對于高速可壓流動,由強體積力(如浮力或者旋轉(zhuǎn)力)導致的強耦合流動,或者在非常精細的網(wǎng)格上求解的流動,需要考慮耦合式求解器。耦合式求解器耦合了流動和能量方程,常常很快便可以收斂。耦合式求解器所需要的內(nèi)存約是分離式求解器的1.5??倍,選擇時可以根據(jù)這一情況來權衡利弊。在需要耦合隱式的時候,如果計算機內(nèi)存不夠,就可以采用分離式或耦合顯式。耦合顯式雖然也耦合了流動和能量方程,但是它還是比耦合隱式需要的內(nèi)存少,當然它的收斂性也相應差一些。 需要注意的是,在分離式求解器中提供的幾個物理模型,在耦合式求解器中是沒有的。這些物理模型包括:流體體積模型(VOF),多項混合模型,歐拉混合模型,PDF燃燒模型,預混合燃燒模型,部分預混合燃燒模型,煙灰和NOx模型,Rosseland輻射模型,熔化和凝固等相變模型,指定質(zhì)量 流量的周期流動模型,周期性熱傳導模型和殼傳導模型等。 而下列物理模型只在耦合式求解器中有效,在分離式求解器中無效:理想氣體模型,用戶定義的理想氣體模型,NIST理想氣體模型,非反射邊界條件和用于層流火焰的化學模型。 43 FLUENT中常用的文件格式類型:dbs,msh,cas,dat,trn,jou,profile等有什么用處? 在Gambit目錄中,有三個文件,分別是default_id.dbs,jou,trn文件,對Gambit運行save,將會在工作目錄下保存這三個文件:default_id.dbs,default_id.jou,default_id.trn。 jou文件是gambit命令記錄文件,可以通過運行jou文件來批處理gambit命令; dbs文件是gambit默認的儲存幾何體和網(wǎng)格數(shù)據(jù)的文件; trn文件是記錄gambit命令顯示窗(transcript)信息的文件; msh文件可以在gambit劃分網(wǎng)格和設置好邊界條件之后export中選擇msh文件輸出格式,該文件可以被fluent求解器讀取。 Case文件包括網(wǎng)格,邊界條件,解的參數(shù),用戶界面和圖形環(huán)境。 Data文件包含每個網(wǎng)格單元的流動值以及收斂的歷史紀錄(殘差值)。Fluent自動保存文件類型,默認為date和case文件 Profile文件邊界輪廓用于指定求解域的邊界區(qū)域的流動條件。例如,它們可以用于指定入口平面的速度場。 讀入輪廓文件,點擊菜單File/Read/Profile...彈出選擇文件對話框,你就可以讀入邊界輪廓文件了。 寫入輪廓文件,你也可以在指定邊界或者表面的條件上創(chuàng)建輪廓文件。例如:你可以在一個算例的出口條件中創(chuàng)建一個輪廓文件,然后在其它算例中讀入該輪廓文件,并使用出口輪廓作為新算例的入口輪廓。要寫一個輪廓文件,你需要使用Write Profile面板(Figure 1),菜單:File/Write/Profile 44 在計算區(qū)域內(nèi)的某一個面(2D)或一個體(3D)內(nèi)定義體積熱源或組分質(zhì)量源。如何把這個zone定義出來?而且這個zone仍然是流體流動的。 在gambit中先將需要的zone定義出來,對于要隨流體流動我覺得這個可以用動網(wǎng)格來處理 在動網(wǎng)格設置界面 將這個隨流體流動的zone設置成剛體 這樣既可以作為zone不影響流體流通 也可以隨流體流動 只是其運動的udf不好定義 最好根據(jù)其流動規(guī)律編動網(wǎng)格udf 46 如何選擇單、雙精度解算器的選擇 Fluent的單雙精度求解器適合于所有的計算平臺,在大多數(shù)情況下,單精度求解器就能很好地滿足計算精度要求,且計算量小。 但在有些情況下推薦使用雙精度求解器: 1, 如果幾何體包含完全不同的尺度特征(如一個長而壁薄的管),用雙精度的; 2, 如果模型中存在通過小直徑管道相連的多個封閉區(qū)域,不同區(qū)域之間存在很大的壓差,用雙精度。 3, 對于有較高的熱傳導率的問題或?qū)τ谟休^大的長寬比的網(wǎng)格,用雙精度。 47 求解器為flunet5/6在設置邊界條件時,specify boundary types下的types中有三項關于interior,interface,internal設置,在什么情況下設置相應的條件?它們之間的區(qū)別是什么?interior好像是把邊界設置為內(nèi)容默認的一部分;interface是兩個不同區(qū)域的邊界區(qū),比如說離心泵的葉輪旋轉(zhuǎn)區(qū)和葉輪出口的交界面;internal;請問以上三種 每個的功能?最好能舉一兩個例子說明一下,因為這三個都是內(nèi)部條件吧,好像用的很多。interface,interior,internal boundary區(qū)別?在Fluent中,Interface意思為“交接面”,主要用途有三個:多重坐標系模型中靜態(tài)區(qū)域與運動區(qū)域之間的交接面的定義;滑移網(wǎng)格交接處的交接面定義,例如:兩車交會,轉(zhuǎn)子與定子葉柵模型,等等,在Fluent中,interface的交接重合處默認為interior,非重合處默認 為wall;非一致網(wǎng)格交接處,例如:上下網(wǎng)格網(wǎng)格間距不同等。Interior意思為“內(nèi)部的”,在Fluent中指計算區(qū)域。Internal意思為“內(nèi)部的”,比如說內(nèi)能,內(nèi)部放射率等,具體應用不太清楚。 48 FLUENT并行計算中Flexlm如何對多個License的管理?在FLEXlm LMTOOLS Utility-〉config services->service name里選好你要啟動的軟件的配備的service name,然后配置好下邊的path to the lmgrd.exe file和path to the license file,然后save service,轉(zhuǎn)到FLEXlm LMTOOLS Utility->config services-〉start/stop/reread下,選中要啟動的license,start server即可 49 在“solver”中2D 、axisymmetric和axisymmetric swirl如何區(qū)別?對于2D和3D各有什么適用范圍?從字面的意思很好理解axisymmetric和axisymmetric swirl的差別:axisymmetric:是軸對稱的意思,也就是關于一個坐標軸對稱,2D的axisymmetric問題仍為2D問題。而axisymmetric swirl:是軸對稱旋轉(zhuǎn)的意思,就是一個區(qū)域關于一條坐標軸回轉(zhuǎn)所產(chǎn)生的區(qū)域,這產(chǎn)生的將是一個回轉(zhuǎn)體,是3D的問題。在Fluent中使用這個,是將一個3D的問題簡化為2D問題,減少計算量,需要注意的是,在Fluent中,回轉(zhuǎn)軸必須是x軸。 50 在設置速度邊界條件時,提到了“Velocity formulation(Absolute和Relative)”都是指的動量方程的相對速度表示和絕對速度表示,這兩個速度如何理解? 在定義速度入口邊界條件時,Reference Frame中有Absolute和Relative to Adjacent Cell Zone的選項,關于這個,Fluent用戶手冊上是這樣寫的:“ If the cell zone adjacent to the velocity inlet is moving, you can choose to specify relative or absolute velocities by selecting Relative to Adjacent Cell Zone or Absolute in the Reference Frame drop-down list. If the adjacent cell zone is not moving, Absolute and Relative to Adjacent Cell Zone will be equivalent, so you need not visit the list. ” 如果速度入口處的單元在計算的過程中有運動發(fā)生的情況(如果你使用了運動參考系或者滑移網(wǎng)格),你可以選擇使用指定相對于鄰近單元區(qū)域的速度或在參考坐標系中的絕對速度來定于入口處的速度;如果速度入口處的相鄰單元在計算過程中沒有發(fā)生運動,那么這兩種方法所定義的速度 是等價的。 Specifying Relative or Absolute Velocity If the cell zone adjacent to the wall is moving (e.g., if you are using a moving reference frame or a sliding mesh), you can choose to specify velocities relative to the zone motion by enabling the Relative to Adjacent Cell Zone option. If you choose to specify relative velocities, a velocity of zero means that the wall is stationary in the relative frame, and therefore moving at the speed of the adjacent cell zone in the absolute frame. If you choose to specify absolute velocities (by enabling the Absolute option), a velocity of zero means that the wall is stationary in the absolute frame, and therefore moving at the speed of the adjacent cell zone--but in the opposite direction--in the relative reference frame. If you are using one or more moving reference frames, sliding meshes, or mixing planes, and you want the wall to be fixed in the moving frame, it is recommended that you specify relative velocities (the default) rather than absolute velocities. Then, if you modify the speed of the adjacent cell zone, you will not need to make any changes to the wall velocities, as you would if you specified absolute velocities. Note that if the adjacent cell zone is not moving, the absolute and relative options are equivalent. 這個問題好像問的不是特別清楚,在Fluent6.3中,問題出現(xiàn)的這個Velocity formulation(Absolute和Relative)設置,應該是設置求解器時出現(xiàn)的選項,在使用Pressure-based的求解器時,Fluent允許用戶定義的速度形式有絕對的和相對的,使用相對的速度形式是為了在Fluent中使用運動參考系以及滑移網(wǎng)格方便定義速度,關于這兩個速度的理解很簡單,可以參?忌廈嫻乃得鰨蝗綣褂肈ensity-based的求解器,這個求解器的算法只允許統(tǒng)一使用絕對的速度形式。 51 對于出口有回流的問題,在出口應該選用什么樣的邊界條件(壓力出口邊界條件、質(zhì)量出口邊界條件等)計算效果會更好?答:給定流動出口的靜壓。對于有回流的出口,壓力出口邊界條件比質(zhì)量出口邊界條件邊界條件更容易收斂。 壓力出口邊界條件壓力根據(jù)內(nèi)部流動計算結(jié)果給定。其它量都是根據(jù)內(nèi)部流動外推出邊界條件。該邊界條件可以處理出口有回流問題,合理的給定出口回流條件,有利于解決有回流出口問題的收斂困難問題。 出口回流條件需要給定:回流總溫(如果有能量方程),湍流參數(shù)(湍流計算),回流組分質(zhì)量分數(shù)(有限速率模型模擬組分輸運),混合物質(zhì)量分數(shù)及其方差(PDF 計算燃燒)。如果有回流出現(xiàn),給的表壓將視為總壓,所以不必給出回流壓力?；亓髁鲃臃较蚺c出口邊界垂直。 52 對于不同求解器,離散格式的選擇應注意哪些細節(jié)?實際計算中一階迎風差分與二階迎風差分有什么異同?離散格式對求解器性能的影響 控制方程的擴散項一般采用中心差分格式離散,而對流項則可采用多種不同的格式進行離散。Fluent允許用戶為對流項選擇不同的離散格式(注意:粘性項總是自動地使用二階精度的離散格式)。默認情況下,當使用分離式求解器時,所有方程中的對流項均用一階迎風格式離散;當使用耦合 式求解器時,流動方程使用二階精度格式,其他方程使用一階精度格式進行離散。此外,當選擇分離式求解器時,用戶還可為壓力選擇插值方式。 當流動與網(wǎng)格對齊時,如使用四邊形或六面體網(wǎng)格模擬層流流動,使用一階精度離散格式是可以接受的。但當流動斜穿網(wǎng)格線時,一階精度格式將產(chǎn)生明顯的離散誤差(數(shù)值擴散)。因此,對于2D三角形及3D四面體網(wǎng)格,注意使用二階精度格式,特別是對復雜流動更是如此。一般來講,在一 階精度格式下容易收斂,但精度較差。有時,為了加快計算速度,可先在一階精度格式下計算,然后再轉(zhuǎn)到二階精度格式下計算。如果使用二階精度格式遇到難于收斂的情況,則可考慮改換一階精度格式。 對于轉(zhuǎn)動及有旋流的計算,在使用四邊形及六面體網(wǎng)格式,具有三階精度的QUICK格式可能產(chǎn)生比二階精度更好的結(jié)果。但是,一般情況下,用二階精度就已足夠,即使使用QUICK格式,結(jié)果也不一定好。乘方格式(Power-law Scheme)一般產(chǎn)生與一階精度格式相同精度的結(jié)果。中心差分格式一般只用于大渦模擬,而且要求網(wǎng)格很細的情況。 53 對于FLUENT的耦合解算器,對時間步進格式的主要控制是Courant數(shù)(CFL),那么Courant數(shù)對計算結(jié)果有何影響?courant number實際上是指時間步長和空間步長的相對關系,系統(tǒng)自動減小courant數(shù),這種情況一般出現(xiàn)在存在尖銳外形的計算域,當局部的流速過大或者壓差過大時出錯,把局部的網(wǎng)格加密再試一下。 在Fluent中,用courant number來調(diào)節(jié)計算的穩(wěn)定性與收斂性。一般來說,隨著courant number的從小到大的變化,收斂速度逐漸加快,但是穩(wěn)定性逐漸降低。所以具體的問題,在計算的過程中,最好是把courant number從小開始設置,看看迭代殘差的收斂情況,如果收斂速度較慢而且比較穩(wěn)定的話,可以適當?shù)脑黾觕ourant number的大小,根據(jù)自己具體的問題,找出一個比較合適的courant number,讓收斂速度能夠足夠的快,而且能夠保持它的穩(wěn)定性。 54 在分離求解器中,FLUENT提供了壓力速度耦和的三種方法:SIMPLE,SIMPLEC及PISO,它們的應用有什么不同在FLUENT中,可以使用標準SIMPLE算法和SIMPLEC(SIMPLE-Consistent)算法,默認是SIMPLE算法,但是對于許多問題如果使用SIMPLEC可能會得到更好的結(jié)果,尤其是可以應用增加 的亞松馳迭代時,具體介紹如下: 對于相對簡單的問題(如:沒有附加模型激活的層流流動),其收斂性已經(jīng)被壓力速度耦合所限制,你通?？梢杂肧IMPLEC算法很快得到收斂解。在SIMPLEC中,壓力校正亞松馳因子通常設為1.0,它有助于收斂。但是,在有些問題中,將壓力校正松弛因子增加到1.0可能會導致不穩(wěn)定。 對于所有的過渡流動計算,強烈推薦使用PISO算法鄰近校正。它允許你使用大的時間步,而且對于動量和壓力都可以使用亞松馳因子1.0。對于定常狀態(tài)問題,具有鄰近校正的PISO并不會比具有較好的亞松馳因子的SIMPLE或SIMPLEC好。對于具有較大扭曲網(wǎng)格上的定常狀態(tài)和過渡計算推薦使 用PISO傾斜校正。當你使用PISO鄰近校正時,對所有方程都推薦使用亞松馳因子為1.0或者接近1.0。如果你只對高度扭曲的網(wǎng)格使用PISO傾斜校正,請設定動量和壓力的亞松馳因子之和為1.0比如:壓力亞松馳因子0.3,動量亞松馳因子0.7)。如果你同時使用PISO的兩種校正方法,推薦參?腜ISO鄰近校正中所用的方法 55 對于大多數(shù)情況,在選擇選擇壓力插值格式時,標準格式已經(jīng)足夠了,但是對于特定的某些模型使用其它格式有什么特別的要求?壓力插值方式的列表只在使用Pressure-based求解器中出現(xiàn)。一般情況下可選擇Standard;對于含有高回旋數(shù)的流動,高Rayleigh數(shù)的自然對流,高速旋轉(zhuǎn)流動 ,多孔介質(zhì)流動,高曲率計算區(qū)域等流動情況,選擇PRESTO格式;對于可壓縮流動,選擇Second Order;當然也可以選擇Second Order以提高精度;對于含有大體力的流動,選擇Body Force Weighted。注意:Second Order格式不可以用于多孔介質(zhì);在使用VOF和Mixture多相流模型時,只能使用PRESTO或Body Force Weighted格式。關于壓力插值格式的詳細內(nèi)容,請參考Fluent用戶手冊 57 討論在數(shù)值模擬過程中采用四面體網(wǎng)格計算效果好,還是采用六面體網(wǎng)格更妙呢?在2D中,FLUENT 可以使用三角形和四邊形單元以及它們的混合單元所構(gòu)成的網(wǎng)格。在3D中,它可以使用四面體,六面體,棱錐,和楔形單元所構(gòu)成的網(wǎng)格。選擇那種類型的單元取決于你的應用。當選擇網(wǎng)格類型的時候,應當考慮以下問題: 設置時間(setup time) 計算成本(computational expense) 數(shù)值耗散(numerical diffusion ) 1.設置時間 在工程實踐中,許多流動問題都涉及到比較復雜的幾何形狀。一般來說,對于這樣的問題,建立結(jié)構(gòu)或多塊(是由四邊形或六面體元素組成的)網(wǎng)格是極其耗費時間的。所以對于復雜幾何形狀的問題,設置網(wǎng)格的時間是使用三角形或四面體單元的非結(jié)構(gòu)網(wǎng)格的主要動機。然而,如果所使用 的幾何相對比較簡單,那么使用哪種網(wǎng)格在設置時間方面可能不會有明顯的節(jié)省。 如果你已經(jīng)有了一個建立好的結(jié)構(gòu)代碼的網(wǎng)格,例如FLUENT 4,很明顯,在FLUENT中使用這個網(wǎng)格比重新再生成一個網(wǎng)格要節(jié)省時間。這也許是你在FLUENT 模擬中使用四邊形或六面體單元的一個非常強的動機。注意,對于從其它代碼導入結(jié)構(gòu)網(wǎng)格,包括FLUENT 4,FLUENT 有一個篩選的范圍。 2.計算成本 當幾何比較復雜或流程的長度尺度的范圍比較大的時候,可以創(chuàng)建是一個三角形/四面體網(wǎng)格,因為它與由四邊形/六面體元素所組成的且與之等價的網(wǎng)格比較起來,單元要少的多。這是因為一個三角形/ 四面體網(wǎng)格允許單元群集在被選擇的流動區(qū)域中,而結(jié)構(gòu)四邊形/六面體網(wǎng)格一般會把單元強加到所不需要的區(qū)域中。對于中等復雜幾何,非結(jié)構(gòu)四邊形/六面體網(wǎng)格能構(gòu)提供許多三角形/ 四面體網(wǎng)格所能提供的優(yōu)越條件。 在一些情形下使用四邊形/六面體元素是比較經(jīng)濟的,四邊形/六面體元素的一個特點是它們允許一個比三角形/四面體單元大的多的縱橫比。一個三角形/ 四面體單元中的一個大的縱橫比總是會影響單元的偏斜(skewness),而這不是所希望的,因為它可能妨礙計算的精確與收斂。所以,如果你有一個相對簡單的幾何,在這個幾何中流動與幾何形狀吻合的很好,例如一個瘦長管道,你可以運用一個高縱橫比的四邊形/六面體單元的網(wǎng)格。這?鐾裼滌械牡ピ贍鼙熱切? 四面體少的多。 3.數(shù)值耗散 在多維情形中,一個錯誤的主要來源是數(shù)值耗散,術語也為偽耗散(false diffusion)。之所以稱為“偽耗散”是因為耗散不是一個真實現(xiàn)象,而是它對一個流動計算的影響近似于增加真實耗散系數(shù)的影響。 關于數(shù)值耗散的觀點有: 當真實耗散小,即情形出現(xiàn)對流受控時(即本身物理耗散比較小時),數(shù)值的耗散是最值得注意的。 關于流體流動的所有實際的數(shù)值設計包括有限數(shù)量的數(shù)值耗散。這是因為數(shù)值耗散起于切斷錯誤,而切斷錯誤是一個表達離散形式的流體流動方程的結(jié)果。 用于FLUENT 中的二階離散方案有助于減小數(shù)值耗散對解的影響。 數(shù)值耗散的總數(shù)反過來與網(wǎng)格的分解有關。因此,處理數(shù)值耗散的一個方法是改進網(wǎng)格。 當流動與網(wǎng)格相吻一致時,數(shù)值耗散減到最小。 最后這一點與網(wǎng)格的選擇非常有關。很明顯,如果你選擇一個三角形/ 四面體網(wǎng)格,那么流動與網(wǎng)格總不能一致。另一方面,如果你使用一個四邊形/六面體網(wǎng)格,這種情況也可能會發(fā)生,但對于復雜的流動則不會。在一個簡單流動中,例如過一長管道的流動,你可以依靠一個四邊形/六面體網(wǎng)格以盡可能的降低數(shù)值的耗散。在這種情形,使用一個四邊形/六?嫣逋窨贍苡行┯欣跫蛭朧褂靡桓鋈切? 四面體單元比起來,你將能夠使用比較少的單元而得到一個更好的解。 59 在UDF中compiled型的執(zhí)行方式和interpreted型的執(zhí)行方式有什么不同編譯型UDF: 采用與FLUENT 本身執(zhí)行命令相同的方式構(gòu)建的。采用一個稱為Makefile的腳本來引導c 編譯器構(gòu)造一個當?shù)啬繕司幋a庫(目標編碼庫包含有將高級c 語言源代碼轉(zhuǎn)換為機器語言。)這個共享庫在運行時通過“動態(tài)加載”過程載入到FLUENT 中。目標庫特指那些使用的計算機體系結(jié)構(gòu),和運行的特殊FLUENT 版本。因此,FLUENT 版本升級,計算機操作系統(tǒng)改變以及在另一臺不同類型的計算機上運行時,這個庫必須進行重構(gòu)。 編譯型UDF 通過用戶界面將原代碼進行編譯,分為兩個過程。這兩個過程是:訪問編譯UDF 面板,從源文件第一次構(gòu)建共享庫的目標文件中;然后加載共享庫到FLUENT 中。 采用與FLUENT 本身執(zhí)行命令相同的方式構(gòu)建的。采用一個稱為Makefile的腳本來引導c 編譯器構(gòu)造一個當?shù)啬繕司幋a庫(目標編碼庫包含有將高級c 語言源代碼轉(zhuǎn)換為機器語言。)這個共享庫在運行時通過“動態(tài)加載”過程載入到FLUENT 中。目標庫特指那些使用的計算機體系結(jié)構(gòu),和運行的特殊FLUENT 版本。因此,FLUENT 版本升級,計算機操作系統(tǒng)改變以及在另一臺不同類型的計算機上運行時,這個庫必須進行重構(gòu)。 編譯型UDF 通過用戶界面將原代碼進行編譯,分為兩個過程。這兩個過程是:訪問編譯UDF 面板,從源文件第一次構(gòu)建共享庫的目標文件中;然后加載共享庫到FLUENT 中。 解釋型UDF: 解釋型UDF 同樣也是通過圖形用戶界面解釋原代碼,卻只有單一過程。這一過程伴隨著運行,包含對解釋型UDF 面板的訪問,這一面板位于源文件中的解釋函數(shù)。 在FLUENT內(nèi)部,源代碼通過c 編譯器被編譯為即時的、體系結(jié)構(gòu)獨立的機器語言。UDF 調(diào)用時,機器編碼通過內(nèi)部模擬器或者解釋器執(zhí)行。額外層次的代碼導致操作不利,但是允許解釋型UDF 在不同計算結(jié)構(gòu),操作系統(tǒng)和FLUENT 版本上很容易實現(xiàn)共享。如果迭代速度成為焦點時,解釋型UDF 可以不用修改就用編譯編碼直接運行。 解釋型UDF 使用的解釋器不需要有標準的c 編譯器的所有功能。特別是解釋型UDF 不含有下列C 程序語言部分: goto 語句聲明;無ANSI-C 語法原形;沒有直接數(shù)據(jù)結(jié)構(gòu)引用;局部結(jié)構(gòu)的聲明;聯(lián)合函數(shù)指針;函數(shù)陣列; 解釋型UDF與編譯型UDF的區(qū)別: 在解釋型與編譯型UDF 之間的主要的不同之處是很重要的,例如當你想在UDF 中引進新的數(shù)據(jù)結(jié)構(gòu)時。解釋型不能通過直接數(shù)據(jù)引用獲得FLUENT 解算器的數(shù)據(jù);只能間接的通過FLUENT 預先提供的宏來獲取數(shù)據(jù)。具體請參考第7 章。 在解釋型與編譯型UDF 之間的主要的不同之處是很重要的,例如當你想在UDF 中引進新的數(shù)據(jù)結(jié)構(gòu)時。解釋型不能通過直接數(shù)據(jù)引用獲得FLUENT 解算器的數(shù)據(jù);只能間接的通過FLUENT 預先提供的宏來獲取數(shù)據(jù)。具體請參考第7 章。 總結(jié)一下,當選擇寫解釋型或者編譯型UDF時,記住以下幾條: 解釋型UDF:對別的運行系統(tǒng)是可移植的,可以作為編譯型運行,不需要c 編譯器,比編譯型的要慢,在使用C 程序語言時有限制,不能鏈接到編譯系統(tǒng)或者用戶庫,只能通過預先提供的宏訪問FLUENT 中存儲的數(shù)據(jù)。 編譯型UDF:運行要快于解釋型UDF,對C 程序語言沒有限制,可以使用任何ANSI-compliant c 編譯器進行編譯,可以調(diào)用其他語言寫的函數(shù)(特別是獨立于系統(tǒng)和編譯器的),如果包含某些解釋器不能處理的c 語言部分時用解釋型UDF 是不行的。 總之,當決定哪一類型的udf 應用到你的模型時: 對小的,直接的函數(shù)用解釋型;對復雜函數(shù)使用編譯型61 FLUENT help和GAMBIT help能教會我們(特別是剛?cè)腴T的新手)學習什么基本知識?答:可以了解其基本原理和基本的操作。不過我覺得對于新手熟悉軟件最好的還是tutorial guide63 FLUENT模擬飛行器外部流場,最高MA多少時就不準確了?MA達到一定的程度做模擬需注意哪些問題?答:不準確的標準是什么?沒有判斷標準就沒辦法判斷。一般來說fluent計算馬赫數(shù)大于3~5之后就不是很理想了(不過相信版本越新結(jié)果越好)。計算的時候應該從低馬赫數(shù)慢慢往上算。?熱縊等綣撲懵硨帳?的話,就在馬赫數(shù)4的計算結(jié)果上算。另外,求解器需選擇耦合和顯式的。(對于6.3來說,選擇基于密度的求解器) 68 做飛機設計時,經(jīng)常計算一些翼型,可是經(jīng)常出現(xiàn)計算出來的阻力是負值,出現(xiàn)負值究竟是什么原因,是網(wǎng)格的問題還是計算參數(shù)設置的問題?如果這個問題對于某個人經(jīng)常出現(xiàn)的話,那就比較奇怪了,阻力是負值,難道就是傳說中的前緣吸力現(xiàn)身?呵呵,只是開個玩笑:),估計肯定是 計算錯了或者是設置錯了。在飛機翼型氣動里面,阻力主要有兩種成份:壓差阻力和摩擦阻力。應該是正值的。排除是計算過程的其他問題,我覺得在使用Fluent進行這方面的計算時,需要注意兩個方面:1.參考值的設置,也就是Report->Reference Values...這些參考值,是用來計算Re,以及升力,阻力,力矩系數(shù)所要用到的。如果設置不當,即使計算過程是對的,所得到的升阻力等系數(shù)也是不對的。對于2D翼型仿真計算,比較容易出錯的就是里面的Area該寫什么,單位是平方米,這里應該填寫翼型的弦長(Chord Length),The area here is actually area per unit depth;就是每單位展長的面積。2.在監(jiān)視力的時候,關于力的矢量方向設置,Solve->Monitor->Force...這個矢量方向千萬不要小看,不能填錯,填錯了就可能出現(xiàn)阻力是負值的錯誤,Fluent之前的版本所附帶的例子,關于NACA0012翼型的計算中,這里的矢量就設置錯了,受錯誤例子的?跋歟賈夷潛臼櫓腥切我硇偷哪歉隼右采柚么砦螅謔櫚牡?12頁的第6步的第(7)小步就設置錯誤,升力系數(shù)的力方向矢量,應該是X=-0.087155,Y=0.996195;前面他也寫到要注意:要確保阻力和升力分別與來流平行和垂直,那么這兩個力矢量肯定是垂直的了,那么這兩個矢量?牡慍司涂隙ǖ扔諏懔耍凰業(yè)氖牽贔luent6.3版本的例子中,這個錯誤已經(jīng)改正過來了。 74 大概需要劃分100萬個左右的單元,且只計算穩(wěn)態(tài)流動,請問這樣的問題PC機上算的了嗎?如果能算至少需要怎樣的計算機配置呢?答:一般來說,按照1000個節(jié)點對1MB內(nèi)存這樣預估就差不多了,只計算穩(wěn)態(tài)流動,pc機應該差不多了,不過因為一般的pc機可能在連續(xù)計算5、6天之后就出現(xiàn) 浮點運算錯誤,所以如果計算不是很復雜,采用的求解器和湍流模型不是太好計算資源,應該還是可以的。如果使用pc機計算,建議至少采用2GB內(nèi)存,主板最好固態(tài)電容,不易爆漿,電源最好功率大典,應該差不多了,現(xiàn)在流行四核cpu的,可以考慮使用四核的,這樣的配置下來也不比服 務器差多少。 76 GAMBIT劃分三維網(wǎng)格后,怎樣知道結(jié)點數(shù)?如何知道總生成多少網(wǎng)格(整個模型)? 答:個人一般是將網(wǎng)格讀入fluent后,通過grid->info->size來看:77 在FLUENT的后處理中可以顯示一個管道的。某個標量的。圓截面平均值沿管道軸線(中心線)的變化曲線嗎?何顯示空間某一點的數(shù)值呀(比如某一點溫度)? 正確的方法應該是輸入命令畫曲線命令輸入狀態(tài)下直接按回車>plot>c-a-a (就是circum-average-axial)再空按回車顯示可以選擇的值(從溫度到nusselt數(shù)應有盡有)比如輸入>temp (溫度)>100 (軸向數(shù)據(jù)點個數(shù))>filename.txt (文件名,隨便取)>no (不知道什么,order point)然后在plot-file里選擇輸出就可以了另定義空間點的方法為surface-point,輸入點的坐標或者直接在網(wǎng)格上標記,然后就能在后處理時看到這個點的選項了。80 如何在gambit中輸入cad和Pro/e的圖形?如何將FLUNET的結(jié)果EXPORT成ANSYS的文件?答:autocad需要將圖形轉(zhuǎn)化為sat格式,pro/e可以將文件轉(zhuǎn)化為igse或者stp格式。在fluent的flie/export 中可以選擇導出ansys格式的文件 87 courant數(shù):在模擬高壓的流場的時候,迭代的時候總是自動減小其數(shù)值,這是什么原因造成的,為什么?怎么修改?這是流場的壓力梯度較大,Fluent自身逐步降低時間步長,防止計算發(fā)散。我一般的處理辦法是:先將邊界條件上的壓力設置較低點,使得壓力梯度較小一點,等到收斂的感 覺差不多,在這個基礎上,逐漸把壓力增大,這樣就不容易發(fā)散。94 把帶網(wǎng)格的幾個volume,copy到另一處,但原來split的界面,現(xiàn)在都變成了wall,怎么才能把wall變成內(nèi)部流體呢?直接邊界面定義為interior即可 97 在udf中,U,V,W代表的速度,分別代表什么方向的,直角坐標還是柱坐標? 應該是在直角坐標系中,U、V、W分別代表是X、Y、Z的方向98 Gambit的網(wǎng)格相連問題:如果物體是由兩個相連的模型所結(jié)合,一個的網(wǎng)格劃分比較密、另一個比較稀疏,用Gambit有辦法將兩個網(wǎng)格密度不同的物體,相連在一起嗎? 請參考第16題答案。將兩種網(wǎng)格交界的地方設置成一對interface即可 100 在FLUENT里定義流體的密度時,定義為不可壓理想流體是用在什么地方呀,講義上說是用于可變密度的不可壓流動,不知如何理解?答:define/matirial 中定義?？勺兠芏鹊牟豢蓧嚎s流動,就是說在該流動下,流體介質(zhì)的密度可以認為不變。比如說空氣在流速在0.3馬赫的情況下都可以認為是密度不變的101 已經(jīng)建好的模型,想修改一些尺寸,但不知道頂點的座標,請問如何在gambit中顯示點的座標?答:在gambit中的geometr-〉vortex->summarize vortices即可顯示點的坐標。102 在FLUENT模擬以后用display下的操作都無法顯示,不過剛開始用的是好的,然后就不行了,為什么?答: DirectX 控制面板中的“加速”功能禁用即可 103 能否同時設置進口和出口都為壓力的邊界條件?在這樣的邊界條件設置情況下發(fā)現(xiàn)沒有收斂,研究的物理模型只是知道進口和出口的壓力,不知道怎么修改才能使其收斂? 當然可以同時設置進口和出口都為壓力的邊界條件。如果沒有收斂,需要首先看看求解器、湍流模型、氣體性質(zhì)和邊界條件時有沒有出現(xiàn)warning;其次,還是我上邊的帖子所說的,對于可壓流動,采用壓力邊界條件,不能一下把壓力和溫度加到所需值,應該首先設置較低的壓力或溫度,?緩籩鸞ピ齟?，最后达惦A約核璧鬧怠? 104 在FLUENT計算時,有時候計算時間會特別長,為了避免斷電或其它情況影響計算,應設置自動保存功能,如何設置自動保存功能?在非定常計算中讀入自動保存文件時如下出現(xiàn)問題: Writing "F:\propane\16\160575.cas"... Error: sopenoutputfile&: unable to open file for output Error Object: "F:\propane\16\160575.cas" Error: Error writing "F:\propane\16\160575.cas". Error Object: #f 非定常的,算了一段之后停下來,改天繼續(xù)算的時候,自動保存的時候出現(xiàn)問題,請問如何解決?答:File->write->Autosave就可以實現(xiàn)自動保存。只要你在寫自動保存文件的時候,文件名另取一個就行,比如Writing "F:\propane\16\160575_1.cas 105 gambit劃分時運動部分與靜止部分交接面:一個系統(tǒng)的兩塊,運動部分與靜止部分交接部分近似認為沒有空隙(無限小,雖然實際上是不可能的),假設考慮做成一個實體,那么似乎要一起運動或靜止;假設分開做成兩個實體,那么交接處的兩個不完全重合的面要設為WALL還是什么呢,設 成WALL不就不能過流了嗎?將這一對接觸面設置成Interface就行了,具體請參考第47題的解答 106 在計算模擬中,continuity總不收斂,除了加密網(wǎng)格,還有別的辦法嗎?別的條件都已經(jīng)收斂了,就差它自己了,還有收斂的標準是什么?是不是到了一定的尺度就能收斂了,比如10-e5具體的數(shù)量級就收斂了continuity 是質(zhì)量殘差,具體是表示本次計算結(jié)果與上次計算結(jié)果的差別,如果別的條件收斂了,就差它?？梢渣creport, 打開里面FLUX選項,算出進口與出口的質(zhì)量流量差,看它是否小于0.5%.如果小于,可以判斷它收斂. 108 想把gambit的圖形保存成圖片,可是底色總是黑色,怎么改為白色呀。用windows中畫圖板的反色,好像失真很多。如何處理?答:首先點開GAMBIT的EDIT菜單,其次點GRAPHICS,在下拉列表中點到 WINDOWS BACKGROUND COLOR BLACK 一項 在下面VALUE 中填寫WHITE,再點左面的MODIFY,就可以了.110 在分析一個轉(zhuǎn)輪時,想求得轉(zhuǎn)輪的轉(zhuǎn)矩,不知道fluent中有什么方法可以提供該數(shù)據(jù)。本來想到用葉片上面的壓力乘半徑,然后做積分運算,但是由于葉片正反壁面統(tǒng)一定義的,即全部定義為wall-rn1,所以分不出方向來了答:report/force/moment 定義需要計算的面和旋轉(zhuǎn)中心就ok了111 如何在gambit中實現(xiàn)坐標軸的變換:有一個三維的網(wǎng)格,想在柱坐標中實現(xiàn),可是gambit中一直顯示直角坐標?答:對于這個問題,你可以嘗試一下:Operation--->Tools--->Coordinate System--->Activate Coordinate System 113 利用vof非穩(wěn)態(tài)求解,結(jié)果明顯沒有收斂的情況下,為什么就開始提示收斂,雖然可以不管它,繼續(xù)算下去達到收斂。但是求解怎么會提前收斂?可以吧殘差圖的k和E 改小點,就好點了。另外 vof 中 殘差圖一般是波動的吧所以最好 設置檢測面 比如說進出口流量 來確定你的計算的準確性116 在Gambit中如何將兩個dbs文件到入:把爐膛分成了三個dbs文件,現(xiàn)在想導入兩個dbs文件,在Gambit中進行操作,但好象使用open命令就只能open一個dbs文件,請問這要怎么處理?答:將其中一個導出成iges或者別的格式,然后就能和dbs一起導入了。 119 用GAMBIT生成網(wǎng)格時要是出現(xiàn)負值怎么辦啊?有什么辦法可以改正嗎,只能將網(wǎng)格重新畫嗎?答:好像只能重新生成。也就是要在線上重新布點,重新生成面網(wǎng)格等等。120 scale是把你所畫模型中的單位轉(zhuǎn)化為Fluent默認的m,而unite是根據(jù)你自己的需要轉(zhuǎn)化單位,也就是把Fluent中默認的m轉(zhuǎn)畫為其他的單位,兩中方法對計算沒有什么影響嗎?答:scale是對幾何進行比例縮放,而unit只是改變單位,不改變幾何外形的大小。比如,一個是1m的幾何外形,?ü齭cale將m變?yōu)閙m,那么幾何外形就變成了1mm。如果通過unite將m改為mm,那么幾何外形不變,還是1000mm,只是表示的單位變成mm了 121 GAMBIT處理技巧:兩個圓內(nèi)切產(chǎn)生的尖角那個面如何生成網(wǎng)格質(zhì)量才比較好?答:可以采用劃分結(jié)構(gòu)網(wǎng)格的方法(對于狹縫的一般處理都是生成長寬比很大的結(jié)構(gòu)網(wǎng)格);或者將這個尖角導個圓弧之后再劃分網(wǎng)格(也就是進行幾何簡化)128 在gambit中對一體積成功的進行了體網(wǎng)格,網(wǎng)格進行了examine mesh,也沒有什么問題,可當要進行邊界類型(boundary type)的設定時,卻發(fā)現(xiàn)type 只有node, element_side兩項,沒有什么wall,pressure_outlet等。為何無法定義邊界?答:因為沒有選擇求解器為fluent 5/6132 網(wǎng)格數(shù)量和內(nèi)存之間的關系是什么?大概有這樣一個估計:“1k網(wǎng)格=1M內(nèi)存。”對于一臺有1G內(nèi)存的計算機,你能接受的計算網(wǎng)格數(shù)最好少于100萬,當然這只是一個粗略的說法,影響計算速度的因素還有Fluent計算的設置等。 133 如何在FLUENT中進行密度的選擇?FLUENT中的密度我個人認為是指流體的狀態(tài)方程,包括不可壓縮(常數(shù)),不可壓縮理想流體;可壓縮理想流體;或者定義成溫度的線性關系;也可以通過udf定義成壓力和溫度的復雜函數(shù)(但是我親自實踐過)定義成壓力的函數(shù)要么收斂后結(jié)果不對,要?錘靜皇樟玻ㄎ葉ㄒ宓氖莍f97標準水蒸汽性質(zhì))。因此,如果說問題中的密度要寫成壓力的復雜函數(shù),會對收斂帶來一定困難。設置,在所選材料的density里面選 142 什么是多孔介質(zhì);在那些方面應用?多孔介質(zhì)由固體物質(zhì)組成的骨架和由骨架分隔成大量密集成群的微小空隙構(gòu)成的介質(zhì)。 多孔介質(zhì)在很多方面都有運用,像醫(yī)學,工業(yè),環(huán)境,軍事等! 我個人曾經(jīng)做過這方面的實驗。主要是多孔陶瓷,運用它來進行除塵的! 143 有沒有介紹多孔介質(zhì)的專業(yè)書籍? 《多孔介質(zhì)流體動力學》貝爾 著;李競生,陳崇希 譯 北京:中國建筑工業(yè)出版社,1983.8《 多孔介質(zhì)污染物遷移動力學》仵彥卿編著上海:上海交通大學出版社,2007《多孔介質(zhì)——流體滲移與孔隙結(jié)構(gòu)》 FAL DULLIEN 著 石油工業(yè)出版 《多孔介質(zhì)傳熱傳質(zhì)理論與應用》作者:劉偉 范愛武 黃曉明 科學出版社《多孔介質(zhì)溶質(zhì)運移動力學 》張永祥,陳鴻漢著 北京 地震出版社 2000 155 如何區(qū)分層流和紊流? 以什么為標準來區(qū)分呢?從層流過渡到紊流的標準是什么? 答:自然界中的流體流動狀態(tài)主要有兩種形式,即層流laminar和湍流(就是問題中所說的紊流)turbulence.層流是指流體在流動過程中兩層之間沒有相互混滲,而湍流是指流體不是處于分層流動狀態(tài)。對于圓管內(nèi)流動,雷諾數(shù)小于等于2300,管流一定為層流,雷諾數(shù)大于等于8000到12000??，管流覛gㄎ牧鰨著凳笥?300而小于8000時,流動處于層流與湍流的過渡區(qū)。對于一般流動,在計算雷諾數(shù)時,可以用水力半徑代替管徑。 157 gambit不響應鍵盤的原因。要把輸入法調(diào)到英文輸入狀態(tài),, 158 gambit不能正常啟動的原因有哪些1、在Gambit建模過程中出現(xiàn)界面突然跳出,并且下次運行Gambit時,界面調(diào)不出來,這時只需刪去gambit工作目錄下的(默認的工作目錄為\FLUENT.INC\ntbin\ntx86)后綴為*.lok的文件,就會恢復正常。 2、出錯信息“IDENTIFIER "default_id" CURRENTLY OPEN”,Gambit的缺省文件已經(jīng)打開,gambit運行失敗,到用戶默認目錄刪default_id.*等文件。 3、出錯信息“unable find Exceed X Server ” ,GAMBIT需要裝EXCEED才能用,推薦EXCEED 6.2。 159 在fluent中如何設置工作目錄?在Gambit中如何設置工作目錄?找到桌面上的Fluent或者Gambit圖標,右鍵圖標,“屬性”->"起始位置”...將起始位置設置為你想要的文件夾目錄就可以了。 160 在計算過程中其它指數(shù)都收斂了,就continuity不收斂是怎么回事?在初始化設置中,那些項影響continuity的收斂?在計算過程中其他指數(shù)都收斂了,就continuity不收斂,這種情況一般出現(xiàn)在多相流中,在初始化設置中,可能把上次計算結(jié)果的進口參數(shù)作為初始化設置,可以加快conti nuity的收斂,不過更重要的是改進網(wǎng)格質(zhì)量

相關標簽搜索：CFD常見問題(轉(zhuǎn)載） CFD培訓 CFD流體分析培訓 cfd視頻 fluent cfx pumplinx軟件培訓 Fluent、CFX流體分析 HFSS電磁分析 Ansys培訓 Abaqus培訓 Autoform培訓 有限元培訓