什么是CFD中的"矩陣預(yù)條件"!
2017-03-30 by:CAE仿真在線 來源:互聯(lián)網(wǎng)
任何一個矩陣方程組都可以通過直接求解(高斯消去法或LU分解法)來計算。但是,LU分解之后會把原本稀疏的矩陣變成密集矩陣。所以直接求解法會耗費非常多的計算資源。同時,CFD計算中離散的誤差遠大于矩陣求解的誤差,因此我們沒必要那么精確的計算矩陣方程。只要比離散精度高就可以。
對于稀疏矩陣,通常我們選用迭代求解技術(shù),通常我們估計一個初始值,然后系統(tǒng)的使這個值逼近真實解。如果每次迭代很經(jīng)濟并且需要很少的迭代次數(shù)的話,那么迭代求解就比直接求解經(jīng)濟的多了。在CFD中,我們通常使用迭代求解技術(shù)。
然而,迭代求解器一個眾所周知的弱點就是不是很穩(wěn)健。雖然迭代求解器可以快速的求解矩陣系統(tǒng),但是這個缺點導致其在工業(yè)應(yīng)用上收到了嚴重的限制。穩(wěn)健性和效率可以通過預(yù)條件(pre-conditioning)來提高。預(yù)條件矩陣M可以有很多不同的選擇,但是通常其滿足一些最低的要求。比如,求解Mx=b需要是非常經(jīng)濟的,同時,M應(yīng)該和A接近并且應(yīng)該是非奇異的。
舉例:我們有矩陣系統(tǒng):
其中phi是我們需要求的。如果我們另:A=M-N,我們有:
移項有(注:不要被下面的n和n-1搞暈了,在迭代收斂的時候,n和n-1的phi是相等的,因此下面這個方程就是上面那個方程):
方程左右乘以M的逆,最后有:
在CFD中,我們稱M即為預(yù)條件(預(yù)處理)矩陣(Pre-conditioning Matrix)。
一種非常受歡迎的預(yù)條件矩陣建立于不完全LU分解上,有關(guān)不完全LU分解,coming soon!
Computational Methods for Fluid dynamics. J.H. Ferziger. Page.67.
相關(guān)標簽搜索:什么是CFD中的"矩陣預(yù)條件"! Fluent培訓 Fluent流體培訓 Fluent軟件培訓 fluent技術(shù)教程 fluent在線視頻教程 fluent資料下載 fluent分析理論 fluent化學反應(yīng) fluent軟件下載 UDF編程代做 Fluent、CFX流體分析 HFSS電磁分析