特征线理论及应用.ppt

《特征线理论及应用.ppt》由会员分享，可在线阅读，更多相关《特征线理论及应用.ppt（151页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第二章第二章 特征线理论及应用特征线理论及应用气体动力学中，有大量问题是用双曲型偏微分方程来描述的，气体动力学中，有大量问题是用双曲型偏微分方程来描述的，很难得到解析结果，在这种情况下，有两种数值解法：很难得到解析结果，在这种情况下，有两种数值解法：1）特征线数值解法特征线数值解法：求解域用特征线网格进行离散，求各网：求解域用特征线网格进行离散，求各网格结点上的解；气体动力学中，有大量流动问题是用双曲型格结点上的解；气体动力学中，有大量流动问题是用双曲型偏微分方程来描述的，宜于用特征线方法求解。偏微分方程来描述的，宜于用特征线方法求解。2）有限差分法有限差分法：求解域的有限差分网格一般是正交的

2、，根据：求解域的有限差分网格一般是正交的，根据由偏微分方程构造的差分格式来求各网格结点上的解。由偏微分方程构造的差分格式来求各网格结点上的解。2.1 特征线理论特征线理论 特征线的数学定义特征线的数学定义考虑一个一般的一阶双曲型偏微分方程：考虑一个一般的一阶双曲型偏微分方程：x,y 是两个自变量，是两个自变量，u(x,y)是因变量。系数是因变量。系数A1、A2及及非齐次项非齐次项F1可以是可以是 x，y，u 的函数。的函数。（1）将偏微分方程改写为：将偏微分方程改写为：设未知函数设未知函数u(x,y)连续，连续，u 的一阶导数可以写作：的一阶导数可以写作：【注注:u的一阶导数可以不连续的一阶导

3、数可以不连续】偏微分方程的特征线定义为：xy平面内具有斜率为平面内具有斜率为 的曲线。的曲线。（2）（3）沿着特征线或：偏微分方程可化简为：代入代入 式式（4）得到偏微分方程的相容方程【是平面上这样一族曲线：沿着此族中任一曲线(a)，可以把待求物理量的一阶偏微分控制方程变换成等价的常微分控制方程(b)，称为原偏微分方程或偏微分方程组的相容方程】u 特征线的第一个数学意义：特征线的第一个数学意义：(a)(b)u特征线的第二个数学意义：特征线的第二个数学意义：上两式表明：上两式表明：沿着特征线，分母和分子均为零。沿着特征线，分母和分子均为零。即沿着特征线，即沿着特征线，表明：表明：1 1）沿特征线

4、因变量的一阶导数具有不定值，可以是不连续的，在这种情）沿特征线因变量的一阶导数具有不定值，可以是不连续的，在这种情况下，特征线是弱间断（第一类间断线）。况下，特征线是弱间断（第一类间断线）。2 2）在气体动力学中，特征线可以是弱扰动波传播的迹线，或者说弱扰动）在气体动力学中，特征线可以是弱扰动波传播的迹线，或者说弱扰动传播的迹线就是特征线。传播的迹线就是特征线。因此，因变量的一阶导数只允许有弱间断，如果在物理平面上有激因此，因变量的一阶导数只允许有弱间断，如果在物理平面上有激波出现，在强间断面上便无法建立因变量的全微分式，也就不能用特征线波出现，在强间断面上便无法建立因变量的全微分式，也就不能

5、用特征线方法求解。方法求解。例：一阶偏微分方程的初始条件是的初始条件是2）沿此特征线的相容方程）沿此特征线的相容方程3）u(2,4)的值的值用特征线法确定：用特征线法确定：1）通过点）通过点(2,4)的特征线的特征线解解：（1）对照一般形式的双曲型偏微分方程对照一般形式的双曲型偏微分方程该方程对应的系数：该方程对应的系数：A1=1,A2=2x,F1=3x2则特征线方程为：则特征线方程为：积分得：积分得：为确定过点为确定过点(2,4)的特征线，将的特征线，将x=2,y=4，代入上式得：，代入上式得：所以，所求的特征线方程是：所以，所求的特征线方程是：对上式积分，得：对上式积分，得：（2）偏微分方

6、程的相容方程为：偏微分方程的相容方程为：如何确定 C2？初始条件初始条件 u(0,y)=5y+10及特征线方程及特征线方程u(0,0)=10因此相容方程为：因此相容方程为：2.2 2.2 一维等熵流动的特征线数值解法一维等熵流动的特征线数值解法 基本方程与黎曼不变量基本方程与黎曼不变量（连续方程）（连续方程）（动量方程）（动量方程）（以一维等直截面管为例）（以一维等直截面管为例）基本方程基本方程等熵流动中只有一个等熵流动中只有一个状态参量状态参量独立：独立：将基本方程中的将基本方程中的 用用 代替，得：代替，得：基本方程可化为：基本方程可化为：两式相加减两式相加减合并，基本方程可写作：合并，基

