原创力文档- 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）