LBM-matlab代码(Red-and-Blue-Fluids)

%FLUIDPROPERTIES(2IMISCIBLEFLUIDS)physicalproperties%dPdL=Pressure/dL;%Externalpressuregradient[atm/cm]cP_visco=1;%[cP]1CPwaterviscosity20Cdensity=1;%fluiddensity%latticeprop.somega=1,cs2=1/3;%omega:relaxationfrequencyomegaR=0.9;%RedomegaB=0.8;%Blue%colourcompositioni.e.amountofRedandBluecolorsatstartfct_Red=0.40;%volumetricfractionofRedfluidfct_Blue=1-fct_Red;%fractionBluefluidSigma=0.0002;%surfacetensionparametermin_M_Grad=0.0001;%mingradientforre-colouringAdes_RNB=-1;%Ades_RNB=1ifredisonthesurfaceofthedryarea%N_c=9;%declarativestatemetsrelatedtothe2fluidsf=zeros([Nr,Mc,N_c]);%arrayoffluiddensitydistribution%f=zeros(Nr,Mc,N_c);Red=zeros(Nr,Mc,N_c);%e.g.oilBlue=zeros(Nr,Mc,N_c);%e.g.waterrho=zeros(Nr,Mc);%macro-scopicdensityrhoR=zeros(Nr,Mc);rhoB=zeros(Nr,Mc);%alsoforthe2fluidsrhoR_B=zeros(Nr,Mc);%Conc_fct=zeros(Nr,Mc);%if(0)isOpen=matlabpool('size')0ifisOpen==1,matlabpoolclose,endmatlabpoolopenlocal4,f=zeros([Nr,Mc,N_c],'single',codistributor());%arrayend