diff options
| -rw-r--r-- | apps/adrian/cylinder2d/cylinder2d.cpp | 14 |
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); |