《理想不可压缩流体的有旋流动和无旋流动.pptx》由会员分享，可在线阅读，更多相关《理想不可压缩流体的有旋流动和无旋流动.pptx（81页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第七章 理想不可压缩流体的有旋流动和无旋流动本章内容：在许多工程实际问题中，流动参数不仅在流动方向上发生变化，而且在垂直于流动方向的横截面上也要发生变化。要研究此类问题，就要用多维流的分析方法。本章主要讨论理想流体多维流动的基本规律，为解决工程实际中类似的问题提供理论依据，也为进一步研究粘性流体多维流动奠定必要的基础。第1页/共81页第一节 流体流动的连续性方程 当把流体的流动看作是连续介质的流动，它必然遵守质量守恒定律。对于一定的控制体，必须满足式（322）。它表示在控制体内由于流体密度变化所引起的流体质量随时间的变化率等于单位时间内通过控制体的流体质量的净通量。首先推导在笛卡儿坐标系中微分
2、形式的连续性方程。图7-1 微元六面体 第2页/共81页 设该微元六面体中心点O（x,y,z）上流体质点的速度为 、，密度为 ，于是和 轴垂直的两个平面上的质量流量如图所示。在 方向上，单位时间通过EFGH面流入的流体质量为：（a）单位时间通过ABCD面流出的流体质量：（b）则在 方向单位时间内通过微元体表面的净通量为（b）-（a），即（c1）第3页/共81页同理可得 和 方向单位时间通过微元体表面的净通量分别为:（c2）（c3）因此，单位时间流过微元体控制面的总净通量为:（c）微元六面体内由于密度随时间的变化而引起的质量的变化率为：（d d）第4页/共81页 将式（c），（d）代入式（7-1
3、），取 0，则可得到流场中任一点的连续性方程的一般表达式为:或（7-17-1）（7-1a）连续性方程表示了单位时间内控制体内流体质量的增量等于流体在控制体表面上的净通量。它适用于理想流体和粘性流体、定常流动和非定常流动。在定常流动中，由于 （7-27-2）（7-37-3）或（7-3a）对于不可压缩流体（=常数）第5页/共81页在其它正交坐标系中流场中任一点的连续性方程和柱坐标系中的表示式为：(7-4)对于不可压缩流体(7-4a)式中 为极径；为极角。球坐标系中的表示式为:(7-5)（7-5a）式中 为径矩；为纬度；为径度。第6页/共81页【例7-1】已知不可压缩流体运动速度已知不可压缩流体运动
4、速度 在在 ，两个轴方向的分量为两个轴方向的分量为 ，。且在。且在 处，有处，有 。试求。试求 轴方向的速度分量轴方向的速度分量 。【解】对不可压缩流体连续性方程为：【解】对不可压缩流体连续性方程为：将已知条件代入上式，有将已知条件代入上式，有 又由已知条件对任何又由已知条件对任何 ，当，当 时，时，。故有。故有 第7页/共81页第二节 流体微团的运动分析 流体与刚体的主要不同在于它具有流动性，极易变形。因此，流体微团在运动过程中不但象刚体那样可以有移动和转动，而且还会发生变形运动。一般情况下，流体微团的运动可以分解为移动，转动和变形运动。图7-2 流体微团运动速度分量 第8页/共81页 如图
5、7-2所示，在流场中任取一微元平行六面体，其边长分别为 dx、dy、dz，微元体中心点沿三个坐标轴的速度分量为 、。顶点E的速度分量可按照泰勒级数展开，略去二阶以上无穷小项求得，如图。为了简化讨论，先分析流体微团的平面运动，如图7-3。该平面经过微元平行六面体的中心点且平行于xoy面。由于流体微团各个点的速度不一样，在dt时间间隔中经过移动、转动和变形运动（包括角变形运动和线变形运动），流体微团的位置和形状都发生了变化。具体分析如下：图7-3 流体微团的平面运动第9页/共81页（1）移动：由图7-3看出，A、B、C、D各点速度分量中都含有 、项，如果只考虑这两项，则经过时间dt，矩形ABCD向
6、右移动 的距离，向上移动 的距离。移动到新位置后，形状保持不变，如图7-4(a)所示。（2）线变形运动：如果只考虑AB边和CD边在x轴方向上的速度差 ，则经过时间dt，AD边和BC边在x轴方向上伸长了 的距离；如果只考虑AD边和BC边在y轴方向上的速度差 ，则经过时间dt，根据连续性条件，AB边和CD边在y轴方向上缩短了 的距离，这就是流体微团的线变形，如图7-4(b)。每秒钟单位长度的伸长或缩短量称为线应变速度，在x轴方向的线应变速度分量为：同样可得在y轴方向和z轴方向的分量分别为 、。第10页/共81页图7-4 流体微团平面运动的分析 第11页/共81页 （3）角变形运动和旋转运动：如图7
7、-4（c）、（d）所示，取图7-3中的 来分析。如果只考虑B点和A在y轴方向上的速度差 ，则经过时间dt，B点运动到B点，运动距离为 ，使AB边产生了角变形运动，变形角度为 ；如果只考虑D点和A点在x轴方向上的速度差 ，则经过时间dt，D点运动到D点，运动距离为 ,使AD边产生了角变形运动，变形角度为 。变形角可按下列公式求得。变形角速度为：第12页/共81页 上面只考虑了角变形运动，实际上流体微团在运动中变形和旋转是同时完成的。设流体微团旋转角度为 ，变形角度为 ，如图7-4(d)所示 由 、式可得：如果 ，则 ，也就是只发生了角变形运动，矩形变成了平行四边形。如果 ，则 ，矩形ABCD各边
8、都向逆时针旋转了同一微元角度 ，矩形只发生旋转运动，形状不变。一般情况是 ，即 ，矩形ABCD在发生旋转运动的同时，还要发生角变形运动，结果也变成了平行四边形。第13页/共81页 在旋转运动中，流体微团的旋转角速度定义为每秒内绕同一转轴的两条互相垂直的微元线段旋转角度的平均值。于是流体微团沿z轴的旋转角速度分量：同理，可求得流体微团沿x轴和y轴的旋转角速度分量和 。于是，流体微团的旋转角速度分量为：写成矢量形式为：(7-6)（7-8）第14页/共81页 在角变形运动中，流体微团的角变形速度定义为每秒内一个直角的角度变化量，则在xoy面内的角变形是 。于是流体微团在垂直于z轴的平面上的角变形速度
9、分量 ，即 同样可求得在垂直于x轴和y轴的平面上的角变形速度分量之半 和 。于是，流体微团的角变形速度之半的分量是:(7-9)如果在式（7-10）的第一式右端加入两组等于零的项和，其第15页/共81页值不变。经过简单组合，可将该式写成:同理，有:将式（7-6），（7-9）代入以上三式，得 第16页/共81页上式中，各速度分量的第一项是移动速度分量，第二、三、四项分别是由线变形运动、角变形运动和旋转运动所引起的线速度分量。此关系也称为亥姆霍兹(Helmholtz)速度分解定理，该定理可简述为：在某流场O点邻近的任意点A上的速度可以分成三个部分，与O点相同的平移速度（移运）；绕O点转动在A点引起的
10、速度（旋转运动）；由于变形（包括线变形和角变形）在A点引起的速度（变形运动）。亥姆霍兹速度分解定理对于流体力学的发展有深远的影响。由于把旋转运动从一般运动中分离出来，才使我们有可能把运动分成无旋运动和有旋运动。正是由于把流体的变形运动从一般运动中分离出来，才使我们有可能将流体变形速度与流体应力联系起来，这对于粘性流体运动规律的研究有重大的影响。第17页/共81页第三节 有旋流动和无旋流动 根据流体微团在流动中是否旋转，可将流体的流动分为两类：有旋流动和无旋流动。数学条件：当当无旋流动有旋流动 通常以 是否等于零作为判别流动是否有旋或无旋的判别条件。在笛卡儿坐标系中：（7-11）第18页/共81
11、页即当流场速度同时满足：时流动无旋 需要指出的是，有旋流动和无旋流动仅由流体微团本身是否发生旋转来决定，而与流体微团本身的运动轨迹无关。如图7-5（a），流体微团的运动为旋转的圆周运动，其微团自身不旋转，流场为无旋流动；图7-5（b）流体微团的运动尽管为直线运动，但流体微团在运动过程中自身在旋转，所以，该流动为有旋流动。（a）（b）图7-5 流体微团运动轨迹 第19页/共81页【例7-2】某一流动速度场为 ，其中 是不为零的常数，流线是平行于 轴的直线。试判别该流动是有旋流动还是无旋流动。【解】由于所以该流动是有旋运动。第20页/共81页第四节 理想流体运动微分方程式欧拉积分和伯努里积分 一、
12、运动微分方程 理想流体运动微分方程式是研究流体运动学的重要理论基础。可以用牛顿第二定律加以推导。在流场中取一平行六面体，如图76所示。其边长分别为dx，dy，dz，中心点为A(x,y,z)。中心点的压强为p=p(x,y,z)，密度为=(x,y,z)。因为研究的对象为理想流体，作用于六个面上的表面力只有压力，作用于微元体上的单位质量力 沿三个坐标轴的分量分别为。图76 理想流体运动微分方程用图 第21页/共81页 微元体在质量力和表面力的作用下产生的加速度 ，根据牛顿第二定律：两端同除以微元体的质量 ,并整理有：(7-12)写成矢量式：(7-13)第22页/共81页将加速度的表达式代入（712）
13、有：（714）其矢量式为：（715）公式（714）为理想流体运动微分方程式，物理上表示了作用在单位质量流体上的质量力、表面力和惯性力相平衡。该式推导过程中对流体的压缩性没加限制，故可适用于理想的可压流体和不可压缩流体，适用于有旋流动和无旋流动。第23页/共81页 将（714）作恒等变形，便可以直接由运动微分方程判定流动是有旋还是无旋流动，在式（7-14）的第一式右端同时加减 、，得:由式（7-8）得:（7-16）写成矢量形式（7-17）第24页/共81页 如果流体是在有势的质量力作用下，流场是正压性的，则：此时存在一压强函数:（718）将压强函数对坐标的偏导数有:将上述关系代入式（7-16），
14、得:（7-19）写成矢量形关系式（7-20）理想正压性流体在有势的质量力作用下的运动微分关系理想正压性流体在有势的质量力作用下的运动微分关系第25页/共81页二、欧拉积分 当理想正压性流体在有势的质量力作用下作定常无旋流动时，式（7-19）右端为零。若在流场中任取一有向微元线段 ，其在三个坐标轴的投影分别为dx、dy、dz，将它们分别依次乘式（7-19）并相加，得:积分（7-21）上式为欧拉积分的结果，表明理想正压性流体在有势的质量力作用下作定常无旋流动时，单位质量流体的总机械能在流场中保持不变。第26页/共81页三、伯努里积分 当理想正压性流体在有势的质量力作用下作定常有旋流动时，式（7-1
15、9）右端第一项等于零。由流线的特性知，此时流线与迹线重合，在流场中沿流线取一有向微元线段 ，其在三个坐标轴上的投影分别为 ，将它们的左、右端分别依次乘式（7-19）的左、右端，相加有 积分有（7-22）该积分为伯努里积分。表明理想正压性流体在有势的质量力作用下作定常有旋流动时，单位质量流体的总机械能沿流线保持不变。通常沿不同流线积分常数值有所不同。第27页/共81页第五节 理想流体的旋涡运动 本节主要讲述理想流体有旋运动的理论基础，重点是速度环量及其表征环量和旋涡强度间关系的斯托克斯定理。一、涡线、涡管、涡束和旋涡强度 涡量用来描述流体微团的旋转运动。涡量的定义为：(7-23)涡量是点的坐标和
16、时间的函数。它在直角坐标系中的投影为：(7-24)在流场的全部或部分存在角速度的场，称为涡量场。如同在速度场中引入了流线、流管、流束和流量一样。在涡量场中同样也引入涡线、涡管、涡束和旋涡强度的概念。第28页/共81页1涡线:涡线是在给定瞬时和涡量矢量相切的曲线。如图7-7所示。图7-7 涡线 图7-8 涡管根据涡通量矢量与涡线相切的条件，涡线的微分方程为:(7-25)2涡管、涡束：在涡量场中任取一不是涡线的封闭曲线，在同一时刻过该曲线每一点的涡线形成的管状曲面称作涡管。如图7-8所示。截面无限小的涡管称为微元涡管。涡管中充满着的作旋转运动的流体称为涡束，微元涡管中的涡束称为微元涡束或涡丝。第2
17、9页/共81页3旋涡强度（涡通量）在涡量场中取一微元面积dA，见图7-9（a），其上流体微团的涡通量为 ，为dA的外法线方向，定义（7-26）为任意微元面积dA上的旋涡强度，也称涡通量。任意面积A上的旋涡强度为：（7-27）如果面积A是涡束的某一横截面积，就称为涡束旋涡强度，它也是旋转角速度矢量的通量。旋涡强度不仅取决于，而且取决于面积A。第30页/共81页二、二、二、二、速度环量、斯托克斯定理速度环量、斯托克斯定理速度环量、斯托克斯定理速度环量、斯托克斯定理1速度环量:在流场的某封闭周线上，如图7-9（b），流体速度矢量沿周线的线积分，定义为速度环量，用符号 表示，即:(7-28)速度环量是
18、一代数量，它的正负与速度的方向和线积分的绕行方向有关。对非定常流动，速度环量是一个瞬时的概念，应根据同一瞬时曲线上各点的速度计算，积分时为参变量。图7-9微元面积、微元有向线段 第31页/共81页 2斯托克斯（Stokes）定理:在涡量场中，沿任意封闭周线的速度环量等于通过该周线所包围曲面面积的旋涡强度，即:（7-29）这一定理将旋涡强度与速度环量联系起来，给出了通过速度环量计算旋涡强度的方法。【例7-3】已知二维流场的速度分布为 ，试求绕圆 的速度环量。【解】此题用极坐标求解比较方便，坐标变换为:速度变换为,第32页/共81页【例7-4】一二维元涡量场，在一圆心在坐标原点、半径 的圆区域内，
19、流体的涡通量 。若流体微团在半径 处的速度分量 为常数，它的值是多少？【解】由斯托克斯定理得：第33页/共81页三、三、汤姆孙定理、亥姆霍兹定理汤姆孙定理、亥姆霍兹定理1汤姆孙（Thomson）定理 理想正压性流体在有势的质量力作用下，沿任何封闭流体周线的速度环量不随时间变化，即：证明：在流场中任取一由流体质点组成的封闭周线K，它随流体的运动而移动变形，但组成该线的流体质点不变。沿该线的速度环量可表示为式（7-28），它随时间的变化率为：（7-30）第34页/共81页（7-30a）由于质点线K始终由同样的流体质点组成，将其代入式（7-30a）等号右端第一项积分式：第35页/共81页由理想流体的
20、欧拉运动微分方程，式（7-30a）等号右端第二项积分式可表示为:将上面的结果代入式（7-30a），并考虑到 都是单值连续函数，得:（7-30b）或斯托克斯定理和汤姆孙定理表明，理想正压性流体在有势的质量力作用下，涡旋不会自行产生，也不会自行消失。第36页/共81页 2亥姆霍兹（Helmholtz）定理亥姆霍兹关于旋涡的三个定理，解释了涡旋的基本性质，是研究理想流体有旋流动的基本定理。（1）亥姆霍兹第一定理：在理想正压性流体的有旋流场中，同一涡管各截面上的旋涡强度相同。如图710所示，在同一涡管上任取两截面A1、A2，在A1、A2之间的涡管表面上取两条无限靠近的线段a1a2和b1b2。由于图71
21、0 同一涡管上的两截面图711 涡管上的封闭轴线第37页/共81页 封闭周线a1a2b1b2a1所围成的涡管表面无涡线通过，旋涡强度为零。根据斯托克斯定理，沿封闭周线的速度环量等于零，即：由于 而 ，故得 该定理说明，在理想正压性流体中，涡管既不能开始，也不能终止。但可以自成封闭的环形涡管，或开始于边界、终止于边界。第38页/共81页（2）亥姆霍兹第二定理（涡管守恒定理）亥姆霍兹第二定理（涡管守恒定理）理想正压性流体在有势的质量力作用下，流场中的涡管始终由相同的流体质点组成。如图711所示，K为涡管表面上的封闭周线，其包围的面积内涡通量等于零。由斯托克斯定理知，周线K上的速度环量应等于零；又由
22、汤姆孙定理，K上的速度环量将永远为零，即周线K上的流体质点将永远在涡管表面上。换言之，涡管上流体质点将永远在涡管上，即涡管是由相同的流体质点组成的，但其形状可能随时变化。第39页/共81页（3）亥姆霍兹第三定理（涡管强度守恒定理）亥姆霍兹第三定理（涡管强度守恒定理）理想正压性流体在有势的质量力作用下，任一涡管强度不随时间变化。若周线K为包围涡管任意的截面A的边界线。由汤姆孙定理知，该周线上的速度环量为常数。根据斯托克斯定理截面A上的旋涡强度为常数。因为A为任意截面，所以整个涡管各个截面旋涡强度都不瞬时间发生变化，即涡管的旋涡强度不随时间变化。由亥姆霍兹三定理可知，粘性流体的剪切应力将消耗能量，
23、使涡管强度逐渐减弱。第40页/共81页第六节第六节 二维旋涡的速度和压强分布二维旋涡的速度和压强分布假设在理想不可压缩的重力流体中，有一象刚体一样以等角速度 绕自身轴旋转的无限长铅垂直涡束，其涡通量为J。涡束周围的流体在涡束的诱导下绕涡束轴等速圆周运动，由斯托克斯定理知 。由于直线涡束无限长，该问题可作一个平面问题研究。可以证明涡束内的流动为有旋流动，称为涡核区，其半径为 ；涡束外的流动区域为无旋流动，称为环流区。在环流区内，速度分布为:在环流区内，压强分布由伯努里方程式导出。环流区内半径为 的点和无穷远处的伯努里方程:（7-31）第41页/共81页式中的 即为 ，为无穷远处的压强。将 代入上
24、式得:由上式可知，在涡束外部的势流区内，随着环流半径的减小，流速上升而压强降低；在涡束边缘上，流速达该区的最高值，而压强则是该区的最低值，即：涡束内部的速度分布为:（7-32）（7-33）第42页/共81页由于涡束内部为有旋流动，伯努利积分常数随流线变化，故其压强分布可由欧拉运动微分方程导出。对于平面定常流动，欧拉运动微分方程为:将涡核内任意点的速度投影到直角坐标上，则有，代入上式得:将 和 分别乘以以上二式，相加后得:第43页/共81页或积分得:在与环流区交界处，代入上式，得积分常数：得涡核区的压强分布为：（7-30）第44页/共81页由上式可知涡管中心的压强最低，其大小为 ，涡核区边缘至涡
25、核中心的压强差为 由以上讨论可知，涡核区和环流区的压强差相等，其数值均为 。涡核区的压强比环流区的的低。在涡束内部,半径愈小,压强愈低，沿径向存在较大的压强梯度，所以产生向涡核中心的抽吸作用，涡旋越强，抽吸作用越大。自然界中的龙卷风和深水旋涡就具有这种流动特征，具有很大的破坏力。在工程实际中有许多利用涡流流动特性装置，如锅炉中的旋风燃烧室、离心式除尘器、离心式超声波发生器、离心式泵和风机、离心式分选机等。第45页/共81页第七节第七节 速度势和流函数速度势和流函数一 速度势函数对于无旋流场，处处满足：，由矢量分析知，任一标量函数梯度的旋度恒为零，所以速度 一定是某个标量函数 的梯度，即:因 则
26、有:即流场的速度等于势函数 的梯度。因此，称 为速度势函数，简称速度势；称无旋流动为有势流动，简称势流。这与单位质量有势力和有势力场的势函数的关系相类似。（735）（736）第46页/共81页证明：结论：无旋条件是速度有势的充要条件。无旋必然有势，有势必须无旋。所以无旋流场又称为有势流场。速度势的存在与流体是否可压缩、流动是否定常无关。在笛卡儿坐标系中：，由 则 ，代入 或 ，有所以 得证第47页/共81页以上给出了在直角坐标系中速度势函数和速度的关系，在柱坐标系中，有势流动的速度势函数与速度的线积分有密切关系。若势流中有一曲线AB，速度沿该曲线积分为上式表明，有势流动中沿AB曲线的速度线积分
27、等于终点B和起点A的速度势之差。由于速度势是单值的，则该线积分与积分路径无关。这与力做的功和位势的关系相类似。当速度沿封闭轴线积分时即，周线上的速度环量等于零。（7-33）（7-34）（7-35）第48页/共81页 根据无旋条件，速度有势：代入不可压缩连续性条件可得：或上述方程称作不可压无旋流动的基本方程。在笛卡儿坐标系中：在柱坐标系中：式中 为拉普拉斯算子。满足拉普拉斯方程的函数为调和函数，故速度势是调和函数。（7-36）（7-37）（7-38）第49页/共81页二二 流函数流函数在笛卡儿坐标系中，平面、不可压缩流体的连续性方程可写成:若定义某一个函数（流函数）令:平面不可压缩流体流函数的基
28、本性质1、等流函数线为流线当 常数时即：（7-39）,第50页/共81页2、流体通过两流线间单位高度的体积流量等于两条流流体通过两流线间单位高度的体积流量等于两条流线的流函数之差线的流函数之差在xy平面上任取A和B点，AB连线如图7-12所示，则（AB为使 与 同号）不论是理想流体还是粘性流体，不论是有旋的还是无旋的流动，只要是不可压缩（或定常可压缩）流体的平面（或轴对称）流动，就存在流函数。图7-12 流量与流函数的关系第51页/共81页由不可压缩流体、平面、无旋流动条件有：将速度和流函数的关系代入上式得在极坐标系中：故不可压缩流体的平面无旋流动流函数也满足拉普拉斯方程，也是调和函数。（7-
29、40）（7-41）第52页/共81页三三 速度势函数和流函数的关系速度势函数和流函数的关系 对于不可压缩流体的平面无旋流动，速度势函数和流函数都是调和函数，且具有以下关系：该数学关系式称为柯西黎曼（CauchyRiemen）条件。由它可得：两族曲线的正交条件。在平面上它们构成处处正交的网络，称为流网。第53页/共81页【例7-6】已知不可压缩流体平面势流，其速度势，试求速度投影和流函数。【解】由速度势可求得速度分量 ，由速度和流函数的关系 ，将速度代入流函数的关系式积分得 将上式对求偏导数，并考虑速度和流函数的关系则有：上式对积分，得：代入原式有：第54页/共81页第八节第八节 几种简单的平面
30、势流几种简单的平面势流 一、均匀等速流一、均匀等速流二、点二、点源和点汇源和点汇 三、点涡第55页/共81页一一 均匀等速流均匀等速流流速的大小和方向沿流线不变的流动为均匀流；若流线平行且流速相等，则称均匀等速流。例如 ，其中 为常数，便是这样的流动。由于积分得:（7-43a）由于 积分得 （7-43b）在以上二式中均取积分常数为零（下同），这对流动的计算并无影响。第56页/共81页显然，等势线 与流线 是相互垂直的两族直线，如图7-13所示。若已知来流速度 与x 轴的夹角 ，则有:由于流场中各点的速度相同，流动无旋，故处处有 常数，即在流场中各点的总势能保持不变。若是水平面上的均匀等势流，或
31、者不计重力的影响（例如大气），则p=常数，即压强在流场中处处相等。图图7-13 均匀等速流均匀等速流（7-43d）（7-43c）第57页/共81页二 点源和点汇 无限大平面上，流体从一点沿径向直线均匀地向外流出的流动，称为点源，这个点称为源点；如果流体沿径向均匀的流向一点，称为点汇，这个点称为汇点。不论是点源还是点汇，流场中只有径向速度，即图7-14 源流和汇流(a)(a)(b)(b)第58页/共81页根据流体的连续性原理，在极坐标中流体流过任意单位高度圆柱面的体积流量（也称为源流或汇流的强度）都相等，即 上式中点源取正号，点汇取负号。根据上式，只是 的函数，所以 积分得以上讨论表明，当 时，
32、源点和汇点是奇点，以上 和 只有在 0时才有意义。流函数和速度的关系为:（7-44a），第59页/共81页因此，只是 的函数，故有 上式积分得根据以上得到的流函数和势函数可知，等势线为不同半径的同心圆，即 =常数；流线为不同极角的径线，即 =常数。在水平面 面上，对半径 处和无穷远处列伯努利方程 代入速度值后由上式可知，压强随着半径的减小而降低。零压强处的半径为 。以上各式仅适用于 的区域。（7-44b）（7-44c）第60页/共81页三三 点点涡涡若直线涡束的半径 ，则垂直于该涡束的平面内的流动称为点涡或自由涡流，涡流中心称为涡点。涡点以外势流区的速度分布仍为 由以上关系式知，时，所以涡点为
33、奇点，该式仅适用于 区域。由此式可见，只是 的函数。故有 积分得速度和流函数的关系为上式表明 只是 的函数，所以（7-45a）图7-15 点涡 第61页/共81页 上式积分得 由上可知，点涡流场的等势线为不同极角的径线，即 =常数；流线为不同半径的同心圆，即 =常数。与点源（或点汇）相反。点涡的强度即沿围绕点涡轴线上的环量 0时，环流为逆时针方向；0，环流为顺时针方向。由斯托克斯定理知，点涡的强度 取决于旋涡的强度。涡点以外势流区的压强和前述二维涡流流场压强分布相同，其分布关系仍为式（7-32）。零压强处的半径为 上述各式的实际适用范围为 的区域。以上几种简单的平面势流实际中很少应用，但它们是
34、势流的基本单元，若把几种基本单元叠加在一起，可以形成许多有实际意义的复杂流动。（7-45b）第62页/共81页第九节 简单平面势流的叠加 一、汇流和点涡叠加的流动螺旋流 二、源流和汇流叠加的流动偶极子流 几个简单有势流动叠加得到的新的有势流动，其速度势函数和流函数分别等于原有几个有势流动的速度势函数和流函数的代数和，速度分量为原有速度分量的代数和。研究势流叠加原理的意义：将简单的势流叠加起来，得到新的复杂流动的流函数和势函数，可以用来求解复杂流动。第63页/共81页一一 汇流和点涡叠加的流动汇流和点涡叠加的流动螺旋流螺旋流若点源和点涡均位于坐标原点，组成一新的流场，其速度势和流函数为（7-48
35、）（7-49）（7-50）（7-51）图图7-16 螺旋流网螺旋流网 令以上的速度势和流函数为常数，得到的等势线和流线方 程分别为：v其图像为图71 6所示，等势线和流线是两组相v互正交的对数螺旋线，故称汇流和点涡叠加的流动v为螺旋流。其速度分布为:v其适用范围应为：v压强分布可用前述方法导出，表达式为第64页/共81页二二 源流和汇流叠加的流动源流和汇流叠加的流动偶极子流偶极子流（7-52）（7-53）图图7-17 点源和点汇叠加点源和点汇叠加 图图7-18 偶极流偶极流 v组合流动的速度势和流函数为v两个强度 相等的位于点A(-a,0)的点源和位于点B(a,0)的点汇叠加，如图717所示。
36、由于 是AP、BP之间的夹角，在流线上 =常数，=常数。其图像为经过源点和汇点的圆线族第65页/共81页 当 时，源点和汇点无限接近，流量为无限增大，使得 取有限值，称这种流动为偶极流。M为偶极子矩，其方向由源点指向汇点。当 为微量时，故由式（7-52）（7-53）可得偶极流的速度势和流函数分别为 即 即 （7-54）（7-55）第66页/共81页 若令式（7-54）等于常数 ，则得等势线方程即等势线的图像为圆心在（）点上，半径为 并与y轴在原点相切的圆族，如图7-18中虚线所示。令式（7-55）式等于常数 时，可得流线方程：即流线的图像是圆心为（）.半径为 并与x轴在原点相切的圆族，如图7-
37、18中实线所示。对速度势函数求偏导数，得出的偶极流的速度分布为（7-56），第67页/共81页第十节第十节 平行流绕过圆柱体无环流的平面流动平行流绕过圆柱体无环流的平面流动 平行流（均匀等速流）和偶极流叠加,可用来描述流体绕过圆柱体无环流的流动.若均匀等速流的速度为 ,沿x轴正向流动,偶极流的偶极矩为M。一、平行流与偶极流的叠加1.流网 平行流：偶极流：叠加：（757）（758）第68页/共81页 流线方程为：流线方程为：当常数C取不同的数值时，可得如图719所示的流普。当C0时对应的流线，称为零流线。图719流体对圆柱体的无环量绕流 2、零流线 当常数C0时，即零流线的流线方程：由 ，得 。
38、或 即：可见，零流线为以坐标原点为圆心，为半径的圆和x轴。第69页/共81页二、平行绕流圆柱体无环流的流动二、平行绕流圆柱体无环流的流动1、流函数和速度势：2、流场中的速度分析（1）直角坐标系：因为：所以：（759a）（759b）（）（）b:在（r0，0）和（r0，0）处a:当讨论：时，即为平行流。为驻点，即A，B为驻点。第70页/共81页（2）对于极坐标：对于极坐标：讨论：（760）a:半径为r的圆形曲线上的速度环量b：当 时，故平行流绕圆柱体的流动为势流。；时当 时，即C、D点的速度最大（如图720）。图720（761）第71页/共81页三、圆柱面上的压强分布三、圆柱面上的压强分布 圆柱面
39、上的压强分布可由伯努利方程求得。在无穷远处，速度为 ，压强为 。则 工程上为了处理问题方便起见，引入一个无量纲压强系数 ，则 。由其中讨论：1、前、后驻点：2、C、D点：3、在 和 的范围内，圆柱面上的压强作用是对称的，即作用在其上的压力是平衡的。（762）第72页/共81页四、对于理想流体，平行流无环流绕流圆柱体时，既不产四、对于理想流体，平行流无环流绕流圆柱体时，既不产生阻力也无升力。生阻力也无升力。证明：如图，在单位长度的圆柱体上作用在微元弧段上的总压力和阻力分别为 证毕。图721 说明：无升力、无阻力只适用于理想流体，实际流体不适用，上述即为达朗伯疑题。第73页/共81页第十一节第十一
40、节 平行流绕流圆柱体有环流的平面流平行流绕流圆柱体有环流的平面流动动库塔儒可夫斯基公式库塔儒可夫斯基公式一、平行流绕流圆柱体有环流的流动是无环流流动和一个环流 的叠加。1、平行流绕圆柱体无环流的流动其中2、环流第74页/共81页 3、叠加的结果、叠加的结果 流场如图所示。上部和环流方向一致，速度加快，下部方向相反，速度减慢，上部压强降低，下部升高。（763）（764）图722 平行流与纯环流的叠加则可得（765）（766）第75页/共81页4、验证是否满足两个边界条、验证是否满足两个边界条件件（1）是否满足圆柱面为流线的条件由式（763）（765）知：当 时，令又，当 时，故满足以圆柱面为边界
41、的流动。则（2）是否满足来流速度为 的边界条件 由式（763）知：当 时，。故满足无穷远处的条件因此这种叠加是正确的。第76页/共81页二、圆柱面上的速度分布及驻点的位置二、圆柱面上的速度分布及驻点的位置 1、圆柱面上的速度分布 由式（765）（766）知：当 时，该式说明：圆柱面上径向速度为零，即无分离，切向速度为 的正弦函数。2、驻点的位置 （1）时 当驻点在圆柱面上时，此时，讨论：当 时，且 则有两个驻点。因 即 和 。且随着 增大，也增大，驻点向中间移动。如图（723a）（767）（768）第77页/共81页当当 时，时，驻点移动到最下驻点移动到最下方。如图方。如图（723b）当 时，
42、。可令 和 为零，得（）（）两个驻点，一个在圆柱体内，如图（723c）当 时，和上述的情况类似，只是驻点的位置在上部。（c）（b）（a）图723 圆柱面上的驻点位置第78页/共81页三、圆柱面上的压强分布三、圆柱面上的压强分布1、压强分布在圆柱面上 列无穷远处和圆柱面上的伯努利方程（769）则可得当 、为常数时，2、阻力、升力 库塔儒可夫斯基公式作用在单位长度圆柱体上的阻力和升力为：第79页/共81页意义：理想流体绕圆柱体有环流的流动中，在垂直于来流意义：理想流体绕圆柱体有环流的流动中，在垂直于来流方向上，流体作用在单位长度的圆柱体上的升力等于流体方向上，流体作用在单位长度的圆柱体上的升力等于流体的密度、来流速度和环量的乘积。升力的方向为的密度、来流速度和环量的乘积。升力的方向为 的方向的方向反环流转反环流转 。式771即为库塔（美国）儒可夫斯基（俄罗斯）公式。（771）（770）图724 升力的方向第80页/共81页感谢您的观看！第81页/共81页