From ac6ac9383dfcc3a688579dc1fb56874bc6ff392e Mon Sep 17 00:00:00 2001
From: Adrian Kummerlaender
Date: Thu, 17 Jan 2019 13:15:00 +0100
Subject: Increase refined cylinder2d Reynolds number to 1000
---
apps/adrian/cylinder2d/cylinder2d.cpp | 14 +++++++-------
1 file 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);
--
cgit v1.2.3