空气动力学课件第3章

空气动力学基础第三章理想不可压缩流体平面位流2010年版本北京航空航天大学《空气动力学》北京市精品课第3章理想不可压缩流体平面位流3．1理想不可压缩流体平面位流的基本方程3．2几种简单的二维位流3．2．1直匀流3．2．2点源3．2．3偶极子3．2．4点涡3．3一些简单的流动迭加举例3．3．1直匀流加点源3．3．2直匀流加偶极子3．3．3直匀流加偶极子加点涡3．4二维对称物体绕流的数值解2010年版本北京航空航天大学《空气动力学》北京市精品课3.1、理想不可压缩流体平面位流的基本方程

对于理想不可压缩流体，流动的基本方程是连续方程和欧拉运动方程组。在第二章中已给出这些方程的推导过程，本章应该讨论怎样求解这些方程。但是，要想得到这些偏微分方程的解，并非易事。因为实际飞行器的外形都比较复杂，要在满足这些复杂边界条件下求得基本方程的解，困难是相当大的。为了简化求解问题，本章首先介绍流体力学中一类简单的流动问题，理想不可压缩流体的无旋流动。2010年版本北京航空航天大学《空气动力学》北京市精品课3.1、理想不可压缩流体平面位流的基本方程这是早期流体力学发展的一种理想化近似模型，比求解真实粘性流动问题要容易的多。在粘性作用可忽略的区域，这种理想模型的解还是有相当的可信程度。2010年版本北京航空航天大学《空气动力学》北京市精品课3.1、理想不可压缩流体平面位流的基本方程1、不可压缩理想流体无旋流动的基本方程初始条件和边界条件为在t=t0时刻，在物体的边界上在无穷远处2010年版本北京航空航天大学《空气动力学》北京市精品课3.1、平面不可压位流的基本方程如果没有无旋流条件进一步简化上述方程，求解起来也是很困难的。这是因为方程中的对流项是非线性的，而且方程中的速度V和压强p相互偶合影响，需要一并求出。但是，对于无旋流动，问题的复杂性可进一步简化，特别是可将速度和压力分开求解。这是因为，对于无旋运动情况，流场的速度旋度为零，即存在速度势函数（位函数）为2010年版本北京航空航天大学《空气动力学》北京市精品课3.1、平面不可压位流的基本方程如果将上式代入不可压缩流体的连续方程中，得到

2010年版本北京航空航天大学《空气动力学》北京市精品课3.1、平面不可压位流的基本方程由此可见，利用无旋流动和连续条件所得到的这个方程是大家熟知的二阶线性偏微分方程，拉普拉斯方程，这是一个纯运动学方程。如果对这个方程赋予适定的定解条件，就可以单独解出速度位函数，继而求出速度值。与压强p没有进行偶合求解，那么如何确定压强呢？在这种情况下，可将速度值作为已知量代入运动方程中，解出p值。实际求解并不是直接代入运动方程中，而是利用Bernoulli(或Lagrange)积分得到。2010年版本北京航空航天大学《空气动力学》北京市精品课3.1、平面不可压位流的基本方程对于理想不可压缩流体，在质量力有势条件下，对于无旋流动，运动方程的积分形式为

对于定常流动，质量力只有重力，得到如果忽略质量力（在空气动力学中经常不考虑重力的作用）

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

在流体力学中的边界条件多数属于第二类边界条件，及在边界上给定速度势函数的偏导数。2010年版本北京航空航天大学《空气动力学》北京市精品课3.1、平面不可压位流的基本方程

2、速度势函数的性质

（1）速度势函数沿着某一方向的偏导数等该方向的速度分量，速度势函数沿着流线方向增加。由此可得出，速度势函数允许相差任意常数，而不影响流体的运动。

2010年版本北京航空航天大学《空气动力学》北京市精品课3.1、平面不可压位流的基本方程（2）速度势函数满足拉普拉斯方程，是调和函数。满足解的线性迭加原理。如果速度势函数满足拉普拉斯方程，则它们的线性组合也满拉普拉斯方程。

2010年版本北京航空航天大学《空气动力学》北京市精品课3.1、平面不可压位流的基本方程（3）速度势函数相等的点连成的线称为等势线，速度方向垂直于等势线。

（4）连接任意两点的速度曲线等于该两点的速度势函数之差。速度线积分与路径无关，仅决定于两点的位置。如果是封闭曲线，速度环量为零。2010年版本北京航空航天大学《空气动力学》北京市精品课3.1、平面不可压位流的基本方程

