• 大小: 3.12MB
    文件类型: .rar
    金币: 2
    下载: 1 次
    发布日期: 2023-09-21
  • 语言: 其他
  • 标签: LBM  OpenLB  3D  

资源简介

资源里包含学习格子Boltzmann方法的经典论文《A 3D Lattice Boltzmann Code for Modeling Flow and Multi-Component Dispersion》和《Lattice Boltzmann Method for Thermal Compressible Flows》及OpenLB软件库源代码,以及几篇中文硕士论文推荐

资源截图

代码片段和文件信息

% 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(‘all‘);clear(‘all‘)

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 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件    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个文件信息

评论

共有 条评论