- 1、本文档共3页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
流函数- 涡量法的二维方腔流数值模拟 基本方程: t———时间步长, s; γ ...
流函数- 涡量法的二维方腔流数值模拟
基本方程:
t———时间步长, s;
γ———流体的运动粘度, m2 /s;
———涡量, s- 1;
———流函数, m2 /s;
v———速度矢量, m/s;
差分格式:
采用FTCS格式,对于有:
采用FTCS格式,对于有:
边界条件:
速度的边界条件:固体壁面处设置无滑移边界条件, 即u=0, v=0;
流函数的边界条件:在固体壁面及平板驱动处的流函数=0;
涡量的边界条件:
平板驱动处:
左右壁面上:;
网格划分:
采用等距结构化网格划分(40*40)
编程计算:
本算例采用MATLAB进行编译,其主要优势是语言简单,可以方便地描绘出方腔环流的等值线图等。
主要语句:
while norm(c1)1e-4|norm(c2)1e-4
n=n+1;
O1=O;
E1=E;
t=t+dt;
for i=2:I
for j=2:J
O(i,j)=O(i,j)-dt*(u(i,j)*(O(i+1,j)-O(i-1,j))/(2*dx)+v(i,j)*(O(i,j+1)-O(i,j-1))/(2*dy))+dt/re*((O(i+1,j)-2*O(i,j)+O(i-1,j))/dx^2+(O(i,j+1)-2*O(i,j)+O(i,j-1))/dy^2);
end
end
for i=2:I
for j=2:J
E(i,j)=E(i,j)+0.25*(E(i+1,j)+E(i-1,j)+E(i,j+1)+E(i,j-1)-4*E(i,j)+dx^2*O(i,j));
end
end
% 0为涡量,E流函数%
计算效果图:
分别设置雷诺数为200,500,1000
雷诺数Re=200时的流函数图及速度矢量图
雷诺数Re=500时的流函数图及速度矢量图
雷诺数Re=1000时的流函数图及速度矢量图
文档评论(0)