计算流体.zip_CFD_流体_流体计算_计算流体_计算流体力学

采用 C 型网格，求解无源项的 Laplace 方程，画出网格。

1）求解

其中，

，

，

各项均采用中心差分，可得

，

其中，

，

，

2）结果如图所示

3）源程序（采用 matlab 编程）

M=151;N=241;nx=10;ny=5;

x=zeros(M,N);y=zeros(M,N);x1=zeros(M,N);y1=zeros(M,N);

%第一行

b1=4;

for j=1:N

a=(j-1)/(N-1);

x(1,j)=(exp(a*b1)-1)/(exp(b1)-1)*nx;

end

y(1,i)=0.1781*x(1,i)^0.5-0.0756*x(1,i)-0.2122*x(1,i)^2+0.1705*x(1,i)^3-0.0609*x(1,

i)^4;

end

y(1,j:N)=0;

%第 M 行

for j=101:N

x(M,j)=(j-101)/(N-101)*nx;

end

y(M,101:N)=ny;