北航空气动力学课件第三章

《北航空气动力学课件第三章》由会员分享，可在线阅读，更多相关《北航空气动力学课件第三章（78页珍藏版）》请在金锄头文库上搜索。

1、Folie1空气动力学基础第三章理想不可压缩流体平面位流（6学时）Folie220102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课第第3 3章章 理想不可压缩流体平面位流理想不可压缩流体平面位流31 理想不可压缩流体平面位流的基本方程32 几种简单的二维位流321 直匀流322 点源323 偶极子324 点涡33 一些简单的流动迭加举例331 直匀流加点源332 直匀流加偶极子333 直匀流加偶极子加点涡34 二维对称物体绕流的数值解Folie320102010年版本年版

2、本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、理想不可压缩流体平面位流的基本方程、理想不可压缩流体平面位流的基本方程 对于理想不可压缩流体，流动的基本方程是连续方程和欧拉运动方程组。在第二章中已给出这些方程的推导过程，本章应该讨论怎样求解这些方程。但是，要想得到这些偏微分方程的解，并非易事。因为实际飞行器的外形都比较复杂，要在满足这些复杂边界条件下求得基本方程的解，困难是相当大的。 为了简化求解问题，本章首先介绍流体力学中一类简单的流动问题，理想不可压缩流体的无旋流动。 这是早期流体力学发

3、展的一种理想化近似模型，比求解真实粘性流动问题要容易的多。在粘性作用可忽略的区域，这种理想模型的解还是有相当的可信程度。 Folie420102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、理想不可压缩流体平面位流的基本方程、理想不可压缩流体平面位流的基本方程1、不可压缩理想流体无旋流动的基本方程初始条件和边界条件为在t=t0时刻，在物体的边界上在无穷远处思考：思考： 为什么需要边界条件？为什么需要边界条件？Folie520102010年版本年版本年版本年版本北京航空航天大

4、学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、平面不可压位流的基本方程、平面不可压位流的基本方程如果没有无旋条件进一步简化上述方程，求解起来也是很困难的。这是因为方程中的对流项是非线性的，而且方程中的速度 V 和压强 p 相互耦合影响，需要一并求出。但是，对于无旋流无旋流动，问题的复杂性可进一步简化，特别是可将速度和可将速度和压力分开求解力分开求解。这是因为，对于无旋运动情况，流场的速度旋度为零，即 存在速度势函数（位函数）为思考：思考： 速度和压力需要耦合求解是什么意思？速度和压力需要耦合求解是什么意思？Folie

5、620102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、平面不可压位流的基本方程、平面不可压位流的基本方程如果将上式代入不可压缩流体的连续方程中，得到：由此可见，利用无旋流无旋流动和连续条件条件所得到的这个方程是大家熟知的二二阶线性偏微分方程性偏微分方程，拉普拉斯方程拉普拉斯方程，这是一个纯运运动学方程学方程。如果对这个方程赋予适当的定解条件，就可以单独解出速度位函数，继而求出速度值。与压强p没有进行耦合求解，那么如何确定压强呢？在这种情况下，可将速度值作为已知量代入运动

6、方程中，解出p值。实际求解并不是直接代入运动方程中，而是利用Bernoulli(或Lagrange)积分得到。 Folie720102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、平面不可压位流的基本方程、平面不可压位流的基本方程 对于理想不可压缩流体，在质量力有势条件下，对于无旋流动，运动方程的积分形式为 对于定常流动，质量力只有重力，得到 如果忽略质量力（在空气动力学中经常不考虑重力）Folie820102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北

7、京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、平面不可压位流的基本方程、平面不可压位流的基本方程 由此说明，只要把速度势函数解出，压强p可直接由Bernoulli方程得到。在这种情况下整个求解步骤概括为：（1）根据纯运动学方程求出速度势函数和速度分量；（2）由Bernoulli方程确定流场中各点的压强。这使得速度和压强的求解过程分开进行，从而大大简化了问题的复杂性。综合起来对于理想不可压缩流体无旋流动，控制方程及其初边界条件为初始条件为 边界条件为 在流体力学中的边界条件多数属于第二第二类边界条件界条件，及在边界上给定速度势函数的偏导数

