計(jì)算流體力學(xué)課件：有限體積法1

有限體積法（1）知識(shí)點(diǎn)：

1有限體積法的基本概念

無(wú)粘通量及粘性通量的計(jì)算

多塊網(wǎng)格2知識(shí)回顧：Roe格式平均斜率線性化，以平均增長(zhǎng)率代替瞬時(shí)增長(zhǎng)率[j,j+1]區(qū)間內(nèi)連續(xù)，且可通過(guò)相似變換對(duì)角化

應(yīng)當(dāng)具有的性質(zhì)常系數(shù)方程的Riemann解3知識(shí)回顧2：LU-SGSStep1:求解Step2:求解

§7.1結(jié)構(gòu)網(wǎng)格有限體積法有限體積法主要優(yōu)勢(shì)：處理復(fù)雜網(wǎng)格差分法處理復(fù)雜外形——坐標(biāo)變換無(wú)粘項(xiàng)需要計(jì)算9次導(dǎo)數(shù)，計(jì)算量大；守恒性被破壞（1）兩端同乘無(wú)粘項(xiàng)僅需計(jì)算三次導(dǎo)數(shù)5差分法有限體積法優(yōu)點(diǎn)簡(jiǎn)單、計(jì)算量小、易于提高精度本身包含幾何信息，易處理復(fù)雜網(wǎng)格不足差分離散與幾何解耦，難以處理復(fù)雜網(wǎng)格復(fù)雜、不易提高精度均勻流場(chǎng)數(shù)值解（誤差），未必為0幾何誘導(dǎo)誤差：對(duì)于均勻場(chǎng)實(shí)際問(wèn)題：外形復(fù)雜，光滑的結(jié)構(gòu)網(wǎng)格生成困難不易保障幾何守恒性坐標(biāo)變換函數(shù)必須足夠光滑——否則損失精度61.基本概念1）控制體節(jié)點(diǎn)（中心）型控制體與網(wǎng)格型控制體72）積分型控制方程在控制體上積分物理含義：

控制體內(nèi)總質(zhì)量/動(dòng)量/能量的增加=穿過(guò)控制體邊界流入的凈質(zhì)量/動(dòng)量/能量無(wú)粘通量粘性通量8物理含義：?jiǎn)挝粫r(shí)間內(nèi)，無(wú)粘流動(dòng)流過(guò)垂直于n方向的單位面積的質(zhì)量、動(dòng)量和能量法向速度93）有限體積法中物理量的含義4）殘差殘差=凈通量=右端項(xiàng)含義：控制體內(nèi)的平均量（平均質(zhì)量密度、平均動(dòng)量密度、平均能量密度）控制體幾何中心處的當(dāng)?shù)孛芏?、?dòng)量密度、能量密度二階精度近似102.無(wú)粘通量的計(jì)算常用方法（流過(guò)AB邊的通量）：a.利用周圍點(diǎn)的值，計(jì)算出(I+1/2,J)點(diǎn)處的物理量；b.利用該處的物理量，計(jì)算出流過(guò)AB邊的流通量方法1：中心型有限體積法人工粘性項(xiàng)用中點(diǎn)的值代替AB面上的均值，二階精度11方法2：迎風(fēng)型有限體積法Step1:

利用（偏）左側(cè)點(diǎn)及（偏）右側(cè)點(diǎn)的值，計(jì)算出I+1/2,J點(diǎn)的值計(jì)算方法：與差分法完全相同各種差分格式，均可直接使用也稱為“差分格式”該過(guò)程稱為“重構(gòu)”（很多文獻(xiàn)中稱為“插值”）有限體積與有限差分共通之處,可直接使用差分格式12常見(jiàn)的差分格式：2階NND格式minmod(a,b):a,b符號(hào)相反時(shí)取0，符號(hào)相同時(shí)取絕對(duì)值小的3階迎風(fēng)3階MUSCL格式TVD,WENO,GVC,保單調(diào)格式……13重構(gòu)方式：原始變量、守恒變量及特征變量以NND格式為例：守恒變量重構(gòu)原始變量重構(gòu)特征變量重構(gòu)先算出UI+1/2(可用UI和UI+1的算術(shù)平均或Roe平均），再利用該值算出SI+1/2

是控制體內(nèi)的平均值

（稱為數(shù)值流通量）

的含義14重要概念澄清：重構(gòu)與插值A(chǔ).有限體積法：j+1/2j-1/2確實(shí)為f在xj+1/2點(diǎn)的值！

通常做法：1）用計(jì)算出2）u在xj+1/2點(diǎn)的值！關(guān)鍵：是用計(jì)算（稱為重構(gòu)），而不是用計(jì)算（是標(biāo)準(zhǔn)的插值）；否則最高也只能達(dá)到2階精度。15重要概念澄清：重構(gòu)與插值B.有限差分法：j+1/2

h(x)為f(x)的“重構(gòu)函數(shù)”差分法與有限體積法的數(shù)值格式完全相同16Step2:

利用，計(jì)算出通量

方法1：流通矢量分裂（FVS）方法2：通量差分分裂(FDS)通過(guò)Riemann解，獲得通量可利用近似Riemann階（Roe,HLL,HLLC)Roe格式：xyx’y’方法3：AUSM方法壓力項(xiàng)單獨(dú)處理，其余項(xiàng)類似VanLeer分裂小技巧：利用局部坐標(biāo)系，計(jì)算通量時(shí)簡(jiǎn)化為1維問(wèn)題AByx’y’“擴(kuò)展的”一維問(wèn)題x,y坐標(biāo)系下的Riemann問(wèn)題（二維）x’,y’坐標(biāo)系下的Riemann問(wèn)題（一維）x’,y’坐標(biāo)系下，切向速度v’表現(xiàn)相當(dāng)于被動(dòng)標(biāo)量183.粘性通量的計(jì)算關(guān)鍵問(wèn)題：計(jì)算

