diff --git a/apps/gpu/LBM/ActuatorLine/ActuatorLine.cpp b/apps/gpu/LBM/ActuatorLine/ActuatorLine.cpp index 410e77a83df24050ebfd6a6c894ed3f8f09c0e41..62957069ee279e3ff3e93ea1c4d0203bf5c5a2f2 100644 --- a/apps/gpu/LBM/ActuatorLine/ActuatorLine.cpp +++ b/apps/gpu/LBM/ActuatorLine/ActuatorLine.cpp @@ -68,10 +68,6 @@ //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// -//////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// -//////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// -//////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// - void multipleLevel(const std::string& configPath) { @@ -100,7 +96,7 @@ void multipleLevel(const std::string& configPath) const real L_x = 6*reference_diameter; - const real L_y = 4*reference_diameter; + const real L_y = 3*reference_diameter; const real L_z = 4*reference_diameter; const real viscosity = 1.56e-5; @@ -126,21 +122,21 @@ void multipleLevel(const std::string& configPath) const real dx = reference_diameter/real(nodes_per_diameter); - real turbPos[3] = {1*reference_diameter, 2*reference_diameter, 2*reference_diameter}; + real turbPos[3] = {1*reference_diameter, 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]-0.4*reference_diameter, turbPos[1]-0.8*reference_diameter, turbPos[2]-0.8*reference_diameter, + gridBuilder->addGrid( new Cuboid( turbPos[0]-0.4*reference_diameter, turbPos[1]-1*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); @@ -207,8 +203,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->setVelocityBoundaryCondition(SideType::GEOMETRY, 0.0, 0.0, 0.0); + bcFactory.setGeometryBoundaryCondition(BoundaryConditionFactory::NoSlipBC::NoSlipCompressible); bcFactory.setVelocityBoundaryCondition(BoundaryConditionFactory::VelocityBC::VelocityAndPressureCompressible); bcFactory.setPressureBoundaryCondition(BoundaryConditionFactory::PressureBC::OutflowNonReflective); @@ -216,6 +212,8 @@ void multipleLevel(const std::string& configPath) SPtr<TurbulenceModelFactory> tmFactory = std::make_shared<TurbulenceModelFactory>(para); tmFactory->readConfigFile(config); + // gridBuilder->writeGridsToVtk(para->getOutputPath() + "Grid" + "_"); + //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// int level = 1; // grid level at which the turbine samples velocities and distributes forces diff --git a/apps/gpu/LBM/ActuatorLine/Pole.stl b/apps/gpu/LBM/ActuatorLine/Pole.stl index 9373b1f54e346426bf8e2a6607d4442cbf6d67a4..dae651651f5174393d36b139b08eaf818de9f491 100644 --- a/apps/gpu/LBM/ActuatorLine/Pole.stl +++ b/apps/gpu/LBM/ActuatorLine/Pole.stl @@ -1,534 +1,254 @@ solid Visualization Toolkit generated SLA File - facet normal 0.987689 0 -0.156432 + facet normal 0.951056 0 -0.309018 outer loop - vertex 148 252 252 - vertex 148 0 252 - vertex 147.608 252 249.528 + vertex 150 126 252 + vertex 150 0 252 + vertex 148.09 126 246.122 endloop endfacet - facet normal -0.987689 0 0.156432 + facet normal -0.951056 0 0.309018 outer loop - vertex 147.608 0 249.528 - vertex 147.608 252 249.528 - vertex 148 0 252 + vertex 148.09 0 246.122 + vertex 148.09 126 246.122 + vertex 150 0 252 endloop endfacet - facet normal 0.987689 0 -0.156432 + facet normal 0.951056 0 -0.309018 outer loop - vertex 147.608 252 249.528 - vertex 147.608 0 249.528 - vertex 146.472 252 247.298 + vertex 148.09 126 246.122 + vertex 148.09 0 246.122 + vertex 143.09 126 242.489 endloop endfacet - facet normal -0.987689 0 0.156432 + facet normal -0.951056 0 0.309018 outer loop - vertex 146.472 0 247.298 - vertex 146.472 252 247.298 - vertex 147.608 0 249.528 + vertex 143.09 0 242.489 + vertex 143.09 126 242.489 + vertex 148.09 0 246.122 endloop endfacet - facet normal 0.987689 0 -0.156432 + facet normal 0.951056 0 -0.309018 outer loop - vertex 146.472 252 247.298 - vertex 146.472 0 247.298 - vertex 144.702 252 245.528 + vertex 143.09 126 242.489 + vertex 143.09 0 242.489 + vertex 136.91 126 242.489 endloop endfacet - facet normal -0.987689 0 0.156432 + facet normal -0.951056 0 0.309018 outer loop - vertex 144.702 0 245.528 - vertex 144.702 252 245.528 - vertex 146.472 0 247.298 + vertex 136.91 0 242.489 + vertex 136.91 126 242.489 + vertex 143.09 0 242.489 endloop endfacet - facet normal 0.987689 0 -0.156432 + facet normal 0.951056 0 -0.309018 outer loop - vertex 144.702 252 245.528 - vertex 144.702 0 245.528 - vertex 142.472 252 244.392 + vertex 136.91 126 242.489 + vertex 136.91 0 242.489 + vertex 131.91 126 246.122 endloop endfacet - facet normal -0.987689 0 0.156432 + facet normal -0.951056 0 0.309018 outer loop - vertex 142.472 0 244.392 - vertex 142.472 252 244.392 - vertex 144.702 0 245.528 + vertex 131.91 0 246.122 + vertex 131.91 126 246.122 + vertex 136.91 0 242.489 endloop endfacet - facet normal 0.987689 0 -0.156432 + facet normal 0.951056 0 -0.309018 outer loop - vertex 142.472 252 244.392 - vertex 142.472 0 244.392 - vertex 140 252 244 + vertex 131.91 126 246.122 + vertex 131.91 0 246.122 + vertex 130 126 252 endloop endfacet - facet normal -0.987689 0 0.156432 + facet normal -0.951056 0 0.309018 outer loop - vertex 140 0 244 - vertex 140 252 244 - vertex 142.472 0 244.392 + vertex 130 0 252 + vertex 130 126 252 + vertex 131.91 0 246.122 endloop endfacet - facet normal 0.987689 0 -0.156432 + facet normal 0.951056 0 -0.309018 outer loop - vertex 140 252 244 - vertex 140 0 244 - vertex 137.528 252 244.392 + vertex 130 126 252 + vertex 130 0 252 + vertex 131.91 126 257.878 endloop endfacet - facet normal -0.987689 0 0.156432 + facet normal -0.951056 0 0.309018 outer loop - vertex 137.528 0 244.392 - vertex 137.528 252 244.392 - vertex 140 0 244 + vertex 131.91 0 257.878 + vertex 131.91 126 257.878 + vertex 130 0 252 endloop endfacet - facet normal 0.987689 0 -0.156432 + facet normal 0.951056 0 -0.309018 outer loop - vertex 137.528 252 244.392 - vertex 137.528 0 244.392 - vertex 135.298 252 245.528 + vertex 131.91 126 257.878 + vertex 131.91 0 257.878 + vertex 136.91 126 261.511 endloop endfacet - facet normal -0.987689 0 0.156432 + facet normal -0.951056 0 0.309018 outer loop - vertex 135.298 0 245.528 - vertex 135.298 252 245.528 - vertex 137.528 0 244.392 + vertex 136.91 0 261.511 + vertex 136.91 126 261.511 + vertex 131.91 0 257.878 endloop endfacet - facet normal 0.987689 0 -0.156432 + facet normal 0.951056 0 -0.309018 outer loop - vertex 135.298 252 245.528 - vertex 135.298 0 245.528 - vertex 133.528 252 247.298 + vertex 136.91 126 261.511 + vertex 136.91 0 261.511 + vertex 143.09 126 261.511 endloop endfacet - facet normal -0.987689 0 0.156432 + facet normal -0.951056 0 0.309018 outer loop - vertex 133.528 0 247.298 - vertex 133.528 252 247.298 - vertex 135.298 0 245.528 + vertex 143.09 0 261.511 + vertex 143.09 126 261.511 + vertex 136.91 0 261.511 endloop endfacet - facet normal 0.987689 0 -0.156432 + facet normal 0.951056 0 -0.309018 outer loop - vertex 133.528 252 247.298 - vertex 133.528 0 247.298 - vertex 132.392 252 249.528 + vertex 143.09 126 261.511 + vertex 143.09 0 261.511 + vertex 148.09 126 257.878 endloop endfacet - facet normal -0.987689 0 0.156432 + facet normal -0.951056 0 0.309018 outer loop - vertex 132.392 0 249.528 - vertex 132.392 252 249.528 - vertex 133.528 0 247.298 + vertex 148.09 0 257.878 + vertex 148.09 126 257.878 + vertex 143.09 0 261.511 endloop endfacet - facet normal 0.987689 0 -0.156432 + facet normal 0.951056 0 -0.309018 outer loop - vertex 132.392 252 249.528 - vertex 132.392 0 249.528 - vertex 132 252 252 + vertex 148.09 126 257.878 + vertex 148.09 0 257.878 + vertex 150 126 252 endloop endfacet - facet normal -0.987689 0 0.156432 + facet normal -0.951056 0 0.309018 outer loop - vertex 132 0 252 - vertex 132 252 252 - vertex 132.392 0 249.528 + vertex 150 0 252 + vertex 150 126 252 + vertex 148.09 0 257.878 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 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 + facet normal 0.309017 0 -0.951056 outer loop - vertex 142.472 252 259.608 - vertex 144.702 252 258.472 - vertex 144.702 252 245.528 + vertex 150 126 252 + vertex 148.09 126 246.122 + vertex 148.09 126 257.878 endloop endfacet - facet normal 0.156428 0.00944581 -0.987644 + facet normal -0.308774 0.0396513 0.950309 outer loop - vertex 144.702 252 258.472 - vertex 146.472 252 256.702 - vertex 144.702 252 245.528 + vertex 148.09 126 246.122 + vertex 143.09 126 242.489 + vertex 148.09 126 257.878 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 + vertex 148.09 126 257.878 + vertex 143.09 126 242.489 + vertex 143.09 126 261.511 endloop endfacet - facet normal -0.156428 0.00944581 0.987644 + facet normal 0 0.028819 -0.999585 outer loop - vertex 146.472 252 247.298 - vertex 147.608 252 254.472 - vertex 147.608 252 249.528 + vertex 143.09 126 261.511 + vertex 143.09 126 242.489 + vertex 136.91 126 261.511 endloop endfacet - facet normal 0.156434 0 -0.987688 + facet normal -0.309017 0 0.951057 outer loop - vertex 148 252 252 - vertex 147.608 252 249.528 - vertex 147.608 252 254.472 + vertex 136.91 126 261.511 + vertex 143.09 126 242.489 + vertex 131.91 126 257.878 endloop endfacet - facet normal 0.156434 0 -0.987688 + facet normal 0.308981 0.0151556 -0.950947 outer loop - vertex 147.608 0 254.472 - vertex 146.472 0 256.702 - vertex 148 0 252 + vertex 131.91 126 257.878 + vertex 143.09 126 242.489 + vertex 130 126 252 endloop endfacet - facet normal -0.156434 0 0.987688 + facet normal -0.587785 0 0.809018 outer loop - vertex 148 0 252 - vertex 146.472 0 256.702 - vertex 147.608 0 249.528 + vertex 130 126 252 + vertex 143.09 126 242.489 + vertex 131.91 126 246.122 endloop endfacet - facet normal 0.309002 0.00980961 -0.951011 + facet normal 0.587785 0 -0.809018 outer loop - vertex 147.608 0 249.528 - vertex 146.472 0 256.702 - vertex 146.472 0 247.298 + vertex 136.91 126 242.489 + vertex 131.91 126 246.122 + vertex 143.09 126 242.489 endloop endfacet - facet normal -0.309017 0 0.951056 + facet normal -0.587785 0 0.809018 outer loop - vertex 146.472 0 247.298 - vertex 146.472 0 256.702 - vertex 144.702 0 245.528 + vertex 130 0 252 + vertex 131.91 0 246.122 + vertex 131.91 0 257.878 endloop endfacet - facet normal 0.45397 0.00944579 -0.890967 + facet normal 0.58754 0.028819 -0.808682 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 + vertex 131.91 0 246.122 + vertex 136.91 0 242.489 + vertex 131.91 0 257.878 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 + vertex 131.91 0 257.878 + vertex 136.91 0 242.489 + vertex 136.91 0 261.511 endloop endfacet - facet normal 0.156435 0 -0.987688 + facet normal -0.308774 0.0396513 0.950309 outer loop - vertex 135.298 0 258.472 - vertex 133.528 0 256.702 - vertex 135.298 0 245.528 + vertex 136.91 0 261.511 + vertex 136.91 0 242.489 + vertex 143.09 0 261.511 endloop endfacet - facet normal 0 0.00450917 0.99999 + facet normal 0.587785 0 -0.809017 outer loop - vertex 133.528 0 256.702 - vertex 132.392 0 254.472 - vertex 135.298 0 245.528 + vertex 143.09 0 261.511 + vertex 136.91 0 242.489 + vertex 148.09 0 257.878 endloop endfacet - facet normal 0 0.00155373 -0.999999 + facet normal -0.587147 0.0465989 0.808138 outer loop - vertex 135.298 0 245.528 - vertex 132.392 0 254.472 - vertex 133.528 0 247.298 + vertex 148.09 0 257.878 + vertex 136.91 0 242.489 + vertex 150 0 252 endloop endfacet - facet normal -0.156432 0 0.987689 + facet normal 0.587785 0 0.809018 outer loop - vertex 133.528 0 247.298 - vertex 132.392 0 254.472 - vertex 132.392 0 249.528 + vertex 150 0 252 + vertex 136.91 0 242.489 + vertex 148.09 0 246.122 endloop endfacet - facet normal 0.156432 0 -0.987689 + facet normal -0.587785 0 -0.809018 outer loop - vertex 132 0 252 - vertex 132.392 0 249.528 - vertex 132.392 0 254.472 + vertex 143.09 0 242.489 + vertex 148.09 0 246.122 + vertex 136.91 0 242.489 endloop endfacet endsolid diff --git a/apps/gpu/LBM/ActuatorLine/configActuatorLine.txt b/apps/gpu/LBM/ActuatorLine/configActuatorLine.txt index f28c2ff3c6254c033d037200f8588becb3f73f16..ae4ea633a164b6d8e9b09b7cb754cda4f18a1e60 100644 --- a/apps/gpu/LBM/ActuatorLine/configActuatorLine.txt +++ b/apps/gpu/LBM/ActuatorLine/configActuatorLine.txt @@ -8,7 +8,7 @@ GridPath=. ################################################## ReferenceDiameter=126 -NodesPerDiameter=40 #20=800mb # 32=3200mb +NodesPerDiameter=40 Velocity=9 ################################################## tStartOut=0