8、。Folie9 边界条件是在流场边界上规定的条件，边界通常分为内边界和外边界条件是在流场边界上规定的条件，边界通常分为内边界和外边界。对飞行器或物体而言，内边界即飞行器或物体表面，外边界为无边界。对飞行器或物体而言，内边界即飞行器或物体表面，外边界为无穷远。穷远。、平面不可压位流的基本方程、平面不可压位流的基本方程 （边界条件）（边界条件）按照在边界上所给条件是针对位函数自身还是位函数的法向导数，按照在边界上所给条件是针对位函数自身还是位函数的法向导数，边界条件分为三种类型：边界条件分为三种类型：（1）第一边值问题（狄利希特问题）：给出边界上位函数自身值）第一边值问题（狄利希特问题）：给出边界

9、上位函数自身值（2）第二边值问题（诺曼问题）：给出边界上位函数的法向导数值）第二边值问题（诺曼问题）：给出边界上位函数的法向导数值（3）第三边值问题（庞卡莱问题）：给出部分边界上位函数自身值，）第三边值问题（庞卡莱问题）：给出部分边界上位函数自身值，部分边界上位函数的法向导数值部分边界上位函数的法向导数值空气动力问题大多数属于第二边值问题空气动力问题大多数属于第二边值问题20102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课Folie10将坐标系与飞行器或物体固连，则外边界

10、在远离物体处，速度为将坐标系与飞行器或物体固连，则外边界在远离物体处，速度为 V ，内，内边边界是物体表面，不允界是物体表面，不允许许流体穿流体穿过过或表面法向速度或表面法向速度为为零零外外边边界界内内边边界界 n为为物面法向物面法向可以可以证证明，拉普拉斯方程的解若在明，拉普拉斯方程的解若在给给定定边边界上能界上能满满足上述条件，足上述条件，则则解是唯一的。解是唯一的。求不可求不可压压理想无旋流理想无旋流绕绕物体的流物体的流动问题动问题就就转转化化为为求解拉求解拉普拉斯方程的普拉斯方程的满满足足给给定定边边条的特解条的特解这这一数学一数学问题问题20102010年版本年版本年版本年版本北京航

11、空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课Folie1120102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、平面不可压位流的基本方程、平面不可压位流的基本方程2、速度势函数的性质 （1）速度势函数沿着某一方向的偏导数等于该方向的速度分量，速度速度势函数沿着流函数沿着流线方向增加方向增加。由此可得出，速度势函数允许相差任意常数，而不影响流体的运动。Folie1220102010年版本年

12、版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、平面不可压位流的基本方程、平面不可压位流的基本方程（2）速度势函数满足拉普拉斯方程，是调和函数。满足解的线性迭加原理。如果速度势函数满足拉普拉斯方程，则它们的线性性组合合也满足拉普拉斯方程。（3）速度势函数相等的点连成的线称为等等势线，速度方向垂直于等势线（4）连接任意两点的速度线积分等于该两点速度势函数之差。速度线积分与路径无关，仅决定于两点的位置。如果是封闭曲线，速度环量为零。Folie1320102010年版本年版本年版本年版本北京航空

13、航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、平面不可压位流的基本方程、平面不可压位流的基本方程3、流函数及其性质根据高等数学中，格林公式可知（平面问题的线积分与面积分的关系）如果令 由此可见，下列线积分与路径无关（围绕封闭曲线的线积分为零） 存在的充分必要条件是Folie1420102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、平面不可压位流的基本方程、平面不可压位流的基本方程这是不可压

14、缩流体平面流动的连续方程。这样下列微分一定是某个函数的全微分，即 这个函数称为流函数。由此可见，对于不可压缩流体的平面流平面流动（二维问题），无论是理想流体还是粘性流体，无论是有涡流动还是无涡流动，均存在流函数。思考：思考： 为什么二维流动一定存在流函数？为什么二维流动一定存在流函数？Folie1520102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、平面不可压位流的基本方程、平面不可压位流的基本方程流函数的概念是1781年Lagrange首先引进的。流函数具有下列性质（

15、1）流函数）流函数值可以差任意常数而不影响流可以差任意常数而不影响流动（2）流函数）流函数值相等的点的相等的点的连线是流是流线。即等流函数。即等流函数线的切的切线方向与方向与速度矢量方向重合速度矢量方向重合在流函数相等的线上，有上式即为平面流动的流线方程。Folie1620102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、平面不可压位流的基本方程、平面不可压位流的基本方程 （3）流函数在某一方向的偏）流函数在某一方向的偏导数等于数等于顺时针旋旋转90度度方向的速度分方向的

