From 6372c796c62d60d2642ffc81365ba9e7a87fdc68 Mon Sep 17 00:00:00 2001 From: Adrian Kummerlaender Date: Mon, 21 Jan 2019 20:39:38 +0100 Subject: Add Knudsen functor to locate refinement areas --- apps/adrian/cylinder2d/cylinder2d.cpp | 12 +++++++----- 1 file changed, 7 insertions(+), 5 deletions(-) (limited to 'apps/adrian/cylinder2d') diff --git a/apps/adrian/cylinder2d/cylinder2d.cpp b/apps/adrian/cylinder2d/cylinder2d.cpp index 3baa4ef..2e18ebd 100644 --- a/apps/adrian/cylinder2d/cylinder2d.cpp +++ b/apps/adrian/cylinder2d/cylinder2d.cpp @@ -171,12 +171,14 @@ void getResults(const std::string& prefix, auto& sGeometry = grid.getSuperGeometry(); SuperVTMwriter2D vtmWriter(prefix + "cylinder2d"); - SuperLatticePhysVelocity2D velocity(sLattice, converter); - SuperLatticePhysPressure2D pressure(sLattice, converter); - SuperLatticeGeometry2D geometry(sLattice, sGeometry); + SuperLatticePhysVelocity2D velocity(sLattice, converter); + SuperLatticePhysPressure2D pressure(sLattice, converter); + SuperLatticeGeometry2D geometry(sLattice, sGeometry); + SuperLatticeKnudsen2D knudsen(sLattice); vtmWriter.addFunctor(geometry); vtmWriter.addFunctor(velocity); vtmWriter.addFunctor(pressure); + vtmWriter.addFunctor(knudsen); const int statIter = converter.getLatticeTime(maxPhysT/10.); @@ -235,8 +237,8 @@ int main(int argc, char* argv[]) sOnLatticeBoundaryCondition2D fineSBoundaryCondition(fineGrid.getSuperLattice()); createLocalBoundaryCondition2D(fineSBoundaryCondition); - const Vector fineExtend2 {5*cylinderD, 1.5*cylinderD}; - const Vector fineOrigin2 {1.25*cylinderD, (ly-fineExtend2[1])/2}; + const Vector fineExtend2 {1.75*cylinderD, 1.75*cylinderD}; + const Vector fineOrigin2 {2*cylinderD-fineExtend2[0]/2, 2*cylinderD-fineExtend2[1]/2}; auto& fineGrid2 = fineGrid.refine(fineOrigin2, fineExtend2); prepareGeometry(fineGrid2); -- cgit v1.2.3