fluent收斂終極方案_第1頁
fluent收斂終極方案_第2頁
fluent收斂終極方案_第3頁
fluent收斂終極方案_第4頁
fluent收斂終極方案_第5頁
已閱讀5頁,還剩4頁未讀 繼續免費閱讀

下載本文檔

版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領

文檔簡介

1、Continuity收斂問題 FLUENT求解器設置FLUENT求解器設置主要包括:1、壓力-速度耦合方程格式選擇2、對流插值 3、梯度插值 4、壓力插值下面對這幾種設置做詳細說明。一、壓力-速度耦合方程求解算法FLUENT中主要有四種算法:SIMPLE,SIMPLEC,PISO,FSM(1)SIMPLE(semi-implicit method for pressure-linked equations)半隱式連接壓力方程方法,是FLUENT的默認格式。(2)SIMPLEC(SIMPLE-consistent)。對于簡單的問題收斂非常快速,不對壓力進行修正,所以壓力松弛因子可以設置為1(3)

2、Pressure-Implicit with Splitting of Operators (PISO)。對非定常流動問題或者包含比平均網格傾斜度更高的網格適用(4)Fractional Step Method (FSM)對非定常流的分步方法。用于NITA格式,與PISO具有相同的特性。二、對流插值(動量方程)FLUENT有五種方法:一階迎風格式、冪率格式、二階迎風格式、MUSL三階格式、QUICK格式(1)FLUENT默認采用一階格式。容易收斂,但精度較差,主要用于初值計算。(2)Power Lar.冪率格式,當雷諾數低于5時,計算精度比一階格式要高。(3)二階迎風格式。二階迎風格式相對于一

3、階格式來說,使用更小的截斷誤差,適用于三角形、四面體網格或流動與網格不在同一直線上;二階格式收斂可能比較慢。(4)MUSL(monotone upstream-centered schemes for conservation laws).當地3階離散格式。主要用于非結構網格,在預測二次流,漩渦,力等時更精確。(5)QUICK(Quadratic upwind interpolation)格式。此格式用于四邊形/六面體時具有三階精度,用于雜交網格或三角形/四面體時只具有二階精度。三、梯度插值梯度插值主要是針對擴散項。FLUENT有三種梯度插值方案:green-gauss cell-based,

4、Green-gauss node-based,least-quares cell based.(1)格林-高斯基于單元體。求解方法可能會出現偽擴散。(2)格林-高斯基于節點。 求解更精確,最小化偽擴散,推薦用于三角形網格上(3)基于單元體的最小二乘法插值。推薦用于多面體網格,與基于節點的格林-高斯格式具有相同的精度和格式。四、壓力插值壓力基分離求解器主要有五種壓力插值算法。(1)標準格式(Standard)。為FLUENT缺省格式,對大表妹邊界層附近的曲線發現壓力梯度流動求解精度會降低(但不能用于流動中壓力急劇變化的地方此時應該使用PRESTO!格式代替)(2)PRESTO!主要用于高旋流,壓

5、力急劇變化流(如多孔介質、風扇模型等),或劇烈彎曲的區域。(3)Linear(線性格式)。當其他選項導致收斂困難或出現非物理解時使用此格式。(4)second order(二階格式)。用于可壓縮流動,不能用于多孔介質、階躍、風扇、VOF/MIXTURE多相流。(5)Body Force Weighted體積力。當體積力很大時,如高雷諾數自然對流或高回旋流動中采用此格式。(1)連續性方程不收斂是怎么回事? 在計算過程中其它指數都收斂了,就continuity不收斂是怎么回事。 這和fluent程序的求解方法SIMPLE有關。SIMPLE根據連續方程推導出壓力修正方法求解壓力。由于連續方程中流場耦

6、合項被過渡簡化,使得壓力修正方程不能準確反映流場的變化,從而導致該方程收斂緩慢。 你可以試驗SIMPLEC方法,應該會收斂快些。 在計算模擬中,continuity總不收斂,除了加密網格,還有別的辦法嗎?別的條件都已經收斂了,就差它自己了,還有收斂的標準是什么?是不是到了一定的尺度就能收斂了,比如10e5具體的數量級就收斂了continuity 是質量殘差,具體是表示本次計算結果與上次計算結果的差別,如果別的條件收斂了,就差它。可以點report,打開里面FLUX選項,算出進口與出口的質量流量差,看它是否小于0.5%.如果小于,可以判斷它收斂. (2) fluent殘差曲線圖中continui