16、速度分量量 根据流函数这一性质，如果沿着流线取s，反时针旋转90度取n方向，则有 （4）理想不可）理想不可压缩流体平面流体平面势流，流函数流，流函数满足拉普拉斯方程足拉普拉斯方程Folie1720102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、平面不可压位流的基本方程、平面不可压位流的基本方程（5）过同一点的等速度同一点的等速度势函数函数线与等流函数与等流函数线正交（等正交（等势线与流与流线正交）正交）等流函数线是流线，有另一方面，过该点的等势函数线方程为在同一点处，流

17、线与等势线的斜率乘积为说明流线与等势线在同一点正交。Folie1820102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、平面不可压位流的基本方程、平面不可压位流的基本方程（6）流网及其特征）流网及其特征 在理想不可压缩流体定常平面势流中，每一点均存在速度势函数和流函数值。这样在流场中存在两族曲线，一族为流线，另一族为等势线，且彼此相互正交。把由这种正交曲线构成的网格叫做流网。在流网中，每一个网格的边长之比等于势函数和流函数的增值之比。 Folie1920102010年版本

18、年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、平面不可压位流的基本方程、平面不可压位流的基本方程 流网不仅可以显示流速的分布情况流速的分布情况(方向）方向），也可以反映速度的大反映速度的大小小。如流线密的地方流速大，流线稀疏的地方流速小。 如果相邻流线之间的流函数差为常数，等于单宽流量增量。即 表示流速与相邻流线的间距成反比，因此流线的疏密程度反映了速度的大小。 Folie2020102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空

19、气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、平面不可压位流的基本方程、平面不可压位流的基本方程（1）以速度势函数为未知函数的提法（2）以流函数为未知函数的提法（3）以复位势w(z)为未知函数提法 理想不可理想不可压缩流体平面定常无旋流流体平面定常无旋流动数学数学问题的提法的提法共有三种数学提法。设给定一平面物体C，无穷远为直均流，在绕流物体不脱体的情况下，求这个绕流问题。需要求解满足一定定解条件定解条件的在C外区域内的解析函数解析函数。Folie21位函数位函数的性质小结的性质小结(1)速度位函数由无旋条件无旋条件定义，位函数值可以差任意常数而不影响流动。(2) 速度位

20、函数沿着某一方向的偏导数等于该方向的速度分量，速度位函数的数值沿着流线方向增加。(3) 对于理想不可压缩无旋流动，从从连续方程出方程出发，速度位函数满足拉普拉斯方程，是调和函数，满足解的线性迭加原理。(4) 速度位函数相等的点连成的线称为等位等位线，速度方向垂直于等位线。(5) 连接任意两点的速度线积分等于该两点的速度位函数之差。速度线积分与路径无关，仅决定于两点的位置。对封闭曲线，速度环量为零。Folie22流函数流函数的性质小结的性质小结(1) 流函数由平面不可平面不可压缩流流动的连续条件条件定义，流函数值可以差任意常数而不影响流动。(2) 等流函数线是流流线。即等流函数线的切线方向与速度

21、矢量方向重合。(3) 对于理想不可压缩平面无旋流动，流函数满足拉普拉斯方程，是调和函数，解也满足叠加原理。存在条件存在条件物理意物理意义拉普拉斯方程拉普拉斯方程获得途径得途径势(位)函数无旋与流线正交始终满足连续方程流函数二维流线无旋时才满足无旋条件Folie23 位函数位函数 和流函数和流函数 之间满足之间满足柯西柯西-黎曼条件黎曼条件： 速度分量与位函数和流函数之间的关系是：速度分量与位函数和流函数之间的关系是：Folie2420102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国

22、家精品课课、几种简单的二维位流、几种简单的二维位流1、直匀流、直匀流直匀流是一种速度不变的最简单的平行流动。其流速为位函数为 常用平行于 x 轴的直匀流，从左面远方流来，流速为 。相应的流函数和势函数为Folie2520102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、几种简单的二维位流、几种简单的二维位流2、点源、点源 源可以有正负。正源是从流场上某一点有一定的流量向四面八方流开去的一种流动。负源（又名汇）是一种与正源流向相反的向心流动。如果把源放在坐在坐标原点原点上，

23、那末这流动便只有 ，而没有 。 设半径为r处的流速是 vr，那么这个源的总流量是 xy流量是常数，故流速 vr 与半径成反比。流函数的表达式是 或 Folie2620102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、几种简单的二维位流、几种简单的二维位流位函数从 的式子积分得到在极坐标系中，速度分量与流函数和势函数偏导数关系式为Folie2720102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学

24、力学力学国家精品国家精品国家精品国家精品课课、几种简单的二维位流、几种简单的二维位流如果源的位置不在坐标原点，而在A点（,）处Folie2820102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、几种简单的二维位流、几种简单的二维位流3、偶极子、偶极子 等强度的一个源和一个汇，放在 x 轴线上，源放在（-h，0）处，汇放在（0，0）处。从源出来的流量都进入汇。Folie2920102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大

25、学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、几种简单的二维位流、几种简单的二维位流应用叠加原理，位函数和流函数如下 其中 表示流场点 P 分别与源和汇连线与x轴之间的夹角。现在我们考虑一种极限情况，当h0，但同时Q增大，使Folie3020102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、几种简单的二维位流、几种简单的二维位流保持不变的极限情况。这时位函数变成： 对偶极子而言，等位线是一些圆心在 x 轴上的圆，且都过原点。Folie3120

26、102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、几种简单的二维位流、几种简单的二维位流 偶极子的流函数偶极子的流函数: 取 h0而 Qh/2=M 保持不变的极限结果，是 流线也是一些圆，圆心都在 y 轴上，且都过源点O： 两个分速的表达式是： 合速度为：Folie3220102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、几种简单的二维位流、几种简单

27、的二维位流 要注意，偶极子是源汇无限靠近的极限情况，它是有有轴线方向方向的，原来的源和汇放在哪条直线上，那条直线就是它的轴线。前面表示的偶极子是以 x 轴为轴线的，其正向为轴线上的流线方向，前面的偶极子是指向负 x 方向的。如果偶极子轴线和 x 轴成 角，正向指向第三象限，则势函数和对应流函数分别为： 如果偶极子位于（，），轴线和 x 轴成 角，正向指向第三象限，则势函数和流函数分别为： Folie3320102010年版本年版本年版本年版本、几种简单的二维位流、几种简单的二维位流4、点、点涡 点涡是位于原点的一个点涡的流动，流线是一些同心圆。流速只有 V，而没有Vr 。 式中的 是个常数，称

28、为点涡的强度，逆时针 方向为正。分速 V 和离中心点的距离 r 成反比 ，指向逆时针方向。其位函数和流函数分别为 等势线是射线，流线是圆Folie34由几何条件可立刻写出由几何条件可立刻写出 u 、 v 分量：分量：xyuvV位函数可由上式代入位函数可由上式代入 然后积分求出，但方便的还是利用然后积分求出，但方便的还是利用极座标关系：极座标关系：积分后得：积分后得：显然等位线显然等位线=C=C是是一系列射线一系列射线怎么求位函数？怎么求位函数？Folie35求流函数可由极座标下流函数与位函数的柯西黎曼关系：求流函数可由极座标下流函数与位函数的柯西黎曼关系：积分得：积分得：显然流线显然流线 =

29、C = C 是一系列同心圆，可见点涡与点源的位函数与流函数只是一系列同心圆，可见点涡与点源的位函数与流函数只是对调了一下（上述负号只是代表涡转向）。是对调了一下（上述负号只是代表涡转向）。怎么求流函数？怎么求流函数？Folie3620102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、几种简单的二维位流、几种简单的二维位流 如点涡位置不在原点，而在（，），则点涡的位函数和流函数为： 沿任意形状的围线计算环量，只要这个围线把点涡包围在内，环量值都是 ，但不包含点涡在内的围线，

30、其环量等于零。Folie3720102010年版本年版本年版本年版本、几种简单的二维位流、几种简单的二维位流 这种点涡其实应该看作是一根在 z 方向无限长的直涡线。涡本来是有旋流动，但像这样一根单独的涡线所产生的流场，除真正的涡心那一条线（在平面里就是一点）之外，其余的地方仍是无旋流动。当 r 0时，速度趋近于无穷大，相应的压强也趋于负无限大，这是不现实的。按这个速度分布规律，速度在半径方向的变化率是： * 当当 r 很小之后，很小之后，这这个个变变化率极大，化率极大，这时这时粘性力必然要粘性力必然要 起作用（粘性力与速度的法向起作用（粘性力与速度的法向变变化率成正比）。化率成正比）。结结果，

31、果， 实际涡总实际涡总是有一个核，核内流体的不是与是有一个核，核内流体的不是与 r 成反比，而成反比，而 是与是与r 成正比。但核外的流速是与成正比。但核外的流速是与 r 成反比的，如成反比的，如图图所所 示。核内是有旋流，核外是无旋流。示。核内是有旋流，核外是无旋流。这这个核的尺寸究竟个核的尺寸究竟 有多大？它是因流体的粘性大小及有多大？它是因流体的粘性大小及涡涡强强大小而不同的大小而不同的。 Folie3820102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、几种简单

32、的二维位流、几种简单的二维位流 一一般般地地说说，这这个个尺尺寸寸不不大大，我我们们作作外外部部流流场场的的计计算算时时，可可以以不不管管它它，把把它它看看作作很很微微小小就就行行了了。这这里里要要说说明明的的一一个个事事实实是是，涡涡对对于于外外部部流流场场是是产产生生诱诱导导速速度度的的（即即扰扰动动），其），其值值与至中心的距离成反比，但与至中心的距离成反比，但对对它自己的核心是没有它自己的核心是没有诱导诱导速度的。速度的。rV rVk / rr0pFolie39直匀流：xy基本解位函数、流函数小结：基本解位函数、流函数小结：abFolie4020102010年版本年版本年版本年版本北京

33、航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、 一些简单的迭加举例一些简单的迭加举例1 、直匀流加点源、直匀流加点源 在一个平行于 x 轴由左向右流去的直匀流里，加一个强度为Q的源，把坐标原点放在源所在的地方，迭加得到的位函数是： 两个分速是： 在 x 轴线上有一个合速为零的点，即驻点 A。Folie41、 一些简单的迭加举例一些简单的迭加举例令 uA=0, vA=0 即得驻点 xA 坐标为：Folie4220102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大

34、学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、 一些简单的迭加举例一些简单的迭加举例流动的流函数是：对于零流线 是一条通过坐标原点的水平线。对于 的流线方程为：得到解为： 说明是通过驻点的一条水平流线。对于非水平流线，半径 r：如对于 相应的半径 r 为：yFolie43 流线BAB的形状可以根据流函数=c 画出来，也可以从流量关系推算出来。由流函数表达： 由驻点坐标（y=0,=) 定常数c，得 cQ / 2 ，从而得流线BAB的方程为：用直角坐标表达，注意到反正切的值域为-/2,/2：Folie44 该流线与 y 轴交于 处，当 即流线在无穷远处趋于宽度为 的直

