流函数-涡量法的二维方腔流数值模拟.docVIP

流函数-涡量法的二维方腔流数值模拟 流函数-涡量法的二维方腔流数值模拟 PAGE 流函数-涡量法的二维方腔流数值模拟 流函数- 涡量法的二维方腔流数值模拟 基本方程： 在直角坐标系下，不可压非定常流体所满足的流函数涡量形式的N-S方程为 其中 为雷诺数 差分格式： 采用FTCS格式有： 对于本问题，将方腔四边同时分为等分，则有 故 在直角坐标系下，不可压定常流体所满足的流函数涡量形式的N-S方程为 其中 为雷诺数 差分格式： 采用FTCS格式有： 对于本问题，将方腔四边同时分为等分，则有，则有即 边界条件： 在腔体的两侧和顶边， （第二式由泰勒级数展开得到） 在底边 （第二式由泰勒级数展开得到） 其中代表边界，代表与边界相邻的节点。 而 即 Matlab程序为： 不可压非定常流体 clear; %参数设置 Re=10; %雷诺数取10,100，500，1000 L=1; %空穴几何尺寸 n=100; dh=L/n;%delta h dt=1e-4; %时间步长 psi=zeros(n+1,n+1); xi=zeros(n+1,n+1); rho=1; for k=1:1000000 err=0; %边界条件