7、本方程可写作：定义定义则基本方程化为以基本方程化为以v G 为新的未知函数的偏微分方程：为新的未知函数的偏微分方程：基本方程基本方程偏微分方程偏微分方程特征线特征线？相容方程？相容方程？在x-t 平面上，把dx/dt=v c 曲线称为偏微分方程的特征线。CCxtC表示第一族特征线表示第一族特征线；C 表示第二族特征线。表示第二族特征线。解相容方程解相容方程:对多方气体：对多方气体：其相容方程的解为：其相容方程的解为：由由声速：声速：沿着特征线沿着特征线结论：结论：特征线的基本性质特征线的基本性质1）一维非定常流动中，平面一维非定常流动中，平面x-t上任一点，都有两条不同上任一点，都有两条不同族

8、的特征线，沿各特征线有各自不同的黎曼不变量；族的特征线，沿各特征线有各自不同的黎曼不变量；2）特征线上参量特征线上参量v,c,p,的一阶导数可以不连续，但这的一阶导数可以不连续，但这些参量本身是连续的，称因变量的一阶导数不连续的点叫些参量本身是连续的，称因变量的一阶导数不连续的点叫做弱间断。如果初始某一点有弱间断，那么这个弱间断必做弱间断。如果初始某一点有弱间断，那么这个弱间断必定会沿着过该点的特征线向外传播定会沿着过该点的特征线向外传播。3）两个相邻的，不同类型流动区域的分界线，必定是特两个相邻的，不同类型流动区域的分界线，必定是特征线。征线。三类流态中的特征线三类流态中的特征线vcc0v0

9、定常均匀流动定常均匀流动相容关系描述的状态特征线相容关系描述的状态特征线xt特征线特征线（不代表波的传播迹线）（不代表波的传播迹线）xt(0)(I)(II)v/c0c/c0(0)(1)(2)(3)(4)简单波流动简单波流动特征线特征线相容关系描述的状态特征线相容关系描述的状态特征线活塞运动迹线活塞运动迹线复合波流动复合波流动特征线特征线相容关系描述的状态特征线相容关系描述的状态特征线xtc/c1v/c1C+C-7652348910234 5678910 依赖区和影响区依赖区和影响区由于沿着两族特征线，分别有：由于沿着两族特征线，分别有：可以把可以把 J 和和 J 看作是两个新的函数，则看作是两

10、个新的函数，则利用 J 和 J表示的特征线方程为：第 I 族特征线斜率仅由 J-决定；第 II 族特征线斜率仅由 J+决定。xtDABD点的依赖区点的依赖区MCC CC 在平面运动中，沿着特征线黎曼不变量保持不变，这在平面运动中，沿着特征线黎曼不变量保持不变，这一重要性质清楚地揭示出流体动力学中的一些依赖关系。一重要性质清楚地揭示出流体动力学中的一些依赖关系。设设 t=0 时各量沿时各量沿x轴的分布为轴的分布为v0(x),c0(x),于是可知黎曼不变量的于是可知黎曼不变量的相应分布为相应分布为 。则则(x,t)平面上任意一点平面上任意一点D(x,t)上上的状态，将直接由的状态，将直接由x轴上点

11、轴上点A(xA,0),B(xB,0)两点上的状态决定。两点上的状态决定。所以，点所以，点 所处的状态将完全由且所处的状态将完全由且只由线段只由线段ABAB上的值决定，上的值决定，线段线段AB就称为点就称为点D的依赖区的依赖区。DMM点的影响区点的影响区xtAB C C 同样，能够受到同样，能够受到ABAB线线段间某点段间某点MM的初始值影响的初始值影响的区域，是由发自的区域，是由发自MM点的点的 与发自与发自MM点的点的 所包所包围的区域，而这个区域围的区域，而这个区域之外的地方，都不受之外的地方，都不受MM点点的影响。这个区域称为的影响。这个区域称为MM点的影响区点的影响区。PQ例例：已知初

12、始时刻：已知初始时刻 v(x,0),c(x,0)D(x3,t)A(x1,0)MB(x2,0)xt,求求D点的点的v(x,t),c(x,t)C+C-解：在D(x3,t)点，有根据：根据：得：得：由由2.3 两个偏微分方程的特征线法两个偏微分方程的特征线法考虑下面两个偏微分方程组成的方程组：考虑下面两个偏微分方程组成的方程组：x,y是自变量，是自变量，u(x,y)和和v(x,y)是两个因变量。系数是两个因变量。系数A、B及非齐次项及非齐次项F可以是可以是 x、y、u和和v的函数，方程组是准线性的。的函数，方程组是准线性的。以上两个方程进行线性组合：以上两个方程进行线性组合：假设待求函数假设待求函数