35、线。 从物理上这个结果很好理解，从源流出的流量只能限制在围线中，由速度分布知：而源的流量为Q，以速度 V 流过时将占据宽度 D=Q / V Folie4520102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、 一些简单的迭加举例一些简单的迭加举例 全部流线谱中，经过驻点A的流线BAB是一条特殊的流线， 。它像一道围墙一样，把流场划分成为两部分。外面的是直匀流绕此围墙的流动，里面的是源流在此围墙限制之内的流动。流线是气流不可逾越的线。一个物体放在气流里，它的边界也是气流不可

36、逾越的界线，气流只能与物体边界相切着流过去。所以，我们可以把外部流动看作是在直匀流中放了一个BAB那样形状的物体所造成的流动。不过这个物体后面是不封口的，称半无限体。这个半无限体在 +x 无限远处，其宽度（y 向尺寸）趋向一个渐近值 D 为：Folie4620102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、 一些简单的迭加举例一些简单的迭加举例 通常将压强表为无量纲的压强系数系数 Cp ，其定义是当地静压减去来流静压再除以来流的动压头。 不可压无粘流时： 沿这个半无限体

37、的外表面，压强系数是：why?* 压强系数与来流参数具体值压强系数与来流参数具体值 p 、V 无关，具有通用性无关，具有通用性Folie4720102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、 一些简单的迭加举例一些简单的迭加举例 首先，A点是驻点，这一点的 Cp 一定等于 +1。从驻点往后， Cp迅速下降，在距A不很远的地方， Cp 降到零，该点流速已达远前方的来流速度。此后气流继沿物面加速，走了一段之后，流速达最大值， Cp 达最小值。这一点称最大速度点，或最低压强

