理想流体的有旋流动和无旋流动.

1、第八章 理想流体的有旋 流动和无旋流动主要内容 理想流体理想流体微分微分形式的基本方程形式的基本方程( (连续方程、连续方程、运动方程运动方程) ) 流体微团运动分析（平移运动、变形运动）流体微团运动分析（平移运动、变形运动） 二维势流以及叶栅、叶型绕流的升力计算二维势流以及叶栅、叶型绕流的升力计算为工程实践提供理论依据；同时是研究黏性流体多维流动的基础为工程实践提供理论依据；同时是研究黏性流体多维流动的基础 两种方法 （1）微元控制体分析法）微元控制体分析法 （2）有限控制体分析法）有限控制体分析法一、微分形式的连续方程一、微分形式的连续方程oyxzdmxdmxdxdydzdtdt时间内时间

2、内x x方向：方向：流入质量流入质量流出质量流出质量净流出质量净流出质量dydzdtvdmxxdydzdtdxx)v(vdmxxxdxdydzdtx)v(dmdmMxxxx同理：dxdydzdty)v(Myydxdydzdtz)v(Mzzdtdt时间内，控制体总净流出质量：时间内，控制体总净流出质量：zyxMMMMdxdydzdt)v(divdxdydzdtv由质量守恒：控制体总净流出质量，必等于控制体内由于密度变化而减少的质量，即dxdydzdttdxdydzdt)v(divdxdydzdtz)v(y)v(x)v(zyx连续性微分方程连续性微分方程连续性方程的微分形式连续性方程的微分形式不可

3、压缩流体不可压缩流体0vdiv常数0zvyvxvzyx0vt)v(divt定常流动定常流动即即二、 流体微团运动分解 流体微团：指大量流体质点组成的具有线流体微团：指大量流体质点组成的具有线性尺度效应的微小流体团。性尺度效应的微小流体团。 流体在运动过程中可能发生变形或旋转，流体在运动过程中可能发生变形或旋转，只要微团的运动分析清楚了，流场的运动只要微团的运动分析清楚了，流场的运动就知道了。就知道了。zxyMM0dxdydz tttM0M流体微团的体积在单位时间的相对变化。流体微团的体积在单位时间的相对变化。平移速度平移速度v vx x,v,vy y 代表微团平移运动。代表微团平移运动。 xv

4、xxxyvyyyzvzzzv1zvyvxvtd)V(dVzyxxyyvdtdx1xvdtdy2)(21xvyvyxzvkji21zyxyvdtdx1xvdtdy2xyyxxvyvdtdd)(21)(2121)xvxv(ijjiij2121zzyyxx,zxyzxy,tttijijv2（有旋）0v（无旋）0vTaylorTaylor展开并略去高阶小量，有展开并略去高阶小量，有 zzyyxxtzyxtzzyyxxMvvvvv),(),(zzvyyvxxvvvzzvyyvxxvvvzzvyyvxxvvvzzzzMzyyyyMyxxxxMx3, 2, 1ixxvvvjjiiiMt t 时刻：流体微团时

5、刻：流体微团 ),(0zyxMkjivwvutzyx),(kjirzyx),(zzyyxxMkjivMMMMwvutzzyyxx),()(21)(21ijjiijjiijijjixvxvxvxvxv)(21ijjiijxvxv)(21ijjiijxvxvjjiiiMxxvvv变形运动变形运动旋转运动旋转运动平移运动平移运动亥姆霍兹运动分解定理三、理想流体运动方程 牛顿第二定律maF 流体平衡的欧拉方程：01pf 欧拉方法中加速度的表达形式：zvvyvvxvvtvaxzxyxxxxvvtva流体运动欧拉方程zpfzvvyvvxvvtvypfzvvyvvxvvtvxpfzvvyvvxvvtvzzz

6、zyzxzyyzyyyxyxxzxyxxx111兰姆方程兰姆方程pfvwvtv1222当当 流动是无旋的；流动是无旋的；否则，流动是有旋的。否则，流动是有旋的。0zyxwww四、定解条件 1、起始条件：起始瞬时流场中的流动分布、起始条件：起始瞬时流场中的流动分布z , y, x,z , y, xpp,z , y, xvv ,t 0是研究非定常流动必不可少的定解条件是研究非定常流动必不可少的定解条件 2、边界条件：方程组的解在流场边界上应当满足的条件。、边界条件：方程组的解在流场边界上应当满足的条件。A、固体壁面：壁面上流体质点的法向速度应等于对应点上壁面的法向速度、固体壁面：壁面上流体质点的法

7、向速度应等于对应点上壁面的法向速度 流体与固体壁面的作用力也必沿壁面法线方向；流体与固体壁面的作用力也必沿壁面法线方向；B、流体交界面：在交界面同一点，两种流体法向速度相等，对于平面，压力相、流体交界面：在交界面同一点，两种流体法向速度相等，对于平面，压力相等；等；C、无穷远处：一般给定参数；、无穷远处：一般给定参数；D、流道进、出口处，可根据具体情况确定。、流道进、出口处，可根据具体情况确定。0v积分时时间变量积分时时间变量t 作常数处理。作常数处理。 任一时刻，涡线上每一点的任一时刻，涡线上每一点的切向量都与该点的涡向量相切。涡线微分方程切向量都与该点的涡向量相切。涡线微分方程0 r d)

