diff --git a/src/VirtualFluids_GPU/DataStructureInitializer/GridProvider.cpp b/src/VirtualFluids_GPU/DataStructureInitializer/GridProvider.cpp index 7f6a2a15e55d4963c33669422bb6325af8221a42..f3e83c03cf099c4a804f52671a1049353c70b641 100644 --- a/src/VirtualFluids_GPU/DataStructureInitializer/GridProvider.cpp +++ b/src/VirtualFluids_GPU/DataStructureInitializer/GridProvider.cpp @@ -48,11 +48,11 @@ void GridProvider::setInitalNodeValues(const int numberOfNodes, const int level) const real coordZ = para->getParH(level)->coordZ_SP[j]; const real velocity = para->getVelocity(); - para->getParH(level)->rho_SP[j] = real((velocity * velocity) * 3.0 / 4.0 * (cos(coordX * 4.0*PI / gridX) + cos(coordZ * 4.0*PI / gridZ))) * gridZ / gridX; + para->getParH(level)->rho_SP[j] = real(0.0); //real((velocity * velocity) * 3.0 / 4.0 * (cos(coordX * 4.0*PI / gridX) + cos(coordZ * 4.0*PI / gridZ))) * gridZ / gridX; para->getParH(level)->vy_SP[j] = real(0.0); - para->getParH(level)->vx_SP[j] = real( velocity * sin(coordX * 2.0*PI / gridX) * cos(coordZ * 2.0*PI / gridZ)) + uAdvect * (1.0 + para->getParH(level)->rho_SP[j]); - para->getParH(level)->vz_SP[j] = real(-velocity * cos(coordX * 2.0*PI / gridX) * sin(coordZ * 2.0*PI / gridZ)); // *(real)(gridZ) / (real)(gridX); + para->getParH(level)->vx_SP[j] = real(0.0); //real( velocity * sin(coordX * 2.0*PI / gridX) * cos(coordZ * 2.0*PI / gridZ)) + uAdvect * (1.0 + para->getParH(level)->rho_SP[j]); + para->getParH(level)->vz_SP[j] = real(0.0); //real(-velocity * cos(coordX * 2.0*PI / gridX) * sin(coordZ * 2.0*PI / gridZ)); // *(real)(gridZ) / (real)(gridX); //para->getParH(level)->vx_SP[j] = 0.001; //para->getParH(level)->vy_SP[j] = 0.0; diff --git a/src/VirtualFluids_GPU/LBM/Simulation.cpp b/src/VirtualFluids_GPU/LBM/Simulation.cpp index 7d2a5d00fb8f90a5cdfaed8404fe5922fe963482..062ead70a3805706eaacf71dfd23ba3af9cbdc38 100644 --- a/src/VirtualFluids_GPU/LBM/Simulation.cpp +++ b/src/VirtualFluids_GPU/LBM/Simulation.cpp @@ -807,13 +807,13 @@ void Simulation::run() //para->getParD(0)->neighborX_SP, para->getParD(0)->neighborY_SP, para->getParD(0)->neighborZ_SP, //para->getParD(0)->size_Mat_SP, para->getParD(0)->evenOrOdd); // getLastCudaError("QVelDevComp27 execution failed"); - //QVelDevCompZeroPress27(para->getParD(0)->numberofthreads, para->getParD(0)->nx, para->getParD(0)->ny, - // para->getParD(0)->Qinflow.Vx, para->getParD(0)->Qinflow.Vy, para->getParD(0)->Qinflow.Vz, - // para->getParD(0)->d0SP.f[0], para->getParD(0)->Qinflow.k, para->getParD(0)->Qinflow.q27[0], - // para->getParD(0)->kInflowQ, para->getParD(0)->Qinflow.kArray, para->getParD(0)->omega, - // para->getParD(0)->neighborX_SP, para->getParD(0)->neighborY_SP, para->getParD(0)->neighborZ_SP, - // para->getParD(0)->size_Mat_SP, para->getParD(0)->evenOrOdd); - //getLastCudaError("QVelDevComp27 execution failed"); + QVelDevCompZeroPress27(para->getParD(0)->numberofthreads, para->getParD(0)->nx, para->getParD(0)->ny, + para->getParD(0)->Qinflow.Vx, para->getParD(0)->Qinflow.Vy, para->getParD(0)->Qinflow.Vz, + para->getParD(0)->d0SP.f[0], para->getParD(0)->Qinflow.k, para->getParD(0)->Qinflow.q27[0], + para->getParD(0)->kInflowQ, para->getParD(0)->Qinflow.kArray, para->getParD(0)->omega, + para->getParD(0)->neighborX_SP, para->getParD(0)->neighborY_SP, para->getParD(0)->neighborZ_SP, + para->getParD(0)->size_Mat_SP, para->getParD(0)->evenOrOdd); + getLastCudaError("QVelDevComp27 execution failed"); //////////////////////////////////////////////////////////////////////////// //QVeloDevEQ27(para->getParD(0)->numberofthreads, diff --git a/targets/apps/LBM/DrivenCavity/configDrivenCavity.txt b/targets/apps/LBM/DrivenCavity/configDrivenCavity.txt index 5b0aa6de2da4a59a2fc0cf7a3b3f42e2ae8d305e..a6793f0b389330ba814b5bc183e1d85e8323aacf 100644 --- a/targets/apps/LBM/DrivenCavity/configDrivenCavity.txt +++ b/targets/apps/LBM/DrivenCavity/configDrivenCavity.txt @@ -24,11 +24,11 @@ NOGL=1 #LBM Version ################################################## D3Qxx=27 -MainKernelName=CumulantAA2016CompSP27 +#MainKernelName=CumulantAA2016CompSP27 ################################################## #simulation parameter ################################################## -TimeEnd=10000 +TimeEnd=100000 TimeOut=1000 TimeStartOut=0 \ No newline at end of file