summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rw-r--r--apps/adrian/poiseuille2d/poiseuille2d.cpp15
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};