3、流函数及其性质流函数的概念是1781年Lagrange首先引进的。流函数具有下列性质（1）流函数值可以差任意常数而不影响流动（2）流函数值相等的点的连线是流线。即等流函数线的切线方向与速度矢量方向重合。在流函数相等的线上，有上式即为平面流动的流线方程。（3）流函数在某一方向的偏导数顺时针旋转90度方向的速度分量。2010年版本北京航空航天大学《空气动力学》北京市精品课3.1、平面不可压位流的基本方程

根据流函数这一性质，如果沿着流线取s，反时针旋转90度取n方向，则有（流函数增值方向沿速度方向反时针旋转90度方向）（4）理想不可压缩流体平面势流，流函数满足拉普拉斯方程。即2010年版本北京航空航天大学《空气动力学》北京市精品课3.1、平面不可压位流的基本方程2010年版本北京航空航天大学《空气动力学》北京市精品课3.1、平面不可压位流的基本方程2010年版本北京航空航天大学《空气动力学》北京市精品课3.1、平面不可压位流的基本方程（3）流网及其特征在理想不可压缩流体定常平面势流中，每一点均存在速度势函数和流函数值。这样在流场中存在两族曲线，一族为流线，另一族为等势线且彼此相互正交。把由这种正交曲线构成的网格叫做流网。在流网中，每一个网格的边长之比等于势函数和流函数的增值之比。如网格正方形。2010年版本北京航空航天大学《空气动力学》北京市精品课3.1、平面不可压位流的基本方程流网不仅可以显示流速的分布情况，也可以反映速度的大小。如流线密的地方流速大，流线稀疏的地方流速小。

2010年版本北京航空航天大学《空气动力学》北京市精品课3.1、平面不可压位流的基本方程如果相邻流线之间的流函数差为常数，等于单宽流量增量。即表示流速与网格间距成反比，因此流线的疏密程度反映了速度的大小。5、理想不可压缩流体平面定常无旋流动数学问题的提法对于理想不可压缩平面定常无旋流动问题的数学提法共有三种。设给定一平面物体C，无穷远为直均流，在绕流物体不脱体的情况下，求这个绕流问题。2010年版本北京航空航天大学《空气动力学》北京市精品课3.1、平面不可压位流的基本方程（1）以速度势函数为未知函数的提法

（2）以流函数为未知函数的提法（3）以复位势w(z)为未知函数提法需要求解满足一定定解条件的在C外区域内的解析函数。2010年版本北京航空航天大学《空气动力学》北京市精品课3.2、几种简单的二维位流

1、直匀流直匀流是一种速度不变的最简单的平行流动。其流速为位函数为常用平行于x轴的直匀流，从左面远方流来，流速为。相应的流函数和势函数为2010年版本北京航空航天大学《空气动力学》北京市精品课3.2、几种简单的二维位流

2、点源源可以有正负。正源是从流场上某一点有一定的流量向四面八方流开去的一种流动。负源（又名汇）是一种与正源流向相反的向心流动。如果把源放在坐标原点上，那末这流动便只有υr，而没有。设半径为r处的流速是υr，那末这个源的总流量是2010年版本北京航空航天大学《空气动力学》北京市精品课3.2、几种简单的二维位流流量是常数，故流速υr与半径成反比。流函数的表达式是或2010年版本北京航空航天大学《空气动力学》北京市精品课3.2、几种简单的二维位流位函数从的式子积分得到在极坐标系中，速度分量与流函数和势函数偏导数关系式为2010年版本北京航空航天大学《空气动力学》北京市精品课3.2、几种简单的二维位流如果源的位置不在坐标原点，而在A（ξ,η）处2010年版本北京航空航天大学《空气动力学》北京市精品课3.2、几种简单的二维位流

3、偶极子等强度的一个源和一个汇，放在x轴线上，源放在（-h，0）处，汇放在（0，0）处。从源出来的流量都进入汇。2010年版本北京航空航天大学《空气动力学》北京市精品课3.2、几种简单的二维位流应用叠加原理，位函数和流函数如下