對(duì)的導(dǎo)數(shù)同樣計(jì)算19方法1.利用Jacobian變換方法2.利用Green積分公式計(jì)算最終公式與方法1相同20

§7.2邊界條件及網(wǎng)格分塊1.邊界條件的處理方法常用方法：虛網(wǎng)格設(shè)置1或2層虛網(wǎng)格（如果采用高精度格式，需設(shè)置更多重虛網(wǎng)格）虛網(wǎng)格點(diǎn)上的值根據(jù)邊界條件給定

絕熱固壁邊界條件超聲速入口邊界條件超聲速出口邊界條件21亞聲速入口邊界條件外部計(jì)算域內(nèi)部邊界特征分析：內(nèi)點(diǎn)提供一個(gè)邊界條件，外部提供其余邊界條件調(diào)節(jié)作用，內(nèi)點(diǎn)速度降低時(shí)，增加入口壓力；

內(nèi)點(diǎn)速度增加時(shí)，降低入口壓力；22亞聲速出口邊界條件外部計(jì)算域內(nèi)部邊界特征分析：外部提供一個(gè)邊界條件；

通常情況下，指定背壓如果邊界處壓力比內(nèi)點(diǎn)高，則邊界處密度增加，速度降低；

如果邊界壓力比內(nèi)點(diǎn)低，則邊界處密度降低，速度升高；232.多塊網(wǎng)格復(fù)雜外形情況下，單塊結(jié)構(gòu)網(wǎng)格很難實(shí)現(xiàn)，需要由多塊網(wǎng)格實(shí)現(xiàn)241）對(duì)接網(wǎng)格對(duì)接網(wǎng)格示意圖相鄰的網(wǎng)格塊共享交界網(wǎng)格點(diǎn)處理方法：

利用虛網(wǎng)格傳遞信息；

利用連接信息，將內(nèi)點(diǎn)的物理量復(fù)制給對(duì)應(yīng)的虛網(wǎng)格點(diǎn)；Block2的值Block1虛網(wǎng)格的值252）重疊網(wǎng)格重疊網(wǎng)格示意圖處理方法：

使用虛網(wǎng)格交換信息；

利用插值，用網(wǎng)格內(nèi)點(diǎn)的信息計(jì)算出另一套網(wǎng)格虛網(wǎng)格上的信息。26

§7.2非結(jié)構(gòu)網(wǎng)格有限體積法簡(jiǎn)介1.結(jié)構(gòu)網(wǎng)格與非結(jié)構(gòu)網(wǎng)格結(jié)構(gòu)網(wǎng)格：網(wǎng)格以二維（或三維）表格的拓?fù)湫问脚帕蟹墙Y(jié)構(gòu)網(wǎng)格：網(wǎng)格點(diǎn)無(wú)法以二維（或三維）表格的拓?fù)渑帕卸S非結(jié)構(gòu)網(wǎng)格單元：三角形、四邊形單元結(jié)構(gòu)網(wǎng)格是特殊的非結(jié)構(gòu)網(wǎng)格27非結(jié)構(gòu)網(wǎng)格示意圖（多塊）結(jié)構(gòu)網(wǎng)格示意圖結(jié)構(gòu)網(wǎng)格非結(jié)構(gòu)網(wǎng)格優(yōu)點(diǎn)計(jì)算效率高；易于提高精度；易于處理邊界層

易于網(wǎng)格生成；尤其適用于復(fù)雜外形不足復(fù)雜外形網(wǎng)格生成困難；處理復(fù)雜外形網(wǎng)格質(zhì)量不易保證，精度降低

計(jì)算效率較低；不易提高精度；不易處理邊界層（可采用混合網(wǎng)格）網(wǎng)格長(zhǎng)寬比：1000邊界層內(nèi)的網(wǎng)格282.二階精度非結(jié)構(gòu)網(wǎng)格有限體積法Step1:利用插值，計(jì)算出P點(diǎn)的物理量值

ABC梯度可以用周圍點(diǎn)的值計(jì)算插值過(guò)程中，可使用限制器（例如梯度型限制器等）例如，利用四個(gè)點(diǎn)處的值，計(jì)算P為AB面中心點(diǎn)29ABCStep2:

利用通量技術(shù)（FVS或FDS），

計(jì)算出P點(diǎn)的流通量

Step3:將P點(diǎn)的流通量乘以面AB的面積，得到穿過(guò)面AB的流通量（以中點(diǎn)的流通量代替整個(gè)面的平均流通量，二階精度）Step4:

時(shí)間推進(jìn)303.高階精度有限體積法（k-exact）需要高精度計(jì)算面AB的通量需要計(jì)算該面上多個(gè)點(diǎn)的值（通常為高斯點(diǎn)）Step1:通過(guò)高階插值，計(jì)算出AB面上多個(gè)點(diǎn)的值：二維、三維的高階插值比較復(fù)雜ABC假設(shè)物理量為多項(xiàng)式分布k=1,2,3,4,5定出系數(shù)高精度插值，需要非常多的單元以確定系數(shù)，難度很大31Step2:

利用通量技術(shù)（FVS或FDS），

計(jì)算出Pa,Pb,Pc點(diǎn)的流通量

ABCStep3:

利用數(shù)值積分（例如高斯積分），計(jì)算面AB的通量Step4:

時(shí)間推進(jìn)32（選作題）作業(yè)7.1編制有限體積程序,并計(jì)算圖示鈍楔繞流問(wèn)題r=1幾何及流動(dòng)參數(shù)