13、u(x,y)和和v(x,y)在在x,y平面上是连平面上是连续的，则连续函数的全微分为：续的，则连续函数的全微分为：上式作对比，可以发现，若存在一条斜率为下式的平面曲线：上式作对比，可以发现，若存在一条斜率为下式的平面曲线：沿着该曲线，偏微分方程就化为全微分方程：沿着该曲线，偏微分方程就化为全微分方程：1 1）特征线方程）特征线方程式式 可化为：可化为：为使得关于为使得关于1 2的方程组有非零解，的方程组有非零解，系数行列式为零，即：系数行列式为零，即：化简行列式得：化简行列式得：其中，其中，对应于椭圆型方程，没有实数解对应于抛物型方程，过每一点有一条特征线对应于抛物型方程，过每一点有一条特征线

14、对应于双曲型方程，过每一点有两条特征线对应于双曲型方程，过每一点有两条特征线物理特征线方程物理特征线方程 对于一个一阶的偏微分方程总是可以用对于一个一阶的偏微分方程总是可以用特征线法求解特征线法求解;但是对于两个一阶的偏微分方程组来说，但是对于两个一阶的偏微分方程组来说，只有双曲型方程才能利用两条特征线求出只有双曲型方程才能利用两条特征线求出两个因变量的数值解。两个因变量的数值解。2 2）相容方程）相容方程由：由：解出解出代入全微分方程代入全微分方程求得相容方程：求得相容方程：MN y x 2.4 初初 值值(Cauchy)问问 题：题：两个偏微分方程的特征线数值解法两个偏微分方程的特征线数值

15、解法p(xp,yp)C+C-FGMN是物理平面上一条不是特是物理平面上一条不是特征线的曲线，沿着该线各点征线的曲线，沿着该线各点的的x,y和和u,v都是已知的，求此都是已知的，求此曲线邻域内的解。曲线邻域内的解。1）先确定）先确定F点的位置点的位置由由C-和和C+的特征线方程：的特征线方程：求得求得F点的位置点的位置2）求）求F点处的因变量点处的因变量 值值上式中包含的上式中包含的3）由）由F点处的因变量点处的因变量 值，将值，将 代代 入特征线方程，重新计算特征线方程中的系数，重复入特征线方程，重新计算特征线方程中的系数，重复1)2)过程，重新计算过过程，重新计算过M点的点的C和过和过N点的

16、点的C两特征线的坐两特征线的坐 标，反复迭代，一直计算到满足精度为止。标，反复迭代，一直计算到满足精度为止。上述过程重复进行，从而得到一条新的初值线上述过程重复进行，从而得到一条新的初值线 ，再沿，再沿着新初值线重复下一轮运算过程，一直可以计算到初值线着新初值线重复下一轮运算过程，一直可以计算到初值线的的AB与与K包围的区域。包围的区域。AN y x C+C-FGMBKGH 当已知函数的初值线当已知函数的初值线 ，则可以沿着，则可以沿着x,y平面上平面上M点的点的特征线的方向，用常微分方程组求解特征线的方向，用常微分方程组求解u,v两个曲面的函数值，两个曲面的函数值，而不是沿任意方向用偏微分方

17、程组求解而不是沿任意方向用偏微分方程组求解u,v两个曲面。两个曲面。由于在每条特征线上各有自己的相容性方程，而每个相容由于在每条特征线上各有自己的相容性方程，而每个相容性方程中又有性方程中又有du,dv两个函数的微分，所以单个相容性方程无两个函数的微分，所以单个相容性方程无法求解；但任意点法求解；但任意点p(xp,yp)上有两条特征线到达，其上的函数上有两条特征线到达，其上的函数全微分全微分du+,du-,dv+,dv-虽然沿着不同的特征线发展，其终值虽然沿着不同的特征线发展，其终值up,vp却是同一个，因此经过却是同一个，因此经过P点的两条特征线上的各一个相点的两条特征线上的各一个相容方程可

18、以联立求解。容方程可以联立求解。总总 结结AMNA A 点的依赖区点的依赖区BB B 点的影响区点的影响区xyyx初值线初值线初值线初值线D平面二维的依赖区和影响区平面二维的依赖区和影响区1、扰动、扰动：当流场中的一个区域，由于物体运动、物面转折或炸药爆炸等原因使气流参数发生变化，破坏了原来的平衡状态时，即为气体受到了扰动。2、波：、波：气体的扰动都是以波的形式向流场各处传播的。在超声速流场 中，在某处使气体膨胀或者压缩的任何扰动都是通过等熵波(连续波)或激波(间断波)传播到流场一定范围内。2.5 二维定常超音速无旋流动的特征线解法二维定常超音速无旋流动的特征线解法3、弱扰动波、弱扰动波：压缩

