diff --git a/apps/gpu/LBM/DrivenCavity/DrivenCavity.cpp b/apps/gpu/LBM/DrivenCavity/DrivenCavity.cpp index 869989462593edf14b198efa598f95c1b7051ba9..9f31c9358620331dd45dd1ff1a32f29b677a6d13 100644 --- a/apps/gpu/LBM/DrivenCavity/DrivenCavity.cpp +++ b/apps/gpu/LBM/DrivenCavity/DrivenCavity.cpp @@ -95,7 +95,7 @@ const real dt = (real)1.0e-3; //0.5e-3; const uint nx = 64; -std::string path("."); +std::string path("D:/out/DrivenCavity/new"); std::string simulationName("DrivenCavityChim"); @@ -123,8 +123,18 @@ void multipleLevel(const std::string& configPath) real dx = L / real(nx); - gridBuilder->addCoarseGrid(-0.5 * L, -0.5 * L, -0.5 * L, - 0.5 * L, 0.5 * L, 0.5 * L, dx); + //gridBuilder->addCoarseGrid(-0.5 * L, -0.5 * L, -0.5 * L, + // 0.5 * L, 0.5 * L, 0.5 * L, dx); + + gridBuilder->addCoarseGrid(-2.0 * dx, -0.5 * L, -0.5 * L, + 2.0 * dx, 0.5 * L, 0.5 * L, dx); + + auto refBox = new Cuboid(-0.1 * L, -0.1 * L, -0.1 * L, + 0.1 * L, 0.1 * L, 0.1 * L); + + gridBuilder->addGrid(refBox, 1); + + gridBuilder->setNumberOfLayers(0, 0); gridBuilder->setPeriodicBoundaryCondition(false, false, false); @@ -166,13 +176,15 @@ void multipleLevel(const std::string& configPath) para->setDevices(std::vector<uint>{(uint)0}); + para->setOutputPath( path ); + para->setOutputPrefix( simulationName ); para->setFName(para->getOutputPath() + "/" + para->getOutputPrefix()); para->setPrintFiles(true); - para->setMaxLevel(1); + para->setMaxLevel(2); para->setVelocity(velocityLB); para->setViscosity(viscosityLB); @@ -202,6 +214,10 @@ void multipleLevel(const std::string& configPath) //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// + gridBuilder->writeGridsToVtk(path + "/grid/"); + + //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// + SPtr<CudaMemoryManager> cudaMemoryManager = CudaMemoryManager::make(para); SPtr<GridProvider> gridGenerator = GridProvider::makeGridGenerator(gridBuilder, para, cudaMemoryManager); diff --git a/src/gpu/GridGenerator/grid/GridImp.h b/src/gpu/GridGenerator/grid/GridImp.h index b096f5ff85dcd725ff065dbb6fc31d75c016c869..9a8c209d2e6d3113c44482473ef591e0af4cb44d 100644 --- a/src/gpu/GridGenerator/grid/GridImp.h +++ b/src/gpu/GridGenerator/grid/GridImp.h @@ -28,7 +28,7 @@ // //! \file GridImp.h //! \ingroup grid -//! \author Soeren Peters, Stephan Lenz, Martin Sch�nherr +//! \author Soeren Peters, Stephan Lenz, Martin Schoenherr //======================================================================================= #ifndef GRID_IMP_H #define GRID_IMP_H @@ -197,7 +197,6 @@ public: void fixRefinementIntoWall(uint xIndex, uint yIndex, uint zIndex, int dir); void findStopperNode(uint index); void findEndOfGridStopperNode(uint index); - void findEndOfGridStopperPeriodicNode(uint index); void findSolidStopperNode(uint index); void findBoundarySolidNode(uint index);