第五講高分辨率格式初步目的為什么那么重視高分辨率格式？是為了提高捕獲激波和間斷面的質量對于光滑流動，無論是中心差分格式，還是矢通量分裂格式，其結果都是令人滿意的，彼此之間無優劣之分。對于含激波這樣強間斷的流動，格式的優劣就顯現出來了。中心差分是頻散格式，求解過程中，激波前后會產生強烈的數值振蕩。加入人工粘性，能夠抑制數值振蕩但是過大的耗散使激波變得很寬（5個網格點左右），從而降低了激波的分辯率。矢通量分裂格式是耗散格式，可以消除激波前后的數值振蕩，但由于精度較低（一階）、耗散較大使得激波的梯度被抹平，激波寬度仍然較大（3個網格點左右）。與中心差分格式不同，矢通量分裂格式屬于迎風格式，與后面要講的高分辯率格式屬同類，它為什么分辯率不高呢？解釋如下：遷移方程ututaux0令(11-1)(11-2)式中,-aa0,a2-a|a|02-aa0,a2-a|a|02(11-3)結合式(11-2),(11-3),式(11-1)成為:迎風格式n1Uiut1

u迎風格式n1Uiut1

ux—ax2ux0

x(11-4)（與矢通量格式相當）為nnuiui11oax2nnui1uix(11-5)經過整理，可將迎風格式寫成如下形式,n1nuiuitnnn1nuiuitnnui1ui1an

ui12uinuin12

x(11-6)從上式可見，當a0時，二階耗散消失了，這時差分格式會得到非物理解。實際上a0的點在實際流動中就是音速點。一維Euler方程擾動傳播速度為ua和ua（這里a表示音速）而音速點恰恰是激波區不可避免的！

（一）激波管（黎曼問題）初始條件為：膈膜左側：uuL,ppL,L，當xx0,t0膈膜右側：uuR,ppR,R,當xx0,t0假定：plpr,lr，而uluR0膈膜破裂之后，t0,將會發生什么？膈膜的初始位置膨脹波接觸間斷激波流動狀態（t0）接觸間斷：p3p2,u3u2,G.OQ7.5O10,B9G.GO245G5.00x(m)7.50LS.QCO-Ofl。?加2.5054。x(m)lOOxLO2O.0Q2,505.勃7.5910,OQ1czb-o-9GD**12101.30M數赫馬o,ocx(m)0<002,505.007.5010.00x(m)x(m)x(m)激波管內流場情況（膈膜破裂之后6.1msec）初始條件為：Pl105N/m2,

L1.0Kg/m3,uL0m/sPR104N/m2,R0.125Kg/m3,uR0m/s.②區：設p/P,1,根據Rankine-Hugoniot條件（簡稱PR1R-H條件）有:u2P1u2P111P工12aRP2PPRfP.③區根據間斷面的定義，可得P3P2PPRfPu3u2fpPR;而p滿足下列關系式:2P11、.2P11、.ip2包1aR1民丁pWPR可用牛頓迭代法從上式中求出P,于是②區和③區的氣流狀態完全確定

設膈膜位置為.膨脹波掃過的區域⑤設膈膜位置為x0,膈膜破裂時間點為t0。則根據Riemann不變量的規律。沿特征線Cdt

dx1、力),有沿特征線Cdt

dx1、力),有ua沿特征線Cdtdx1、力),有ua沿特征線C0dt

dx可得:U5常數3可得:U5常數3l1一Ul2P5PP5PluaLP5Pl可見，當膈膜破裂之后，激波管內氣流參數完全由膈膜兩側的狀態參數所決定，當然考察的位置x和時間t也是決定因素。（二）Riemann問題的理論分析1.對于一維問題，非線性守恒方程系統：ptFx0FFU（11-7a）Ux,0U0（x）x（11-7b）古典解：如果方程（11-7a）的解U連續且U的一階偏導數分段連續，則稱U為古典解。弱解：對于方程（11-7a）,如果它的解U為間斷的不連續函數，則稱其為弱解（也可稱為廣義解）。弱解滿足下列條件：1）解U在所有連續點滿足初始條件（11-7b）;2）在解U的光滑處，方程（11-7a）局部成立；3）在解U的間斷處滿足Rankine-Hugoniot條件F（Ur）FUlC（UrUl）（11-8）式中，C為間斷面的移動速度；下標“R”表示間斷面右側變量，“L”表水左側變量。也可寫成：FUCU式中，f表示函數f在間斷左右的跳躍。

2.Riemann(黎曼)問題:對于Euler方程：(11-9)UFx0(11-9)式中uUuFu2p

eepu如果初始條件為:Ux,0UUx,0ULx0URx0當t0之后，由初始間斷的左右狀態確定波態、波的強度和波之間的流動特性的問題稱為Riemann問題。3.間斷面的分類：1）弱間斷：參數連續，而參數的導數不連續，例如膨脹波（稀疏波）；2）強間斷：參數不連續（產生間斷），例如激波和接觸間斷。激波：滿足R-H條件的間斷；接觸間斷：流體不穿過波面，即滿足：p0,u0即，間斷前后壓力和速度不變，只有密度發生跳躍0可以證明，Riemann問題的解是由一個向左移動的波，一個處于中間的接觸間斷和一個向右移動的波構成，向左向右移動的波或為

激波、或為膨脹波（稀疏波）。三個波之間的參數呈定常狀態。顯然，t0之后，Riemann問題的解只由初始間斷的左右狀態UR、UL所確定，所以Riemann解可寫成：RXUx,tU-,Ul,Ur

t上標R表示Riemann解。Godunov格式Godunov格式的基本原理：在離散點的界面上求解Riemann問題。右圖表示以單元1i—,i—為單位TOC\o"1-5"\h\z2分片平均的幾何意義。由圖可見，單元1,,一i-,i—的平均2值就是Ui。Godunov認為從n1t的未知解時間層tnt的已知解Un1t的未知解時間層tnt的已知解U：求得下一時間層t11...一一第1步：將已知解U：在單兀i—,i-內進行平均i22并得到相鄰單元的平均值,Uin1解:L\ndxuinxr并得到相鄰單元的平均值,Uin1解:L\ndxuinxrUnnni1，Ui1Ui1，第2步：根據激波管原理,在相鄰單元界面上，求得(11-10)RiemannURUR0,Uin1,UinuRuR00nlUn(11-11)i弓界面步：下一時間層n1t的未知解由Riemann解(11-11),、，1在單兀i一,i2內積分獲得:_1步：下一時間層n1t的未知解由Riemann解(11-11),、，1在單兀i一,i2內積分獲得:_1TUin1—URxTx-，Ul,Urdx(11-12)由于界面iRiemann解是不同的，需要分段積分，于是（11-12)式可寫為Uin1x-Ux0R-,Uin1,Uin0UR__x2t,Uin,Uin1d(11-13)在以上三步中,第在以上三步中,第1和第3步均為在單元內的積分，與方程的物理本質無關，而第2步則利用了激波管的物理特性，Godunov求解控制方程的獨特之處就是這第2步：將離散的數值求解化為求