diff --git a/apps/gpu/LBM/BoundaryLayer/BoundaryLayer.cpp b/apps/gpu/LBM/BoundaryLayer/BoundaryLayer.cpp index 51b890d205c52650343c9aa2608f9bc3e237da56..7d2d79d8c21f1966a5b2369f2649df2efc3b4bcb 100644 --- a/apps/gpu/LBM/BoundaryLayer/BoundaryLayer.cpp +++ b/apps/gpu/LBM/BoundaryLayer/BoundaryLayer.cpp @@ -238,11 +238,11 @@ void multipleLevel(const std::string& configPath) if(isFirstSubDomain || isMidSubDomain) { xGridMax += overlap; - xGridMin -= overlap; + // xGridMin -= overlap; } if(isLastSubDomain || isMidSubDomain) { - xGridMax += overlap; + // xGridMax += overlap; xGridMin -= overlap; } @@ -252,8 +252,8 @@ void multipleLevel(const std::string& configPath) { gridBuilder->setNumberOfLayers(4,0); real xMaxRefinement = readPrecursor? xGridMax-H: xGridMax; //Stop refinement some distance before outlet if domain ist not periodic - gridBuilder->addGrid( new Cuboid( xGridMin+dx, 0.f, 0.f, xMaxRefinement, L_y, 0.5*L_z) , 1 ); - + // gridBuilder->addGrid( new Cuboid( xGridMin+dx, 0.f, 0.f, xMaxRefinement, L_y, 0.5*L_z) , 1 ); + gridBuilder->addGrid( new Cuboid( 0.f, 0.f, 0.f, 5000.0, L_y, 0.5*L_z) , 1 ); para->setMaxLevel(2); scalingFactory.setScalingFactory(GridScalingFactory::GridScaling::ScaleCompressible); } @@ -302,8 +302,9 @@ void multipleLevel(const std::string& configPath) { if(isFirstSubDomain || nProcs == 1) { - auto precursor = createFileCollection(precursorDirectory + "/precursor", FileType::VTK); - gridBuilder->setPrecursorBoundaryCondition(SideType::MX, precursor, nTReadPrecursor); + // auto precursor = createFileCollection(precursorDirectory + "/precursor", FileType::VTK); + // gridBuilder->setPrecursorBoundaryCondition(SideType::MX, precursor, nTReadPrecursor); + gridBuilder->setVelocityBoundaryCondition(SideType::MX, velocityLB, 0.0, 0.0); } if(isLastSubDomain || nProcs == 1) diff --git a/src/gpu/VirtualFluids_GPU/Parameter/Parameter.cpp b/src/gpu/VirtualFluids_GPU/Parameter/Parameter.cpp index df53a0d9300ead34ad96b5257b1c8b5debad3a70..d1568932d9802d1060e37d26ed37fa48de1abb05 100644 --- a/src/gpu/VirtualFluids_GPU/Parameter/Parameter.cpp +++ b/src/gpu/VirtualFluids_GPU/Parameter/Parameter.cpp @@ -861,6 +861,10 @@ real Parameter::getScaledForceRatio(int level) { return this->getForceRatio()*(level+1); } +real Parameter::getScaledStressRatio(int level) +{ + return this->getVelocityRatio()*this->getVelocityRatio(); +} void Parameter::setRealX(real RealX) { ic.RealX = RealX; diff --git a/src/gpu/VirtualFluids_GPU/Parameter/Parameter.h b/src/gpu/VirtualFluids_GPU/Parameter/Parameter.h index cf71b20bb22f35948f567b5501c066f1732efdbf..de5e2b0a4d983efdd272649962bfc9cfa60c7f19 100644 --- a/src/gpu/VirtualFluids_GPU/Parameter/Parameter.h +++ b/src/gpu/VirtualFluids_GPU/Parameter/Parameter.h @@ -778,6 +778,8 @@ public: real getScaledDensityRatio(int level); //! \returns the pressure ratio in SI/LB units scaled to the respective level real getScaledPressureRatio(int level); + //! \returns the stress ratio in SI/LB units scaled to the respective level + real getScaledStressRatio(int level) //! \returns the time ratio in SI/LB units scaled to the respective level real getScaledTimeRatio(int level); //! \returns the length ratio in SI/LB units scaled to the respective level diff --git a/src/gpu/VirtualFluids_GPU/PreCollisionInteractor/Probes/Probe.cu b/src/gpu/VirtualFluids_GPU/PreCollisionInteractor/Probes/Probe.cu index 3bde6e8d80d947b2870935c967a00385f6ecfb24..44f76220f77db01d63a153cccec40bad7cca5449 100644 --- a/src/gpu/VirtualFluids_GPU/PreCollisionInteractor/Probes/Probe.cu +++ b/src/gpu/VirtualFluids_GPU/PreCollisionInteractor/Probes/Probe.cu @@ -187,7 +187,7 @@ void Probe::init(Parameter* para, GridProvider* gridProvider, CudaMemoryManager* this->velocityRatio = std::bind(&Parameter::getScaledVelocityRatio, para, _1); this->densityRatio = std::bind(&Parameter::getScaledDensityRatio, para, _1); this->forceRatio = std::bind(&Parameter::getScaledForceRatio, para, _1); - this->stressRatio = std::bind(&Parameter::getScaledPressureRatio, para, _1); + this->stressRatio = std::bind(&Parameter::getScaledStressRatio, para, _1); this->viscosityRatio = std::bind(&Parameter::getScaledViscosityRatio, para, _1); this->nondimensional = std::bind(&Probe::getNondimensionalConversionFactor, this, _1);