19、扰动(p0);膨胀扰动(p0)；活塞右移形成压缩波活塞右移形成压缩波 活塞左移形成膨胀波活塞左移形成膨胀波p1+dpp1-dp2.5.1 弱扰动波的一维传播弱扰动波的一维传播定常问题？非定常问题？定常问题？非定常问题？参考坐标系：选取与弱扰动波一起运动的坐标系音音 速速非定常流动 定常流动：弱扰动波相对于波前气体的传播速度为音速。cc-dvp,Tp+dp,+d,T+dTx 正方向控制体扰动区未扰动区连续方程连续方程：（1）动量方程：动量方程：（2）由（1）和（2），得：证明：弱扰动的传播过程为等熵过程。由于弱扰动的传播过程很快，可以认为是绝热过程。由绝热可压流体的能量方程，有：去掉髙阶小量，得

20、：根据由比热焓表示的热力学第一定律，得：由（2）式，得：因此，弱扰动的传播过程是等熵过程。由完全气体的等熵方程得到：对T T288K288K的空气，流体中的音速是气体介质状态参数的函数流体中的音速是气体介质状态参数的函数。在相同的温度下，不同介质有不同的音速。在相同的温度下，不同介质有不同的音速。在同一气体中，音速随着气体温度升高而升高，并在同一气体中，音速随着气体温度升高而升高，并 与气体的热力学温度的平方根成比例。与气体的热力学温度的平方根成比例。音速是弱扰动波相对于波前气体的传播速度。音速是弱扰动波相对于波前气体的传播速度。音 速 的 特 性：马 赫 数气体在某点的流速与当地音速之比。M

