有限體積法介紹

1、有限體積法1 有限體積法基本原理上一章講到的有限差分法將數(shù)值網(wǎng)格的節(jié)點上定義為計算節(jié)點，并在網(wǎng)格節(jié)點上對微分形式的流體基本方程進行離散，用網(wǎng)格節(jié)點上的物理量的代數(shù)方程作為原PDE的近似。在本章所要學習的有限體積法則采用了不同的離散形式。首先，有限體積法離散的是積分形式的流體力學基本方程：(1)計算域用數(shù)值網(wǎng)格劃分成若干小控制體。和有限差分法不同的是，有限體積法的網(wǎng)格定義了控制體的邊界，而不是計算節(jié)點。有限體積法的計算節(jié)點定義在小控制體內部。一般有限體積法的計算節(jié)點有兩種定義方法，一種是將網(wǎng)格節(jié)點定義在控制體的中心，另一種方法中，相鄰兩個控制體的計算節(jié)點到公共邊界的距離相等。第一種方法的優(yōu)點在于

2、用計算節(jié)點的值作為控制體上物理量的平均值具有二階的精度；第二種方法的好處是在控制體邊界上的中心差分格式具有較高的精度。積分形式的守恒方程在小控制體和計算域上都是成立的。為了獲得每一個控制體上的代數(shù)方程，面積分和體積分需要用求面積公式來近似。2 面積分的近似采用結構化網(wǎng)格，在二維情況下，每一個控制體有4個面，二維情況，每一個控制體有6個表面。計算節(jié)點用大寫字母表示，控制體邊界和節(jié)點用小寫字母表示。為了保證守恒性，控制體不能重疊，每一個面都是相鄰兩個控制體的唯一公共邊界。控制體邊界上的積分等于控制體個表面的積分的和：(2)上式中，f可以表示或。PENWSNNEEWWSSNWNESWSEewnsne

3、nwsesw顯然，為了獲得邊界上的積分，必須知道f在邊界上的詳細分布情況，這是不可能實現(xiàn)的，由于只是計算節(jié)點上的函數(shù)值，因此必須采用近似的方法來計算積分。整個近似過程分成兩步第一步：用邊界上幾個點的近似積分公式第二步：邊界點上的函數(shù)值用計算節(jié)點函數(shù)值的插值函數(shù)近似面積分可采用以下不同精度的積分公式：二階精度積分：(3)上式中為邊界中點出的函數(shù)值。近似為方格中心點的值乘以方格的面積。三階精度積分：(4)四階精度積分：(5)應該注意的是，采用不同精度的積分公式，在相應的邊界點的插值時也應采用相應精度的插值函數(shù)。積分公式的精度越高，近似公式就越復雜。3 體積分的近似和面積分相似，體積分也有不同精度的

4、近似公式二階精度積分公式(6)采用雙二次樣條函數(shù)(7)可以得到四階精度的積分公式：(8)4 函數(shù)的插值在上節(jié)講到的積分的近似公式中用到了非計算節(jié)點上的函數(shù)值，被積函數(shù)f中包含了多個物理量及其偏微分，如對流項,擴散項，在源項中也有類似情況，這里假定流場和流體的物性參數(shù)是已知的，物理量及其偏導數(shù)在控制面上的值需要通過計算節(jié)點上物理量的插值得到。下面已e面為例進行討論。4.1 迎風插值（UDS）用上游計算節(jié)點的函數(shù)值近似相當于對一階偏導數(shù)采用迎風格式，因此用UDS來表示這種近似方法，在UDS中：(9)UDS是唯一無條件滿足有界性要求的近似格式，在數(shù)值過程中不會產生數(shù)值振蕩。UDS存在數(shù)值粘性。根據(jù)T

5、aylor公式，該格式具有一階精度，并具有數(shù)值粘性：(10)在多維問題中，如果流動方向和網(wǎng)格是斜交的，截斷誤差會在垂直于流動方向以及流線方向產生擴散，這是一種非常嚴重的誤差，函數(shù)的峰值或函數(shù)值的快速變化會被抹平，為了得到高精度結果需要采用非常精細的網(wǎng)格。4.2 線性插值（CDS）(11)(12)線性插值具有二階精度，線性插值相當于FDM中的CDS格式，因此用CDS表示。CDS格式會產生數(shù)值振蕩。對于擴散項(13)4.3 三階迎風格式（QUICK）和UDS類似，QUICK格式也和流動方向有關(14)其中：；(15a)；(15b)4.4 高階格式（4階精度CDS）采用三次曲線可擬合出四階精度的中心

6、插值公式，在均勻網(wǎng)格中，四階公式為：(16)(17)5 邊界的處理對于對流項，在入口處一般給出了流量或函數(shù)值，在邊界和對稱面上流量為零，在出口處假設和出口的法向坐標無關，因此可采用迎風格式。對于擴散項則可能需要采用偏心格式。6 有限體積法應用舉例例：考慮一標量在已知流場中的輸運過程（如圖4.4所示），輸運方程為：(18)邊界條件：；北部入口邊界；西部壁面邊界對稱條件；南部邊界梯度為0；東部出口條件，流線方程f(y)壁面對稱邊界入口，f0出口流線，xy=c對流項：(19)為質量通量。(20)若采用UDS格式，代數(shù)方程組中各項系數(shù)為：；(21)若采用CDS格式，代數(shù)方程組中各項系數(shù)為；(22)根據(jù)

7、連續(xù)性方程：(23)相鄰CV之間的關系：；(24)其余相鄰CV有類似關系擴散項采用CDS格式(25)代數(shù)方程組中擴散項系數(shù)為：；(26)對于任意控制體(27)，l為任意指標P，E，W，S，N。(28)邊界條件的處理：對于西部和北部邊界，由于給定了函數(shù)值，對流項可直接代入函數(shù)值而無需插值，擴散項則采用一側差分(29)這里，W點和P的w邊中點重合。南邊和西邊的梯度為零，以南邊為例，由于梯度為零，代數(shù)方程變?yōu)椋?30)6 SIMPLE方法考慮定常不可壓流動問題，控制方程為：連續(xù)性方程：(31)動量方程：(32)不可壓縮問題求解的困難在于壓力場的求解。主要原因在于壓力p沒有獨立的方程組。先考慮一維問題

8、：對于動量方程：(33)若采用CDS格式簡化后得：(34)根據(jù)連續(xù)性方程，則有，由于相鄰節(jié)點之間的壓力沒有聯(lián)系方程，容易造成壓力交錯現(xiàn)象。為了解決這一問題，可采用交錯網(wǎng)格技術，即速度場和壓力場采用不同的網(wǎng)格。以二維問題為例，交錯網(wǎng)格的布置如下圖所示：PEWNSueuwvnvs主控制體為壓力控制體（黑色實線網(wǎng)格），u的控制體（紅色虛線網(wǎng)格）的計算節(jié)點在主控制體的e邊，控制體的e，w邊界通過主控制體的計算節(jié)點，v控制體（藍色雙點劃線網(wǎng)格）的計算節(jié)點在主控制體的n邊，該控制體的n，s面經過主控制體的計算節(jié)點。在u的控制體中，采用有限體積法離散可得u的代數(shù)方程：(35)壓力場的求解采用壓力校正方法。即采用預估的壓力場求速度，再用連續(xù)性方程校正壓力場。當連續(xù)性方程得到滿足時，壓力場就是真實的壓力場。具體步驟如下1 預測壓力場2 將預測壓力場代入動量方程，分別求解速度場(36)3 用連續(xù)性方程校正壓力設方程的精確解為u，v，p；(37)其中