summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorAdrian Kummerlaender2019-01-17 13:15:00 +0100
committerAdrian Kummerlaender2019-06-24 15:17:09 +0200
commitac6ac9383dfcc3a688579dc1fb56874bc6ff392e (patch)
treec0dc9a74f7621f47c490e128142050dfa260dfc6
parent5b8bad24546c0c01b11113fcd9b212bc633b6757 (diff)
downloadgrid_refinement_openlb-ac6ac9383dfcc3a688579dc1fb56874bc6ff392e.tar
grid_refinement_openlb-ac6ac9383dfcc3a688579dc1fb56874bc6ff392e.tar.gz
grid_refinement_openlb-ac6ac9383dfcc3a688579dc1fb56874bc6ff392e.tar.bz2
grid_refinement_openlb-ac6ac9383dfcc3a688579dc1fb56874bc6ff392e.tar.xz
grid_refinement_openlb-ac6ac9383dfcc3a688579dc1fb56874bc6ff392e.zip
Increase refined cylinder2d Reynolds number to 1000
-rw-r--r--apps/adrian/cylinder2d/cylinder2d.cpp14
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);