21、 1 超音速流弱扰动波传播的绝对速度：v=0 (两道弱扰动波向上游和下游传播速度均为两道弱扰动波向上游和下游传播速度均为c)vc (两道弱扰动波均向下游传播两道弱扰动波均向下游传播)马 赫 数n 流体力学中表征流体压缩性影响的相似准数。为纪念E.马赫而命名。n 马赫数表示作用于流体微团的惯性力与弹性力之比。n 在不可压缩流动中，流体密度不变，声速为无限大，马赫数为零。在可压缩流动中，马赫数越大，流体的密度变化越大，即流体表现出的可压缩性越大。n 通常，按不同的马赫数范围，工程上常把流动划分为低速流动(M0.3)、亚声速流动(0.3M0.8)、跨声速流动(0.8M1.2)、超声速流动(1.2Ma

22、 cv=cctvt作用区寂静区马赫锥半顶角 马赫角当扰动源和气体间的相对速度不同时，波面的传播有以下四种情况:1）无相对运动（无相对运动（v=0）:扰动源静止（即扰动源运动速度扰动源静止（即扰动源运动速度v=0），波），波面为一系列的同心球面，球心就是扰动波源所在的位置，一定时间后，面为一系列的同心球面，球心就是扰动波源所在的位置，一定时间后，将传播到整个空间。将传播到整个空间。2）扰动源以亚音速（扰动源以亚音速（vc）运动）运动:扰动源始终处于其发出的扰动波阵扰动源始终处于其发出的扰动波阵面之前，传播范围为一圆锥形空间（马赫锥），马赫锥顶角之半为马面之前，传播范围为一圆锥形空间（马赫锥），马

23、赫锥顶角之半为马赫角，锥面即为马赫波赫角，锥面即为马赫波。1 1）二维定常等熵流动基本方程二维定常等熵流动基本方程连续方程：连续方程：2.5.3 特征特征线法求解法求解二维定常超音速无旋流动二维定常超音速无旋流动动量方程：动量方程：能量方程：能量方程：等熵状态方程等熵状态方程连续方程中的两项可化简为：连续方程中的两项可化简为：动量方程动量方程动量方程动量方程两式相加，以及无旋流条件，整理得：两式相加，以及无旋流条件，整理得：与标准线性偏微分方程进行比较，有：与标准线性偏微分方程进行比较，有：则得：则得：2 2）特征线方程）特征线方程当当即即所以特征线方法能用于解所以特征线方法能用于解超音速超音

24、速条条件下的平面二维定常无旋流动。件下的平面二维定常无旋流动。,特征线存在；特征线存在；,可用特征线方法求解。可用特征线方法求解。3 3）相容方程相容方程4 4）速度平面上的特征线）速度平面上的特征线超音速定常无旋流动，设速度为超音速定常无旋流动，设速度为V，气流方向，气流方向角为角为，则速度分量速度分量为：表示速度和表示速度和x x轴正方向的夹角轴正方向的夹角则特征线方程为：则相容方程为：沿沿C 特征线特征线沿沿C 特征线特征线3.4 变截面等熵管流变截面等熵管流2.5.4 气流速度与通道截面的关系气流速度与通道截面的关系(1)基本方程基本方程微分形式的连续方程微分形式的连续方程微分形式的动

25、量方程微分形式的动量方程微分形式的气体状态方程微分形式的气体状态方程积分形式的能量方程积分形式的能量方程(2)(2)截面变化造成的影响截面变化造成的影响由动量方程，得：结合连续方程，得：截面变化与速度变化的关系(3)(3)三种流动情况三种流动情况a.亚音速流动(M1):dv 和dA 的符号相同。截面积扩大，速度增加；截面积减小，速度减小。V 增V 减V 增V 减截面积变化对流动参数的影响截面积变化对流动参数的影响dA0气流参数比气流参数比M1M1dv/v000dM/M000dp/p000d/000dT/T001 LO2.5.5 2.5.5 超音速气流遇壁面外折超音速气流遇壁面外折 引起引起膨胀

26、马赫波膨胀马赫波B 是极微小的角度，则O点相当于一个弱扰动源弱扰动源，扰动的传播范围是在由O点发出的马赫波OL的下游，扰动的影响是使气流外折气流外折 。OL与原始气流夹角是马赫角 。原始气流到了原始气流到了OL处感受到壁面外折的影响，方向折转处感受到壁面外折的影响，方向折转 角，沿着角，沿着OB壁面的壁面的方向流动，相当于放宽气流的通道方向流动，相当于放宽气流的通道,dA 0。经过膨胀波以后，气流参数的变化趋势怎么样？经过膨胀波以后，气流参数的变化趋势怎么样？首先，流速首先，流速V V 是不断增大的，即是不断增大的，即由微分形式动量方程：由微分形式动量方程：由绝能流的能量方程：由绝能流的能量方

27、程：由状态方程：由状态方程：1d1M11 L1 3 4 2 L2 L4 L3在在 点处点处，气流受到O1L1的扰动，气流折转角度为：速度变为：气流在 处受外折微小角度 以后，又在 和 继续外折角度 和 。在在 点处点处，气流受到O1L1的扰动，气流折转角度为：速度变为：由于在在 点处点处，气流受到O1L1的扰动，气流折转角度为：速度变为：所以：后产生的膨胀波相对于原始气流的倾斜角都比前一道的小，膨胀波不可能彼此相交，因而形成一个连续的膨胀区域。根据极限的概念，曲线可以看作是无数条微元折线的极限，因而，超声速气流绕外凸曲壁膨胀加速的情况与上面分析完全一样，只是单个的膨胀波连成连续的膨胀波了,也称

28、作“膨胀马赫波”壁面从 连续外折到 点，经过无限多次折转 角后，总折角为：1 2 3 4 5xy为壁面和水平方向的夹角；1为马赫波与壁面的夹角，注意正负号。M1M3M4M2M5左伸膨胀马赫波左伸膨胀马赫波超音速气流沿连续下弯壁面的流动正负号规定：由正负号规定：由x正方向逆时针旋转为正，顺时针旋转为负。正方向逆时针旋转为正，顺时针旋转为负。波线方程：波线方程：外凸的壁面上方，形成膨胀波，即：M1M2M3M4M5随着流线折转，连续膨胀马赫波是发散的。随着流线折转，连续膨胀马赫波是发散的。xy 1 2 3 4超音速气流沿连续上弯壁面的流动右伸膨胀马赫波右伸膨胀马赫波波线方程：波线方程：左伸马赫波波线

29、方程为：跨波线：(特征线特征线)波线方程为：右伸马赫波：跨波线：(特征线特征线)Prandtl-Meyer(普朗陀迈耶普朗陀迈耶)流动流动 如果壁面弯曲段缩成一个点，气流流过如图所示的外凸壁时，可以看作由一系列折转无限小的外凸壁的流动，气流每折转一个角度，就产生一道膨胀波，而气流每经过一道膨胀波，马赫数增大，马赫角减小。因此，这些膨胀波发散，如果壁面的几个折转点都无限接近于点O1，就形成了普朗特迈耶流动。普朗陀迈耶流动的形成M1O2O1O3O12.6 简单波简单波1 1）齐次可约方程组）齐次可约方程组当上面偏微分方程组中的当上面偏微分方程组中的，称方程组为，称方程组为齐次方程组齐次方程组；当齐

30、次方程组的系数都只依赖于当齐次方程组的系数都只依赖于u,v，则方程组称为是，则方程组称为是可约方程组可约方程组。齐次可约方程组齐次可约方程组的相容方程为：的相容方程为：由于右端各项均与由于右端各项均与x,y无关，只与无关，只与u,v有关，上式可以独立积分，有关，上式可以独立积分，从而得到从而得到方程组的特解方程组的特解。2)简单波解的一般形式为了更好地了解简单波的性质，我们需要求简单波的为了更好地了解简单波的性质，我们需要求简单波的一般形式的解：一般形式的解：由一维等熵流的基本方程组：由一维等熵流的基本方程组：对于基本方程组中的第二个方程自动满足，第一个方程化为：的简单波流动，由偏微分方程形式

