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

fix bug in TestPoisuille

Browse Source
This commit is contained in:
James McClure 2019-10-25 12:07:01 +02:00
parent ac5ec4492e
commit f3f08722d3
1 changed files with 17 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
tests/TestPoiseuille.cpp
View File

@ -16,6 +16,8 @@ void ParallelPlates(ScaLBL_MRTModel &MRT){
int Nx = MRT.Nx;
int Ny = MRT.Ny;
int Nz = MRT.Nz;
Array<char> id_solid(Nx,Ny,Nz);
for (k=0;k<Nz;k++){
for (j=0;j<Ny;j++){
for (i=0;i<Nx;i++){
@ -23,9 +25,24 @@ void ParallelPlates(ScaLBL_MRTModel &MRT){
if (i<2) MRT.Mask->id[n] = 0;
else if (i>Nx-3) MRT.Mask->id[n] = 0;
else MRT.Mask->id[n]=1;
if (MRT.Mask->id[n] == 0) id_solid(i,j,k) = 1;
else id_solid(i,j,k) = 0;
}
}
}
// Initialize the signed distance function
for (int k=0;k<Nz;k++){
for (int j=0;j<Ny;j++){
for (int i=0;i<Nx;i++){
int n=k*Nx*Ny+j*Nx+i;
// Initialize distance to +/- 1
MRT.Distance(i,j,k) = 2.0*double(id_solid(i,j,k))-1.0;
}
}
}
if (rank==0) printf("Initialized solid phase -- Converting to Signed Distance function \n");
CalcDist(MRT.Distance,id_solid,*MRT.Dm);
if (rank == 0) cout << "Domain set." << endl;
}
//***************************************************************************************