沖擊波第一講

沖擊波理論

——研究生課程主講人：彭金華Email：pengjh@教學目的

本課程旨在比較深入、系統(tǒng)地介紹氣體中運動的定常、非定常沖擊波傳播及與其它間斷面的相互作用，使學生掌握基本物理概念和計算方法，以便為開展科學研究和解決有關工程技術問題奠定基礎。課程大綱1基本概念和方程（6學時）1.1守恒方程1.2介質狀態(tài)方程1.3理想流體運動方程組1.4伯努力方程1.5不可壓縮流體運動方程組1.6流體力學方程組的積分形式1.7間斷面及間斷關系式第一講第二講

課程大綱（續(xù)）2正沖擊波（15學時）2.1沖擊波基本概念和關系式2.2多方氣體沖擊波關系式2.3凝聚介質沖擊波關系式2.4雨貢紐曲線及瑞利曲線2.5沖擊波基本性質2.6沖擊波熵增及耗散過程2.7弱沖擊波的聲學近似2.8沖擊波的相互作用2.9沖擊波與稀疏波的相互作用2.10沖擊波與交界面的相互作用2.11初始間斷分解第三講第四講第五講第六講第七講課程大綱（續(xù)）3斜沖擊波（6學時）3.1斜沖擊波極曲線3.2斜沖擊波在固壁上的正規(guī)反射3.3斜沖擊波在固壁上的馬赫反射3.4斜沖擊波在自由面上的正規(guī)反射3.5斜沖擊波在物質界面上的正規(guī)折射3.6兩沖擊波斜碰撞4非定常沖擊波傳播（3學時）4.1二維Whitham方法4.2沖擊波的繞射4.3點爆炸問題的自模擬解4.4球面沖擊波的聚心運動5沖擊波技術的應用第八講第九講第十講教材選用1)李維新.一維不定常流與沖擊波.北京：國防工業(yè)出版社.20032）周毓麟.一維非定常流體力學.北京：科學出版社.19983）王繼海.二維非定常流和激波.北京：科學出版社.1994考核上課出勤率，回答問題及聽課情況，占總成績10%；學期中，每人寫一篇讀書報告或準備一節(jié)課的教學內容，上講臺交流，占總成績20%；學期末，開卷考試，考試時間2小時，試卷分100分，占總成績70%。第一章基本概念和方程

1.1守恒方程質點：介質的微元叫作“流體質點”或“質點”。當說質點速度時，指的并非各分子的速度，而是微元整體的速度，當說到質點密度、壓力等狀態(tài)量時，指的則是該微元體現(xiàn)的平衡態(tài)宏觀量。宏觀小、微觀大

守恒方程的一般形式

強度量：單位體積的量，例如密度、動量密度、能量密度、壓力等，這類量不隨體積的增加而增加；廣延量：強度量對體積積分的結果，例如質量、動量、能量、熵等，這類量對體積是可加的。設L(r，t)是所討論宏觀系統(tǒng)中介質的某一強度量，它是空間坐標r=r(x，y，z)和時間t的函數(shù)。在系統(tǒng)中任取一個體積V，則L(r，t)對應的廣延量是當L是一守恒量時，對于非孤立系統(tǒng)，的變化由兩項組成：一項是單位時間內在體積V內ψ的產(chǎn)生項，即源項，把它記作P（ψ）；另一項是單位時間內通過體積V的表面積S流走ψ的流項，將它記作J（ψ），即（1.1）這里ψ（t）只是t的函數(shù)，故與的含義相同。對P(ψ)和J(ψ)也可用其相應的強度量表出其中σ是單位時間單位體積內ψ的源，而其中j是單位時間內通過表面單位面積的ψ的流，這里j和面積dS都是矢量，定義表面積的外法線方向為正。一般形式守恒方程的積分形式（1.2）再利用格林(Green)公式把式中最后一項的面積分化為體積分，上式可化為（1.3）其中▽是符號算子，在直角坐標系(x，y，z)中因(1.3)式對任意的體積V都成立，當所有的量在V內是連續(xù)變量時，該式就意味著積分號內整個被積函數(shù)應等于零，故得守恒方程的微分形式（1.4）對于孤立系統(tǒng)，不存在與外界的交換，也無源，這時ψ的守恒方程為這里和以后都用表示當?shù)氐臅r間微商，以表示隨體微商，它們的關系是其中u=u（u，v，w）是介質的速度矢量。質量守恒方程

質量對應的強度量，即單位體積的質量是密度ρ，現(xiàn)令L=ρ。因質量不產(chǎn)生也不消亡，故源項σ=0。ρ的流只有運流，故流項j=ρu，這里u是介質的宏觀速度。于是，代人(1.4)式得（1.5）這就是熟知的質量守恒方程，也稱為連續(xù)性方程。展開上式中的散度所以或者若在運動過程中介質的ρ始終保持不變，即dρ／dt=0，則這種介質稱為不可壓縮介質。對不可壓縮介質，連續(xù)性方程特別簡單，為動量守恒方程

動量的強度量是動量密度ρu，即現(xiàn)在令L=ρu。當存在外力場的作用時，根據(jù)牛頓定律，外力對介質的作用將導致介質動量增加，故外力是產(chǎn)生動量的源。設F是作用于介質單位質量的外力，則ρF為作用于單位體積的外力，于是動量密度的源σ=ρF。動量密度本身是一個矢量，它的流則應是個張量。其中運流即隨質點運動帶走的動量密度流是ρuu，這里ρuu是并矢張量，例如分量ρuux就代表動量ρu在x方向的流量。另外是擴散流，因為介質中的應力張量∏要導致動量的擴散，所以在所討論系統(tǒng)的表面積上將產(chǎn)生流過該面積的擴散流－∏，這里取負號是因為應力朝表面積外法向為正，故應力給外界產(chǎn)生的動量為正，而給本系統(tǒng)產(chǎn)生的動量則為負。所以，動量密度的流，j=ρuu－∏。于是，根據(jù)(1.4)式得動量守恒方程（1.6）所以（1.7）納維—斯托克斯(Navier—Stokes)方程