其中表示流场点P分别与源和汇连线与x轴之间的夹角。现在我们考虑一种极限情况，当h→0，但同时Q增大，使2010年版本北京航空航天大学《空气动力学》北京市精品课3.2、几种简单的二维位流保持不变的极限情况。这时位函数变成偶极子。等位线是一些圆心在x轴上的圆，且都过原点。2010年版本北京航空航天大学《空气动力学》北京市精品课3.2、几种简单的二维位流流函数的式子，取h→0而Qh/2π=M保持不变的极限结果，是2010年版本北京航空航天大学《空气动力学》北京市精品课3.2、几种简单的二维位流流线也是一些圆，圆心都在y轴上，且都过源点O。两个分速的表达式是:合速度为2010年版本北京航空航天大学《空气动力学》北京市精品课3.2、几种简单的二维位流要注意，偶极子是源汇无限靠近的极限情况，它是有轴线方向的，原来的源和汇放在哪条直线上，那条直线就是它的轴线。前面表示的偶极子是以x轴为轴线的，其正向为轴线上的流线方向，前面的偶极子是指向负x方向的。如果偶极子轴线和x轴成θ角，正向指向第三象限，则势函数为

2010年版本北京航空航天大学《空气动力学》北京市精品课3.2、几种简单的二维位流相应的流函数为如果偶极子位于（ξ，η），轴线和x轴成θ角，正向指向第三象限，则势函数和流函数分别为2010年版本北京航空航天大学《空气动力学》北京市精品课3.2、几种简单的二维位流

4、点涡点涡是位于原点的一个点涡的流动，流线是一些同心圆。流速只有，而没有。式中的是个常数，称为点涡的强度，反时针方向为正。分速和离中心点的距离r成反比，指向是反时针方向的。其位函数和流函数分别为（等势线是射线，流线是圆）2010年版本北京航空航天大学《空气动力学》北京市精品课3.2、几种简单的二维位流如果点涡的位置不在原点，而在（ξ，η），则点涡的位函数和流函数分别是沿任意形状的围线计算环量，值都是，只要这个围线把点涡包围在内，但不包含点涡在内的围线，其环量等于零。2010年版本北京航空航天大学《空气动力学》北京市精品课3.2、几种简单的二维位流

这种点涡其实应该看作是一根在z方向无限长的直涡线。涡本来是有旋流动，但像这样一根单独的涡线所产生的流场，除真正的涡心那一条线（在平面里就是一点）之外，其余的地方仍是无旋流动。当r→0时，速度趋近于无穷大，相应的压强也趋于负无限大，这是不现实的。按这个速度分布规律，速度在半径方向的变化率是

2010年版本北京航空航天大学《空气动力学》北京市精品课3.2、几种简单的二维位流当r很小之后，这个变化率极大，这时粘性力必然要起作用（粘性力与速度的法向变化率成正比）。结果，实际涡总是有一个核，核内流体的不是与r成反比，而是与r成正比。但核外的流速是与r成反比的，如图所示。核内是有旋流，核外是无旋流。这个核的尺寸究竟有多大？它是因流体的粘性大小及涡强大小而不同的。2010年版本北京航空航天大学《空气动力学》北京市精品课3.2、几种简单的二维位流一般地说，这个尺寸不大，我们作外部流场的计算时，可以不管它，把它看作很微小就行了。这里要说明的一个事实是，涡对于外部流场是产生诱导速度的（即扰动），其值与至中心的距离成反比，但对它自己的核心是没有诱导速度的。2010年版本北京航空航天大学《空气动力学》北京市精品课小结2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例

1、直匀流加点源在一个平行于x轴由左向右流去的直匀流里，加一个强度为Q的源，把坐标原点放在源所在的地方，迭加得到的位函数是

两个分速是在x轴线上有一个合速为零的点，即驻点A。2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例令

即得驻点xA坐标为2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例2、直匀流加偶极子（无环量的圆柱绕流）只有当正源和负源的总强度等于零时，物形才是封闭的。设直匀流平行于x轴，由左向右流。再把一个轴线指向负x的偶极子放在坐标原点处。这时，流动的位函数是流动是直匀流流过一个圆。圆的半径可以从驻点A的坐标定出来。令2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例

得到

a就是圆半径。这样位函数可以写成为流函数方程为2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例

ψ=0是一条特殊的流线。容易证明，该流线通过驻点的x轴线；另外还有是半径为a的圆。两个速度分量为2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例在圆周上，r=a，速度分量为相应的压强系数为2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例在圆周前后驻点，θ=0°，θ=π，压强系数等于1.0。从前驻点往后流，在θ=150°处流速加快到和来流的流速一样大了。以后继续加速，在θ=π/2处达最大速度，其值二倍于来流的速度，Cp是（–3.0）。过了最大速度点以后，气流减速，在θ=0°处降为零，这一点称为后驻点。这个流动不仅上下是对称的，而且左右也是对称的，物面上的压强分布也是对称的，结果哪个方向的合力也没有。2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例不过实际流动左右是不对称的，由于实际流体是有粘性的缘故，气流过了最大速度点以后，不可能始终贴着物体流下去，不可能进行完全的减速结果水平方向是有一个阻力的。达朗培尔疑题 达朗培尔（D’Alembert）18世纪法国著名数学家，他提出，在理想不可压流中，任何一个封闭物体的绕流，其阻力都是零。这个结论不符合事实。这个矛盾多少耽误了一点流体力学的发展，那时人们以为用无粘的位流去处理实际流动是没有什么价值的。

