diff options
Diffstat (limited to 'apps/adrian/poiseuille2d')
-rw-r--r-- | apps/adrian/poiseuille2d/poiseuille2d.cpp | 6 |
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(); |