38、点。过了最大速度点之后气流开始减速，到无限远的右方，流速减到和远前方来流一样大 。这是大多钝头物体低速流动的特点。头部附近形成一个低速高压区，随后速度迅速上升，压强急剧下降。Folie48直均流加变强度点源直均流加等强度点源、点汇用分布的点源、点汇构造物面直均流加偶极子实验演示的直均流加点源和点汇的其他例子实验演示的直均流加点源和点汇的其他例子Folie4920102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、 一些简单的迭加举例一些简单的迭加举例 2、直匀流加偶极子（无

39、、直匀流加偶极子（无环量的量的圆柱柱绕流）流） 只有当正源和负源的总强度等于零时，物形才是封闭的。设直匀流平行于 x 轴，由左向右流。再把一个轴线指向负 x 的偶极子放在坐标原点处。这时，流动的位函数和流函数分别是： 流流动是直匀流流是直匀流流过一个一个圆。圆的半径可以从驻点点A的坐标定出来。令：why?Folie5020102010年版本年版本年版本年版本、 一些简单的迭加举例一些简单的迭加举例 得到a 就是圆半径。这样位函数和流函数可以写为： =0是一条特殊的流是一条特殊的流线。容易证明，该流线通过驻点的 x 轴线；另外还有 是半径为 a 的圆。两个速度分量为：Folie512010201

