绕圆柱无环量和有环量仿真设计

《绕圆柱无环量和有环量仿真设计》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、绕圆柱无环量和有环量仿真设计由流体力学的相关知识可知，所有的真实流体都具有粘性和一定程度的可压缩性，但是在一•些情况下，粘性及压缩性的影响较小，可以忽略，采用简化的理想不川■压缩模型能很好的近似实际流动，另一些情况则是为了揭示出流动的特性的冇价值的规律，采用简化的流动模型能使研究变得简中.。基本的平面势流包括均匀直线流动、点源和点汇、点涡和偶极流。根据流函数和速度势函数的定义以及将流体近似看作是理想不吋压缩流体的条件下，可以写出各个基本平面势流的流函数和速度势函数，得到流函数后，讨以利用数学编程MATLAB软件将函数曲线

2、画在坐标系下，更加直观的观察各种流动的特点以及不同基本势流叠加后的流动特性。本文即对相关内容进行了处理。柱的无环量流动由理论分析可知，均匀直线流动与偶极流叠加后可以表示均匀来流绕圆柱的无环量流动。设均匀直线流动的速度为错误！未找到引用源。，沿x轴正方向，偶极中心位于坐标原点、强度为M、偶极轴沿负x方句。由均匀直线流动和偶极流的流函数可叠加得绕岡柱的无环量流动的流函数方程为：屮=错误！未找到引用源。rsin(0)_错误！未找到引用源。因此符合流动的流线方程为：错误！未找到引用源。rsin(e)-错误！未找到引用源。=C(常

3、量)不同的C值对应不同的流函数曲线，对方程中的常熟错误！未找到引用源。、M以及C取不同的值编程副出一族流函数曲线，上诉方程是用极坐标给出的，为了编程方便，现将其变为宵角坐标表示的形式：因为r=错误！未找到引用源。；sin(e)=错误！未找到引用源。=错误！未找到引用源。，所以方程化为：错误！未找到引用源。y-错误！未找到引用源。=C(常量)MATLAB程序为ezplot(’23=50*y-557*y/(2*pi*(xl+yl))’)holdonezplot(’-23=50*y-557*y/(2*pi*(x"2+y"2))

4、’)holdonezplot(’0=50*y-557*y/(2*pi*(x'2+yD)’)holdonezplot(’46=50*y-557*y/(2*pi*(x*2+y*2))>)holdonezplot(’-46=50*y-557*y/(2*pi*(xl+yl))’)得到的绕圆柱的无环量流动的流线分布图为：二、绕x柱的有环量流动如果在绕岡柱的无环量流动之上再在岡心叠加一个点涡流动，由于点涡流动的流线为同心圆族，故圆柱面仍为流线，则可组成绕圆柱的有环量流动，设点涡为顺时针方向、强度为-r。由均匀直线流动、偶极流、点涡三

5、者叠加后的复合流动，其流函数为：屮=错误!未找到引用源。（r-错误！未找到引用源。）sin（0）+错误！未找到引用源。ln（r）其流线方程为：错误！未找到引用源。（r-错误！未找到引用源。）sin（e）+错误！未找到引用源。ln（r）=C（常量）同样将其化为直角坐标的形式为：错误！未找到引用源。（错误！未找到引用源。-错误！未找到引用源。）错误！未找到引用源。+错误！未找到引用源。In（错误!未找到引用源。）=C（常量）对其中的各个常数错误！未找到引用源。、R、r以及c取不同的值可以得到不同的流线分布图：X相应程序为:e

6、zplotC0=50*(sqrt(x2+y^2)-9/sqrt(x2+y^2))*y/sqrt(x^2+y^2)+0.6/(2*pi)*exp(sqrt(x^2+y2))，,[-66,-55])holdonezplot(’23=50*(sqrt(x2+y"2)-9/sqrt(x2+y"2))*y/sqrt(x2+y2)+0.6/(2郎i)*exp(sqrt(x2+y2))’，[-66,-55])holdonezplot(’40=50*(sqrt(x"2+y"2)-9/sqrt(x"2+y"2))*y/sqrt(x"2+y

7、’2)+0•6/(2*pi)*exp(sqrt(x2+yD)’，[-66,-55])holdonezplot(’60=50*(sqrt(x2+y"2)-9/sqrt(x2+y"2))*y/sqrt(x2+y2)+0•6/(2*pi)*exp(sqrt(x"2+y"2))’，[-66,-55])holdonezplot(’80=50*(sqrt(x2+y2)-9/sqrt(x2+y2))*y/sqrt(x2+y"2)+0.6/(2*pi)*exp(sqrt(x"2+y"2))’，[-66，-55])holdon当点涡的强度

8、减小时，可得流线分布图为：X对应MATLAB程序为:ezplot(’0=50*(sqrt(x^2+y^2)-9/sqrt(x^2+y^2))*y/sqrt(x"2+y"2)+0.3/(2*pi)*exp(sqrt(x"2+y"2))’，[-66,-55])holdonezplot(’23=50*(sqrt(x"2+y