matlab计算涡度的函数_流函数涡量法的二维方腔流数值模拟matlab编程.doc

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

基本方程：

在直角坐标系下，不可压非定常流体所满足的流函数涡量形式的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;

%边界