40、0年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、 一些简单的迭加举例一些简单的迭加举例在圆周上，r = a，速度分量为：相应的压强系数为：* 绕绕圆圆流流动动在在表表面面上上只只有有周周向向速速度度，没没有有径径向向速速度度Folie5220102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、 一些简单的迭加举例一些简单的迭加举例 在圆周前后驻点，=0， =，

41、压强系数等于。从前驻点往后流，在 =150 处流速加快到和来流的流速一样大了。以后继续加速，在=/2 处达最大速度，其值二倍于来流的速度，Cp 是（）。过了最大速度点以后，气流减速，在 =0处降为零，这一点称为后驻点。这个流动不仅上下是对称的，而且左右也是对称的，物面上的压强分布也是对称的，结果哪个方向的合力也没有。 不过实际流动左右是不对称的，由于实际流体是有粘性的缘故，气流过了最大速度点以后，不可能始终贴着物体流下去，不可能进行完全的减速结果水平方向是有一个阻力的 。Folie5320102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空

42、气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、 一些简单的迭加举例一些简单的迭加举例达朗培达朗培尔疑疑题 达朗培尔（DAlembert）18世纪法国著名数学家，他提出，在理想不可压流中，任何一个封闭物体的绕流，其阻力都是零。这个结论不符合事实。这个矛盾多少耽误了一点流体力学的发展，那时人们以为用无粘的位流去处理实际流动是没有什么价值的。 后来才知道，这样撇开粘性来处理问题，是一种很有价值的合乎逻辑的抽象，它能使我们把影响流动的各种因素分开来看清楚。譬如，早期由经验得出来的良好翼型，最大的升阻比不过是几十比一，后来在位流理论指导下，设计出来的翼型的最大升阻比竟达三百比一。这

43、就是无粘抽象的指导意义 。Folie5420102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课粘性流体绕圆柱的流动显示实验粘性流体绕圆柱的流动显示实验二维圆柱扰流的卡门涡街二维圆柱扰流的卡门涡街有攻角机翼绕流尾流场中的旋涡有攻角机翼绕流尾流场中的旋涡Folie5520102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、 一些简单的迭加举例一些简单的迭加举

44、例3、直匀流加偶极子加点、直匀流加偶极子加点涡（有（有环量的量的圆柱柱绕流）流） 在直匀流加偶极子直匀流加偶极子的流动之上，再在圆心处加一个强度为（ ）的点点涡（顺时针转为负）。 这时的流函数和位函数为Folie56、 一些简单的迭加举例一些简单的迭加举例在极坐标下，两个分速度为r=a仍是一条流仍是一条流线。在这个圆上 Vr =0，圆周速度为：驻点现在不在其位置可以从： ，定出来：在第三和第四象限内，前后驻点对 y 轴是对称的。这个角度离开和0的多少决定于环量对速度乘半径a之比值；比值越大，驻点越往下移。xy* 此时流函数数值不为零此时流函数数值不为零Folie57下图给出几种不同点涡强度下驻

45、点位置图画：显然，有环量的绕圆流动其左右仍是对称的，但上下已不对称了，因此在垂直于来流的 y 方向合力就不会为零。垂直于来流方向的空气动力分力称为升力，可以通过沿圆柱表面压强积分（利用伯努利方程将压强表为速度分布后积分求得）。Folie5820102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、 一些简单的迭加举例一些简单的迭加举例 用动量定理来计算绕圆柱的有环量流动的升力。以原点为中心，画一个半径为 r1很大的控制面 S，整个的控制面还包括圆的表面 S1 及连接 S 和

