6第六章 不可压缩流体的平面势流
第六章 不可压缩流体的平面势流
?6-1 有势流动的速度势函数
,一、速度势函数
对于无旋流动,有
,,w,,,,,y,z,
,u,w,,,,z,x,
,,,u,,,,x,y, (1)
udx,,dy,wdz,(x,y,z,t)根据数学分析可知:上式成立是成为某一函数的
,全微分的充要条件。称为速度势函数,简称速度势。
d,,udx,,dy,wdz即:
,,,,,,d,dx,dy,dz,,x,y,z又有:
,,,,,,,,w,?u,,y,x,z,,
又由矢量分析:
,,,,,,,,,,,,,V,ui,i,wk,i,i,k,,grad,,,,,x,y,z (2) 即速度势的梯度等于流场的速度。
,,,,r,r在柱坐标中:径向速度:
,,,,,,,,,sr,,切向速度:
,,,,z,z轴向速度:
,,,V由此可见,对任意方向的偏导数,就是速度在该方向的投影,这是的
一个重要性质。
,(rotV,0),(x,y,z,t)函数称为速度势函数,简称速度势,对无旋流动,总
有速度势存在,所以,无旋流动也称为有势流动。
,Γ在有势流动中,和的关系为:
,BBB,,,Γ,V,ds,udx,,dy,wdz,d,,,,,BAAB,,,AAA (3)
即在有势流动中,沿AB曲线的切向速度线积分(速度环量)等于终点B与起点A的速度势之差。
又:在有势流动中,沿任一封闭周线K的速度环量
,,,,Γ,V,ds,udx,,dy,wdz,d,,,,KKK
d,,0,,K若是单值或由斯托克斯定理,则
二、势函数方程
,,,,,,,,w,u,,y,x,z将,,代入不可压流体连续方程:
,,u,,w,,,0,x,y,z
222,,,,,,2,,,,,,0222xyz,,,则有: (4)
222,,,2,,,,,,222,x,y,z(其中称为拉普拉斯算子)
,即在不可压流体的有势流动中,速度势满足拉普拉斯方程。凡是满足拉普拉斯方程的函数,数学上称为调和函数,所以,速度势点数是一个调和函数。
,对柱面坐标,的拉普拉斯方程为:
222,,1,,1,,,,2,,,,,,222r,rr,r,,z, (5)
,,,,,,,,,,,,zr,,rr,,,z〔推导过程为:将,,代入柱面坐标的连续方程,即可〕 根据以上讨论可知:只要流体流动无旋。则必然存在单值的速度势函数,
,反之,若流场中存在单值的势函数,则此流动必为无旋流动。此外,流动无旋,流场中沿封闭曲线的速度环量为零。
,,,0,,即:
,,,0,Γ,0
因此,求解不可压流体的无旋流动问
题
快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题
,便可归结为求解速度势问题,以
2pVz,,,constg2g,此求得速度场,再由无旋流动的伯努利方程,求得压力分布。
?6-2 流函数
,一、不可压缩流体的流函数
以上引进的势函数虽然能使问题简化,但它仅限于有势流动(即无旋流动),
对于有旋流动,我们必须根据定常,不可压二元流动的连续方程,引出流函数的
概念。
推导如下:由二元不可压流体的连续方程
,u,,,,0,x,y
,,u,,,,x,y则: (1) 又:平面流动的流线微分方程为:
dxdy,u,
,udy,,dx,0 (2)
,(x,y)由数学分析可知:(1)是(2)成为某一函数的全微分的充要条件。即:
d,,udy,,dx (3)
,,,,d,dx,dy,,x,y又:
,,,,u,,,,,y,x所以, , (4)
22,,,,u,,,,,,,x,y,x,y,x,y则, ,
22,,,,u,,,,,,,0,x,y,y,x,x,y
,说明满足连续性方程
二、流函数的基本性质
(1) 等流函数线为流线
udy,,dx,d,,0,显然,在流线上,,即
,,c,,c即: 即,的曲线为流线。
,在每条流线上的常数值各不相同。
(2) 即:平面流动中两条流线间通过的流体流量等于两条流线上的流函数之差。
,,,12ψ 证:取两道流成,再取曲线AB垂2 ψ ψ 1 直于各流线,假定垂直纸面的尺寸为1,在AB Y , ,υ V Vdl曲线上取微元线段,其上速度为,则通过 dl 曲线AB的体积流量为: B u ,,dy BB ,,q,Vdl,ucos(V,x),,cos(V,y)dln-dx ,, AA A -dx B,dydx,, ,u,,()dl,,,Adldl,,X O BB ,,,udy,,dx,d,,,,,,,,, BA21,,AA
dx,cos(V,y)dlX指向减小方向,为负。为使为正,所以在dx前加负号。
证毕。
,由此可见:两根流线之间的流量等于两流函数的差值。同时,由于在引出这个概念时,没有涉及流体是有粘性还是无粘性(即理想或实际),有旋或无旋。所以,不论是有粘性还是无粘性,有旋还是无旋,只要是不可压流体的平面流动,就存在流函数。
三、流函数方程
222,,,,,,2,,,,,,0222xyz,,, (5)
四、边界条件
若无穷远处均匀来流绕流一物体时,在不分离的情况下,对于固定不动的边界,在壁面上流体的法向速度为0,而壁面必然是流线,通常令壁面上的流函数值为0,因此,壁面上的边界条件可写作:
,,V,,0n,,n或=0 (6a)
对于无穷远处均匀来流,当取X轴与来流方向一致时,则有
,,,,,,,,uV,,,,xy,,,,,,,,,,,0,,,,yx, (6b)
,,与之间的关系 五、
1(满足柯西黎曼条件
,,对不可压流体的平面无旋(有势)流动,则必然同时存在和,而对平面无
,,0z旋流动,由,可推出
,,,u,,0,x,y则 (7)
,,,,u,,,,,y,x再将 ,代入上式
22,,,,2,,,,0,22,x,y得: (8)
22,,,,1,1,2,,,,,,02222rrrr,,,,对极坐标: (9)
,所以不可压流体平面无旋流动的流函数,满足拉普拉斯方程,也是调和函数。
又:对平面无旋流动,必然存在
,,,,,u,,,,x,y,,,,,,,,,,,,,y,x,,由 柯西黎曼条件
2(流线与等势线正交
,,,,,,0是流线与等势线正交的条件
,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,i,j,i,j,,,,,0,,,,,x,y,x,y,x,x,y,y,,,,
,,c,,c,,c式就是等势线簇和流线簇互相正交的条件。所以说明(等势
,,c线簇) 和流线簇正交。
,,c,,c在XY坐标平面上由及画图,构成正交网格,称为流网。如下图所示。
?6-3 几种简单的势流流动
,,不可压流体平面无旋流动的流函数和势函数,满足拉普拉斯方程,而且,拉普拉斯方程是线性齐次方程,其解具有叠加性。
,,,12设是两个有势流动,均满足:
2222,,,,,,,,222211,,,,,0,,,,,0122222xyxy,,,,,
叠加后,可得;
22,,,,,,,,,,,,21212,,,,,,,,,01222xy,, (1)
同理:
22,,,,,,,,,,,,21212,,,,,,,,,01222xy,, (2)
,,
V,,,V,,,1122,
,,,
,,V,,,,,,,,,,,,,,,V,V121212
一、 一、
y
平行流(均匀直线流动,无旋,参看下图)
0 x V υ ,
θ u
设流体作等速直线流动,流场中各点速度大小,方向均相同。即
u,Vcos,,,Vsin,,,,
,,xVcos,,yVsin,,, (3a)
d,,,,dx,udy
,,,xVsin,,yVcos,,, (3b)
,,0,,0当取x轴与来流方向一致时,则有,,
,,xV,,yV,,,
,,c,,c显然,与互相垂直(斜率互为负倒数),并且都满足拉普拉斯方程。 又:由位流伯努利方程
2pVz,,,const,g2g
pz,,const,gV,const由,则
p,const若平行流在水平面上进行(即z=常数),或流体重度可忽略不计,则 即流场中压力处处相等。
二、点源与点汇
设在无限大平面上流体从一点沿径向直线均匀地向各方流出。(如图所示)这种流动称为点源。这个点称为源点。
反之,若流体沿径向直线匀地从各方流入一点,这种流动称为点汇,这个点称为汇点。
y y
x x
,r显然,点源和点汇流动都只有径向速度 点汇 点源 ,, ,,r,r且
,,0,切向速度
又:对半径为r,单位长度的圆柱面,由质量守恒,则流体通过统一圆柱面的流量Q应相等。
Q,,,rQ,2,r,1,,,常数2,rr则: (4)
式中,Q是点源(或点汇)单位时间流入(或流出)的流量。称为点源或点汇的
,,0,,0,,rr强度。对点源Q,0,我们取+Q,对点汇,Q,0,,取-Q。
,,Q,,d,dr,dr,,drr,,r2r则
QQ22,,,,lnr,,lnx,y,,22 (5)
,,r当r=0时,与都变成无穷大,所以源点(或汇点)是奇点。所以(4),(5)仅仅在源点(汇点)以外才适用。
又由柯西黎曼条件:
,,,,,,,,,,,,x,y,y,x
,,,,d,dx,dy,,x,y 则
Q,,,,2,积分: (6)
或由极坐标的柯西黎曼条件:
,,,,,,,,,,,,r,s,s,r
QQ,,,,,drdsrd,,,,,,,,,,rs2r2,,,,则:
,,constr,const等势线,即是半径不同得同心圆,(由
2,cQQ,,lnr,C,r,e,,,const2,,,const圆的方程),与流线(即)互相正交。(参图虚线为等势线,实线为流线)。
,,又将与代入极坐标的拉普拉斯算子
22,1,1,2,,,,222,rr,rr,,
22,,,0,,,0均满足拉普拉斯方程 ,
由此可见,点源与点汇均是无旋流动。
Q,,,0r,,r2对无限水平面的无穷远处
三、涡流与点涡
1、速度分布
设有一旋涡强度为J的直线涡束,该涡束半径为r沿Z轴方向为无限长(如0
,图),且该涡束好像刚体一样以等角度绕自身轴旋转,由于假设直线涡束沿Z轴方向无限长,即认为在与Z轴垂直的所有平面上流动情况都一样。所以,此种流动可视为平面运动处理。而涡束周围的流体将被带动着做旋转运动。如图所示,这种运动称为涡流。
设涡束轴为Z轴,则由涡束所诱导的环流的流线就是以坐标原点为圆心的
圆心园。
,由于涡束以等角速旋转,则速度沿半径的变化规律由斯托克斯定理求得: 由该定理:沿任何圆周流线的速度环量等于涡束的旋涡强度。
2,Γ,,ds,,rd,,,,,0即:
2,,r2,,J,const,2,r,,0
Γ,,,,2,r (7)
r,,,,,,0,r,,,,, 显然,由等环量,则在涡束外,
,,r,,,,涡束内,流体如同刚体一样以等角速绕自身轴旋转,即显然常r,,,,r,,,,,,数,则,
通常称涡束外为势流旋转区,涡束内为涡核区。
证:除原点外,涡流为无旋流动
Γ,,,2,r由(1)
ΓyΓy,,,,,u,,sin,,22,2,x,y2,rr而
ΓxΓx,,,,cos,,,,222,rr2,x,y
2222Γ2(x,y),2(x,y),,1,u,,0,(,),z2224,(x,y),,2xy则 =,
因此,除原点外,涡流为无旋流动证毕。
Y 2、涡核区的压力分布
推导如下:
由平面定常流动的欧拉微分方程 υ υ θ ,u,u1,p ,u,,,u x,y ,x,y,,x r ,,,,1,p ,u,,,r 0 ,x,y,,g
涡核内任一点的速度为:
,u,,y,,,,r,,,,,,x, 则
,,u,,,u,,0,,,,0,,,y,y,x,x,;,
代入欧拉方程,则:
1,p,1,p2,,,0(-),x,dxx,dx,,,,x,x,,,1,p1,p2,,,()0,y,,,,dyydy,,,,y,,y,
u,,,sin,,,r,sin,,,y,, 其中
由于u与x轴反向所以加负号
,,,cos,,r,cos,,x,,
两式则:
2,11,p,p222,,d(x,y),dp(xdx,ydy),(dx,dy)2,,,x,y 积分得:
22,,p21,2,r,C,,C,p,,,,C,22,2
r,r,p,p,,,,00,0边界条件:
1122p,,,,C,C,p,,,000022则 代入通解
则涡核区任一点压力涡束边界压力
1122p,,,,p,,,,00,022 ,为强制旋涡(即涡束)的边界速度 由于涡束外是势流,由伯努利方程:则
1122p,p,,,,p,p,,,,,000022
代入上式,得涡核区的压力分布为:
122p,p,,,,,,,,02
r,r,,,p,p0,00又:当,时,
12p,p,,,,002则
,,0,而当 r=0时 p=p=p,代入涡核区任一点压力涡束边界压力 中心C
12p,p,,,c002则
12?p,p,p,p,,,,c0002
由此可见,涡核内外压力降相等,都等于以涡核边缘的速度计算的动压头。
涡核内外的压力分布如图所示。由上图可见,由于涡核区的压力比涡核外势流旋转区时压力低,所以涡流有抽吸作用。如自然界中的龙卷风,旋风等。就是
,
Vr() r0
rr
rr pp,,
p(r)因为旋转中心的压力低,产生真空吸力,都可以把碎纸,树叶,甚至房屋等卷入(吸入)空中。另外,流体机械中,离心泵,离心式风机外壳中流体的流动也是涡流的实例。
r,00又:当涡束半径时,涡束成为一条涡线,这样的涡流称为点涡。
Γ,,,,,,r,0,02,r当时,,则,所以涡点是奇点。
,,3、点涡的及
下面求点涡的速度势及流点数
,1,Γ,,,,,,,0,r,r,,2,r,r 由于 ,
d,,,dr,,rd,r,
,,ΓΓΓy,,,,,d,d,,arctg,,,,,,,222x 再由柯西黎曼条件:
,,,,,,,,,,,,y,x,x,y
Γ,,,lnr,2则
Γ,0Γ,0对逆时针,;顺时针,。
四、点源,点汇?偶极流
如图所示:位于A点(-a,0)的点源和位于B点(a,0)的点汇迭加后,其速度势
y
P(x,y)
r A r Br
为:
QQAB,lnr,lnrAB,,,,,,,22AB
22r,y,(x,a)A其中:
22r,y,(x,a)B
Q,Q,QAB若,则
2222y,(x,a)QQy,(x,a)rQQA,ln,ln,,,,(lnrlnr)ln22AB22,,24y,(x,a)y,(x,a)22r,,B
同理,流函数为:
QQ(),,,,,,,,ABP22,,
由几何关系:一三角形外角等于两不相邻内角之和,则
,,,,,,,,,,,,BPAPAB
,P其中为动点P与AB两点连接线之间的夹角。
,,const,,,constP由流线方程,这就是说,流线是经过源点A与汇点B的圆线簇。如图所示
若点源、点汇无限接近,即a?0时,假定a?0时,Q??,并且使lim2aQ,Ma,0Q,,则这种流动称为偶极流?无旋流动,M为一个有限的常数值,称为偶极矩。
下面推导偶极流的速度势及流函数
如图所示,由(6)则:
rrra,2cos,QQQAABA,ln,ln(1,),ln(1,),,rrr222,,,BBB
r,r,r,,,,,Q,,2aQ,MABAB且当2a?0时,,,,,再由
23,,ln(1,),,,?,,23级数展开式,。当ε为无穷小时,略去高阶项,则ln(1,,),,
,,,2acosQ,A,limln(1,),,,a,202rBQ,,,,
,2acosQA,lim,a,202rQ,,B
,,McosMrcos,,2,,2r2r
Mx,222,x,y那么: (8)
同理:
Q,,(,,,)AB2,
再由:2a?0,Q??,2aQ?M
MyMy,,,,,2222,2,xyr,得: (9)
再令上式等于常数C,并改写方程,得流线方程 1
MM222x,(y,),(),,4C4C11
MM,4,C4,C11此为圆的方程:即流线是半径为,圆心在(0,)且与X轴在原
点相切的圆周簇。如图实线所示。
同理:令(9)等于常数C,得 2
等势线方程:
MM222(x,),y,(),,4C4C22
单独的偶极流没有什么实际意义,但这种复合流动对讨论下一节绕圆柱体的
流动是有用的。
流线
等势线
?6,4 平面势流的叠加流动
一、点汇+点涡?螺旋流
QQ,22,rxy,,ln,,ln,1,,,,22,Q,,,,,1,2,,对点汇
ΓΓy,,,,,arctg2,,,,22x,Γ,,,,lnr2,2,,对点涡
迭加,则:
1,,,,,,,(Qlnr,Γ,)122,
(1)
1,,,,,,,(Q,,Γlnr)122, (2) 令(1),(2)两式等于常数,则得到
Γ,Qr,Ce1等势线 (3)
Q,,Γr,Ce2流 线 (4) 式中C、C是两个常数,这是两组相互正交的对数螺旋公式,称为螺旋流。 12
,,Γ,,,,r,,2πr则 切向速度:
,,Q,,,,r,r2,r 径向速度:
22ΓQ,222,,V,,,r,224,r速度代入伯努利方程
2pV(,z,,const,其中,z,const),g2g
2222Γ,QΓ,Q
2222p4,rp4,r1122,,,,g2gg2g,,则得流场中压力分布:
r1122()()p,p,Γ,Q,122228,grr12 (5)
在锅炉设备中,离心式除尘器,离心式喷油嘴等均是点汇,点涡的实例。而
水泵、风机外壳中的流动则是点源,点涡的实例。
二、绕圆柱体无环量的流动:,均匀直线流动+偶极流
V,
在
工程
路基工程安全技术交底工程项目施工成本控制工程量增项单年度零星工程技术标正投影法基本原理
实践中常常遇到流体横向绕过圆柱体的流动,(如图所示),各种冷却及加热设备中都来用这种流动方式。(如锅炉设备中的过热器省煤器,汽轮机设备中的冷凝器等)。因此,分析及掌握绕圆柱体的流动规律是很重要的。
V,设在无穷远处有一速度为的平行流,方向与圆柱体垂直,并绕过一半径为r的无限长圆柱体。由图可见,流体在流近柱体前,不受任何干扰,流线是一0
组平行直线,流近柱体时,由于柱体的阻碍作用,流线发生弯曲,在柱体后一定
y,,距离后又恢复为一组平行直线,显然,在处,流体没有受到干扰,流线仍是平行直线。
当流体沿X轴流向圆柱体时,在A点处与柱体发生碰撞,显然,A点是驻点,该处流体速度为0,然后折转方向,贴圆柱体表面流动,因此,半径为r的0圆就是一条流线。在B点处流体要急转成水平方向,所以B点为后驻点。然后,流体又沿X轴方向流动。所以,紧贴圆柱的那一条流线以半径为r的圆周以及0与X轴所组成。
再分析稍离圆柱体的流线形状,发现在Y轴附近,有相当长的一段区域,与偶极流的流线极为相似。这说明理想流体绕流圆柱体的流动状况与偶极流的流动有着内在联系。
经前人研究认为:理想流体绕流圆柱体的流动,正是以平行于X轴的均匀直线流动与以坐标原点为偶极点的偶极流两者叠加的结果。事实上,这种设想也符合绕流圆柱体的流动规律。
则这种组合流动的流函数
,,M1M1,,,V1,rsin,,,,,,,,Vrsin,,rsin,2,,12,22Vr,,,,2,r
相应的速度场为为:
,,M1,,,V,,(1,)cos,r2,,,rV2r,,,,,,,M1,,,,,,,V1,sin,,,2,,,,r2Vr,,,,,,
令,,0,,,0r,可得两个驻点坐标,即
,,0,r,M/2,V;,,,,r,M/2,V,,。
r,,00r在圆柱体表面上有,因此可以求得圆柱体物面方程为一半径为的圆周,
Mr,02,V,。
因此,平行流绕圆柱体无环量流动的流函数可以写成:
22rr00,Vy,,(1),Vr,sin,(1),22,2r,rx,y0r (6)式中 同理势函数可以写成:
2Mxr0,,Vx,,Vr,cos,(1),22,2r,r2,x,y0r (7)式中
r,r0因为,在圆柱体内,无实际意义。
对式(7)求导,则流场中任一点的速度分量为
2,r0,,V,,(1)cos,,r2r,,2,,r0,,,V,,,,1sin,,,2,,,r,,,(8)
r,r0上式中。
V,iVr,,,由上式可知,当时,
沿包围圆柱体圆形周线的速度环量为
错误~不能通过编辑域代码创建对象。错误~不能通过编辑域代码创
建对象。
说明平行流绕圆柱体的平面流动没有速度环量。
r,r0当时,即在圆柱面上
,,0,,,,2Vsin,,,,,, (9)
,,,,0,,,0,,180在(B点)或(A点)处,(后、前驻点)
,,,,,,,,2V,,max,2 时,,即等于无穷远处来流速度的两倍。 圆柱面上任一点的压力,由伯努利方程
22,pVp,,,,,,,g2g,g2g 得:
2222222V,,pp4VsinppV(14sin),,,,,,,,,,,,,,,,,,,g2ggg2gg2g,,,,
,22,p,p,V(1,4sin,),,2 (10)
工程上常用无因次压力系数来表示流体作用在柱体上任一点处的压力。
无因次压力系数定义为:
p,p,C,p12,V,2 (11)
将(10)代入(11)得:
2C,1,4sin,P (12)
Vp,,由此可见,对绕流圆柱体,C既与r无关,也与和无关,这就是利用P0
C的方便所在。 P
根据式(12)算出的理论无因次压力系数曲线如图中实线所示。其中纵坐标为
C,横坐标为θ,θ角是从前驻点A沿顺时针方向增加。由图可见: P
12p,p,,V,,A,,,02,速度等于零,C,1,压力达到最大值 P
32P,P,,V,,00,2V,,902,,速度最大,C,-3 压力达最小值 ,P
,,,180,后驻点 C,1 P
,,,,180,,,3600,,,180与完全对称,因此,流体作用在圆柱面上的压力全部对称,合力为零。下面进行
证明
住所证明下载场所使用证明下载诊断证明下载住所证明下载爱问住所证明下载爱问
:
如图所示:
dF d, , V ,r o
ds,rd,0 垂直纸面方向取单位长圆柱,在其上取微元弧段,作用在上面的微
dF,prd,,1小是压力,
,,,,dFprdcos,x0,,,dFprd,sin,y0,则:
sin,cos,式中的负号是考虑到,当,为正值时,dF,dF的方向分别与xxy轴y轴相反。
122p,p,,V(1,4sin,),,2将代入上式积分
2,1,,22F,F,,rp,,V(1,4sin,)cos,d,,0x,,D0,,,02,,则:阻力:
2,1,,22F,F,,rp,,V(1,4sin,)sin,d,,0,,Ly0,,,02,, 升力:
上式说明:运动流体作用在圆柱面上的x方向,y方向的压力合力均为零。那么根据运动的相对性原理,亦即圆柱体在流体中作匀速运动时,将不受阻力。这个结论可推广到其它任意形状的物体。显然,这个结论与实验观察完全违背。这就是著名的达朗伯疑题。
其发生的原因是对整个流场作了位流(无旋流动)的假设,直到普朗特提出了边界层的概念,才解决了这个问题。因为根据边界层理论,即使对低粘流体(如空气和水),在运动流体紧靠物体表面的范围内,也不能忽略摩擦影响。在这个范围内,流动有旋,并且在柱体后半部分将发生边界层分离,形成旋涡区,使压力下降。从而破坏了圆柱面上压力分布的对称性,形成阻力。
小结:
平行流绕圆柱体的平面流动有以下几个特点:
V,V,(1) (1) 在无穷远处,速度为来流速度,即;
,,0,,,0r,(2) (2) 在圆柱体表面上有
Mr,0r2,V,0(3) (3) 圆柱体物面方程,这是一个半径为的圆,因为
,,0r该圆表面上有,也称驻圆。
,,0(4) (4) 该流动的速度环量。
?6,5 绕圆柱体有环量的流动
库塔,儒可夫斯基公式
平行流绕圆柱体有环流的平面流动相当于平行流绕圆柱体的无环流平面流
,,,0动+纯环流。而点涡外的速度场就是纯环流,而纯环流除原点外,,所以,纯环流是无旋流动,又称为势涡。则组合流动的速度势上数及流函数为:
2rΓ0,,,V(1,)rsin,lnr,2,r2 (1)
2rΓ0(1)cos,,,,V,r,,22r, (2) 其中我们讨论的纯环流是顺时针方向的,即Γ,0,故Γ前加了负号。则:
2,,r,0,,,,V(1,)cos,,,2,,rrr,r要求,02,r,Γ,0,,,V(1,)sin,,,,,2,r,2πrr,, (3)
r,r0因为当,在圆柱体内无实际意议
,,,,无环量,有环量 从上式可以看到:。
即径向分速不因有环量而变。
r,r0 在圆柱表面上,即时,由(3)
,,0,,(4) ,Γ, ,,,,2Vsin,,,,2,r (5)0,
下面求驻点位置:
,,0,再由(4)令(驻点)
Γsin,,,,0πrV40,则: 其中Γ前已加了负号,所以Γ为正值(5)
说明θ在三、四象限,即驻点位置在三、四象限。θ为驻点位置角。
Γ,sin,,4,rV0,将(5)变成
,sin,,sin(,,),sin(,,,)又:由三角函数:
则:
Γ,1Γ,4,rV,rV40,0,? 当时,即,圆柱表面上有两个驻点,位置角分别为
Γ,,,sin,,,,,,,,,,V,const,,如上图所示,若,则随故A、B两个驻点向下移动,并互相靠拢。
Γ,1Γ,4,rV,rV40,0,sin,,,1?当Γ?且为,即时,,即
,,,,(sin(,,),sin(,,,),1)2,则,两驻点合二为一。如图所示
Γ,1,sin,,1Γ,4,rV,rV40,0,?当Γ继续?,且,即时,则这时圆柱面,上不存在驻点,驻点沿y轴向下移到相应位置。如图所示。
驻点半径r求法如下:
令(3)式:
(a)
2 ,,r,0,,V,,,,(1)cos0,,,2 (b),rr,,2,rΓ,1,0,,,,V,,,(1)0,,2,r,,,rr2,
3,,,,或,,,22由(a)式
2r0,Γ,,2,rV(1,)sin,,2r由(b)式,因Γ前面已考虑了符号,故Γ恒为“+”,
3,,,sin,2故只有为“-”,才能满足Γ为“+”,所以θ在三、四象限,受(a)约束,故代入(b)式,求得驻点半径:
1
2,,Γ1Γ22r,,(),4r0,,,4V22πV,,,,
讨论:
通过上面的分析,可看出:驻点的位置依赖于Γ的数值。 ?Γ,0时,即圆柱无环量绕流时,驻点两个。均在x轴与圆柱表面的两个
0,,,,,,,交点上。
Γ,4,rV0,?Γ继续增加,为,两个驻点对称地下移。在三、四象限;
,,,,Γ,4,rV0,sin,,,12?Γ继续增加,驻点继续下移,当时,,,两驻点合二为一,但仍在圆柱表面;
Γ,4,rV0,?Γ继续增加,为时,驻点下移脱离圆柱表面,沿y轴进入流体之内。
?Γ继续增加,驻点位置沿y方向?,还可以从流体力学直观加以解释。由于Γ为顺时针方向的点涡运动产生的,故与圆柱体的上半部分速度方向相同,复合后流速增加,而下半部分,速度方向相反,互相抵销,故迭加后流速减小。
下面求圆柱表面上的压力公布:
22pVpV,,,,,,,g2g,g2g 由位流伯努利方程:
1122p,,V,p,,V,,22
Γ,,V2Vsin,,,,,,2,r,,00r这里,则,代入伯努利方程,得圆柱表面上的压力分布为:
,,Γ12,,ppVV,,,(2sin,),,,,,,r220,, (6)
则流体作用在单位长度圆柱表面上总压力在x、y方向的分量为:
,2F,F,,pcos,rd,xD,0阻力:
,2F,F,,psin,rd,yL0,0升力:
将(6)代入上式进行积分,得:
阻力: F,0 (7) D
F,,V,L,升力: (8)
(8)就是库塔-儒柯夫斯基升力公式。
V,F(或F)的方向:为的方向逆环量方向转90?,即为F(升力)的方向。LyL如上图所示。
升力F产生的原因:由于加入环量,上下速度场发生变化,引起上下压力L
差。产生升力。
自然界中,如鸟在空中飞翔,打乒乓球中的旋转球,削球用力方向不同,即加的环量方向不同,产生的旋转方向不同,由于摩擦效应,(考虑粘性),球周围的空气主要带动旋转,(如拉弧圈球)即升力方向就不同。球呈曲线运动,升力向下,相当于球重加大,升力向上,相当于球重减小。再如飞机的飞行,透平机械的工作原理,都可用升力理论进行解释,也是儒柯夫斯基升力定理的具体应用。
,V,,上述结论,对弧立具型也适用F,,其中F为流体作用在孤立翼型上LL
的升力F,F。 Lyy
p
d ds= rd ,,o , x r o