summaryrefslogtreecommitdiff
path: root/apps/adrian/poiseuille2d/poiseuille2d.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'apps/adrian/poiseuille2d/poiseuille2d.cpp')
-rw-r--r--apps/adrian/poiseuille2d/poiseuille2d.cpp6
1 files changed, 3 insertions, 3 deletions
diff --git a/apps/adrian/poiseuille2d/poiseuille2d.cpp b/apps/adrian/poiseuille2d/poiseuille2d.cpp
index d92934b..19dc4d7 100644
--- a/apps/adrian/poiseuille2d/poiseuille2d.cpp
+++ b/apps/adrian/poiseuille2d/poiseuille2d.cpp
@@ -38,9 +38,9 @@ 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 int N = 50; // resolution of the model
const T Re = 100.; // Reynolds number
-const T baseTau = 0.8; // Relaxation time of coarsest grid
+const T uMax = 0.02; // Max lattice speed
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
@@ -184,7 +184,7 @@ int main(int argc, char* argv[])
const Vector<T,2> coarseExtend {lx/2, ly};
IndicatorCuboid2D<T> coarseCuboid(coarseExtend, coarseOrigin);
- Grid2D<T,DESCRIPTOR> coarseGrid(coarseCuboid, N, baseTau, Re);
+ Grid2D<T,DESCRIPTOR> coarseGrid(coarseCuboid, N, LatticeVelocity<T>(uMax), Re);
prepareGeometry(coarseGrid);
const T coarseDeltaX = coarseGrid.getConverter().getPhysDeltaX();