%?A?short?and?simple?gravity-driven?LBM?solver?based?on?the?code?snippets?
%?in?Sukop?and?Thorne‘s?‘Lattice?Boltzmann?Modeling‘

%?Note?indexing?differences?between?book‘s?C?code?and?MATLAB:?
%?C?uses?0?for?the?first?index?value?while?MATLAB?starts?at?one.
%?Numerous?changes?are?needed.?In?some?places?I?have?just?
%?explicitly?added?one?to?the?C?index.

close（‘a(chǎn)ll‘）;clear（‘a(chǎn)ll‘）

display（‘initialize‘）

LY=100
LX=300
tau?=?1
g=0.001

%set?solid?nodes?at?walls?on?top?and?bottom
is_solid_node=zeros（LYLX）;
for?i=1:LX
????is_solid_node（1i）=1;
????is_solid_node（10i）=1;
????is_solid_node（20i）=1;
????is_solid_node（60i）=1;
????is_solid_node（80i）=1;
????is_solid_node（LYi）=1;
end
for?i=5:（LY-5）
????is_solid_node（i40）=1;
????is_solid_node（i80）=1;
????is_solid_node（i120）=1;
????is_solid_node（i160）=1;
????is_solid_node（i180）=1;
????is_solid_node（i220）=1;
end

display（‘solid?nodes‘）
is_solid_node

%define?initial?density?and?fs
rho=ones（LYLX）;

f（::1）?=?（4./9.?）*rho;
f（::2）?=?（1./9.?）*rho;
f（::3）?=?（1./9.?）*rho;
f（::4）?=?（1./9.?）*rho;
f（::5）?=?（1./9.?）*rho;
f（::6）?=?（1./36.）*rho;
f（::7）?=?（1./36.）*rho;
f（::8）?=?（1./36.）*rho;
f（::9）?=?（1./36.）*rho;

display（‘intitial?f‘）

f;

%define?lattice?velocity?vectors

ex（0+1）=?0;?ey（0+1）=?0;
ex（1+1）=?1;?ey（1+1）=?0;
ex（2+1）=?0;?ey（2+1）=?1;
ex（3+1）=-1;?ey（3+1）=?0;
ex（4+1）=?0;?ey（4+1）=-1;
ex（5+1）=?1;?ey（5+1）=?1;
ex（6+1）=-1;?ey（6+1）=?1;
ex（7+1）=-1;?ey（7+1）=-1;
ex（8+1）=?1;?ey（8+1）=-1;


for?ts=1:300?%Time?loop
????
????ts
????
????%?Computing?macroscopic?density?rho?and?velocity?u=（uxuy）.
????for?j=1:LY?
????????
????????for?i=1:LX
????????????
????????????u_x（ji）?=?0.0;
????????????u_y（ji）?=?0.0;
????????????rho（ji）?=?0.0;
????????????
????????????if?~is_solid_node（ji）
????????????????
????????????????for?a=0:8
????????????????????
????????????????????rho（ji）?=?rho（ji）?+?f（jia+1）;
???????????????????????????????????????
????????????????????u_x（ji）?=?u_x（ji）?+?ex（a+1）*f（jia+1）;
????????????????????u_y（ji）?=?u_y（ji）?+?ey（a+1）*f（jia+1）;
????????????????????
????????????????end
????????????????
????????????????u_x（ji）?=?u_x（ji）/rho（ji）;
????????????????u_y（ji）?=?u_y（ji）/rho（ji）;
????????????????
????????????end
????????????
????????????%add?space?matricies?for?plotting
????????????x（ji）=i;
????????????y（ji）=j;
????????????
????????end
????end
????
???????
????%?Compute?the?equilibrium?distribution?function?feq.
????f1=3.;
????f2=9./2.;
????f3=3./2.;
????
????for?j=1:LY
????????
????????for?i=1:LX
????????????
????????????if?~is_solid_node（ji）
????????????????
????????????????rt0?=?（4./9.?）*rho（ji）;
????????????????rt1?=?（1./9.?）*rho（ji）;
????????????????rt2?=?（1./36.）*rho（ji）;
????????????????ueqxij?=??u_x（ji）+tau*g;?%add?forcing
????????????????ueqyij?=??u_y（ji）;
????????????????uxsq???=??ueqxij?*?ueqxij;%changes?from?book?here!?See?Book‘s?

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

?????文件????1420974??2009-07-03?14:33??格子Boltzmann\A?3D?Lattice?Boltzmann?Code?for?Modeling?Flow?and?Multi-Component?Dispersion.pdf

?????文件??????10372??2009-07-03?14:13??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\ajaxtabs.js

?????文件????????591??2009-07-03?14:13??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\application_pdf.png

?????文件????????295??2009-07-03?14:13??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\at.png

?????文件????????964??2009-07-03?14:13??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\autocomplete.js

?????文件??????12917??2009-07-03?14:13??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\autosuggest.js

?????文件????????382??2009-07-03?14:13??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\br.png

?????文件??????10115??2009-07-03?14:12??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\browserSniffer.js

?????文件????????356??2009-07-03?14:13??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\ca.png

?????文件????????211??2009-07-03?14:13??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\ch.png

?????文件????????370??2009-07-03?14:13??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\cn.png

?????文件???????2370??2009-07-03?14:26??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\common.css

?????文件????????180??2009-07-03?14:13??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\de.png

?????文件???????1008??2009-07-03?14:26??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\dynCalendar.css

?????文件????????190??2009-07-03?14:26??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\eth-ie6.css

?????文件????????623??2009-07-03?14:26??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\eth-my.css

?????文件??????35698??2009-07-03?14:26??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\eth.css

?????文件????????505??2009-07-03?14:13??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\fi.png

?????文件????????276??2009-07-03?14:13??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\fr.png

?????文件????????644??2009-07-03?14:13??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\gb.png

?????文件??????34788??2009-07-03?14:12??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\global.js

?????文件????????541??2009-07-03?14:13??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\gr.png

?????文件???????2816??2009-07-03?14:13??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\hh.png

?????文件????????677??2009-07-03?14:13??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\hk.png

?????文件????????269??2009-07-03?14:25??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\hp.png

?????文件???????3521??2009-07-03?14:13??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\html_optimized.js

?????文件????????260??2009-07-03?14:13??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\in.png

?????文件????????778??2009-07-03?14:13??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\information.png

?????文件????????300??2009-07-03?14:13??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\ir.png

?????文件????????214??2009-07-03?14:13??格子Boltzmann\ETH?ETH?E-Collection?Lattice?Boltzmann?method?for?thermal?compressible?flows?-?ETH?E-Collection.files\it.png

............此處省略63個(gè)文件信息