第六章 势流理论本章内容:1. 势流问题求解的思路2.库塔----儒可夫斯基条件3. 势流的迭加法绕圆柱的无环绕流,绕圆柱的有环绕流4.布拉休斯公式5.库塔----儒可夫斯基定理学习这部分内容的目的有二:其一,获得解决势流问题的入门知识,即关键问题是求解速度势求出速度势之后,可按一定的步骤解出速度分布、压力分布,以及流体和固体之间的作用力其二,明确两点重要结论:1)园柱体在理想流体中作等速直线运动时,阻力为零(达朗贝尔疑题) ;升力也为零2)园柱本身转动同时作等速直线运动时,则受到升力作用(麦格鲁斯效应) 本章重点:1、 平面势流问题求解的基本思想2、 势流迭加法3、 物面条件,无穷远处条件4、 绕圆柱有环流,无环流流动的结论,即速度分布,压力分布,压力系数分布,驻点位置,流线图谱,升力,阻力,环流方向等5、 四个简单势流的速度势函数,流函数及其流线图谱6、 麦马格鲁斯效应的概念7、 计算任意形状柱体受流体作用力的卜拉修斯定理8、 附加惯性力,附加质量的概念本章难点:1.绕圆柱有环流,无环流流动的结论,即速度分布,压力分布,压力系数分布,驻点位置,流线图谱,升力,阻力,环流方向等2.任意形状柱体受流体作用力的卜拉修斯定理3.附加惯性力,附加质量的概念§6-1 几种简单的平面势流平面流动:平面上任何一点的速度、加速度都平行于所在平面,无垂直于该平面的分量;与该平面相平行的所有其它平面上的流动情况完全一样。
例如:1)绕一个无穷长机翼的流动, 2)船舶在水面上的垂直振荡问题,由于船长比宽度及吃水大得多,且船型纵向变化比较缓慢,可以近似认为流体只在垂直于船长方向的平面内流动,如图6-2所示如果我们在船长方向将船分割成许多薄片,并且假定绕各薄片的流动互不影响的话,则这一问题就可以按平面问题处理这一近似方法在船舶流体力学领域内称为切片理论一、均匀流 流体质点沿x轴平行的均匀速度 Vo,如图6-5 所示, Vx =V o, Vy=0平面流动速度势的全微分为 dxVdxdxy0积分:φ=Vox (6-4 ) 如图6-3流函数的全微分为, dydxydxo积分:ψ=Voy (6-5 ) 如图6-4由(6-4)和(6-5)可得: 流线:y=const,一组平行于x轴的直线,如图6 -3 中的实线等势线:x=const,一组平行于y轴的直线,如图6-3 中的虚线。
均匀流的速度势还可用来表示平行平壁间的流动或薄平板的均匀纵向绕流,如图6-4 所示二、源或汇平面源:流体由坐标原点出发沿射线流出,反之,流体从各个方向流过来汇聚于一点,谓之平面汇:与源的流动方向相反设源的体积流量为Q,速度以源为中心,沿矢径方向向外,沿圆周切线方向速度分量为零现以原点为中心,任一半径r作一圆,则根据不可压缩流体的连续性方程,体积流量Q 2πrvr=Q∴vr=Q/2πr (6-6)在直角坐标中,有 xyVx在极坐标中有: (6-7) 图6-6 点源和点汇rrsVr 1极坐标中 φ 和 ψ 的全微分:(6-8)2ln2QrdQrVddr rssr流线:为 θ=const,从原点引出的一组射线;等势线为r=const,就是和流线正交的一组同心圆由(6-6)式可看出,当Q>0,则vr>0,坐标原点为源点;如果Q<0,则vr<0,流体向原点汇合, 图6-7 扩大壁面和源的互换性乃是汇点源(汇)的速度势,还适用于扩大(收缩) ,渠道中理想流体的流动,如图6-7 所示三、偶极子 图6-8 偶极偶极流:流量相等的源和汇无限靠近,且随着其间距 δx→0,其流量Q→∞,且Qδx→M(δx→0) (6-9 )则这种流动的极限状态称为偶极子,M称为偶极矩。
用迭加法求 φ 和 ψ由图6-8 ( a)所示:)ln(2211rQ2cosxr因此 )cos1ln(2ln2)l(n222 111rxQrrr式中z=δxcosθ1r2是个小量,我们利用泰劳展开式将 φ 展开并略去 δx二阶以上小量得32)1ln(zz当 δx→0时,Qδx→M,θ1→θ,r2→r其中r,θ 为A点的极坐标,这样便可从上式得到偶极子的速度势为 (6-10)直角坐标有(6-11)2yxM对于流函数: )()(2121 Q图6-8(a)三角形BCD:r2δθ=δxsinθ1,有 21sinrx所以 2sinrxMn θ r 2 当 δx→0时,Qδx→M,r 2→r,θ 1→θ,所以(6-12)rsin直角坐标有 (6-13)2yx令 ψ=C 即得流线族: cM2或 12yx即 012cyx配方后得 (6-14)21124)(流线:圆心在y轴上与 x 轴相切的一组圆,如图6-10(b)中的实线。
流体是沿着上述的圆周,由坐标原点流出,重新又流入原点等势线:中心在x轴上与y轴相切的一组圆,并与 ψ=const 正交,如图6-8(b)中的虚rcos221cosrxQ线应当注意的是,偶极子是有轴线和有方向的源和汇所在的直线就是偶极子的轴线,由汇指向源的方向,就是偶极轴的方向如图6-8 所示的偶极子的方向是x轴的负向四、点涡(环流)流场中坐标原点处有一根无穷长直涡索,方向垂直于平面xy平面,与xy平面的交点为一个点涡点涡在平面上的诱导速度沿着以点涡为中心的圆周的切线方向,大小与半径成反比,即(6-15)02rsvv极坐标下: ddsr2积分得: (6-16) 流函数 rrvrs积分: (6-17) ln2流线:ψ=const 就是r=C,即一组以涡点为中心的同心圆, 如图6-9 所示注意:Γ>0对应于反时针的转动,Γ<0对应于顺时针的涡旋。
§6-3 绕圆柱体的无环量流动,达朗贝尔谬理势流迭加法:均匀流、源汇、偶极子、点涡这样一些几种简单的势流,具有可迭加性将它们之中的两个或两个以上迭加起来,在用物面边界条件来控制,会获得有实际意义的结果绕圆柱体的无环流流动就是一个典型的实例理想流体的边界条件:1) 无穷远条件(远场条件)r=∞, 0yxv或r=∞, sincor2)物面条件(近场条件): r=r 0 ,v n =v r =0 称为不可穿透条件零流线: r=r0处 ψ=0是一条流线圆柱在静止无界流体中作等速直线运动 = 均匀流动+ 偶极子流动均匀流和偶极子迭加后的速度势和流函数为: ( 6-18) 1202MCosvrr(6-19) 120Sini观察 ψ=0这条流线,由(6-28)式,我们有:0)2(0rvSin若sinθ=0,有 θ=0或 π,因此 ψ=0的流线中有一部分是x轴;若v 0r -M2πr=0, 0rMv即r2=M2πv0, 令 , 就有r=r0, 即r=r0的圆周02r20rv也是 ψ=0的流线的一部分, 如图6-10 所示。
验证边界条件,将 代入 φ,有20rv(6-20))(cos0速度(6-21))1(sin1c2002rvrvrr 当r→∞,从上式可得 sico0vr当r=r0时,vr=0这样就证明了均匀流和偶极子迭加的速度势,完全满足绕圆柱体无环流流动的远场和近场的边界条件,它在r≥r0的流动情况与均匀流绕圆柱的流动情况完全一样想象把均匀流加偶极子的流动图案中r<r0的那一部分去掉(不感兴趣) ,而在其中充实以一个r=r0的圆柱体,对流场流动不会有任何影响因此,绕圆柱体无环流流动的速度势就是均匀流加偶极子的速度势圆柱表面的速度分布圆柱表面上速度分布:r=r0时:(6-22)sin20vr负号表示其方向与 s 坐标轴方向相反, 如图6-10驻点位置:A,C两点 θ=π 或0,vs=0称为驻点或分流点对B,D两点:(6-23)022vB,D两点:速度达到最大值,与圆柱体半径无关,恰等于来流速度v0之两倍流体从较远处以流速v0流向圆柱,当接近圆柱时,流速逐渐减小,到达A点时速度降至零然后分为二支向两侧流去,同时速度逐渐增大。
B,D两点:速度增至2v0,达最大值然后又逐渐减小,在C点汇合时,速度又降至零离开C点后,又逐渐加速,流向后方的无限远处时,再恢复为v0圆柱表面上压力分布:运动是定常,设无穷远均匀流中的压力为p0,忽略了质量力,拉格朗日方程 202vpv将园柱表面上速度分布代入,即得圆柱表面上压力分布 (6-24))sin41(220vp物面上的压力分布往往用下式定义的无因次系数来表示:(6-25)201vCp由(6-24)式可得 (6-26)4sinp压力系数见图6-11(a)中压力分布既对称于x轴也对称于y轴A,C两点压力最大cp=1B,D两点压力最小cp=-3 (6-27)沿 ψ=0这条流线压力变化为:左方无限远处,cp=0,当流体流向圆柱体时,压力逐渐增大,流到A点时压力为极大值cp=1由A点分为两支分别流向B,D点,压力逐渐减小,到达这两点时压力为极小值,cp=-3由B,D点流向C点时,压力逐渐增大,到达C点,恢复到极大值,cp=1由C点流向右方无限远处,压力又再次减小,最后压力重新降至p0,cp=0。
(a)理想流体;(b)真实流体图6-11 圆柱体表面上的压力系数分布理想流体对圆柱体的作用力:因为其压力分布对称于x轴,显然合力在y轴上的分力L(升力)为零;同样,因其压力分布对称于y轴,故合力在x轴上的分力R(阻力)为零,即升力L=0阻力 R=0 (6-28)这一结果与实验结果有严重矛盾,称为达朗贝尔谬理如图6-11(b)所示为圆柱表面压力分布的实测结果和图6- 11(a)相比较,其前半部保持不变,后半部发生很大变化,C点处由正压变为负压,破坏了压力分布对y轴的对称性,从而引起了作用于物体的阻力其原理,边界层理论一章再详细讨论达朗贝尔谬理在理论上仍然是很有意义达朗贝尔谬理成立的条件可归纳为下列五点;1)在理想流体中;2)物体周围的流场是无边际的;3)物体周围的流场中没有源、汇、涡等奇点存在;4)物体作等速直线运动;5)流动在物体表面上没有分离如果上述条件全部成立,那么任何物体的的确确将不会上面任一条件被破坏,则物体即将遭受到流体的作用力(阻力或升力) 因此,根据达。