LBM绕柱Matlab代码

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% cylinder.m: Channel flow past a cylinderical        
%             obstacle using a LB method            
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Lattice Boltzmann sample written in C
% Copyright （C） 2006 Jonas Latt
% Address: Rue General Dufour 24  1211 Geneva 4 Switzerland 
% E-mail: Jonas.Latt@cui.unige.ch
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% This program is free software; you can redistribute it and/or
% modify it under the terms of the GNU General Public License
% as published by the Free Software Foundation; either version 2
% of the License or （at your option） any later version.
% This program is distributed in the hope that it will be useful
% but WITHOUT ANY WARRANTY; without even t