LBM boiling

Lattice Boltzmann Method (LBM) Shan-Chen 单组份伪势模型模拟沸腾。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% boiling.m: the main program for simulating droplet motion on surfaces
%            with Shan Chen lattice Boltzmann model
%                        
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Shan Chen Lattice Boltzmann sample in Matlab
% Copyright Wei Gong
% Address: Nottingham NG7 2RD UK
% E-mail: wei.gong@nottingham.ac.uk
% Reference: Li Qing et al. “Lattice Boltzmann modeling of boiling heat 
%            transfer: The boiling curve and the effects of wettability.“ 
%            International Journal of Heat and Mass Transfer 85 （2015）: 
%            787-796.

clear clc

global g;

format long;

t_max = single（1000000）;    % maximum iteration
Nwri  = single（100）;        % output data frequency

constant;         % constant setting
initialization;   % initial condition setting

% load time1;

% iteration （first 100 iters without gravity）
g = single（0）;      
for time = single（1:100）    
    disp（time）;
    % streaming
    for i = 1:9
        ff（i : :） = ...  
            circshift（fe（i : :） [0 ex（i） ey（i）]）;
    end
    % Zou-He boundary conditions
    ff = zouhe（ff）; 
    % macro parameters
    macrop;
    % interaction forces
    [S Fx Fy] = forces（ p ff ）; 
    % collision
    collision; 
    % visualization
    visua;
end

% iteration （after 100 iters with gravity）
g = single（3e-5）;
for time = single（101:t_max）
    disp（time）;
    % streaming
    for i = 1:9
        ff（i : :） = ...  
            circshift（fe（i : :） [0 ex（i） ey（i）]）;
    end
    % Zou-He boundary conditions
    ff = zouhe（ff）;
    % macro parameters
    macrop;
    % interaction forces
    [S Fx Fy] = forces（ p ff ）; 
    % collision
    collision; 
    % visualization
    visua;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-02-20 23:11  LBM boiling\
     文件        1910  2017-08-16 20:33  LBM boiling\boiling.m
     文件        1088  2017-08-16 20:33  LBM boiling\collision.m
     文件        3336  2017-08-16 20:33  LBM boiling\constant.m
     文件        1806  2017-08-16 20:33  LBM boiling\forces.m
     文件        2523  2017-08-16 20:33  LBM boiling\initialization.m
     文件         896  2017-08-16 20:33  LBM boiling\macrop.m
     文件        2227  2017-08-16 20:33  LBM boiling\RK.m
     文件        1834  2017-08-16 20:33  LBM boiling\solveEOS.m
     文件        1224  2017-08-16 20:33  LBM boiling\visua.m
     文件        1008  2017-08-16 20:33  LBM boiling\zouhe.m