31、，猜想：因变量可能为两个自变量某组合形式的函数：且且且且发现方程左边为：=0=0猜想正确。表明因变量表明因变量表明因变量表明因变量v v可看做可看做可看做可看做X X的函数。的函数。的函数。的函数。则得通解为：是自变量的任意函数，由问题的的任意函数，由问题的边界条件边界条件决定决定。或？对于的简单波流动，其通解为：或以上解表明，简单波就是向一个方向传播的波，或者说简单波是单向行波。？沿着特征线沿着特征线前面由特征线方法求解曾得到：前面由特征线方法求解曾得到：前面由特征线方法求解曾得到：前面由特征线方法求解曾得到：(C+(C+族特征线）族特征线）族特征线）族特征线）沿着特征线沿着特征线(C-(C

32、-族特征线）族特征线）族特征线）族特征线）由简单波通解得到：由简单波通解得到：由简单波通解得到：由简单波通解得到：当当当当 迹线与迹线与迹线与迹线与C+C+族特征线重合，族特征线重合，族特征线重合，族特征线重合，或，或，或，或，C+C+族特征线是迹线。族特征线是迹线。族特征线是迹线。族特征线是迹线。时，时，时，时，由简单波通解得到：由简单波通解得到：由简单波通解得到：由简单波通解得到：当当当当 迹线与迹线与迹线与迹线与C-C-族特征线重合，族特征线重合，族特征线重合，族特征线重合，或，或，或，或，C-C-族特征线是质点迹线。族特征线是质点迹线。族特征线是质点迹线。族特征线是质点迹线。I）中心简

33、单波：当边界条件使得 时，得到的解为：或中心简单波中心简单波为特征线的共同起为特征线的共同起点，称为中心点。点，称为中心点。II）向前（右传）和向后（左传）简单波：由于简单波时单向行波，按传播方向分为向前和向后简单波。若波的传播速度大于质点速度，则流体质点将从右侧进入波动区，这种简单波称为向前（右传）简单波；反之，称为向后（左传）简单波。的波是向前(右传)简单波；的波是向后(左传)简单波；传播速度为传播速度为III）压缩波和稀疏波：穿过简单波后，若流体的密度和压强增大，称这波为压缩简单波，简称压缩波；若流体的密度和压强减小，则称为稀疏波。xt(0)(I)(II)xt活塞朝气体推进活塞朝气体推进

34、出现间断以前，简单波解成立出现间断以前，简单波解成立活塞背离气体移动活塞背离气体移动活塞背离气体移动活塞背离气体移动例1：设一无限长管道的左半段充满高压气体，右半段是真空，两段间 用一薄膜相隔。高压气体是完全气体，初始状态为:求拆除薄膜后气体的运动。0解解：将气体与真空的交界面初始位置取为坐标原点：将气体与真空的交界面初始位置取为坐标原点x=0=0，当，当t=0=0 时拆除薄膜的一瞬间，气体界面将被加速到某一个速度向右时拆除薄膜的一瞬间，气体界面将被加速到某一个速度向右 运动，同时，向高压气体内传入一个向后稀疏波。运动，同时，向高压气体内传入一个向后稀疏波。则波的传播迹线方程与则波的传播迹线方

35、程与 族特征线重合，其方程为：族特征线重合，其方程为：结合初始条件结合初始条件 ，上式积分为：，上式积分为：沿沿 特征线，满足：特征线，满足：1)2)3)联立两式联立两式得波内的流场分布：得波内的流场分布：界面气体质点的飞散速度？4)5)界面与真空相毗邻的是一个自由面，该自由面上满足：由得：将上式代入中，得界面的飞散速度为：6)7)例例2 2：稀疏波的解：稀疏波的解 由充满气体的管道中抽离活塞产生的管内气体运动，是稀由充满气体的管道中抽离活塞产生的管内气体运动，是稀疏波的一个典型例子。疏波的一个典型例子。t=0=0时刻向左抽动活塞，活塞速度为时刻向左抽动活塞，活塞速度为w(t)=-)=-at,

