流体力学读书笔记

高等流体力学读书笔记

论文题目: 特征线法读书笔记

姓名: 杨志伟

学号: 113108000839

专业:兵器发射理论与技术

指导教师: 周建伟

日期: 2013年12月

1 特征线法

1.1 理论的引出

在考虑了两对面管壁都外折使得两束膨胀波相交，以及膨胀波束在自由边界上反射等问题时，单有绕外钝角流动的公式就不够使用了，需要一种使用于解复杂问题的方法，这就是特征线法。在概况性上说，特征线法确实比绕外钝角的解法进了一步，只要是两个自变数的双曲线型偏微分方程都能用。定常超声速流（包括平面及轴对称的无旋和有旋流）与一维费定常流（不论亚声速还是超声速）的运动方程都是双曲型的。这几种流动能在数学上归在一起，正是反映了在物理上这几种流动都是以波的形式进行变化这样一个事实。

定常亚声速流场上，流动的变化不是以波的形式进行的，任何扰动都没有界线可言，扰动遍及全场，变化都是连续的，任何流动参数（速度、密度和压强等）不仅本身连续，而且它对空间坐标的导数也都连续。与此相反，在定常超声速流场上，扰动都是有界的，像激波在流场中是以突跃面的形式存在的，流动参数本身在突跃面上有突跃的变化，称为强突跃；另一种扰动也是有界的，例如膨胀波（或微弱压缩），界线是马赫波，流动参数本身在波上是连续的，但它的导数在波上可以不连续。如图1所以，定常超声速气流流过外钝角，在第一道膨胀波

O L

11的上游，各流动参数都是均一的，对 的导数到处都是零，但一到

O L上便开始

11

变化了，虽然流动参数本身在

O L还是连续的，但无变化的直匀流区突然在这条

11

线上有变化的扇形膨胀地带相接，诸流动参数的导数在

O L上必是突然从零变为

11

某一定值。在最后一道波

O L上，导数从一定值跃变为零。在中间各道波上，流

12

动参数的导数取了特殊的突跃值——零。

图1 外钝角绕流

流动参数的导数在其上发生突跃的线（或面）称为弱突跃线（或面），以区别与流动参数本身发生突跃的强突跃面（如激波）。马赫线（膨胀波或微弱压缩波）就是弱突跃线。当然，突跃之中可以包括突跃值为零（即不发生突跃）的情况。这样就可以把绕外钝角流动的扇形地带中的每一道膨胀波都包括在内，都是弱突跃线。

上述几种流动的运动方程都是双曲型二阶偏微分方程。解双曲型二阶偏微分方程时，存在有一类特殊的曲线，即函数的二阶导数在这种曲线上所以不连续，因而知道了因数在这种曲线一侧的数值时，不可能靠级数展开去推得函数在曲线另一侧的数值来。双曲型二阶偏微分方程所独有的这种函数导数可以不连续的曲线，在数学上称为特征线。流场上的弱突跃线就是数学上的特征线。根据数学上的特征线理论，用数值解法解决定常超声速流动的问题，就是气体动力学中的所谓特征线法。

特征线理论除了定义特征线、并根据定义找出特征线方程之外，还确定；沿特征线上的各参数的变化之间必须遵守一定的关系，这是特征线所具有的另一方面的性质。所以说，特征线有两种性质：一是跨过这种曲线，函数(例如位函数)的二阶导数可以不连续；二是沿着这种曲线，函数的一阶导数(例如流速)的变化之间又必须遵守一定的规律。就是说，既有突跃的一面、又有规律性的一面。在做计算时，要同时利用这两种性质，尤其是利用沿特征线的规律性。这个规律性在绕外钝角流动的问题中已经用过了，那就是气流折角与流速有一一对应的关系。不过，这个关系在绕外钝角流动的问题里，没有从数学上强调它是沿特征线的关

系，表面上看来反而像是跨特征线似的。原因是，在那个简单问题上，没有全面讨论流场，只是画了要用的马赫线。事实上，与流线成μ角的马赫线还有另一族存在，如图2上的PQR 线。后面可以证明，在这个具体问题上，所谓沿特征线的变化关系正是沿PQR 那一族马赫线的变化关系。

图2 外钝角绕流的马赫线

本节先把定常平面超声速流和轴对称超声速流的运动方程列一下，以说明数学上的共同点，然后导出这些方程列的特征线来。接下去详细讲解特征线法在乎面无旋定常超声速流问题上的应用及解题法。后面再讲超声流流过圆锥体的稻确解及轴对称特征线法。

1.2 两个自变数的运动方程

平面不可压位流的运动方程是

22220x y

φφ

∂∂+=∂∂ 而在平面定常无旋可压流的条件下，此式府写为如下形式

2

2222222221111120c x x c y y

c x y x y φφφφφφφ

⎡⎤⎡⎤⎛⎫⎛⎫∂∂∂∂∂∂∂⎛⎫-∙+-∙-∙∙∙

=⎢⎥⎢⎥ ⎪ ⎪ ⎪∂∂∂∂∂∂∂∂⎝⎭⎢⎥⎝⎭⎝⎭⎢⎥⎣⎦⎣⎦ 通乘以2c ，并引用符号

x x φφ∂=∂，y y φ

φ∂=∂，22xx x φφ∂=∂，22yy y φφ∂=∂，2xy x y

φφ∂=∂∂

则上式化为

()()2

22220x xx x y xy y yy c

c φφφφφφφ-∙-∙∙+-∙= （1-1）

定常无旋轴对称流的运动方程直接写为

()()2

2

22

220x

xx x r xr r

rr r c c c r φφφφφφφφ-∙-∙∙+-∙+= （1-2）

这两个方程是同一类型的，都是具有两个独立变数的二阶非线性偏微分方程。它们有时也称为拟线性方程，因为对最高阶导数而言是线性的。此二式可写成一个共同的形式

20xx xy yy A B C D φφφ+++= （1-3）

一般说来，式中的系数A ，B ，C ，D 是x ，y ，x φ，y φ的函数。这样一个二阶偏微分方程，究竟属于哪种类型，要看判别式()2B AC -大于、等于、还是小于零。现在式（1-1）中的

()()()()22222222222

222x y x y x y B AC c c c c c V c φφφφφφ⎡⎤-=∙--∙-=+-=-⎣⎦

式（1-2）中的

()()22222

222x y B AC c c c V c φφ⎡⎤-=+-=-⎣⎦

在超声速流(V c >)中，这两个判别式都大于零，所以式（1-1）与（1-2）都是双曲型的方程。现在总概括成式（1-3）来研究，目的是要得到解决实际问题的具体办法。不过，并不能单纯从数学角度来研究式（1-3）的解法，而是要与气体流动的物理情况密切结合起来。函数φ就是速度位，x φ就是x 分速，y φ就是y 分速。也常把x φ和y φ组成的平面称为“速度面”。

1.3 哥西问题及特征线方程

在一定的边界条件下直接积分式（1-3）是很困难的，所以就希望用数值解的办法求出需要的答案。以平面流动来说，假设知道了待求函数在某一条曲线上的函数值，如果能够设法求出函数在该曲线附近的数值，就能把全流场上的解答一步步地找出来，这在数学上叫做哥西问题。哥西问题的数学提法是这样：给定函数的偏微分方程，并在xy 平面上沿某一条曲线L 给定该函数的数值，问能否推出L 曲线附近的函数值。