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

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

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

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

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

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

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

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