粘性流體的動量方程，其標量形式（1.8）其中μ是粘性系數(shù)，稱為運動粘性系數(shù)歐拉(Euler)方程

對于不可壓縮粘性流體，(1.8)式化簡為（1.9）對于非粘性流體，在無外力作用情況下，動量守恒方程就化為（1.10）這個方程也叫作歐拉(Euler)方程。

能量守恒方程

單位體積的總能為ρE，即令L=ρE（1.11）總能的源有兩部分，一是介質本身釋放的能量，二是外力F對介質做的功，即總能的流包括：①隨介質運動帶走的能量，即運流ρEu；②因熱傳導而在單位時間內流過單位面積的能量流q；③應力單位時間內在單位面積上所做的功。于是能量流項為將以上各項代人(1.4)式，就得到總能守恒方程為(1.12)或寫為(1.13)并利用到質量守恒方程(1.5)，則(1.12)式可化為(1.14)內能守恒方程(1.15)常用的內能守恒方程(1.16)也稱為內能平衡方程。它表明，介質內能的增量等于如下幾項之和：①周圍介質對本介質做的壓縮功，即；②外界向介質輸入的熱量；③介質表面上應力做的功；④介質本身釋放的能量。當無能源、無耗散應力時，內能守恒方程則為(1.17)這表明，外界向介質輸入的熱量，將用于增加介質的內能和使介質對外做功。這就是大家熟知的熱力學第一定律。在無能源、無熱傳導、無耗散作用的腈況下，內能守恒方程非常簡單，即(1.18)守恒方程小結最一般形式的流體動力學方程組：（1.19）非守恒形式的流體動力學方程組：（1.20）在無能源、無外力、無熱傳導的情況下，粘性流體動力學方程組為（1.21）同上情況下，非粘性流體動力學方程組是（1.22）

若把這組方程寫為隨體微商，即拉格朗日(Lagrange)時間微商的形式則為

（1.23）

以上得到的流體動力學方程組，其方程個數(shù)是五個(其中動量方程是三個)，而方程中待求物理量為ρ、p、e、u（u，v，w）共六個，比方程的個數(shù)多一個。一維運動情況也如此，方程是三個，待求量共四個。為了對問題求解．還需再補充一個方程，這就要給出一個表達狀態(tài)量ρ、p和e之間關系的方程，即狀態(tài)方程。所以，求解流體動力學問題，除流體動力學方程組外，還需再加一個狀態(tài)方程，才能構成封閉方程組。1.2介質狀態(tài)方程四個熱力學關系式是

（1.24）

第一式是熱力學第一定律，其余各式是由第一式及如下定義導出的。（1.25）hpτepτFTsGTs狀態(tài)方程是涉及介質具體性質的熱力學量之間的關系式，通常是指介質的p，τ，T之間的關系式，并常用τ，T或p，T為自變量，即或有時也把內能函數(shù)視為狀態(tài)方程。在英文文獻中狀態(tài)方程equationofstate（EOS）一詞通常指p=p(τ,T)，有時也稱它為溫態(tài)方程thermalEOS,e=e(τ,T)稱為能態(tài)方程caloricEOS。根據(jù)熱力學理論，有了以上兩個方程，介質的熱力學性質就全部知道了。在流體動力學中通常多采用p、ρ、e之間關系的狀態(tài)方程，即或一般流體(氣體)的性質

（1）熵不變時，壓力總是隨密度的增加(比容的減小)而增加。（1.26）當ρ=0時，=0。所以永遠為正，于是可以定義一個如下的恒正的量：（1.27）c稱為聲速，是一個很重要的量。(2)熵不變時，聲速將隨密度的增加而增加。即有（1.28）(3)比容不變時，壓力隨熵的增加而增加。即有（1.29）(4)對于氣體，其密度可趨近于零，再作以下假定：當ρ→0時（1.30）

理想氣體

pτ=RT（1.31）式中R是常數(shù)．它等于氣體普適常數(shù)除以氣體的摩爾質量。理想氣體有時也叫作完全氣體。（1.32）（1.33）（1.33）（1.34）（1.33）多方氣體

當理想氣體的比定容熱容cV為常數(shù)時，則由(1.32)式積分可得（1.35）由熱力學第一定律對多方氣體得常用的多方氣體的狀態(tài)方程（1.36）對于等熵過程，A(S)為常數(shù)，故多方氣體的等熵狀態(tài)方程為（1.37）指數(shù)γ是比熱比，也稱為多方指數(shù)或絕熱指數(shù)。從γ的定義(1.34)看到，總有γ>1，根據(jù)統(tǒng)計力學和熱力學得知，其中l(wèi)是氣體粒子的內部自由度。對于單原子氣體(如氫、氖)，內部自由度l=0，故γ=5／3；對雙原子氣體(如氧、氮、空氣等)，在溫度不高時有兩個轉動自由度，l=2，故γ=7／5，當溫度較高時，振動自由度被激發(fā)，自由度增加到l=4，這時γ=9／7。對多方氣體．還容易得到(1.38)(1.39)(1.40)凝聚介質

通?？砂涯劢橘|的壓力和內能分別表示為如下一般形式(1.41)(1.42)式中pc和ec，分別是冷壓和冷能，它們只是比容的函數(shù)，pT和eT是熱壓和熱能，它們同時依賴于比容和溫度。

格留納森