Pretreatment for counter-current spontaneous imbibition

grid_preprocess.m

%use Pretreatment (testTreatment.m) to get file 'Proceded_LBM_DVM_original_from_online_577_388.raw';
%And use the following command to get 'decimalFile_577_388.dat':
fileID=fopen('Proceded_LBM_DVM_original_from_online_577_388.raw');
A=fread(fileID,[577,388],'uint8=>uint');
dlmwrite('decimalFile_577_388.dat', A);

%then add the inlet layers
load decimalFile_577_388.dat
B=decimalFile_577_388;
finalmatrix=B(2:577,:);
dlmwrite('cut_spontaneous.dat',finalmatrix); % cut 577 by 1 to 576 in order to be divided by 24
load cut_spontaneous.dat
C=cut_spontaneous;
A=zeros(24,388);% add 24 lines on the x direction to make 600X388
A(:,1:352)=1;
addleft=A;
finalfinalmatrix=cat(1,addleft,C);
dlmwrite('Grid_try_3.dat', finalfinalmatrix);

Fracture 6 has the same dimension as fracture 5. The only difference is that fracture 6 added some channels to make it look less symmetrical.