46、S1 的两条割线。不过这两条割线上的压力和动量进出 都对消了，不必管它受力情况左右对称，不 会有 X合力。我们只只计算算 Y 方向合力方向合力就行了。 彻体力略去不计；流流动是定常是定常的。Folie5920102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、 一些简单的迭加举例一些简单的迭加举例动量积分方程变为：在 r1 的大圆上： LpLvFolie6020102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空

47、气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、 一些简单的迭加举例一些简单的迭加举例在上述表达式中，奇函数积分为零，只有偶函数积分。 对于单位时间动量的净流出量计算如下：Folie61、 一些简单的迭加举例一些简单的迭加举例在 y 方向的速度分量是：单位长度圆柱所受的总升力为：Folie6220102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、 一些简单的迭加举例一些简单的迭加举例 只要是一个封闭物体，代表这个物体作用的正负源的强度总和必须等于零。这种正负源

48、放在一起的情况，在远离物体的地方（我们可以取r1 很大），其作用和一个偶极子没什么区别。这就说明了物形对升力没有直接的关系，关键的问题在于必须有一个绕物体的环量存在。有了环量又有一个直匀流，便有一个升力。 库塔塔-儒可夫斯基定理儒可夫斯基定理 一个封闭物体所受升力L等于来流的密度 乘速度再乘以气流方向逆着环流旋转90。Folie6320102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课、 一些简单的迭加举例一些简单的迭加举例 环量之所以能产生一个 Y 向的合力，也可以从圆柱

49、体上的压力分布直接看到。其中有环量和无环量绕流情况作了对比。无环量时，上半圆（由至0）上的压力分布和下半圆（由至2）上的压力分布对称，结果是合力为零。有环量时，上半圆上的负压远远超过下半圆上的负压，所以有一个向上的合力，即升力。这个力的来源主要靠上半圆上的吸力 。Folie64从“野渡无人舟自横”到“香蕉球”技术浅谈“香蕉球”的力学原理Folie6520102010年版本年版本年版本年版本3 34 4 二维对称物体绕流的数值解二维对称物体绕流的数值解 把直匀流和分布的偶极子（或总强度为零的分布点源和点汇）叠加起来，所得到的组合流动为对称封闭物体绕流。设直匀流沿 x 轴正向流来，其速度为V，在

50、x 轴上 x =a和 x =b 范围内连续分布一系列的偶极子，单位长度内偶极子的强度设为m（偶极子密度偶极子密度）。 如果偶极子密度的分布形式已知，则离原点距离为 的小区间内由偶极子产生的流函数为： 总流函数为： 。物体的外形可以用零流线来表示。改变不同的偶极子密度分布，可以获得不同形状的封闭物体，由流函数和速度以及速度与压强的关系确定流场中各点及物体表面的速度分布和压强分布。 奇点叠加数奇点叠加数值解法解法Folie6620102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课

51、课3 34 4 二维对称物体绕流的数值解二维对称物体绕流的数值解Folie6720102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课3 34 4 二维对称物体绕流的数值解二维对称物体绕流的数值解 对于实际问题，往往是给定物体的外形来确定其流动的特性。在这种情况下，偶极子密度分布函数的确定需要由流函数求解。对偶极子密度来说，流函数是一个积分方程，求它的解是比较困难的。但是随着计算机技术的发展，可以用数值方法比较迅速地获得这种方程的有一定准确度的数值解。 下面简单地叙述用数值方

52、法求解已知物体形状确定绕物体流动特性的过程。Folie6820102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课3 34 4 二维对称物体绕流的数值解二维对称物体绕流的数值解1、数、数值解法步解法步骤 首先，我们把偶极子分布区域分成等宽度的 n 段，设每段的宽度为，段数 n 可根据计算机容量及结果的准确度要求而确定。流场中某一定点P处的流函数为： 式中 为第 j 段中点离原点的距离； 为第 j 段内偶极子密度的平均值； 表示第 j 段内偶极子的强度。Folie6920102

53、010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课3 34 4 二维对称物体绕流的数值解二维对称物体绕流的数值解用物面用物面边界条件来确定待求的偶极子密度界条件来确定待求的偶极子密度对于给定物体外形上的n个已知点（xi，yi），就可以得到一个对未知函数的 n 元一次联立代数方程组： 其中 为影响系数，表示 处的单位偶极子 密度对物体表面某点Pi (xi, yi) 处的流函数贡献。Folie7020102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学

