From 5eb4f8c195fa685284d0c154a5bbfcfa8c1a7a19 Mon Sep 17 00:00:00 2001
From: HenrikAsmuth <henrik.asmuth@geo.uu.se>
Date: Tue, 13 Dec 2022 14:06:37 +0100
Subject: [PATCH] Fix wrong scaling factor for stress in Probes
---
apps/gpu/LBM/BoundaryLayer/BoundaryLayer.cpp | 13 +++++++------
src/gpu/VirtualFluids_GPU/Parameter/Parameter.cpp | 4 ++++
src/gpu/VirtualFluids_GPU/Parameter/Parameter.h | 2 ++
.../PreCollisionInteractor/Probes/Probe.cu | 2 +-
4 files changed, 14 insertions(+), 7 deletions(-)
diff --git a/apps/gpu/LBM/BoundaryLayer/BoundaryLayer.cpp b/apps/gpu/LBM/BoundaryLayer/BoundaryLayer.cpp
index 51b890d20..7d2d79d8c 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 df53a0d93..d1568932d 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 cf71b20bb..de5e2b0a4 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 3bde6e8d8..44f76220f 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);
--
GitLab