From 946cb3702e6baacec1faaeab9ad66dfd99320b55 Mon Sep 17 00:00:00 2001 From: Adrian Kummerlaender Date: Thu, 10 Jan 2019 21:55:43 +0100 Subject: Fix grid alignment of refinement area --- apps/adrian/poiseuille2d/poiseuille2d.cpp | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) (limited to 'apps/adrian') diff --git a/apps/adrian/poiseuille2d/poiseuille2d.cpp b/apps/adrian/poiseuille2d/poiseuille2d.cpp index df166e1..b16a19b 100644 --- a/apps/adrian/poiseuille2d/poiseuille2d.cpp +++ b/apps/adrian/poiseuille2d/poiseuille2d.cpp @@ -37,9 +37,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 = 64; // resolution of the model -const T Re = 500.; // Reynolds number -const T baseTau = 0.57; // Relaxation time of coarsest grid +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 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 @@ -188,8 +188,8 @@ int main(int argc, char* argv[]) prepareGeometry(coarseGrid->getConverter(), coarseGrid->getSuperGeometry()); const Vector wantedFineExtend {3.0, 1.5}; - const Vector fineOrigin = coarseGrid->alignToGrid({0.8, (ly-wantedFineExtend[1])/2}); - const Vector fineExtend = coarseGrid->alignToGrid(fineOrigin + wantedFineExtend) - fineOrigin; + const Vector fineOrigin = coarseGrid->alignLocationToGrid({0.8, (ly-wantedFineExtend[1])/2}); + const Vector fineExtend = coarseGrid->alignExtendToGrid(wantedFineExtend); auto fineGrid = &coarseGrid->refine(fineOrigin, fineExtend); prepareGeometry(fineGrid->getConverter(), fineGrid->getSuperGeometry()); @@ -222,8 +222,8 @@ int main(int argc, char* argv[]) createLocalBoundaryCondition2D(fineSBoundaryCondition); const Vector wantedFineExtend2 {0.6, 0.4}; - const Vector fineOrigin2 = fineGrid->alignToGrid({1.05, (ly-wantedFineExtend2[1])/2}); - const Vector fineExtend2 = fineGrid->alignToGrid(fineOrigin2 + wantedFineExtend2) - fineOrigin2; + const Vector fineOrigin2 = fineGrid->alignLocationToGrid({1.05, (ly-wantedFineExtend2[1])/2}); + const Vector fineExtend2 = fineGrid->alignExtendToGrid(wantedFineExtend2); auto fineGrid2 = &fineGrid->refine(fineOrigin2, fineExtend2); prepareGeometry(fineGrid2->getConverter(), fineGrid2->getSuperGeometry()); -- cgit v1.2.3