7、ty是什么含義?是質量守恒方程的反映,也就是連續性的殘差。這個收斂的快并不能說明你的計算就一定正確,還要看動量方程的迭代計算。表示某次迭代與上一次迭代在所有cells積分的差值,continuty表示連續性方程的殘差(3)要加速continuity收斂該設置那些參數? 感覺需要調整courant number FLUENT 中courant number是在耦合求解的時候才出現的。正確的調整,可以更好地加速收斂和解的增強穩定性。 courant number 實際上是指時間步長和空間步長的相對關系,系統自動減小courant 數,這種情況一般出現在存在尖銳外形的計算域,當局部的流速過大或者壓差

8、過大時出錯,把局部的網格加密再試一下。 在fluent 中,用courant number 來調節計算的穩定性與收斂性。一般來說,隨著courantnumber 的從小到大的變化,收斂速度逐漸加快,但是穩定性逐漸降低。所以具體的問題,在計算的過程中,最好是把ourant number 從小開始設置,看看迭代殘差的收斂情況,如果收斂速度較慢而且比較穩定的話,可以適當的增加courant number 的大小,根據自己具體的問題,找出一個比較合適的courant number,讓收斂速度能夠足夠的快,而且能夠保持它的穩定性。個人認為這應該和你采用的算法有關SIMPLE算法是根據連續方程推導出壓力修

9、正方法求解壓力。由于連續方程中流場耦合項被過渡簡化,使得壓力修正方程不能準確反映流場的變化,從而導致該方程收斂緩慢。試著用SIMPLEC算法看看。(4) 正在學習Fluent,模擬圓管內的流動,速度入口,出口outflow運行后xy的速度很快就到1e-06了,但是continuity老是降不下去,維持在1e-00和1e-03之間,減小松弛因子好像也沒什么變化大家有什么建議嗎? 你查看了流量是否平衡嗎?在report->flux里面操作,mass flow rate,把所有進出口都選上,compute一下,看看nut flux是什么水平,如果它的值小于總進口流量的1%,并且其他檢測量在繼續

10、迭代之后不會發生波動,也可以認為你的解是收斂的。 造成連續方程高殘差不收斂的原因主要有以下幾點:1.網格質量,主要可能是相鄰單元的尺寸大小相差較大,它們的尺寸之比最好控制在1.2以內,不能超過1.4.2.離散格式及壓力速度耦合方法,如果是結構網格,建議使用高階格式,如2階迎風格式等,如果是非結構網格,除pressure保持standard格式不變外,其他格式改用高階格式;壓力速度耦合關系,如果使用SIMPLE,SIMPLEC,PISO等segerated solver對聯系方程收斂沒有提高的話,可以嘗試使用coupled solver。另外,對于梯度的計算,不論使用結構或非結構網格,都可以改用

11、node-based來提高計算精度。一些情況:1.監測流場某個變量來判斷收斂更合理一些.2.網格質量.3.Velocity inlet boundary conditions are not appropriate for compressible flowproblems. (4)要加速continuity收斂該設置那些參數? 感覺需要調整courant number FLUENT 中courant number是在耦合求解的時候才出現的。正確的調整,可以更好地加速收斂和解的增強穩定性。 courant number 實際上是指時間步長和空間步長的相對關系,系統自動減小courant 數,這

12、種情況一般出現在存在尖銳外形的計算域,當局部的流速過大或者壓差過大時出錯,把局部的網格加密再試一下。 在fluent 中,用courant number 來調節計算的穩定性與收斂性。一般來說,隨著courantnumber 的從小到大的變化,收斂速度逐漸加快,但是穩定性逐漸降低。所以具體的問題,在計算的過程中,最好是把ourant number 從小開始設置,看看迭代殘差的收斂情況,如果收斂速度較慢而且比較穩定的話,可以適當的增加courant number 的大小,根據自己具體的問題,找出一個比較合適的courant number,讓收斂速度能夠足夠的快,而且能夠保持它的穩定性。個人認為這應

13、該和你采用的算法有關SIMPLE算法是根據連續方程推導出壓力修正方法求解壓力。由于連續方程中流場耦合項被過渡簡化,使得壓力修正方程不能準確反映流場的變化,從而導致該方程收斂緩慢。試著用SIMPLEC算法看看。ANSYS仿真計算代做:模態分析,瞬態動力學,諧響應分析和譜分析、械結構的疲勞、損傷,CFD流體;結構的強度評估和優化; 企鵝:690294845利用FLUENT不收斂通常怎么解決?、一般首先是改變初值,嘗試不同的初始化,事實上好像初始化很關鍵,對于收斂。、FLUENT的收斂最基礎的是網格的質量,計算的時候看怎樣選擇CFL數,這個靠經驗、首先查找網格問題,如果問題復雜比如多相流問題,與模型