2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例

3、直匀流加偶极子加点涡（有环量的圆柱绕流）在直匀流加偶极子的流动之上，再在圆心处加一个强度为（–）的点涡（顺时针转为负）。这时的流函数和位函数为2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例在极坐标下，两个分速度为r=a仍是一条流线。在这个圆上Vr=0，圆周速度为

驻点现在不在其位置可以从2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例定出来在第三和第四象限内，前后驻点对y轴是对称的。这个角度离开π和0°的多少决定于环量对速度乘半径a之比值；比值越大，驻点越往下移。

2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例2010年版本北京航空航天大学《空气动力学》北京市精品课3.3、一些简单的迭加举例2010年版本北京航空航天大学《空气动力学》北京市精品课3．4二维对称物体绕流的数值解

把直匀流和分布的偶极子（或总强度为零的分布点源和点汇）叠加起来，所得到的组合流动为对称封闭物体绕流。设直匀流沿x轴正向流来，其速度为V∞，在x轴上x=a和x=b范围内连续分布一系列的偶极子，单位长度内偶极子的强度设为m（偶极子密度）。这样组合的流函数为如果偶极子密度的分布形式已知，则离原点距离为ξ的小区间内由偶极子产生的流函数为2010年版本北京航空航天大学《空气动力学》北京市精品课3．4二维对称物体绕流的数值解总流函数为物体的外形可以用零流线来表示。改变不同的偶极子密度分布，可以获得不同形状的封闭物体，由流函数和速度以及速度与压强的关系确定流场中各点及物体表面的速度分布和压强分布。2010年版本北京航空航天大学《空气动力学》北京市精品课3．4二维对称物体绕流的数值解对于实际问题，往往是给定物体的外形来确定其流动的特性。在这种情况下，偶极子密度分布函数的确定需要由流函数求解。对偶极子密度来说，流函数是一个积分方程，求它的解是比较困难的。但是随着计算机技术的发展，可以用数值方法比较迅速地获得这种方程的有一定准确度的数值解。

2010年版本北京航空航天大学《空气动力学》北京市精品课3．4二维对称物体绕流的数值解下面简单地叙述用数值方法求解已知物体形状确定绕物体流动特性的过程。首先，我们把偶极子分布区域分成等宽度的n段，设每段的宽度为△ξ，段数n可根据计算机容量及结果的准确度要求而确定。2010年版本北京航空航天大学《空气动力学》北京市精品课3．4二维对称物体绕流的数值解1、数值解法步骤首先，我们把偶极子分布区域分成等宽度的n段，设每段的宽度为△ξ，段数n可根据计算机容量及结果的准确度要求而确定。流场中某一定点P处的流函数为式中为第j段中点离原点的距离；为第j段内偶极子密度的平均值；表示第j段内偶极子的强度。2010年版本北京航空航天大学《空气动力学》北京市精品课3．4二维对称物体绕流的数值解用物面边界条件来确定待求的偶极子密度 对于给定物体外形上的n个已知点（xi，yi），就可以得到一个对未知函数的n元一次联立代数方程组 其中

为影响系数，表示（xi，yi）处的单位偶极子密度对物体表面某点Pi处的流函数贡献。2010年版本北京航空航天大学《空气动力学》北京市精品课3．4二维对称物体绕流的数值解展开上式，即

利用解一次方程组的各种计算方法，求解上面方程组，确定偶极子密度。2010年版本北京航空航天大学《空气动力学》北京市精品课3．4二维对称物体绕流的数值解一旦所给定物体外形的偶极子密度分布已经解得，则可以确定流场内任意点处的流函数。此后即可由流函数与速度的关系式及伯努利方程，确定流场内各点处的速度及压强值。在上述过程中，我们实际上是把第j段中分布的偶极子用集中在该段中点处的等强度偶极子来代替了。显然，如果分段数量较多，这种近似表示才有一定的准确性。理论上，当段数n趋于无限大时，