Modifiy params to promote vortex street formation

1 files changed, 7 insertions, 8 deletions

diff --git a/apps/adrian/poiseuille2d/poiseuille2d.cpp b/apps/adrian/poiseuille2d/poiseuille2d.cpp
index f4ebffa..20d2585 100644
--- a/apps/adrian/poiseuille2d/poiseuille2d.cpp
+++ b/apps/adrian/poiseuille2d/poiseuille2d.cpp
@@ -37,8 +37,8 @@ typedef double T;
 
 const T lx  = 4.0;            // length of the channel
 const T ly  = 1.0;            // height of the channel
-const int N = 30;             // resolution of the model
-const T Re  = 100.;           // Reynolds number
+const int N = 60;             // resolution of the model
+const T Re  = 500.;           // Reynolds number
 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
@@ -81,9 +81,8 @@ void prepareGeometry(UnitConverter<T,DESCRIPTOR> const& converter,
 
   // Set material number for vertically centered obstacle
   {
-    const Vector<T,2> extend {0.1, 0.1};
-    const Vector<T,2> origin {1.25, (ly-extend[1])/2};
-    IndicatorCuboid2D<T> obstacle(extend, origin);
+    const Vector<T,2> origin {1.25, ly/2};
+    IndicatorCircle2D<T> obstacle(origin, 0.1);
     superGeometry.rename(1,2,obstacle);
   }
 
@@ -113,7 +112,7 @@ void prepareLattice(UnitConverter<T,DESCRIPTOR> const& converter,
   sLattice.defineDynamics(superGeometry, 4, &bulkDynamics); // outflow
 
   sBoundaryCondition.addVelocityBoundary(superGeometry, 3, omega);
-  sBoundaryCondition.addPressureBoundary(superGeometry, 4, omega);
+  sBoundaryCondition.addVelocityBoundary(superGeometry, 4, omega);
 
   const T Lx = converter.getLatticeLength(lx);
   const T Ly = converter.getLatticeLength(ly);
@@ -162,7 +161,7 @@ void getResults(const std::string& prefix,
     vtmWriter.createMasterFile();
   }
 
-  if (iT%100==0) {
+  if (iT%10==0) {
     vtmWriter.write(iT);
   }
 
@@ -184,7 +183,7 @@ int main(int argc, char* argv[])
   const Vector<T,2> coarseExtend {lx, ly};
   IndicatorCuboid2D<T> coarseCuboid(coarseExtend, coarseOrigin);
 
-  auto coarseGrid = Grid2D<T,DESCRIPTOR>::make(coarseCuboid, N, 0.8, Re);
+  auto coarseGrid = Grid2D<T,DESCRIPTOR>::make(coarseCuboid, N, 0.57, Re);
   prepareGeometry(coarseGrid->getConverter(), coarseGrid->getSuperGeometry());
 
   const Vector<T,2> wantedFineExtend {2.0, 0.75};