流体力学上机.docVIP

能动硕21 王赵波 2120306014 第六章 6.27 解：根据题意， 点源强度Q = 24 , x正向速度 u = 10 m/s.用极坐标()表示速度流函数： ; 复合叠加得： 求滞止点 速度分量： 令得： 即滞止点坐标为（），故距原点距离为 绕流物体的形状 流函数 Matlab对不同C 值作图 子程序： function [x,y,streamfun]=simpyfuild(ftype,strense,x0,y0) close all nx=1000;xmin=-1.5;xmax=1.5; ny=1000;ymin=-1.5;ymax=1.5; [x,y]=meshgrid(linspace(xmin,xmax,nx),linspace(ymin,ymax,ny)); radius=inline(sqrt((x-x0).^2+(y-y0).^2),x,y,x0,y0); r=radius(x,y,x0,y0); theta=atan2(y-y0,x-x0); switch ftype case uniform %均匀来流 streamfun=strense*y; case source %源汇 streamfun=strense*theta/(2*pi); end figure(name,ftype,numbertitle,off) levmax=max(max(streamfun)); levmin=min(min(streamfun)); lev=linspace(levmin,levmax,100); contour(x,y,streamfun,lev) title([the figure of ftype]); axis equal xlabel(x) ylabel(y) 主程序： clc clear [x,y,z1]=simpyfuild(source,24,0,0); [x,y,z2]=simpyfuild(uniform,10,0,0); contour(x,y,z1+z2,-18:3:18) hold on plot(-6/5/pi,0,o) 绕物体的流动，物体的外轮廓线必是一条流线，流动沿物体表面分流时，分流线只能在滞止点A与物体外轮廓线相交，因此滞止点所在的流线为物体外轮廓线。即为A点所在的曲线。 （3）在和处物体外廓线之间的距离 将滞止点坐标代入流线方程，则，外廓线方程为： 由上式，可得时， 则在和处物体外廓线之间的距离 （4）处物体外廓线上的速度 将代入速度场方程得到：。 第七章 7.36 解:对两自由液面列伯努利方程得 其中，，为大气压 化简得： 沿程阻力损失： 综合上两式，可得： 雷诺数： 用科尔布鲁克公式求解： 其中，，故 迭代计算得: a=1; for i=1:100 f=(1/(-2.0*log10(1.80*10^-3+2.85*10^-6*sqrt(70+2.8/a))))^2 if abs(a-f)10^-8 break end a=f end f= 0.0334; 所以带入得： 流量 A