FLUENT中求解方法的區(qū)別與選擇

2016-10-07  by:CAE仿真在線  來源:互聯網

1.非耦合求解 ( Segregated );2.耦合隱式求解 ( Coupled Implicit );3.耦合顯式求解 ( Coupled Explicit )

非耦合求解方法主要用于不可壓縮或壓縮性不強的流體流動。耦合求解則可以用在高速可壓縮流動。FLUENT默認設置是非耦合求解,但對于高速可壓流動,有強的體積力(浮力或離心力)的流動,求解問題時網格要比較密,建議采用耦合隱式求解方法,可以耦合求解能量和動量方程,能比較快地得到收斂解。缺點是需要的內存比較大(是非耦合求解迭代時間的1.5-2倍)。如果必須要耦合求解,但是你的機器內存不夠,這時候可以考慮用耦合顯式解法器求解問題。該解法器也耦合了動量,能量及組分方程,但內存卻比隱式求解方法小。缺點是收斂時間比較長。

這里需要指出的是非耦合求解的一些模型在耦合求解解法器里并不都有。耦合解法器沒有的模型包括:多相流模型,混合分數/PDF燃燒模型,預混燃燒模型,污染物生成模型,相變模型,Rosseland輻射模型,確定質量流率的周期性流動模型及周期性換熱模型等。

隱式( Implicit ):對于給定變量,單元內的未知值用鄰近單元的已知和未知值計算得出。因此,每一個未知值會在不止一個方程中出現,這些方程必須同時解來給出未知量。

顯式( Explicit ):對于給定變量,每一個單元內的未知量用只包含已知量的關系式計算得到。因此未知量只在一個方程中出現,而且每一個單元內的未知量的方程只需解一次就可以給出未知量的值。

一階迎風格式( First Order Upwind ):當需要一階精度時,我們假定描述單元內變量平均值的單元中心變量就是整個單元內各個變量的值,而且單元表面的量等于單元內的量。因此,當選擇一階迎風格式時,表面值被設定等于迎風單元的單元中心值。

二階迎風格式( Second Order Upwind ):當需要二階精度時,使用多維線性重建方法來計算單元表面處的值。在這種方法中,通過單元中心解在單元中心處的泰勒展開來實現單元表面的二階精度值。因此,當使用二階迎風格式時,用下面的方程來計算表面值;

QUICK格式:對于四邊形和六面體網格,我們可以確定它們唯一的上游和下游表面以及單元。FLUENT還提供了計算對流變量在表面處高階值的QUICK格式。QUICK類型的格式是通過變量的二階迎風與中心插值加上適當的權因子得到的;

亞松馳( Under-Relaxation ):由于FLUENT所解方程組的非線性,我們有必要控制的變化。一般用亞松馳方法來實現控制,該方法在每一部迭代中減少了f的變化量。亞松馳最簡單的形式為:單元內變量f等于原來的值f_old加上亞松馳因子a與f變化的積.

SIMPLE:SIMPLE算法使用壓力和速度之間的相互校正關系來強制質量守恒并獲取壓力場。

一階與二階的比較

當流動和網格成一條線時(如:矩形網格或者六面體網格模擬矩形導管的層流流動),可以使用一階迎風離散格式。但是,當流動和網格不在一條線上時(即:流動斜穿網格線)一階對流離散增加了對流離散的誤差(數值耗散)。對于三角形和四面體網格,流動從來就不會和網格成一條線,此時一般要使用二階離散來獲取更高精度的結果。對于四邊形或者六面體網格,如果使用二階離散格式,尤其是對于復雜流動來說,你可以獲取更好的結果。

總而言之,一階離散一般會比二階離散收斂得好,但是精度要差,尤其是對于三角形或者四面體網格精度更差。

對于大多數情況,你可以在計算的開始使用二階格式。對于有些情況,你應該以一階離散開始計算,在進行了初步迭代之后再轉到二階格式。例如,如果你解高馬赫數流動問題,初始解科所預期的解相差較大,你就應該先用一階格式迭代幾步然后打開二階格式繼續(xù)計算直至收斂。

對于與網格成一條線的簡單流動(如:劃分為矩形網格或者六面體網格的矩形導管的層流流動),數值耗散自然會很低,所以一般使用一階格式替代二階格式而不損失精度。

最后,如果你使用二階格式遇到收斂性問題,你就應該嘗試使用一階格式。

選擇壓力插值格式

如壓力插值格式所述,當使用分離求解器時我們可以采用很多壓力插值格式。對于大多數情況,標準格式已經足夠了,但是對于特定的某些模型使用其它格式可能會更好:

l        對于具有較大體積力的問題,推薦使用體積力加權格式。

l        對于具有高渦流數,高Rayleigh數自然對流,高速旋轉流動,包含多孔介質的流動和高度扭曲區(qū)域的流動,使用PRESTO!格式。

注意:PRESTO!只能用于四邊形或者六面體網格。

l        對于可壓流動推薦使用二階格式。

當其它格式不適用時,使用二階格式來提高精度(如:對于流過具有非六面體或者非四邊形網格的曲面邊界的流動。)

選擇壓力速度耦合方法

在分離求解器中,FLUENT提供了壓力速度耦合的三種方法:SIMPLE,SIMPLEC以及PISO。定常狀態(tài)計算一般使用SIMPLE或者SIMPLEC方法,對于過渡計算推薦使用PISO方法。PISO方法還可以用于高度傾斜網格的定常狀態(tài)計算和過渡計算。需要注意的是壓力速度耦合只用于分離求解器,對于耦合求解器你不可以使用它。

SIMPLE與SIMPLEC比較

在FLUENT中,可以使用標準SIMPLE算法和SIMPLEC(SIMPLE-Consistent)算法,默認是SIMPLE算法,但是對于許多問題如果使用SIMPLEC可能會得到更好的結果,尤其是可以應用增加的亞松馳迭代時,具體介紹如下:

對于相對簡單的問題(如:沒有附加模型激活的層流流動),其收斂性已經被壓力速度耦合所限制,你通?？梢杂?/b>SIMPLEC算法很快得到收斂解。在SIMPLEC中,壓力校正亞松馳因子通常設為1.0,它有助于收斂。但是,在有些問題中,將壓力校正松弛因子增加到1.0可能會導致不穩(wěn)定。對于這種情況,你需要使用更為保守的亞松馳或者使用SIMPLE算法。對于包含湍流和/或附加物理模型的復雜流動,只要用壓力速度耦合做限制,SIMPLEC會提高收斂性。它通常是一種限制收斂性的附加模擬參數,在這種情況下,SIMPLE和SIMPLEC會給出相似的收斂速度。