diff --git a/source/Applications/pChannel/configBombadilpChannel.cfg b/source/Applications/pChannel/configBombadilpChannel.cfg index 905962b71783cfc90870ae21e53e7017241bebc9..69f20a844bec121fb29b5059fb35a49b21c81925 100644 --- a/source/Applications/pChannel/configBombadilpChannel.cfg +++ b/source/Applications/pChannel/configBombadilpChannel.cfg @@ -2,9 +2,9 @@ #Simulation parameters for porous channel # -pathOut = d:/temp/ChannelFlow-test +pathOut = d:/temp/ChannelFlow pathGeo = d:/Projects/SFB880/GeometrienPoroeseMedien/isotrop/PA80-110 -numOfThreads = 1 +numOfThreads = 4 availMem = 14e9 logToFile = false @@ -83,19 +83,20 @@ Re = 15000 #real velocity is 54.95 m/s u_LB = 0.1 -newStart = true -restartStep = 10 +newStart = false +restartStep = 230000 -cpStep = 10 -cpStart = 10 +cpStep = 100 +cpStart = 100 averaging = false averagingReset = false timeAvStart = 21000000 timeAvStop = 2100010000 -endTime = 100 -outTime = 10000 +outTime = 100 +endTime = 230000 + -nupsStep = 100 100 10000000 +nupsStep = 10 10 10000000 diff --git a/source/Applications/pChannel/pChannel.cpp b/source/Applications/pChannel/pChannel.cpp index 62a2d405a2be28c5465060805063797fb0db21df..2dad903437de6554ccc7b69aae1486cf22c035b1 100644 --- a/source/Applications/pChannel/pChannel.cpp +++ b/source/Applications/pChannel/pChannel.cpp @@ -505,6 +505,7 @@ void run(string configname) calculator->addCoProcessor(nupsCoProcessor); calculator->addCoProcessor(AdjForcCoProcessor); calculator->addCoProcessor(migCoProcessor); + //calculator->addCoProcessor(restartCoProcessor); calculator->addCoProcessor(writeMQSelectCoProcessor); calculator->addCoProcessor(writeMQCoProcessor); diff --git a/source/VirtualFluidsCore/Visitors/SetForcingBlockVisitor.cpp b/source/VirtualFluidsCore/Visitors/SetForcingBlockVisitor.cpp index 88ca54934de4fc7d7ea3ce21f696edd365cd445f..100ee2fe499f691484017f60654ad922be3f9d5d 100644 --- a/source/VirtualFluidsCore/Visitors/SetForcingBlockVisitor.cpp +++ b/source/VirtualFluidsCore/Visitors/SetForcingBlockVisitor.cpp @@ -32,12 +32,12 @@ SetForcingBlockVisitor::SetForcingBlockVisitor(const std::string& sForcingX1, co ////////////////////////////////////////////////////////////////////////// void SetForcingBlockVisitor::visit(SPtr<Grid3D> grid, SPtr<Block3D> block) { - SPtr<LBMKernel> kernel = dynamicPointerCast<LBMKernel>(block->getKernel()); - if (!kernel) - throw std::runtime_error("SetForcingBlockVisitor: Kernel is not a LBMKernel"); - if(block->getRank() == grid->getRank()) { + SPtr<LBMKernel> kernel = dynamicPointerCast<LBMKernel>(block->getKernel()); + if (!kernel) + throw UbException(UB_EXARGS, "LBMKernel is not exist"); + switch (ftype) { case 0: