资源描述：
《特征线理论及应用课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、第二章特征线理论及应用气体动力学中，有大量问题是用双曲型偏微分方程来描述的，很难得到解析结果，在这种情况下，有两种数值解法：1）特征线数值解法：求解域用特征线网格进行离散，求各网格结点上的解；气体动力学中，有大量流动问题是用双曲型偏微分方程来描述的，宜于用特征线方法求解。2）有限差分法：求解域的有限差分网格一般是正交的，根据由偏微分方程构造的差分格式来求各网格结点上的解。§2.1特征线理论特征线的数学定义考虑一个一般的一阶双曲型偏微分方程：x,y是两个自变量，u(x,y)是因变量。系数A1、A2及非齐次项F1可以是x，y，u的函数。（1）将偏微分方程改写为：设未知函数u(x,y)
2、连续，u的一阶导数可以写作：【注:u的一阶导数可以不连续】偏微分方程的特征线定义为：xy平面内具有斜率为的曲线。（2）（3）沿着特征线或：偏微分方程可化简为：代入式（4）得到偏微分方程的相容方程【是平面上这样一族曲线：沿着此族中任一曲线(a)，可以把待求物理量的一阶偏微分控制方程变换成等价的常微分控制方程(b)，称为原偏微分方程或偏微分方程组的相容方程】特征线的第一个数学意义：(a)(b)特征线的第二个数学意义：上两式表明：沿着特征线，分母和分子均为零。即沿着特征线，表明：1）沿特征线因变量的一阶导数具有不定值，可以是不连续的，在这种情况下，特征线是弱间断（第一类间断线）。2）在
3、气体动力学中，特征线可以是弱扰动波传播的迹线，或者说弱扰动传播的迹线就是特征线。因此，因变量的一阶导数只允许有弱间断，如果在物理平面上有激波出现，在强间断面上便无法建立因变量的全微分式，也就不能用特征线方法求解。例：一阶偏微分方程的初始条件是2）沿此特征线的相容方程3）u(2,4)的值用特征线法确定：1）通过点(2,4)的特征线解：（1）对照一般形式的双曲型偏微分方程该方程对应的系数：A1=1,A2=2x,F1=－3x2则特征线方程为：积分得：为确定过点(2,4)的特征线，将x=2,y=4，代入上式得：所以，所求的特征线方程是：对上式积分，得：（2）偏微分方程的相容方程为：如何确
4、定C2？初始条件u(0,y)=5y+10及特征线方程u(0,0)=10因此相容方程为：§2.2一维等熵流动的特征线数值解法基本方程与黎曼不变量（连续方程）（动量方程）（以一维等直截面管为例）基本方程等熵流动中只有一个状态参量独立：将基本方程中的用代替，得：基本方程可化为：两式相加减合并，基本方程可写作：定义则基本方程化为以vG为新的未知函数的偏微分方程：基本方程——偏微分方程特征线？？相容方程？？在x-t平面上，把dx/dt=vc曲线称为偏微分方程的特征线。C＋C－xtC＋表示第一族特征线；C－表示第二族特征线。解相容方程:对多方气体：其相容方程的解为：由声速：沿着特征线沿着
5、特征线结论：特征线的基本性质1）一维非定常流动中，平面x-t上任一点，都有两条不同族的特征线，沿各特征线有各自不同的黎曼不变量；2）特征线上参量v,c,p,…的一阶导数可以不连续，但这些参量本身是连续的，称因变量的一阶导数不连续的点叫做弱间断。如果初始某一点有弱间断，那么这个弱间断必定会沿着过该点的特征线向外传播。3）两个相邻的，不同类型流动区域的分界线，必定是特征线。三类流态中的特征线vcc0v0定常均匀流动相容关系描述的状态特征线xt特征线（不代表波的传播迹线）xt(0)(I)(II)v/c0c/c0(0)(1)(2)(3)(4)简单波流动特征线相容关系描述的状态特征线活塞
6、运动迹线复合波流动特征线相容关系描述的状态特征线xtc/c1v/c1C+C-76523489102345678910依赖区和影响区由于沿着两族特征线，分别有：可以把J＋和J－看作是两个新的函数，则利用J＋和J－表示的特征线方程为：第I族特征线斜率仅由J-决定；第II族特征线斜率仅由J+决定。xtDABD点的依赖区MC＋C－C＋C－在平面运动中，沿着特征线黎曼不变量保持不变，这一重要性质清楚地揭示出流体动力学中的一些依赖关系。设t=0时各量沿x轴的分布为v0(x),c0(x),于是可知黎曼不变量的相应分布为。则(x,t)平面上任意一点D(x,t)上的状态，将直接由x轴上点A(xA,
7、0),B(xB,0)两点上的状态决定。所以，点所处的状态将完全由且只由线段AB上的值决定，线段AB就称为点D的依赖区。DMM点的影响区xtABC＋C－同样，能够受到AB线段间某点M的初始值影响的区域，是由发自M点的与发自M点的所包围的区域，而这个区域之外的地方，都不受M点的影响。这个区域称为M点的影响区。PQ例：已知初始时刻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)点，有根据：得：由§
