绕圆柱无环量流动和有环量流动流线分布图.docx

1、绕圆柱无环量流动和有环量流动流线分布图绕圆柱无环量流动和有环量流动流线分布图一、绕圆柱无环量流动由均匀直线流动和偶极流的流函数可叠加得绕圆柱无环量流动的流函数方程为:=8 r sin( 8) 兀相应流动的流线方程为：8 r sin( 0 ) - =C (常量)兀不同C值对应不同的流函数曲线，将极坐标形式变成直角坐标形式：8 y =C (常量)兀1.8=10, M=500, C=10,-10,0,20,-20ezplot( 10*y-500*y/(2*pi*(xA2+yA2)=10 )hold onezplot( 10*y-500*y/(2*pi*(xA2+yA2)=-10 )hold onez

2、plot( 10*y-500*y/(2*pi*(xA2+yA2)=0 )hold onezplot( 10*y-500*y/(2*pi*(xA2+yA2)=20 )hold onezplot( 10*y-500*y/(2*pi*(xA2+yA2)=-20 )hold on2.8=20, M=520, C=20,-20,0,40,-40ezplot( 20*y-520*y/(2*pi*(xA2+yA2)=20 )hold onezplot( 20*y-520*y/(2*pi*(xA2+yA2)=-20 )hold onezplot( 20*y-520*y/(2*pi*(xA2+yA2)=0 )h

3、old onezplot( 20*y-520*y/(2*pi*(xA2+yA2)=40 )hold onezplot( 20*y-520*y/(2*pi*(xA2+yA2)=-40 )hold on20 六520 y/(2 tt (x2+y2)=-40* i l I h i I i-S -4 -2 D 2 4 63.8=30, M=540, C=30,-30,0,60,-60ezplot( 30*y-540*y/(2*pi*(xA2+yA2)=30 )hold onezplot( 30*y-540*y/(2*pi*(xA2+yA2)=-30 )hold onezplot( 30*y-540*y

4、/(2*pi*(xA2+yA2)=0 )hold onezplot( 30*y-540*y/(2*pi*(xA2+yA2)=60 )hold onezplot( 30*y-540*y/(2*pi*(xA2+yA2)=-60 )hold on30 y-540 y/(2 打(x2+y2)-60-6 4 -2 D 2 4 64.oo=10O, M=1000, C=30,-30,0,60,-60ezplot( 100*y-1000*y/(2*pi*(xA2+yA2)=30 )hold onezplot( 100*y-1000*y/(2*pi*(xA2+yA2)=-30 )hold onezplot(

5、100*y-1000*y/(2*pi*(xA2+yA2)=0 )hold onezplot( 100*y-1000*y/(2*pi*(xA2+yA2)=60 )hold onezplot( 100*y-1000*y/(2*pi*(xA2+yA2)=-60 )hold on6100 y-1000 W(2 77 (k，/)=-60-G -4 *2 0 2 4 6X5.8=10, M=100, C=30,-30,0,60,-60ezplot( 10*y-100*y/(2*pi*(xA2+yA2)=30 )hold onezplot( 10*y-100*y/(2*pi*(xA2+yA2)=-30 )h

6、old onezplot( 10*y-100*y/(2*pi*(xA2+yA2)=0 )hold onezplot( 10*y-100*y/(2*pi*(xA2+yA2)=60 )hold onezplot( 10*y-100*y/(2*pi*(xA2+yA2)=-60 )hold on二、绕圆柱有环量流动由均匀直线流动、偶极流、点涡三者叠加后的复合流动，流函数方程为:， 、 ， r .w= oo (r)sin( 0)+ln(r)兀相应流动的流线方程为：F _ .、一 一 .(r-)sin( 9 )+一 ln(r)=C (吊重)兀不同C值对应不同的流函数曲线，将极坐标形式变成直角坐标形式： r

7、 、一一( - ) + 一 ln =C (吊室)兀8=10, R=3, r=0.1, C=0,10,20,30,40ezplot( 10*y*(1-9/(xA2+yA2)+0.1*exp(sqrt(xA2+yA2)/(2*pi)=0 )hold onezplot( 10*y*(1-9/(xA2+yA2)+0.1*exp(sqrt(xA2+yA2)/(2*pi)=10 )hold onezplot( 10*y*(1-9/(xA2+yA2)+0.1*exp(sqrt(xA2+yA2)/(2*pi)=20 )hold onezplot( 10*y*(1-9/(xA2+yA2)+0.1*exp(sqr

8、t(xA2+yA2)/(2*pi)=30 )hold onezplot( 10*y*(1-9/(xA2+yA2)+0.1*exp(sqrt(xA2+yA2)/(2*pi)=40 )hold on10 y (1 -9/(A2+y2)0.1 exp(sqrt(x- 2)/(2 开)=40-6 -4 -2 0 2 4 68=20, R=2, r=0.1, C=0,20,-20,40,-40ezplot( 20*y*(1-4/(xA2+yA2)+0.1*exp(sqrt(xA2+yA2)/(2*pi)=0 )hold onezplot( 20*y*(1-4/(xA2+yA2)+0.1*exp(sqrt

9、(xA2+yA2)/(2*pi)=20 )hold onezplot( 20*y*(1-4/(xA2+yA2)+0.1*exp(sqrt(xA2+yA2)/(2*pi)=-20 )hold onezplot( 20*y*(1-4/(xA2+yA2)+0.1*exp(sqrt(xA2+yA2)/(2*pi)=40 )hold onezplot( 20*y*(1-4/(xA2+yA2)+0.1*exp(sqrt(xA2+yA2)/(2*pi)=-40 )hold on20 y (14/(1+/)+。J *p 佰qrl仪)=-402.4-6 4-2 0 2 4 68=50, R=4, r=0.1,

10、C=0,10,-10,20,-20ezplot( 50*y*(1-16/(xA2+yA2)+0.1*exp(sqrt(xA2+yA2)/(2*pi)=0 )hold on ezplot( 50*y*(1-16/(xA2+yA2)+0.1*exp(sqrt(xA2+yA2)/(2*pi)=10 )hold onezplot( 50*y*(1-16/(xA2+yA2)+0.1*exp(sqrt(xA2+yA2)/(2*pi)=-10 )hold on ezplot( 50*y*(1-16/(xA2+yA2)+0.1*exp(sqrt(xA2+yA2)/(2*pi)=20 )hold on ezpl

11、ot( 50*y*(1-16/(xA2+yA2)+0.1*exp(sqrt(xA2+yA2)/(2*pi)=-20 )hold on50 y (M6/(A2+y2)+0.1 expsqrt(xy2)/(2 兀)=-208=50, R=3, r=-0.3, C=0,20,-20ezplot( 50*y*(1-9/(xA2+yA2)-0.3*exp(sqrt(xA2+yA2)/(2*pi)=0 ，-66,-5 5) hold onezplot( 50*y*(1-9/(xA2+yA2)-0.3*exp(sqrt(xA2+yA2)/(2*pi)=-20 ,-66,-5 5) hold on ezplo

12、t( 50*y*(1-9/(xA2+yA2)-0.3*exp(sqrt(xA2+yA2)/(2*pi)=20 ,-66,-5 5) hold on50 y (V9/(x2+y3)-0.3 exp(sqrt(x2+/)y(2 开户208=1000, R=5, r =-1, C=0,20,-20,40 , -40ezplot( 1000*y*(1-25/(xA2+yA2)-1*exp(sqrt(xA2+yA2)/(2*pi)=0 ，-66,-5 5)hold onezplot( 1000*y*(1-25/(xA2+yA2)-1*exp(sqrt(xA2+yA2)/(2*pi)=-20 ,-66,-5 5)hold onezplot( 1000*y*(1-25/(xA2+yA2)-1*exp(sqrt(xA2+yA2)/(2*pi)=20 ,-66,-5 5)hold onezplot( 1000*y*(1-25/(xA2+yA2)-1*exp(sqrt(xA2+yA2)/(2*pi)=40 ,-66,-5 5)hold onezplot( 1000*y*(1-25/(xA2+yA2)-1*exp(sqrt(xA2+yA2)/(2*pi)=-40 ,-66,-5 5)hold on1000 y (1-25/(x2+y2)-1 exp(sqrt(x2+y2)/(2 tt)=-40-6-2