From 4b4887769564bff24783eceaaa69a1d72ebcde32 Mon Sep 17 00:00:00 2001 From: Adrian Kummerlaender Date: Tue, 8 Jan 2019 19:32:53 +0100 Subject: Modifiy params to promote vortex street formation --- apps/adrian/poiseuille2d/poiseuille2d.cpp | 15 +++++++-------- 1 file changed, 7 insertions(+), 8 deletions(-) (limited to 'apps/adrian/poiseuille2d') 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 const& converter, // Set material number for vertically centered obstacle { - const Vector extend {0.1, 0.1}; - const Vector origin {1.25, (ly-extend[1])/2}; - IndicatorCuboid2D obstacle(extend, origin); + const Vector origin {1.25, ly/2}; + IndicatorCircle2D obstacle(origin, 0.1); superGeometry.rename(1,2,obstacle); } @@ -113,7 +112,7 @@ void prepareLattice(UnitConverter 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 coarseExtend {lx, ly}; IndicatorCuboid2D coarseCuboid(coarseExtend, coarseOrigin); - auto coarseGrid = Grid2D::make(coarseCuboid, N, 0.8, Re); + auto coarseGrid = Grid2D::make(coarseCuboid, N, 0.57, Re); prepareGeometry(coarseGrid->getConverter(), coarseGrid->getSuperGeometry()); const Vector wantedFineExtend {2.0, 0.75}; -- cgit v1.2.3