36、设活塞运动到设活塞运动到A点处时，点处时，速度达到速度达到uA时，活塞不再加速，而是以常速时，活塞不再加速，而是以常速uA继续向左运动继续向左运动xt(0)(I)(II)A 0分析：图中分析：图中(0)(0)区域是尚未受扰动的常态区，区域区域是尚未受扰动的常态区，区域(I)(I)是简单波区，在活塞轨迹是简单波区，在活塞轨迹 上上A A点之后活塞速度为常数，该段的特征线是一族平行线，出现对应的点之后活塞速度为常数，该段的特征线是一族平行线，出现对应的 一个区域一个区域(II)(II)，它与区域，它与区域(I)(I)的分界线是发自的分界线是发自A A点的一条点的一条 特征线。特征线。1)1)区域（

37、区域（I I）内的解：）内的解：对完全气体：对完全气体：沿着沿着x轴，轴，1)根据什么？根据什么？根据什么？根据什么？有：有：如何求如何求?思考：思考：2）联立联立联立联立边界条件？边界条件？根据活塞的边界条件确定任意函数根据活塞的边界条件确定任意函数活塞的轨迹为：活塞的轨迹为：在气体边界上满足：紧靠活塞的气体其位置和速度与活塞相同，在气体边界上满足：紧靠活塞的气体其位置和速度与活塞相同，即：即：3)4)将将 边界条件以及边界条件以及 代入（代入（2 2）式）式5)的函数形式将将代入（代入（2 2）式，得：）式，得：经整理，得：经整理，得：6)7)1）由于在波头 x=c0 t 上气体的速度应v

38、=0，所以取了根号前为 号的根；2）区域I中的坐标x c0t，因此v0，即波内气体都随活塞向同 一方向运动；求解关于求解关于 的二次代数方程，解得：的二次代数方程，解得：8)区域区域(I)(I)中的中的C C+族特征线是一族直线，其方程为：族特征线是一族直线，其方程为：对上式求积分，积分的起点在活塞上对上式求积分，积分的起点在活塞上9)10)由初始条件：由初始条件：代入代入(11)(11)式，于是得到对应于初始时刻的式，于是得到对应于初始时刻的 的方程为：的方程为：11)当活塞在当活塞在 到达到达A A点，速度达到点，速度达到 ,注：指右活塞开始加速运动时为起点发出的第一条特征线；指右活塞开始

39、匀速运动时发出的一条特征线，是II区和I区的交界线。2)2)区域（区域（IIII）内的解：）内的解：区域区域II中的中的 族特征线来自区域族特征线来自区域0，故整个区域，故整个区域II内有内有，12）上式在活塞轨迹线上也成立。上式在活塞轨迹线上也成立。IIII区活塞轨迹线上有区活塞轨迹线上有 ,则得活塞轨迹线上气流声速为：则得活塞轨迹线上气流声速为：13）则沿活塞轨迹线上的黎曼不变量则沿活塞轨迹线上的黎曼不变量 为：为：则整个区域则整个区域IIII中的黎曼不变量中的黎曼不变量 相等，则联立两式可得：相等，则联立两式可得：14）15）从而解出区域从而解出区域IIII中的解：中的解：16）为活塞匀

40、速运动的速度。为活塞匀速运动的速度。17）区域区域IIII的的 族特征线即为区域族特征线即为区域I I中的中的 特征线的平行线特征线的平行线,3)3)逃逸速度：逃逸速度：活塞迹线气体微团迹线亚音速流动区亚音速流动区超音速流动区超音速流动区当活塞由静止连续地向左加速，产生一族向右稀疏波，跨过右传波，满足：在活塞和气体不分离的条件下，与活塞毗邻的气体速度 等于活塞运动速度 ，与活塞毗邻的气体音速18）19）随着活塞向左的速度增加，音速增大还是减小？随着活塞向左的速度增加，音速增大还是减小？族特征线（波线）方程为当即从活塞上会产生一道平行于t轴的驻波20）21）在驻波上当地气流速度(绝对值)等于当地

41、音速，记做：驻波将简单稀疏波区分成两个区域：亚音速流动区和超音速流动区。继续增大活塞速度22）使音速 降到零值，对应的活塞速度为：与之毗邻的气体速度为：23）逃逸速度逃逸速度是气体通过稀疏波膨胀所能达到的最大极限速度，是气体通过稀疏波膨胀所能达到的最大极限速度，之后即使继续加大活塞速度，气体不可能继续加速，因为之后即使继续加大活塞速度，气体不可能继续加速，因为气体与活塞从此分离，它们之间形成真空。实际上逃逸速气体与活塞从此分离，它们之间形成真空。实际上逃逸速度不可能达到。度不可能达到。例例3 3：无限长管道内高压气体推动刚体运动的解：无限长管道内高压气体推动刚体运动的解设在截面积为A管道内x=

