OilfieldToolsNet/LBPM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Found the (a) reason why the velocity BC were misbehaving

Browse Source
This commit is contained in:
James E McClure
2015-07-23 15:31:11 -04:00
parent ac7e410cba
commit 90cbc71a45
1 changed files with 18 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

28
tests/lbpm_color_simulator.cpp
View File

@@ -364,13 +364,6 @@ int main(int argc, char **argv)
if (Averages.SDs(n) > 0.0){ if (Averages.SDs(n) > 0.0){
sum++; sum++;
} }
// The following turns off communication if external BC are being set
if (BoundaryCondition > 0){
if (kproc==0 && k==0) id[n]=0;
if (kproc==0 && k==1) id[n]=0;
if (kproc==nprocz-1 && k==Nz-2) id[n]=0;
if (kproc==nprocz-1 && k==Nz-1) id[n]=0;
}
} }
} }
} }
@@ -382,6 +375,21 @@ int main(int argc, char **argv)
fread(id,1,N,IDFILE); fread(id,1,N,IDFILE);
fclose(IDFILE); fclose(IDFILE);
for ( k=0;k<Nz;k++){
for ( j=0;j<Ny;j++){
for ( i=0;i<Nx;i++){
n = k*Nx*Ny+j*Nx+i;
// The following turns off communication if external BC are being set
if (BoundaryCondition > 0){
if (kproc==0 && k==0) id[n]=0;
if (kproc==0 && k==1) id[n]=0;
if (kproc==nprocz-1 && k==Nz-2) id[n]=0;
if (kproc==nprocz-1 && k==Nz-1) id[n]=0;
}
}
}
}
// Set up kstart, kfinish so that the reservoirs are excluded from averaging // Set up kstart, kfinish so that the reservoirs are excluded from averaging
int kstart,kfinish; int kstart,kfinish;
kstart = 1; kstart = 1;
@@ -822,13 +830,13 @@ int main(int argc, char **argv)
Averages.ComponentAverages(); Averages.ComponentAverages();
Averages.SortBlobs(); Averages.SortBlobs();
Averages.PrintComponents(timestep); Averages.PrintComponents(timestep);
*/ //************************************************************************/ //************************************************************************/
/*
int NumberComponents_NWP = ComputeGlobalPhaseComponent(Dm.Nx-2,Dm.Ny-2,Dm.Nz-2,Dm.rank_info,Averages.PhaseID,1,Averages.Label_NWP); int NumberComponents_NWP = ComputeGlobalPhaseComponent(Dm.Nx-2,Dm.Ny-2,Dm.Nz-2,Dm.rank_info,Averages.PhaseID,1,Averages.Label_NWP);
printf("Number of non-wetting phase components: %i \n ",NumberComponents_NWP); printf("Number of non-wetting phase components: %i \n ",NumberComponents_NWP);
DeviceBarrier(); DeviceBarrier();
CopyToHost(Averages.Phase.get(),Phi,N*sizeof(double)); CopyToHost(Averages.Phase.get(),Phi,N*sizeof(double));
*/
// Create the MeshDataStruct // Create the MeshDataStruct
fillHalo<double> fillData(Dm.rank_info,Nx-2,Ny-2,Nz-2,1,1,1,0,1); fillHalo<double> fillData(Dm.rank_info,Nx-2,Ny-2,Nz-2,1,1,1,0,1);
std::vector<IO::MeshDataStruct> meshData(1); std::vector<IO::MeshDataStruct> meshData(1);