《北航空气动力学课件》由会员分享，可在线阅读，更多相关《北航空气动力学课件（76页珍藏版）》请在金锄头文库上搜索。
1、Folie1空气动力学基础空气动力学基础 第三章理想不可压缩流体平面位流第三章理想不可压缩流体平面位流（6学时）Folie220102010年版本年版本北京航空航天大学北京航空航天大学空气动力学空气动力学国家精品课国家精品课第3章 理想不可压缩流体平面位流31 理想不可压缩流体平面位流的基本方程 32 几种简单的二维位流 321 直匀流 322 点源 323 偶极子 324 点涡 33 一些简单的流动迭加举例 331 直匀流加点源 332 直匀流加偶极子 333 直匀流加偶极子加点涡 34 二维对称物体绕流的数值解Folie320102010年版本年版本北京航空航天大学北京航空航天大学空气动力
2、学空气动力学国家精品课国家精品课3.1、理想不可压缩流体平面位流的基本方程对于理想不可压缩流体，流动的基本方程是连续方程和欧拉运动方程组。在第二章中已给出这些方程的推导过程，本章应该讨论怎样求解这些方程。但是，要想得到这些偏微分方程的解，并非易事。因为实际飞行器的外形都比较复杂，要在满足这些复杂边界条件下求得基本方程的解，困难是相当大的。为了简化求解问题，本章首先介绍流体力学中一类简单的流动问题，理想不可压缩流体的无旋流动。 这是早期流体力学发展的一种理想化近似模型，比求解真实粘性流动问题要容易的多。在粘性作用可忽略的区域，这种理想模型的解还是有相当的可信程度。 Folie420102010年
3、版本年版本北京航空航天大学北京航空航天大学空气动力学空气动力学国家精品课国家精品课3.1、理想不可压缩流体平面位流的基本方程1、不可压缩理想流体无旋流动的基本方程初始条件和边界条件为在t=t0时刻，在物体的边界上在无穷远处思考： 为什么需要边界条件？Folie520102010年版本年版本北京航空航天大学北京航空航天大学空气动力学空气动力学国家精品课国家精品课3.1、平面不可压位流的基本方程如果没有无旋条件进一步简化上述方程，求解起来也是很 困难的。这是因为方程中的对流项是非线性的，而且方程中的 速度 V 和压强 p 相互耦合影响，需要一并求出。但是，对于无旋流动，问题的复杂性可进一步简化，特
4、别是可将速度和压力 分开求解。这是因为，对于无旋运动情况，流场的速度旋度为 零，即存在速度势函数（位函数）为思考： 速度和压力需要耦合求解是什么意思？Folie620102010年版本年版本北京航空航天大学北京航空航天大学空气动力学空气动力学国家精品课国家精品课3.1、平面不可压位流的基本方程如果将上式代入不可压缩流体的连续方程中，得到：由此可见，利用无旋流动和连续条件所得到的这个方程是大家熟知的二阶线性偏微分方程，拉普拉斯方程，这是一个纯运动学方程。如果对这个方程赋予适当的定解条件，就可以单独解出速度位函数，继而求出速度值。与压强p没有进行耦合求解，那么如何确定压强呢？在这种情况下，可将速度
5、值作为已知量代入运动方程中，解出p值。实际求解并不是直接代入运动方程中，而是利用Bernoulli(或Lagrange)积分得到。 Folie720102010年版本年版本北京航空航天大学北京航空航天大学空气动力学空气动力学国家精品课国家精品课3.1、平面不可压位流的基本方程由此说明，只要把速度势函数解出，压强p可直接由Bernoulli方程得到。在这种情况下整个求解步骤概括为：（1）根据纯运动学方程求出速度势函数和速度分量；（2）由Bernoulli方程确定流场中各点的压强。这使得速度和压强的求解过程分开进行，从而大大简化了问题的复杂性。综合起来对于理想不可压 缩流体无旋流动，控制方程及其初
6、边界条件为初始条件为边界条件为在流体力学中的边界条件多数属于第二类边界条件，及在边界上给定速度势函数的偏导数。Folie8边界条件是在流场边界上规定的条件，边界通常分为内边界和外边界。对飞行器或物体而言，内边界即飞行器或物体表面，外边界为无 穷远。3.1、平面不可压位流的基本方程 （边界条件）按照在边界上所给条件是针对位函数自身还是位函数的法向导数，边界条件分为三种类型：（1）第一边值问题（狄利希特问题）：给出边界上位函数自身值（2）第二边值问题（诺曼问题）：给出边界上位函数的法向导数值（3）第三边值问题（庞卡莱问题）：给出部分边界上位函数自身值，部分边界上位函数的法向导数值空气动力问题大多数
7、属于第二边值问题20102010年版本年版本北京航空航天大学北京航空航天大学空气动力学空气动力学国家精品课国家精品课Folie9将坐标系与飞行器或物体固连，则外边界在远离物体处，速度为 V ，内边边界是物体表面，不允许许流体穿过过或表面法向速度为为零外边边界内边边界 n为为物面法向可以证证明，拉普拉斯方程的解若在给给定边边界上能满满足上述条件，则则 解是唯一的。求不可压压理想无旋流绕绕物体的流动问题动问题 就转转化为为求解拉 普拉斯方程的满满足给给定边边条的特解这这一数学问题问题20102010年版本年版本北京航空航天大学北京航空航天大学空气动力学空气动力学国家精品课国家精品课Folie102
8、0102010年版本年版本北京航空航天大学北京航空航天大学空气动力学空气动力学国家精品课国家精品课3.1、平面不可压位流的基本方程2、速度势函数的性质（1）速度势函数沿着某一方向的偏导数等于该方向的速度分量，速度势函数沿着流线方向增加。由此可得出，速度势函数允许相差任意常数，而不影响流体的运动。Folie1120102010年版本年版本北京航空航天大学北京航空航天大学空气动力学空气动力学国家精品课国家精品课3.1、平面不可压位流的基本方程（2）速度势函数满足拉普拉斯方程，是调和函数。满足解的线性迭加原理。如果速度势函数满足拉普拉斯方程，则它们的线性组合也满足拉普拉斯方程。（3）速度势函数相等的
9、点连成的线称为等势线，速度方向垂直于等势线（4）连接任意两点的速度线积分等于该两点速度势函数之差。速度线积分与路径无关，仅决定于两点的位置。如果是封闭曲线，速度环量为零。Folie1220102010年版本年版本北京航空航天大学北京航空航天大学空气动力学空气动力学国家精品课国家精品课3.1、平面不可压位流的基本方程3、流函数及其性质根据高等数学中，格林公式可知（平面问题的线积分与面积分的关系）如果令由此可见，下列线积分与路径无关（围绕封闭曲线的线积分为零）存在的充分必要条件是Folie1320102010年版本年版本北京航空航天大学北京航空航天大学空气动力学空气动力学国家精品课国家精品课3.1
10、、平面不可压位流的基本方程这是不可压缩流体平面流动的连续方程。这样下列微分一定是某个函数的全微分，即这个函数称为流函数。由此可见，对于不可压缩流体的平面流动（二维问题），无论是理想流体还是粘性流体，无论是有涡流动还是无涡流动，均存在流函数。思考： 为什么二维流动一定存在流函数？Folie1420102010年版本年版本北京航空航天大学北京航空航天大学空气动力学空气动力学国家精品课国家精品课3.1、平面不可压位流的基本方程流函数的概念是1781年Lagrange首先引进的。流函数具有下列性质（1）流函数值可以差任意常数而不影响流动（2）流函数值相等的点的连线是流线。即等流函数线的切线方向与速度矢
11、量方向重合在流函数相等的线上，有上式即为平面流动的流线方程。Folie1520102010年版本年版本北京航空航天大学北京航空航天大学空气动力学空气动力学国家精品课国家精品课3.1、平面不可压位流的基本方程（3）流函数在某一方向的偏导数等于顺时针旋转90度方向的速度分量根据流函数这一性质，如果沿着流线取s，反时针旋转90度取n方向，则有（4）理想不可压缩流体平面势流，流函数满足拉普拉斯方程Folie1620102010年版本年版本北京航空航天大学北京航空航天大学空气动力学空气动力学国家精品课国家精品课3.1、平面不可压位流的基本方程（5）过同一点的等速度势函数线与等流函数线正交（等势线与流线正
12、交）等流函数线是流线，有另一方面，过该点的等势函数线方程为在同一点处，流线与等势线的斜率乘积为说明流线与等势线在同一点正交。Folie1720102010年版本年版本北京航空航天大学北京航空航天大学空气动力学空气动力学国家精品课国家精品课3.1、平面不可压位流的基本方程（6）流网及其特征在理想不可压缩流体定常平面势流中，每一点均存在速度势函数和流函数值。这样在流场中存在两族曲线，一族为流线，另一族为等势线，且彼此相互正交。把由这种正交曲线构成的网格叫做流网。在流网中，每一个网格的边长之比等于势函数和流函数的增值之比。 Folie1820102010年版本年版本北京航空航天大学北京航空航天大学空
13、气动力学空气动力学国家精品课国家精品课3.1、平面不可压位流的基本方程流网不仅可以显示流速的分布情况(方向），也可以反映速度的大小。如流线密的地方流速大，流线稀疏的地方流速小。如果相邻流线之间的流函数差为常数，等于单宽流量增量。即表示流速与相邻流线的间距成反比，因此流线的疏密程度反映了速度的大小。Folie1920102010年版本年版本北京航空航天大学北京航空航天大学空气动力学空气动力学国家精品课国家精品课3.1、平面不可压位流的基本方程（1）以速度势函数为未知函数的提法（2）以流函数为未知函数的提法（3）以复位势w(z)为未知函数提法理想不可压缩流体平面定常无旋流动数学问题的提法共有三种数
14、学提法。设给定一平面物体C，无穷远为直均流，在绕流物体不脱体的情况下，求这个绕流问题。需要求解满足一定定解条件的在C外区域内的解析函数 。Folie20位函数的性质小结(1)速度位函数由无旋条件定义，位函数值可以差任意常数而不影响流动。(2) 速度位函数沿着某一方向的偏导数等于该方向的速度分量，速度位函数的数值沿着流线方向增加。(3) 对于理想不可压缩无旋流动，从连续方程出发，速度位函数满足拉普拉斯方程，是调和函数，满足解的线性迭加原理。(4) 速度位函数相等的点连成的线称为等位线，速度方向垂直于等位线。(5) 连接任意两点的速度线积分等于该两点的速度位函数之差。速度线积分与路径无关，仅决定于
15、两点的位置。对封闭曲线，速度环量为零。Folie21位函数 和流函数 之间满足柯西-黎曼条件：速度分量与位函数和流函数之间的关系是：Folie2220102010年版本年版本北京航空航天大学北京航空航天大学空气动力学空气动力学国家精品课国家精品课3.2、几种简单的二维位流1、直匀流直匀流是一种速度不变的最简单的平行流动。其流速为位函数为常用平行于 x 轴的直匀流，从左面远方流来，流速为 。相应的流函数和势函数为Folie2320102010年版本年版本北京航空航天大学北京航空航天大学空气动力学空气动力学国家精品课国家精品课3.2、几种简单的二维位流2、点源源可以有正负。正源是从流场上某一点有一定的流量向四面八方流开去的一种流动。负源（又名汇）是一种与正源流向相反的向心流动。如果把源放在坐标原点上，那末这流动便只有 ，而没有 。设半径为r处的流速是 vr，那么这个源的总流量是 xy流量是常数，故流速 vr 与半径成反比。流函数的表达式是 或 Folie2420102010年版本年版本北京航空航天大学北京航空航天大学空气动力学空气动力学国家精品课国家精品课3.2、几种简单的二维位流位函数从 的式子积分得到在极坐标系中，速度分量与流函数和势函数偏导数关系式为Folie2520102010年版本年版本北京航空航天大学北京航空航天大
