特征線理論及應(yīng)用課件.ppt

第二章特征線理論及應(yīng)用 氣體動(dòng)力學(xué)中 有大量問題是用雙曲型偏微分方程來描述的 很難得到解析結(jié)果 在這種情況下 有兩種數(shù)值解法 1 特征線數(shù)值解法 求解域用特征線網(wǎng)格進(jìn)行離散 求各網(wǎng)格結(jié)點(diǎn)上的解 氣體動(dòng)力學(xué)中 有大量流動(dòng)問題是用雙曲型偏微分方程來描述的 宜于用特征線方法求解 2 有限差分法 求解域的有限差分網(wǎng)格一般是正交的 根據(jù)由偏微分方程構(gòu)造的差分格式來求各網(wǎng)格結(jié)點(diǎn)上的解 2 1特征線理論 特征線的數(shù)學(xué)定義 考慮一個(gè)一般的一階雙曲型偏微分方程 x y是兩個(gè)自變量 u x y 是因變量 系數(shù)A1 A2及非齊次項(xiàng)F1可以是x y u的函數(shù) 1 將偏微分方程改寫為 設(shè)未知函數(shù)u x y 連續(xù) u的一階導(dǎo)數(shù)可以寫作 注 u的一階導(dǎo)數(shù)可以不連續(xù) 偏微分方程的特征線定義為 xy平面內(nèi)具有斜率為的曲線 2 3 沿著特征線 或 偏微分方程可化簡(jiǎn)為 代入式 4 得到偏微分方程的相容方程 是平面上這樣一族曲線 沿著此族中任一曲線 a 可以把待求物理量的一階偏微分控制方程變換成等價(jià)的常微分控制方程 b 稱為原偏微分方程或偏微分方程組的相容方程 特征線的第一個(gè)數(shù)學(xué)意義 a b 特征線的第二個(gè)數(shù)學(xué)意義 上兩式表明 沿著特征線 分母和分子均為零 即沿著特征線 表明 1 沿特征線因變量的一階導(dǎo)數(shù)具有不定值 可以是不連續(xù)的 在這種情況下 特征線是弱間斷 第一類間斷線 2 在氣體動(dòng)力學(xué)中 特征線可以是弱擾動(dòng)波傳播的跡線 或者說弱擾動(dòng)傳播的跡線就是特征線 因此 因變量的一階導(dǎo)數(shù)只允許有弱間斷 如果在物理平面上有激波出現(xiàn) 在強(qiáng)間斷面上便無法建立因變量的全微分式 也就不能用特征線方法求解 例 一階偏微分方程 的初始條件是 2 沿此特征線的相容方程 3 u 2 4 的值 用特征線法確定 1 通過點(diǎn) 2 4 的特征線 解 1 對(duì)照一般形式的雙曲型偏微分方程 該方程對(duì)應(yīng)的系數(shù) A1 1 A2 2x F1 3x2 則特征線方程為 積分得 為確定過點(diǎn) 2 4 的特征線 將x 2 y 4 代入上式得 所以 所求的特征線方程是 對(duì)上式積分 得 2 偏微分方程的相容方程為 如何確定C2 初始條件u 0 y 5y 10 及特征線方程 u 0 0 10 因此相容方程為 2 2一維等熵流動(dòng)的特征線數(shù)值解法 基本方程與黎曼不變量 連續(xù)方程 動(dòng)量方程 以一維等直截面管為例 基本方程 等熵流動(dòng)中只有一個(gè)狀態(tài)參量獨(dú)立 將基本方程中的用代替 得 基本方程可化為 定義 則 基本方程化為以v G為新的未知函數(shù)的偏微分方程 基本方程 偏微分方程 特征線 相容方程 在x t平面上 把dx dt v c曲線稱為偏微分方程的特征線 C C x t C 表示第一族特征線 C 表示第二族特征線 解相容方程 對(duì)多方氣體 其相容方程的解為 由 聲速 沿著特征線 沿著特征線 結(jié)論 特征線的基本性質(zhì) 1 一維非定常流動(dòng)中 平面x t上任一點(diǎn) 都有兩條不同族的特征線 沿各特征線有各自不同的黎曼不變量 2 特征線上參量v c p 的一階導(dǎo)數(shù)可以不連續(xù) 但這些參量本身是連續(xù)的 稱因變量的一階導(dǎo)數(shù)不連續(xù)的點(diǎn)叫做弱間斷 如果初始某一點(diǎn)有弱間斷 那么這個(gè)弱間斷必定會(huì)沿著過該點(diǎn)的特征線向外傳播 3 兩個(gè)相鄰的 不同類型流動(dòng)區(qū)域的分界線 必定是特征線 三類流態(tài)中的特征線 定常均勻流動(dòng) 相容關(guān)系描述的狀態(tài)特征線 特征線 不代表波的傳播跡線 v c0 c c0 簡(jiǎn)單波流動(dòng) 特征線 相容關(guān)系描述的狀態(tài)特征線 活塞運(yùn)動(dòng)跡線 復(fù)合波流動(dòng) 特征線 相容關(guān)系描述的狀態(tài)特征線 x t c c1 v c1 C C 7 6 5 2 3 4 8 9 10 2 3 4 5 6 7 8 9 10 依賴區(qū)和影響區(qū) 由于沿著兩族特征線 分別有 可以把J 和J 看作是兩個(gè)新的函數(shù) 則 利用J 和J 表示的特征線方程為 第I族特征線斜率僅由J 決定 第II族特征線斜率僅由J 決定 x t D A B D點(diǎn)的依賴區(qū) M C C C C 在平面運(yùn)動(dòng)中 沿著特征線黎曼不變量保持不變 這一重要性質(zhì)清楚地揭示出流體動(dòng)力學(xué)中的一些依賴關(guān)系 設(shè)t 0時(shí)各量沿x軸的分布為v0 x c0 x 于是可知黎曼不變量的相應(yīng)分布為 則 x t 平面上任意一點(diǎn)D x t 上的狀態(tài) 將直接由x軸上點(diǎn)A xA 0 B xB 0 兩點(diǎn)上的狀態(tài)決定 所以 點(diǎn)所處的狀態(tài)將完全由且只由線段AB上的值決定 線段AB就稱為點(diǎn)D的依賴區(qū) 同樣 能夠受到AB線段間某點(diǎn)M的初始值影響的區(qū)域 是由發(fā)自M點(diǎn)的與發(fā)自M點(diǎn)的所包圍的區(qū)域 而這個(gè)區(qū)域之外的地方 都不受M點(diǎn)的影響 這個(gè)區(qū)域稱為M點(diǎn)的影響區(qū) P Q 例 已知初始時(shí)刻v x 0 c x 0 求D點(diǎn)的v x t c x t C C 解 在D x3 t 點(diǎn) 有 根據(jù) 得 由 2 3兩個(gè)偏微分方程的特征線法 考慮下面兩個(gè)偏微分方程組成的方程組 x y是自變量 u x y 和v x y 是兩個(gè)因變量 系數(shù)A B及非齊次項(xiàng)F可以是x y u和v的函數(shù) 方程組是準(zhǔn)線性的 以上兩個(gè)方程進(jìn)行線性組合 假設(shè)待求函數(shù)u x y 和v x y 在x y平面上是連續(xù)的 則連續(xù)函數(shù)的全微分為 上式作對(duì)比 可以發(fā)現(xiàn) 若存在一條斜率為下式的平面曲線 沿著該曲線 偏微分方程就化為全微分方程 1 特征線方程 式可化為 為使得關(guān)于 1 2的方程組有非零解 系數(shù)行列式為零 即 化簡(jiǎn)行列式得 其中 對(duì)應(yīng)于橢圓型方程 沒有實(shí)數(shù)解 對(duì)應(yīng)于拋物型方程 過每一點(diǎn)有一條特征線 對(duì)應(yīng)于雙曲型方程 過每一點(diǎn)有兩條特征線 物理特征線方程 對(duì)于一個(gè)一階的偏微分方程總是可以用特征線法求解 但是對(duì)于兩個(gè)一階的偏微分方程組來說 只有雙曲型方程才能利用兩條特征線求出兩個(gè)因變量的數(shù)值解 2 相容方程 由 解出 代入全微分方程 求得相容方程 2 4初值 Cauchy 問題 兩個(gè)偏微分方程的特征線數(shù)值解法 p xp yp C C F G MN是物理平面上一條不是特征線的曲線 沿著該線各點(diǎn)的x y和u v都是已知的 求此曲線鄰域內(nèi)的解 1 先確定F點(diǎn)的位置 由C 和C 的特征線方程 求得F點(diǎn)的位置 2 求F點(diǎn)處的因變量值 上式中包含的 3 由F點(diǎn)處的因變量值 將代入特征線方程 重新計(jì)算特征線方程中的系數(shù) 重復(fù)1 2 過程 重新計(jì)算過M點(diǎn)的C 和過N點(diǎn)的C 兩特征線的坐標(biāo) 反復(fù)迭代 一直計(jì)算到滿足精度為止 上述過程重復(fù)進(jìn)行 從而得到一條新的初值線 再沿著新初值線重復(fù)下一輪運(yùn)算過程 一直可以計(jì)算到初值線的AB與K包圍的區(qū)域 A N y x C C F G M B K G H 當(dāng)已知函數(shù)的初值線 則可以沿著x y平面上M點(diǎn)的特征線的方向 用常微分方程組求解u v兩個(gè)曲面的函數(shù)值 而不是沿任意方向用偏微分方程組求解u v兩個(gè)曲面 由于在每條特征線上各有自己的相容性方程 而每個(gè)相容性方程中又有du dv兩個(gè)函數(shù)的微分 所以單個(gè)相容性方程無法求解 但任意點(diǎn)p xp yp 上有兩條特征線到達(dá) 其上的函數(shù)全微分du du dv dv 雖然沿著不同的特征線發(fā)展 其終值up vp卻是同一個(gè) 因此經(jīng)過P點(diǎn)的兩條特征線上的各一個(gè)相容方程可以聯(lián)立求解 總結(jié) A M N A點(diǎn)的依賴區(qū) B B點(diǎn)的影響區(qū) x y y x 初值線 初值線 D 平面二維的依賴區(qū)和影響區(qū) 1 擾動(dòng) 當(dāng)流場(chǎng)中的一個(gè)區(qū)域 由于物體運(yùn)動(dòng) 物面轉(zhuǎn)折或炸藥爆炸等原因使氣流參數(shù)發(fā)生變化 破壞了原來的平衡狀態(tài)時(shí) 即為氣體受到了擾動(dòng) 2 波 氣體的擾動(dòng)都是以波的形式向流場(chǎng)各處傳播的 在超聲速流場(chǎng)中 在某處使氣體膨脹或者壓縮的任何擾動(dòng)都是通過等熵波 連續(xù)波 或激波 間斷波 傳播到流場(chǎng)一定范圍內(nèi) 2 5二維定常超音速無旋流動(dòng)的特征線解法 3 弱擾動(dòng)波 壓縮擾動(dòng) p 0 膨脹擾動(dòng) p 0 2 5 1弱擾動(dòng)波的一維傳播 定常問題 非定常問題 參考坐標(biāo)系 選取與弱擾動(dòng)波一起運(yùn)動(dòng)的坐標(biāo)系 音速 非定常流動(dòng) 定常流動(dòng) 弱擾動(dòng)波相對(duì)于波前氣體的傳播速度為音速 x正方向 控制體 擾動(dòng)區(qū) 未擾動(dòng)區(qū) 連續(xù)方程 1 動(dòng)量方程 2 由 1 和 2 得 證明 弱擾動(dòng)的傳播過程為等熵過程 由于弱擾動(dòng)的傳播過程很快 可以認(rèn)為是絕熱過程 由絕熱可壓流體的能量方程 有 去掉髙階小量 得 根據(jù)由比熱焓表示的熱力學(xué)第一定律 得 由 2 式 得 因此 弱擾動(dòng)的傳播過程是等熵過程 由完全氣體的等熵方程 得到 對(duì)T 288K的空氣 流體中的音速是氣體介質(zhì)狀態(tài)參數(shù)的函數(shù) 在相同的溫度下 不同介質(zhì)有不同的音速 在同一氣體中 音速隨著氣體溫度升高而升高 并與氣體的熱力學(xué)溫度的平方根成比例 音速是弱擾動(dòng)波相對(duì)于波前氣體的傳播速度 音速的特性 馬赫數(shù) 氣體在某點(diǎn)的流速與當(dāng)?shù)匾羲僦?M1超音速流 弱擾動(dòng)波傳播的絕對(duì)速度 v 0 兩道弱擾動(dòng)波向上游和下游傳播速度均為c vc 兩道弱擾動(dòng)波均向下游傳播 馬赫數(shù) 流體力學(xué)中表征流體壓縮性影響的相似準(zhǔn)數(shù) 為紀(jì)念E 馬赫而命名 馬赫數(shù)表示作用于流體微團(tuán)的慣性力與彈性力之比 在不可壓縮流動(dòng)中 流體密度不變 聲速為無限大 馬赫數(shù)為零 在可壓縮流動(dòng)中 馬赫數(shù)越大 流體的密度變化越大 即流體表現(xiàn)出的可壓縮性越大 通常 按不同的馬赫數(shù)范圍 工程上常把流動(dòng)劃分為低速流動(dòng) M 0 3 亞聲速流動(dòng) 0 3 M 0 8 跨聲速流動(dòng) 0 8 M 1 2 超聲速流動(dòng) 1 2 Ma 5 和高超聲速流動(dòng) M 5 等 馬赫數(shù)的性質(zhì) 2 5 2微擾動(dòng)在空間的傳播特征 擾動(dòng)源靜止 氣流速度對(duì)擾動(dòng)傳播特性的影響 擾動(dòng)波波形 擾動(dòng)中心 ct ct 擾動(dòng)不能超越擾動(dòng)源向前傳播 擾動(dòng)波集中在線的一側(cè) 超音速氣流中擾動(dòng)集中在馬赫錐內(nèi) 馬赫角 vt 馬赫錐半頂角 馬赫角 當(dāng)擾動(dòng)源和氣體間的相對(duì)速度不同時(shí) 波面的傳播有以下四種情況 1 無相對(duì)運(yùn)動(dòng) v 0 擾動(dòng)源靜止 即擾動(dòng)源運(yùn)動(dòng)速度v 0 波面為一系列的同心球面 球心就是擾動(dòng)波源所在的位置 一定時(shí)間后 將傳播到整個(gè)空間 2 擾動(dòng)源以亞音速 vc 運(yùn)動(dòng) 擾動(dòng)源始終處于其發(fā)出的擾動(dòng)波陣面之前 傳播范圍為一圓錐形空間 馬赫錐 馬赫錐頂角之半為馬赫角 錐面即為馬赫波 1 二維定常等熵流動(dòng)基本方程 連續(xù)方程 2 5 3特征線法求解二維定常超音速無旋流動(dòng) 動(dòng)量方程 能量方程 等熵狀態(tài)方程 連續(xù)方程中的兩項(xiàng)可化簡(jiǎn)為 動(dòng)量方程 動(dòng)量方程 兩式相加 以及無旋流條件 整理得 與標(biāo)準(zhǔn)線性偏微分方程進(jìn)行比較 有 則得 2 特征線方程 當(dāng) 即 所以特征線方法能用于解超音速條件下的平面二維定常無旋流動(dòng) 特征線存在 可用特征線方法求解 3 相容方程 4 速度平面上的特征線 超音速定常無旋流動(dòng) 設(shè)速度為V 氣流方向角為 則速度分量為 表示速度和x軸正方向的夾角 則特征線方程為 2020 1 30 77 可編輯 則相容方程為 沿C 特征線 沿C 特征線 3 4變截面等熵管流 2 5 4氣流速度與通道截面的關(guān)系 1 基本方程 微分形式的連續(xù)方程 微分形式的動(dòng)量方程 微分形式的氣體狀態(tài)方程 積分形式的能量方程 2 截面變化造成的影響 由動(dòng)量方程 得 結(jié)合連續(xù)方程 得 截面變化與速度變化的關(guān)系 3 三種流動(dòng)情況 a 亞音速流動(dòng) M 1 dv和dA的符號(hào)相反 截面積縮小 速度增加 截面積擴(kuò)大 速度減小 c 等音速流動(dòng) M 1 無論何種類型的流動(dòng) M 1處的截面積具有極小值 該截面為臨界截面 臨界截面一定是管道的最小截面 但最小截面不一定是臨界截面 b 超音速流動(dòng) M 1 dv和dA的符號(hào)相同 截面積擴(kuò)大 速度增加 截面積減小 速度減小 截面積變化對(duì)流動(dòng)參數(shù)的影響 2 5 5超音速氣流遇壁面外折引起膨脹馬赫波 B 是極微小的角度 則O點(diǎn)相當(dāng)于一個(gè)弱擾動(dòng)源 擾動(dòng)的傳播范圍是在由O點(diǎn)發(fā)出的馬赫波OL的下游 擾動(dòng)的影響是使氣流外折 OL與原始?xì)饬鲓A角是馬赫角 原始?xì)饬鞯搅薕L處感受到壁面外折的影響 方向折轉(zhuǎn)角 沿著OB壁面的方向流動(dòng) 相當(dāng)于放寬氣流的通道 dA 0 經(jīng)過膨脹波以后 氣流參數(shù)的變化趨勢(shì)怎么樣 首先 流速V是不斷增大的 即 由微分形式動(dòng)量方程 由絕能流的能量方程 由狀態(tài)方程 1 d 1 M1 1 L1 3 4 2 L2 L4 L3 在點(diǎn)處 氣流受到O1L1的擾動(dòng) 氣流折轉(zhuǎn)角度為 速度變?yōu)?氣流在處受外折微小角度以后 又在和繼續(xù)外折角度和 在點(diǎn)處 氣流受到O1L1的擾動(dòng) 氣流折轉(zhuǎn)角度為 速度變?yōu)?由于 在點(diǎn)處 氣流受到O1L1的擾動(dòng) 氣流折轉(zhuǎn)角度為 速度變?yōu)?所以 后產(chǎn)生的膨脹波相對(duì)于原始?xì)饬鞯膬A斜角都比前一道的小 膨脹波不可能彼此相交 因而形成一個(gè)連續(xù)的膨脹區(qū)域 根據(jù)極限的概念 曲線可以看作是無數(shù)條微元折線的極限 因而 超聲速氣流繞外凸曲壁膨脹加速的情況與上面分析完全一樣 只是單個(gè)的膨脹波連成連續(xù)的膨脹波了 也稱作 膨脹馬赫波 壁面從連續(xù)外折到點(diǎn) 經(jīng)過無限多次折轉(zhuǎn)角后 總折角為 1 2 3 4 5 x y 為壁面和水平方向的夾角 1為馬赫波與壁面的夾角 注意 正負(fù)號(hào) M1 M3 M4 M2 M5 左伸膨脹馬赫波 超音速氣流沿連續(xù)下彎壁面的流動(dòng) 正負(fù)號(hào)規(guī)定 由x正方向逆時(shí)針旋轉(zhuǎn)為正 順時(shí)針旋轉(zhuǎn)為負(fù) 波線方程 外凸的壁面上方 形成膨脹波 即 M1 M2 M3 M4 M5 隨著流線折轉(zhuǎn) 連續(xù)膨脹馬赫波是發(fā)散的 x y 1 2 3 4 超音速氣流沿連續(xù)上彎壁面的流動(dòng) 右伸膨脹馬赫波 波線方程 左伸馬赫波 波線方程為 跨波線 特征線 波線方程為 右伸馬赫波 跨波線 特征線 Prandtl Meyer 普朗陀 邁耶 流動(dòng) 如果壁面彎曲段縮成一個(gè)點(diǎn) 氣流流過如圖所示的外凸壁時(shí) 可以看作由一系列折轉(zhuǎn)無限小的外凸壁的流動(dòng) 氣流每折轉(zhuǎn)一個(gè)角度 就產(chǎn)生一道膨脹波 而氣流每經(jīng)過一道膨脹波 馬赫數(shù)增大 馬赫角減小 因此 這些膨脹波發(fā)散 如果壁面的幾個(gè)折轉(zhuǎn)點(diǎn)都無限接近于點(diǎn)O1 就形成了普朗特 邁耶流動(dòng) 普朗陀 邁耶流動(dòng)的形成 O1 2 6簡(jiǎn)單波 1 齊次可約方程組 當(dāng)上面偏微分方程組中的 稱方程組為齊次方程組 當(dāng)齊次方程組的系數(shù)都只依賴于u v 則方程組稱為是可約方程組 齊次可約方程組的相容方程為 由于右端各項(xiàng)均與x y無關(guān) 只與u v有關(guān) 上式可以獨(dú)立積分 從而得到方程組的特解 2 簡(jiǎn)單波解的一般形式 為了更好地了解簡(jiǎn)單波的性質(zhì) 我們需要求簡(jiǎn)單波的一般形式的解 由一維等熵流的基本方程組 對(duì)于 基本方程組中的第二個(gè)方程自動(dòng)滿足 第一個(gè)方程化為 的簡(jiǎn)單波流動(dòng) 由偏微分方程形式 猜想 因變量可能為兩個(gè)自變量某組合形式的函數(shù) 且 發(fā)現(xiàn)方程左邊為 0 猜想正確 表明因變量v可看做X的函數(shù) 則得通解為 是自變量 的任意函數(shù) 由問題的邊界條件決定 或 對(duì)于 的簡(jiǎn)單波流動(dòng) 其通解為 或 以上解表明 簡(jiǎn)單波就是向一個(gè)方向傳播的波 或者說簡(jiǎn)單波是單向行波 沿著特征線 前面由特征線方法求解曾得到 C 族特征線 沿著特征線 C 族特征線 由簡(jiǎn)單波通解得到 當(dāng) 跡線與C 族特征線重合 或 C 族特征線是跡線 時(shí) 由簡(jiǎn)單波通解得到 當(dāng) 跡線與C 族特征線重合 或 C 族特征線是質(zhì)點(diǎn)跡線 I 中心簡(jiǎn)單波 當(dāng)邊界條件使得時(shí) 得到的解為 或 中心簡(jiǎn)單波 為特征線的共同起點(diǎn) 稱為中心點(diǎn) II 向前 右傳 和向后 左傳 簡(jiǎn)單波 由于簡(jiǎn)單波時(shí)單向行波 按傳播方向分為向前和向后簡(jiǎn)單波 若波的傳播速度大于質(zhì)點(diǎn)速度 則流體質(zhì)點(diǎn)將從右側(cè)進(jìn)入波動(dòng)區(qū) 這種簡(jiǎn)單波稱為向前 右傳 簡(jiǎn)單波 反之 稱為向后 左傳 簡(jiǎn)單波 的波是向前 右傳 簡(jiǎn)單波 的波是向后 左傳 簡(jiǎn)單波 傳播速度為 傳播速度為 III 壓縮波和稀疏波 穿過簡(jiǎn)單波后 若流體的密度和壓強(qiáng)增大 稱這波為壓縮簡(jiǎn)單波 簡(jiǎn)稱壓縮波 若流體的密度和壓強(qiáng)減小 則稱為稀疏波 活塞背離氣體移動(dòng) 例1 設(shè)一無限長(zhǎng)管道的左半段充滿高壓氣體 右半段是真空 兩段間用一薄膜相隔 高壓氣體是完全氣體 初始狀態(tài)為 求拆除薄膜后氣體的運(yùn)動(dòng) 0 解 將氣體與真空的交界面初始位置取為坐標(biāo)原點(diǎn)x 0 當(dāng)t 0時(shí)拆除薄膜的一瞬間 氣體界面將被加速到某一個(gè)速度向右運(yùn)動(dòng) 同時(shí) 向高壓氣體內(nèi)傳入一個(gè)向后稀疏波 則波的傳播跡線方程與族特征線重合 其方程為 結(jié)合初始條件 上式積分為 沿特征線 滿足 1 2 3 聯(lián)立兩式 得波內(nèi)的流場(chǎng)分布 界面氣體質(zhì)點(diǎn)的飛散速度 4 5 界面與真空相毗鄰的是一個(gè)自由面 該自由面上滿足 由 得 將上式代入 中 得界面的飛散速度為 6 7 例2 稀疏波的解 由充滿氣體的管道中抽離活塞產(chǎn)生的管內(nèi)氣體運(yùn)動(dòng) 是稀疏波的一個(gè)典型例子 t 0時(shí)刻向左抽動(dòng)活塞 活塞速度為w t at 設(shè)活塞運(yùn)動(dòng)到A點(diǎn)處時(shí) 速度達(dá)到uA時(shí) 活塞不再加速 而是以常速uA繼續(xù)向左運(yùn)動(dòng) 分析 圖中 0 區(qū)域是尚未受擾動(dòng)的常態(tài)區(qū) 區(qū)域 I 是簡(jiǎn)單波區(qū) 在活塞軌跡上A點(diǎn)之后活塞速度為常數(shù) 該段的特征線是一族平行線 出現(xiàn)對(duì)應(yīng)的一個(gè)區(qū)域 II 它與區(qū)域 I 的分界線是發(fā)自A點(diǎn)的一條特征線 1 區(qū)域 I 內(nèi)的解 對(duì)完全氣體 沿著x軸 1 根據(jù)什么 根據(jù)什么 有 如何求 思考 2 聯(lián)立 邊界條件 根據(jù)活塞的邊界條件確定任意函數(shù) 活塞的軌跡為 在氣體邊界上滿足 緊靠活塞的氣體其位置和速度與活塞相同 即 3 4 將邊界條件以及代入 2 式 5 的函數(shù)形式 將 代入 2 式 得 經(jīng)整理 得 6 7 1 由于在波頭x c0t上氣體的速度應(yīng)v 0 所以取了根號(hào)前為 號(hào)的根 2 區(qū)域I中的坐標(biāo)x c0t 因此v 0 即波內(nèi)氣體都隨活塞向同一方向運(yùn)動(dòng) 求解關(guān)于的二次代數(shù)方程 解得 8 區(qū)域 I 中的C 族特征線是一族直線 其方程為 對(duì)上式求積分 積分的起點(diǎn)在活塞上 9 10 由初始條件 代入 11 式 于是得到對(duì)應(yīng)于初始時(shí)刻的的方程為 11 當(dāng)活塞在到達(dá)A點(diǎn) 速度達(dá)到 注 指右活塞開始加速運(yùn)動(dòng)時(shí)為起點(diǎn)發(fā)出的第一條特征線 指右活塞開始勻速運(yùn)動(dòng)時(shí)發(fā)出的一條特征線 是II區(qū)和I區(qū)的交界線 2 區(qū)域 II 內(nèi)的解 區(qū)域II中的族特征線來自區(qū)域0 故整個(gè)區(qū)域II內(nèi)有 12 上式在活塞軌跡線上也成立 II區(qū)活塞軌跡線上有 則得活塞軌跡線上氣流聲速為 13 則沿活塞軌跡線上的黎曼不變量為 則整個(gè)區(qū)域II中的黎曼不變量相等 則聯(lián)立兩式可得 14 15 從而解出區(qū)域II中的解 16 為活塞勻速運(yùn)動(dòng)的速度 17 區(qū)域II的族特征線即為區(qū)域I中的特征線的平行線 3 逃逸速度 活塞跡線 氣體微團(tuán)跡線 亞音速流動(dòng)區(qū) 超音速流動(dòng)區(qū) 當(dāng)活塞由靜止連續(xù)地向左加速 產(chǎn)生一族向右稀疏波 跨過右傳波 滿足 在活塞和氣體不分離的條件下 與活塞毗鄰的氣體速度等于活塞運(yùn)動(dòng)速度 與活塞毗鄰的氣體音速 18 19 隨著活塞向左的速度增加 音速增大還是減小 族特征線 波線 方程為 當(dāng) 即從活塞上會(huì)產(chǎn)生一道平行于t軸的駐波 20 21 在駐波上當(dāng)?shù)貧饬魉俣?絕對(duì)值 等于當(dāng)?shù)匾羲?記做 駐波將簡(jiǎn)單稀疏波區(qū)分成兩個(gè)區(qū)域 亞音速流動(dòng)區(qū)和超音速流動(dòng)區(qū) 繼續(xù)增大活塞速度 22 使音速降到零值 對(duì)應(yīng)的活塞速度為 與之毗鄰的氣體速度為 23 逃逸速度是氣體通過稀疏波膨脹所能達(dá)到的最大極限速度 之后即使繼續(xù)加大活塞速度 氣體不可能繼續(xù)加速 因?yàn)闅怏w與活塞從此分離 它們之間形成真空 實(shí)際上逃逸速度不可能達(dá)到 例3 無限長(zhǎng)管道內(nèi)高壓氣體推動(dòng)剛體運(yùn)動(dòng)的解 設(shè)在截面積為A管道內(nèi)x 0截面處有一剛體 其質(zhì)量為M0 剛體右側(cè)是真空 左側(cè)充滿高壓氣體 其初始狀態(tài)為 剛體在高壓氣體的推動(dòng)下在時(shí)開始向右運(yùn)動(dòng) 求剛體松開后的運(yùn)動(dòng)軌跡 解 剛體松開后 在高壓氣體的作用下向右運(yùn)動(dòng) 從而有一左傳稀疏波 使得剛體上的壓力隨時(shí)間發(fā)生變