diff --git a/apps/gpu/LBM/ActuatorLine/ActuatorLine.cpp b/apps/gpu/LBM/ActuatorLine/ActuatorLine.cpp index bb5572cacaeb6947a82d5b396157bab8acfede76..b19557ac2c86a596a7bcb2392ba8e2cb2e58e5ba 100644 --- a/apps/gpu/LBM/ActuatorLine/ActuatorLine.cpp +++ b/apps/gpu/LBM/ActuatorLine/ActuatorLine.cpp @@ -149,6 +149,8 @@ void multipleLevel(const std::string& configPath) const real viscosityLB = viscosity * dt / (dx * dx); // LB units + VF_LOG_INFO("Knoten pro Turbinendurchmesser = {}", reference_diameter/dx); + VF_LOG_INFO("velocity [m/s] = {}", velocity); VF_LOG_INFO("velocity [dx/dt] = {}", velocityLB); VF_LOG_INFO("viscosity [10^8 dx^2/dt] = {}", viscosityLB*1e8); @@ -174,9 +176,13 @@ void multipleLevel(const std::string& configPath) vz = (real)0.0; }); - para->setTimestepStartOut( uint(tStartOut/dt) ); - para->setTimestepOut( uint(tOut/dt) ); - para->setTimestepEnd( uint(tEnd/dt) ); + para->setTimestepStartOut( uint(tStartOut) ); + para->setTimestepOut( uint(tOut) ); + para->setTimestepEnd( uint(tEnd) ); + + // para->setTimestepStartOut( uint(tStartOut/dt) ); + // para->setTimestepOut( uint(tOut/dt) ); + // para->setTimestepEnd( uint(tEnd/dt) ); para->setIsBodyForce( true ); para->setUseStreams( true ); @@ -184,7 +190,7 @@ void multipleLevel(const std::string& configPath) ///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// gridBuilder->setVelocityBoundaryCondition(SideType::MX, velocityLB, 0.0, 0.0); - gridBuilder->setVelocityBoundaryCondition(SideType::MY, 0.0 , 0.0, 0.0); + gridBuilder->setVelocityBoundaryCondition(SideType::MY, velocityLB, 0.0, 0.0); gridBuilder->setVelocityBoundaryCondition(SideType::PY, velocityLB, 0.0, 0.0); gridBuilder->setVelocityBoundaryCondition(SideType::MZ, velocityLB, 0.0, 0.0); gridBuilder->setVelocityBoundaryCondition(SideType::PZ, velocityLB, 0.0, 0.0); @@ -214,7 +220,6 @@ void multipleLevel(const std::string& configPath) actuator_farm->addTurbine(turbPos[0], turbPos[1], turbPos[2], reference_diameter, omega, 0, 0, bladeRadii); para->addActuator( actuator_farm ); - // SPtr<PointProbe> pointProbe = std::make_shared<PointProbe>("pointProbe", para->getOutputPath(), 100, 1, 500, 100); // std::vector<real> probeCoordsX = {reference_diameter,2*reference_diameter,5*reference_diameter}; // std::vector<real> probeCoordsY = {3*reference_diameter,3*reference_diameter,3*reference_diameter}; diff --git a/apps/gpu/LBM/ActuatorLine/configActuatorLine.txt b/apps/gpu/LBM/ActuatorLine/configActuatorLine.txt index 3f94f328a34638750ff142e66d161da0f241fe4b..e6475bbae26ba816ed8d4bbd1d489613e4747f6d 100644 --- a/apps/gpu/LBM/ActuatorLine/configActuatorLine.txt +++ b/apps/gpu/LBM/ActuatorLine/configActuatorLine.txt @@ -11,7 +11,7 @@ ReferenceDiameter=126 NodesPerDiameter=20 Velocity=9 ################################################## -tStartOut=100 +tStartOut=0 tOut=50 tEnd=500 ##################################################