42、0截面处有一刚体，其质量为M0，刚体右侧是真空，左侧充满高压气体，其初始状态为：刚体在高压气体的推动下在 时开始向右运动。求刚体松开后的运动轨迹。解：刚体松开后，在高压气体的作用下向右运动，从而有一左传稀疏刚体松开后，在高压气体的作用下向右运动，从而有一左传稀疏 波，使得刚体上的压力随时间发生变化，从而影响刚体的运动速度。波，使得刚体上的压力随时间发生变化，从而影响刚体的运动速度。设刚体运动轨迹为：设刚体运动轨迹为：与刚体毗邻的与刚体毗邻的气体的速度与刚体相同，压气体的速度与刚体相同，压力为作用于刚体上的压力力为作用于刚体上的压力：1）2）由于沿由于沿 族特征线的稀疏波内满足：族特征线的稀疏波

43、内满足：则与刚体毗邻的气体音速为：则与刚体毗邻的气体音速为：3）气体的压力为：求刚体的运动轨迹：4）5）积分式5）：首先令：代入初始条件：5）式化为：得刚体运动速度：对刚体运动速度进行积分，得刚体的运动轨迹：极限速度5 5）压缩波的解）压缩波的解 当活塞朝气体推进时，将导致气体的压力和密度增当活塞朝气体推进时，将导致气体的压力和密度增大，这时产生的简单波就是压缩波。大，这时产生的简单波就是压缩波。若活塞作加速运动，则活塞速度将不断增大，活塞面若活塞作加速运动，则活塞速度将不断增大，活塞面上的气体的速度上的气体的速度 和和 随之增大，于是活塞轨迹随之增大，于是活塞轨迹上发出的上发出的 族特征线式

44、聚拢的，在一定时候将出现同族族特征线式聚拢的，在一定时候将出现同族特征线相交，在相交处出现间断解。特征线相交，在相交处出现间断解。xt活塞朝气体推进活塞朝气体推进出现间断以前，简单波解成立出现间断以前，简单波解成立在出现间断之前，简单波解成立。在出现间断之前，简单波解成立。(I)简单波区的解为：简单波区的解为：1）2）求求的函数形式的函数形式设活塞速度为设活塞速度为 ，则活塞轨迹为：，则活塞轨迹为：3）将活塞轨迹将活塞轨迹及活塞运动时间及活塞运动时间代入 2）式，得得 的函数形式：的函数形式：4）将代入 2)式，得简单波区内的解：5）压缩波在不同时刻的速度分布示意图速度出现多值，不合理形成间断

45、的时间和地点形成间断的时间和地点设时刻 于 处形成间断，这就意味着 时的速度分布曲线成为垂直的，即：在活塞是均匀加速的情况下，间断出现在波头上，波前气体是静止的，即：由则：可求得间断发生的时刻和位置。可求得间断发生的时刻和位置。由得：得：vc12vc21vc21vc21xtx txtxt 四类简单波四类简单波右传膨胀波右传压缩波左传膨胀波左传压缩波 实线为波传播迹线，虚线为气体质点传播迹线n 设管内有流体以速度设管内有流体以速度300 m/s向右运动，设想一小锤向右运动，设想一小锤敲击管壁，于是流体中有疏密相间的音波在传播，传敲击管壁，于是流体中有疏密相间的音波在传播，传播速度为播速度为280

46、 m/s，该弱扰动波传播向上游和下游传播，该弱扰动波传播向上游和下游传播的绝对速度分别为的绝对速度分别为 和和 和和 。n 偏微分方程过点偏微分方程过点(1，1)的特征线方程为的特征线方程为 ，若初始条件为若初始条件为u(0,x)=2x+1,则沿上述特征线的相容方则沿上述特征线的相容方程为：程为：。作业：作业：n已知一流场的两族特征线的黎曼不变量分别为，用表已知一流场的两族特征线的黎曼不变量分别为，用表示的特征线的方程为：示的特征线的方程为：。n充满气体的管道中向充满气体的管道中向右右抽离活塞，气流初始时刻的音抽离活塞，气流初始时刻的音速为速为c0，速度为零。活塞从零时刻开始以速度，速度为零。活塞从零时刻开始以速度 运动，当运动到运动，当运动到 位置时，活位置时，活塞开始以恒速度匀速向右推动。求活塞在加速运动阶塞开始以恒速度匀速向右推动。求活塞在加速运动阶段的速度和音速的解。段的速度和音速的解。n 已知下列线性偏微分方程组，导出特征线方程已知下列线性偏微分方程组，导出特征线方程和相容方程：和相容方程：