工程流体力学第四章不可压缩流体的有旋流动和二维无旋流动

1、第四章不可压缩流体的有旋流 动和二维无旋流动 第一节 流体微团运动分析 第二节 有旋流动和无旋流动 第三节 无旋流动的速度势函数 第四节 二维平面流动的流函数 第五节 基本的平面有势流动 第六节 平面势流的叠加流动 欢迎进入第四章的学习 流体由于具有易变形的特性（易流动性），因此流体 的运动要比工程力学中的刚体的运动复杂得多。在流体运 动中，有旋流动和无旋流动是流体运动的两种类型。由流 体微团运动分析可知，有旋流动是指流体微团旋转角速度 的流动，无旋流动是指 的流动。 实际上，黏性流体的流动大多数是有旋流动，而且有 时是以明显的旋涡形式出现的，如桥墩背流面的旋涡区， 船只运动时船尾后形成的旋涡

2、，大气中形成的龙卷风等等。 但在更多的情况下，流体运动的有旋性并不是一眼就能看 得出来的，如当流体绕流物体时，在物体表面附近形成的 速度梯度很大的薄层内，每一点都有旋涡，而这些旋涡肉 眼却是观察不到的。至于工程中大量存在着的紊流运动， 更是充满着尺度不同的大小旋涡。 0 0 流体的无旋流动虽然在工程上出现得较少，但 无旋流动比有旋流动在数学处理上简单 得多，因 此，对二维平面势流在理论研究方面较成熟。对 工程中的某些问题，在特定条件下对黏性较小的 流体运动进行无旋处理，用势流理论去研究其运 动规律，特别是绕流物体的流动规律，对工程实 践具有指导意义和应用价值。因此，本章先阐述 有旋流动的基本概

3、念及基本性质，然后再介绍二 维平面势流理论。 第一节 流体微团运动分析 刚体的一般运动可以分解为移动和转动 两部分。流体与刚体的主要不同在于它具 有流 动性，极易变形。因此，任一流体微 团在运动过程中不但与刚体一样可以移动 和转动，而且还会发生变形运动。所以， 在一般情况下流体微团的运动可以分解为 移动、转动和变形运动三部分。 一、表示流体微团运动特征的速度表达式一、表示流体微团运动特征的速度表达式 z z u y y u x x u uu c ddd z z v y y v x x v vvcddd z z w y y w x x w wwcddd 图 4-1 分析流体微团运动用图 y y

4、u x v z x w z u z x w z u y x v y u x x u uucd 2 1 d 2 1 d 2 1 d 2 1 d z z v y w x y u x v z y w z v x y u x v y v vvcd 2 1 d 2 1 d 2 1 d 2 1 dy x x w z u z z v y w y z v y w x z u x w z z w wwcd 2 1 d 2 1 d 2 1 d 2 1 d 剪切变形速率 、 、 、 、 、 ， 引入记号，并赋予运动特征名称： 线变形速率 、 、 ， xx 、 yy 、 zz ， z w y v x u zzyyxx

5、 , xy yx yz zy xz zx x w z u z v y w y u x v xzzx zyyz yxxy 2 1 2 1 2 1 （4-1） （4-2） 于是可得到表示流体微团运动特征的速度表达式为 旋转角速度 、 、 ， x y z y u x v x w z u z v y w z y x 2 1 2 1 2 1 （4-3） xyyxzww zxzxyvv yzzyxuu yxzyzxzzc xzyzyxyyc zyxzxyxxc ddddd ddddd ddddd （4-4） 二、流体微团运动的分解二、流体微团运动的分解 为进一步分析流体微团的分解运动及其几何特 征，对式(

6、4-4)有较深刻的理解，现在分别说明流 体微团在运动过程中所呈现出的平移运动、线变 形运动、角变形运动和旋转运动。 为简化分析，仅讨论在 平面上流体微团的运 动。假设在时刻 ，流体微团abcd为矩形，其上 各点的速度分量如图4-2所示。由于微团上各点的 速度不同，经过时间 ，势必发生不同的运动， 微团的位置和形状都将发生变化，现分析如下。 xoy t t d 1平移运动 图 4-2 分析流体微团平面运动用图 a 2线变形运动 b 图4-3 流体微团平面运动的分解(a) 图4-3 流体微团平面运动的分解(b) 图4-3 流体微团平面运动的分解(c) 图4-3 流体微团平面运动的分解(d) 3角变

