fixed bugs in some cases

apps/gpu/DiluteGravityCurrents_Case2/DiluteGravityCurrents_Case2.cpp

@@ -109,9 +109,9 @@ int main(int argc, char *argv[])
 
         real dx = L / real(nx);
 
-        gridBuilder->addCoarseGrid(-12.5 * L, 0 * L, 0 * L, 12.5 * L, 1.5 * L, 1 * L, dx);
+        gridBuilder->addCoarseGrid(0 * L, 0 * L, 0 * L, 15 * L, 1 * L, 1 * L, dx);
 
-        gridBuilder->setPeriodicBoundaryCondition(true, true, false);
+        gridBuilder->setPeriodicBoundaryCondition(false, false, false);
 
         gridBuilder->buildGrids(lbmOrGks, false);
 
@@ -122,7 +122,7 @@ int main(int argc, char *argv[])
         if (lbmOrGks == LBM) {
             SPtr<Parameter> para = Parameter::make();
 
-            real DeltaT = 2.5 * (real)nx / para->getVelocityLB();
+            real DeltaT = 0.5 * (real)nx / para->getVelocityLB();
 
             const uint timeStepOut = (int)DeltaT;
             const uint timeStepEnd = (int)DeltaT * 12;
@@ -192,6 +192,10 @@ int main(int argc, char *argv[])
             // set boundary conditions
             //////////////////////////////////////////////////////////////////////////
 
+            gridBuilder->setNoSlipBoundaryCondition(SideType::PX);
+            gridBuilder->setNoSlipBoundaryCondition(SideType::MX);
+            gridBuilder->setNoSlipBoundaryCondition(SideType::PY);
+            gridBuilder->setNoSlipBoundaryCondition(SideType::MY);
             gridBuilder->setNoSlipBoundaryCondition(SideType::PZ);
             gridBuilder->setNoSlipBoundaryCondition(SideType::MZ);
 

apps/gpu/DiluteGravityCurrents_Case4/DiluteGravityCurrents_Case4.cpp

@@ -109,9 +109,9 @@ int main(int argc, char *argv[])
 
         real dx = L / real(nx);
 
-        gridBuilder->addCoarseGrid(-12.5 * L, 0 * L, 0 * L, 12.5 * L, 1.5 * L, 1 * L, dx);
+        gridBuilder->addCoarseGrid(0 * L, 0 * L, 0 * L, 15 * L, 1 * L, 1 * L, dx);
 
-        gridBuilder->setPeriodicBoundaryCondition(true, true, false);
+        gridBuilder->setPeriodicBoundaryCondition(false, false, false);
 
         gridBuilder->buildGrids(lbmOrGks, false);
 
@@ -122,7 +122,7 @@ int main(int argc, char *argv[])
         if (lbmOrGks == LBM) {
             SPtr<Parameter> para = Parameter::make();
 
-            real DeltaT = 2.5 * (real)nx / para->getVelocityLB();
+            real DeltaT = 0.5 * (real)nx / para->getVelocityLB();
 
             const uint timeStepOut = (int)DeltaT;
             const uint timeStepEnd = (int)DeltaT * 12;
@@ -192,6 +192,10 @@ int main(int argc, char *argv[])
             // set boundary conditions
             //////////////////////////////////////////////////////////////////////////
 
+            gridBuilder->setNoSlipBoundaryCondition(SideType::PX);
+            gridBuilder->setNoSlipBoundaryCondition(SideType::MX);
+            gridBuilder->setNoSlipBoundaryCondition(SideType::PY);
+            gridBuilder->setNoSlipBoundaryCondition(SideType::MY);
             gridBuilder->setNoSlipBoundaryCondition(SideType::PZ);
             gridBuilder->setNoSlipBoundaryCondition(SideType::MZ);
 

apps/gpu/DiluteGravityCurrents_Case6/DiluteGravityCurrents_Case6.cpp

@@ -109,9 +109,9 @@ int main(int argc, char *argv[])
 
         real dx = L / real(nx);
 
-        gridBuilder->addCoarseGrid(-12.5 * L, 0 * L, 0 * L, 12.5 * L, 1.5 * L, 1 * L, dx);
+        gridBuilder->addCoarseGrid(0 * L, 0 * L, 0 * L, 15 * L, 1 * L, 1 * L, dx);
 
-        gridBuilder->setPeriodicBoundaryCondition(true, true, false);
+        gridBuilder->setPeriodicBoundaryCondition(false, false, false);
 
         gridBuilder->buildGrids(lbmOrGks, false);
 
@@ -122,7 +122,7 @@ int main(int argc, char *argv[])
         if (lbmOrGks == LBM) {
             SPtr<Parameter> para = Parameter::make();
 
-            real DeltaT = 2.5 * (real)nx / para->getVelocityLB();
+            real DeltaT = 0.5 * (real)nx / para->getVelocityLB();
 
             const uint timeStepOut = (int)DeltaT;
             const uint timeStepEnd = (int)DeltaT * 12;
@@ -192,6 +192,10 @@ int main(int argc, char *argv[])
             // set boundary conditions
             //////////////////////////////////////////////////////////////////////////
 
+            gridBuilder->setNoSlipBoundaryCondition(SideType::PX);
+            gridBuilder->setNoSlipBoundaryCondition(SideType::MX);
+            gridBuilder->setNoSlipBoundaryCondition(SideType::PY);
+            gridBuilder->setNoSlipBoundaryCondition(SideType::MY);
             gridBuilder->setNoSlipBoundaryCondition(SideType::PZ);
             gridBuilder->setNoSlipBoundaryCondition(SideType::MZ);
 

apps/gpu/DiluteGravityCurrents_Case7/DiluteGravityCurrents_Case7.cpp

@@ -80,7 +80,7 @@ int main(int argc, char *argv[])
         const real Re   = 15000;// Set Reynolds number for given case
         const real Sc   = 1;    // Set Schmidt number
         const real L    = 1;    // Set height of channel (H)
-        const uint nx   = 140;  // Set the number of nodes discretizing channel depth
+        const uint nx   = 104;  // Set the number of nodes discretizing channel depth
         const real Gate = L;    // Set distance to gate (x0)
 
         // this testcase is implemented for LBM(GPU) only

apps/gpu/DiluteGravityCurrents_Case8/DiluteGravityCurrents_Case8.cpp

@@ -109,9 +109,9 @@ int main(int argc, char *argv[])
 
         real dx = L / real(nx);
 
-        gridBuilder->addCoarseGrid(-12.5 * L, 0 * L, 0 * L, 12.5 * L, 1.5 * L, 1 * L, dx);
+        gridBuilder->addCoarseGrid(0 * L, 0 * L, 0 * L, 15 * L, 1 * L, 1 * L, dx);
 
-        gridBuilder->setPeriodicBoundaryCondition(true, true, false);
+        gridBuilder->setPeriodicBoundaryCondition(false, false, false);
 
         gridBuilder->buildGrids(lbmOrGks, false);
 
@@ -122,7 +122,7 @@ int main(int argc, char *argv[])
         if (lbmOrGks == LBM) {
             SPtr<Parameter> para = Parameter::make();
 
-            real DeltaT = 2.5 * (real)nx / para->getVelocityLB();
+            real DeltaT = 0.5 * (real)nx / para->getVelocityLB();
 
             const uint timeStepOut = (int)DeltaT;
             const uint timeStepEnd = (int)DeltaT * 12;
@@ -192,6 +192,10 @@ int main(int argc, char *argv[])
             // set boundary conditions
             //////////////////////////////////////////////////////////////////////////
 
+            gridBuilder->setNoSlipBoundaryCondition(SideType::PX);
+            gridBuilder->setNoSlipBoundaryCondition(SideType::MX);
+            gridBuilder->setNoSlipBoundaryCondition(SideType::PY);
+            gridBuilder->setNoSlipBoundaryCondition(SideType::MY);
             gridBuilder->setNoSlipBoundaryCondition(SideType::PZ);
             gridBuilder->setNoSlipBoundaryCondition(SideType::MZ);