LBM入门小案例matlab格式书写
2016-12-25
3 0 0
4.0 分
其他
如何获取积分?
代码注释清晰,适合初学者学习使用
%给定初始的密度,流速,设定粒子分布函数的初始值
rho=ones(nx,ny); %rho=1u=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
相关源码推荐
使用Java开发Android AOA Android开放式附件
0
0
暂无评分
Golang AOA Android开放式附件HID控制
0
0
暂无评分
pclite
0
0
暂无评分
GMSK线性接收机
0
0
暂无评分
传世AFT中变打金版版本--以前很火爆的,收藏很久
0
0
暂无评分
暂无评论