LBM入门小案例matlab格式书写
于 2022-01-27 发布
文件大小:1.35 kB
0 131
下载积分: 2
下载次数: 1
代码说明：
代码注释清晰，适合初学者学习使用 %给定初始的密度，流速，设定粒子分布函数的初始值 rho=ones(nx,ny); %rho=1 u=zeros(nx,ny); %x方向速度分量 v=zeros(nx,ny); %y方向速度分量 f=zeros(nx,ny,9); %nx*ny个节点，每个节点有9个方向 for m=1:9 %粒子分布函数f(:,:,m)表示同一页的数值，f1=1/9 f2=1/9...1/36, 4/9 f(:,:,m)=w(m)*rho.*(1 ... +3*(e(1,m)*u+e(2,m)*v)... +9/2*(e(1,m)*u+e(2,m)*v).^2 ... -3/2*(u.^2+v.^2)); end ftemp=zeros(nx,ny,9); flag=ones(nx,ny); %nx*ny的单位矩阵Flag，用于区分上下边界的位置 flag(:,1)=0; %flag矩阵的第一列赋值为0 flag(:,ny)=0; %矩阵的ny列赋值为0
下载说明：请别用迅雷下载，失败请重下，重下不扣分！
发表评论
