Increase refined cylinder2d Reynolds number to 1000

1 files changed, 7 insertions, 7 deletions

diff --git a/apps/adrian/cylinder2d/cylinder2d.cpp b/apps/adrian/cylinder2d/cylinder2d.cpp
index 0cdfa7b..cb0cfc6 100644
--- a/apps/adrian/cylinder2d/cylinder2d.cpp
+++ b/apps/adrian/cylinder2d/cylinder2d.cpp
@@ -39,9 +39,9 @@ typedef double T;
 
 const T lx  = 8.0;            // length of the channel
 const T ly  = 2.0;            // height of the channel
-const int N = 50;             // resolution of the model
-const T Re  = 200.;           // Reynolds number
-const T baseTau = 0.8;        // Relaxation time of coarsest grid
+const int N = 100;            // resolution of the model
+const T Re  = 1000.;          // Reynolds number
+const T baseTau = 0.52;       // Relaxation time of coarsest grid
 const T maxPhysT = 60.;       // max. simulation time in s, SI unit
 const T physInterval = 0.25;  // interval for the convergence check in s
 const T residuum = 1e-5;      // residuum for the convergence check
@@ -86,7 +86,7 @@ void prepareGeometry(Grid2D<T,DESCRIPTOR>& grid)
   // Set material number for vertically centered obstacle
   {
     const Vector<T,2> origin {1.25, ly/2};
-    IndicatorCircle2D<T> obstacle(origin, 0.1);
+    IndicatorCircle2D<T> obstacle(origin, 0.15);
     sGeometry.rename(1,2,obstacle);
   }
 
@@ -118,7 +118,7 @@ void prepareLattice(Grid2D<T,DESCRIPTOR>& grid,
   sLattice.defineDynamics(sGeometry, 4, &bulkDynamics); // outflow
 
   sBoundaryCondition.addVelocityBoundary(sGeometry, 3, omega);
-  sBoundaryCondition.addVelocityBoundary(sGeometry, 4, omega);
+  sBoundaryCondition.addPressureBoundary(sGeometry, 4, omega);
 
   const T Lx = converter.getLatticeLength(lx);
   const T Ly = converter.getLatticeLength(ly);
@@ -221,8 +221,8 @@ int main(int argc, char* argv[])
   sOnLatticeBoundaryCondition2D<T, DESCRIPTOR> fineSBoundaryCondition(fineGrid.getSuperLattice());
   createLocalBoundaryCondition2D<T, DESCRIPTOR>(fineSBoundaryCondition);
 
-  const Vector<T,2> fineExtend2 {0.6, 0.4};
-  const Vector<T,2> fineOrigin2 {1.05, (ly-fineExtend2[1])/2};
+  const Vector<T,2> fineExtend2 {0.8, 0.6};
+  const Vector<T,2> fineOrigin2 {0.95, (ly-fineExtend2[1])/2};
 
   auto& fineGrid2 = fineGrid.refine(fineOrigin2, fineExtend2);
   prepareGeometry(fineGrid2);