基本信息

---

下载

   友情提示：（无需注册或充值，赞助后即可获取资源下载链接）

---

     嘿，亲！知识可是无价之宝呢，但咱这精心整理的资料也耗费了不少心血呀。小小地破费一下，绝对物超所值哦！如有下载和支付问题，请联系我们QQ(微信同号)：813200300

本次赞助数额为： 3 元　

支付宝支付

微信支付

   源码介绍

---

学习3DLBM的好例子

% D3Q19 LATTICE CONSTANTS  
tNS = [1/3, 1/18,1/18,1/18,1/18,1/18,1/18, 1/36,1/36,1/36,1/36,1/36,1/36,1/36,1/36,1/36,1/36,1/36,1/36];  
cxNS = [ 0, 1,-1, 0, 0, 0, 0, 1, 1,-1,-1, 1,-1, 1,-1, 0, 0, 0, 0];  
cyNS = [ 0, 0, 0, 1,-1, 0, 0, 1,-1, 1,-1, 0 ,0, 0, 0, 1, 1,-1,-1];  
czNS = [ 0, 0, 0, 0, 0, 1,-1, 0, 0, 0, 0, 1, 1,-1,-1, 1,-1 ,1,-1];  
 
[y,x,z] = meshgrid(1:ly,1:lx,lz);  
 
drho = 0.001;  
delta_rho = -drho*(1-2.0*rand(lx,ly,lz));  
 
cent_x1= round(lx/3); cent_y1= round(ly/2); cent_z1= round(lz/2);  
cent_x2= round(2*lx/3); cent_y2= round(ly/2); cent_z2= round(lz/2);  
rad= round(lx/10);  
 
circ_in= zeros(lx,ly,lz); circ_out= ones(lx,ly,lz);  
 
circ_in(cent_x1-rad:cent_x1 rad,cent_y1-rad:cent_y1 rad,cent_z1-rad:cent_z1 rad)= 1;  
circ_in(cent_x2-rad:cent_x2 rad,cent_y2-rad:cent_y2 rad,cent_z2-rad:cent_z2 rad)= 1;  
circ_out(cent_x1-rad:cent_x1 rad,cent_y1-rad:cent_y1 rad,cent_z1-rad:cent_z1 rad)= 0;  
circ_out(cent_x2-rad:cent_x2 rad,cent_y2-rad:cent_y2 rad,cent_z2-rad:cent_z2 rad)= 0;  
 
% INITIAL CONDITION FOR BOTH DISTRIBUTION FUNCTIONS: (T=0) ==> TIn(i) = t(i)  
for i=1:19  
% fIn(i,1:lx,1:ly,1:lz) = tNS(i)*circ_in*2;  
% gIn(i,1:lx,1:ly,1:lz) = tNS(i)*circ_out*2;  
fIn(i,1:lx,1:ly,1:lz) = tNS(i).*(1.0 delta_rho);  
gIn(i,1:lx,1:ly,1:lz) = tNS(i).*(1.0 - delta_rho);  
end