7、形运动 c y u x v t yxxy 2 1 d )/2dd( y u x v yxxy 2 1 z v y w zyyz 2 1 x w z u xzzx 2 1 4旋转运动 d y u x v t z 2 1 d d-d 2 1 tdd/ )-d( 2 1 y u x v x w z u z v y w z y x 2 1 2 1 2 1 222 zyx )( 2 1 vkji zyx x y 综上所述，在一般情况下，流体微团的运动总是可 以分解成：整体平移运动、旋转运动、线变形运动及角变 形运动，与此相对应的是平移速度、旋转角速度、线变形 速率和剪切变形速率。 第二节 有旋流动和无旋

8、流动 一、有旋流动和无旋流动的定义一、有旋流动和无旋流动的定义 二、速度环量和旋涡强度二、速度环量和旋涡强度 一、有旋流动和无旋流动的定义一、有旋流动和无旋流动的定义 流体的流动是有旋还是无旋，是由流体微团本身是否旋转 来决定的。流体在流动中，如果流场中有若干处流体微团具 有绕通过其自身轴线的旋转运动，则称为有旋流动。如果在 整个流场中各处的流体微团均不绕自身轴线的旋转运动，则 称为无旋流动。这里需要说明的是，判断流体流动是有旋流 动还是无旋流动，仅仅由流体微团本身是否绕自身轴线的旋 转运动来决定，而与流体微团的运动轨迹无关，在图4-4(a)中， 虽然流体微团运动轨迹是圆形，但由于微团本身不旋

9、转，故 它是无旋流动；在图4-4(b)中，虽然流体微团运动轨迹是直线， 但微团绕自身轴线旋转，故它是有旋流动。在日常生活中也 有类似的例子，例如儿童玩的活动转椅，当转轮绕水平轴旋 转时，每个儿童坐的椅子都绕水平轴作圆周运动，但是每个 儿童始终是头向上，脸朝着一个方向，即儿童对地来说没有 旋转。 图4-4 流体微团运动 无旋流动 有旋流动 判断流体微团无旋流动的条件是：流体中每一个流体微团都满足 根据式（4-3），则有 0 zyx , z v y w , x w z u y u x v （4-8） 二、速度环量和旋涡强度二、速度环量和旋涡强度 1速度环量 为了进一步了解流场的运动性质，引入流体力

10、 学中重要的基本概念之一速度环量。 在流场中任取封闭曲线k，如图4-5所示。速度 沿该封闭曲线的线积分称为速度沿封闭曲线k的环 量，简称速度环量，用 表示，即 式中 在封闭曲线上的速度矢量； 速度与该点上切线之间的夹角。 速度环量是个标量，但具有正负号。 v kk svsvdcosd v 图4-5 沿封闭曲线的速度环量 在封闭曲线k 上的速度矢 量 速度 与该点 上切线之间 的夹角 v 速度环量的正负不仅与速度方向有关，而且与积分时 所取的绕行方向有关。通常规定逆时针方向为k的正方向， 即封闭曲线所包围的面积总在前进方向的左侧，如图4-5 所示。当沿顺时针方向绕行时，式（4-9）应加一负号。

11、实际上，速度环量所表征的是流体质点沿封闭曲线k运动 的总的趋势的大小，或者说所反映的是流体的有旋性。 由于和，则 kwj vi uv kzj yi xs dddd zwyvxusvdddd 代入式（4-9），得 kk zwyvxusv)ddd(d （4-10） 2旋涡强度 沿封闭曲线的速度环量与有旋流动之间有一个重要 的关系，现仅以平面流动为例找出这个关系。如图4- 6所示，在平面上取一微元矩形封闭曲线，其面 积，流体在a点的速度分量为和，则b、 c和d点的速度分量分别为： xoy yxaddd uv x x u uud b x x v vvd b y y u x x u uudd c y y

12、 v x x v vvdd c y y u uud d y y v vvd d 图4-6 沿微元矩形的速度环量 x x u ud x x v vd y y u x x u udd y y v x x v vdd y y u ud y y v vd 于是，沿封闭曲线反时针方向abcda的速度环量 将 、 、 、 和 、 、 、 各值代入上式，略去高于一阶 的无穷小各项，再将式（4-3）的第三式代入后，得 然后将式（4-11）对面积积分，得 y vv x uu y vv x uu d 2 d 2 d 2 d 2 d addccbba a u b u c u d u a v b v c v d v

13、ayx y u x v z d2ddd （4-11） a z d2 （4-12） 于是得到速度环量与旋转角速度之间关系的斯托克斯定理： 沿封闭曲线的速度环量等于该封闭周线内所有的旋转角速 度的面积积分的二倍，称之为旋涡强度i，即 和 式中 在微元面积 的外法线 上的分量。 ai nd 2d ai n d2 （4-13） n adn 由式（4-11）可导出另一个表示有旋流动的量， 称为涡量，以 表示之。它定义为单位面积上的速 度环量，是一个矢量。它在z轴方向的分量为 对于流体的空间流动，同样可求得x和y轴方向 涡量的分量 和 。于是得 即 zz y u x v a 2 d d zz yy xx

14、y u x v x w z u z v y w 2 2 2 v 2 （4-14） （4-15） 也就是说，在有旋流动中，流体运动速度 的旋度称为 涡量。 由此可见，在流体流动中，如果涡量的三个分量中有一 个不等于零，即为有旋流动。如果在一个流动区域内各处 的涡量或它的分量都等于零，也就是沿任何封闭曲线的速 度环量都等于零，则在这个区域内的流动一定是无旋流动。 下面举两个简单的例子来说明速度环量和旋涡强度的物 理意义，以及有旋流动和无旋流动的区别。 v 【例例4-1】 一个以角速度 按反时针方向作像刚体 一样的旋转的流动，如图4-7所示。试求在这个流 场中沿封闭曲线的速度环量，并证明它是有旋流

15、动 . (解) 【例例4-2】 一个流体绕o点作同心圆的平面流动， 流场中各点的圆周速度的大小与该 点半径成反 比，即 ，其中c为常数，如图4-8所示。试 求在流场中沿封闭曲线的速度环量，并分析它的 流动情况。(解) rcv 【解解】 在流场中对应于任意两个半径 和 的圆周 速度各为 和 ，沿图中画斜线扇形部分 的周界abcda的速度环量 可见，在这个区域内是有旋流动。又由于扇形面积 于是 上式正是斯托克斯定理的一个例证。 以上结论可推广适用于圆内任意区域内。 1 r 2 r 11 rv 22 rv )()( 2 1 2 211221122dacdbcababcda rrrvrvrvrv )(

16、 2 d 2 1 2 2 2 1 rrrra r r a2 abcda 图4-7 有旋流动中速度环量的计算图4-8 无旋流动中速度环量的计算 【解解】 沿扇形面积周界的速度环量 可见，在这区域内是无旋流动。这结论可推广适用于任何 不包围圆心o的区域内，例如 。若包有圆心( )， 该处速度等于无限大，应作例外来处理。现在求沿半径 的圆周封闭曲线的速度环量 上式说明，绕任何一个圆周的流场中，速度环量都不 等于零，并保持一个常数，所以是有 旋流动。但凡是绕 不包括圆心在内的任何圆周的速度环量必等于零，故在圆 心o点处必有旋涡存在，圆心是一个孤立涡点，称为奇点。 0 1 1 2 2 dacdbcaba

17、bcda r r c r r c adcba 0r 2 0 2d常数cr r c 第三节 无旋流动的速度势函数 如前所述，在流场中流体微团的旋转角速度 在任意 时刻处处为零，即满足 的流动为无旋流动，无旋流 动也称为有势流动。 一、速度势函数引入一、速度势函数引入 二二、速度势函数的性质、速度势函数的性质 0 v 一、速度势函数引入一、速度势函数引入 由数学分析可知， 是 成为某一标量 函数 全微分的充分必要条件。则函数 称为速度 势函数。因此，也可以说，存在速度势函数 的流动为有 势流动，简称势流。根据全微分理论，势函数 的全微分 可写成 于是得 0 v zwyvxuddd )(tzyx，

18、z z y y x x dddd z w y v x u ，（4-16） 按矢量分析 对于圆柱坐标系，则有 于是 从以上分析可知，不论是可压缩流体还是不可 压缩流体，也不论是定常流动还是非定常流动， 只要满足无旋流动条件，必然存在速度势函数。 grad k z j y i x kwj viuv z v r v r v zr ， 1 zvrvrv zr dddd （4-17） （4-18） 二、速度势函数的性质二、速度势函数的性质 （1）不可压缩流体的有势流动中，势函数 满足拉普拉斯 方程，势函数 是调和函数。 将式（4-16）代入到不可压缩流体的连续性方程（3-28） 中，则有 式中 为拉普拉

19、斯算子，式（4-19）称为拉普拉 斯方程，所以在不可压流体的有势流动中，速度势必定满 足拉普拉斯方程，而凡是满足拉普拉斯方程的函数，在数 学分析中称为调和函数，所以速度势函数是一个调和函数。 0 2 2 2 2 2 2 2 zyx 2 2 2 2 2 2 2 zyx 0 z w y v x u 从上可见，在不可压流体的有势流动中，拉普拉斯 方程实质是连续方程的一种特殊形 式，这样把求解无 旋流动的问题，就变为求解满足一定边界条件下的拉普 拉斯方程的问题。 （2）任意曲线上的速度环量等于曲线两端点上速度势 函数 值之差。而与曲线的形状无关。 根据速度环量的定义，沿任意曲线ab的线积分 这样，将求

20、环量问题，变为求速度势函数值之差的问题。 对于任意封闭曲线，若a点和b点重合，速度势函数是单 值且连续的，则流场中沿任一条封闭曲线的速度环量等于 零，即 。 ab b a b a b a ab dwdzvdyudxsdv )( 0ab 第四节 二维平面流动的流函数 一、流函数的引入一、流函数的引入 对于流体的平面流动，其流线的微分方程为 ,将其改写成下列形式 （4-20） 在不可压缩流体的平面流动中，速度场必须满足不可压缩流体的连续性方程， 即 或 （4-21） 由数学分析可知，式（4-21）是（ ）成为某函数全微分的充分必要 条件，以 表示该函数，则有 （4-22） 函数称为流场的流函数。由

21、式（4-22）可得 （4-23） vyuxdd 0ddyuxv 0 y v x u y v x u yuxvdd yuxvy y x x ddddd ），（yx x v y u ， 由式(4-22)，令 ，即 常数，可得流线微分方程式(4-20)。 由此可见, 常数的曲线即为流线，若给定一组常数值，就可得 到流线簇。或者说，只要给定流场中某一固定点的坐标（ ）代入 流函数 ，便可得到一条过该点的确定的流线。因此，借助流函数可 以形象地描述不可压缩平面流场。 对于极坐标系，可写成 （4-24） （4-25） 在已知速度分布的情况下，流函数的求法与速度势函数一样，可由曲 线积分得出。 至此可看到，

22、在不可压缩平面流动中，只要求出了流函数 ， 由式（4-23）或式（4-24）就可求出速度分布。反之，只要流动满足 不可压缩流体的连续性方程，不论流场是否有旋，流动是否定常，流 体是理想流体还是黏性流体，必然存在流函数 。 这里需说明，等流函数线与流线等同，仅在平面流动时成立。对 于三维流动，不存在流函数，也就不存在等流函数线，但流线还是存 在的。 0d ），（yx 00 yx ， r v r 1 r v dddrvrv r ），（yx 二、流函数的性质二、流函数的性质 （1）对于不可压缩流体的平面流动，流函数 永远 满足 连续性方程。 将式（4-23）代入式（4-21）得 即流函数永远满足连续

23、性方程。 （2）对于不可压缩流体的平面势流，流函数 满足拉普 拉斯方程，流函数也是调和函数。 对于平面无旋流动， ，则 将式（4-23）代入上式 因此，不可压缩流体平面无旋流动的流函数也满足拉普拉 斯方程，也是一个调和函数。 因此，在平面不可压缩流体的有势流场中的求解问题，可 以转化为求解一个满足边界条件的 的拉普拉斯方程. yxxy 22 0 z 0 y u x v 0 2 2 2 2 2 yx （3）平面流动中，通过两条流线间任一 曲线单位厚度的体积流量等于两条流线的 流函数之差。这就是流函数 的物理意义。 如图4-9所示，在两流线间任一曲线 ab，则通过单位厚度的体积流量为 （4-26）

24、 由式（4-26）可知，平面流动中两条流线 间通过的流量等于这两条流线上的流函数 之差。 图4-9 说明流函数物理意义用图 2 1 2 1 2 1 2 1 dd)d(d x x y y x x y y v x x y y xvyuq 12 , ),( 22 11 d yx yx 三、三、 和和 的关系的关系 （1）满足柯西-黎曼条件 如果是不可压缩流体的平面无旋流动，必然同时存在着速 度势和流函数，比较式（4-16）和式（4-23），可得到速度 势函数和流函数之间存在的如下关系 （4-27） （4-28） 这是一对非常重要的关系式，在高等数学中称作柯西-黎曼 条件。因此， 和 互为共轭调和函数

25、，这就有可能使我们利 用复变函数这样一种有力的工具求解此类问题。 当势函数 和 流函数二者知其一时，另一个则可利用 式（4-27）的关系求出，而至多相差一任意常数。 xyyx ， 0 yyxx （2）流线与等势线正交。 式（4-28）是等势线簇 常数和流线簇 常数互相正交的条 件，若在同一流场中绘出相应的一 系列流线和等势线，则它们必然构 成正交网格，称为流网，如图4-10 所示。 ），（yx ），（yx 图4-10 流网 0 yyxx 【例例4-3】 有一不可压流体平面流动的速度分布 为 。该平面流动是否存在流函数和速度 势函数；若存在，试求出其表达式；若在流场中a （1m，1m）处的绝对压

26、强为1.4105pa，流体的密度 1.2kg/m3，则b（2m，5m）处的绝对压强是多少？ 【解解】 （1）由不可压流体平面流动的连续性方程 该流动满足连续性方程，流动是存在的，存在流函数。 由于是平面流动 该流动无旋，存在速度势函数。 yvxu44， 0)4()4( y y x xy v x u 0 yx 0 44 2 1 2 1 y x x y y u x v z （2）由流函数的全微分得： 积分 由速度势函数的全微分得： 积分 （3）由于 ，因此，a和b处的速度分别为 由伯努里方程 可得 yxxyyuxvy y x x d4d4ddddd cxy 4 yyxxyvxuy y x x d4

27、d4ddddd cyx)( 2 22 222 vuv )(32) 14() 14( 22222 a smv )(464) 54() 24( 22222 b smv 2 bb 2 aa 2 1 2 1 vpvp )(8 .139740)46432(2 . 1 2 1 104 . 1)( 2 1 52 b 2 aab pavvpp 第五节 基本的平面有势流动 流体的平面有势流动是相当复杂的，很多复 杂的平面有势流动可以由一些简单的有势 流动叠 加而成。所以，我们首先介绍几种基本的平面有 势流动，它包括均匀直线流动，点源和点汇、点 涡等 一、均匀直线流动一、均匀直线流动 流体作均匀直线流动时，流场中

28、各点速度的大小相等，方 向相同，即 和 。由 式（4-16）和式（4-23），得 于是速度势和流函数各为 以上两式中的积分常数 和 可以任意选取，而不影响流体 的流动图形（称为流谱）。 0 uu 0 vv 00, v xy vu yx u 10000 ddddcyvxuyvxuy y x x 20000 d)d(ddcyuxvyuxvy y x x 1 c 2 c 若令 ，即得均匀直线流动的速度势和流函数各为 （4-29） （4-30） 由式（4-29）和式（4-30）可知，等势线簇（ 常 数）和流线簇（ =常数）互相垂直，如图4-11 所示。各流线与轴的夹 角等于 。 由于流场中各点的速度都

29、相等，根据伯努里方程（3-41）， 得 常数 如果均匀直线流动在水平面上，或流体为气体，一般可以忽 略重力的影响，于是 常数 即流场中压强处处相等。 0 21 cc yvxu 00 yuxv 00 yvxu 00 yuxv 00 0 01 - tg u v g p z p 图4-11 均匀直线流的流谱 二、平面点源和点汇二、平面点源和点汇 如果在无限平面上流体不断从一点沿径向直线均匀地向 各方流出，则这种流动称为点源，这个点称为源点（图4-12， a）；若流体不断沿径向直线均匀地从各方流入一点，则这 种流动称为点汇，这个点称为汇点（图4-12，b）。显然，这两 种流动的流线都是从原点 o发出的

30、放射线，即从源点流出和 向汇点流入都只有径向速度 。现将极坐标的原点作为源点 或汇点，则 r v r r v 0 v rv r dd 图4-12 点源和点汇的流谱 点源 点汇 back 根据流动的连续性条件，流体每秒通过任一半径为 的单位长度圆柱 面上的流量 都应该相等，即 常数 由此得 （4-31） 式中 是点源或点汇在每秒内流出或流入的流量，称为点源强 度或点汇强度。对于点源， 与 同向， 取正号；对于点汇， 与 异向， 取负号，于是 积分得 式中积分常数 是任意给定的，现令 。又由于 ，于是得 速度势 （4-32） 当 时，速度势 和 速度都变成无穷大，源点和汇点都是奇点。所 以速度势

31、和速度 的表达式（4-31）和式（4-32）只有在源点和汇点 以外才能应用。 r v q vr qrv12 r q v v r 2 v q r vr v rrv q v q r rqvd 2 d cr qv ln 2 c0c 22 yxr 22 ln 2 ln 2 yx q r q vv 0r r v r v 现在求流函数，由式（4-25） 积分得(令式中的积分常数为零) （4-33） 等势线簇（ 常数，即 常数）是同心圆簇（在图4-12中 用虚线表示）与流线簇（ 常数，即 常数）成正交。而 且除源点或汇点外，整个平面上都是有势流动。 如果 平面是无限水平面，则根据伯努里方程（341） 式中

32、为 在处的流体压强，该处的速度为零。 将式（4-31）代入上式，得 （4-34） 由式（4-34）可知，压强 随着半径 的减小而降低。当 时， 。图4-13表示当 时，点汇沿 半径 的压强分布。 d 2 dddd v rr q rvrvrv x yqq vv1- tg 22 r xoy g p g v g p r 2 2 pr 22 2 1 8r q pp v p r 2/ 1 22 0 )8/( pqrr v rr0 0p r 图4-13 点汇沿半径的压强分布 三、点涡三、点涡 设有一旋涡强度为 的无限长直线涡束，该涡束以等角 速度 绕自身轴旋转，并带动涡束周围的流体绕其环流。由 于直线涡束

33、为无限长，所以可以认为与涡束垂直的所有平面 上的流动情况都一样。也就是说，这种绕无限长直线涡束的 流动可以作为平面流动来处理。由涡束所诱导出的环流的流 线是许多同心圆，如图4-14所示。根据斯托克斯定理可知， 沿任一同心圆周流线的速度环量等于涡束的旋涡强度，即 常数 于是 （4-35） 因此涡束外的速度与半径成反比。若涡束的半径 ，则 成为一条涡线，这样的流动称为点涡，又称为纯环流。但当 时， ，所以涡点是一个奇点。 i irv 2 0 2 r v r v， 0 0 r 0 0 r v 图4-14 点涡的流谱 现在求点涡的速度势和流函数。由于 由 积分后得速度势 （4-36） 又由于 由 积分

34、后得流函数 （4-37） 当 时，环流为反时针方向，如图4-14所示；当 时，环流 为顺时针方向。 由式（4-36）和式（4-37）可知，点涡的等势线簇是经过涡点的放射 线，而流线簇是同心圆。而且除涡点外，整个平面上都是有势流动。 rr v r vr 2 1 0 ， d 2 d 1 dd r r r r x y 1- tg 22 rr v r vr 2 0 1 ， r r r r r r d 2 d 1 dd rln 2 0 0 设涡束的半径为 ，涡束边缘上的速度为 ，压强 为 ； 时的速度显然为零，而压强为 。代入伯努里 方程（3-41），得涡束外区域内的压强 分布为 （4-38） 由式（4

35、-38）可知，在涡束外区域内的压强随着半径的减小 而降低，涡束外缘上的压强为 或 （4-39） 所以涡束外区域内从涡束边缘到无穷远处的压强降是一个常 数。又由式（4-38）可知，在 处，压强 ，显然这 是不可能的。所以在涡束内确实存在如同刚体一样以等角速 度旋转的旋涡区域，称为涡核区。由式（4-39）可得涡核的 半径 0 r 0 0 2 r v 0 p r p 22 22 1 82r p v pp 2 0 2 22 0 0 1 82r p v pp 2 0 2 2 2 00 1 82 1 r vpp 0r p 常数 g v g p z 2 2 由于涡核内是有旋流动，故流体的压强可以根据欧拉运动

36、微 分方程求得。平面定常流动的欧拉运动微分方程为 将涡核内任一点的速度 和 代入上两式，得 以 和 分别乘以上两式，然后相加，得 或 积分得 x p y u v x u u 1 y p y v v x v u 1 yuxv x p x 1 2 y p y 1 2 xdyd y y p x x p yyxxdd 1 )dd( 2 pyxd 1 d 2 22 2 ）（ cvcrcyxp 222222 2 1 2 1 2 1 ）（ 在 处， ，代入上式，得 最后得涡核区域内的压强分布为 （4-40） 或 （4-40a） 于是涡核中心的压强 而涡核边缘的压强 所以 可见，涡核内、外的压强降相等，都等于

37、用涡核边缘速度计 算的动压头。涡核内、外的速度分布和压强分布如图4-15所 示。 0 rr 00 vvpp 、 2 0 2 0 2 00 2 00 2 1 2 1 vpvvpvpc 22 0 2 1 vvpp 222 0 2 2 1 rrpp 2 0 22 0 rpvppc 2 0 22 00 2 1 2 1 rpvpp 2 0 22 000 2 1 2 1 2 1 rvpppppp cc ）（ 图5-14 涡流中涡核内、外的速度和压强分布 第六节 平面势流的叠加流动 从上节可以看到，只有对一些简单的有势流动， 才能求出它们流函数和势函数，但当流动较复杂时， 根据流动直接求解流函数和势函数往往

38、十分困难。 我们可以将一些简单有势流动进行叠加，得到较复 杂的流动，这样一来，为求解流动复杂的流场提供 了一个有力的工具。因此，本节先介绍势流的叠加 原理，然后再介绍几种典型的有实际意义的叠加流 动。 一、势流叠加原理一、势流叠加原理 前面我们知道，速度势函数和流函数都满足拉普拉斯方 程。凡是满足拉普拉斯方程的函数，在数学分析上都称为调 和函数，所以速度势函数和流函数都是调和函数。根据调和 函数的性质，即若干个调和函数的线性组合仍然是调和函 数，可将若干个速度势函数(或流函数)线性组合成一个代表 某一有势流动的速度势函数(或流函数)。现将若干个速度势 函数 、 、 、叠加，得 （4-41） 而

39、 （4-42） 显然，叠加后新的速度势函数也满足拉普拉斯方程。同样， 叠加后新的流函数也满足拉普拉斯方程，即 （4-43） 1 2 3 321 0)( 3 2 2 2 1 2 321 22 0 3 2 2 2 1 22 这个叠加原理方法简单，在实际应用上有很大意义，可 以应用这个原理把上一节所讨论的几个简单的基本平面有势 流动叠加成所需要的复杂有势流动。 将新的速度势函数 分别对 、 和 取偏导数，就等于 新的有势流动的速度分别在 、 和 轴方向上的分量： （4-44） 或 （4-45） 即 （4-46） x y z zzzz yyyy xxxx 321 321 321 321 321 321

40、 wwww vvvv uuuu xyz 321 vvvv 由此可见，叠加后所得的复杂有势流动的 速度为叠加前原来的有势流动速度的矢量 和。 由此，可得出一个重要结论：叠加两个或多个不 可压平面势流流动组成一个新的复合流动，只要把 各原始流动的势函数或流函数简单地代数相加，就 可得到该复合流动的势函数或流函数。该结论称为 势流的叠加原理。 二、螺旋流二、螺旋流 螺旋流是点涡和点汇的叠加。将式（4-36）和式（4-32） 相加以及将式（4-37）和式（4-33）相加即得新的有势流动 的速度势和流函数 （4-47） （4-48） 式中 取反时针方向为正。于是得等势线方程 常数 或 （4-49） 流线

41、方程为 常数 或 （4-50） 显然，等势线簇和流线簇是两组互相正交的对数螺旋线簇 （图4-16），称为螺旋流。流体从四周向中心流动。 ）（rq v ln 2 1 ）（ v qrln 2 1 rq v ln v q cre 1 v qrln qv cr e 2 图4-16 螺旋流的流谱 研究螺旋流在工程上有重要意义。例如旋流燃烧室、旋 风除尘设备及多级离心泵反导叶中的旋转气流即可看成是这 种螺旋流。 螺旋流的速度分布为 （4-51） （4-52） （4-53） 代入伯努里方程（3-41），得流场的压强分布 （4-54） r r v 2 1 r q r v v r 2 22 22 222 4r

42、q vvv v r 2 2 2 1 22 2 21 11 8rr qpp v） （ 三、偶极流三、偶极流 将流量各为 的点源和 的点汇相距2a距离放在x轴 上，叠加后的流动图形如图4-17所示，它的速度势和流函数 各为 （4-55） （4-56） 由流线方程（4-56） 常数，得 常数，所以流线是经 过源点a和汇点b的圆簇，而且从源点流出的流量全部流入汇 点。 2222 21 lnln 2 lnln 2 yaxyax q rr q vv ）（）（）（ 22 22 ln 4yax yaxqv ）（ ）（ 22 21 vv qq ）（ v q v q 图4-17 点源和点汇的叠加 常数 现在分析一

43、种在点源和点汇无限接近的同时，流量无限增 大（即 ），以至使 保持一个有限常数值 的 极限情况。在这种极限情况下的流动称为偶极流， 称为偶 极矩或偶极强度。偶极流是有方向的，一般规定由点源指向 点汇的方向为正向。如图4-18所示，偶极流指向 轴方向， 这时的偶极矩 取正值。 偶极流的速度势可由式（4-55）根据上述极限条件求得， 将式（4-55）改写成 v qa，02 v aq2 m m m x 2 21 2 1 21 1ln 2 ln 2 lnln 2r rrq r rq rr q vvv ）（ 常数 常数 图4-18 偶极流的流谱 从图4-19中可知，当a点和b点向原点o无限接近 时， ，

44、而且当 ， 时 ， ， ，又由于 当 为无穷小时，可以略去高阶项，得 。因 此，偶极流的速度势 或 （4-57） 121 cos2arr 02 a v qmaqv2 rrr 21 021 432 1ln 432 ）（ )1ln( 2 1 02 2 1 02 cos2 2 lim cos2 1ln 2 lim r aq r aq v q a v q a vv 2 cos 22 cos r rm r m 222 22yx xm r xm 图4-19 推导偶极流用图 在图4-19中，bc为从b点向ap所作的垂线，则 又当 ， ， ，所以 ，代入式（4-56） 得偶极流的流函数 或 （4-58） 令式

45、（4-58）等于常数 ，于是得流线方程 （4-59） 即流线簇是半径为 、圆心为（0， ），且与轴在原 点相切的圆簇，如图4-18中实线所示。 又令式（4-57）等于常数，得等势线方程 （4-60） 即等势线簇是半径为 、圆心为（ ，0）且与轴在原 点相切的圆簇，如图4-18中虚线所示。 12 sin2sinbcar 02 a0a sin sin2ar 2 0202 sin 2 sin2 2 lim 2 lim r rm r aqq v q a v q a vv 222 22yx ym r ym 1 c 2 1 2 1 2 44 c m c m yx 1 4 c m 1 4 c m 2 2 2

46、 2 2 44 c m y c m x 2 4 c m 2 4 c m 四、绕圆柱体无环量流动四、绕圆柱体无环量流动 将均匀直线流与偶极流叠加，可以得到绕圆柱体无环量 流动。设有一在无穷远处速度 为 、平行于x轴、由左向右 流的均匀直线流，与在坐标原点o上偶极矩为m、方向与x轴 相反的偶极流叠加，如图4-20所示，组合流动的流函数为 （4-61） 流线方程 （4-62） 选取不同的常数值 ，可得到如图4-20所示的流动图形。对 的所谓零流线的方程为 或 ， v 2222 1 2 1 2yxv m yv yx ym yv c yx ym yv 22 2 c 0 c 0 1 2 1 22 yxv

47、m yv 0y v m yx 2 22 图4-20 均匀流绕圆柱体无环量流动 由此可知，零流线是一个以坐标原点为圆心、半径 的 圆周与正负x轴 和 所构成的图形。该流线到a点处分为 两段，沿上、下两个半圆周流到b点，又重新汇合。这个平 面组合流动的流函数为 (4-63) 同样，也可得到它的速度势 (4-64) 以上两式中， ，这是因为 的圆柱体内的流动没有实 际意义。 v m r 2 0 b b a a sin11 2 2 0 22 2 0 r r r v yx r yv 22 2 0 22 1 2yx r xv yx xm xv cos1 2 2 0 r r r v r 0 r 0 rr 流场中任一点的速度分量为 (4-65) 在 ， 处， ， 。这表示，在离开圆柱体无穷 远处是速度为 的均匀直线流动。在图4-20中的a点（ ， 0）和b点( ，0)处， ，a点为前驻点，b点为后驻 点。 用极坐标表示的