From c4b77485b61650a0b93a7898803f6f760ed3f459 Mon Sep 17 00:00:00 2001
From: Anna Wellmann <a.wellmann@tu-bs.de>
Date: Thu, 15 Dec 2022 10:33:30 +0000
Subject: [PATCH] Finer grid, but slower
---
apps/gpu/LBM/ActuatorLine/ActuatorLine.cpp | 21 +-
apps/gpu/LBM/ActuatorLine/Pole.stl | 550 +++++++++++++++---
.../LBM/ActuatorLine/configActuatorLine.txt | 2 +-
.../CumulantK17Sponge_Device.cu | 4 +-
4 files changed, 484 insertions(+), 93 deletions(-)
diff --git a/apps/gpu/LBM/ActuatorLine/ActuatorLine.cpp b/apps/gpu/LBM/ActuatorLine/ActuatorLine.cpp
index 4aa89fe67..410e77a83 100644
--- a/apps/gpu/LBM/ActuatorLine/ActuatorLine.cpp
+++ b/apps/gpu/LBM/ActuatorLine/ActuatorLine.cpp
@@ -99,9 +99,9 @@ void multipleLevel(const std::string& configPath)
const real velocity = config.getValue<real>("Velocity");
- const real L_x = 10*reference_diameter;
- const real L_y = 5*reference_diameter;
- const real L_z = 6*reference_diameter;
+ const real L_x = 6*reference_diameter;
+ const real L_y = 4*reference_diameter;
+ const real L_z = 4*reference_diameter;
const real viscosity = 1.56e-5;
@@ -126,22 +126,21 @@ void multipleLevel(const std::string& configPath)
const real dx = reference_diameter/real(nodes_per_diameter);
- real turbPos[3] = {static_cast<real>(1.5*reference_diameter), 2*reference_diameter, 3*reference_diameter};
+ real turbPos[3] = {1*reference_diameter, 2*reference_diameter, 2*reference_diameter};
gridBuilder->addCoarseGrid(0.0, 0.0, 0.0,
L_x, L_y, L_z, dx);
gridBuilder->setNumberOfLayers(4,0);
- gridBuilder->addGrid( new Cuboid( turbPos[0]-1.0*reference_diameter, turbPos[1]-1.1*reference_diameter, turbPos[2]-1.1*reference_diameter,
- turbPos[0]+5.0*reference_diameter, turbPos[1]+1.1*reference_diameter, turbPos[2]+1.1*reference_diameter) , 1 );
+ gridBuilder->addGrid( new Cuboid( turbPos[0]-0.4*reference_diameter, turbPos[1]-0.8*reference_diameter, turbPos[2]-0.8*reference_diameter,
+ turbPos[0]+3.0*reference_diameter, turbPos[1]+0.8*reference_diameter, turbPos[2]+0.8*reference_diameter) , 1 );
para->setMaxLevel(2);
scalingFactory.setScalingFactory(GridScalingFactory::GridScaling::ScaleCompressible);
std::string stlPath = "./apps/gpu/LBM/ActuatorLine/Pole.stl";
- Object *sphere = TriangularMesh::make(stlPath);
-
- gridBuilder->addGeometry(sphere);
+ // Object *sphere = TriangularMesh::make(stlPath);
+ // gridBuilder->addGeometry(sphere);
gridBuilder->setPeriodicBoundaryCondition(false, false, false);
@@ -208,8 +207,8 @@ void multipleLevel(const std::string& configPath)
gridBuilder->setVelocityBoundaryCondition(SideType::PZ, velocityLB, 0.0, 0.0);
gridBuilder->setPressureBoundaryCondition(SideType::PX, 0.0);
- gridBuilder->setNoSlipBoundaryCondition(SideType::GEOMETRY);
- bcFactory.setGeometryBoundaryCondition(BoundaryConditionFactory::NoSlipBC::NoSlipCompressible);
+ // gridBuilder->setNoSlipBoundaryCondition(SideType::GEOMETRY);
+ // bcFactory.setGeometryBoundaryCondition(BoundaryConditionFactory::NoSlipBC::NoSlipCompressible);
bcFactory.setVelocityBoundaryCondition(BoundaryConditionFactory::VelocityBC::VelocityAndPressureCompressible);
bcFactory.setPressureBoundaryCondition(BoundaryConditionFactory::PressureBC::OutflowNonReflective);
diff --git a/apps/gpu/LBM/ActuatorLine/Pole.stl b/apps/gpu/LBM/ActuatorLine/Pole.stl
index 5b6453f68..9373b1f54 100644
--- a/apps/gpu/LBM/ActuatorLine/Pole.stl
+++ b/apps/gpu/LBM/ActuatorLine/Pole.stl
@@ -1,142 +1,534 @@
solid Visualization Toolkit generated SLA File
- facet normal 0.866025 0 -0.5
+ facet normal 0.987689 0 -0.156432
outer loop
- vertex 214 256 378
- vertex 214 -4 378
- vertex 209 256 369.34
+ vertex 148 252 252
+ vertex 148 0 252
+ vertex 147.608 252 249.528
endloop
endfacet
- facet normal -0.866025 0 0.5
+ facet normal -0.987689 0 0.156432
outer loop
- vertex 209 -4 369.34
- vertex 209 256 369.34
- vertex 214 -4 378
+ vertex 147.608 0 249.528
+ vertex 147.608 252 249.528
+ vertex 148 0 252
endloop
endfacet
- facet normal 0.866025 0 -0.5
+ facet normal 0.987689 0 -0.156432
outer loop
- vertex 209 256 369.34
- vertex 209 -4 369.34
- vertex 199 256 369.34
+ vertex 147.608 252 249.528
+ vertex 147.608 0 249.528
+ vertex 146.472 252 247.298
endloop
endfacet
- facet normal -0.866025 0 0.5
+ facet normal -0.987689 0 0.156432
outer loop
- vertex 199 -4 369.34
- vertex 199 256 369.34
- vertex 209 -4 369.34
+ vertex 146.472 0 247.298
+ vertex 146.472 252 247.298
+ vertex 147.608 0 249.528
endloop
endfacet
- facet normal 0.866025 0 -0.5
+ facet normal 0.987689 0 -0.156432
outer loop
- vertex 199 256 369.34
- vertex 199 -4 369.34
- vertex 194 256 378
+ vertex 146.472 252 247.298
+ vertex 146.472 0 247.298
+ vertex 144.702 252 245.528
endloop
endfacet
- facet normal -0.866025 0 0.5
+ facet normal -0.987689 0 0.156432
outer loop
- vertex 194 -4 378
- vertex 194 256 378
- vertex 199 -4 369.34
+ vertex 144.702 0 245.528
+ vertex 144.702 252 245.528
+ vertex 146.472 0 247.298
endloop
endfacet
- facet normal 0.866025 0 -0.5
+ facet normal 0.987689 0 -0.156432
outer loop
- vertex 194 256 378
- vertex 194 -4 378
- vertex 199 256 386.66
+ vertex 144.702 252 245.528
+ vertex 144.702 0 245.528
+ vertex 142.472 252 244.392
endloop
endfacet
- facet normal -0.866025 0 0.5
+ facet normal -0.987689 0 0.156432
outer loop
- vertex 199 -4 386.66
- vertex 199 256 386.66
- vertex 194 -4 378
+ vertex 142.472 0 244.392
+ vertex 142.472 252 244.392
+ vertex 144.702 0 245.528
endloop
endfacet
- facet normal 0.866025 0 -0.5
+ facet normal 0.987689 0 -0.156432
outer loop
- vertex 199 256 386.66
- vertex 199 -4 386.66
- vertex 209 256 386.66
+ vertex 142.472 252 244.392
+ vertex 142.472 0 244.392
+ vertex 140 252 244
endloop
endfacet
- facet normal -0.866025 0 0.5
+ facet normal -0.987689 0 0.156432
outer loop
- vertex 209 -4 386.66
- vertex 209 256 386.66
- vertex 199 -4 386.66
+ vertex 140 0 244
+ vertex 140 252 244
+ vertex 142.472 0 244.392
endloop
endfacet
- facet normal 0.866025 0 -0.5
+ facet normal 0.987689 0 -0.156432
outer loop
- vertex 209 256 386.66
- vertex 209 -4 386.66
- vertex 214 256 378
+ vertex 140 252 244
+ vertex 140 0 244
+ vertex 137.528 252 244.392
endloop
endfacet
- facet normal -0.866025 0 0.5
+ facet normal -0.987689 0 0.156432
outer loop
- vertex 214 -4 378
- vertex 214 256 378
- vertex 209 -4 386.66
+ vertex 137.528 0 244.392
+ vertex 137.528 252 244.392
+ vertex 140 0 244
endloop
endfacet
- facet normal 0.5 0 -0.866026
+ facet normal 0.987689 0 -0.156432
outer loop
- vertex 214 256 378
- vertex 209 256 369.34
- vertex 209 256 386.66
+ vertex 137.528 252 244.392
+ vertex 137.528 0 244.392
+ vertex 135.298 252 245.528
endloop
endfacet
- facet normal -0.5 0 0.866026
+ facet normal -0.987689 0 0.156432
outer loop
- vertex 209 256 386.66
- vertex 209 256 369.34
- vertex 199 256 386.66
+ vertex 135.298 0 245.528
+ vertex 135.298 252 245.528
+ vertex 137.528 0 244.392
endloop
endfacet
- facet normal 0 0.0332902 0.999446
+ facet normal 0.987689 0 -0.156432
outer loop
- vertex 209 256 369.34
- vertex 199 256 369.34
- vertex 199 256 386.66
+ vertex 135.298 252 245.528
+ vertex 135.298 0 245.528
+ vertex 133.528 252 247.298
+ endloop
+ endfacet
+ facet normal -0.987689 0 0.156432
+ outer loop
+ vertex 133.528 0 247.298
+ vertex 133.528 252 247.298
+ vertex 135.298 0 245.528
+ endloop
+ endfacet
+ facet normal 0.987689 0 -0.156432
+ outer loop
+ vertex 133.528 252 247.298
+ vertex 133.528 0 247.298
+ vertex 132.392 252 249.528
+ endloop
+ endfacet
+ facet normal -0.987689 0 0.156432
+ outer loop
+ vertex 132.392 0 249.528
+ vertex 132.392 252 249.528
+ vertex 133.528 0 247.298
+ endloop
+ endfacet
+ facet normal 0.987689 0 -0.156432
+ outer loop
+ vertex 132.392 252 249.528
+ vertex 132.392 0 249.528
+ vertex 132 252 252
+ endloop
+ endfacet
+ facet normal -0.987689 0 0.156432
+ outer loop
+ vertex 132 0 252
+ vertex 132 252 252
+ vertex 132.392 0 249.528
+ endloop
+ endfacet
+ facet normal 0.987689 0 -0.156432
+ outer loop
+ vertex 132 252 252
+ vertex 132 0 252
+ vertex 132.392 252 254.472
+ endloop
+ endfacet
+ facet normal -0.987689 0 0.156432
+ outer loop
+ vertex 132.392 0 254.472
+ vertex 132.392 252 254.472
+ vertex 132 0 252
+ endloop
+ endfacet
+ facet normal 0.987689 0 -0.156432
+ outer loop
+ vertex 132.392 252 254.472
+ vertex 132.392 0 254.472
+ vertex 133.528 252 256.702
+ endloop
+ endfacet
+ facet normal -0.987689 0 0.156432
+ outer loop
+ vertex 133.528 252 256.702
+ vertex 132.392 0 254.472
+ vertex 133.528 0 256.702
+ endloop
+ endfacet
+ facet normal 0.987689 0 -0.156432
+ outer loop
+ vertex 133.528 252 256.702
+ vertex 133.528 0 256.702
+ vertex 135.298 252 258.472
+ endloop
+ endfacet
+ facet normal -0.987689 0 0.156432
+ outer loop
+ vertex 135.298 0 258.472
+ vertex 135.298 252 258.472
+ vertex 133.528 0 256.702
+ endloop
+ endfacet
+ facet normal 0.987689 0 -0.156432
+ outer loop
+ vertex 135.298 252 258.472
+ vertex 135.298 0 258.472
+ vertex 137.528 252 259.608
+ endloop
+ endfacet
+ facet normal -0.987689 0 0.156432
+ outer loop
+ vertex 137.528 0 259.608
+ vertex 137.528 252 259.608
+ vertex 135.298 0 258.472
+ endloop
+ endfacet
+ facet normal 0.987689 0 -0.156432
+ outer loop
+ vertex 137.528 252 259.608
+ vertex 137.528 0 259.608
+ vertex 140 252 260
+ endloop
+ endfacet
+ facet normal -0.987689 0 0.156432
+ outer loop
+ vertex 140 0 260
+ vertex 140 252 260
+ vertex 137.528 0 259.608
+ endloop
+ endfacet
+ facet normal 0.987689 0 -0.156432
+ outer loop
+ vertex 140 252 260
+ vertex 140 0 260
+ vertex 142.472 252 259.608
+ endloop
+ endfacet
+ facet normal -0.987689 0 0.156432
+ outer loop
+ vertex 142.472 0 259.608
+ vertex 142.472 252 259.608
+ vertex 140 0 260
+ endloop
+ endfacet
+ facet normal 0.987689 0 -0.156432
+ outer loop
+ vertex 142.472 252 259.608
+ vertex 142.472 0 259.608
+ vertex 144.702 252 258.472
+ endloop
+ endfacet
+ facet normal -0.987689 0 0.156432
+ outer loop
+ vertex 144.702 0 258.472
+ vertex 144.702 252 258.472
+ vertex 142.472 0 259.608
+ endloop
+ endfacet
+ facet normal 0.987689 0 -0.156432
+ outer loop
+ vertex 144.702 252 258.472
+ vertex 144.702 0 258.472
+ vertex 146.472 252 256.702
+ endloop
+ endfacet
+ facet normal -0.987689 0 0.156432
+ outer loop
+ vertex 146.472 0 256.702
+ vertex 146.472 252 256.702
+ vertex 144.702 0 258.472
+ endloop
+ endfacet
+ facet normal 0.987689 0 -0.156432
+ outer loop
+ vertex 146.472 252 256.702
+ vertex 146.472 0 256.702
+ vertex 147.608 252 254.472
+ endloop
+ endfacet
+ facet normal -0.987689 0 0.156432
+ outer loop
+ vertex 147.608 252 254.472
+ vertex 146.472 0 256.702
+ vertex 147.608 0 254.472
+ endloop
+ endfacet
+ facet normal 0.987689 0 -0.156432
+ outer loop
+ vertex 147.608 252 254.472
+ vertex 147.608 0 254.472
+ vertex 148 252 252
+ endloop
+ endfacet
+ facet normal -0.987689 0 0.156432
+ outer loop
+ vertex 148 0 252
+ vertex 148 252 252
+ vertex 147.608 0 254.472
+ endloop
+ endfacet
+ facet normal -0.156432 0 -0.987689
+ outer loop
+ vertex 132.392 252 254.472
+ vertex 133.528 252 256.702
+ vertex 132 252 252
+ endloop
+ endfacet
+ facet normal 0.156431 0.00306919 0.987684
+ outer loop
+ vertex 132 252 252
+ vertex 133.528 252 256.702
+ vertex 132.392 252 249.528
endloop
endfacet
facet normal 0 0 -1
outer loop
- vertex 194 256 378
- vertex 199 256 386.66
- vertex 199 256 369.34
+ vertex 132.392 252 249.528
+ vertex 133.528 252 256.702
+ vertex 133.528 252 247.298
+ endloop
+ endfacet
+ facet normal -0 0.00450917 0.99999
+ outer loop
+ vertex 133.528 252 247.298
+ vertex 133.528 252 256.702
+ vertex 135.298 252 245.528
+ endloop
+ endfacet
+ facet normal 0.156435 0 -0.987688
+ outer loop
+ vertex 135.298 252 245.528
+ vertex 133.528 252 256.702
+ vertex 137.528 252 244.392
+ endloop
+ endfacet
+ facet normal -0.156432 0.00583797 0.987671
+ outer loop
+ vertex 137.528 252 244.392
+ vertex 133.528 252 256.702
+ vertex 140 252 244
+ endloop
+ endfacet
+ facet normal -0.309017 0 0.951056
+ outer loop
+ vertex 133.528 252 256.702
+ vertex 135.298 252 258.472
+ vertex 140 252 244
+ endloop
+ endfacet
+ facet normal 0.30901 0.00702308 -0.951033
+ outer loop
+ vertex 135.298 252 258.472
+ vertex 137.528 252 259.608
+ vertex 140 252 244
+ endloop
+ endfacet
+ facet normal -0.156434 0 0.987688
+ outer loop
+ vertex 137.528 252 259.608
+ vertex 140 252 260
+ vertex 140 252 244
+ endloop
+ endfacet
+ facet normal 0.156434 0 -0.987688
+ outer loop
+ vertex 140 252 244
+ vertex 140 252 260
+ vertex 142.472 252 244.392
+ endloop
+ endfacet
+ facet normal -0.309009 0.00702306 0.951033
+ outer loop
+ vertex 142.472 252 244.392
+ vertex 140 252 260
+ vertex 144.702 252 245.528
+ endloop
+ endfacet
+ facet normal 0.309005 0.00884943 -0.951019
+ outer loop
+ vertex 140 252 260
+ vertex 142.472 252 259.608
+ vertex 144.702 252 245.528
+ endloop
+ endfacet
+ facet normal -0.156435 0 0.987688
+ outer loop
+ vertex 142.472 252 259.608
+ vertex 144.702 252 258.472
+ vertex 144.702 252 245.528
+ endloop
+ endfacet
+ facet normal 0.156428 0.00944581 -0.987644
+ outer loop
+ vertex 144.702 252 258.472
+ vertex 146.472 252 256.702
+ vertex 144.702 252 245.528
+ endloop
+ endfacet
+ facet normal 0 0 1
+ outer loop
+ vertex 146.472 252 256.702
+ vertex 147.608 252 254.472
+ vertex 144.702 252 245.528
+ endloop
+ endfacet
+ facet normal 0 0 -1
+ outer loop
+ vertex 144.702 252 245.528
+ vertex 147.608 252 254.472
+ vertex 146.472 252 247.298
+ endloop
+ endfacet
+ facet normal -0.156428 0.00944581 0.987644
+ outer loop
+ vertex 146.472 252 247.298
+ vertex 147.608 252 254.472
+ vertex 147.608 252 249.528
+ endloop
+ endfacet
+ facet normal 0.156434 0 -0.987688
+ outer loop
+ vertex 148 252 252
+ vertex 147.608 252 249.528
+ vertex 147.608 252 254.472
+ endloop
+ endfacet
+ facet normal 0.156434 0 -0.987688
+ outer loop
+ vertex 147.608 0 254.472
+ vertex 146.472 0 256.702
+ vertex 148 0 252
+ endloop
+ endfacet
+ facet normal -0.156434 0 0.987688
+ outer loop
+ vertex 148 0 252
+ vertex 146.472 0 256.702
+ vertex 147.608 0 249.528
+ endloop
+ endfacet
+ facet normal 0.309002 0.00980961 -0.951011
+ outer loop
+ vertex 147.608 0 249.528
+ vertex 146.472 0 256.702
+ vertex 146.472 0 247.298
+ endloop
+ endfacet
+ facet normal -0.309017 0 0.951056
+ outer loop
+ vertex 146.472 0 247.298
+ vertex 146.472 0 256.702
+ vertex 144.702 0 245.528
+ endloop
+ endfacet
+ facet normal 0.45397 0.00944579 -0.890967
+ outer loop
+ vertex 144.702 0 245.528
+ vertex 146.472 0 256.702
+ vertex 142.472 0 244.392
+ endloop
+ endfacet
+ facet normal -0.45399 0 0.891007
+ outer loop
+ vertex 142.472 0 244.392
+ vertex 146.472 0 256.702
+ vertex 140 0 244
+ endloop
+ endfacet
+ facet normal -0.309005 0.00884946 0.951019
+ outer loop
+ vertex 146.472 0 256.702
+ vertex 144.702 0 258.472
+ vertex 140 0 244
+ endloop
+ endfacet
+ facet normal 0.309017 0 -0.951057
+ outer loop
+ vertex 144.702 0 258.472
+ vertex 142.472 0 259.608
+ vertex 140 0 244
+ endloop
+ endfacet
+ facet normal -0.156429 0.00803518 0.987656
+ outer loop
+ vertex 142.472 0 259.608
+ vertex 140 0 260
+ vertex 140 0 244
+ endloop
+ endfacet
+ facet normal 0.156432 0.00583801 -0.987672
+ outer loop
+ vertex 140 0 244
+ vertex 140 0 260
+ vertex 137.528 0 244.392
+ endloop
+ endfacet
+ facet normal -0.309017 0 0.951056
+ outer loop
+ vertex 137.528 0 244.392
+ vertex 140 0 260
+ vertex 135.298 0 245.528
+ endloop
+ endfacet
+ facet normal 0.309017 0 -0.951056
+ outer loop
+ vertex 140 0 260
+ vertex 137.528 0 259.608
+ vertex 135.298 0 245.528
+ endloop
+ endfacet
+ facet normal -0.156432 0.00583797 0.987671
+ outer loop
+ vertex 137.528 0 259.608
+ vertex 135.298 0 258.472
+ vertex 135.298 0 245.528
+ endloop
+ endfacet
+ facet normal 0.156435 0 -0.987688
+ outer loop
+ vertex 135.298 0 258.472
+ vertex 133.528 0 256.702
+ vertex 135.298 0 245.528
endloop
endfacet
- facet normal -0.866025 0 0.5
+ facet normal 0 0.00450917 0.99999
outer loop
- vertex 194 -4 378
- vertex 199 -4 369.34
- vertex 199 -4 386.66
+ vertex 133.528 0 256.702
+ vertex 132.392 0 254.472
+ vertex 135.298 0 245.528
endloop
endfacet
- facet normal 0.866025 0 -0.5
+ facet normal 0 0.00155373 -0.999999
outer loop
- vertex 199 -4 386.66
- vertex 199 -4 369.34
- vertex 209 -4 386.66
+ vertex 135.298 0 245.528
+ vertex 132.392 0 254.472
+ vertex 133.528 0 247.298
endloop
endfacet
- facet normal 0.499907 0.0192272 -0.865865
+ facet normal -0.156432 0 0.987689
outer loop
- vertex 199 -4 369.34
- vertex 209 -4 369.34
- vertex 209 -4 386.66
+ vertex 133.528 0 247.298
+ vertex 132.392 0 254.472
+ vertex 132.392 0 249.528
endloop
endfacet
- facet normal -0.5 0 0.866026
+ facet normal 0.156432 0 -0.987689
outer loop
- vertex 214 -4 378
- vertex 209 -4 386.66
- vertex 209 -4 369.34
+ vertex 132 0 252
+ vertex 132.392 0 249.528
+ vertex 132.392 0 254.472
endloop
endfacet
endsolid
diff --git a/apps/gpu/LBM/ActuatorLine/configActuatorLine.txt b/apps/gpu/LBM/ActuatorLine/configActuatorLine.txt
index 11f0d471d..f28c2ff3c 100644
--- a/apps/gpu/LBM/ActuatorLine/configActuatorLine.txt
+++ b/apps/gpu/LBM/ActuatorLine/configActuatorLine.txt
@@ -8,7 +8,7 @@
GridPath=.
##################################################
ReferenceDiameter=126
-NodesPerDiameter=32 #20=800mb # 32=3200mb
+NodesPerDiameter=40 #20=800mb # 32=3200mb
Velocity=9
##################################################
tStartOut=0
diff --git a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/FluidFlow/Compressible/CumulantK17Sponge/CumulantK17Sponge_Device.cu b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/FluidFlow/Compressible/CumulantK17Sponge/CumulantK17Sponge_Device.cu
index a6bf00fb0..ce6bbce68 100644
--- a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/FluidFlow/Compressible/CumulantK17Sponge/CumulantK17Sponge_Device.cu
+++ b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/FluidFlow/Compressible/CumulantK17Sponge/CumulantK17Sponge_Device.cu
@@ -358,8 +358,8 @@ __global__ void LB_Kernel_CumulantK17Sponge(
// 3.
// Calculate modified omega for sponge bob layer
- real startXsponge = 1060.0f;
- real endXsponge = 1260.0f;
+ real startXsponge = 656.0f;
+ real endXsponge = 756.0f;
real sizeSponge = endXsponge - startXsponge;
real OxyyPxzz = c8o1 * (-c2o1 + omega) * (c1o1 + c2o1 * omega) / (-c8o1 - c14o1 * omega + c7o1 * omega * omega);
--
GitLab