LBM Midrange Repulsion

基于 MATLAB 的 Midrange repulsion 伪势多相流模型。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% collision.m: cillision step
%                        
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Shan Chen Lattice Boltzmann sample in Matlab
% Copyright Wei Gong
% Address: Nottingham NG7 2RD UK
% E-mail: wei.gong@nottingham.ac.uk

u_n   = ee‘*reshape（u 2 lxy）; 
u_squ = u（1 : :）.^2 + u（2 : :）.^2;
fequi = t_k‘*reshape（rho 1 lxy）.*（1 + u_n/c_squ + u_n.^2/2/c_squ^2 - ...
    ones（9 1）*reshape（u_squ 1 lxy）/2/c_squ）; % equilibrium function
fn    = reshape（fequi 9 lx ly） + ...
    （1 - 1/tau）*（ff - reshape（fequi 9 lx ly））;
% fn    = reshape（fequi 9 lx ly） + ...
%     （1 - 1/tau）*（ff - reshape（fequi 9 lx ly）） + ...
%     reshape（ee‘*reshape（F 2 lxy）.*（t_k‘*ones（1 lxy））/c_squ 9 lx ly）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-02-20 23:27  Two phase flow - midrange repulsion\
     文件         841  2017-08-16 20:35  Two phase flow - midrange repulsion\collision.m
     文件        2080  2017-08-16 20:35  Two phase flow - midrange repulsion\constant.m
     文件        1375  2017-08-16 20:35  Two phase flow - midrange repulsion\forces.m
     文件         888  2017-08-16 20:35  Two phase flow - midrange repulsion\initialization.m
     文件         992  2017-08-16 20:35  Two phase flow - midrange repulsion\tflow.m
     文件        1260  2017-08-16 20:35  Two phase flow - midrange repulsion\visua.m