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

Fixed the setup for the square tube test

Browse Source
This commit is contained in:
rexzheli
2016-09-02 14:38:34 -04:00
parent b34d1f1f29
commit e7987c9cc9
1 changed files with 3 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
tests/lbpm_squaretube_pp.cpp
View File

@@ -149,13 +149,14 @@ int main(int argc, char **argv)
double pore_vol;
sum=0;
//double wall_width = (Nx - TubeWidth)/2;
for (k=0;k<Nz;k++){
for (j=0;j<Ny;j++){
for (i=0;i<Nx;i++){
n = k*Nx*Ny + j*Nz + i;
// square capillary tube aligned with the z direction
Averages.SDs(i,j,k) = TubeWidth - fabs(i-0.5*Nx);
Averages.SDs(i,j,k) = min(Averages.SDs(i,j,k),TubeWidth-fabs(j-0.5*Ny));
Averages.SDs(i,j,k) = TubeWidth/2 - fabs(i-0.5*Nx);
Averages.SDs(i,j,k) = min(Averages.SDs(i,j,k),TubeWidth/2-fabs(j-0.5*Ny));
// Initialize phase positions
if (Averages.SDs(i,j,k) < 0.0){