黎曼间断.docVIP

一．Riemann间断解 1.问题的提出 一维非定常理想气体守恒方程 取t＝0时刻的初始值： 求方程（1）满足以上初值条件且在间断线上满足间断条件的解。 间断条件： 2.解的分类 一般情况下，给定的初始间断并不满足关系式（2），因而它是不稳定的，当后，立刻分解为若干满足间断条件的间断解，并以各自的速度分别运动 初始间断的所分解的解，在考虑熵条件下，是唯一确定的。这些准确解，称为Riemann间断解。一般分为五类 图 虚线表示接触间断，实线表示激波，一簇中心直线表示稀疏波(膨胀波)，（e）表示两个稀疏波之间为真空。 时，左右波尚未达到的波前区域内，流动状态仍保持常数分布： 接触间断两边的速度和压力值相同，由表示，密度值不同分别为。 中心稀疏波是弱间断，他反映了气流的等熵膨胀过程。在此区域内，物理量由波前状态连续过渡到波后状态且有，只有物理量的导数在波头和波尾处是间断的。 ④ 对于激波，物理量是由波前状态跳跃到波后状态，熵增，且。 3.求解法 理想气体状态方程变为： 为比热比。 等熵关系式： 声速： 假设左波为一向左运动激波，速度为，激波左侧和右侧气流的参数分别为，。若认为激波是相对静止的，那么左、右侧气 流的相对速度：…………………… (3) 间断条件： ………………………………………………… (4) ………………………………………… (5) ………………………………………… (6) 根据(4)式引入质量流量 ……………………………………………….. (7) 那么(5)式为 ………………………………………. (8) 将(7)式代入(8)式 ………………. (9) 得 ……………………………………….. (10) 将代入(9)式 ………………………………………………(11) 激波前后的密度关系（朗金—雨贡纽方程） ，代入(11)式，得 ……………………………………(12) 同理，对于右激波可得： ……………………………………(13) 假设左波为一中心稀疏波，满足等熵关系式 ，又，……….(14) 在中心稀疏波内，由于Riemann不变量沿特征线保持常数，故有 ….…….(15) 同理，对于右稀疏波也有： …………………………….(16) 将（10）（13）（15）（16）统一写成： 左波： 右波： 其中： 合并(17)式和(18)式，消去U，得： 为了对给定的初始条件判断其间断分解的类型，下面讨论一下函数的性质。 函数在处是连续的，并且有连续的一阶导数。 当时，，那么 当时，是单调上升的凹函数。即有，（19）式，。因此在时，也是单调上升，导数连续的凹函数。由这一性质，可以给出间断分解类型的判别式。 先假设 1. 由可推得。 即左波、右波都是激波，左激波向左边低压区域传播()，右激波向右边低压区域传播()。 左右激波之间的区域，压力连续，速度为U，密度是间断的——该区域内有一接触间断（切向间断）。 2. 可以推出，左为激波，右为稀疏波。 激波与稀疏波之间区域内，由于气流激波后压力升高，而稀疏波后压力降低，压力和速度U连续，密度间断。 3. 可推得：，左波和右波都是稀疏波。 4. 左右均为稀疏波，且在稀疏波之间出现真空区。真空区内各参数值为零。 若假设，与上面的推导类似。 时， ，右波为激波，左波为稀疏波。 参数确定：由（10）式可以推出，此式为超越方程，需要迭代求解，求得后，即可求得U。 由此，一维非定常理想气体在平面上的间断解，，可根据解的类型及的位置分别由以下公式给出。 1.左波区域 当左波为激波时，激波速度和波后密度 ， 左激波的波后区域：，取为。 区域为未扰动的区域，取为。 当左波为稀疏波时，波头速度为和波尾速度可由特征线的表达式给出： 波后的密度可由等熵关系式推得，而可由Riemann不变量求得。 波前和波后：为 为 稀疏波区：根据Riemann不变量关系式： 由特征线，则可以确定声速c为： 由此可以得到： 2. 右波区域 当右波为激波时，激波速度和波后密度 右激波波前区：区域，取为。 右激波波后区：，取为。 当右波为稀疏波时，波头速度和波尾速度可由特征线的表达式给出： 波后的密度可由等熵关系式推得，而可由Riemann不变量求得。 波前区： 为 波后区： 为 稀疏波区：根据Riemann不变量关系式： 由特征线，则可以确定声速c为： 由此可以得到： 3. 出现真空区时 在右稀疏波后压力，密度，且星形区域内速度U已经没有意义，稀疏波波尾速度计算公式不再成立，改写为： 稀疏波区内的流动参数计算方法不变。