D2Q9模型--基于MATLAB的LBM代码

哇咔咔又扒代码啦，感谢我的难兄老田同志贡献了自己的账号。这个代码上面没署名，不知道是哪个朋友写的，侵权立刻删哈。

还是老规矩先宣传一下QQ群群： 格子玻尔兹曼救星：293267908。

%    ------------------
%    ->                   二维通道，左进右出，上下为壁面
%    ------------------
% D2Q9模型
% C6  C2  C5    
%   \  |  /     
% C3  C9  C1    %原C0在此表示为C9
%   /  |  \     
% C7  C4  C8   

clear all
clc
m=41;
n=1001;
dt=1;
dx=1;
dy=1;
x=0:dx:1000;
y=0:dy:40;
alpha=0.02;%粘度
omega = 1/(alpha*3+0.5);
uo=0.2;%初速度
rhoo=5;%初始密度
mstep=40000;

%初始化
rho=ones(m