diff options
| -rw-r--r-- | apps/adrian/poiseuille2d/poiseuille2d.cpp | 15 |
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}; |