西安交通大学流体力学实验.docVIP

西安交通大学 流体力学实验 能动硕71 严冬 1 一、实验目的: 通过编写和调试程序，使用画图软件画出绕圆柱有环量和无环量的流线分布图， 其中来流速度、圆柱直径、环量值Γ自定。 二、实验原理及内容: 1、绕圆柱的无环量流动 绕圆柱的无环量流动可以用均匀直线流动和偶极流叠加来表示。设圆柱半径为R，均匀直线流动的速度为，沿x轴正方向;偶极中心位于原点，强度为M，偶极轴沿负x方向。 V, 上述均匀直线流动中u=，v=0,流函数为= Vrsin= Vy， 势函数=Vrcos= ,,,,V,,,, Vx。 , 上述偶极流又可等效于位于点(-a,0)处强度为+Q的点源与点(a,0)处强度为-Q的点汇 Qx，aQy的叠加(其中)。对于点源,，,，流函数u,v,2222,,2(x，a)，y2(x，a)，y QyQ22，势函数;对于点汇,,arctan(),,ln(x，a)，y，，,,22x，a QyQx,aQy，流函数，势函数,,,arctan()u,,v,,,2222,2x,a,,2(x,a)，y2(x,a)，y Q22。 ,,,ln(x,a)，y,,2 令 lim2aQ,Ma,0Q,, 22y,x2xyMMMy相加得:u= V+ v=-流函数 ,,Vy,,222222,222,2,(x，y)(x，y)2,(x，y) Mx势函数。 ,,Vy，,222,(x，y) 令=0，该流线称为零流线。得(1)=0，=所以。当=0,,,, ,或=时V=0，这两点分别称为前后滞止点;而当时，V取最大值，这两点称为舷,,,,2 点。 2VRy2,给定值，则得到一条流线的方程。分离变量得: ,=c为定值x,,,y(Vy,,), 极坐标下 ,r,r或 ,,arcsin,,arcsin,,2222,rVVR,rVVR,,,, 2 用数学软件Matlab绘制不同值时x、y的关系曲线，就能得出绕圆柱无环量流动的流 线分布图。同时还可绘制均匀直线流动和偶极流的流线、等势线分布图，比较叠加前后它们 的区别。 以下是编程程序(其中V=20m/s, r=1m, Ψ=5/10/15/20): , a=20; r=1:0.01:6; t1=r.*r*a-a; t2=1./t1; t3=5*t2; t4=asin(t3); polar(t4,r) hold on polar(pi-t4,r) hold on polar(-t4,r) hold on polar(pi+t4,r) hold on t5=10*t2; t6=asin(t5); polar(t6,r) hold on polar(pi-t6,r) hold on polar(-t6,r) hold on polar(pi+t6,r) hold on t7=15*t2; t8=asin(t7); polar(t8,r) polar(pi-t8,r) polar(-t8,r) polar(pi+t8,r) t9=20*t2; t10=asin(t9); polar(t10,r) polar(pi-t10,r) polar(-t10,r) polar(pi+t10,r) 运行结果如下: 3 90 6 12060 4 15030 2 1800 210330 240300 270 2、绕圆柱的有环量流动 如果在1的基础上再在圆心处叠加一个点涡流动，由于点涡的流线为同心圆簇，则可组 成绕圆柱有环量流动。 ,22 设点涡为顺时针方向，强度为，由定义得流函数，势函数,,,,lnx，y0,2 y, arctan,,,02,x 2VRy,22, 把1中的流动和上述点涡流动叠加起来，得，势函数,,Vy,，lnx，y,22,,2(x，y)2 2VRxy,,。给定值，则得到一条流线的方程。分离变量得: ,=C 为Vyarctan,,，,,222,x(xy)， ,2r,,,rlnr,2r,,,rlnr或 ,,arcsin,,,,arcsin22222,V(r,R)2,V(r,R),, 22 对求偏导，得速度场(1)(2)V,V(1,R/r)cos,, r, 22V,,V(1，R/r)sin,,,/2,r ,, 4 在圆柱表面上，，所以沿圆柱表面流动无分离。现在确定圆V,0V,,2Vsin,,,/2,Rr,, 柱表面上滞止点的位置。 ,,,得,当时, 有两个解,分别是在滞止点, ,V,0,sin,,,1,s4,RV,,4RV, ,,,,,。当时，有一个解，当,,,,arcsin()arcsin(),1,,,,,,,,s1s2s,R,,R,,24RV,4sin4sin ,时无解，说明此时滞止点不在圆柱表面上。这时由V的偏导数式可以确定滞止点,1,4RV, ,,1,22