14、、邊界、初始條件都有關系。、有時初始條件和邊界條件嚴重影響收斂性,曾經作過一個計算反反復復,通過修改網格,重新定義初始條件,包括具體的選擇的模型, 還有老師經常用的方法就是看看哪個因素不收斂,然后尋找和它有關的條件,改變相應參數。就收斂了、A.檢查是否哪里設定有誤:比方用mm的unit建構的mesh,忘了scale;比方給定的邊界條件不合理。B從算至發散前幾步,看presure分布,看不出來的話,再算幾步, 看看問題大概出在那個區域。 C網格,配合第二點作修正,就重建個更漂亮的,或是更粗略的來處理。D再找不出來的話,換個solver。、解決的辦法是設幾個監測點,比如出流或參數變化較大的地方,若

15、這些地方的參數變化很小,就可以認為是收斂了,盡管此時殘值曲線還沒有降下來。、調節松弛因子也能影響收斂,不過代價是收斂速度。亞松弛因子對收斂的影響所謂亞松馳就是將本層次計算結果與上一層次結果的差值作適當縮減,以避免由于差值過大而引起非線性迭代過程的發散。用通用變量來寫出時,為松馳因子(Relaxation Factors)。數值傳熱學-214FLUENT中的亞松馳:由于FLUENT所解方程組的非線性,我們有必要控制變化。一般用亞松馳方法來實現控制,該方法在每一部迭代中減少了變化量。亞松馳最簡單的形式為:單元內變量等于原來的值加上亞松馳因子a與變化的積: 分離解算器使用亞松馳來控制每一步迭代中的計

16、算變量的更新。這就意味著使用分離解算器解的方程,包括耦合解算器所解的非耦合方程(湍流和其他標量)都會有一個相關的亞松馳因子。在FLUENT中,所有變量的默認亞松馳因子都是對大多數問題的最優值。這個值適合于很多問題,但是對于一些特殊的非線性問題(如:某些湍流或者高Rayleigh數自然對流問題),在計算開始時要慎重減小亞松馳因子。使用默認的亞松馳因子開始計算是很好的習慣。如果經過4到5步的迭代殘差仍然增長,你就需要減小亞松馳因子。有時候,如果發現殘差開始增加,你可以改變亞松馳因子重新計算。在亞松馳因子過大時通常會出現這種情況。最為安全的方法就是在對亞松馳因子做任何修改之前先保存數據文件,并對解的

17、算法做幾步迭代以調節到新的參數。最典型的情況是,亞松馳因子的增加會使殘差有少量的增加,但是隨著解的進行殘差的增加又消失了。如果殘差變化有幾個量級你就需要考慮停止計算并回到最后保存的較好的數據文件。注意:粘性和密度的亞松馳是在每一次迭代之間的。而且,如果直接解焓方程而不是溫度方程(即:對PDF計算),基于焓的溫度的更新是要進行亞松馳的。要查看默認的亞松弛因子的值,你可以在解控制面板點擊默認按鈕。對于大多數流動,不需要修改默認亞松弛因子。但是,如果出現不穩定或者發散你就需要減小默認的亞松弛因子了,其中壓力、動量、k和e的亞松弛因子默認值分別為0.2,0.5,0.5和0.5。對于SIMPLEC格式一

18、般不需要減小壓力的亞松弛因子。在密度和溫度強烈耦合的問題中,如相當高的Rayleigh數的自然或混合對流流動,應該對溫度和/或密度(所用的亞松弛因子小于1.0)進行亞松弛。相反,當溫度和動量方程沒有耦合或者耦合較弱時,流動密度是常數,溫度的亞松弛因子可以設為1.0。 對于其它的標量方程,如漩渦,組分,PDF變量,對于某些問題默認的亞更松弛可能過大,尤其是對于初始計算。你可以將松弛因子設為0.8以使得收斂容易。看了流量是否平衡在report->flux里面操作,mass flow rate,把所有進出口都選上,compute一下,看看nut flux是什么水平,如果它的值小于總進口流量的1