54、北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课3 34 4 二维对称物体绕流的数值解二维对称物体绕流的数值解展开上式，即 利用解一次方程组的各种计算方法，求解上面方程组，确定偶极子密度。Folie7120102010年版本年版本年版本年版本北京航空航天大学北京航空航天大学北京航空航天大学北京航空航天大学空气空气空气空气动动力学力学力学力学国家精品国家精品国家精品国家精品课课3 34 4 二维对称物体绕流的数值解二维对称物体绕流的数值解 一旦所给定物体外形的偶极子密度分布已经解得，则可以确定流场内任意点处的流函数。此后即可由流函数与速度的关系式及伯努利方

55、程，确定流场内各点处的速度及压强值。 在上述过程中，我们实际上是把第j段中分布的偶极子用集中在该段中点处的等强度偶极子来代替了。显然，如果分段数量较多，这种近似表示才有一定的准确性。理论上，当段数 n 趋于无限大时，偶极子密度分布的数值结果趋近于精确解。在实际应用时，由于计算机容量和计算机机时的限制，以及多元一次联立方程组解的不稳定性。分段的数目不宜太多。 也可以由位函数出发，用位函数对应的物面条件来解决实际流动问题。这两种方法是等价的。在实际应用中，用位函数叠加法比用流函数法更广泛。 Folie72本章基本要求本章基本要求掌握平面不可压位流中位函数与流函数的性质与关系。掌握平面不可压位流的基

56、本方程即拉普拉斯方程的特点、叠加原理和边界条件；掌握四种基本而重要的位流流动即：直匀流，点源（点汇）、偶极子和点涡的表达；重点掌握直匀流与偶极子和点涡的叠加；掌握儒可夫斯基升力定律；了解二维对称物体绕流数值解法步骤Folie73 本章本章讨论讨论怎怎样样求解不可求解不可压压理想流体无旋运理想流体无旋运动动的的规规律律在理想不可压条件下欧拉方程和连续方程包括四个方程和四个未知函数（u,v,w,p），理论上是可解的由于飞行器的外形都比较复杂，要在满足如此复杂的边界条件下求该偏微分方程组的解析解是非常困难的，原因在于方程包含非线性项，而且方程中速度与压强相互耦合，需要一并求出人们发现在无旋条件下问题

57、可以得到大大简化，尤其是可以将速度和压强分开求解，这是因为无旋条件可使关于速度位的方程化为线性方程，从而便于单独求得速度位即求出速度，而压强可利用伯努利方程求解本章的思路是，先针对理想不可压无旋流求得一些典型的速度位基本解，将这些基本解进行叠加得到满足非常简单边界条 件的流动。对复杂外形的绕流，介绍用基本解进行叠加的数值解法大意Folie74Folie75小测验（15分钟）1.a. 试写出从 y 流向 +y ，速度值为 V 的直匀流的位函数。 b. 试写出位于原点的点汇的位函数。 c. 试写出位于原点，顺时针旋转的点涡的位函数。 d. 试写出位于原点，轴线指向y轴的偶极子的位函数。 提示：坐标

58、旋转变换关系为： 2.在 x 轴上距原点为b处放置点源（源强度为Q），y 轴为璧面，试求x=b/2 处的速度表达。3.同上图，在 x 轴上距原点为b处放置强度为顺时针旋转的点涡，求该点涡此时的速度大小与方向。xyyxxyo.b1.d 图2题.3题图Folie76解解1a. V Vy y b. c. d. xyyxFolie77解2. 在 x 轴上距原点为b处放置点源（源强度为Q），y 轴为璧面，试求x=b/2 处的速度表达。相当于将璧面去掉在b处放置强度同样为Q的点源：xyo.b2题.3题图.Folie78解解2：题目相当于在：题目相当于在b处叠加一个环量相等方向相反的点处叠加一个环量相等方向相反的点涡，求该点涡速度相当于求涡，求该点涡速度相当于求b点速度。由于点涡对自身没点速度。由于点涡对自身没有诱导速度，因此有诱导速度，因此b处速度只是由镜像涡产生，镜像涡的处速度只是由镜像涡产生，镜像涡的流函数（方便求导）：流函数（方便求导）：xyo.b2题.3题图.如果求如果求b/2处速度，则需写出全部流函数：处速度，则需写出全部流函数：