8、,(),(),(tzyxdztzyxdytzyxdxzyxrd 旋转角速度的值与垂旋转角速度的值与垂直于角速度方向的微源涡管横截面积的乘积直于角速度方向的微源涡管横截面积的乘积的两倍。的两倍。涡量场的通量（涡强）。涡量场的通量（涡强）。AndAdAJ2dzvdyvdxvsdzyxKKvnCS速度环量是标量，其正负号与速度和线积分绕行方向有关，速度环量是标量，其正负号与速度和线积分绕行方向有关，规定：其绕行正方向为逆时针方向，面积的法线与正方向规定：其绕行正方向为逆时针方向，面积的法线与正方向形成右手螺旋系统。形成右手螺旋系统。ABCDxvdyyvvxxdxyvvyydyyvdxxvvyyydy

9、yvdxxvvxxxdxxvvxxdxxvvyyyvxy0从从A点起逆时针方向积分，可以得到点起逆时针方向积分，可以得到微分形式的速度环量为微分形式的速度环量为dyvvdxvvdyvvdxvvdyAyDxDxCyCyBxBxA2222dxdyyvxvdxyzxyyvxvdxdyd2dAwsdvAn 2将各点速度代入，并忽略高阶小量，得到将各点速度代入，并忽略高阶小量，得到JKAddAs dvKvJC例题例题 已知理想流体定常流动的速度分布公式为已知理想流体定常流动的速度分布公式为 试求涡线方程与沿封闭周线试求涡线方程与沿封闭周线 的速度环量，的速度环量，a，b为常数。为常数。02122zyxv

10、v ,xyav0222zbyx例题例题 已知平面流动的流速为已知平面流动的流速为:（1）检查是否连续；（）检查是否连续；（2）是否无旋；）是否无旋； yxxvx422yxyvy22 有旋流动特点特点: :（有旋）0v1、),(),(),(tzyxdztzyxdytzyxdxzyx2、JdzvdyvdxvdAJzyxAn23、第七节 汤姆孙定理 亥姆霍兹定理1、汤姆孙定理：理想正压流体在有势质量力的作用下，速度环量和旋涡都是不能自行产生或消失。理想流体没有粘性。2、亥姆霍兹定理：A、在同一瞬间涡管各截面上的涡通量都相同。涡管不能在流体中开始或者终止，自能形自封闭管圈，或者在边界上开始或终止。B、

11、理想正压流体在有势的质量力作用下，涡管永远保持为由相同流体质点组成的涡管。C、涡管强度不随时间变化，永远保持为定值。CdAAn2涡管不可能在流体中开始或终止，它只能自成封闭形，涡管不可能在流体中开始或终止，它只能自成封闭形，或开始、终止于边界面或伸展到无穷远。或开始、终止于边界面或伸展到无穷远。龙卷风开始和终止于地面与云层。龙卷风开始和终止于地面与云层。烟圈呈环形。烟圈呈环形。第八节 平面涡流 前提：重力作用，理想不可压流体。前提：重力作用，理想不可压流体。 一无限长，涡通量为一无限长，涡通量为J的铅直涡束，象刚体的铅直涡束，象刚体一样以等角速度一样以等角速度绕自身轴旋转。涡束周围绕自身轴旋转

12、。涡束周围的流体受涡束的诱导将绕涡束轴做对应的的流体受涡束的诱导将绕涡束轴做对应的等速圆周运动。等速圆周运动。 涡束内的流动为有旋流动，称为涡核区；涡束内的流动为有旋流动，称为涡核区； 涡束外的流动为无旋流动，称为环流区涡束外的流动为无旋流动，称为环流区。第八节 平面涡流vxy环流区速度分布环流区速度分布:rvv ,vrrrvrb202伯努利方程伯努利方程:2822/222222bbvpprpvpppvp涡核区涡核区:222221, 0bbrrrpprrrvvv由欧拉方程：由欧拉方程：涡核边沿：涡核边沿：第八节 平面涡流 涡核中心区，流速为涡核中心区，流速为0， 压强为：压强为： 涡核边沿至涡

13、核中心的压降为：涡核边沿至涡核中心的压降为：2bcvppbbcbppvpp221结论：在环流区随着半径减小，流速升高，压强结论：在环流区随着半径减小，流速升高，压强降低；涡核区和环流区的压强降相等；涡核区的降低；涡核区和环流区的压强降相等；涡核区的压强比环流区的低，涡核区本身很小，使得径向压强比环流区的低，涡核区本身很小，使得径向压强梯度大，故有压强梯度大，故有向涡核中心抽吸作用向涡核中心抽吸作用。应用：离心泵，旋风燃烧室，离心式除尘器等。应用：离心泵，旋风燃烧室，离心式除尘器等。主要内容：1、速度势；2、流函数；3、简单平面势流及其叠加。第九节 速度势、流函数、流网第九节 速度势、流函数、流

14、网 一、速度势一、速度势yvxv,xvzv,zvyvwxyzxyz 0t , z , y, x为函数全微分的充要条件dzvdyvdxvzyxgradkzjyixvzv ,yv ,xvdzzdyydxxdzyx一、速度势一、速度势BAABBAzyxABddzvdyvdxv22222220zyxAB积分与曲线形状无关不可压缩流体连续性方程：0zvyvxvzyxl不可压缩流体的有势流动中，速度满足拉普拉斯方程；l势函数为调和函数；l势流为不可压缩流体的无旋流。二、流函数二、流函数 对于不可压平面流动yxyxvdyvdxyvxvdyvdxvdxv ,yvxyyxl每条流线上，每条流线上，d=0，=常数

15、，所以常数，所以（x，y）为流函数）为流函数对于无旋流动：对于无旋流动：222220yx流函数流函数为调和函数为调和函数P206，qv=2-1，沿流线全长，两流线间的流量保持不变。三、流网三、流网 对于不可压缩流体平面无旋流动中：对于不可压缩流体平面无旋流动中：xyvyxvyx0yyxx可知流线和等势线是正交的，等势线可知流线和等势线是正交的，等势线iC和流线和流线形成了由相互垂直的交叉线组成的网，形成了由相互垂直的交叉线组成的网，称为流网。称为流网。iK流线等势线1K2K3K1C2C3C4Cyvyxvx第十节 几种简单的平面势流 势函数的作用：对于不可压缩流体的无旋流动问题归结为根据起始条件

16、和边界条件求解拉普拉斯方程问题。求解：均匀等速流，源流和汇流，势涡m0V一、均匀等速流 流线平行，流速相等流线平行，流速相等jvivjvivvyxyx00其中vx0，vy0为常数yvxvdyvdxvdyydxxdyxyx0000yvxvdyvdxvdyydxxdxyxy0000P=常数常数二、源流和汇流 x y =const =const r源流 x y =const =const r汇流定义：在无限平面上，流体从一点沿径向直线均匀向各方流出， 称为源流。 若流体沿径向直线均匀从各方流入一点，称为汇流。二、源流和汇流01rv ,rvrCqrvvr12rqvvr2rdrqdrvdrrdvr2rl

17、nqv201rv ,rvr22vvrqdqrdvd伯努利方程伯努利方程:2222182rqpppvpvr压强随半径减小而降低压强随半径减小而降低三、势涡v，若逆时针，上式加负号。，若逆时针，上式加负号。2rln2 rrvrvr21, 0222282rpvpp第十一节 简单平面势流的叠加nnnnnnnvvvv212121221222212221222212)(0)(0)(结论：调和函数可以叠加，叠加成的新函数仍然是调结论：调和函数可以叠加，叠加成的新函数仍然是调和函数，流动仍然无旋；叠加成新平面势流的速度矢和函数，流动仍然无旋；叠加成新平面势流的速度矢等于原先平面势流速度的矢量和。等于原先平面势

18、流速度的矢量和。（1）汇流与势涡叠加）汇流与势涡叠加螺旋流螺旋流第十一节 简单平面势流的叠加2222821,2ln21ln21rqpprrvrqrvrqrqVVrVV流体自外沿圆周切向进入，又从中央不断流流体自外沿圆周切向进入，又从中央不断流出。离心泵、离心风机等外壳中的流动。出。离心泵、离心风机等外壳中的流动。（2）源流与汇流的叠加）源流与汇流的叠加 偶极子流偶极子流 流量均为流量均为qV，第十一节 简单平面势流的叠加aa2a222122222tan2)(2)()(ln4ln2)ln(ln2ayxayqqyaxyaxqrrqrrqVBAVVBAVBAV当源点和汇点无限接近时，流量同当源点和汇

19、点无限接近时，流量同时为无穷大，以使时为无穷大，以使 ，称为偶极子流。称为偶极子流。M为偶极子矩，也称偶极子流强度。为偶极子矩，也称偶极子流强度。MaqVqaV2lim0等势线方程（等势线方程（8-61），流线方程（），流线方程（8-62）22sin2;cos2sin2;cos2rMvrMvrMrMr第十二节第十二节 均匀等速绕过圆柱体的均匀等速绕过圆柱体的平面流动平面流动 沿沿x轴正向速度为轴正向速度为 v 的均匀等速流与沿的均匀等速流与沿x轴正向轴正向强度为强度为M的偶极子流的叠加的偶极子流的叠加。xrcosMcosrV2rsinMsinrV2流线方程：流线方程：零流线零流线2100202

20、20202vMr ,sinrrMv,CCsinrrMv=0，流体不能穿过流线，零流线的圆代表圆柱的横截面。流体不能穿过流线，零流线的圆代表圆柱的横截面。零流线以原点为圆心，以零流线以原点为圆心，以r0为半径的圆为半径的圆第十二节 均匀等速绕过圆柱体的平面流动)1 (sin1)1 (cos220220rrVrvrrVrvrsinvv ,cosvv ,rr在无穷远处，为均匀等速流：在无穷远处，为均匀等速流：在零流线上：在零流线上：sinvv ,vr20流体不能穿过柱体，圆柱绕流时，流体不能穿过柱体，圆柱绕流时，vrMrr2002,在在=180。，和，和=0。 ，速度为，速度为0，称为驻点，在，称为

21、驻点，在=+90。 速度达到速度达到最大值最大值第十二节 均匀等速绕过圆柱体的平面流动无穷远处无穷远处V，p ）)sin41 (222Vpp2221sin41Vppcp2222VpVp 得全流场压力分布。得全流场压力分布。柱面上（柱面上（r =r0）：）： 压强系数压强系数圆柱的前后驻点压强达到最大值，圆柱的上下顶点压强达到最小圆柱的前后驻点压强达到最大值，圆柱的上下顶点压强达到最小值。压强对称分布，流体作用在圆柱面上总压力为值。压强对称分布，流体作用在圆柱面上总压力为0。 阻力：总压力平行于来流方向的分力。阻力：总压力平行于来流方向的分力。FD 升力：总压力垂直来流方向的分力。升力：总压力垂

22、直来流方向的分力。FL第十二节 均匀等速绕过圆柱体的平面流动第十二节 均匀等速绕过圆柱体的平面流动 圆柱体在理想流体中作等速直线运圆柱体在理想流体中作等速直线运动时，受到流体作用的阻力等于零，动时，受到流体作用的阻力等于零，没有考虑流体的粘性没有考虑流体的粘性。 01d)sin(cosd200rpSpSjinP00LDF;FNo liftNo drag 圆柱受力圆柱受力第十三节 绕圆柱有环流的平面运动点涡偶级均匀流圆柱r)rr(sinVrv)rr(cosVrvr21112202202020)rrr(cosVrln)rrr(sinV220Voyxrr00220rsinVvvrr=r0流线是圆周，

23、可以流线是圆周，可以代替圆柱面代替圆柱面vvvvvvrr,sin,cos,满足无穷远处条件满足无穷远处条件第十三节 绕圆柱有环流的平面运动若若 1，柱面上有两个驻点，柱面上有两个驻点: 和和 ；若若 =1，柱面上只有一个驻点，柱面上只有一个驻点: ；若若 1，柱面上无驻点，柱面上无驻点: 。04rv22,01 rr12环量对流场的影响：环量对流场的影响：04rv04rv8-33（a）8-33（b）8-33（c）第十三节 绕圆柱有环流的平面运动 柱面上（柱面上（r=r0）：）： )2sin2(1 21202rvvpp2)2sin2(1aVcpjinnvdSrvvpdSpFSS0)2sin2(1

24、2122lVFL000DFV LF阻力：阻力：升力：升力：等于密度等于密度 、流速、流速V、环量、环量、和柱体长度的乘积。、和柱体长度的乘积。沿沿V方向逆速度环量旋转方向逆速度环量旋转90所对应的方向。图所对应的方向。图8-34Pressure coefficient 压力分布压力分布库塔库塔- -儒可夫斯基公式儒可夫斯基公式升力产生的原因（Magnus effect）：绕圆柱上下表面流动不对称、环量（旋转）、粘性。绕圆柱上下表面流动不对称、环量（旋转）、粘性。机翼周围流场不对称、环量（机翼几何形状、攻角）、粘性。机翼周围流场不对称、环量（机翼几何形状、攻角）、粘性。鸟的飞翔、球在空中旋转、飞机飞行、流体机械叶片的升力等鸟的飞翔、球在空中旋转、飞机飞行、流体机械叶片的升力等000第十三节 绕圆柱有环流的平面运动第十四节 叶栅的库塔翼型的翼弦b、翼展l、冲角vF第十四节 叶栅的库塔问题： 绕翼型流动的升力是如何产生的？绕翼型流动的升力是如何产生的？主要内容 微分形式的连续性方程微分形式的连续性方程 流体微团的运动分解流体微团的运动分解