19、%,并且其他檢測量在繼續迭代之后不會發生波動,也可以認為你的解是收斂的。造成連續方程高殘差不收斂的原因主要有以下幾點:1.網格質量,主要可能是相鄰單元的尺寸大小相差較大,它們的尺寸之比最好控制在1.2以內,不能超過1.4.2.離散格式及壓力速度耦合方法,如果是結構網格,建議使用高階格式,如2階迎風格式等,如果是非結構網格,除pressure保持standard格式不變外,其他格式改用高階格式;壓力速度耦合關系,如果使用SIMPLE,SIMPLEC,PISO等segerated solver對聯系方程收斂沒有提高的話,可以嘗試使用coupled solver。另外,對于梯度的計算,不論使用結構或

20、非結構網格,都可以改用node-based來提高計算精度。在FLUENT的中文幫助里,對收斂有比較詳細地描述,建議去看看。收斂性可能會受到很多因素影響。大量的計算單元,過于保守的亞松馳因子和復雜的流動物理性質常常是主要的原因。有時候很難確定你是否得到了收斂階。 沒有判斷收斂性的普遍準則。殘差定義對于一類問題是有用的,但是有時候對其它類型問題會造成誤導。因此,最好的方法就是不僅用殘差來判斷收斂性而且還要監視諸如阻力、熱傳導系數等相關的積分量。 在FLUENT的幫助文件里給出了下面幾種典型的情況。 1、如果你對流場的初始猜測很好,初始的連續性殘差會很小從而導致連續性方程的標度殘差很大。在這種情況下

21、,檢查未標度的殘差并與適當的標度如入口的質量流速相比較是很有用的。什么是標度殘差?就是選作用來標準化的殘差值,一般是取第五步吧,所以,一開是殘差就很小,那么,后面的殘差和它一比,值也很難收斂到很小數。 2、對于某些方程,如湍流量,較差的初始猜測可能會造成較高的標度因子。在這種情況下,標度的殘差最開始會很小,隨后會呈非線性增長,最后減小。因此,最好是從殘差變化的行為來判斷收斂性而不僅僅是殘差的本身值來判斷收斂性。你應該確認在幾步迭代(比如說50步)之后殘差繼續減小或者仍然保持較低值,才能得出收斂的結論。 另一個判斷收斂性的流行方法就是要求未標度的殘差減小到三階量級。為了實現這一方法,FLUENT

22、提供了殘差標準化,有關殘差標準化的信息請參閱分離求解器殘差定義和耦合求解器殘差定義兩節。在這種方法中,要求標準化的未標度殘差降到10-3。但是這種要求在很多情況下可能是不合適的。   1、如果你提供了較好的初始猜測,殘差可能不會降到三階量級。比方說,在等溫流動中,如果溫度的初始猜測非常接近最終值,那么能量殘差根本就不會降到三階量級。   2、如果控制方程中包括的非線性源項在計算開始時是零,但是在計算過程中緩慢增加,殘差是不會降到三階量級的。例如,在封閉區域內部的自然對流問題,由于初始的均一溫度猜測不會產生浮力,所以初始的動量殘差可能非常接近零。在這種情況下,初始的接近零的殘

23、差就不適合作為殘差的較好的標度。    3、如果所感興趣的變量在所有的地方都接近零,殘差不會降到三階量級。例如,在完全發展的管流中,截面上的速度為零。如果這些速度初始化為零,那么初始的和最終的殘差都接近零,因此也就不能期待降三階量級。 在這種情況下,最好監視諸如阻力、總熱傳導系數等積分量來判斷解的收斂。檢查非標準化未標度的殘差來確定這個殘差和適當的標度相比是不是很小也是很有用的。 相反,如果初始猜測很差,初始的殘差過大以至于殘差下降三階量級也不能保證收斂。這種情況對于初始猜測很難的k和e方程尤其常見。在這里,檢查你所感興趣的所有積分量就很有用了。如果解是不收斂的,你可以減少收斂公差FLUENT運行過程中,出現殘差曲線震蕩原因?如何解決殘差震蕩的問題?殘差震蕩對計算收斂性和計算結果有什么影響?一. 殘差波動的主要原因:1、高精度格式; 2、網格太粗;3、網格質量差;4、流場本身邊界復雜,流動復雜;5、模

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網頁內容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
  • 4. 未經權益所有人同意不得將文件中的內容挪作商業或盈利用途。
  • 5. 人人文庫網僅提供信息存儲空間,僅對用戶上傳內容的表現方式做保護處理,對用戶上傳分享的文檔內容本身不做任何修改或編輯,并不能對任何下載內容負責。
  • 6. 下載文件中如有侵權或不適當內容,請與我們聯系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論